24小时热门版块排行榜

返回列表

鲤鱼飞

铁虫 (正式写手)

应助: 4 (幼儿园)
金币: 710.4
散金: 15
红花: 4
帖子: 432
在线: 88.1小时
虫号: 698635
注册: 2009-02-09
性别: GG
专业: 金属结构材料

[求助] 求助：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 ]

回复此楼

CAEP

1楼2013-06-17 20:59:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

鲤鱼飞

铁虫 (正式写手)

应助: 4 (幼儿园)
金币: 710.4
散金: 15
红花: 4
帖子: 432
在线: 88.1小时
虫号: 698635
注册: 2009-02-09
性别: GG
专业: 金属结构材料

这是我的INCAR （未作ISIF的，计算也超慢....）：

SYSTEM= metal
ENCUT=400
ISTART=0
ICHARG=2
ISMEAR=-5
SIGMA=0.2
EDIFF=1E-5
PREC=Accurate
ALGO=Very_Fast
IALGO=48
因为我认为体系大，所有加了后两个。不写也可以吧...