24小时热门版块排行榜    

查看: 2029  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

liuhengjin

银虫 (小有名气)


[交流] 【求助】我对表面驰豫,算完了原子却一个都没有动

这是我的INCAR文件

SYSTEM = Ge(100)_(2x1)
ENCUT =400
ISMEAR = 0; SIGMA = 0.05
ISTART = 0; ICHARG = 2
GGA=PE
ALGO =Fast
EDIFF = 1E-5; EDIFFG = -1.0E-2
NELMIN=6
ISIF=3
IBRION=2
LREAL=AUTO

POSCAR为
Ge(110)_(2x1)                           
   1.02200000000000     
     4.0004566145628928    0.0000000000000000    0.0000000000000000
     0.0000000000000000    8.0009132291257874    0.0000000000000000
     0.0000000000000000    0.0000000000000000   29.1437500000000007
  22
Selective dynamics
Direct
  0.0000000000000000  0.0000000000000000  0.0970619772679981   F   F   F
  0.0000000000000000  0.0019822828384809  0.2967260685518983   T   T   T
  0.0000000000000000  0.9276824881117479  0.4932595612465049   T   T   T
  0.5000000000000000  0.2500000000000000  0.0000000000000000   F   F   F
  0.5000000000000000  0.2500000000000000  0.1971038563673539   T   T   T
  0.5000000000000000  0.2500000000000000  0.3964169710606655   T   T   T
  0.5000000000000000  0.0000000000000000  0.0485309886339991   F   F   F
  0.5000000000000000 -0.0029335760595449  0.2466868371805478   T   T   T
  0.5000000000000000  0.0058225771951400  0.4521603543234448   T   T   T
  0.0000000000000000  0.2500000000000000  0.1455929659000006   F   F   F
  0.0000000000000000  0.2500000000000000  0.3481028386587142   T   T   T
  0.0000000000000000  0.5000000000000000  0.0970619772679981   F   F   F
  0.0000000000000000  0.4980177171615190  0.2967260685518983   T   T   T
  0.0000000000000000  0.5723175118862555  0.4932595612465049   T   T   T
  0.5000000000000000  0.7500000000000000  0.0000000000000000   F   F   F
  0.5000000000000000  0.7500000000000000  0.1964311093813827   T   T   T
  0.5000000000000000  0.7500000000000000  0.3954745052361157   T   T   T
  0.5000000000000000  0.5000000000000000  0.0485309886339991   F   F   F
  0.5000000000000000  0.5029335760595451  0.2466868371805478   T   T   T
  0.5000000000000000  0.4941774228048601  0.4521603543234448   T   T   T
  0.0000000000000000  0.7500000000000000  0.1455929659000006   F   F   F
  0.0000000000000000  0.7500000000000000  0.3478604452255454   T   T   T



请大家指导,应该怎么设置?
我要驰豫的表面是Ge(100)_(2x1)
谢谢了

[ Last edited by liuhengjin on 2010-11-21 at 19:43 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

liuhengjin

银虫 (小有名气)


感谢各位的回答,谢谢。不过就6楼的观点,我个人认为,不应该先用ISIF=3来全驰豫,再用ISIF=2,这样会花掉很多时间,因为做全驰豫时底下的原子还是没有变,只是上层的原子改变了

[ Last edited by liuhengjin on 2010-11-23 at 13:39 ]
7楼2010-11-23 13:26:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

bingmou

金虫 (著名写手)


★ ★
mazuju028(金币+2):谢谢提示 2010-11-21 22:36:43
liuhengjin(金币+1):谢谢提醒 2010-11-22 13:58:28
你的NSW打开了吗?默认是NSW=0,也就是离子步为0,所以肯定是计算一步之后立即停止计算,所以离子不会发生移动
2楼2010-11-21 19:53:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)


★ ★
mazuju028(金币+2):谢谢交流 2010-11-21 22:37:46
表面优化你还用ISIF=3 难道你有优化晶胞的体积和形状?


本文来自: 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=2618313&fpage=1
3楼2010-11-21 20:27:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhengjin

银虫 (小有名气)


2楼是对的,不过我还想就3楼的疑问再问下,表面优化不可以设置isif=3吗?这个应该比isif=2更好吧?

[ Last edited by liuhengjin on 2010-11-22 at 13:57 ]
4楼2010-11-22 13:46:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见