本章教程将继续带领大家解读VASP手册。大家可以跟着前面的教程多加练习,循序渐进,夯实基础。解读VASP手册系列教程之后,笔者会做一些实操性比较强的教程,供大家练习。
这里面提到了DOSCAR,我们知道,DOSCAR文件包含的是态密度和积分态密度,因此对于弛豫过程一般不需要生成DOSCAR文件,如果需要得到精确的DOSCAR文件,需要以CONTCAR中的原子坐标为结构文件并在INCAR中设置ISTART=1;NSW=0进行静态计算即可。这是一条关于PROCAR文件的提示,PROCAR文件包含了每个能带的投射波函数信息,必须在INCAR文件中指定RWIGS才能获取该文件。在并行版本中,如果NPAR≠1,则不会评估每个能带的投射波函数信息。
关于LOCPOT文件,该文件包含了总局域势,这里面提到一点,需要确认main.F文件中LEXCHG=-1是否被注释掉,如果LEXCHG=-1被注释掉,则需要在INCAR中设置PREC=High来避免“wrap around errors”。
NGX,NGY,NGZ控制FFT网格中沿三个点阵矢量方向的网格点数。NGXF,NGYF,NGZF控制实空间函数的网格密度,这段文字是对NGXF,NGYF,NGZF参数设置的说明。
这里面讲到了INCAR中的ISTART参数。如果WAVECAR文件丢失或者WAVECAR文件包含了不准确的信息,则需要设置ISTART=0。如果超胞的尺寸、形状发生改变,截断能也和前面的计算不同,并且希望根据新的设置重新定义一组平面波,此种情况设置ISTART=1。如果希望以上一次任务中相同基组重新启动计算任务,通常使用ISTART=2。这里面有个提示,如果ISTART=3,WAVECAR或TMPCAR文件不存在或任何输入数据不一致都将立即停止任务。
这里提到了ISTART的值对ICHARG的影响。ICHARG控制如何产生初始电荷密度,ICHARG=0表示从初始波函数计算电荷密度,ICHARG=1表示从CHGCAR文件读入电荷密度,ICHARG=2则表示用原子电荷密度组装,ICHARG=11表示读入自洽的CHGCAR,并进行能带计算或态密度的非自洽计算,ICHARG=12则表示非自洽的原子密度计算。如果ISTART=0 则ICHARG=2。
以上就是本章教程需要讲解的内容,大家遇到县官问题可以直接回到教程查找解决方法,后续教程会继续讲解VASP手册。
参考资料:https://www.vasp.at/wiki/index.php/The_VASP_Manual