磁性计算求助,设置初始磁矩优化后磁矩为0,想得到反铁磁态的能量如何设置呢?

林家清染

最近在做AFM的计算,在优化结构这一步遇到如题问题,优化的INCAR磁性tag设置如下
IBRION=2
NSW=150
EDIFF=1E-7
EDIFFG=-1E-5
ISIF=2
ISPIN=2
MAGMOM=-4 4 -4 4 4*0 2*0
AMIX = 0.2
BMIX = 0.0001 ! almost zero, but 0 will crash some versions
AMIX_MAG = 0.8
BMIX_MAG = 0.0001 ! almost zero, but 0 will crash some versions
计算以后得到的结果是
magnetization (x)

# of ion       s       p       d       tot
------------------------------------------
    1       -0.000  -0.000  -0.000  -0.000
    2        0.000   0.000   0.000   0.000
    3       -0.000  -0.000  -0.000  -0.000
    4        0.000   0.000   0.000   0.000
    5       -0.000  -0.000  -0.000  -0.000
    6        0.000   0.000   0.000   0.000
    7       -0.000  -0.000  -0.000  -0.000
    8        0.000   0.000   0.000   0.000
    9       -0.000   0.000   0.000  -0.000
   10       -0.000  -0.000  -0.000  -0.000
--------------------------------------------------
tot         -0.000  -0.000  -0.000  -0.000
最后得到的结果是无磁态,想要得到AFM的能量要如何设置参数呢?

林家清染

其实有一点小想法,就是晶格常数.
无磁态的晶格常数我也计算了要大概是3.755左右,这里反铁磁态偷懒了直接用ISIF=3算的结果优化出来的晶格常数也在3.75左右,所以其实结构是自己跑到了NM态附近了吗?(自问无答).还有就是,如果真是这样的话,我还是想得到当前结构的反铁磁态能量要怎么做呢?