来源论文: https://arxiv.org/abs/2606.02917v1 生成时间: Jun 03, 2026 13:21

0. 执行摘要

过渡金属配合物的激发态动力学是现代无机光化学、太阳能转换及光电材料开发的核心基石。其中，三价铬（Cr(III)）配合物因其独特的 $3 ext{d}^3$ 电子排布，展现出极富研究价值的超快系间窜越（Intersystem Crossing, ISC）行为。实验观测表明，Cr(III) 配合物在光激发至 $^4 ext{T}_{2 ext{g}}$ 激发态后，能够在百飞秒（$< 100 ext{ fs}$）内以接近统一的量子产率窜越到亚稳态的 $^2 ext{E}_{ ext{g}}$ 态。然而，由于该自旋翻转过程极快且伴随着强烈的电子-振动双重非绝热耦合，如何在分子层面上精细描绘振动模式如何调制电子能级、驱动这一非辐射弛豫过程，一直是理论与计算化学领域的难点。

本研究采用 Born-Oppenheimer 算符起步分子动力学（BO-AIMD）模拟，结合高精度的多参考状态平均完全活性空间自洽场（SA-CASSCF）方法以及多参考二阶微扰理论（XDW-CASPT2），对典型的 Cr(III) 配合物——三(乙酰丙酮)合铬(III)（$\text{Cr(acac)}_3$）及其简化模型三(1,3-丙二酮)合铬(III)（$\text{Cr(PDO)}_3$）的激发态动力学进行了系统探究。结果表明，配合物在 $^4 ext{T}_{2 ext{g}}$ 激发态上表现出持久的、非局域化的 Jahn-Teller（JT）畸变；通过傅里叶变换动力学轨迹，识别出频率分别为 $219 ext{ cm}^{-1}$（剪切/扭转模式）和 $465 ext{ cm}^{-1}$（金属-配体对称伸缩模式）的两个关键振动。其中，低频的剪切/扭转模式能够强烈调制 $^4 ext{T}_{2 ext{g}}$ 与 $^2 ext{E}_{ ext{g}}$ 之间的瞬时能量差，使其频繁降至零（即达到势能面交叉区域），而自旋-轨道耦合（SOC）矩阵元在振动过程中保持在 $60\sim 80 ext{ cm}^{-1}$ 的高位且基本恒定。这一发现确立了“能隙调制”而非“自旋-轨道耦合调制”作为驱动超快 ISC 的主导物理机制。此外，研究还通过计算从 $^2 ext{E}_{ ext{g}}$ 态到高能配体到金属电荷转移态（$^2\text{LMCT}$）的瞬态吸收光谱，证实了实验中观测到的飞秒相干振动信号起源于 $^2 ext{E}_{ ext{g}}$ 势能面上的相干核运动。本工作为理性设计具有特定激发态寿命与发光性能的 $3 ext{d}$ 过渡金属配合物提供了稳固的理论微观机理支撑。

1. 核心科学问题、理论基础、技术难点与方法细节

1.1 核心科学问题

过渡金属配合物的自旋态转换（如系间窜越，ISC）通常被认为是由自旋-轨道耦合（SOC）作用驱动的自旋禁阻过程。在轻过渡金属（如 $3 ext{d}$ 过渡金属）中，由于相对论效应较弱，其原子自旋-轨道耦合常数显著小于 $4 ext{d}$ 和 $5 ext{d}$ 金属。然而，$\text{Cr(acac)}_3$ 等 $3 ext{d}^3$ 体系却能在 $100 ext{ fs}$ 内完成 $^4 ext{T}_{2 ext{g}} \rightarrow ^2 ext{E}_{ ext{g}}$ 的超快自旋转化。这就提出了一个根本性的科学问题：在弱自旋-轨道耦合作用下，究竟是什么物理机制消除了自旋禁阻的限制，并在飞秒时间尺度上驱动了如此高效的系间窜越？

具体而言，研究聚焦于以下子问题：

1. 激发态分子在向目标自旋态过渡时，经历了怎样的核空间演化？Jahn-Teller 畸变在其中扮演了什么角色？
2. 分子的哪些特定分子振动模式（Vibrational Modes）与该电子能态跃迁存在强耦合？
3. 振动运动是通过周期性地改变电能隙（Energy Gap Modulation），还是通过周期性增强自旋-轨道耦合常数（SOC Modulation）来加快 ISC 速率？

1.2 理论基础

解析该问题的理论大厦建立在非绝热动力学和多参考电子结构理论的基础之上。

1.2.1 Jahn-Teller 效应与多态势能面

