Analytic first order nonadiabatic coupling matrix elements of spin-adapted open-shell time-dependent density functional theory

本文推导、实现并基准测试了自旋适配 X-TDDFT 方法的解析一阶非绝热耦合矩阵元，证明其相较于标准 U-TDDFT 显著降低了误差，并为铜(II) 卟啉等开壳层体系的光物理过程提供了定性正确的见解。

要点

想象一下，你试图预测一个由磁铁制成的复杂机器（比如一个旋转的陀螺）在受到摇晃时会如何反应。在化学世界中，这个“机器”是一个具有未配对电子的分子（如有机自由基或过渡金属配合物），而“摇晃”则是光照射到它身上，使其跃迁到激发态。

科学家使用一种名为**含时密度泛函理论（TDDFT）**的工具来模拟这些反应。可以将 TDDFT 想象为分子的精密天气预报。它预测分子如何运动以及能量如何变化。

然而，这里存在一个问题。标准的 TDDFT（我们称之为U-TDDFT）就像是一个假设风总是沿直线吹的天气预报。它对简单分子尚且有效，但对于具有“未配对”电子的复杂分子（就像我们那个旋转的磁铁陀螺）来说，它就会感到困惑。它将电子的两个“自旋”（我们称之为自旋 A 和自旋 B）视为相互独立，从而导致误差。这就像试图描述两个手牵手的舞伴的舞蹈，但天气预报却假设他们是独自跳舞。

新方案：X-TDDFT

本文作者开发了一种名为X-TDDFT的升级版。这就像是一个新版的天气模型，它理解舞伴们是手牵手的。它强制数学运算遵守量子力学的“自旋”规则。他们此前已利用该方法更准确地预测了这些分子的能量和形状，但还缺少一个关键部分：非绝热耦合矩阵元（NACMEs）。

什么是 NACME？
想象分子是一辆在颠簸道路上行驶的汽车。

• 能量告诉你汽车行驶得有多快。
• 梯度告诉你道路向哪个方向倾斜。
• NACMEs告诉你汽车变道或撞入不同状态的可能性有多大。

在化学中，这种“变道”被称为内转换（IC）。这是分子吸收能量、被激发，然后迅速将能量释放回基态的过程，通常以热的形式而非光的形式释放。如果你的 NACME 计算有误，你可能会认为汽车会保持在车道内，而实际上它即将猛烈地拐进沟里。

他们做了什么？

该团队推导出了数学公式，利用他们新的、具有自旋感知能力的X-TDDFT方法来计算这些“变道”概率（NACMEs）。随后，他们通过两种方式进行了测试：

1. 小型测试（甲醛自由基）： 他们将新方法与“金标准”的超精确计算进行了对比（就像用新 GPS 核对卫星地图）。他们发现，旧方法（U-TDDFT）经常存在巨大偏差——有时误差高达三分之一到三分之二。新方法（X-TDDFT）修正了大部分这些误差，使得预测分子“冷却”速度（内转换速率）要准确得多。在某些情况下，新方法预测的冷却速度比旧方法预测的慢了 100 倍。

2. 大型测试（铜卟啉）： 他们研究了复杂的铜基分子（类似于血液中的血红素，但含有铜）。

   • 旧观点（U-TDDFT）： 预测当分子被激发时，它直接冷却或通过中间态绕行的几率是相等的。
   • 新观点（X-TDDFT）： 预测分子几乎从不直接冷却。它几乎总是选择绕行。
   • 结果： 这彻底改变了这些分子行为的故事。旧方法不仅仅是数字上略有偏差，它搞错了整个故事。它还搞乱了不同版本分子（具有不同化学修饰）之间的比较，使得一种版本看起来比另一种更快，而事实恰恰相反。

结论

该论文得出结论：对于具有未配对电子的分子（如自由基或过渡金属），你不能信任旧的“直线”数学（U-TDDFT）来预测它们如何切换能量状态。

正如你不会用平面地图来导航山脉一样，你也不应该对这些复杂分子使用旧的 TDDFT 方法。新的X-TDDFT方法就像一张三维地形图，揭示了“道路”（能量路径）与以往认为的截然不同。这对于试图设计更好的太阳能电池、LED 或催化剂的科学家至关重要，因为如果你不知道分子会跳入哪条“车道”，你就无法控制其行为。

简而言之： 作者构建了一把更好的尺子，用来测量分子在能量状态之间“跳跃”的情况。他们证明，旧的尺子如此不准确，以至于它讲述的关于这些分子如何工作的故事完全不同，特别是对于那些涉及铜和其他过渡金属的分子。

技术摘要

技术摘要：自旋适配开壳层含时密度泛函理论的解析一阶非绝热耦合矩阵元

问题陈述
自旋适配含时密度泛函理论（X-TDDFT）最近已被确立为计算开壳层分子（如自由基和过渡金属配合物）激发能及解析梯度的稳健方法。通过正确处理自旋污染并隐式包含特定的双激发，X-TDDFT 相较于非限制 TDDFT（U-TDDFT）提供了显著改进。然而，自旋适配对非绝热耦合矩阵元（NACMEs）的影响——这对模拟内转换（IC）速率和非绝热动力学至关重要——此前尚属未知。虽然存在严格的 U-TDDFT NACME 理论，但由于自旋污染，特别是对于以闭壳层到空壳层（CV）激发为主的状态，这些理论往往产生巨大误差。现有的自旋适配 NACME 理论要么局限于半经验哈密顿量，要么仅针对特定的自旋翻转激发，要么缺乏适用于一般开壳层系统所需的理论框架。因此，由于缺乏一致的 X-TDDFT NACME 理论，利用自旋守恒框架对开壳层系统的光物理和光化学进行精确模拟受到了阻碍。

