24小时热门版块排行榜    

查看: 2080  |  回复: 15
当前主题已经存档。

ypf1449

银虫 (初入文坛)

K点取1X1X1是没有办法收敛的,我计算的另外一个半金属开始设置为1X1X1到2000步都没有收敛,K点增大以后很快就收敛了!!
芳圆时态
11楼2008-07-02 10:05:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leongxf

银虫 (初入文坛)

计算结果

采取各位的建议,模拟了很多次,终于收敛了,也找到了不收敛的原因,分析一下,请各位多指教。
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,所以是不是我的计算还是有问题?还是我的结合能计算方法有问题?
12楼2008-07-03 11:49:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mermaid1128

木虫 (著名写手)


wuchenwf(金币+1):thank you 1-27 16:05
k点太小 可以先用所有默认值计算
13楼2008-07-03 14:25:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)


wuli8(金币+1,VIP+0):(*^__^*) 2-19 19:47
我以前计算Ti原子能量遇到不收敛情况,后来我是设置IBRION=2解决的。此外,Ni有磁性,采用自旋极化。
14楼2009-02-19 14:35:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangyc

木虫 (正式写手)

应该和K点选取没有关系吧,增加K点数目只能改进描述原子间的相互作用,单原子计算中需要吗?
15楼2009-02-19 16:16:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilyx8050

禁虫 (正式写手)

本帖内容被屏蔽

16楼2009-02-26 16:40:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leongxf 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见