说明:本文华算科技主要介绍中间体与过渡态在 DFT 中的结构定义、振动特征、路径连接,以及计算模型改变后的归属。
DFT如何定义中间体和过渡态?
反应中暂时形成、随后继续转化的物种称为中间体。在常规静态 DFT 描述中,中间体的平衡构型对应势能面的局部极小值:原子受到小幅扰动后,邻近构型的能量升高。反应物和产物也可具有这种局部结构;“中间体”这个名称还包含它在整条反应中的位置。

DFT 按给定的核坐标求解电子结构,得到能量随原子位置变化的函数 E(R)。沿某个基元步骤连接相邻极小值时,过渡态构型通常对应一阶鞍点:沿一个独立位移方向,离开该点会降能;沿其余内部方向,小幅位移会升能。下文的“过渡态”均指这种静态构型,有限温度下的过渡态集合涉及额外的统计采样。
酸酐交换剖面中的 INT2 与 TS2a 分别标为 9.9 和 42.4 kcal mol−1。两者都高于图中的反应物零点,但 INT2 仍是中间体。局部极小值描述邻近构型的曲率,无需比所有反应物都低能。该图采用 DFT 优化结构,并以耦合簇单点能、溶剂及热化学校正得到自由能;水平短线也只是状态标记。
在驻点附近,力 F=−∇E 接近零,能量可写为 E(R+u)≈E(R)+uTHu/2,H 是能量对坐标的二阶导数矩阵,即 Hessian 矩阵。去除平移、转动等零模后,普通极小值的内部曲率为正,一阶鞍点有且只有一个负曲率本征方向。这个方向可以涉及多个原子的协同位移。

酸酐交换的过渡构型中,多个 C—O 距离分别变化,形成与断裂的键未必同步改变。“旧键断了一半、新键形成了一半”没有统一的键长标准,过渡构型与反应物或产物的几何相似程度随具体反应而变。结构图中的虚线仅突出正在变化的原子间距;外观相似的构型仍可能分别属于极小值、鞍点或尚未收敛的路径点。
力、振动和反应路径怎样确认两类构型?
局部曲率与原子位移
几何优化报告“收敛”,通常只说明规定的力或位移阈值已经满足。对质量加权 Hessian 求本征值 λ,可得到 ω2=λ。局部极小值的内部振动频率应为实数;一阶鞍点有一个 λ<0,相应频率写成虚数。程序常用负数显示虚频的模长,原子并没有以负频率做周期振动。
虚频本征矢给出具体的原子位移方向。若待研究步骤是氢转移,负模应包含供体—H 距离增大和受体—H 距离减小;若负模仅对应远端基团扭转,当前鞍点可能连接两个构象。接近零的微小虚频还可能受残余力、积分网格或有限差分误差影响,需比较加严收敛后频率的符号和位移图样。
相邻状态与路径采样点
从候选过渡态沿负模的正、反方向作微小位移,再跟踪降能路径,可以确定它连接哪些极小值。内禀反应坐标 IRC 在质量加权坐标中沿最陡下降方向延伸;两端继续优化后,应得到目标反应物、产物或相邻中间体。碳酸酐酶簇模型的示例以 TS2 连接 INT1 和 INT2,展示了这种归属方式。

IRC 横轴记录构型变化,没有赋予分子运动时间。端点的键连接、吸附位点与配位状态若与预定反应不同,即使虚频数量正确,也应把鞍点归属到实际找到的反应。对表面模型,沿负模两侧弛豫或追踪最小能量路径同样能核查连接;原子数、电荷、自旋处理和固定原子范围须保持一致。
NEB 用一串离散构型近似连接两个端点的路径,其中 intermediate image 是中间图像。多数图像仍有沿路径的真实受力,并未达到驻点。最高图像只是过渡态候选;爬山图像 NEB 可逼近鞍点,随后仍需确认收敛、负模与路径连接。曲线低点则须独立弛豫并检验曲率,才可认定为新增中间体。

三类路径示例分别出现转变区内的额外低点、化学转变前的较高峰,以及最高点与目标键变化重合。一个端点对可能包含几个基元步骤,新增中间体会要求分别寻找前后两段的鞍点。这里的图用于说明路径分类,出自机器学习势搜索研究;它没有对每个黑点逐一完成中间体或过渡态认证。
表面 slab 常固定底层原子;若振动计算仅移动吸附物和少量表面原子,所得 Hessian 只覆盖这部分坐标。受限模型里的极小值或单负模,尚未检验冻结原子的位移。释放相关表面原子后,吸附物转动、金属重排或质子迁移可能产生新的降能方向。
反应环境与原子运动改变后,原有归属还成立吗?
中间体与过渡态属于指定的计算模型
电荷态、自旋态、质子化程度或显式溶剂改变后,E(R) 本身随之改变。原模型中的中间体可能继续弛豫成另一物种,原有鞍点可能移动或消失。在旧构型上计算一次新条件下的能量,只得到新势能面上的一个数值;要沿用中间体或过渡态名称,须重新核查力、曲率和连接。
酸催化酸酐交换增加了质子化反应物,过渡构型中出现质子参与的键重排,并形成新的中间体。带电模型与中性模型各有反应物参考态,图中的高度应在各自反应定义内比较;评价酸浓度带来的速率变化还涉及质子化平衡。该案例是分子模型,聚合物链段约束和周围基团并未全部进入模型。

对某个中间体 I,通向相邻过渡态的局部活化自由能为 ΔG‡=GTS−GI。提高某个中间体的稳定程度,而过渡态能量变化较小时,离开它的能垒会增大。能量差只含反应物和产物的自由能图,描述各步热力学变化;若没有定位过渡态,连线中间的活化自由能仍未知。
相同物种之间仍可能存在多条路径
相同的吸附物初末态可以通过不同的转动、扩散和断键顺序连接。NbSe2(111) 上的氢转移示例中,改变原子映射后,两条路径包含的分子转动不同,最高点也不同。这里的曲线采用 DFT 数据训练的机器学习势;若用于定量 DFT 机理判断,需要在目标泛函下重新验证相关驻点和能垒。

编号相同的原子在整条路径中的运动应连续;转移的是哪一个 H、吸附物是否先转动、表面是否重构,会决定鞍点连接的实际构型。较低的一个鞍点只描述对应路径。若另一条路径经过额外中间体,反应网络中就需保留该物种及其进入、离开的步骤,分别计算各段能垒。
中间体在指定模型里形成局部极小值后,实验中能否积累还受生成速率和消耗速率控制。浅能阱中的物种可能迅速转化,谱学信号也可能低于检测限;持续形成而缓慢消耗的表面物种则可能占据大量位点。比较实验时,温度、电位、覆盖度和观测时间决定了应检验哪一种布居,而静态驻点计算本身没有给出这段时间内的浓度变化。
