说明：本文华算科技主要介绍反应能垒在第一性原理计算中的计算对象、路径能量来源、自由能修正、模型条件和材料机制判断。

反应能垒到底算哪一段能量？

反应能垒对应同一势能面上从初态到过渡区域的能量上坡。初态和终态是两个经过弛豫的局部稳定构型，二者能量差给出反应前后的热力学高低；ΔE^‡ 描述初态到高能构型之间的动力学阻力。在静态 DFT 中，常用写法为 ΔE^‡=E_TS-E_IS，E_TS 来自过渡态或路径最高候选构型，E_IS 来自相同表面模型中的初态。

能量差的数值归属由晶胞尺寸、slab 厚度、表面覆盖度、电荷态、磁态、溶剂近似和参考态共同限定，任一条件改变后，E_IS 与 E_TS 所在的势能面都会改变。同一模型条件限定了能量基准，反应路径对象限定了高能构型的化学含义；缺少其中任一项，能垒数字很难对应到具体反应步骤。

表面反应中的能垒通常来自键长、配位和吸附位点同时变化的区域。以 Ni/ZrO₂ 氢溢流为例，H₂ 解离、H* 扩散和氧空位形成分别对应不同初态、过渡构型和终态；H-H 键拉长、H-金属键生成、O-H 键形成会把体系推到较高能量位置。E_a 标的是这段上坡高度，E_r 标的是两端稳定态的能量差。

自由能垒把静态电子能量换成有限温度和环境基准下的 ΔG^‡。电子总能差仍是起点，零点能、振动热容、熵、气相分压、液相标准态和电化学势会改变过渡区域与初态的相对高度。ΔE^‡、ΔH^‡ 和 ΔG^‡对应的物理量不同，跨论文比较时应先看作者使用的能量定义。

路径能量怎样从计算里得到？

路径能量来自连接两个稳定态的构型序列。表面扩散、离子迁移和吸附物重排常采用 NEB 或 CI-NEB，把初态和终态之间的中间图像放到势能面上弛豫；分子反应或表面键转化常结合 dimer 方法、约束扫描、过渡态优化和频率分析。路径搜索定位候选高能区域，过渡态优化和振动模式再检验该构型是否接近一阶鞍点。

频率分析给过渡态命名提供直接证据。目标反应坐标方向出现一个虚频，且振动模式对应正在断裂或形成的化学键，才能把候选高能构型和目标步骤联系起来；多个虚频、表面整体滑移虚频或吸附物转动虚频会改变构型归属。虚频方向比“最高图像”这个标签更接近反应坐标本身。

反应坐标通常取自关键键长、原子迁移距离、质子转移距离或吸附物构型变化。N-H 断裂、O-H 生成、Li⁺ 跳迁和 *COOH 形成对应的坐标不同，过渡构型附近的力和曲率也不同。坐标选择决定了路径搜索是否沿着目标化学事件推进。

覆盖度、局部配位和反应气氛改变时，同一反应式对应的势能面会同步改变。相邻 H*、OH* 或 CO* 之间存在排斥、氢键或金属团簇重排时，过渡构型的键长和局域电荷分布随之改变；高覆盖度下的活化能可能低于低覆盖度，也可能因位点阻塞而升高。覆盖度依赖的 E_a比单个稀覆盖度数字更接近工作表面的反应图像。

在反应路径图中，E_a、E_r 和 ΔG 分别对应上坡高度、两端能差和热力学校正。E_a 量化初态到过渡区域的上坡高度，E_r 量化两端稳定态的能量差，ΔG 把热力学修正和环境化学势纳入反应坐标。终态稳定只能说明反应热力学驱动力较大；正向能垒和逆向能垒仍由各自初态到同一过渡区域的高度决定。

初态和终态的几何构型、吸附位、质子化态、磁态和电荷态直接进入能垒数值。初态弛豫不充分时，E_IS 偏高或偏低；终态吸附位、质子化态、磁态或电荷态选错时，路径会连接到另一类反应对象。表面位点、覆盖度、电子态和反应式配平共同限定能垒的可比范围，路径算法只是这个化学定义之后的数值步骤。