方法论
作者推导并实现了 X-TDDFT 方法的解析一阶 NACME，具体针对基态到激发态（ge）以及激发态到激发态（ee）的跃迁。该推导遵循运动方程（EOM）方法，用其自旋适配对应项替换标准的随机相位近似（RPA）激发和退激发算符。

关键理论组成部分包括：

1. 自旋适配算符：该形式体系利用限制开壳层 Kohn-Sham（ROKS）参考态。激发算符 $O^\dagger_I$ 使用自旋适配的 CV(0) 和 CV(1) 激发构建，其中后者与特定的双激发（$\Psi^{\bar{t}a}_{\bar{i}t}$）混合，以确保正确的自旋对称性。
2. 拉格朗日表述：为了处理激发矢量和分子轨道（MO）系数对核坐标的隐式依赖，采用了拉格朗日方法。这使得能够在不显式微分 MO 系数的情况下计算解析导数，并通过拉格朗日乘子（$Z$ 和 $W$ 向量）引入 Pulay 项。
3. 修正项：X-TDDFT NACME 表示为标准的 U-TDDFT NACME 加上特定的修正项，这些修正项源于：
   • 使用自旋适配参考态（ROKS）而非非限制参考态（UKS）。
   • 修改 CV(1) 激发算符以包含双激发。
   • 由此导致的跃迁密度矩阵（TDMs）和 Casida 矩阵（$A$ 和 $B$ 块）的变化。
4. 实现策略：该实现旨在构建于现有的 U-TDDFT NACME 代码之上。作者证明，额外的计算开销是可以接受的，与解析梯度计算相当，因为它主要涉及求解修正的 Z 向量方程和计算中间密度矩阵。

主要贡献

• 理论推导：本文提供了自旋守恒 X-TDDFT 解析 NACME 的首个严格推导，涵盖了 ge 和 ee 跃迁。
• 算法实现：该方法已在 BDF 软件包中实现。它利用了现有的 U-TDDFT 基础设施，仅需在跃迁密度矩阵和拉格朗日乘子中添加特定的修正项。
• 基准测试：作者针对甲醛自由基阳离子（$CH_2O^+$），将该方法与高水平多参考计算（XMS-CASPT2）进行了基准测试。
• 过渡金属配合物的应用：该方法被应用于铜 (II) 卟啉（CuP、CuOEP、CuTPP），以模拟激发态弛豫路径和内转换速率。

结果

• 基准精度（$CH_2O^+$）：与 XMS-CASPT2 参考数据相比，X-TDDFT NACME 将 U-TDDFT NACME 范数的平均绝对偏差（MAD）降低了约 32% 至 65%（具体取决于泛函）。对于具有显著 CV(1) 特征的状态，观察到了最显著的改进，在这些状态下，U-TDDFT 常因自旋污染而失效。然而，作者指出，对于某些混合态（例如 $1^2A_2$），X-TDDFT 偶尔会产生比 U-TDDFT 更大的误差，但这归因于 U-TDDFT 结果中偶然的误差抵消，而非自旋适配理论的失败。
• 铜 (II) 卟啉：
  • NACME 幅值：与 U-TDDFT 相比，X-TDDFT 预测 $2T_1 \to 2S_0$ 跃迁的 NACME 幅值显著更小（高达一个数量级）。这种减小归因于在自旋适配框架内，纯 CV(1) 态中 $\alpha$ 和 $\beta$ 贡献的近乎完美抵消，而在 U-TDDFT 中，由于轨道形状的差异，这种抵消是不完全的。
  • 内转换速率：NACME 范数的减小导致计算出的 IC 速率急剧下降（高达两个数量级）。
  • 弛豫路径：对于未取代的 CuP，U-TDDFT 预测主要的弛豫路径涉及大致相等的直接 IC 至基态和经由配体场态的间接 IC。相比之下，X-TDDFT 预测经由配体场态的间接路径占主导地位。
  • 取代基效应：两种方法预测的取代卟啉（CuP、CuOEP、CuTPP）IC 速率的相对趋势不同。U-TDDFT 预测 CuOEP < CuTPP < CuP 的趋势，而 X-TDDFT 预测 CuTPP < CuOEP < CuP。这表明 U-TDDFT 的误差不仅很大，而且缺乏系统性，甚至无法捕捉到密切相关分子的相对趋势。

意义与主张
本文主张，X-TDDFT NACME 的开发完成了利用自旋适配 TDDFT 研究开壳层分子光物理和光化学过程所需的基本工具包。作者强调，自旋适配不仅仅是对能量和梯度的修正，对于 NACME 以及由此决定的模拟非绝热动力学和 IC 速率的准确性至关重要。

研究结论指出，在开壳层系统中依赖 U-TDDFT NACME 可能导致对弛豫机制的定性错误预测，以及速率常数上既大又缺乏系统性的定量误差。因此，作者建议普遍在开壳层系统中使用 X-TDDFT NACME，即使在这种情况下 U-TDDFT 和 X-TDDFT 的激发能看起来相似。这项工作为未来在自由基和过渡金属配合物中应用非绝热分子动力学及研究超快光物理奠定了基础。