说明：本文华算科技主要介绍零点能校正的物理来源、频率求和方式、在轻元素材料体系中的能量作用，以及它怎样进入吸附、缺陷和反应路径判断。

在 DFT 计算里，结构优化把原子放到势能面附近的极小区域，原子核并没有真的停在一个没有动量的数学点上。量子力学把每个正常振动模式近似成谐振子，最低能级仍然保留 1/2hν 的基态振动能，这部分能量就是零点能。

零点能校正把静态电子总能 EDFT 与振动基态能连到同一能量标尺上。科研计算里真正使用的通常不是某个结构的绝对零点能，而是两个态之间的差值：ΔEZPE=EZPE,product-EZPE,reactant。同一个反应式两侧如果含有相近的重原子骨架，大部分低频模式会抵消；H、D、O-H、N-H、金属-H 这类高频模式的差异会留下更清楚的能量位移。

零点能校正的基本概念是什么？

零点能校正的起点不是温度，也不是软件输出里的附加标签，而是原子核振动的量子化。经典图像中，粒子在势能最低点可以没有动能；量子图像中，位置和动量无法同时压到零，最低振动态仍有有限能量。对材料计算来说，这个能量来自力常数矩阵给出的振动频率，和电子自洽总能属于不同来源。

把一个吸附态、缺陷态或过渡态候选构型优化好以后，DFT 总能给出电子基态下的势能面高度。零点能校正相当于在这个高度上加一层振动基态台阶；台阶高低由各个模式的频率决定。低频金属骨架模式每个只贡献很小的 hν/2，高频 X-H 伸缩振动却能带来几十到几百 meV 的差别。

图1. 谐振子势阱中的零点能位置。

把“校正”理解成修补计算错误时，ZPE 的物理含义会被放到错误位置。它更接近把模型从纯静态电子能量扩展到振动基态能量：结构坐标没有重新优化，电子密度没有重新求解，能量表达多了由频率带来的项。同一份频率输出里，实频数量和模式方向 应和结构身份对应，再谈 ZPE 数值。若一个结构本身有虚频，零点能并不会把它变成稳定构型；虚频模式反映势能面方向没有停在极小点，应回到结构和路径本身处理。

这种校正只在比较里显出意义。单个 slab 的零点能数值可以很大，但两种相似表面的重原子振动会互相抵掉；同一个位点上 H*、*OH、*OOH 的成键环境不同，频率分布不再相同，剩下的 ΔEZPE 就会改变吸附自由能或反应能垒。

频率怎样变成零点能校正项？

频率计算把一个优化结构附近的势能面局部曲率写成 Hessian 矩阵，再把原子质量带入质量加权坐标，得到正常振动模式。非线性分子有 3N-6 个振动模式，线性分子有 3N-5 个；表面吸附体系常把吸附物和邻近可动原子的模式取出来近似处理，固定的底层原子不应被拿来解释吸附物的高频振动。

图2. 不同振动模式对零点能求和的贡献。

在谐振近似下，每个真实频率 νi 贡献 1/2hνi。若频率用 cm-1 表示，换算到 eV 后求和即可得到该构型的 EZPE。反应能校正不是把反应物和产物各自的 ZPE 单独展示，而是计算 ΔErxn,ZPE=ΣEZPE(产物)-ΣEZPE(反应物)；能垒校正则比较过渡态候选构型和初态构型的零点能差。

把频率项放进反应自由能时，平动、转动、熵项和热焓项经常和 ZPE 一并出现。气相小分子通常按分子热力学处理，吸附态则多把平动、转动自由度冻结成低频振动；低频模式如果接近数十 cm-1，谐振近似会放大熵项的不确定，而 ZPE 本身主要受高频模式支配。零点能校正不等同于完整热力学校正，它只是自由能表达式中的振动基态部分。

频率里的质量项让同位素校正有了直接来源：在近似相同的力常数下，ν 与 1/√m 成比例，D 相关伸缩振动低于 H 相关伸缩振动，EZPE 随之降低。把 H2 换成 D2 时，静态电子总能几乎不变，零点能差和振动熵 才是同位素自由能差的主要来源；氢迁移、质子穿膜和同位素筛分都要面对这种质量效应。

