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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】我对表面驰豫,算完了原子却一个都没有动
51/1返回列表
查看: 2027  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-11-21 19:42:14
已阅   回复此楼   关注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
赞一下(1人) 回复此楼
3楼2010-11-21 20:27:15
已阅   回复此楼   关注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，所以肯定是计算一步之后立即停止计算，所以离子不会发生移动
赞一下(1人) 回复此楼
2楼2010-11-21 19:53:17
已阅   回复此楼   关注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的回帖

stractor

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指教 2010-11-22 18:31:17
表面优化ISIF=2合理一些吧。
我是从全驰豫的体材料（ISIF=3)构建表面的。也可以这样理解：表面是个二维结构，只有一个方向不具有周期性，其它两个方向（最下面几层）都还是和体材料的环境一样的，它们还应该和体材料的晶格常数一致，所以驰豫的时候就ISIF=2了。

[ Last edited by stractor on 2010-11-22 at 15:12 ]
赞一下(2人) 回复此楼
5楼2010-11-22 15:07:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