VASP计算分子能量：从孤立分子建模到BSSE校正的全流程

用VASP算分子能量这件事，入门门槛比周期性体系低很多——不需要切表面、选slab模型、调控真空层。但也正因为”简单”，很多细微的陷阱反而更容易被忽视。

第一个绕不开的问题是盒子大小。VASP本质上是平面波基组加周期性边界条件的代码，算分子相当于把一个孤立分子放进一个足够大的真空盒子。”足够大”三个字是操作层面的核心。经验是分子沿各方向的尺寸加至少12Å真空层，总盒子边长不小于18Å。曾经算过一条含苯环的有机分子，初始给了10Å真空层，相邻周期镜像之间的π-π堆叠给总能量贡献了近0.3eV的伪相互作用——看能量曲线觉得不对劲，把真空加到15Å之后才收敛。

K点只需用Gamma点，这一点理解上容易接受，实际操作中却有人忘掉——自动K点网格对分子盒子来说会产生大量冗余计算，而且不会提升任何精度。直接用Gamma-centered 1×1×1。

赝势选择上，对于轻元素分子体系（主要是H、C、N、O），PAW-PBE就够了。但涉及到卤素、稀土、或者对极化敏感的体系时，需要权衡要不要上杂化泛函。HSE06做分子能量的精度确实比PBE好，但代价也是明显的——一个60原子以内的分子，PBE 10分钟搞定，HSE06可能要90分钟。性价比的判断标准很简单：如果你的分子里有电荷转移比较紧张的区域（比如推拉电子对），HSE06值得等；如果只是常规的共价键合分子，PBE加DFT-D3色散修正已经够用。

BSSE校正是真正拉开”硕士算”和”博士算”差距的地方。做过一个二聚体结合能的项目，没做BSSE校正时结合能差了将近0.4eV，做完之后基本跟MP2参考值对上了。VASP的BSSE校正走的是counterpoise方法，需要设置ghost原子的flag，然后用四组单点能相减。步骤多了两步，但时间是值得花的。

有个容易犯的错值得一提：很多人在算分子能量时顺手开了ISIF=3做结构优化，然后忘了在单点能阶段把应力释放掉。分子盒子里的应力是一个纯数值产物——没有任何物理意义——如果INCAR里ISIF没有切回2或者完全不设，CPU会花大量时间在无意义的晶格优化迭代上，而且偶尔会把盒子压到一个偏小的尺寸，间接导致真空层不够的问题。

最后，算分子的自由能需要频率计算。很多人在这一步偷懒不加零点能修正，给出的分子能量其实缺了一个重要的物理项。VASP做频率分析之后在OUTCAR里直接可以查到零点能，把这个加到电子能量上才是可跟实验热力学数据比较的总能。频率计算还有一个副产品——虚频检测。如果虚频数值超过50 cm⁻¹，说明优化得到的那个结构不是势能面上的真极小点，可能是鞍点或者某种中间过渡态，得回到结构优化阶段重新收敛。
还有一个值得强调的收敛性验证步骤：对不同盒子大小做能量外推。具体做法是取四个不同的盒子边长（比如18、20、22、24Å），分别做结构优化加单点能，然后把总能量对盒子大小的倒数做线性外推，取无穷大盒子极限。这个步骤的成本是四倍的单点能计算量，但它能给你一个完全消除了周期性镜像相互作用的分子能量。对于做高精度结合能、吸附能或者反应能垒的项目，这种外推是审稿人眼中的加分项而非多余步骤。

https://www.keyanxueshu.com/category/dft/