求助：磁性体系如何进行SCF收敛

guyingbing

大家好，我是学习理论计算的新手。我现在进行一个磁性体系的结构优化。在体系中一共有三种原子：Co（II）,Mo（IV），O。根据实验和单胞优化结果，设定Co为反铁磁构型，Mo初始磁矩为O，O的初始磁矩0。我首先进行了没有磁性的SCF计算，以便获得磁性计算中使用的波函数和电荷密度。无磁计算顺利收敛，但是进行磁性计算时，我发现体系非常难以收敛。我一共跑了3次scf计算，每一次计算都读取上一次的波函数和电荷密度。计算时发现，第一次scf在425时收敛，第二次scf在185时收敛，但是第三次scf就不收敛。首先，按理来说每一次读取波函数和电荷密度，那么每次收敛应该一次比一次迅速，但是我是一次比一次困难。而且每次收敛的时候，发现能带能快速收敛，但是能量不能收敛。此外，每次scf收敛时，能量是跳跃的。与以往的不一样，以往的是逐渐收敛的。下图1是INCAR，图2到图四分别是第一次，第二次以及第三次scf收敛时的截图。