2.3.4 VASP计算输出总结

最后，我们对输出文件简单做一个总结：

（1） 确认计算顺利结束：result文件最后没有warning，OUTCAR最后输出了“aborting loop because EDIFF is reached”以及能量信息和末尾的计算时间。

（2） 查看电子步收敛情况，看result或者OSZICAR文件；

（3） 查看体系的压强和每个原子受力情况，看OUTCAR文件

（4） 查看体系的总能量，OSZICAR中的E0

目前为止，你已经走完了一次完整的、从结构到输入文件，再到输出文件的计算过程。这个计算是关于金刚石单质Si的静态自洽计算。静态的含义就是离子不进行移动，即INCAR中NSW=0。现在我们需要回过头来理解，自洽计算是什么，以及我们能从一次自洽计算中得到什么。

这里的自洽含义体现在计算过程中会不断更新电荷密度，参见VASPwiki介绍：

https://www.vasp.at/wiki/index.php/Self-consistency_cycle

与自洽相对的是非自洽计算，对于非自洽计算在计算过程中是不更新电荷密度的：要么只停留在初猜电荷密度（ICHARG=12），要么就是使用已有的电荷密度文件（ICHARG=11）。非自洽计算一般在计算能带和更密k点的dos时会用到。

那么一次SCF可以得到什么信息？总的来说，自洽计算会输出体系当前原子构型下的能量（OSZICAR/OUTCAR）、原子此时受力以及体系压强（OUTCAR）、电子态密度信息和费米能级（DOSCAR）、波函数文件（WAVECAR）和电荷密度文件（CHGCAR）……这些是在后续其他计算需要用到的。比如有了能量就可以比较两个结构之间的相对（能量）稳定性；能量、力、位力等信息，可用于后续拟合机器学习势函数。有了电子态密度信息可以直接判断体系是否为金属/半导体等。

某些情况下我们需要计算孤立原子的能量，此时，需要构建一个边长超过10Å的大胞，里面只放置一个孤立原子。K点只取1x1x1即可，然后做静态自洽计算。单个原子的计算一般需要开磁性参数，INCAR中设置ISPIN=2。详细内容可以查看“大师兄”的教程：https://www.bigbrosci.com/2017/10/30/ex08/

注意，由于晶格的周期性，原子旁边仍有原子，但是该距离超过10Å，一般认为足够。原则上需要测试孤立原子能量随尺寸的收敛情况。