VASP计算形成焓：方法原理与参考能选取的关键判断

材料的热力学稳定性是计算材料学中最基础的问题之一，而形成焓是判断稳定性的核心量。VASP计算形成焓的流程并不复杂，但参考能的选取逻辑常常让人困惑——同样一个体系，不同人计算出来的形成焓可能差出几十到几百 meV/atom，这背后不全是误差，而是物理约定的不同。

形成焓的计算定义

形成焓 ΔHf 定义为：化合物从其组成单质直接反应生成时的焓变，归一化到每个原子或每个化学式单元。

以二元化合物 AₘBₙ 为例：

ΔHf = [E(AₘBₙ) – m·E(A) – n·E(B)] / (m+n)

其中：

• E(AₘBₙ) 是化合物的总能（VASP优化后的TOTEN）
• E(A)、E(B) 是单质参考态的总能，归一化到每个原子

这个公式看起来简单，真正的讲究全在 E(A) 和 E(B) 的选取上。

参考态选择：最常见的分歧来源

气态非金属元素是形成焓参考能最容易出问题的地方。氧、氮、氢在标准状态下是双原子分子（O₂、N₂、H₂），VASP计算这些分子时的精度与固体体系存在系统性偏差。

以氧化物为例，标准做法是计算 O₂ 分子在超胞中的总能作为氧的参考能。但DFT（PBE泛函）对 O₂ 分子的描述存在已知误差，导致氧化物形成焓系统性偏低约0.1-0.3 eV/O。解决方案有两种：

1. 实验校正：用已知形成焓的参考氧化物（如 MgO、Al₂O₃）对 O₂ 参考能做化学势校正，补偿DFT的系统误差
2. 采用拟合参考能：Materials Project 数据库使用了一套经过实验形成焓拟合的元素参考能，能显著降低误差

对于与 Materials Project 数据库对比或用于相图计算的形成焓，强烈建议使用与 MP 一致的参考能体系，否则计算结果无法直接用于相平衡分析。

金属单质的参考能相对简单，直接使用对应体相的稳定晶型（Fm-3m铜、BCC铁等）的VASP总能即可。注意每个原子的能量需要对超胞中的原子数做归一化。

DFT+U对形成焓的影响

含有 3d（Fe、Co、Ni、Mn、Cu等）或 4f（稀土）过渡金属的化合物，通常需要引入 Hubbard U 校正来修正 DFT 对局域化 d/f 电子的描述。开启 DFT+U 后，化合物的总能会发生变化，但相应的单质参考能也必须在 相同的 U 值设置下 重新计算。

一个常见错误是：化合物用了 U=4 eV，但单质铁的参考能却用的是纯 DFT（U=0）结果，两者的总能在系统论上不一致，形成焓的误差可能高达0.5 eV/atom。

注意：铁等金属单质的 DFT+U 计算通常需要特别处理，因为 U 值对金属中巡游电子的描述并不总是合适，可能导致单质的磁序或晶格常数出现偏差。这是一个没有完美解法的权衡问题，通常的处理是参考文献中对特定体系的推荐方案。

van der Waals 校正的一致性要求

层状材料或含有分子吸附的体系往往需要 vdW 校正（DFT-D3、DFT-D3BJ、optPBE-vdW 等）。与 DFT+U 的原则一致：化合物和所有单质参考的计算都需要使用 相同的 vdW 校正方案。

石墨烯、层状 MoS₂ 等体系的层间结合主要靠 vdW 相互作用，如果化合物层间开启了 D3 校正，石墨参考能（碳的参考态）也必须用 D3 校正，否则形成焓结果会系统性偏差。

结果验证与对比

计算形成焓的最直接验证是与实验数据（JANAF热化学数据库、FactSage、Materials Project实验值）对比。对于过渡金属氧化物，DFT（PBE）的计算误差通常在 ±0.1-0.3 eV/atom 范围内；引入 U 校正后，对特定氧化物体系可以压缩到 ±0.05-0.1 eV/atom。

如果计算结果偏离实验值超过 0.5 eV/atom，首先排查的是：参考能选取是否一致、收敛是否充分（特别是磁性构型的收敛）、以及晶体结构是否使用了正确的稳定多型。

科研学术网提供热力学稳定性全套计算，包括形成焓、相图构建与缺陷形成能，欢迎访问 https://www.keyanxueshu.com 获取技术支持。