Oscillator Strengths and Transition Dipole Moments from a Simplified Equation-of-Motion Coupled Cluster Formalism within the Frozen-Pair Approximation

本文在 EOM-冻结对耦合簇（EOM-fpCCSD 和 EOM-ptCCSD）框架下，利用避免求解 $\Lambda$ 方程和计算左特征向量的近似方法，推导出了跃迁密度矩阵、偶极矩和振子强度的工作方程，证明了这些模型相比于标准的 EOM-CCSD 能够提供更优的激发态性质。

要点

想象一下，你正试图预测一个分子在受到光照时会如何反应。在化学领域，这就像是在还没调配出新油漆之前，就试图猜出它的颜色。为了准确做到这一点，科学家们使用了一种被称为“耦合簇”（Coupled Cluster）理论的复杂数学方法。它是准确性的“金标准”，但同时也极其昂贵且缓慢——就像是在跑马拉松的同时试图解开一个魔方。

这篇论文介绍了一种解决同样难题的新型、更快速的方法，特别针对那些处于“卡住”状态（即电子以一种棘手的方式配对）的分子。以下是他们所做工作的拆解，使用了简单的类比：

1. 问题所在：“完美”食谱的成本太高

标准方法（称为 EOM-CCSD）就像一位大师级厨师，会对每一种食材进行单独品尝，以获得完美的风味。效果很好，但耗时太长。对于大型分子而言，这种方法慢到无法用于日常实验。

另一方面，较便宜的方法（如 TD-DFT）则像是使用食物加工机：虽然快，但有时会把食材处理错，导致味道不对（结果不准确），尤其是在处理复杂的菜肴时。

2. 解决方案：“冻结对”捷径

作者开发了一种名为 EOM-fpCCSD 和 EOM-ptCCSD 的新方法。

• 类比： 想象一个舞池，一对对舞伴（电子对）正在跳舞。在标准方法中，你必须完美地追踪每一位舞者的脚步动作。而在这种新的“冻结对”方法中，作者说：“让我们把舞伴的手锁在一起，只观察这些‘对’作为一个整体是如何移动的。”
• 通过将这些“对”视为一个单一的、冻结的单元，他们可以忽略大量的冗余数学计算。这使得计算速度大幅提升，同时又没有损失“大师级厨师”级别的准确度。

3. 新的技巧：猜测“左侧”

为了计算分子的发光强度（振子强度）或吸收光的能力（跃迁偶极矩），你通常需要求解方程的两侧：即“右侧”（发生了什么）和“左侧”（输入了什么）。

• 旧方法： 计算“左侧”就像是逐帧倒带电影，以看清演员是如何就位的。这既缓慢又耗费计算资源。
• 新方法： 作者使用了一个聪明的数学捷径（“矩阵逆近似”）。他们不再通过逐帧倒带，而是观察最后一帧，并利用一种聪明的推测来重构开头。
• 结果： 他们完全避免了求解“左侧”方程这一繁重的计算任务，从而节省了更多时间。

4. 测试：水与呋喃

为了验证这个新捷径是否有效，他们对两个分子进行了测试：水（简单分子）和呋喃（一种常出现在有机材料中的环状分子）。

• 他们将“冻结对”的结果与“金标准”（LR-CCSD）进行了对比。
• 结果： 他们的这种新方法与“金标准”几乎完全一致。事实上，对于某些困难的激发态类型（即电子发生双重激发的类型），他们的方法甚至比标准方法更出色、更稳定。
• 他们还测试了两种不同的“地图”（轨道基组）来导航分子：一种是标准地图 (HF)，另一种是优化后的地图 (pCCD)。他们发现，新方法在两种地图上表现同样出色，这意味着它具有很强的灵活性。

5. 核心结论

该论文声称，他们成功为计算分子如何与光相互作用建立了一条“快速通道”。

• 速度： 它避开了计算中最昂贵的部分（求解“左侧”方程和“Lambda”方程）。
• 准确度： 它产生的结果与当今最准确的方法非常接近。
• 可靠性： 即使在标准方法难以收敛（陷入困境）的情况下，它依然表现良好。

简而言之，他们找到了一种方法，可以用更高效的食谱获得超级计算机级别的优质结果，使得研究复杂的电子材料不再需要等待漫长的数学计算过程。

技术摘要

技术摘要：基于方程-运动冻结对耦合集群方法的激发态性质研究

