ASE的Vasp函数（VASP计算杂化能带详细步骤教程）

ASE的Vasp函数（VASP计算杂化能带详细步骤教程）Procedure 1: Standard procedure (suitable for DFT calculations)7. 用 vaspkit-252提取杂化密度泛函能带数据。4. 调用 vaspkit-251命令产生用于能带计算的 KPOINTS文件， KPOINTS文件包含两部分为：第一部分为不可约布里渊区权重不为零的K点，用于自洽计算得到正确的费米能；第二部分为权重为零的沿着特定能带的K点，用于能带计算;5. 先用 PBE计算一次并保存波函数。切记计算结束后，用 vaspkit-252命令检查一下能带是否合理。概念一定要清楚， PBE能带也可以采用这种方法得到，和常规分两步计算得到能带是等价的。6. 修改 INCAR关于 HSE计算参数，然后继续执行一次 VASP计算；

本文前面部分为VASPKIT开发者Vei Wang老师所写，后面部分为笔者自己所加。

一、VASP结合VASPKIT软件HSE06计算步骤:

1. 准备 POSCAR，调用 vaspkit-303(体相材料)或 -302(二维材料)得到 KPATH.in和 PRIMCELL.in文件；对于二维体系，需要检查 PRIMCELL.in文件的真空层是否沿 z方向，如果没有，可调用 vaspkit-923或者 vaspkit-407强制真空层沿 z方向。另外，如果需要改变识别结构对称性精度(symprec默认值为1E-5)，可通过 vaspkit-task302-symprec1E-6实现。

2. cp PRIMCELL.vasp POSCAR后，用 VASP-PBE优化结构。 KPATH.in能带只针对于 PRIMITIVE CELL，缺少这一步，你可能得到错误的结果。如果有必要，比较 KPATH.in文件中的能带路径是否与在线能带路径产生工具 SeeK-Path产生的一致，包括比较 PRIMCELL.vasp和 HIGH_SYMMETRY_POINTS文件。需要指出的是 SeeK-Path只用于体相结构能带路径的产生。

3. cp CONTCAR POSCAR，用 VASP-HSE继续优化结构；严格来说 HSE优化的结构与 PBE方法得到的有差别，只是不同的体系差别程度有所不同(可选)

4. 调用 vaspkit-251命令产生用于能带计算的 KPOINTS文件， KPOINTS文件包含两部分为：第一部分为不可约布里渊区权重不为零的K点，用于自洽计算得到正确的费米能；第二部分为权重为零的沿着特定能带的K点，用于能带计算;

5. 先用 PBE计算一次并保存波函数。切记计算结束后，用 vaspkit-252命令检查一下能带是否合理。概念一定要清楚， PBE能带也可以采用这种方法得到，和常规分两步计算得到能带是等价的。

6. 修改 INCAR关于 HSE计算参数，然后继续执行一次 VASP计算；

7. 用 vaspkit-252提取杂化密度泛函能带数据。

二、下面附上VASP官网HSE06计算教程：

Procedure 1: Standard procedure (suitable for DFT calculations)

Only possible within DFT. Described in Fcc Si bandstructure example:

Standard self-consistent (SC) run

POSCAR

INCAR (see INCAR.dft)

KPOINTS (see KPOINTS.6)

Non-SC calculation (ICHARG=11)

Use preconverged CHGCAR file and a suitable KPOINTS file

INCAR

KPOINTS (see KPOINTSPBEbands)

Procedure 2: 0-weight (Fake) SC procedure (works DFT & hybrid functionals)

This procedure can be applied to compute bandstructure at hybrid functionals and DFT level (see the HSE_bandstructure.sh script).

Standard DFT run

Just as before

INCAR (see INCAR.dft)

KPOINTS (see KPOINST.6)

Hybrid calculation using a suitably modified KPOINTS file

INCAR (see INCAR.hse)

KPOINTS (see KPOINTSHSEbands.6 and README.txt)

Please note that step two requires a WAVECAR obtained from a standard DFT run (not an HSE calculation) otherwise the resulting conduction bands often have a zig-zag structure.

三、网上其它经典HSE06计算教程：

1) http://muchong.com/bbs/viewthread.php?tid=4232787&fpage=1

2) http://muchong.com/bbs/viewthread.php?tid=6105659

3) http://blog.sciencenet.cn/blog-567091-732988.html