Excited-state Properties Beyond the Excitation Energy from Orbital-Optimized Density Functional Calculations II: Absorption Spectra

本研究将 Löwdin 关于非正交行列式的形式理论扩展到了投影 augmented-wave (PAW) 框架下的轨道优化密度泛函计算中，证明了尽管该方法在定性上能够重现单行列式激发态的吸收光谱和峰强，但在处理多构型态时表现欠佳，且并未从精确交换或自相互作用修正中获得系统性的改进。

要点

大局观：捕捉跳舞分子的快照

想象一下，一个分子就像是一个微小而复杂的舞蹈团。当你用光照射它时，舞者（电子）会跃迁到更高的能量层，改变他们的舞步。科学家想要预测关于这次跃迁的两件事：

1. 能量： 进行这次跃迁需要多少能量？（就像跳跃的高度）。
2. 亮度： 跳跃时发出的光有多亮？（这被称为“振子强度”或“吸收强度”）。

长期以来，计算机程序擅长预测跳跃的“高度”，但经常在“亮度”上出错。这篇论文介绍了一种计算这些跃迁的新方法，称为轨道优化（Orbital-Optimized, OO）密度泛函理论。

问题所在：“非正交”的舞池

在标准的计算机模型中，科学家通常假设“基态”（休息时的舞蹈）和“激发态”（跳跃时的舞蹈）是完全分离且互不重叠的，就像两个人站在不同的房间里。这使得数学计算变得简单。

然而，这种新的“轨道优化”方法更加真实。它允许电子专门为这次跳跃重新排列自己。问题在于，这种重新排列意味着“休息”状态和“跳跃”状态不再处于不同的房间，而是现在处于同一个房间，并且略有重叠。

类比： 想象你正在测量两个拥抱在一起的人之间的距离。如果你假设他们是分开站立的，你的测量结果就会出错。因为这些状态在“拥抱”（重叠），计算光闪烁的亮度变得非常棘手。以往的研究主要检查跳跃高度是否正确，但并没有检查当状态发生重叠时，亮度计算是否准确。

他们做了什么：测试新方法

研究人员采用了这种新方法，并在一小组分子（水、氨、甲醛、甲醇和乙烯）上进行了测试。他们想看看：

1. 该方法能否正确预测光闪烁的亮度？
2. 使用不同的“工具”（数学基组）是否会对计算产生影响？
3. 改变“游戏规则”（用于描述电子行为的数学公式）是否能修复亮度误差？

结果：喜忧参半

以下是他们的发现，通过简单的语言进行了解析：

1. “工具箱”很重要（基组）
把“基组”想象成相机的分辨率。

• 低分辨率（简单的工具）： 如果你使用简单的相机，你会错过那些“弥散”电子（即那些飘得离分子很远的电子，如里德堡态）的精细细节。此时亮度计算会有偏差。
• 高分辨率（复杂的工具）： 当他们使用非常精细的工具（如平面波或双扩增基组）时，亮度的计算结果变得更加一致。
• 结论： 要想得到正确的亮度，你需要一台高分辨率的相机，特别是对于那些飘得很远的电子。

2. “规则”无法解决所有问题（交换相关泛函）
科学家尝试改变模拟的“规则”（使用不同的数学公式，如 PBE、PBE0，或添加自相互作用修正）。

• 结果： 改变规则有助于在某些情况下修复跳跃的“高度”（能量），但并不能一致地修复“亮度”。
• 类比： 想象你试图修复一张模糊的照片。你尝试了更换镜头、调整光线和添加滤镜。有时照片变清晰了，但通常亮度仍然是错的。并没有一种单一的“神奇规则”能修复所有分子的亮度。

3. 真正的元凶：是“独舞”还是“群舞”（单构型 vs 多构型）
这是最重要的发现。当激发态像是一个独舞者（单一、清晰的构型）时，该方法表现出色。

