DFT吉布斯自由能理论计算

我第一次做 DFT吉布斯自由能理论计算是算一个析氧中间体的反应自由能，只在 0 K 下比总能，得出某步是放热 0.3 eV。等我把温度拉到 298 K 加上熵项，那一步翻成吸热 0.1 eV，整个限速步判断全反了。那次让我记住，DFT吉布斯自由能理论计算要的是”能拿去和实验温度对标”的数，不是 0 K 的电子游戏。

一、为什么 DFT吉布斯自由能理论计算先做频率

吉布斯自由能 G = E_DFT + ZPE + ΔH_vib + ΔG_trans/rot − TS。除了 E_DFT 来自自洽，其余全靠频率计算。一个气相分子算一次频率就能拿全部分量，表面吸附态则只算与吸附相关的低频模式。曾经我偷懒不给气相分子算频率，用经验熵代替，结果一个分解反应的 ΔG 差出 0.4 eV。做 DFT吉布斯自由能理论计算时，气相分子的精确频率是地基，不能省。

二、核心原理：各分量怎么来

E_DFT 是基态总能；ZPE 是 ½hν 求和；振动对焓和熵的贡献用谐振子公式，高温下还要做高温修正；气相分子还有平动和转动熵，用标准统计力学公式，按 1 bar、指定温度算。pv 项对气相是 RT，对凝聚相可忽略。我用 VASPKIT 或自己脚本把频率文件批量转成热力学量，默认输出 100~800 K 的 G(T) 曲线，因为很多反应的温度依赖恰恰藏在熵项里。

三、关键技术要点：低频模式的误差放大

谐振近似在低频（< 200 cm⁻¹）模式上最不准，而表面吸附的阻碍平移/旋转模式正好落在这里，对熵贡献巨大。我处理时会对这些模式做合理截断或用近似平动模型，否则熵会被谐振公式高估，ΔG 虚低。另一件事是泛函：PBE 的总能系统性偏低，但频率一般还行；做 DFT吉布斯自由能理论计算我默认 PBE+D3，必要时用 HSE06 校总能、PBE 频率，二者混用反而更稳。

四、实操流程：从总能到 G(T)

我的标准动作：结构充分弛豫 → 高精度单点 → 频率计算（气相分子完整、吸附态取相关模式）→ 脚本汇总 ZPE、焓、熵、G 随温度 → 套反应循环算 ΔG_rxn(T) → 和实验平衡常数或电势对标。最容易被砍的是温度扫描，但很多反应在 298 K 和 500 K 的限速步会换人，DFT吉布斯自由能理论计算的价值恰恰在温度维度。

补一个实算：某固相析氧路径，0 K 下第二步 ΔG 是 −0.25 eV（放热），加上 298 K 振动熵和气相 O₂ 平动熵后变成 +0.32 eV（吸热），和实验过电位指示的限速步一致。如果只报 0 K 总能，你会把限速步判成别的一步，后续催化剂改性方向全错。这之后我做机理分析，默认把 G(T) 曲线画出来，不再只贴 0 K 数字。

五、常见踩坑：把总能当自由能

DFT吉布斯自由能理论计算最普遍的错就是只比 E_DFT 就下结论，尤其对气相分子数变化的反应，熵项动辄 ±0.3 eV。另一个坑是忽略溶剂，电化学体系不加隐式溶剂，电势意义全失。还有人把表面态的”假频”（数值噪声产生的虚频）当真频率加进去，ZPE 被污染。

六、复盘：自由能是机理的裁判

回过头看，DFT吉布斯自由能理论计算真正的价值是让机理判断经得起温度考验，而不是给你一个好看的总能差。我接项目默认先问反应条件和温度区间，再决定频率怎么算、熵怎么取。把总能、频率、温度三件事串好，ΔG 才敢写进机理讨论。被证明有用的，是愿意为一根 G(T) 曲线多等两小时算频率的耐心。下次有人甩给我一个”放热”结论，我会先问——你的熵加了吗，温度对吗。