vasp入门到精通[29]自旋轨道耦合SOC
1、比一个6个骒貉缭塄原子的Mn原胞(),做SOC优化计算POSCAR 和 INCAR参数如下图,INCAR如下:SYSTEM=Mn6ISTART=0ICHARG=2GGA=PEENCUT=600ISPI绿覆冗猩N=2MAGMOM=3 3 0 -3 3 0 3 3 0 -3 3 0 5 1 0 5 1 0PREC=AccurateISYM=0EDIFF=1.0E-7EDIFFG=-0.002GGA_COMPAT=.FALSE.VOSKOWN=1ISMEAR=0LORBIT=11NEDOS=2001EMIN=-8EMAX=8LMAXMIX=4LREAL=.TRUE.NSIM=4NPAR=4NELM=200IBRION=2NSW=200LSORBIT=.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-7EDIFFG=-0.002GGA_COMPAT=.FALSE.VOSKOWN=1SAXIS = 0 0 1(量子化轴,这个讲起来有点费时,小伙伴设置成这样子就可以)
3、计算结果如下:运行tail -200 OUTCAR,看一下最后的数据:
4、tail -400 OUTCAR可以再往上看到每个轨道的矩阵元(小编省略了)L=1。。。L=2。。。L=3。。。
5、从上图可以看出:最后给出了三个方向的每个原子的磁矩,也就是能看出磁性非共线结果小编画出这6个原子磁矩的矢量图:下面的是初始结构,设置磁矩方向
6、再下面的是优化后磁矩的俯视图
7、今天就到这里了,下期见!