DFT计算反应的吉布斯自由能，是连接电子结构计算与实验热力学数据的核心桥梁。很多团队在用VASP或QE算催化反应路径时，只算了电子自由能（E_DFT），就直接拿去和实验的ΔG对比——这个跳跃通常会导致0.2-0.5eV的系统偏差，在某些体系里甚至能颠覆整个反应机理的判断。
裸的DFT总能量E_DFT，对应的是T=0K、无振动、无熵贡献的电子基态能量。要逼近实验条件下的吉布斯自由能G(T)，至少要把三样东西加进去：零点能（ZPE）、振动熵（TS_vib）、平动/转动/电子熵（对气相分子）。
具体公式是：
G(T)=EDFT+ZPE+∫0TCvdT−TStrans−TSrot−TSvib−TSelec
听起来繁琐，但VASP的用户可以借助vibrational.py（来自VASPKIT或Phonopy后处理脚本）把这套流程半自动化。关键是：频率计算必须在收敛精度更高的标准下做——EDIFF至少设1e-6，EDIFFG设-0.01，否则算出来的振动频率噪声会直接污染熵修正。
固-气反应（比如HER、OER、CO2RR在电极表面的反应）里，气相反应物的熵贡献往往被低估或错误处理。一个典型的坑：直接用标准状态下的实验熵值（比如S°_H2 = 130.7 J/mol·K）而不做温度插值，在低温（<200K）或高温（>800K）条件下会引入不可忽视的误差。
更严谨的做法是：用DFT算气相分子的频率，自己算S_vib(T)，然后和NIST-JANAF表里的实验熵做交叉验证。我通常在300K先对比一下，偏差能控制在5 J/mol·K以内，说明频率算得靠谱，再往其它温度推。
如果反应发生在溶液相，G(T)还得加上溶剂化自由能ΔG_solv。VASP用户有两个主流选择：
LSOL=.TRUE.，适合均相体系，计算成本低，但对特定离子-溶剂相互作用描述粗糙。实际项目中，我倾向于先用VASPsol快速扫一遍整个反应网络，把ΔG超过1.0eV的能垒点挑出来，再对这几个关键态做explicit solvent的精细计算。这样在精度和效率之间能取得一个合理的平衡。
很多DFT热力学研究只做了300K一个点的修正，但真实催化反应的工作温度往往在室温到600K之间变化。热容C_v(T)的贡献（公式里的∫C_v dT项）在高温下不可忽略——对气体分子，这项在600K时能达到0.1-0.2eV的量级。
幸运的是，Ideal gas的热容有解析表达式（Einstein模型和Debye模型的组合），很多后处理脚本（比如ASE的Thermochemistry模块）已经内置了这部分计算。核心输入只有三个：振动频率列表、分子质量、对称性数。
本文涉及的脚本和后处理工具，可根据具体体系的需求进一步定制。吉布斯自由能修正不是”可选项”，而是任何声称能与实验对标的DFT热力学计算的必要步骤。
