05_vasp态密度计算 - Yiwei666/13_vasp GitHub Wiki

计算流程

  1. 自洽计算获取CHGCAR文件
  2. vaspkit获取KPOINTS文件
  3. 能带非自洽计算（非必要，可选）
  4. 态密度非自洽计算
  5. vaspkit后处理态密度数据
  6. origin绘图

自洽计算获取CHGCAR文件

参考能带计算

vaspkit获取KPOINTS文件

参考能带计算

态密度非自洽计算

态密度计算需要如下文件：INCAR、POSCAR、POTCAR、KPOINTS、CHGCAR，上述文件中除INCAR文件外，其余文件与能带计算均相同。

  • INCAR输入文件
ALGO = Normal
EDIFF = 1e-05
ENCUT = 520
IBRION = -1
ISIF = 3
ISMEAR = 0
ISPIN = 2

ICHARG = 11

KPAR = 8
LAECHG = True
LASPH = True

# LCHARG = True
LCHARG = .FALSE.

LORBIT = 11
LREAL = Auto
LVHAR = True
LWAVE = False
MAGMOM = 4*0.003 4*-0.0
NELM = 100
NPAR = 4
NSW = 0
PREC = Accurate
SIGMA = 0.05
SYSTEM = Mp-8044

EMIN = -10     #此为DOS图的能量范围，根据能带的能量范围来决定min和max是多少。
EMAX = 10
NEDOS  = 801   #这个应该是DOS的采点个数，弄高一点无所谓。

注意：上述态密度计算的INCAR文件相比于能带计算的INCAR文件，主要新增了如下部分

EMIN = -10     #此为DOS图的能量范围，根据能带的能量范围来决定min和max是多少。
EMAX = 10
NEDOS  = 801   #这个应该是DOS的采点个数，弄高一点无所谓。
  • 态密度非自洽计算输出文件
.
├── AECCAR0
├── AECCAR1
├── AECCAR2
├── CHG
├── CHGCAR
├── CONTCAR
├── DOSCAR
├── EIGENVAL
├── INCAR
├── KPOINTS
├── LOCPOT
├── OSZICAR
├── OUTCAR
├── PCDAT
├── POSCAR
├── POTCAR
├── PROCAR
├── REPORT
├── slurm-2678854.out
├── sub.sh
├── vasprun.xml
├── WAVECAR
└── XDATCAR

vaspkit后处理态密度数据

  • 操作流程
vaspkit --> 11 --> 113
  • vaspkit读取文件和输出文件

读取文件：INCAR，DOSCAR，PROCAR 等文件

 -->> (01) Reading Input Parameters From INCAR File...
 +---------------------------------------------------------------+
 |       >>> The Fermi Energy will be set to zero eV <<<         |
 +---------------------------------------------------------------+
 -->> (02) Reading Fermi-Energy from DOSCAR File...
 -->> (03) Reading DOS Data From DOSCAR File...
 -->> (04) Reading Structure from POSCAR File...

输出文件：各种元素的 PDOS_V_UP/DW.dat 和 IPDOS_V_UP/DW.dat

 -->> (05) Written PDOS_V_UP/DW.dat Files!
 -->> (06) Written IPDOS_V_UP/DW.dat Files!
 -->> (07) Written PDOS_P_UP/DW.dat Files!
 -->> (08) Written IPDOS_P_UP/DW.dat Files!
 -->> (09) Written PDOS_C_UP/DW.dat Files!
 -->> (10) Written IPDOS_C_UP/DW.dat Files!
  • vaspkit导出的态密度相关数据文件
.
├── CHGCAR
├── DOSCAR
├── FERMI_ENERGY
├── INCAR
├── IPDOS_C_DW.dat
├── IPDOS_C_UP.dat
├── IPDOS_S_DW.dat
├── IPDOS_S_UP.dat
├── IPDOS_V_DW.dat
├── IPDOS_V_UP.dat
├── KPOINTS
├── OUTCAR
├── PDOS_C_DW.dat
├── PDOS_C_UP.dat
├── PDOS_S_DW.dat
├── PDOS_S_UP.dat
├── PDOS_V_DW.dat
├── PDOS_V_UP.dat
├── POSCAR
├── POTCAR
└── PROCAR

origin绘图

将上述vaspkit导出的PDOS文件用于绘图，输出文件包含了每种元素的上/下自旋的态密度，对于闭壳层，绘制一种即可

├── PDOS_C_DW.dat
├── PDOS_C_UP.dat
├── PDOS_S_DW.dat
├── PDOS_S_UP.dat
├── PDOS_V_DW.dat
├── PDOS_V_UP.dat

参考资料