应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 如何让计算离子从晶体内部扩散到表面的能量？
152/2返回列表上一页12
查看: 3106  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

【答案】应助回帖

★ ★ ★ ★ ★
franch: 金币+5, 应助指数+1, 谢谢回帖交流。。。 2014-03-31 22:26:03
引用回帖:
10楼: Originally posted by 1164028932 at 2014-03-30 14:16:14
麻烦你看一下，我根据你的意思建立的反应物与生成物的结构如下图所示，你看一下合理吗?对不对？还有就是优化的时候要固定晶格优化吗?优化时固定其他原子只优化表面原子还是都优化？

reaction.png

production.png
...

OB轴向取的太少了。建议把OB轴向取得再大一些，否则相当于OB向的一条线上的Ba都脱出，这不合理。另外真空层建议至少10A，如果严格的话建议做测试，CASTEP的计算量和晶格大小有关，所以要看能承受多大的计算量了，如果是Dmol3就无所谓，真空层取个三五十A问题都不大。

优化应当作不opt cell（但是因此也有可能造成不收敛）的全几何优化，不要固定，这样拖出前后晶格弛豫的能量会包含在脱出的能量中，也很容易理解，离子脱出后晶格需要重组、再构。不过也有可能会不容易收敛。

不过我估计CASTEP计算这个体系有点危险，速度相当慢（我估计单节点每个structure都要几个小时），因为大致上从尺寸上看K-point基本是gamma了，这样CASTEP就算跨节点也不会有好的效率（CASTEP计算划分是依据K-point的）。速度如果比较慢就不要用Ultrafine，可以降低到fine。另外可以尝试Dmol3（也有可能收敛造成问题）甚至开源的SIESTA或Openmx。差不多的精确程度来看，VASP速度比CASTEP好得多，这个你就要好好请教楼上前辈了。

老菜鸟的一点建议，你先试试吧。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
11楼2014-03-30 16:02:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1164028932

金虫 (初入文坛)
引用回帖:
11楼: Originally posted by 卡开发发 at 2014-03-30 16:02:13
OB轴向取的太少了。建议把OB轴向取得再大一些，否则相当于OB向的一条线上的Ba都脱出，这不合理。另外真空层建议至少10A，如果严格的话建议做测试，CASTEP的计算量和晶格大小有关，所以要看能承受多大的计算量了，如 ...

恩好的，不要谦虚，在我眼里你已经是大神了，我还是请教一下，我这个能用Dmol算过渡态吗?我前几天就是这个结构分别用CASTEP和Dmol计算了同一个路径的扩散能，参数设置基本一样，为什么结果差距很大呀？我看文献没发现有用Dmol算晶体内部离子扩散的?麻烦大神指导一下。谢谢。
赞一下 回复此楼
12楼2014-03-31 21:28:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
12楼: Originally posted by 1164028932 at 2014-03-31 21:28:52
恩好的，不要谦虚，在我眼里你已经是大神了，我还是请教一下，我这个能用Dmol算过渡态吗?我前几天就是这个结构分别用CASTEP和Dmol计算了同一个路径的扩散能，参数设置基本一样，为什么结果差距很大呀？我看文献没发 ...

原子轨道的方法计算固体是不太准确，优化晶格结构都有点困难，castep对于PBC体系肯定好些的。你指的结果差距指的是扩散能垒？除了平面波方法中的pseudopot不一样，另外k-point和Ecut一般需要测试我们是知道的，Dmol3中的k-point以及原子轨道的截断也是要测试的，这些都不一样结果不一样也很正常。不过我觉得你这个体系偏大，CASTEP算起来有点危险，完全准确不可能，趋势正确，计算参数有力有据即可。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
13楼2014-04-01 09:37:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1164028932

金虫 (初入文坛)
引用回帖:
13楼: Originally posted by 卡开发发 at 2014-04-01 09:37:04
原子轨道的方法计算固体是不太准确，优化晶格结构都有点困难，castep对于PBC体系肯定好些的。你指的结果差距指的是扩散能垒？除了平面波方法中的pseudopot不一样，另外k-point和Ecut一般需要测试我们是知道的，Dmo ...

恩就是的，扩散能垒差距很大，你说的PBC体系是什么体系？
赞一下 回复此楼
14楼2014-04-02 22:35:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
14楼: Originally posted by 1164028932 at 2014-04-02 22:35:31
恩就是的，扩散能垒差距很大，你说的PBC体系是什么体系？...

周期边界条件（periodic boundary condition）。我觉得你先看看趋势（和试验方面），比较纯粹的绝对值意义不大，这个本身也依赖泛函的选取，很显然的就是即便同一种泛函采用不同基函数展开，特性也很不同。PBC下Dmol3计算某些性质我听说不是很靠谱。

castep大家很注意Ecut的选取，但是Dmol3大家对Rcut的选取似乎并不是很重视，两者的可比性我觉得都只有这些相关的参数极限外推才能够作比较。

另外，我不确定的是你现在过渡态构型上两者做出来的定性是否正确？Dmol3可以在TS search后计算Freq，然后进行TS opt，经过优化后的构型是否比较过？

靠谱点的话还是VASP。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
15楼2014-04-02 23:00:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1164028932 的主题更新
152/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 今天基金会出结果吗？20260819 +13 kkkl_v 2026-08-19 13/650 2026-08-19 22:37 by 宝贝虫子
[基金申请] 时间戳变了，能看出什么问题？ +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 重要消息，中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
关闭
请填处理意见
关闭