VASP自旋（NM、AM、AFM）

【背景介绍】

          添加自旋（设置ISPIN=2），就是在哈密顿量上加了一个磁学项A(alpha,beta)，它是一个关于不同自旋的电子数的函数。我们在MAGMON上只是初始猜测，实际计算中SCF会不断调整alpha和beta的数值，从而使得体系的能量逐渐下降最终得到体系真正的电子自旋极化程度，也就是（alpha-beta）的大小。

[举例] 在Si中某个原子加入自旋，然后计算，会发现这个加入的磁矩会vanish，因为Si中本来就不可能有自旋自发极化。

【计算对象】

优化时加不加磁性，不同体系可以不同处理，看优化出来的结果差异大小及跟实验结构对比。

【参数设置】

• ISPIN: 1（默认无自旋）2（自旋）
• MAGMOM: 初始磁距，默认值为Nion*1
            （1）计算FM：MAGMON = 0 0 0 0 0 0 0 0 0 1 0 0 1 0

            备注：此处仅仅是计算铁磁性【FM】，而不考虑非共线，或各向异性或自旋轨道耦合，是不需要SAXIS参数和每个原子设置三个分量的，只需要一个值就好。前面的9个0可以缩写为9*0.初始磁矩还能再大一些。（graphene初始磁矩可以是1，其他 的d轨道应该是4-5）2*0 7*0 1 2*0 1 0。不论你把初始自旋设为1还是4，最终结果都可能一样，例如可能收敛到0.3。而有时候会出现即使开始猜测是AFM，但是体系最终弛豫到了FM的configuration，那就是可能不小心从局部最小跳到了全局最小点了，或者这个体系AFM就是不稳定。
  
  [举例] 我不是所有原子都有磁矩，就是纳米带两头上的原子设了磁矩。对啊，你就两个C原子有磁矩。就算这样，你设置AFM算出来总磁矩也该为0。你要算AFM的话就是把那两个MAGMOM为1的C分别设置成1和-1，磁矩不为0的系统就不是AFM了。一模一样确实比较诡异，按理说FM和AFM的原子位置一样是有可能的，这取决于楼主的优化方式。但是能量应该有几个meV的差别。然后磁矩的话，印象中graphene边界上的C的话如果考虑成FM的话体系的总磁矩应该是在0.0几到0.几的。而AFM体系的总磁矩肯定是0，最多有点误差是0.000几（具体几位记不清了，就是VASP输出的最后一位有数字）。建议楼主可以这么处理AFM：在b方向扩展一下，建一个1*2*1的超胞，然后就得到了4个边界C原子。接着将相邻两个边界C的磁矩分别设置成+1和-1。

          （2）非共线设置：3*NIONS*1.0
          ISPIN = 2
          MAGMOM = 0 0 4 0 0 4 0 0 4 0 0 4 0 0 4
          SAXIS = 0 0 1
          备注：五个原子的磁矩都向z轴的设置。为什么每个原子都设置为4呢？