• 独舞者： 对于简单的跃迁（如氨分子），该方法完美地预测了亮度。
• 群舞者： 对于复杂的跃迁——即电子舞蹈是许多种不同可能性混合在一起的情况（多构型）——该方法无法正确预测亮度。
• 具体的失败案例： 他们发现水、甲醛和乙烯存在巨大误差。在这些情况下，激发态是多种不同舞步的混乱混合体。由于计算机模型强行让状态看起来像是一个单一、干净的“独舞”动作，因此导致亮度预测错误。
• 重叠问题： 他们检查了基态和激发态之间的“拥抱”（重叠）是否导致了误差。他们发现，即使改变了重叠程度，亮度误差依然保持不变。因此，重叠并不是主要问题；“独舞 vs 群舞”的本质才是。

与旧方法（LR-TDDFT）的对比

他们将这种新方法与标准方法（LR-TDDFT）进行了对比：

• 标准方法： 擅长预测复杂“群舞”（如乙烯中的亮闪）的亮度，但对那些飘得很远的“独舞者”（里德堡态）的能量预测较差。
• 新方法 (OO)： 非常擅长预测“独舞者”（里德堡态）的能量，但在处理复杂的“群舞”时，在亮度预测上表现挣扎。

核心结论

这篇论文表明，这种新的“轨道优化”方法是理解分子如何吸收光线的强大工具，但有一个前提条件：

• 它非常适用于简单的、单一的“舞步”（里德堡态）。
• 当舞蹈是多种舞步的复杂混合时（多构型状态），它会表现不佳。
• 仅仅改变数学公式或“相机分辨率”并不能修复复杂舞蹈中的误差。要解决这些问题，我们需要一种能够处理“群舞”而非仅仅是“独舞”的方法。

简而言之，该方法在理解电子跳跃的“高度”和简单跳跃的亮度方面迈出了重要一步，但要准确预测那些最复杂、最混乱的电子舞蹈，它仍需要进一步的改进。

技术摘要

技术摘要：超越激发能的轨道优化密度泛函计算 II：吸收光谱

问题陈述
轨道优化（OO）密度泛函计算通过为每个激发态采用变分优化的轨道，而非将其视为基态的线性响应，提供了一种描述激发态的状态特定方法。虽然这种方法已被证明可以改善激发能，特别是在标准线性响应时间依赖密度泛函理论（LR-TDDFT）由于错误的渐近势行为而在描述里德堡（Rydberg）态方面表现不佳的情况下，但它也为计算跃迁性质引入了显著的复杂性。由于每个电子态由一组不同的分子轨道描述，OO 态通常是非正交的。这种非正交性使得计算跃迁偶极矩（TDM）和振子强度变得复杂。以往的基准测试主要集中在激发能或受限的第一个激发态（HOMO–LUMO 跃迁）的振子强度计算上。因此，OO 方法在更广泛的价层激发和里德堡激发光谱中的振子强度表现，以及非正交性对这些性质的影响，仍不为人所熟知。

方法论
本研究扩展了轨道优化密度泛函计算在计算一组小分子（H₂O, CH₂O, NH₃, CH₃OH, 和 C₂H₄）吸收光谱中的应用，涵盖了高达 10 eV 的价层态和里德堡态。

• 理论框架： 作者采用直接轨道优化方案将激发态定位为能量泛函的驻点（鞍点），以避免变分坍缩。为了处理由此产生的非正交态，研究中将 Löwdin 的非正交行列式形式扩展到了投影自治波（PAW）框架内。这使得在非正交 Slater 行列式之间进行矩阵元计算成为可能。
• 实现： 该方法在 GPAW 软件中实现。跃迁矩阵元在 PAW 形式内推导，考虑了平滑伪波函数与全电子波函数之间的变换。该形式包括了对增广球内重叠矩阵和单电子算符（特别是偶极算符）的修正。
• 计算细节： 计算使用了多种交换相关（xc）泛函：广义梯度近似（PBE）、杂化泛函（PBE0），以及使用全局缩放因子（$\alpha = 1$ 和 $\alpha = 1/2$）带有显式 Perdew–Zunger 自相互作用修正（SIC）的 PBE。在 PAW 框架内，使用了平面波（PW）和线性原子轨道组合（LCAO）（aug-cc-pVDZ 和 d-aug-cc-pVDZ）基组。结果与高水平多参考参考数据（MS-CASPT2, CCSD, CCSD(T), 和 CCSDT）以及 LR-TDDFT 计算进行了对比。
• 自旋纯化： 对于开壳层单线态，通过从混合自旋和三线态计算中推导出的自旋纯化公式获得能量和 TDM。

