过渡态最容易被误解成“反应过程中的某个中间结构”。这个说法只说对了一半:过渡态确实位于反应从初态走向终态的途中,但它不是可以稳定停留的中间体,而是势能面上的一阶鞍点。过渡态是连接初态和终态的最高能鞍点,不是能被单独分离出来的稳定物种。
在催化计算里,是否需要过渡态,取决于你要回答的问题。自由能台阶能说明热力学上哪一步上坡,过渡态能垒则说明跨过这一步需要多大活化代价。一个步骤热力学上有利,也可能因为过渡态太高而很慢。

图1 活化能、反应坐标和过渡态位置示意,展示过渡态位于初态到终态之间的能量峰附近。DOI: 10.3390/ijms26104670。


过渡态到底是什么?
为什么它是鞍点?
势能面可以想象成一个多维地形。初态和终态通常在局部低谷,过渡态位于连接两个低谷的山口。沿反应坐标方向,它是最高点;沿其他垂直方向,它又像局部低点。这种“一个方向向下、其余方向尽量稳定”的点,就叫一阶鞍点。
这也是过渡态不同于普通优化结构的地方。普通结构优化会把体系推向附近低谷;过渡态搜索则要找到通往终态的山口。如果一个构型沿所有方向都稳定,它更像中间体;如果只有一个方向不稳定,它才可能是过渡态。
Ea = ETS − EIS 式(1)
式(1)里 Ea 是从初态 IS 跨到过渡态 TS 的活化能。若反应要从终态反向走回初态,也可以用 ETS − EFS 表示反向能垒。能垒大小直接影响反应速率,因此它和自由能差不是同一个概念。
为什么它不能当成中间体?
中间体对应势能面上的局部最低点,轻微扰动后还会回到附近;过渡态在反应坐标方向不稳定,轻微推一下就会滑向初态或终态。实验上常说“捕捉过渡态”也通常指间接证据、类似结构或过渡态类似物,而不是像稳定吸附物那样长期存在。
在 VASP 里,如果你把过渡态结构拿去做普通弛豫,它往往会掉向初态或终态。这个现象不是计算坏了,而是过渡态本身就不该是稳定低谷。过渡态的身份来自反应路径和曲率,而不是来自一张静态结构图。

图2 反应路径中初态、过渡态和终态结构及自由能示例,说明过渡态需要和具体反应坐标绑定。DOI: 10.1021/jacsau.2c00582。


NEB为什么能找过渡态?
图像串在找什么?
NEB 的思路是先给出初态和终态,再在两者之间插入一串图像。每个图像代表反应路径上的一个近似结构,通过弹簧约束和真实力投影,让整串图像逐渐贴近最低能路径。最高能附近的图像,就是过渡态候选。
CI-NEB 会把最高能图像推向更准确的鞍点。它不只是把某个结构能量算高,而是让最高能图像沿路径方向上爬、沿垂直方向下滑。NEB 找的是连接初态和终态的最低能路径,而不是随手插值的一排结构。
初态和终态为什么要先优化好?
如果初态或终态本身没有优化到局部低谷,NEB 路径会被错误端点牵着走。比如吸附物在终态其实应该转到桥位,但终态只放在顶位没有弛豫,NEB 得到的能垒可能对应一条不存在的路径。
端点还要保持同一超胞、同一原子数、同一固定层和相同原子排序。水解离、H 迁移、O-O 成键、C-H 断裂这类反应,都要先确认初末态是同一反应的两端,而不是两个无关构型。端点不清,过渡态再漂亮也没有反应意义。

图3 约束动力学得到的反应坐标和平均力示例,可帮助理解能垒与反应坐标的关系。DOI: 10.1073/pnas.1604590113。