$\text{Cr(III)}$ 的基态为 $^4 ext{A}_{2 ext{g}}$（$O_{\text{h}}$ 对称性下电子排布为 $t_{2 ext{g}}^3$）。光激发（$d\rightarrow d$ 跃迁）将一个电子由 $t_{2 ext{g}}$ 轨道提升至 $e_{\text{g}}$ 轨道，形成 $^4 ext{T}_{2 ext{g}}$ 激发态（配置为 $t_{2 ext{g}}^2 e_{\text{g}}^1$）。由于 $^4 ext{T}_{2 ext{g}}$ 是一个轨道三重简并态，根据 Jahn-Teller 定理，非共线性分子在简并电子态下必然发生自发对称性破缺畸变，以降低几何对称性并消除简并度。这就导致 $^4 ext{T}_{2 ext{g}}$ 势能面（PES）呈现出极其复杂的、具有多个极小值点和势垒的“流沙状”拓扑结构。相比之下，目标状态 $^2 ext{E}_{ ext{g}}$（排布同样为 $t_{2 ext{g}}^3$，但自旋翻转）由于属于壳层内配对转变，其平衡几何结构与基态 $^4 ext{A}_{2 ext{g}}$ 非常相似，受到的 Jahn-Teller 畸变较小。因此，$^4 ext{T}_{2 ext{g}}$ 和 $^2 ext{E}_{ ext{g}}$ 势能面是非平行的，它们在核配位空间的某处必然存在交叉面（Crossing Seams），这为非绝热跃迁提供了天然通道。

1.2.2 费米黄金定则与自旋-轨道核振动耦合

在非绝热理论框架下，初态 $|\text{I}\rangle$ 到末态 $|\text{F}\rangle$ 的系间窜越速率 $k_{\text{ISC}}$ 在弱耦合极限下可近似由费米黄金定则（Fermi’s Golden Rule）描述：

其中 $\hat{H}_{\text{SOC}}$ 是自旋-轨道耦合算符，$\rho(E)$ 是 Franck-Condon 加权态密度。核运动（振动）对该公式中两项都存在潜在的调制作用：

1. 自旋-轨道振动耦合（Spin-Vibronic Coupling）：核运动通过改变化学键长度或角度，直接调制自旋-轨道耦合矩阵元 $\langle \Psi_{\text{I}}(q) | \hat{H}_{\text{SOC}} | \Psi_{\text{F}}(q) \rangle$ 的大小。
2. 能隙调制（Energy Gap Modulation）：核振动引起初末态势能面瞬时能隙 $\Delta E(q) = E_{\text{I}}(q) - E_{\text{F}}(q)$ 的剧烈振荡。当能隙减小至零时， Franck-Condon 重叠积分极大化，且非绝热效应暴增，促进跃迁。

1.3 技术难点

尽管物理图像清晰，但要在定量层面模拟这一过程面临着巨大的技术瓶颈：

1. 强多参考特征（Multireference Character）：过渡金属的 $3 ext{d}$ 轨道能量非常接近，电子相关作用（特别是静态相关）极强。普通的单参考方法（如压倒性应用的密度泛函理论 TD-DFT）无法可靠描述简并或近简并的 $d\text{-}d$ 激发态，甚至可能给出错误的势能面拓扑。因此，必须使用多参考方法（如 CASSCF）。
2. 动力学模拟的计算开销（Computational Bottleneck）：在 Born-Oppenheimer 分子动力学模拟中，需要在每一个时间步（通常 $< 0.5 ext{ fs}$）实时（On-the-fly）计算电子结构梯度。多参考方法（如 CASSCF）的计算开销极大，直接用于包含数十个原子的动力学模拟在过去是不可想象的。
3. 高阶激发态与电荷转移态的精确计算：实验探测 ISC 常用的是瞬态吸收（TA）光谱，其探测信号源于激发态吸收（ESA），通常涉及配体到金属电荷转移态（$^2\text{LMCT}$）等高能状态。这类电荷转移态包含极强的动态电子相关，必须引入昂贵的多参考二阶微扰理论（CASPT2）才能定性正确地进行描述。

1.4 方法细节

为攻克上述瓶颈，本研究设计了一套多尺度、高精度的量子化学联合计算协议，其核心工作流和参数细节如下：

+--------------------------------------------------------+
| 建立简化模型 Cr(PDO)3 并利用 Molpro 优化多状态几何结构   |
| - 基/激发四重态: SA-4-CAS(3,5)SCF/LANL2DZ             |
| - 双重态: SA-2-CAS(3,5)SCF/LANL2DZ                     |
+--------------------------------------------------------+
                           | 
                           v
+--------------------------------------------------------+
| 初始状态抽样 (Wigner Sampling)                           |
| - 基于 SA-4-CAS(3,5) 基态谐振振动频率，抽取 50 组初始坐标 |
+--------------------------------------------------------+
                           | 
                           v
+--------------------------------------------------------+
| 激发态 Born-Oppenheimer 分子动力学 (BO-AIMD)            |
| - 平台: PySpawn (动力学框架) + TeraChem (GPU加速 CASSCF)|
| - 电子结构: SA-4-CAS(3,5)SCF/LANL2DZ                    |
| - 时间步长: 15 a.u. (~0.36 fs), 总步长: 25200 a.u. (~610 fs)|
+--------------------------------------------------------+
                           | 
                           v
+--------------------------------------------------------+
| 动力学轨迹后处理与物理量分析                            |
| - 键长动力学分解 (对称与非对称伸缩模式)                    |
| - 对称伸缩轨迹的傅里叶变换 (FT) 获取振动频率             |
| - 沿轨迹的 4T2g / 2Eg 瞬时能隙与 SOC 矩阵元定量跟踪        |
+--------------------------------------------------------+
                           | 
                           v
+--------------------------------------------------------+
| 瞬态激发态吸收 (ESA) 性质精细表征                       |
| - 软件: OpenMolcas                                     |
| - 方法: XDW-CASPT2/def2-SVP (包含 X2C 标量相对论校正)     |
| - 活性空间: CAS(9,8) (包含 Cr 3d 轨道及 3 个配体 pi 轨道) |
| - 状态平均: SA-41, 虚轨道移动 (Imaginary Shift): 0.25 a.u.|
+--------------------------------------------------------+

