DFT计算能量：相对能量比绝对能量重要，前提是参考系选对了

一个CO在Pt(111)表面吸附能的计算项目，最终提交的数据和实验值差了8 kcal/mol——不是可以忽略的偏差。回过头查原因，三条陷阱全踩了：BSSE没校、零点能用了谐振近似硬套、O₂参考态的自旋多重度没设成3。DFT计算能量的输出文件里，”SCF Done”后面的数字动辄几百到几千Hartree，这个绝对值本身没有直接化学意义——化学上所有有用的结论，结合能、反应能、能垒，全来自能量差。而能量差的可靠性，取决于相减的两项是否在同一精度水平上、参考态是否可比。这个项目恰恰在这两件事上都出了问题。

BSSE校正是弱吸附体系的硬性门槛

结合能E_ads = E_slab+molecule – E_slab – E_gasmolecule，方程看着直白，第一道陷阱藏在基组里。孤立分子在气相优化时只用自身基组，到了表面上却额外”借用”了表面原子的基函数——基函数重叠让表面上分子的能量被人为算低，吸附能随之高估。Counterpoise校正（Boys & Bernardi方法）用表面+分子的全基组重新计算孤立分子的能量，把基组不一致的缝隙堵住。弱吸附体系（物理吸附，结合能<0.5 eV）不做BSSE校正，结论几乎不可信——偏差能占到总结合能的30-50%；强化学吸附（>1 eV）的校正量通常不到0.2 eV，可以省下这步计算时间。这个项目计算的是CO在Pt上的化学吸附，BSSE校正量只有0.08 eV——单独看不算致命，但叠加另外两个偏差后，总误差就不可接受了。认定BSSE校正在这个体系里”可以不做”，是被化学吸附的经验结论误导了——叠加效应让每一条小偏差都变得不可忽略。

零点能的谐振近似是个粗口径

第二个陷阱是零点能（ZPE）。振动频率分析输出的零点能直接加在电子能量上，但谐振近似对低频模式存在系统性高估。表面吸附体系中，分子的平动和转动自由度转变为振动模式，新增一堆低于100 cm⁻¹的软模——这些模式的DFT频率误差可达50%，零点能跟着偏。结合能计算中ZPE贡献量级在0.1-2 kcal/mol，定性判断不校正也能过关，但定量比较不同金属表面的吸附能排序时，ZPE差异足以反转排名。这个项目比较了CO在Pt和Pd上的吸附偏好——不做ZPE校正时Pt比Pd强0.3 kcal/mol，校正后Pd反超Pt 0.1 kcal/mol。排序直接翻转，差距不会说谎。谐振近似在软模密集的体系里是个必须被标注的已知缺陷，不是可以默默跳过的技术细节。

开壳层参考态的自旋多重度，杀伤力最大的隐蔽陷阱

第三个陷阱杀伤力最大，也最容易遗漏。气态H₂O、CO₂这些参考试剂的DFT能量是基态闭壳层，没问题。但O₂基态是开壳层三重态——用闭壳层（自旋多重度=1）算出来的参考能量差了5-10 kcal/mol。算氧化反应能时不指定O₂的自旋多重度为3，得到的反应焓在物理上是错误的。这个项目在计算CO氧化反应能时，O₂的spin multiplicity默认设成了1，差值4.8 kcal/mol直接灌进了反应能。这类错误在文献里反复出现，不是个案而是通病，值得警醒。每一条涉及O₂、NO、NO₂等开壳层分子参与的反应能计算，自旋多重度是第一项需要核验的设置。

热力学校正：从电子能量差到自由能差

反应能和能垒的计算中，热力学校正把”电子能量之差”变成”自由能之差”。Gaussian的freq关键词做完频率分析后输出热力学校正值，但默认针对气相298 K和1 atm标态。如果反应在溶液里跑的、温度是400 K、浓度不是1 M——这些默认值需要手动调整或用脚本处理。VASP和Gaussian的绝对能量没有可比较性（PAW势 vs 全电子/赝势），但同一计算设置内部的能量差通常能对齐到0.1-0.2 eV——前提是k点密度和平面波截断能都收敛到位，这份收敛性检查本身是计算成本的必要组成部分。

回过头看这个项目，三条陷阱叠加后总偏差8 kcal/mol，每条单独看都不致命——BSSE只差0.08 eV，ZPE只差0.3 kcal/mol，自旋只差4.8 kcal/mol——合在一起就把结论推翻了。DFT计算能量的每一步减法，背后都有一个参考系的假定；踩准了，差值才有意义，踩偏了，数据再精确也是精确的错误。