24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

leongxf

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 464.5
散金: 1
帖子: 30
在线: 6.2小时
虫号: 490352
注册: 2008-01-03
专业: 金属材料的合金相、相变及

[交流] [求助]vasp计算单个Ni原子能量不收敛问题

以下是计算使用的POSCAR，和INCAR
POSCAR：
Ni
1
10.0 0.0 0.0
0.0 10.0 0.0
0.0 0.0 10.0 ！周期取10A
1
Cartesian
5.00 5.00 5.00 ！原子在原胞中心

INCAR：
System = fcc Ni
ISTART = 0 ; ICHARG=2
PREC= n
ISMEAR = 0; SIGMA = 0.01;
ISPIN=1
IBRION=-1
NSW=0
NGX=70
NGY=70
NGZ=70

K点取
Automatic mesh
0
Monkhorst Pack
1 1 1
0 0 0

Iteration 60次，没办法收敛。请问大家应该怎样修改文件，新手请教。

[ Last edited by wuchenwf on 2009-6-22 at 22:11 ]

回复此楼

» 猜你喜欢

26年申博自荐-计算机视觉已经有4人回复
考博已经有4人回复
药化及相关博士的申请已经有3人回复
一篇MDPI论文改变了学习工作和生活已经有4人回复
一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇已经有4人回复
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业已经有3人回复
收到国自然专家邀请后几年才会有本子送过来评已经有3人回复

1楼 2008-06-27 19:47:32

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

leongxf

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 464.5
散金: 1
帖子: 30
在线: 6.2小时
虫号: 490352
注册: 2008-01-03
专业: 金属材料的合金相、相变及

计算结果

采取各位的建议，模拟了很多次，终于收敛了，也找到了不收敛的原因，分析一下，请各位多指教。
1. 刚开始的KEPOINTS，取1 1 1，是因为看到论坛上有贴的计算H原子能的例子，我test了这个列子，发现 1 1 1是收敛的，故Ni原子的也取一样的，后来不收敛，我把k点增加到10 10 10，还是不收敛，所以我分析可能是INCAR的参数设置问题。
根据楼上建议，分别尝试以下的INCAR（以下所有k点都为10 10 10 ）
1. 结果不收敛（使用PW91势）
ISTART = 0 ; ICHARG=2
PREC= n
ISMEAR = -5; SIGMA = 0.2;
ISPIN=1
IBRION=-1
NSW=0
2。收敛
ISTART = 0 ; ICHARG=2
PREC= n
ISMEAR = 0; SIGMA = 0.2;
ISPIN=2
VOSKOWN=1
IBRION=-1
NSW=0
3。收敛（使用PBE势）***********因为有人说用此势比较好，故尝试
ISTART = 0 ; ICHARG=2
PREC= n
ISMEAR = 0; SIGMA = 0.2;
ISPIN=2
IBRION=-1
NSW=0
4。收敛（使用PW91势）
ISTART = 0 ; ICHARG=2
PREC= n
ISMEAR = 0; SIGMA = 0.2;
ISPIN=1
IBRION=-1
NSW=0

综上所述，以及我另外做的其他尝试，排除了poscar的影响，以及势函数的影响，发现SIGMA = 0.2是收敛的主要因素。当然我也在以上条件3，4 INCAR中其他条件未变的情况下，仅改变SIGMA = 0.01，发现不收敛了。所以我想
SIGMA = 0.2是收敛的主要因素

2，3，4计算中，4没有计算磁矩，2，3考虑自旋，但计算结果发现
E2=-0.31eV
E3=-0.54ev
2，3区别在于势函数不一样，所以采用不同的势函数计算了体材料Ni的单原子能分别为（我拟和了NI的晶格参数，所以结果都为取平衡晶格参数的值）
PW91 -5.416
PBE -5.454
这样计算的NI的结合能是否为
PW91-E2=-5.106
PBE-E3=-4.914
试验上的数据参考《潜入原子方法理论及其在材料科学中的应用》为4.46ev，所以是不是我的计算还是有问题？还是我的结合能计算方法有问题？