| 查看: 2139 | 回复: 3 | |||
[求助]
求助：VASP计算广义层错能GSFE的方法问题是否有误？
|
|
各位： 看到文献用VASP/CASTEP计算GSFE的文献， 都是人为优化一个结构后构建10-12层的单胞外加一个真空层来处理，一般在40个原子左右.(BCC,FCC )。然后以中间原子基线，让上部原子整体沿滑移方向平移0.1b, 0.2b,0.3b....计算构型的能量，以各个构型的能量差值，除以滑移平面的面积，就是层错能的变化了，极大值 一般在0.5b（BCC。 也有很多几个波浪形式，可以得到稳定、不稳定层错能、孪晶能等） 问题是： （1）请问是计算各个构型的总能量吗？ 如果是单质计算到还可以，如果是做doping，是否要考虑结合能或Mixing能，要扣除掉么？ 我看的文献上没说也，介绍了计算方法，然后直接上图... （2）需要对全部原子进行几何优化（至少在平移层附近两层原子要优化，其它固定？）， 以优化后的结构总能量计算？ 但几何优化后，原子坐标是否变了（比如doping了个头更大的元素作为置换原子），那就不是0.1b, 0.2b 的等距平移了吧？ 如果我只看规律，不做几何优化可否？ 感觉不做几何优化不行.... （3）按照文献方法，我以BCC单质，我建了个40原子的构型，模拟[111](110)滑移系。先以纯单质平移0.1-0.5b，用CASTEP粗略跑了下，曲线关于0.5b对称，而不像很多其它金属的都是有极小、极大值。 但我单机内存不足....于是： 后来用VASP跑。 我以ISIF=2的跑了一下，算了半个月了....居然才算1步（也就是只算了0.1b结构的能量，距离跑到0.5差远了...），我就stop了。 然后我去掉ISIF=2, 直接计算doping一个原子后的总能量，跑了快半个月了，还没算完.，也才1步... 这样下来，做不同元素、doing不同浓度，两年也算不完...... 但我计算纯元素的时候，速度还可以，不到三个星期跑完（汗呀... 相对快吧）。。 为什么doping了一个原子，慢了N倍，而且没做ISIF... （服务器：2节点，共32核） ！！！！！！请斧正！！！！！！！ [ Last edited by 鲤鱼飞 on 2013-6-17 at 21:10 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 科研必备知识 |
» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有87人回复
-
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 六角密排非基面广义层错能计算
已经有12人回复
2楼2013-06-17 21:30:33
ustbwgs
木虫 (正式写手)
- 应助: 14 (小学生)
- 金币: 1449.6
- 散金: 269
- 红花: 12
- 帖子: 720
- 在线: 165.5小时
- 虫号: 1141993
- 注册: 2010-11-08
- 专业: 凝聚态物性 II ：电子结构
3楼2013-10-16 22:09:11
shuoshuohome
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 51.8
- 帖子: 38
- 在线: 42.5小时
- 虫号: 2540055
- 注册: 2013-07-10
- 专业: 金属材料的合金相、相变及
4楼2015-01-28 15:16:03