1.4.1 模型简化策略

实际实验体系 $\text{Cr(acac)}_3$ 含有 43 个原子和 183 个电子（图 2 所示）。为了在保持核心物理特性的同时降低动力学计算成本，研究者用氢原子（$-\text{H}$）替换了乙酰丙酮配体上的甲基（$-\text{CH}_3$），构建了简化模型 $\text{Cr(PDO)}_3$（三(1,3-丙二酮)合铬(III)）。这一简化将体系规模骤降至 25 个原子和 135 个电子，理论上不仅完整保留了铬原子的第一配位层（$\text{CrO}_6$ 八面体），还完好保留了配体的 $\pi$ 共轭电子结构，确保了配体场劈裂能和激发态性质不受本质影响。

1.4.2 几何优化细节

在进行动力学前，使用 Molpro 软件在 $\text{SA-CAS(3,5)SCF/LANL2DZ}$ 理论水平下对两种配合物的各目标态进行了高精度优化：

• 四重态几何优化：采用状态平均 SA-4-CAS(3,5)SCF（对四重基态 $^4 ext{A}_{2 ext{g}}$ 和三个三重简并的四重激发态 $^4 ext{T}_{2 ext{g}}$ 进行等权重平均），该设置能完美平衡基态与各激发分量之间的能量，避免状态翻转导致的优化困难。
• 双重态几何优化：采用 SA-2-CAS(3,5)SCF（等权重平衡两个简并的双重态 $^2 ext{E}_{ ext{g}}$）。为了确保轨道特征的一致性，使用优化好的四重态轨道作为双重态优化的初始猜测。
• 基组与赝势：选用 LANL2DZ 基组。其中对 Cr 原子使用有效核心电荷赝势（ECP），用以替代内部的 10 个核电子，这不仅能够大幅提升计算效率，还能部分消除标量相对论效应产生的误差。

1.4.3 BO-AIMD 模拟细节

激发态动力学模拟基于 PySpawn 动力学驱动程序展开，其底层的电子结构计算调用了以 GPU 加速著称的 TeraChem 软件包。具体参数设定：

• 电子动力学设置：轨迹被约束在最下层的四重激发态 $^4 ext{T}_{2 ext{g}}$ 势能面上单态绝热演化，电子结构方法为 $\text{SA-4-CAS(3,5)SCF}$。
• 初始状态制备：在 $^4 ext{A}_{2 ext{g}}$ 基态极小值结构下，通过简正模式分析，利用 Wigner 分布（Wigner Distribution）进行量子统计抽样，生成 50 组具有量子零点振动动能和坐标偏差的独立初始条件（包括核位置与动量）。之后，分子被瞬时垂直投影（Frank-Condon 激发）至 $^4 ext{T}_{2 ext{g}}$ 激发态，模拟超短激光脉冲的泵浦过程。
• 积分算法与时间分辨率：采用 Velocity Verlet 算法进行经典核动力学积分。时间步长设为 $\Delta t = 15\text{ a.u. } (\approx 0.36\text{ fs})$。每条轨迹传播 $25200\text{ a.u. } (\approx 610\text{ fs})$，总计运行 1700 步。50 条轨迹的统计系综能够有效消除局域热涨落，真实反映核波包的宏观相干演化。

1.4.4 XDW-CASPT2 与高能电荷转移态的确定

由于 $^2\text{E}_{ ext{g}} \rightarrow ^2\text{LMCT}$ 的转换包含强烈的配体-金属电荷转移特征，其描述极度依赖动力学电子相关。研究采用了扩展动力学加权二阶微扰理论（XDW-CASPT2，基于 OpenMolcas 实现）：

• 活性空间扩展：将活性空间从原先仅包含 Cr d 轨道的 $\text{CAS}(3,5)$ 扩大至 $\text{CAS}(9,8)$，即将配体上能量最高的 3 个占有 $\pi$ 轨道也纳入活性空间。这允许 9 个活性电子在 8 个活性轨道中自由排布，从而能够定性正确地描述配体 $\pi \rightarrow$ 铬 $t_{2 ext{g}}$ 的电荷跃迁特征。
• 状态平均与技术参数：轨道优化通过对 SA-41（41 个双重态）进行等权重平均自洽计算，确保所有潜在的 LMCT 态都被包含在内。使用 def2-SVP 基组，并通过精确双组分（X2C）相对论修正方法处理标量相对论效应。在微扰阶段，将电离势-电子亲和势（IPEA）移动设为 0 Hartree，并应用 0.25 a.u. 的虚轨道移动（Imaginary Shift），以彻底消除入侵态（Intruder States）导致的微扰能级发散问题。

2. 关键 Benchmark 体系、计算所得数据与性能分析

2.1 简化模型的合理性验证（Benchmark）

首先，必须证实简化模型 $\text{Cr(PDO)}_3$ 确实可以代表实验真实的 $\text{Cr(acac)}_3$。研究者对比了两者在 $\text{SA-CAS(3,5)SCF}$ 水平下的平衡态几何结构（表 1 与表 2）。