单虚频为什么重要?
虚频在验证什么?
过渡态验证常用振动频率。对真正的一阶鞍点,频率分析应出现且只出现一个虚频。这个虚频对应的振动模式,应沿着目标反应坐标运动,例如 H 从一个位点转移到另一个位点,O-H 键断裂,或 C-H 键被拉长。
如果没有虚频,结构更像局部最低点;如果出现多个虚频,说明结构还没有收敛到一阶鞍点,或者路径包含不止一个不稳定方向。单虚频不是装饰性检查,而是过渡态身份的关键判据。
ν² ∝ k / μ 式(2)
式(2)只是帮助理解频率和曲率的关系:k 可理解为势能面在某个方向的曲率,μ 是相关原子运动的约化质量。虚频对应负曲率,说明沿这个方向结构会向两边滑落。
虚频模式必须和目标反应一致吗?
必须一致。若你想找 H 转移过渡态,虚频却是表面某个无关原子摆动,那这个结构即使只有一个虚频,也不能说明目标反应的过渡态已找到。频率动画或模式矢量要能把结构推向初态和终态。
对表面催化体系,低频摆动、整体平移和软模式容易干扰判断。可以先固定底层,保证吸附物和活性位附近原子参与频率;必要时只对反应区做频率验证。虚频数量和虚频方向要一起看,不能只看日志里有一个负数。

图4 C-H 键断裂反应在不同约束键长下的结构快照,展示反应坐标可对应具体键长变化。DOI: 10.1073/pnas.1604590113。


过渡态和自由能台阶怎样分工?
热力学上坡等于能垒高吗?
不等于。自由能台阶比较的是初态和终态的稳定性差异,过渡态能垒比较的是从初态跨到山口的代价。一个步骤终态更稳定,仍可能因为山口很高而速率慢;一个步骤终态略高,也可能因山口不高而可逆发生。
在 HER、OER、CO₂RR、NRR 等体系中,吸附自由能图常用于快速筛选,但涉及水解离、O-O 成键、N-O 断裂、C-H 活化时,单看台阶往往不够。自由能图回答能不能在热力学上走通,过渡态回答跨过去快不快。
哪些步骤更需要算过渡态?
只涉及质子电子转移且采用 CHE 近似的入门筛选,有时可以先用自由能台阶判断。但如果步骤包含明显的化学键断裂、成键、分子重排、表面迁移或水分子参与,就应考虑 NEB 或其他能垒方法。
比如碱性 HER 的水解离、CO₂ 变成 *COOH、NO₃RR 的 N-O 键断裂、OER 的 O-O 成键,都可能有显著过渡态。此时只说吸附能更优,会漏掉真正控制速率的山口。

图5 O-H 键断裂反应在不同约束键长下的结构快照,说明过渡态分析常要围绕具体断键或成键坐标展开。DOI: 10.1073/pnas.1604590113。


VASP里怎样避免把TS算偏?
第一,初态和终态先充分优化,并确认原子排序一致。第二,插值后检查每个图像是否有原子过近、吸附物穿过表面或路径绕远。第三,NEB 收敛后看能量曲线是否平滑,最高图像是否合理。第四,对最高图像做频率验证,确认单虚频方向对应目标反应。
SYSTEM = neb_ts_check
IMAGES = 5
SPRING = -5
LCLIMB = .TRUE.
IBRION = 3
POTIM = 0
EDIFF = 1E-5
EDIFFG = -0.03
NSW = 200
# 参考示例:端点先优化,CI-NEB 找最高图像,再用频率确认单虚频模式
最终写结果时,不要只报告“找到 TS”。更清楚的表达应包含:初态是什么,终态是什么,TS 对应哪根键形成或断裂,正向和反向能垒是多少,单虚频模式是否沿反应坐标。这样过渡态才真正服务于反应机理判断。
还要区分电子能垒和自由能能垒。NEB 直接给出的通常是 0 K 电子能路径,若要和实验温度下的反应速率联系,还需要考虑 ZPE、熵、溶剂环境和电位修正。对于表面反应,吸附物的低频模式和界面水结构会影响自由能修正,因此能垒数值最好说明采用的是电子能还是自由能口径。
过渡态也不一定只有一条路径。一个反应可能存在顶位迁移、桥位迁移、邻位水辅助和双位点协同等多条路线。若只算一条人工插值路径,就把它写成唯一机制,结论会偏窄。合理的 TS 分析应从多个可能端点和路径中找最低可行山口。
当 NEB 曲线出现多个峰时,还要判断是否包含了多个基本步骤。比如吸附物先旋转、再断键、最后迁移,可能需要拆成两段或三段路径分别计算。把多个步骤硬塞进一条长 NEB,会让最高点难以解释,也不利于虚频验证。
