过渡态与反应中间体辨异:势能面、振动特征、停留时间及反应路径计算

说明:本文华算科技主要介绍过渡态反应中间体势能面位置、振动特征、停留时间和反应路径计算中的差别。

能量图里,过渡态和中间体常被画成相邻的峰与谷。两者都可能对应一套原子坐标,也都能由电子结构计算得到能量,但它们描述的动力学对象完全不同。中间体对应可以容纳概率分布的势能盆地,过渡态对应相邻盆地之间的分界区域。

中间体与过渡态的基本含义

包含 3N 个原子坐标的体系可用一张高维势能面描述。几何优化寻找的是势能梯度接近零的驻点。局部极小值沿各个微小位移方向都升能,周围形成一个盆地;反应物、产物以及反应途中出现的中间体都可占据这样的极小值。

过渡态与反应中间体辨异:势能面、振动特征、停留时间及反应路径计算
图1. 势能面上的局部极小值、过渡态、最小能量路径,以及 NEB 与二聚体搜索示意。DOI:10.3389/fchem.2019.00202

过渡态对应一阶鞍点。它沿成键、断键或原子迁移的反应方向向两侧降能,沿其余独立方向则保持局部最低。梯度为零只能确认驻点,Hessian 矩阵的曲率符号才区分极小值与一阶鞍点。

频率计算把这一区别写成可检查的结果。稳定中间体在去除平动和转动后应具有全实频;合格的一阶鞍点具有唯一一个虚频,其振动箭头应沿着目标键的生成、断裂或迁移方向。两个及以上虚频通常对应高阶鞍点,不能直接充当某个基元步骤的过渡态。

过渡态与反应中间体辨异:势能面、振动特征、停留时间及反应路径计算
图2. Au(111) 表面三步反应的初态、过渡态、中间体、终态结构及相对能量。DOI:10.1038/s41467-019-10407-6

图2中的 Int1、Int2 位于各自能谷,结构可从附近扰动后回到同一极小值;TS1、TS2、TS3 位于相邻能谷之间。“有一套优化坐标”并不会赋予化学物种身份,驻点类型决定这套坐标代表可占据状态还是跨越位置。

反应能量剖面常用短横线画出过渡态能量。横线长度仅用于绘图排版,体系通过峰顶的时间不能从线长读取;中间体横线也只标示该极小值的能量,真实盆地宽度须由多维势能面或自由能面判断。

振动模和停留时间怎样区分两类状态

分子动力学轨迹进入一个势能盆地后,会在极小值附近经历大量振动,直到热涨落把体系送过出口。粗粒化以后,轨迹表现为一个个盆地之间的跳跃。中间体的寿命来自离开该盆地所需的能垒,寿命可长可短,并不要求样品能够分离保存。

过渡态与反应中间体辨异:势能面、振动特征、停留时间及反应路径计算
图3. 分子动力学轨迹在势能盆地内停留,并粗粒化为极小值之间的离散跃迁。DOI:10.3389/fchem.2019.00202

体系到达鞍区后,虚频方向上的微小位移会把轨迹带向反应物侧或产物侧。峰顶缺少恢复振动,热轨迹不会围绕它形成类似中间体的长期占据。过渡态理论把该区域处理成分隔两侧状态的面,并统计单位时间内穿过它的通量。

有限温度下,“过渡态”还可指分界面附近的一组构型。DFT 优化得到的一阶鞍点是这组构型的静态代表,零点能、振动熵、溶剂重排和表面覆盖度会改变活化自由能。单个鞍点结构与过渡态系综属于不同描述,二者也都不构成可积累的反应中间体。

过渡态与反应中间体辨异:势能面、振动特征、停留时间及反应路径计算
图4. 甲醇在 H-ZSM-5 中反应后的 13C 固体核磁谱,硝基甲烷参与时出现 56.5 ppm 的表面甲氧基信号。DOI:10.1038/s41467-025-64524-6

中间体若具有足够浓度和寿命,NMR、红外或瞬态光谱可能捕捉其特征。图4中 56.5 ppm 信号归属于 H-ZSM-5 内的表面甲氧基物种。可观测谱峰支持盆地内存在有限布居;没有谱峰仍可能源于覆盖度低、寿命短或信号重叠,不能据此删除计算路径中的极小值。

