View Categories

常用INCAR参数

全局通用参数

SYSTEM

这个参数是VASP提供给用户用于描述自己此次计算目的、体系、方法等等,用户可以随意填写,也可以不写。
当自洽计算填写此参数(比如SYSTEM=XXX)并输出CHGCAR后,后续非自洽计算如果读取这个CHGCAR,计算刚开始会显示CHGCAR read from XXX

ISTART

ISTART = 0,计算不会读取WAVECAR文件,此时VASP会根据INIWAV的值来确定初始化轨道的方法
ISTART = 1,轨道由WAVECAR文件读取,并重新读取INCAR中的截断能和POSCAR中的晶格矩阵进行计算(即使这些值与WAVECAR中的旧值不同),如果检测到不存在WAVECAR文件,则VASP按ISTART=0执行计算。如果此次计算相比于上次计算晶胞的大小/形状或截断能发生变化,并且希望读取WAVECAR开始计算,则建议将ISTART设置为1
ISTART = 2,轨道由WAVECAR文件读取,即使INCAR中的截断能或者POSCAR中的晶格矩阵与WAVECAR文件上存储的值不同,平面波也不会更改
ISTART = 3,与ISTART=2相同,但是需要额外读取TMPCAR文件,此文件内容是含时间步长的位置和轨道信息,一般不推荐设置!另外如果不存在WAVECAR或TMPCAR文件,计算会立即终止

ICHARG = 0,从初始波函数(WAVECAR)计算电荷密度。若WAVECAR文件无效或缺失,VASP会自动修正为ICHARG = 2。适用于续算场景,但若波函数文件不匹配(如截断能或晶胞尺寸改变),可能导致收敛问题!一般结构优化不推荐输出WAVECAR,因此ICHARG = 0不常用
ICHARG = 1,从CHGCAR文件读取电荷密度,并通过原子电荷密度的线性组合方法从旧位置(CHGCAR中的)外推至新位置。适用于结构微调后的续算(如弛豫或分子动力学),可加速收敛,但是一般结构优化不推荐输出CHGCAR,因此ICHARG = 1不常用
ICHARG = 2,使用原子电荷密度的叠加作为初始电荷密度,当ISTART=0(或者没有WAVECAR文件时),VASP会默认采用此模式。ICHARG = 2适用于全新计算,常用!
ICHARG = 4,从POT文件读取势场(需LVTOT=.TRUE.生成该文件),常用于优化有效势(OEP)方法
ICHARG = 5,外部电荷密度更新模式,允许在SCF过程中通过GAMMA文件修正KS占据数,用于DFT+DMFT等自洽计算,需配合vasp.lock文件控制流程
ICHARG = 10 | 11 | 12 ,通过给ICHARG值加10(也就是0、1、2变为10、11、12),在保持ICHARG=0、1、2的功能前提下可使电荷密度在整个电子最小化过程中保持恒定,用于非自洽计算。适用场景:
ICHARG=10,常用于HSE计算能带结构;
ICHARG=11,常用于基于CHGCAR给定电荷密度计算能带结构和态密度;
ICHARG=12,原子电荷密度叠加的非自洽计算(遵循Harris-Foulkes泛函原理),可正确计算应力和力,支持基于该泛函的从头算分子动力学。
当ICHARG=11/12时,强烈建议将LMAXMIX设为赝势中最大l量子数的两倍(s/p元素设为2,d元素设为4,f元素设为6)

LWAVE

LWAVE = .TRUE.,写入 WAVECAR;LWAVE = .FALSE.,不写入 WAVECAR
WAVECAR文件常常占用大量存储空间,通常来说结构优化、分子动力学等等不需要后续读取WAVECAR计算的任务,都不需要输出WAVECAR
因此只有后续明确会用到WAVECAR的场景,再写LWAVE = .TRUE.,其余情况一律写LWAVE = .FALSE.

LCHARG

LCHARG = .TRUE.,写入 CHGCAR和CHG 文件;LCHARG = .FALSE.,不写入 CHGCAR和CHG 文件
CHGCAR和CHG 文件常常占用大量存储空间,通常来说结构优化、分子动力学等等不需要后续读取CHGCAR和CHG 文件计算的任务,都不需要输出CHGCAR和CHG 文件
因此只有后续明确会用到CHGCAR和CHG 文件的场景,再写LCHARG = .TRUE.,其余情况一律写LCHARG = .FALSE.

ISMEAR

ISMEAR > 0​​:采用Methfessel-Paxton方法,阶数为ISMEAR值,展宽宽度由SIGMA决定。注意:这种方法可能导致绝缘体的部分占据数不物理
ISMEAR = 0​​:采用高斯展宽方法,展宽宽度由SIGMA决定
ISMEAR = -1​​:采用费米展宽方法,展宽宽度由SIGMA决定
ISMEAR = -2​​:从WAVECAR中读取部分占据数并在整个计算过程中保持不变。或者可以通过INCAR文件中的FERWE(和ISPIN=2时的FERDO)标签设置占据数
ISMEAR = -3​​:对INCAR文件中SMEARINGS参数提供的值进行循环计算
​ISMEAR = -4​​:无展宽的四面体方法
​ISMEAR = -5​​:带Blöchl修正的无展宽四面体方法
ISMEAR = -14​​:带费米-狄拉克展宽(SIGMA)的四面体方法
ISMEAR = -15​​:带Blöchl修正和费米-狄拉克展宽(SIGMA)的四面体方法

不同的ISMEAR值适用于不同的体系
CHGCAR和CHG 文件常常占用大量存储空间,通常来说结构优化、分子动力学等等不需要后续读取CHGCAR和CHG 文件计算的任务,都不需要输出CHGCAR和CHG 文件
因此只有后续明确会用到CHGCAR和CHG 文件的场景,再写LCHARG = .TRUE.,其余情况一律写LCHARG = .FALSE.