24小时热门版块排行榜    

查看: 2031  |  回复: 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的回帖

ustbwgs

木虫 (正式写手)



小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by liuhengjin at 2010-11-21 19:42:14:
这是我的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 ...

为什么要每隔两个原子设置一个F  F  F呢
8楼2010-11-23 14:38:19
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见