1、首先用MS建模，把氧原子放在8*8*8的立方体中
2、POSCAR文件为
O atom in a box 1.0 ! universal scaling parameters
8.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 atoms
cart ! positions in cartesian coordinates
0 0 0
3、INCAR文件为
SYSTEM = O atom in a box
ISMEAR = 0 ! 高斯近似
ISPIN = 2 ! 自旋极化计算
4、KPOINTS文件为
Gamma-point only
0
Monkhorst Pack
1 1 1
0 0 0
5、O原子是一个有两个未配对电子的开放外层系统。
可以获得如下的输出代码
6、自旋起始和自旋结束的本征值分别单独计算。在LSDA（局域自旋密度近似）中，它们只能通过有效的自旋向上或自旋向下本征势能相互作用。
在OUTCAR文件中可以看到两个自旋组成，自旋分量1有2个电子对应于2 μ B的磁化