表 1: $\text{Cr(PDO)}_3$ 在各状态下的平衡键长（单位：Å）

表 2: $\text{Cr(acac)}_3$ 在各状态下的平衡键长（单位：Å）

数据分析要点：

1. 极高的晶格相似度：在 $^4\text{A}_{2\text{g}}$ 基态下，无论是 $\text{Cr(PDO)}_3$ 还是 $\text{Cr(acac)}_3$，其优化的 6 个 $\text{Cr-O}$ 键长均在 $1.98\text{ Å}$ 左右，高度简并，与晶体学实验测得的无偏差平均值（约为 $1.94\sim 1.96\text{ Å}$）偏差仅为 $0.02\text{ Å}$。这验证了 LANL2DZ 赝势配比的精确性。
2. Jahn-Teller 畸变的一致性：在 $^4\text{T}_{2 ext{g}}$ 状态下，两者皆表现出经典的 “伸长型”Jahn-Teller 畸变。具体而言，6 个 $\text{Cr-O}$ 键分裂为两短（轴向压缩，$\approx 1.95\text{ Å}$）和四长（赤道面拉伸，最高达 $2.17\text{ Å}$）。这证实了甲基的去除确实没有破坏激发态的配体场分裂特性和轨道力学效应，模型简化完全成立。
3. 双重态的稳定性：在 $^2\text{E}_{ ext{g}}$ 状态下，键长依然保持高度对称（约 $1.98\text{ Å}$），印证了该自旋转换仅涉及排布内部自旋配对，对结构形变贡献较小。

2.2 激发态动力学核心数据分析

2.2.1 瞬态几何变化与 Jahn-Teller 畸变动力学（图 3，图 4）

从 $^4\text{A}_{2 ext{g}}$ 几何垂直投影至 $^4\text{T}_{2 ext{g}}$ 面后：

• 键长相干振荡：系统中的六个 $\text{Cr-O}$ 键立即展现出剧烈的相干拉伸运动，其键长振幅在 $1.9\text{ Å}$ 到 $2.3\text{ Å}$ 之间往复（图 3 所示）。
• 持久的 Jahn-Teller 特征：在 610 fs 的整个动力学演化中，系统始终保持着“四长两短”的 JT 畸变分布，并未退化回高度对称的八面体。为了进一步量化该特性，研究者引入了一个特征概率：“两个最短的 $\text{Cr-O}$ 键处于空间对位（Opposite）的瞬时概率”（图 4 所示）。对于完美、无畸变的对称八面体，此概率应为 $1/5 = 0.2$。动力学演化表明，该概率在 $t=0$ 时处于 0.2 附近，但在 $50\text{ fs}$ 内迅速攀升至 $0.8$ 以上，随后在 $0.4\sim 0.6$ 之间周期性震荡，始终显著高于基准线 0.2。这直接从系综动力学层面证实，Jahn-Teller 畸变在激发态波包演化中是持久存在并处于动态相干振荡中的。

2.2.2 键伸缩振动模式的傅里叶分解（图 5，图 6，图 7）

为了阐明哪些特定核运动参与了相干演化，研究者将这 6 个 $\text{Cr-O}$ 键长动力学分解为全对称运动（Symmetric）与非对称运动（Asymmetric）分量：

• 相干全对称拉伸（图 5）：将每一步中 6 个键长的平均值在 50 条轨迹上进行系综平均。该曲线展示出极其规则、阻尼极小的飞秒尺度相干振荡，波动范围在 $2.00\text{ Å}$ 到 $2.10\text{ Å}$ 之间。其标准差 $\sigma$ 仅为 $0.02\text{ Å}$，表明不同轨迹之间的发散极小，相干性极好。
• 非相干非对称拉伸（图 6）：通过计算轨迹平均的最长与最短 $\text{Cr-O}$ 键长之差 $\Delta(\text{Cr-O})$ 来衡量。该曲线呈现无序的快速涨落，且轨迹间标准差高达 $0.11\text{ Å}$（是对称模式的 5 倍以上），说明非对称畸变在系综平均中迅速去相干并相互抵消。
• 傅里叶变换（FT）功率谱（图 7）：对相干全对称拉伸轨迹进行傅里叶变换，得到了具有两个尖锐特征峰的振动谱，其频率对应为：
  • $219\text{ cm}^{-1}$（对应周期 $152\text{ fs}$）
  • $465\text{ cm}^{-1}$（对应周期 $72\text{ fs}$）

这一计算谱与超快瞬态吸收实验测得的相干振动频率（$193\text{ cm}^{-1}$ 和 $461\text{ cm}^{-1}$，见文献17, 18）具有极高的吻合度（偏差仅在 $\approx 20\text{ cm}^{-1}$ 以内），成功赋予了实验信号明确的结构动力学归属。

2.2.3 简正模式特征分析与能级调制机制（图 8，图 9，图 10）

为了在物理上直观呈现这两个主导振动，研究者通过基态谐振分析，手动定位了两个高度对应的对称简正模式（图 8 所示）：

1. $225\text{ cm}^{-1}$ 模式（对应动力学中 $219\text{ cm}^{-1}$ 峰）：表现为 $\text{CrO}_6$ 八面体中氧原子的非径向切向运动，伴随着乙酰丙酮配体环的周期性扭转/剪切（Twisting/Scissoring）与协同整体呼吸。这一运动极易通过改变受阻配位排斥来影响配体场分裂大小。
2. $487\text{ cm}^{-1}$ 模式（对应动力学中 $465\text{ cm}^{-1}$ 峰）：表现为经典的纯径向共线 $\text{Cr-O}$ 对称伸缩（Collinear Stretching），配体骨架相对静止。

