说明:本文华算科技主要介绍过渡态单虚频的物理含义,以及如何用振动方向、反应路径和同一计算模型下的频率结果核对候选结构。


“一个虚频”具体数的是什么?
在固定组成、电荷态和原子自由度的势能面上,过渡态是梯度接近零的驻点。稳定构型沿允许的微小位移,能量在各独立方向都向上弯曲;基元反应的一阶鞍点则有一个能量向下弯曲的独立方向。沿这个方向朝正、负两侧移动,才可能进入两个相邻的能量盆。
频率计算求的是质量加权 Hessian 的本征值。Hessian 由总能对原子坐标的二阶导数构成;本征值 λi 与角频率平方 ωi2 对应。λii|;有的软件把它显示为负的 cm−1 数值。判定依据是负本征值的个数,频率前的符号只是输出约定。
这里的“一个”针对同一个驻点、同一个势能面和相同的可动原子集合。孤立分子的整体平移与转动、周期模型中的整体平移接近零频;它们与局部负曲率模属于不同对象。若结构尚未达到驻点,频率分解仍可给出数值,负模个数却不再对应严格的鞍点阶数。

碳氮双层的 α-C2N2 向 β-C2N2 转变时,爬山图像弹性带路径出现约 1.90 eV 的能量峰。原研究又在峰顶构型检查虚频振型:图中的双向箭头标出参与 C—N 键断裂、重组的原子位移。路径峰和振型方向相互对应,比单独报告“最高图像”多了一项局部曲率证据。
图上红色曲线是离散路径图像之间的插值。插值峰值本身不输出完整 Hessian,无法给出其余方向的曲率。驻点收敛、一次负本征值与对应振型分别回答“坐标是否停在驻点”“鞍点是几阶”“负曲率对应哪组原子位移”。


怎样核对虚频振型确实属于目标反应?
金红石 TiO2 表面 H2 活化包含两个基元变化:H—H 键断裂形成相邻的 H+ 与 H− 物种,随后 Ti—H 中的氢向邻近 O 位迁移。图中的 TS1 和 TS2 分别服务于这两个步骤。同一反应网络有两个过渡态,每个合格的一阶鞍点各自对应一个负曲率模。

对 TS1,振型应包含 H—H 拉长与 Ti—H、O—H 局部关系的变化;对 TS2,主要关注 H 从 Ti 邻域向 O 位的转移。虚频动画可能同时出现多个原子移动,一个本征矢可以是协同位移。计数依据是独立负本征方向的数目,并非画面中移动的原子数或发生变化的键数。
四个 TiO2 晶面给出了不同的 TS1、TS2 能量位置及键长。图中的曲线按相同反应顺序排列;峰位和能垒随晶面变化,每个晶面上的负本征值个数都应由该晶面的过渡构型计算。把 (001) 面的虚频数值套到 (100) 面,只能得到缺少局部势能面依据的标签。

若唯一负模主要表现为吸附分子整体转动、远处表面原子起伏,而目标 H—H 或 Ti—H 键长几乎不变,这个驻点对应的反应对象就与预设步骤不同。沿本征矢作正、负位移后分别观察键长和局部配位,能确定两侧最初朝哪个结构区域下降;位移方向的正负号本身没有化学优先级。
振动动画中的箭头长度还受绘图时的放大比例影响,不能拿画面的位移幅度估算能垒。对 H2/TiO2,应记录反应中心几条键长随模式坐标的变化方向:TS1 的 H—H 拉长与 TS2 的 H 跨位迁移属于不同几何变化。能垒则来自相同能量参照下初态与鞍点的能量差。


怎样核对路径两端,并理解“唯一”?
金红石 TiO2 晶格中的间隙 H 或 D 可占据两类相邻位点。沿 c 轴观察,同一局部环境里的位点位置和扩散角 γ 有明确几何含义。先确定两端是哪个稳定占位,过渡态虚频振型才有可对照的迁移方向;只看一个原子在图上的箭头,无法判定它通向哪两个结构盆。

在同一体系的爬山图像弹性带结果中,能量随 γ 升高再降低;原研究还比较了鞍点虚频模与迁移轨迹的夹角。原图给出径向和 c 轴坐标变化,说明反应坐标不一定是单一键长。路径峰顶附近应由优化后的驻点和频率结果确认,离散图像或样条曲线的最高点只限定候选区域。

从候选鞍点沿负模的正、负两侧分别作微小位移,再让构型沿势能下降,可以检查所得稳定态;分子体系常用内禀反应坐标,固体表面也可用两侧松弛和最小能量路径互证。若两侧都回到同一构型,或进入与目标步骤无关的吸附位,“单虚频”只完成了局部阶数检查。反应连接关系仍未建立。
两侧下降得到的终点还需与预先设定的初态、终态在原子映射、吸附位和化学计量上对照。对于氢扩散,周期性晶胞可能把看似不同的氢位置映射为等价占位;对于表面加氢,产物若在松弛过程中发生额外质子转移,原先选定的基元步骤已经改变。此时应重新定义端点与对应反应坐标。
“唯一虚频”表示这个候选驻点的可动坐标中只有一个负本征方向。周期表面有多个等价位点,反应物又可能选择不同晶面、吸附构型和迁移路线;每条路线可以各自拥有一阶鞍点。即使一个候选结构通过频率检查,其他位点和构型之间仍可存在各自的过渡态。


额外小虚频和模型改变怎样处理?
碳氮双层 α-C2N2、β-C2N2 与 γ-C4N4 的声子色散没有虚频,描述的是稳定晶相在布里渊区各波矢下的集体振动;图1的局域鞍点虚频来自相变路径上的另一构型。晶相声子谱与过渡态 Hessian 取样的结构、坐标和波矢范围不同,计数不能直接互换。

在大超胞或含柔性吸附物的表面模型中,平移近零模、软扭转、未充分优化的原子坐标和有限差分力误差都可能产生很小的负频率。先核对驻点残余力和额外负模的位移形状,再检验力与能量收敛、位移步长、超胞、固定原子集合。若额外模式在收敛后持续存在,它对应第二个独立下降方向,候选结构应按高阶鞍点重新处理;没有通用的 cm−1 截断值可以代替这些检验。
局部频率计算若只移动吸附物和表面近邻原子,得到的是受限子空间中的 Hessian。未参与频率计算的深层原子或晶格形变,可能具有额外的下降方向;反过来,允许整块 slab 平移又会引入接近零的整体模。报告单虚频时,频率矩阵覆盖哪些原子、哪些原子固定,应与过渡态优化的约束保持一致。
从近似势能面取得候选过渡态时,频率还依赖该势能面的曲率。NewtonNet 的研究将机器学习势优化构型带回 DFT 势能面,比较最负频率的符号,并统计重新优化后连接两端的数量。图中有部分机器学习势上的负频率在 DFT 重新求导时改变符号,不同势能面上的“单虚频”不能自动继承。

对表面反应,固定底层原子、改变电荷态或换用溶剂模型会改变允许位移与局部能量曲率。若研究目标是特定电位和覆盖度下的 H 转移,频率、振型及两侧下降路径应来自该反应环境下同一模型;真空中性表面的一个负模只能界定真空中性表面上的驻点。
