算例01-以O原子为切入理解磁矩信息

本次案例理解自旋极化对轨道占据和能量的影响

氧原子轨道

因此所以氧原子含有两个单电子,具有磁性,且磁矩为2 muB

先看没开磁矩的情况

INCAR设置

SYSTEM = 0
ISMEAR = 0
SIGMA = 0.01
ISPIN = 1 #不考虑自旋极化

查看能带情况

[zsli@c01n01 O]$ grep -A10 "band energies" OUTCAR 
  band No.  band energies     occupation 
      1     -23.9603      2.00000
      2      -9.0203      1.33333
      3      -9.0203      1.33333
      4      -9.0203      1.33333
      5      -0.2407      0.00000
      6       0.2977      0.00000
      7       0.5630      0.00000
      8       0.5630      0.00000
      9       0.5630      0.00000
     10       0.6087      0.00000

可以看到能带1上填充了两个电子，与原子结构模型吻合；
但是在能带2-3上电子均匀分布了，一共4个电子，总数是对的但是分布和模型显然不同
再看能量

因此应该打开自旋极化再试试

打开自旋极化情况

INCAR设置情况

SYSTEM = 0
ISMEAR = 0
SIGMA = 0.01
ISPIN = 2

能带分布情况

[zsli@c01n01 Ospin]$ grep -A10 "band energies" OUTCAR 
  band No.  band energies     occupation 
      1     -25.2168      1.00000
      2     -10.2108      1.00000
      3     -10.2108      1.00000
      4     -10.2108      1.00000
      5      -0.2463      0.00000
      6       0.2707      0.00000
      7       0.5824      0.00000
      8       0.5824      0.00000
      9       0.5824      0.00000
     10       0.5962      0.00000
--
  band No.  band energies     occupation 
      1     -21.9684      1.00000
      2      -7.1826      0.33333
      3      -7.1826      0.33333
      4      -7.1826      0.33333
      5      -0.1495      0.00000
      6       0.3380      0.00000
      7       0.6183      0.00000
      8       0.6183      0.00000
      9       0.6183      0.00000
     10       0.6941      0.00000

这里的轨道任然是均匀分布的不符合实际，这种情况是由于晶胞过于对称造成的，查看POSCAR

Transfer by qvasp
    1.0
    15.000000    0.000000    0.000000
    0.000000    15.000000    0.000000
    0.000000    0.000000    15.000000
  O
  1
  Direct
  0.500000  0.500000  0.500000

三个轴都是公整的15.000，因此只要打破这个高度对称即可 （正常计算的时候不会有这么高对称度因此往往不会遇到这种情况）

对称破缺以后

修改POSCAR

Transfer by qvasp
    1.0
    16.000000    0.000000    0.000000    
    0.000000    15.000000    0.000000    
    0.000000    0.000000    14.000000    
  O   
  1  
  Direct
  0.500000  0.500000  0.500000

对应的能带情况

[zsli@c01n01 Ospin]$ grep -A10 "band energies" OUTCAR 
  band No.  band energies     occupation 
      1     -25.2191      1.00000
      2     -10.8963      1.00000
      3     -10.8942      1.00000
      4      -8.7966      1.00000
      5      -0.2850      0.00000
      6       0.2053      0.00000
      7       0.4635      0.00000
      8       0.4664      0.00000
      9       0.5311      0.00000
     10       0.5986      0.00000
--
  band No.  band energies     occupation 
      1     -21.5541      1.00000
      2      -7.5418      1.00000
      3      -6.4334      0.00000
      4      -6.4330      0.00000
      5      -0.1813      0.00000
      6       0.3336      0.00000
      7       0.4939      0.00000
      8       0.4953      0.00000
      9       0.5581      0.00000
     10       0.5952      0.00000

可以明显看出自旋的极化与理论模型吻合了，
顺便再看看能量

-1.89明显比之前的-0.013小了不少因此更加合理了