图 9 & 图 10: 势能面瞬态演化与能隙周期性关闭

通过沿动力学轨迹实时计算 $^4 ext{T}_{2 ext{g}}$ 与 $^2 ext{E}_{ ext{g}}$ 的状态能量（图 9），发现：

• $^4 ext{T}_{2 ext{g}}$ 的能量在动力学演化中表现平缓（处于势能面谷底附近）。
• 双重态 $^2 ext{E}_{ ext{g}}$（包含近简并的 $^2\text{E}_{\text{g}0}$ 和 $^2\text{E}_{\text{g}1}$）的能量则表现出极强的简谐波动，其震荡周期约为 $154\text{ fs}$。这完美对应于上述 $225\text{ cm}^{-1}$ 配体扭转/剪切模式的振动周期。
• 能隙频繁归零（图 10）：在单条轨迹层面分析 $E(^4 ext{T}_{2 ext{g}}) - E(^2\text{E}_{\text{g}0})$ 的能隙差 $\Delta E$。在绝大多数轨迹中，该瞬时能隙在振动过程中频繁跨越零点（$\Delta E = 0\text{ eV}$）。这确立了非绝热表面交叉（Conical Intersection-like Seamless Crossing）在空间中分布极广，核波包在扭转振动驱动下，每一个振动周期都会多次穿过自旋态转换的“窗口区”。

2.3 自旋-轨道耦合（SOC）与能隙调制（图 11）

为了最终确认动力学转换的控制性物理量，研究者沿动力学轨迹计算了 $^4 ext{T}_{2 ext{g}}$ 与 $^2 ext{E}_{ ext{g}}$ 的瞬时自旋-轨道耦合矩阵元 $H_{\text{SOC}}$（图 11 所示）。

• 对于四重态磁量子数 $M_{\text{S}} = 1.5$ 组分，平均 $H_{\text{SOC}}$ 保持在 $\approx 60\text{ cm}^{-1}$。
• 对于 $M_{\text{S}} = 0.5$ 组分，平均 $H_{\text{SOC}}$ 保持在 $\approx 78\text{ cm}^{-1}$。
• 高度的静态特征：在整个 610 fs 的剧烈化学键拉伸与扭转动力学模拟中，该值波动极小，没有表现出任何可观测的周期性振荡趋势。这一发现具有重大的理论科学价值：它彻底排除了“自旋-轨道振动耦合（自旋-轨道算符被振动强烈调制）”作为加速 ISC 的原因；相反，它雄辩地证明，扭转振动模式是通过周期性地“关闭和开启” $^4 ext{T}_{2 ext{g}}$ / $^2 ext{E}_{ ext{g}}$ 势能面能隙、使 Franck-Condon 因子阶段性最大化来驱动超快 ISC 的。

2.4 ESA 瞬态吸收计算与实验验证（表 3）

为了进一步证明实验 TA 光谱中观测到的飞秒相干振动动力学特征确实发生在 $^2 ext{E}_{ ext{g}}$ 势能面上，研究者在两个关键地标几何点利用高级的 XDW-CASPT2 方法计算了 $^2 ext{E}_{ ext{g}} \rightarrow ^2\text{LMCT}$ 的垂直瞬态激发态吸收（ESA）能级（表 3）。

表 3: $^2\text{E}_{ ext{g}} \rightarrow ^2\text{LMCT}$ 垂直激发能的几何相关性

理论分析： 计算得到的吸收能在 $1.67\sim 1.83\text{ eV}$ 之间。这与实验测定瞬态相干信号时使用的探测光子能量 $2.09\text{ eV}$（波长 $592\text{ nm}$） 处于高度合理的物理邻近区。由于 $^2 ext{E}_{ ext{g}}$ 面的极小值位置与 $^4\text{A}_{2 ext{g}}$ 极为贴近，而核波包在由 $^4 ext{T}_{2 ext{g}}$ 窜越至 $^2 ext{E}_{ ext{g}}$ 时正处于拉伸和扭转畸变极大的位置，这意味着跃迁后的波包将在 $^2 ext{E}_{ ext{g}}$ 面上顺着扭转和伸伸坐标展开大幅度的拉伸振荡，其瞬时跃迁能量差将随核波包运动在 $1.67\text{ eV}$ 到 $> 1.83\text{ eV}$ 之间进行飞秒尺度的周期性摆动。这完美阐明了实验探测到的强瞬态透射率震荡正是起源于核波包在目标自旋面上该相干核振动的延续。

3. 代码实现细节、复现指南与开源工具

为了使量子化学研究人员能够复现本研究中报道的激发态 Born-Oppenheimer 分子动力学模拟及后续能谱分析，以下给出了具体而微的计算环境部署与实现指南。

3.1 软件栈要求及开源仓库

1. PySpawn (2.0.0+)：基于 Python 的多状态非绝热动力学核心框架，负责核经典方程的积分（Velocity Verlet）、初态 Wigner 抽样以及调用外部电子结构软件计算力常数和能量。
   • 开源地址：GitHub - mccoygroup/pyspawn 或 Levine 课题组分支。
