24小时热门版块排行榜

返回列表

当前主题已经存档。

ypf1449

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 323.5
帖子: 16
在线:
虫号: 547775
注册: 2008-04-17
性别: MM
专业: 凝聚态物理

K点取1X1X1是没有办法收敛的，我计算的另外一个半金属开始设置为1X1X1到2000步都没有收敛，K点增大以后很快就收敛了！！

赞一下

回复此楼

芳圆时态

11楼2008-07-02 10:05:43

已阅回复此楼关注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，所以是不是我的计算还是有问题？还是我的结合能计算方法有问题？

赞一下(1人)

回复此楼

12楼2008-07-03 11:49:37

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