前言

一、前期准备

二、计算流程

1.静态自洽（scf）

2.能带计算

3.态密度计算

总结

前言

能带结构是目前采用第一性原理计算所得到的常用信息，可用来结合解释金属、半导体和绝缘体的区别。能带用来定性地阐明了晶体中电子运动的普遍特点。价带，或称价电带，通常指绝对零度时，固体材料里电子的最高能量。在导带中，电子的能量范围高于价带，而所有在传导带中的电子均可经由外在的电场加速而形成电流。对与半导体以及绝缘体而言，价带的上方有一个能隙，能隙上方的能带则是传导带，电子进入传导带后才能在固体材料内自由移动，形成电流。对金属而言，则没有能隙介于价带与传导带之间，因此价带是特指半导体与绝缘体的状况。原则上讲，态密度可以作为能带结构的一个可视化结果。很多分析和能带的分析结果可以一一对应，很多术语也和能带分析相通。但是因为它更直观，因此在结果讨论中用得比能带分析更广泛一些。本文主要讲解使用VASP计算BN能带和态密度。

一、前期准备

需要使用VASP、vaspkit、Originlab

vasp的准备，需要INCAR（控制文件），POSCAR（晶格位置文件或者结构文件），POTCAR（赝势文件），KPOINTS（是k点路径），四个输入文件以及一个运行文件。只要有了POSCAR其他的三个文件可以可以通过vaspkit得到。其中赝势文件是提前导入设置好路径。

二、计算流程

1.静态自洽（scf）

新建一个文件夹scf，INCAR可以简单设置为如下

Global Parameters
ISTART =  1            (Read existing wavefunction; if there)
# ISPIN =  2           (Spin polarised DFT)
# ICHARG =  11         (Non-self-consistent: GGA/LDA band structures)
LREAL  = Auto          (Projection operators: automatic)
# ENCUT  =  400        (Cut-off energy for plane wave basis set, in eV)
PREC   =  Normal       (Precision level)
LWAVE  = .TRUE.        (Write WAVECAR or not)
LCHARG = .TRUE.        (Write CHGCAR or not)
ADDGRID= .TRUE.        (Increase grid; helps GGA convergence)
# LVTOT  = .TRUE.      (Write total electrostatic potential into LOCPOT or not)
# LVHAR  = .TRUE.      (Write ionic + Hartree electrostatic potential into LOCPOT or not)
# NELECT =             (No. of electrons: charged cells; be careful)
# LPLANE = .TRUE.      (Real space distribution; supercells)
# NPAR   = 4           (Max is no. nodes; don't set for hybrids)
# NWRITE = 2           (Medium-level output)
# KPAR   = 2           (Divides k-grid into separate groups)
# NGX    = 500         (FFT grid mesh density for nice charge/potential plots)
# NGY    = 500         (FFT grid mesh density for nice charge/potential plots)
# NGZ    = 500         (FFT grid mesh density for nice charge/potential plots)Static Calculation
ISMEAR =  0            (gaussian smearing method)
SIGMA  =  0.05         (please check the width of the smearing)
LORBIT =  11           (PAW radii for projected DOS)
NEDOS  =  2001         (DOSCAR points)
NELM   =  60           (Max electronic SCF steps)
EDIFF  =  1E-08        (SCF energy convergence; in eV)

按照惯例准备四个输入文件INCAR、POSCAR、KPOINT、POTCAR。

都可以通过vaspkit得到

vaspkit输入103，根据POSCAR自动生成POTCAR

vaspkit输入102，根据POSCAR自动生成KPOINT

vaspkit输入1，根据POSCAR自动生成INCAR

准备好文件以后提交任务进行计算

计算完成后会得到CHGCAR文件

2.能带计算

1、新建一个名为band的文件夹，把scf文件夹中的INCAR、POSCAR、KPOINT、POTCAR、CHGCAR复制过来。

2、修改INCAR的内容

ISTART=1

ICHARG=11

LCHARG = .FALSE.

LORBIT=11 # or 10,算投影能带有用

3、使用vaspkit生成高对称点的KPOINTS文件

呼出vaspkit界面后直接输入302自动生成二维结构的高对称点路径，输入303生成三维块体结构的高对称点路径。

会生成一个KPATH.in文件，直接把它cp成KPOINTS文件，然后可以检查修改一下里面的内容：

完成后提交任务进行计算

4、导出作图数据

再次呼出vaspkit，直接输入211，即可生成可以导入OriginLab中绘图的数据BAND.dat。

5、Origin画图

把BAND.dat导入到Origin中做折线图

3.态密度计算

1、态密度计算和计算能带的做法一样，先新建一个目录，命名为dos，然后复制scf目录里面的INCAR、POSCAR、POTCAR、KPOINTS，CHGCAR过来。要修改一下INCAR，或者直接复制下面的INCAR。