主要贡献

1. 将 Löwdin 规则扩展至 PAW： 本文推导并实现了专门针对 PAW 形式的非正交行列式的 Löwdin 规则，从而能够使用平面波和网格表示法计算非正交态的跃迁性质。
2. 系统性基准测试： 研究对多种分子（超越 HOMO–LUMO 跃迁）的多个激发态的 OO 振子强度进行了全面评估，并对比了不同的基组和 xc 泛函。
3. 非正交性分析： 本工作明确调查了 OO 计算中基态与激发态之间的非零重叠是否会在振子强度中引入系统误差，这是此前文献中提出的疑虑，但在全吸收光谱研究中尚未得到完全解决的问题。

结果

• 基组依赖性： 弥散里德堡态的振子强度对基组表示具有敏感性。相比于双重增广（d-aug）和平面波（PW）表示，单重增广 LCAO 基组（aug）往往会高估激发能，并对高度弥散态产生不准确的振子强度。然而，明亮的价层跃迁对基组并不敏感。
• 泛函依赖性： 改变 xc 泛函（PBE, PBE0, PBE-SIC/2, PBE-SIC）会改变激发能和特定的振子强度，但不会在整个数据集上产生一致的准确性等级。虽然 SIC 改善了渐近势，但它并未一致地纠正振子强度的误差。
• 单构型与多构型特征： 根据电子性质出现了明显的区别。
  • 单构型态： 对于由单一行列式主导的态（例如氨分子中的许多里德堡态），OO 计算在振子强度和峰位置方面与参考值表现出极佳的一致性。
  • 多构型态： 对于具有强多构型特征或显著构型混合的态（例如 H₂O S3/S4, CH₂O S5, C₂H₄ S1），OO 计算表现出显著差异。具体而言，单行列式描述无法捕捉正确的强度分布，导致在所有测试的泛函中都存在振子强度的过高或过低估计。
• 非正交性与重叠： 分析显示，基态/激发态重叠的大小与振子强度的误差之间没有明显的关联。例如，通过杂化泛函或 SIC 减小重叠，并未系统地改善像 H₂O S3 这样存在问题的态的振子强度。因此，非正交性并未被确定为造成最大偏差的主要来源。
• 与 LR-TDDFT 的比较： OO 方法在描述里德堡态的能量和强度方面通常优于 LR-TDDFT。相反，对于具有强构型混合的明亮价层跃迁，LR-TDDFT 通常能提供更好的描述，而单行列式 OO 方法在此类状态下表现挣扎。

意义与结论
本文得出结论，虽然 OO 密度泛函计算是预测单构型里德堡态激发能和吸收光谱的稳健工具，但其对振子强度的定量可靠性受限于 Kohn–Sham 的单行列式性质。当激发态由单一主导构型充分描述时，该方法表现良好；但对于具有强多构型特征的态，无论选择何种 xc 泛函或基组，该方法都会失效。

作者断言，在特定跃迁（如水 S3/S4、甲醛 S5、乙烯 S1）中观察到的巨大振子强度误差，是单行列式近似的内在特性，而非非正交性或基组不完备性的产物。因此，对于多构型特征显著的系统，应将 OO-KS 光谱用于定性解释（例如峰位置），而非关于强度的定量解释。研究表明，未来的改进需要将 OO 框架扩展到显式处理多构型激发态，而不仅仅是依靠泛函或基组的修改。