vasp入门到精通[29]自旋轨道耦合SOC

1、比一个6个原子的Mn原胞（），做SOC优化计算

POSCAR 和 INCAR参数如下图，INCAR如下：

SYSTEM=Mn6

ISTART=0

ICHARG=2

GGA=PE

ENCUT=600

ISPIN=2

MAGMOM=3 3 0 -3 3 0 3 3 0 -3 3 0 5 1 0 5 1 0

PREC=Accurate

ISYM=0

EDIFF=1.0E-7

EDIFFG=-0.002

GGA_COMPAT=.FALSE.

VOSKOWN=1

ISMEAR=0

LORBIT=11

NEDOS=2001

EMIN=-8

EMAX=8

LMAXMIX=4

LREAL=.TRUE.

NSIM=4

NPAR=4

NELM=200

IBRION=2

NSW=200

LSORBIT=.TRUE.

SAXIS = 0 0 1

#I_CONSTRAINED_M=1

2、参数注释：

ISPIN=2 (打开自旋极化)

MAGMOM=3 3 0 -3 3 0 3 3 0 -3 3 0 5 1 0 5 1 0（6个原子，每个给定三个方向的磁矩，顺序为r1x,  r1y, r1z, r2x, r2y, rxz... ...所以是有18个数字）

PREC=Accurate（计算精度要高）

EDIFF=1.0E-7

EDIFFG=-0.002

GGA_COMPAT=.FALSE.

VOSKOWN=1

SAXIS = 0 0 1（量子化轴，这个讲起来有点费时，小伙伴设置成这样子就可以）

3、计算结果如下：

运行tail -200 OUTCAR，看一下最后的数据：

4、tail -400 OUTCAR

可以再往上看到每个轨道的矩阵元（小编省略了）

L=1

。。。

L=2

。。。

L=3

。。。

5、从上图可以看出:

最后给出了三个方向的每个原子的磁矩，也就是能看出磁性非共线结果

小编画出这6个原子磁矩的矢量图：

下面的是初始结构，设置磁矩方向

6、再下面的是优化后磁矩的俯视图

7、今天就到这里了，下期见！