VASP计算单原子自旋极化后的原子磁矩

1、首先用MS建模，把氧原子放在8*8*8的立方体中

2、POSCAR文件为O atom in a box 1.0 ! universal scaling parameters8.0 0.0 0.0 ! lattice vector a(1)0.0 8.0 0.0 ! lattice vector a(2)0.0 0.0 8.0 ! lattice vector a(3)1 ! number of atomscart ! positions in cartesian coordinates0 0 0

3、INCAR文件为SYSTEM = O atom in a boxISMEAR = 0 ! 高斯近似ISPIN = 2 ! 自旋极化计算

4、KPOINTS文件为Gamma-point only0Monkhorst Pack1 1 10 0 0

5、O原子是一个有两个未配对电子的开放外层系统。可以获得如下的输出代码

6、自旋起始和自旋结束的本征值分别单独计算。在LSDA（局域自旋密度近似）中，它们只能通过有效的自旋向上或自旋向下本征势能相互作用。在OUTCAR文件中可以看到两个自旋组成，自旋分量1有2个电子对应于2 μ B的磁化