HSE06计算态密度和能带

misaraty 收录于学术

2021-01-12 约 720 字预计阅读 2 分钟

前言

简述HSE06方法计算态密度（DOS）和能带（BANDS）的流程。

计算流程

PBE进行Relax计算。
PBE进行SCF计算（ICHARG = 2）。
HSE06进行SCF计算（取2.的WAVECAR，ISTART = 1，ICHARG = 0）。
HSE06进行非自洽计算（DOS、BANDS等）（取3.的WAVECAR、CHGCAR，ISTART = 1，ICHARG = 11）。

示例

以CsPbBr₃单胞、G点为例。

若3.中INCAR，LEPSILON = .T.（计算介电参数），则OUTCAR会报错，HF is not implemented in the linear response routines Please use the finite difference versions.。

若计算缺陷形成能，3.和4.中INCAR，应设置ALGO = Normal（速度快）或者All（速度慢），不推荐设置ALGO = VeryFast、Fast、Damped（能量可能不收敛）。

DOS

运行，

1

sumo-dosplot --config colours.conf --format jpg --dpi 300 --ylabel PDOS --columns 1 --xmin -3 --xmax 3

BANDS

能带计算前，需要将3. SCF_HSE06中的IBZKPT，

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11


Automatically generated mesh
       8
Reciprocal lattice
    0.00000000000000    0.00000000000000    0.00000000000000             1
    0.50000000000000   -0.00000000000000   -0.00000000000000             1
   -0.00000000000000    0.50000000000000   -0.00000000000000             1
    0.50000000000000    0.50000000000000   -0.00000000000000             1
   -0.00000000000000   -0.00000000000000    0.50000000000000             1
    0.50000000000000   -0.00000000000000    0.50000000000000             1
   -0.00000000000000    0.50000000000000    0.50000000000000             1
    0.50000000000000    0.50000000000000    0.50000000000000             1

修改为（仅为示例；实际计算中，高对称点要多些），

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17


Automatically generated mesh
       14
Reciprocal lattice
    0.00000000000000    0.00000000000000    0.00000000000000             1
    0.50000000000000   -0.00000000000000   -0.00000000000000             1
   -0.00000000000000    0.50000000000000   -0.00000000000000             1
    0.50000000000000    0.50000000000000   -0.00000000000000             1
   -0.00000000000000   -0.00000000000000    0.50000000000000             1
    0.50000000000000   -0.00000000000000    0.50000000000000             1
   -0.00000000000000    0.50000000000000    0.50000000000000             1
    0.50000000000000    0.50000000000000    0.50000000000000             1
   0.000   0.000   0.000   0  ! \Gamma
   0.500   0.500   0.000   0  ! M
   0.500   0.500   0.500   0  ! A
   0.000   0.000   0.500   0  ! Z   
   0.000   0.500   0.500   0  ! R
   0.000   0.500   0.000   0  ! X

运行，

1

sumo-bandplot --project Br.p,Pb.s,Pb.p --mode stacked --circle-size 90 --format jpg --dpi 300

或根据vaspkit生成BAND.dat和KPATH.dat，使用Origin、Matplotlib、Matlab等作图。

目录

HSE06计算态密度和能带

计算流程

示例

DOS

BANDS

文档