【求助】为什么静态计算和原先的磁矩差异如此之大！
RT，在结构优化后我把NSW等注释掉，进行静态计算，
但是，为什么得出的磁矩根刚才计算时的出入非常大！
结构优化时我的磁矩是4.1，静态计算后变成了6.7,这是怎么回事？
恳请各位帮忙看一下，谢谢了！
我作静态计算的的INCAR文件：
SYSTEM = layer
ISTART = 0
ICHARG = 2
ENCUT = 520
ISPIN = 2
#MAGMOM = 8*0.0 8*1.0 24*0.0
EDIFF = 1E-4
EDFIIG = -1E-2
ISMEAR = 2
NPAR = 4
#NELM=200
#NELMDL=8
#WEIMIN=0
AMIX = 0.2
BMIX = 0.00001
AMIX_MAG = 0.8
BMIX_MAG = 0.00001
LREAL = .TRUE.
Dynamic:
#IBRION = 2
#NSW = 500
#POTIM = 0.2
#LORBIT = 11
#LDAU = .TRUE.
#LDAUTYPE = 2
#LDAUL = 2 -1 -1
#LDAUU = 3.00 0.00 0.00
#LDAUJ = 0.87 0.00 0.00
#LDAUPRINT = 2
LMAXMIX = 4 ! Important: mix paw occupancies up to L=4
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
结构优化时的给出的体系性质是不可信的，应该看结构优化后进行自洽计算的结果
我在ESPRESSO计算graphene ribbon时也遇到类似的情况。
做结构优化时，出来的是个顺磁的
做静态计算出来是反铁磁的：（
关键是文献报道的是饭铁磁的：/
你在做静态计算的时候没有输入MAGMOM值，所以vasp默认按照铁磁序优化，有可能和你原来的磁结构不相同，而恰好体系在铁磁结构也有个稳定态，就得到了不同的结果。vasp说明里面提到可以在INCAR里设置ICHARG=1让程序从前面的CHGCAR文件读取电荷密度，同时也读取磁矩的信息，但是几乎很难以实用，因为改变了k点或者NBANDS等等参数后，可能造成读取不到合适的电荷密度，所以我一直都在静态计算的时候也设置MAGMOM值，这样和优化时得到的磁矩基本上相�
，
那么如果这个静态优化的能量比结构优化时给出的能量低0,可以认为这个磁矩是体系最稳态时的磁矩么？
做静态计算时没有输入MAGMOM值
默认的MAGMON是所有原子均为1的铁磁性排列吧。