中间体与过渡态在 DFT 里怎么区分?

说明:本文华算科技主要介绍中间体与过渡态在 DFT 中的结构定义、振动特征、路径连接,以及计算模型改变后的归属。

DFT如何定义中间体和过渡态?

反应中暂时形成、随后继续转化的物种称为中间体。在常规静态 DFT 描述中,中间体的平衡构型对应势能面的局部极小值:原子受到小幅扰动后,邻近构型的能量升高。反应物和产物也可具有这种局部结构;“中间体”这个名称还包含它在整条反应中的位置。

中间体与过渡态在 DFT 里怎么区分?
图1. 甲基丙烯酸酐与4-戊烯酸酐交换的两条自由能剖面,标出中间体及过渡态。DOI:10.1038/s42004-025-01736-3

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 矩阵。去除平移、转动等零模后,普通极小值的内部曲率为正,一阶鞍点有且只有一个负曲率本征方向。这个方向可以涉及多个原子的协同位移。

中间体与过渡态在 DFT 里怎么区分?
图2. 两条酸酐交换路径中的六个过渡态构型及反应中心的原子间距。DOI:10.1038/s42004-025-01736-3

酸酐交换的过渡构型中,多个 C—O 距离分别变化,形成与断裂的键未必同步改变。“旧键断了一半、新键形成了一半”没有统一的键长标准,过渡构型与反应物或产物的几何相似程度随具体反应而变。结构图中的虚线仅突出正在变化的原子间距;外观相似的构型仍可能分别属于极小值、鞍点或尚未收敛的路径点。

力、振动和反应路径怎样确认两类构型?

局部曲率与原子位移

几何优化报告“收敛”,通常只说明规定的力或位移阈值已经满足。对质量加权 Hessian 求本征值 λ,可得到 ω2=λ。局部极小值的内部振动频率应为实数;一阶鞍点有一个 λ<0,相应频率写成虚数。程序常用负数显示虚频的模长,原子并没有以负频率做周期振动。

虚频本征矢给出具体的原子位移方向。若待研究步骤是氢转移,负模应包含供体—H 距离增大和受体—H 距离减小;若负模仅对应远端基团扭转,当前鞍点可能连接两个构象。接近零的微小虚频还可能受残余力、积分网格或有限差分误差影响,需比较加严收敛后频率的符号和位移图样。

相邻状态与路径采样点

从候选过渡态沿负模的正、反方向作微小位移,再跟踪降能路径,可以确定它连接哪些极小值。内禀反应坐标 IRC 在质量加权坐标中沿最陡下降方向延伸;两端继续优化后,应得到目标反应物、产物或相邻中间体。碳酸酐酶簇模型的示例以 TS2 连接 INT1 和 INT2,展示了这种归属方式。

中间体与过渡态在 DFT 里怎么区分?
图3. 碳酸酐酶簇模型中连接 INT1、TS2 与 INT2 的 IRC 能量剖面及构型。DOI:10.1038/s41598-025-32638-y

IRC 横轴记录构型变化,没有赋予分子运动时间。端点的键连接、吸附位点与配位状态若与预定反应不同,即使虚频数量正确,也应把鞍点归属到实际找到的反应。对表面模型,沿负模两侧弛豫或追踪最小能量路径同样能核查连接;原子数、电荷、自旋处理和固定原子范围须保持一致。

NEB 用一串离散构型近似连接两个端点的路径,其中 intermediate image 是中间图像。多数图像仍有沿路径的真实受力,并未达到驻点。最高图像只是过渡态候选;爬山图像 NEB 可逼近鞍点,随后仍需确认收敛、负模与路径连接。曲线低点则须独立弛豫并检验曲率,才可认定为新增中间体。

中间体与过渡态在 DFT 里怎么区分?
图4. 三类 NEB 路径的能量与图像序号;背景色区分初态类、化学转变区和终态类构型。DOI:10.1038/s41524-026-02298-1

三类路径示例分别出现转变区内的额外低点、化学转变前的较高峰,以及最高点与目标键变化重合。一个端点对可能包含几个基元步骤,新增中间体会要求分别寻找前后两段的鞍点。这里的图用于说明路径分类,出自机器学习势搜索研究;它没有对每个黑点逐一完成中间体或过渡态认证。

表面 slab 常固定底层原子;若振动计算仅移动吸附物和少量表面原子,所得 Hessian 只覆盖这部分坐标。受限模型里的极小值或单负模,尚未检验冻结原子的位移。释放相关表面原子后,吸附物转动、金属重排或质子迁移可能产生新的降能方向。

反应环境与原子运动改变后,原有归属还成立吗?

中间体与过渡态属于指定的计算模型

电荷态、自旋态、质子化程度或显式溶剂改变后,E(R) 本身随之改变。原模型中的中间体可能继续弛豫成另一物种,原有鞍点可能移动或消失。在旧构型上计算一次新条件下的能量,只得到新势能面上的一个数值;要沿用中间体或过渡态名称,须重新核查力、曲率和连接。

酸催化酸酐交换增加了质子化反应物,过渡构型中出现质子参与的键重排,并形成新的中间体。带电模型与中性模型各有反应物参考态,图中的高度应在各自反应定义内比较;评价酸浓度带来的速率变化还涉及质子化平衡。该案例是分子模型,聚合物链段约束和周围基团并未全部进入模型。

中间体与过渡态在 DFT 里怎么区分?
图5. 酸催化酸酐交换的反应结构、过渡构型及自由能剖面。DOI:10.1038/s42004-025-01736-3

对某个中间体 I,通向相邻过渡态的局部活化自由能为 ΔG‡=GTS−GI。提高某个中间体的稳定程度,而过渡态能量变化较小时,离开它的能垒会增大。能量差只含反应物和产物的自由能图,描述各步热力学变化;若没有定位过渡态,连线中间的活化自由能仍未知。

相同物种之间仍可能存在多条路径

相同的吸附物初末态可以通过不同的转动、扩散和断键顺序连接。NbSe2(111) 上的氢转移示例中,改变原子映射后,两条路径包含的分子转动不同,最高点也不同。这里的曲线采用 DFT 数据训练的机器学习势;若用于定量 DFT 机理判断,需要在目标泛函下重新验证相关驻点和能垒。

中间体与过渡态在 DFT 里怎么区分?
图6. 原子映射改变后的过渡态能量差分布,以及 NbSe2(111) 上氢转移的路径与构型。DOI:10.1038/s41524-026-02298-1

编号相同的原子在整条路径中的运动应连续;转移的是哪一个 H、吸附物是否先转动、表面是否重构,会决定鞍点连接的实际构型。较低的一个鞍点只描述对应路径。若另一条路径经过额外中间体,反应网络中就需保留该物种及其进入、离开的步骤,分别计算各段能垒。

中间体在指定模型里形成局部极小值后,实验中能否积累还受生成速率和消耗速率控制。浅能阱中的物种可能迅速转化,谱学信号也可能低于检测限;持续形成而缓慢消耗的表面物种则可能占据大量位点。比较实验时,温度、电位、覆盖度和观测时间决定了应检验哪一种布居,而静态驻点计算本身没有给出这段时间内的浓度变化。

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
赞 (0)
上一篇 2小时前
下一篇 2小时前

相关推荐