VASP计算前的验证

VASP计算前的验证

一、检验赝势的好坏：


（一）方法：对单个原子进行计算；

（二）要求：
1、对称性和自旋极化均采用默认值；
2、ENCUT要足够大；
3、原胞的大小要足够大，一般设置为15 Å足矣，对某些元素还可以取得更小一些。

（三）以计算单个Fe原子为例：

1、INCAR文件：

SYSTEM = Fe atom
ENCUT = 450.00 eV
NELMDL = 5 ! make five delays till charge mixing，详细意义见注释一
ISMEAR = 0
SIGMA=0.1

2、POSCAR文件：

atom
15.00
1.00 0.00 0.00
0.00 1.00 0.00
0.00 0.00 1.00
1
Direct
0 0 0

3、KPOINTS文件：

Automatic
0
Gamma
1 1 1
0 0 0

4、POTCAR文件：（略）

注释一：关键词“NELMDL”:
A）此关键词的用途：指定计算开始时电子非自洽迭代的步数（即NELMDL gives the number of non-selfconsistent steps at the beginning），目的是make calculations faster。“非自洽”指的是保持charge density不变，由于Charge density is used to set up the Hamiltonian, 所以“非自洽”也指保持初始的哈密顿量不变。
B）默认值（default value）:
NELMDL = -5 (当ISTART=0, INIWAV=1, and IALGO=8时)
NELMDL = -12 (当 ISTART=0, INIWAV=1, and IALGO=48时)
NELMDL = 0 (其他情况下)

NELMDL might be positive or negative.
A positive number means that a delay is applied after each ionic movement -- in general not a convenient option. （在每次核运动之后）
A negative value results in a delay only for the start-configuration. （只在第一步核运动之前）

C）关键词“NELMDL”为什么可以减少计算所需的时间？
Charge density is used to set up the Hamiltonian, then the wavefunctions are optimized iteratively so that they get closer to the exact wavefunctions of this Hamiltonian. From the optimized wavefunctions a new charge density is calculated, which is then mixed with the old input-charge density. A brief flowchart is given below.（参自Manual P105页）

一般情况下，the initial guessed wavefunctions是比较离谱的，在前NELMDL次非自洽迭代过程中保持charge density不变、保持初始的哈密顿量不变，只对wavefunctions进行优化，在得到一个与the exact wavefunctions of initial Hamiltonian较为接近的wavefunctions后，再开始同时优化charge density。这样一来，计算时间要比一开始就同时优化charge density 和wavefunctions短得多。

注释二：为什么这里只需要一个k点？
For atoms and molecules, the Bloch theorem does not apply, hence there is no need to use more than one single k-point. When more k-points are used, only the interaction between the atoms (which should be zero) is described more accurately.

（三）计算任务执行方法：
输入：vasp

（四）赝势好的判断标准：计算得到的OUTCAR文件中的“energy without entropy”能量值在 0.001～0.01 eV之间。