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

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的磁化