24小时热门版块排行榜    

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

vaspnew

金虫 (小有名气)

[交流] 【求助】优化晶格常数问题 已有2人参与

我在优化一个12个原子(两种原子)的体系,a=b 不等于c,a和b成60度夹角。
我的思路是这样的:
在INCAR中设置
IBRION=2
ISIF=2
EDIFF=1E-4
EDIFFG=-0.01
用脚本对不同的晶格常数进行计算,能量最低的就是优化后的晶格常数。

我的疑问是:
两类原子的键长我都按照相同了,实际上应该略有差别,那么
上述方法是否正确?
在这里设置ISIF=2 是否合适?
我优化后的数据:
7        -111.46234
7.02        -111.47876
7.04        -111.48719
7.06        -109.99737
7.08        -111.48075
7.1        -111.46591
7.12        -111.43544
7.14        -111.41768
7.16        -111.38364
7.18        -111.34314
7.2        -111.29591

其中7.06那个很怪异!!

[ Last edited by vaspnew on 2010-7-11 at 16:28 ]
回复此楼

» 猜你喜欢

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

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

vaspnew

金虫 (小有名气)

引用回帖:
Originally posted by cenwanglai at 2010-07-11 15:53:33:
可以。
用IBRION=-1也行,优化得到晶格常数后在做一次ISIF-3的优化。

我的结果感觉有问题
不知道怎么回事

7        -111.46234
7.02        -111.47876
7.04        -111.48719
7.06        -109.99737
7.08        -111.48075
7.1        -111.46591
7.12        -111.43544
7.14        -111.41768
7.16        -111.38364
7.18        -111.34314
7.2        -111.29591
3楼2010-07-11 16:29:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主


hedaors(金币+1):谢谢回复! 2010-07-11 16:03:41
vaspnew(金币+2):3q 2010-07-11 16:11:21
可以。
用IBRION=-1也行,优化得到晶格常数后在做一次ISIF-3的优化。
2楼2010-07-11 15:53:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见