Global Parameters
ISTART = 1 (Read existing wavefunction; if there)
ISPIN = 2 (Spin polarised DFT)
ICHARG = 11 (Non-self-consistent: GGA/LDA band structures)
LREAL = Auto (Projection operators: automatic)
# ENCUT = 400 (Cut-off energy for plane wave basis set, in eV)
PREC = Normal (Precision level)
LWAVE = .FALSE. (Write WAVECAR or not)
LCHARG = .TRUE. (Write CHGCAR or not)
ADDGRID= .TRUE. (Increase grid; helps GGA convergence)
LVTOT = .TRUE. (Write total electrostatic potential into LOCPOT or not)
# LVHAR = .TRUE. (Write ionic + Hartree electrostatic potential into LOCPOT or not)
# NELECT = (No. of electrons: charged cells; be careful)
# LPLANE = .TRUE. (Real space distribution; supercells)
# NPAR = 4 (Max is no. nodes; don't set for hybrids)
# NWRITE = 2 (Medium-level output)
# KPAR = 2 (Divides k-grid into separate groups)
# NGX = 500 (FFT grid mesh density for nice charge/potential plots)
# NGY = 500 (FFT grid mesh density for nice charge/potential plots)
# NGZ = 500 (FFT grid mesh density for nice charge/potential plots)

Static Calculation
ISMEAR = 0 (gaussian smearing method)
SIGMA = 0.05 (please check the width of the smearing)
LORBIT = 11 (PAW radii for projected DOS)
NEDOS = 2001 (DOSCAR points)
NELM = 60 (Max electronic SCF steps)
EDIFF = 1E-08 (SCF energy convergence; in eV)
EMIN = -10
EMAX = 10

KPOINTS要设置得比原来做scf的时候要高。这次设成SCF计算的两倍，20-40之间。其他的输入文件不用改。

2、绘制DOS图

打开vaspkit，输入111即可导出TDOS.dat文件，用于绘制DOS图。

总结

本文只是个人刚学习VASP的一点心得体会。在VASP上肯定还有很多不足的地方需要不断改进学习，在此做一点点小小的记录。

使用VASP计算材料能带和态密度相关推荐

VASP计算AgGaS2能带及态密度及光学性质
第一步:构型优化 1．准备四个输入文件 POSCAR INCAR POTCAR KPOINTS POSCAR: 从ms中导入AgGaS2结构,选择CASTEP,file,save,并保存成原包.这样, ...
vasp结构优化英语_VASP态密度、能带、DFT+U实例分析教程
今天主要记录强关联电子体系+U.DOS和能带计算过程中参数的设置,相信大家在读一些计算类文献时会看到文章中计算部分有The GGA + U calculations were performed us ...
使用vasp计算dos（态密度）的一些思考
一:计算前的准备工作 1,得到一个稳定的弛豫过后的POSCR结构 2,对结构进行一次静态计算二:计算态密度的参数准备对于态密度计算INCAR参数的一些思考 ISTAT=1 这里我们选择1,但可以选 ...
CP2K+Gaussian+LAMMPS+ReaxFF计算材料化学系列四大专题。
背景: CP2K是一款较为强大的AIMD计算程序,免费开源,可高效并行.由于CP2K在做DFT时可以速度非常快地计算成千上万个原子的体系,因此在大规模模拟体系中经常被一些学者用到.比如计算原子数达到一 ...
VASP计算弹性常数
记录一下最近计算的不同压力下Cr3C2的弹性常数 1. vasp直接计算 vasp5.1版本后支持设置IBRION=6来直接计算材料的弹性常数,最直接的做法就是在INCAR中设置以下三个参数 IBRI ...
VASP计算新手礼包：基于CentOS 7.6 系统安装VASP 5.4.1
2020-9-24更新:更新VASP5.4.1网盘链接: 链接:https://pan.baidu.com/s/1n9C9CmF0j7iMilo1cPe1cA 提取码:69dj 2019-9-9 更新 ...
VASP计算光学性质
PS:请见文末的打赏选项. VASP计算光学性质的INCAR是根据HSE06计算的自洽INCAR来改的 SCF INCAR: SYSTEM = System ISTART = 1 ICHARG = 1 ...
脑的计算模型能带我们走多远
我们的脑虽然只有1．5公斤重,但是它的结构却极为复杂.这不仅表现为构成它的基本单元--神经细胞的数目达到1俨个,其数量级几乎和我们已知的宇宙中星体的数量级相仿,而且神经细胞和神经细胞之间还存在着极为复 ...
Phonopy-Spectroscopy计算材料红外和Raman光谱
经验分享 Phonopy- spectroscopy是一个为Phonopy代码添加模拟振动光谱能力的项目. 它可以从Phonopy或VASP计算红外(IR)强度,和计算Raman活度张量和远离共振近似 ...

使用VASP计算材料能带和态密度

前言

一、前期准备

二、计算流程

1.静态自洽（scf）

2.能带计算

3.态密度计算

总结

使用VASP计算材料能带和态密度相关推荐

最新文章

热门文章