2. TeraChem (1.9+)：商业 GPU 加速电子结构计算程序，用于在分子动力学模拟的每个时间步高效计算 SA-CASSCF 的梯度。由于其利用了 GPU 上的多中心双电子积分并行算法，CASSCF 计算速度比传统 CPU 软件快 1-2 个数量级。
   • 官方主页：TeraChem - GPU-accelerated Quantum Chemistry
3. OpenMolcas (v21.11+)：开源的多参考量子化学程序包，用于动力学之后，计算高精度的 XDW-CASPT2 能级和瞬态偶极矩。
   • 开源地址：GitHub - OpenMolcas/openmolcas
4. Molpro (2020+)：用于初始基态和激发态的高精度几何优化及谐振频率分析。
   • 官方主页：Molpro Home

3.2 动力学模拟复现关键步骤与输入流示例

步骤 1: 几何优化与频率分析 (Molpro)

在 Molpro 中构建 $\text{Cr(PDO)}_3$ 几何。使用 SA-4-CAS(3,5) 优化 $^4\text{T}_{2 ext{g}}$ 的输入核心代码段如下：

***,Cr_PDO3_Opt
memory,512,m
geometry={
  Cr   0.00000   0.00000   0.00000
  O    1.38500   1.38500   0.00000
  ! ... 其余 23 个原子配位坐标 ...
}
basis=lanl2dz

{rhf;                              ! 首先进行单参考Hartree-Fock计算作为猜测}
{multi;
  occ,15,10;                       ! 定义活性轨道占据}
  closed,12,8;                     ! 冻结内层轨道}
  active,3,2;                      ! 建立活性空间 (3,5)}
  state,4;                         ! 平均 4 个四重态}
  weight,0.25,0.25,0.25,0.25;      ! 等权重平均}
}
optg;                              ! 几何优化命令
frequencies;                       ! 执行频率分析获取 Hessian 用于 Wigner 抽样
--- 

步骤 2: Wigner 初始条件抽样 (PySpawn)

运行 PySpawn 提供的 Wigner 抽样脚本，读取 Molpro 频率分析输出文件（通常为 .out 或 .xml），生成包含 50 个轨迹初始核坐标与核动量数据的基础文件 wigner_samples.dat：

python -m pyspawn.wigner -f cr_pdo3_freq.out -n 50 -o wigner_samples.dat --temp 0

步骤 3: 编写 PySpawn-TeraChem 动力学控制脚本

创建一个名为 run_aimd.py 的动力学控制文件，定义与 TeraChem 的数据交互和时间推进参数。以下为核心配置骨架：

import pyspawn
import pyspawn.aimd
import sys

# 初始化动力学模拟器
sim = pyspawn.aimd.simulation()

# 读取 Wigner 抽样点数据
sim.load_wigner_samples("wigner_samples.dat")

# 配置计算后端：TeraChem 接口
sim.set_quantum_backend("terachem")
sim.set_tc_binary("/usr/local/bin/terachem")

# 设置 TeraChem 运行参数模板
tc_template = {
    "method": "casscf",
    "casscf_active_electrons": 3,
    "casscf_active_orbitals": 5,
    "casscf_states": 4,
    "casscf_state_weights": "0.25 0.25 0.25 0.25",
    "basis": "lanl2dz",
    "purify": "no",
    "precision": "mixed",  # 开启混合精度混合计算，加速GPU
    "gpudevice": "0"       # 指定运行的 GPU 号
}
sim.set_tc_options(tc_template)

# 设置动力学核心参数
sim.set_timestep(15.0)       # 步长 15 a.u. (~0.36 fs)
sim.set_max_steps(1700)      # 最大演化步数 (~610 fs)
sim.set_propagator("verlet")
sim.set_target_state(1)      # 目标激发态：4T2g (第二最低四重态分量)

# 并行或串行分发 50 条独立轨迹运行
sim.run_trajectories()

步骤 4: 沿轨迹进行 XDW-CASPT2 分析 (OpenMolcas)

动力学运行完成后，从 PySpawn 生成的轨迹文件（如 traj_0.xyz）中提取特定步数的核几何结构，转换为 OpenMolcas 格式，并调用以下 .input 输入脚本来精细确定能隙与激发性质：

&GATEWAY
  Title = Cr(PDO)3 XDW-CASPT2 calculation
  Coord = frame_400.xyz
  Basis = def2-SVP
  Relativistic = X2C

&SEWARD

&RASSCF
  Title = SA-41-CAS(9,8) active space run
  Spin = 2
  nActEl = 9 0 0
  Active = 8
  Inactive = 28
  CIsd
  SexE = 1
  MultiState = 41 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41

&CASPT2
  Title = XDW-CASPT2 to include dynamic correlation
  Multistate = RootHomogenous
  Shift = 0.25
  IPEA = 0.0

4. 关键引用文献与本工作的局限性学术评论

4.1 关键引用文献

本工作建立在过渡金属多光激激发态动力学一系列前沿研究成果基础之上，其中最为关键的支撑性科学文献包括：