问题陈述
准确描述大分子的电子激发态及其性质（如振荡强度 OS 和跃迁偶极矩 TDM）仍然是当前量子化学方法面临的重要挑战。虽然时间相关密度泛函理论（TD-DFT）具有计算效率高的优势，但它往往受限于泛函依赖性、在共轭体系中对键交替描述不准确的问题，以及在模拟电荷转移、里德堡态和双激发态方面的失效。相反，标准的波函数方法（如方程-运动耦合集群 EOM-CCSD）提供了高精度和系统性的可改进性，但由于其陡峭的标度增长，对于常规应用到大型有机电子材料而言，计算成本过于昂ê。因此，迫切需要一种能够在计算成本与耦合集群理论的高精度之间取得平衡的方法。

方法论
本工作推导了在方程-运动冻结对耦合集群（EOM-fpCC）框架内，特别是针对 EOM-fpCCSD 和 EOM-ptCCSD 模型，关于跃迁密度矩阵的工作方程。该方法构建于对（pair）耦合集群双激发（pCCD）拟设之上，该拟设将簇算符限制在零陈年度扇区（电子对激发），并将其扩展至激发态。

关键方法论特征包括：

• 基态参考： 这些方法利用 pCCD 基态波函数。在 EOM-ptCCSD 模型中，簇算符被划分为内部（pCCD 对）和外部（破裂对）组分。在 EOM-fpCCSD 模型中，投影方程针对每个激发流形进行强制执行，从而有效地解耦了零陈年度和非零陈年度扇区。
• 近似跃迁矩： 为了避免求解计算量巨大的耦合集群 $\Lambda$ 方程以及显式计算左 EOM 特征向量，作者采用了两种近似方法：
  1. 期望值近似： 左特征向量算符 $\hat{\Lambda}$ 被近似为簇算符 $\hat{T}$（即 $\hat{\Lambda}^\dagger \approx \hat{T}$）。
  2. 矩阵逆近似： 左特征向量（$L^\dagger$）通过关系式 $L^\dagger \approx (R^\dagger R)^{-1}R$ 使用右特征向量（$R$）进行近似。
• 轨道基组： 计算同时使用了正则 Hartree–Fock (HF) 轨道和自然 pCCD 轨道，以评估激发态性质对轨道基组选择的依赖性。

主要贡献

• 工作方程的推导： 作者展示了在 EOM-fpCCSD 和 EOM-ptCCSD 形式体系内，关于跃迁偶极矩和振荡强度的形式化推导。
• 计算效率： 所提出的方案无需求解 $\Lambda$ 方程及显式计算左特征向量，与标准的线性响应（LR）方法相比，显著降低了计算开销。
• 基准测试： 本研究在多种基组（cc-pVDZ, cc-pVTZ, cc-pVQZ, 和 aug-cc-pVTZ）下，针对水分子和呋喃，将新开发的模型与标准 LR-CCSD 方法进行了基准测试。

结果

• 准确性： 结果表明，EOM-fpCCSD 和 EOM-ptCCSD 模型对激发态性质的描述较标准 EOM-CCSD 变体有所改进。具体而言，对于偶极强度（DS）和振荡强度（OS），EOM-fpCCSD 相对于 LR-CCSD 参考数据的百分比误差最小。
• 轨道无关性： 研究发现，在 EOM-fpCCSD 和 EOM-ptCCSD 模型中，轨道基组的选择（正则 HF 与自然 pCCD）对激发态描述的影响微乎其微。这与较简单的 LR-pCCD 模型形成对比，并表明即使在 HF 参考态不稳定时，这些方法依然具有鲁棒性。
• 误差指标： 在大多数情况下，TDM 和 OS 的总体误差通常保持在 5% 以下。这些方法对基组大小或弥散函数包含情况没有显著的依赖性。
• 尺寸扩展性（Size-Intensivity）： 作者指出，由于右跃迁密度矩阵中存在不连通项，按设计而言，近似的激发态性质不具备尺寸扩展性。然而，对于所研究的小分子（水和呋喃），由于 LR-CCSD 和标准 EOM-CCSD 产生的性质相同，这一尺寸扩展性问题可以忽略不计。

意义与主张
本文声称，最近引入的单参考耦合集群方法（EOM-fpCCSD 和 EOM-ptCCSD）可以可靠地扩展到复杂电子结构的激发态光谱和性质的研究中。研究表明，这些模型在保持耦合集群理论系统可改进性的同时，通过显著降低计算标度，实现了成本与精度之间的理想权衡。能够使用自然 pCCD 轨道而不损失精度，为在标准 HF 参考态可能不可靠的复杂系统中模拟激发态性质打开了大门。作者得出结论，这些方法特别适用于描述具有显著双激发特性以及电荷转移态的情况，而这些状态通常是其他方法的挑战所在。