24小时热门版块排行榜    

查看: 1311  |  回复: 2

小小鑫000

铁虫 (小有名气)

[求助] 关于vasp优化Ni3Fe团簇构型并计算其能量遇到的问题

各位前辈,大家好:
最近在用vasp优化Ni3Fe 团簇(此体系有四个原子)构型时,发现比优化Ni4团簇难收敛一点,我当时觉得正常,后来又进一步静态计算Ni3Fe的能量时,就很难收敛了,走了100多步都不收敛,但是在计算Ni4团簇的能量时,走个20多步就收敛了。我贴出我的文件来:
POSCAR:
CIF file
1.0
       15.0000000000         0.0000000000         0.0000000000
        0.0000000000        15.0000000000         0.0000000000
        0.0000000000         0.0000000000        15.0000000000
   Ni   Fe
    3    1
Direct
     0.494969994         0.464340001         0.395099998
     0.419880003         0.433369994         0.534799993
     0.581040025         0.425199986         0.526149988
     0.507070005         0.571850002         0.518890023
INCAR:
SYSTEM = NiFe
ENCUT = 450
ISTART = 0 ; ICHARG = 2
ISMEAR = 0; SIGMA = 0.05
EDIFF = 1E-4; EDIFFG = -1.0E-2
NSW = 100; IBRION = 2
POTIM = 0.1
ISIF = 2
PREC = Accurate
ISPIN = 2
后来我考虑是磁矩的问题,于是我在INCAR里加入了磁矩并考虑了铁磁性和反铁磁性两种情况:
POSCAR不变,
铁磁性的INCAR为
SYSTEM = Ni
ENCUT = 450
ISTART = 0 ; ICHARG = 2
ISMEAR = 0; SIGMA = 0.05
EDIFF = 1E-4; EDIFFG = -1.0E-2
NSW = 100; IBRION = 2
POTIM = 0.1
ISIF = 2
PREC = Accurate
ISPIN = 2
MAGMON = 3 3 3 3

反铁磁性的INCAR为
SYSTEM = Ni
ENCUT = 450
ISTART = 0 ; ICHARG = 2
ISMEAR = 0; SIGMA = 0.05
EDIFF = 1E-4; EDIFFG = -1.0E-2
NSW = 100; IBRION = 2
POTIM = 0.1
ISIF = 2
PREC = Accurate
ISPIN = 2
MAGMON = 3 3 -3 -3
结果这两种磁性优化出来的结果完全相同,甚至走的每一步都一样,并且和没加时磁矩完全相同,磁矩似乎没起任何作用。

不知道是什么问题,还望各位前辈指教。
回复此楼

» 猜你喜欢

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

好好研究,天天向上!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小地雷2014

铜虫 (小有名气)

顶一下
2楼2015-03-19 15:05:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小地雷2014

铜虫 (小有名气)

楼主为什么反铁磁矩为什么这么设置,还有就是为什么磁矩设置成3 ,有什么依据吗
3楼2015-03-19 15:06:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小小鑫000 的主题更新
信息提示
请填处理意见