1. Juban, E. A.; McCusker, J. K. J. Am. Chem. Soc. 2005, 127, 6857–6865. (文献[21])
   • 学术贡献：该工作首次利用高时间分辨率的瞬态吸收光谱，在实验上确立了 $\text{Cr(acac)}_3$ 的 $^4 ext{T}_{2 ext{g}} \rightarrow ^2 ext{E}_{ ext{g}}$ 的系间窜越时间尺度小于 $100\text{ fs}$。这打破了自旋禁阻转换极慢的传统科学认知，提出了飞秒 ISC 的科学假说。
2. Schrauben, J. N.; Dillman, K. L.; Beck, W. F.; McCusker, J. K. Chem. Sci. 2010, 1, 405–410. (文献[22])
   • 学术贡献：实验观测到了超快 ISC 后 $^2 ext{E}_{ ext{g}}$ 激发态面上高达数皮秒的相干振动（振动相干特征），频率约为 $164\text{ cm}^{-1}$。指出这一特定的扭转振动模式直接构成驱动超快自旋跃迁的“反应坐标”。
3. Ando, H.; Iuchi, S.; Sato, H. Chem. Phys. Lett. 2012, 535, 177–181. (文献[23])
   • 学术贡献：先驱性的理论探索，使用一维核波包动力学模拟了该体系。然而受限于计算算力，其仅将核维度退化为单一的 $\text{Cr-O}$ 全对称伸缩坐标，未能揭示出扭转多维运动的协同作用。
4. Paulus, B. C.; McCusker, J. K. Faraday Discuss. 2022, 237, 274–299. (文献[18])
   • 学术贡献：提出并探讨了利用振动退相干时间（Dephasing Times）来反向推导和辨识第一行过渡金属无机色原体超快动力学反应坐标的系统性方法。

4.2 本工作局限性之独家、深刻评论

尽管本工作利用极其精良的量子化学手段优雅地展示了振动如何驱动系间窜越，但从高标准学术审视的角度来看，依然存在以下具有学术争鸣空间的局限：

4.2.1 气相分子模型的非物理简化与溶剂化效应缺失

本动力学模拟完全是在**气相环境（Gas Phase）**中开展的，即不包含任何显式（Explicit）或隐式（Implicit）溶剂分子。然而在化学实验中，所有的超快激光光谱和催化反应都是在溶液环境（如乙腈、甲醇或水）中进行的。

• 科学影响：溶液环境会引起剧烈的振动能量弛豫（Vibrational Energy Relaxation, VER）和动态退相干（Dephasing）。极性或非极性溶剂分子的碰撞、静电相互作用将极大地加速 $219\text{ cm}^{-1}$ 扭转模式和 $465\text{ cm}^{-1}$ 伸缩模式相干波包的能量消散，使其可能无法像本研究所模拟的那样维持长达 $600\text{ fs}$ 以上的低阻尼相干振荡。这种相干周期的缩短是否会显著降低真实的 ISC 速率，在气相模拟中无法得到答复。

4.2.2 忽略实际动力学演化中的非绝热轨跃（Non-adiabatic Surface Hopping）

由于技术限制，本研究的 BO-AIMD 动力学轨迹被严格局限在单一的绝热 $^4 ext{T}_{2 ext{g}}$ 势能面上（未采用多态非绝热表面跃迁方法，如 AIMS 或 Fewest Switches Surface Hopping）。

• 科学影响：虽然研究者沿绝热轨迹监测了能隙频繁接近零，并依此定性判断波包具有极佳的窜越机会，但他们并未在动力学模拟中直接、实时地计算并观察到量子态密度的自发跨态衰减过程和窜越动力学曲线。真正的非绝热动力学不仅受能隙大小控制，还受到非绝热耦合矢量（NACV）和 spin-orbit 耦合方向空间张量重合度的微观调节。未明确进行多面非绝热跳跃模拟，导致无法计算出定量的超快 ISC 半衰期（例如难以给出一个具体的 “50 fs” 还是 “80 fs” 的理论衰减常数），削弱了与实验时间的精确对比。

4.2.3 简化的 $\text{CAS}(3,5)$ 活性空间无法涵盖激发态自洽极化效应

在 50 条长达 $610\text{ fs}$ 动力学轨迹（共计达 8.5 万次 CASSCF 计算）中，出于算力妥协，研究者只能使用最窄的 $\text{CAS}(3,5)$ 活性空间。这一设置仅包含了 Cr 的五个 $\text{3d}$ 轨道。

• 科学影响：配体骨架上的电子排布在动力学演化中被完全约束在“闭壳层配体”近似下。然而在真实的 $^4 ext{T}_{2 ext{g}}$ 激发态中，由于键长的剧烈拉伸（从 $1.9\text{ Å}$ 涨落至 $2.3\text{ Å}$），必然伴随着显著的、动态的配体-金属电荷再分布（Charge Redistribution）。仅将活性空间局限于 Cr 轨道，割裂了金属与配体之间的动态反馈耦合，会导致势能面的陡峭度（梯度）可能偏离真实值，进而使计算得到的振动频率（如 $219\text{ cm}^{-1}$）带有潜在的方法学系统偏差。

5. 补充理论探索：Jahn-Teller 动力学与自旋-相对论效应的深度交叉

为了给致力于无机过渡金属发光设计和光物理研究的读者提供更为深刻的洞察，以下针对本研究涉及的核心物理现象进行深层次的补充理论剖析：

5.1 过渡金属中 $T \otimes (e+t)$ Jahn-Teller 动力学的多维复杂性

