纯干货：解读VASP手册（三）

今天的教程继续解读VASP手册，我们直接进入正题。 这里面提到了一个重要的参数，EDIFF。这个参数是电子自洽迭代的循环中，总能收敛的标准。官网手册给出的建议是设置EDIFF=10-4即可满足大多数情况的要求，因此没有特殊要求一般不需要设置更小的数值。

这里说到了IBRION=0的情况，IBRION=0表示标准的分子动力学模拟。采用Verlet算法来积分原子的牛顿运动方程。通过POTIM来控制时间步长(单位是fs)。SMASS控制系综的设置。 IBRION=3是采用最速下降算法来优化原子的位置。在CI-NEB方法计算过渡态时需要设置IBRION=3。对于IBRION=3，POTIM的默认值为0.5，手册建议必须通过POTIM参数提供合理的时间步长。 这里是对IWAVPR如何设置的建议。IWAVPR确定如何将轨道和/或电荷密度从一种离子构型外推至另一种构型。IWAVPR可以设置为0 | 1 | 2 | 3 | 10 | 11 | 12 | 13。手册建议不要设置为13。因为IWAVPR=13的效果与IWAVPR=1和2之间的混合的效果一样，但是，通常比IWAVPR=2更糟。 这里面是针对ISMEAR的建议。ISMEAR决定了如何确定每个波函数的占有数。ISMEAR可以取-5 | -4 | -3 | -2 | 0 | N，当ISMEAR=N时，表示Methfessel-Paxton N阶smearing方法。这里给出的建议是针对于ISMEAR=N的情况。 这部分对建议ISMEAR的取值给出了建议。对于半导体和绝缘体不要设置ISMEAR>0，绝缘体一般设置ISMEAR=0或ISMEAR=-5。一般情况下，采用ISMEAR=0，并选取一个合适的SIGMA值都是可以得到合理的结果。 这里面说到了计算偶极时的设置问题。计算偶极需要设置“DIPOL = center of cell (in direct coordinates)”，计算偶极需要定义晶胞中心，这里给出了关于晶胞中心如何确定的建议，有需要计算偶极的朋友一定要注意这里的建议。 这里面有一条关于IALGO和LDIAG设置的提示，IALGO用来具体确定电子优化采用何种算法，在4.5以上的版本中其默认值为38。LDIAG是用来确定是否进行子空间对角化，默认值是.TRUE.。 以上就是本次教程解读的VASP手册，后面的教程将继续给大家解读VASP手册。 参考资料： https://www.vasp.at/wiki/index.php/The_VASP_Manual