实验很少直接记录过渡态的稳态浓度,相关结构多由动力学同位素效应、速率温度关系、飞秒谱或计算约束推断。“能否分离”不是驻点分类标准,势能曲率给出几何身份,形成与消耗速率控制实验上能否捕捉。

多步反应中两类状态怎样改变动力学

一个中间体把总反应切成两个或更多基元步骤,每一步各有自己的过渡态。某一步的正向活化能以该步起始极小值为参照,不能始终从整条路径的最初反应物量起。深能谷会提高后续逃逸能垒,并可能积累覆盖位点;浅能谷仍能改变相邻步骤的速率常数。

过渡态与反应中间体辨异:势能面、振动特征、停留时间及反应路径计算
图5. 二聚与单体钒氧位点上的 NH3-SCR 循环结构、各步过渡态及相对能量。DOI:10.1038/s41467-024-47878-1

图5的字母 A 到 O 对应催化循环中的一系列吸附构型,中间标出的 TS1 至 TS4 对应跨越事件。二聚钒氧位点上 NH3 脱氢相关能垒为 1.17 eV,单体位点上 V–OOH 脱水相关能垒为 1.44 eV。过渡态控制相邻盆地间的通行速率,中间体能量与覆盖度则改变通行的起点人口。

反应能图中的最高峰未必对应稳态速率最敏感的步骤。各中间体覆盖度会随温度、分压、溶剂和位点竞争改变,表观速率取决于整套生成与消耗方程。能垒与中间体布居在微观动力学中承担不同变量,只按绝对峰高排序会丢失参照态和覆盖度信息。

过渡态与反应中间体辨异:势能面、振动特征、停留时间及反应路径计算
图6. H-ZSM-5 孔道内甲醇、表面甲氧基与二甲醚转化步骤的自由能及关键结构。DOI:10.1038/s41467-025-64524-6

表面甲氧基既可由甲醇脱水生成,也可与另一分子甲醇生成二甲醚。图6为每个步骤分别标出 ΔG,过渡态结构中的质子转移和 C–O 键重排对应跨越方向。同一个中间体能够连接多条反应支路,每条支路需由各自的鞍点决定动力学权重。

生成速率高于消耗速率时,中间体的盆地占据上升;生成慢、消耗快时,即使极小值真实存在,稳态浓度也可能很低。结构身份与稳态丰度必须分开记录,低丰度不会把极小值改成鞍点。

NEB 中间图像为什么不能命名为反应中间体

NEB 在已知初态和终态之间布置若干图像,图像编号处在端点之间,故常被软件界面称为 intermediate images。这里的 intermediate 是数值序列位置。每个图像都是反应路径的离散采样副本;即使 NEB 有效力已经收敛,普通图像沿路径方向的真实力通常仍不为零,因而没有独立驻点身份。

收敛后的普通 NEB 用有限图像逼近最小能量路径,最高能图像可能仍偏离真正鞍点;攀爬图像方法会把最高能副本推向鞍点。最高能图像需要频率与路径连接检查,相邻的普通图像更不能按文件序号命名成化学中间体。

若一条 NEB 路径在峰间出现清楚的能谷,可把该构型单独解除弹簧约束并进行几何优化与频率计算。优化后得到全实频极小值,且两侧各自存在一阶鞍点时,原来的单步路径应按两个基元步骤处理。能量折线因图像过少产生的小凹陷没有这种含义。

过渡态与反应中间体辨异:势能面、振动特征、停留时间及反应路径计算
图7. Monsanto 反应的两条竞争路径,路径在中间体 m5 处分支并经过不同能垒。DOI:10.1038/s41467-024-47997-9

反应网络把中间体记为独立节点,把连接节点的过渡事件记为边。图7的两条路径在 m5 处分叉,m5 的结构和能量参与后续两种反应分支的起点定义。真正的中间体能够承接多条进入与离开路径,NEB 普通图像只属于某一次路径离散化。

计算记录中可分别使用“第 i 个 NEB 图像”“一阶鞍点 TS”“局部极小值 Int”三类名称。名称应跟随数值对象和曲率证据,这样能量剖面、频率结果、微观动力学状态和实验谱学归属才保持同一套物理定义。

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
(0)
上一篇 3小时前
下一篇 2025年5月21日 上午10:30

相关推荐