为什么静态能垒会转成活化自由能？

静态 DFT 总能描述 0 K 近似下的电子结构能量，真实催化反应处在有限温度、气氛、电位和溶液环境中。气相分子的平动、转动和振动熵使吸附前后自由能差偏离总能差；吸附态的自由度受表面束缚，振动频率对 ZPE 和熵项贡献较大。ΔG^‡ 把过渡构型和初态放到同一温度、压力或电位基准下，反应速率讨论通常使用这个量。

电化学路径中的活化自由能还受电子化学势控制。CO₂ 还原、HER、NO₃^– 还原和 ORR 中，质子-电子转移改变中间体电荷、键级和吸附强度；外加电位移动电子能级，也改变 *COOH、*H、*NO₂ 或 *OOH 的相对稳定性。恒电位自由能图把反应坐标、位点和电位写在同一个能量标尺上。

电极附近的水分子取向、阳离子位置和局域电场会改变过渡区域的稳定程度。水分子氢键网络、碱金属阳离子、水合质子和局域电场可稳定 *OOH、H₂O₂ 脱附态或含氧中间体，也可抬高 4e^– ORR 的竞争路径。界面水结构进入自由能图后，同一催化位点在不同电解质中会呈现不同的有效能垒。

气相小分子、固体表面中间体和电极界面物种对 ZPE、熵项、溶剂化项的敏感来源并不相同。气相小分子的熵项常主导吸附自由能变化，固体表面中间体更多受振动频率和几何约束影响，电极界面叠加双电层、局域 pH 和离子分布。同一电子能垒可对应不同 ΔG^‡，原因来自温度、标准态和环境模型，而非单纯数值噪声。

当两个材料的能垒差只有 0.05-0.10 eV，修正项和模型误差已经足以改变排序。低频振动模式、溶剂化近似、带电中间体校正、表面电荷补偿和吸附物覆盖度都会影响 E_TS-E_IS。能垒排序在小差值区间应同时报告模型条件和修正方式，避免把数值第三位小数写成材料优劣。

含过渡金属位点、氧空位或自由基中间体的路径中，自旋态和电荷局域方式会改变过渡态构型。不同自旋解可能给出不同 E_TS；若初态和过渡态的电子态不在同一分支上，能垒差值会混入电子态切换误差。

能垒数值怎样服务材料机制判断？

反应能垒进入材料比较时，描述的是指定模型中某个步骤的动力学阻力。0.2 eV 的降低在 Arrhenius 因子中可能带来显著速率差异，前因子、位点数量、路径简并度、覆盖度和传质仍会进入实际速率。低 ΔG^‡支持该步骤在给定模型下动力学阻力较小，宏观活性仍由完整反应网络和工作条件共同决定。

催化循环由多个吸附、转化、脱附和位点再生步骤组成，单步能垒低不必然对应整体速率高。C-H 断裂或 O-H 形成若只是循环中的快步骤，前一步中间体生成或后一步产物脱附仍可能控制通量。DRC、SRC 和微观动力学模型把多个 ΔG^‡、覆盖度和电位响应写入速率方程，速率控制权重由整个反应网络分配。

电子结构图谱用于解释能垒变化的来源。d 带中心移动、Bader 电荷变化、差分电荷密度重排或 PDOS 杂化增强，可说明过渡构型中吸附键、金属-吸附物轨道耦合和局域电荷的变化。过渡态稳定化幅度要与初态稳定化幅度同标尺比较；若催化剂同样稳定反应物，E_TS-E_IS 未必下降。

实际材料筛选里，能垒结果通常与反应热、吸附自由能、覆盖度相图、电子结构和稳定性指标共同使用。对于 SnO₂(110) 上 CO₂ 还原这类竞争反应，Sn4、Sn5 位点的自由能图给出候选路径，DRC/SRC 指出速率和选择性敏感步骤，HER 分支提供副反应参照。反应能垒限定具体步骤的动力学权重，表面结构、位点数量和反应网络仍保留各自的判据位置。