赞
踩
我们通过前面的教程带大家了解了VASP理论计算的基础知识,包括准备输入文件,解读输出文件,计算电子结构,分析电荷等方面,建议大家根据教程多加练习。最近笔者在做计算的过程中,遇到了一些棘手的问题,困扰了很久,最后翻阅VASP手册解决了问题。笔者又一次意识到了VASP手册的重要性,因此笔者决定接下来几期教程给大家一起解读一下VASP手册,和大家一起夯实基础。
在VASP手册里面,经常会有“Mind:”、“IMPORTANT:”以及“IMPORTANT NOTE:”开头的段落,这些就是大家常见的疑难点以及易错的地方,笔者主要是从这些点来和大家一起来解读一下VASP手册。在这里按照从前向后的顺序进行解读。
这里面提到了关于POTCAR的问题,需要注意的是POTCAR文件里面的赝势顺序需要和POSCAR里面的原子顺序一致,这个是前面强调过的。同时在这里强调一下,POTCAR里面有默认的截断能数值(ENMAX 和 ENMIN),对于多原子体系,如果INCAR中不设置截断能,则默认取所有截断能默认值的最大值。
关于自动生成k-网格的一个说明,生成的k-网格与倒数单元必须属于同一类Bravais格子,否则VASP将会显示计算错误。这里面还有一点需要注意,即如果网格的原点不是Γ,则k网格的对称性可能导致无法将网格划分为四面体。因此这里给出的建议是如果使用了四面体方法,最好仅使用原点在Γ点的网格。
“Selective dynamics”这项参数不是必须的。在进行结构弛豫时,对于需要固定部分原子进行弛豫的体系,需要在POSCAR里面加上“Selective dynamics”。如果没有添加“Selective dynamics”,VASP默认将所有的原子都驰豫。如果添加了“Selective dynamics”,但是原子后面没有加TTT或者FFF,VASP也会默认所有原子进行驰豫。
这里面提到了CHGCAR文件,dynamic simulation (IBRION=0)得到的CHGCAR文件的电荷率不是CONTCAR文件中原子位置的自洽电荷率,因此不能以该CHGCAR文件为基础进行能带结构计算。对于静态计算和弛豫(IBRION=-1,1,2),写入的电荷密度是最后一步的自洽电荷密度,可用于精确的能带结构计算。
同时,dynamic simulation (IBRION=0)的WAVECAR文件也不能直接用于其他计算,需要以CONTCAR中的原子坐标为结构文件并在INCAR中设置ISTART=1; NSW=0; ICHARG=1进行静态计算即可。
dynamic simulation (IBRION=0)的EIGENVAL文件也需要执行上述相同操作才可以用于其他计算。
以上是本次教程解读的手册内容,有些内容过于基础,笔者认为没有解读的必要就没有写,对于前面教程已经讲过的点,在这篇教程也没有提及,后面教程讲到相关的内容时可能会讲。对于一些英文单词没有翻译,主要是为了避免歧义。
参考内容:
https://www.vasp.at/wiki/index.php/The_VASP_Manual
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。