在群论框架下，$\text{Cr(III)}$ 的 $^4 ext{T}_{2 ext{g}}$ 激发态属于 $O_{\text{h}}$ 点群的三重轨道简并表示。根据哈斯-特勒规则，该三重态将与两种不同对称性的核简正模式发生强烈的线性耦合：

1. $e_{\text{g}}$ 简正模式（对称伸缩/剪切畸变）：诱导经典的正方畸变，将 $O_{\text{h}}$ 对称性拉低至 $D_{4 ext{h}}$。这一过程将产生三个等价的正方畸变能量极小值谷（三个经典的 Jahn-Teller 极小角，常被称为“Wells”）。
2. $t_{2 ext{g}}$ 简正模式（扭转/剪切运动）：诱导三角畸变，将点群对称性拉低至 $D_{3 ext{d}}$。这一多维势能面耦合常被称为 $T \otimes (e+t)$ Jahn-Teller 耦合。

本研究所识别出的 $225\text{ cm}^{-1}$ Twisting/Scissoring 主导振动，正是一种含有强烈 $t_{2 ext{g}}$ 对称性的核简正运动。这意味着分子的核波包在激发后，并非静止在 $D_{4 ext{h}}$ 或 $D_{3 ext{d}}$ 的某一个单一、不变的极小值点处，而是顺着多维 Jahn-Teller “双碗状”（Double-well 或 Triple-well）势能面交叉的环形通道（Mexican-hat type topology 的高维扩展）进行极其复杂的相干漫游。这一多维漫游运动，使得体系能频繁、高效地穿越多条势能面交叉分界线，成为了超快 ISC 能够摆脱热激活阻碍、发生无势垒隧穿的终极物理源泉。

5.2 为什么 SOC 随核振动保持静态？spin-orbit 耦合的原子本征论

在量子化学中，自旋-轨道耦合算符 $\hat{H}_{\text{SOC}}$ 具有强烈的单中心本征特征。通常使用 Breit-Pauli 近似或一电子有效核电荷算符来表述：

其中 $\xi(r) = \frac{e^2}{2m^2c^2} \frac{1}{r} \frac{\partial V(r)}{\partial r}$。由于势能梯度 $\frac{\partial V(r)}{\partial r}$ 极其陡峭且强烈局域在重原子核附近，该积分的数值贡献绝大部分来源于电子在靠近重原子核（即 $\text{Cr}$ 的核心区域）的径向波函数重叠。

在动力学演化中，虽然 $\text{Cr-O}$ 键长在 $1.9\sim 2.3\text{ Å}$ 之间剧烈震荡，且配体发生了明显的扭转，但这一键长变化和配体角度弯曲仅仅轻微地扰动了价层配位轨道的微小轨道杂化比例。而铬原子的 $\text{3d}$ 核心轨道本身并没有发生任何收缩或形变。因此，在整个动力学推进中，自旋-轨道耦合的径向主导项几乎没有发生物理改变，从而在理论上完美解释了为什么图 11 中的 SOC 矩阵元始终恒定地维持在 $60\sim 80\text{ cm}^{-1}$ 之间。这一重要结论警示我们：在设计新型 $3 ext{d}$ 发光配合物时，不需要、也不可能通过改变配体刚性来调制 SOC 算符矩阵元本身的大小；调控的唯一关键在于通过配体工程，精确控制目标简正模式的振动频率与势能面交点的空间相对位置，实现“能隙关闭窗口”的理性调制。

5.3 理性分子设计启示：如何通过配体刚性精确控制无机发光配合物的系间窜越速率？

基于本研究建立的“振动驱动能隙调制超快 ISC”理论，可总结出以下极具实操价值的过渡金属 photocatalyst 结构调控设计指南：

+-------------------------------------------------------------+
|                   目标：调控系间窜越速率                     |
+-------------------------------------------------------------+
                               |
              +----------------+----------------+
              |                                 |
              v                                 v
+---------------------------+     +---------------------------+
|       加速系间窜越         |     |       减缓系间窜越         | 
|  (如: 提高光化学敏化效率)  |     |  (如: 延长荧光寿命和发光)  |
+---------------------------+     +---------------------------+
              |                                 | 
              v                                 v
+---------------------------+     +---------------------------+
| 1. 选用更具柔性的配体主链  |     | 1. 在配体环上引入空间位阻 |
|    (引入更低频的扭转振动)  |     |    基团 (如叔丁基, 苯基)  |
| 2. 减弱配体空间阻碍，使其  |     | 2. 增加配体骨架刚性       |
|    能够自由扭转/剪切以闭合 |     |    (Stiffening the Mode)  |
|    4T2g / 2Eg 能隙。      |     | 3. 从而抑制 Twisting 振动 |
+---------------------------+     +---------------------------+

这一调控机理已被实验部分证实：例如在乙酰丙酮配体环上引入庞大的叔丁基（$t\text{-butyl}$）等阻碍基团后，实验测得其 ISC 发生显著变慢，这正是因为巨大的立体位阻强烈地锚定了八面体扭转，使活性扭转模式发生了严重的频率蓝移或物理受限，切断了波包频繁抵达非绝热交叉 seam 点的通道。本工作用严谨的计算量子化学手段，为这一结构微调策略提供了清晰的动力学演化图像与物理逻辑闭环。