图3. HX 在固相、液相、气相和水合状态之间的热力学循环。

开放文献中常把气相分子的 GiDFT 写成 EiDFT、ZPE、0 K 到 298.15 K 的焓变和 TS 项的组合。对电催化反应，H2、H2O、NH3、CO2 这类参照分子一旦进入反应式，ZPE 与参照态误差会一起移动自由能台阶；若文章只报静态 ΔE，读者很难判断氢化步、脱氢步或质子电子转移步的能量排序。

ZPE 校正影响哪些反应能量？

零点能校正对重元素骨架的相对排序通常较小，对含氢、含质子转移、含轻原子伸缩振动的结构更敏感。HER 的 H* 吸附、ORR/OER 中的 *OH 与 *OOH、CO2RR 中的 *COOH、氨合成中的 *NHx，都把 X-H 频率带进吸附自由能。几十 meV 的 ΔEZPE 足以改变两个位点接近时的排序。

缺陷计算里，氢、锂、质子和小分子插层的局部振动更值得检查。氢在硅、氧化物或二维材料中可能从键中心位、反键位、空位附近位点切换到另一个局域环境，频率谱会随局部配位改变。此时零点能不只是一个统一常数，而是 每个构型自己的振动指纹。

图4. I4-H 复合缺陷中氢原子的不同局部位置和相对能量。

这类氢缺陷图并不是 ZPE 图谱本身，却提醒我们一件事：ZPE 来自结构附近的振动曲率。氢放在 I4 缺陷笼内、键中心位或远离缺陷的晶格位置，局部 Si-H 约束不同，振动频率也会分开。若两个构型静态形成能只差 0.05-0.10 eV，忽略轻元素 ZPE 可能把候选态顺序写反。

图5. I4 缺陷的带隙能级和局域轨道分布。

把零点能校正写进文章时，还要分清它和电子结构图谱的角色。PDOS、带隙缺陷能级、差分电荷密度描述电子态和成键，ZPE 描述原子核在局部势阱中的基态振动。电子结构解释活性位为什么形成，ZPE 修正这些态之间的能量差；两者互相补位，却不能互相替代。

过渡态如何开展 ZPE 校正？

反应路径里最常见的写法是对初态、中间态、终态和过渡态候选构型分别做频率分析，再把 ΔEZPE 加到静态反应能或静态能垒上。对氢迁移、质子耦合电子转移、H2/D2 同位素筛分这类过程，参与反应坐标的轻元素模式 往往决定校正方向。

过渡态频率要格外小心。严格的过渡态应有一个沿反应坐标的虚频，其余模式为实频；ZPE 求和通常排除这一个虚频，只对正频模式求和。若过渡态候选结构出现多个虚频，问题在路径或构型上，不是把虚频取绝对值再求和。NEB 最高图像也只是路径上的候选点，若要给出 ZPE 校正后的活化能，仍应对爬山图像或独立过渡态优化结果做频率检查。

图6. 不同泛函给出的共电解反应自由能与实验参照差异。

自由能图中的每个台阶都像一本能量账。EDFT、ΔEZPE、热焓、熵、电位、pH、溶剂化和气相参照误差分别改变账本里的不同项目。上面的反应能误差图给出一个很实际的提醒：即便泛函层级提高，分子参照和修正项仍会带来 eV 量级偏差；在这种背景下，ZPE 校正应写成可追踪的能量项，而不是一句“已校正”带过。

读一篇材料计算文章时，可以沿着三个问题核对 ZPE：哪些原子参与高频振动，哪些结构之间在比较，校正后的差值相对静态能量改了多少。若 ΔEZPE 只有 0.01 eV，而静态差值是 0.5 eV，它通常不会改写结论；若两个位点或两个路径只差几十 meV，ZPE、同位素质量和低频模式处理 就会进入判断核心。

把校正后的能垒用于速率讨论时，数值尺度要和 kBT 换到同一张能量尺上。室温下 kBT 约为 0.026 eV，0.05 eV 的 ZPE 差已经能让 Arrhenius 因子出现成倍变化；若路径判断停在 0 K 静态能垒，轻元素转移步骤的相对快慢 可能会被低估或高估。