过渡态理论全解析：高维构型空间含义、反应路径分叉、有限温自由能瓶颈及反应坐标检验

说明：本文华算科技主要介绍过渡态在高维构型空间中的物理含义、多条反应路径的分叉、有限温自由能瓶颈，以及反应坐标的检验方法。

过渡态在构型空间里是什么？

在玻恩—奥本海默近似下，电子结构计算为每一组核坐标给出能量，所有原子坐标共同张成构型空间。一个含 N 个原子的体系通常由 3N 个笛卡尔坐标描述；平移、转动、晶胞和约束条件会改变实际自由度。势能面是能量对整组核坐标的函数，常见的一维反应曲线只是它在某条路径上的投影。

局部极小值对应能够停留的反应物、中间体或产物构型。连接两个相邻势阱的最低能路径经过一个驻点，驻点附近有一条方向向两侧降低，其余独立方向向上升高。零温势能面上的一级鞍点便是基元步骤的过渡态构型。原子受力接近零只满足驻点条件；下降方向确实通向指定的两个势阱，才能建立几何连接。

将高维势能面沿选定路径投影到一条反应坐标后，曲线峰顶对应的构型随路径改变，纵向能差必须从该步实际出发的势阱计算。若体系先形成吸附复合物，再越过成键鞍点，活化能应从吸附态起算；直接拿气相反应物能量作零点，会改变数值含义。

图1b中的峰值适合表达单步能垒，图1c则把多个中间态与路径连成网络。同一净反应可穿过不同节点，每条边各有自己的过渡态。反应式只限定起点、终点和计量关系，无法单独指定原子按何种次序移动。

过渡态构型还包含集体运动。水分子的反对称伸缩要求两个 O—H 键协同变化；表面质子转移可能牵动氢键链、吸附物取向与邻近溶剂。参与反应的坐标常由多项原子位移线性或非线性组合而成，单个键长只在其他自由度能快速随动时才够用。

一个基元变化为何会分出多条路径？

原子计量相同，局部环境仍可给出多种重排顺序。分子能选取不同构象，晶体缺陷能沿不等价晶向迁移，催化表面存在台阶、空位与覆盖度差异，质子还能借助不同数目的水分子接力。路径身份由中间体序列和成键变化共同确定，不能只按产物名称归类。

Ga³⁺/SiO₂ 界面上的乙烯低聚反应展示了这种分叉。网络节点包含 Ga-乙基、Ga-乙烯基、Ga-丁基等吸附物，连边颜色记录各路径的活化自由能。低能垒路径形成主循环，较高能垒支路仍可能生成异构化、裂化或结焦前驱物，温度与停留时间会改变其通量。

比较路径时，最高峰的绝对能量不能直接代替控速信息。每个峰要减去其前方能够积累的中间体自由能；若深势阱占据大部分位点，后续一个中等高度的峰也可能控制周转。局部能垒、覆盖度和可逆回流共同决定稳态速率。

图4把低聚、异构化、裂化与炔烃生成路径画在同一能量尺度上。周期模型和簇模型对若干峰高给出相近趋势，长程骨架约束与位阻仍造成数值偏差。过渡态排序会随表面模型变化，活性位配位、晶胞范围及共吸附物均属于路径定义的一部分。

自旋交叉、电荷转移和激发态反应还会涉及多张电子势能面。沿单一基态势能面找到的鞍点只描述绝热路线；若电子态在反应区发生交叉，控制跃迁的结构可能靠近最低能交叉点。电子态身份应与几何路径一并追踪，否则相似构型可能对应不同反应路线。

零温鞍点怎样对应有限温自由能瓶颈？

静态 DFT 优化给出势能面上的单个几何点。实际反应在有限温度下采样许多溶剂排布、表面振动与吸附构象，同一反应进度 s 对应一簇微观构型。由概率分布可写出 F(s) = −k_BT ln P(s) + C，自由能峰包含能量项与构型数目贡献。

水合 γ-Al₂O₃(110) 表面可用 Al_β 与晶格氧、水氧的配位数 CN_a、CN_w 描述溶出早期步骤。二维自由能面出现多组极小值与连接路径，投影曲线包含逐次水合和 Al—O—Al 断裂。溶剂配位数本身参与反应坐标，固定少数水分子的静态结构难以容纳相同的质子与配位重排。

静态能垒经过零点能、振动熵与热力学校正后，可近似得到活化自由能 ΔG^‡。谐振近似适合局部较硬的振动，低频转动、平移、柔性吸附物和液体溶剂会放大误差。温度、压力、覆盖度和标准态必须随 ΔG^‡ 一起报告，不同条件下的峰高不宜直接混排。

酸性沸石中的质子跳跃和 Al—O(H) 断裂给出两类对照：静态路径记录初态、鞍点、中间体与终态的电子能量，增强采样则沿集体变量恢复自由能曲线。图6h、i采用机器学习集体变量描述 FAU 骨架中的反应。静态峰与自由能峰回答的是不同统计条件下的问题，数值接近时也不能省略定义。

经典过渡态理论把体系穿过分界面的平衡通量换算为速率，常写为 k = κ(k_BT/h)exp[−ΔG^‡/(k_BT)]。κ 是传输系数，用来修正穿过后回返的轨迹。自由能峰只给出指数项的主要尺度；回穿、量子隧穿和非绝热跃迁会改写速率。

怎样检查反应坐标有没有漏掉慢变量？

给定构型 x，从相同位置发射许多带有热涨落的短轨迹，先到达产物盆地 B 的比例记作 committor p_B(x)。p_B 接近 0 的构型倾向返回 A，接近 1 的构型倾向进入 B；p_B ≈ 0.5 的构型组成过渡态集合。有限温描述关注一组具有相近反应倾向的构型，而非一幅固定快照。

优良反应坐标应随 p_B 单调变化，在候选分界面附近给出集中于 0.5 的 committor 分布。图7中的 HLDA 坐标与 p_B 呈清晰对应，第三集体模态在跃迁区变化最强。committor 检验直接询问“下一步先到哪一个势阱”，它比只观察坐标曲线的峰形更接近动力学定义。

候选坐标取同一数值时，若 p_B 仍从 0 到 1 大范围散开，体系还记得某个未纳入的慢变量。界面水的取向、质子线连通状态、表面重构、吸附物方位、局部电荷与自旋都可能承担这类记忆。宽而多峰的 committor 分布是坐标缺项的信号，可据此增加集体变量或拆分路径族。

静态鞍点的检查重点不同：力收敛保证驻点位置，振动曲率确认唯一不稳定方向，沿该方向正反位移或路径跟踪检验两端连接。对于周期固体，应排除整体平移、数值噪声和过软晶格模。几何连接证据限定“它连接谁”，自由能采样与 committor 再限定有限温条件下是否形成有效分界面。

计算结果中出现一个虚频、一个 NEB 峰或一条自由能峰，都只完成对应物理量的识别。可用于反应判断的过渡态还应写明电子态、原子路径、能量零点、温度环境与相邻势阱。这些条件变化后，过渡态构型、活化自由能乃至占优路线均可能随之改变。