Simple and efficient computational strategies for calculating orbital energies and pair-orbital energies from pCCD-based methods

本文引入了基于对耦合簇双激发（pCCD）拟设及其轨道优化变体的廉价计算策略，用于计算轨道能和对轨道能，从而以较低的计算成本精确预测电离能、电子亲和能和电荷能隙。

要点

想象一下你正在试图理解一个分子的“个性”——具体来说，是它释放电子的难易程度（就像一个慷慨的人），或者它抓取电子的渴望程度（就像一个囤积者）。在化学世界中，这些特性被称为电离能（移除电子的难度）和电子亲和能（原子对额外电子的渴望程度）。

几十年来，科学家们一直使用一套被称为库普曼斯定理（Koopmans' Theorem）的规则来快速估算这些数值。你可以把库普曼斯定理看作是一个“快速且粗略”的经验法则：“如果你知道一个电子坐在座位上的能量，你就可以猜出把它踢出去需要多少成本。”

然而，这个旧规则存在一个缺陷。它假设电子是互不往来的孤独个体。但在现实中，电子是社交型生物；它们不断地相互作用、排斥并围绕彼此起舞。这种“电子社交”被称为相关性（correlation）。当你忽略它时，你的预测可能会产生巨大的偏差，尤其是在用于太阳能电池的复杂有机分子中。

旧工具的问题

为了获得精确的答案，科学家们使用极其精确但极其昂贵的方法。这就像试图计算一场沙尘暴中每一颗沙粒的精确轨迹一样。对于大型分子来说，这太慢也太贵了。

另一方面，“快速且粗略”的方法虽然很快，但因为忽略了电子社交，所以往往是不准确的。

新方案：“配对”方法

本文作者引入了一种基于 pCCD（对耦合簇双激发，pair Coupled Cluster Doubles）的新型、负担得起的策略。

以下是类比：

• 旧方法（哈特里-福克方法/Hartree-Fock）： 将电子视为房间里互不说话的陌生人。你分别计算每个人的能量。
• 新方法（pCCD）： 承认电子是成对出现的（就像舞伴）。该方法不再忽略它们，而是专门关注这些“对”是如何相互作用的。它是一个中间地带：它比“超精确”的方法快得多，但捕捉电子“社交”的能力又比旧的快速方法要好得多。

他们究竟做了什么？

研究人员将这种“侧重于对”的方法应用于一种“修正后的库普曼斯定理”。

1. 升级： 他们对旧的“快速规则”进行了微调，以包含这些电子对的影响。他们不再仅仅观察单个电子的能量，而是观察“对”的能量以及分子其余部分对它的反应。
2. 测试： 他们在两组对象上测试了这种新方法：
   • 简单原子： 如氦（Helium）、氖（Neon）和锌（Zinc）。他们将新的“快速”猜测与昂贵的、超精确的计算以及现实世界的实验进行了对比。
   • 有机分子： 他们研究了 24 种常用于太阳能电池（作为“受体”部分，即捕捉光的部分）的有机分子。

结果

• 对于原子： 新方法表现得非常好。它能够高精度地预测移除或添加电子的能量成本，其表现往往优于旧的“快速”方法，并接近于那些昂贵的计算方法。
• 对于分子： 这正是变得有趣的地方。
  • 旧的“快速”方法（使用标准数学）在预测分子如何接受电子（即“囤积者”特性）方面表现很差。
  • 使用“配对”方法的这种新方法解决了这个问题。它对给予和获取电子都给出了更平衡的视角。
  • 重大胜利： 他们可以非常可靠地预测“能量间隙”（给予和获取电子之间的差异）。这个间隙对于设计更好的太阳能电池至关重要。

这为什么重要？

论文声称，这种新方法是一种快速、廉价且可靠的筛选新材料的方式。

想象一下你正在设计一座太阳能城市，你是一名建筑师。你有成千上万种潜在的建筑模块（分子）可供选择。

• 超精确方法就像雇佣 100 名工程师来测试每一块砖。这很完美，但耗时且昂贵。
• 旧的快速方法就像通过观察外观来猜测砖块的强度。它很快，但你可能会选到一块弱砖。
• 这种新方法则像是一位经验丰富的精明工头，他只需看一眼砖块，就能以极快的速度以 90% 的准确率判断出其强度。

作者总结道，他们的这种方法是一种提供“平衡处理”能量的“低成本”工具。它允许科学家快速筛选数千种有机分子，以寻找有机电子和太阳能电池的最佳候选材料，而无需等待数周让计算机完成计算。

简而言之： 他们找到了一种方法，让快速运行的计算机程序变得足够“聪明”，能够理解电子是如何成对起舞的，从而能以极低的成本对太阳能电池材料做出准确预测。

技术摘要

技术摘要：基于 pCCD 方法的轨道能与对轨道能

问题陈述
有机光伏（OPV）和有机电子器件的发展在很大程度上依赖于对基本电子性质的准确预测，特别是决定 HOMO-LUMO 能隙的最高占据分子轨道（HOMO）和最低未占据分子轨道（LUMO）能量。虽然标准方法如 Hartree-Fock (HF) 和密度泛函近似 (DFA) 通过 Koopmans 定理和 Janak 定理提供了较低的计算成本，但由于缺乏电子相关性、离域误差以及虚轨道的物理意义缺失，它们往往无法可靠地描述 $\pi$ 扩展系统。相反，高精度的量子化学方法对于大规模建模有机电子器件的基础构建块而言，计算成本往往过高。此外，包括基于双电子对（geminal）的方法（如 pair Coupled Cluster Doubles, pCCD）在内的波函数方法通常可以提供占据数，但缺乏用于计算近似电离能（IP）、电子亲和能（EA）和电荷能隙所需的具有计算性价比的方案。

方法论
作者引入了从 pCCD 波函数中导出轨道能和对轨道能的经济型计算策略，具体利用了标准的 pCCD 拟设及其轨道优化（oo）变体。该方法建立在以下理论框架之上：

1. 修正的 Koopmans 定理： 作者提出了一种对标准 Koopmans 定理的修正，以纳入 pCCD 态中的电子相关效应。该方法并非仅仅依赖于 Fock 矩阵的对角元（如正则 HF），而是通过评估 pCCD 的相似变换哈密顿量 $\hat{H}^{(pCCD)} = e^{-\hat{T}_{pCCD}} \hat{H} e^{\hat{T}_{pCCD}}$ 来计算能量差。

   • 对于单电离/单电子附着，能量表达式包含了涉及 pCCD 簇振幅（$t$）和双电子积分的修正项，有效地减去或加上了特定轨道的相关能贡献。
   • 对于双电离（DIP）和双电子附着（DEA），作者推导了包括同自旋（三重态，$m_s=1$）和反自旋（单重态，$m_s=0$）在内的两种情况的表达式，并纳入了源自 (D)IP/(D)EA-EOM-pCCD 哈密顿量对角部分的关联修正。

2. 轨道集： 本研究使用两组不同的轨道对这些模型进行评估：

   • 正则 Hartree-Fock (HF) 轨道。
   • 天然 pCCD 优化轨道 (pCCD)，这些轨道是变分优化的，代表了 (D)IP/D(EA)-EOM-pCCD 模型的对角近似。

3. 参考数据： 为了基准测试新模型，作者利用 (D)IP/(D)EA-EOM-pCCD 实现作为高水平理论参考。他们还针对八种中性原子（He, Be, Ne, Mg, Ca, Ar, Kr, Zn）以及一组与 OPV 相关的 24 个有机受体分子进行了对比。

主要贡献

• 推导修正方程： 本文提供了 pCCD 框架内轨道能和对轨道能的显式推导，将 Koopmans 定理扩展到多行列图像，且无需进行广义 Fock 矩阵的对角化（不同于扩展 Koopmans 定理）。
• 高效近似： 所提出的方法具有 $O(N^2 V^2)$ 的标度性（其中 $N$ 是占据轨道，$V$ 是虚轨道），这比 pCCD 中的轨道优化或标准 EOM-pCCD 的 $O(N^2 V^4)$ 标度性要低得多。
• 系统基准测试： 研究通过多种基组（cc-pVDZ, cc-pVTZ, cc-pVQZ 以及增广版本）对原子和分子系统进行了系统的性能评估。

结果

• 原子系统：
  • 对于电离能（IPs），使用天然 pCCD 轨道的修正 Koopmans 定理提供了最可靠的结果，其误差围绕实验值分布，且随基组增大而减小。它优于标准的 Koopmans 方法和基于 HF 的变体。
  • 对于电子亲和能（EAs），pCCD 优化轨道一致地给出比 HF 更高的数值，并表现出较低的基组依赖性。尽管并非所有测试的原子都有可靠的实验 EA 值，但基于 pCCD 的模型与 EA-EOM-pCCD 参考值表现出良好的一致性。
  • 对于 DIP 和 DEA，这些方法与实验值表现出定性一致，但在某些惰性气体（Ne, Ar, Kr）和重原子（Zn）的 EOM-pCCD 参考计算中遇到了数值收敛问题。
• 分子系统（有机受体）：
  • IPs： 基于正则 HF 轨道的 Koopmans 方法表现良好，与 CCSD(T) 参考值的误差通常在 0.5 eV 以下。然而，将 Koopmans 定理应用于 pCCD 轨道会导致 IP 的高估。
  • EAs： 相反，基于 pCCD 的轨道相比于经常缺乏物理意义的正则 HF 虚轨道，能提供更可靠的虚轨道和 EA 描述。
  • 电荷能隙： 基于 pCCD 轨道能（Koopmans(pCCD) 和 modified Koopmans(pCCD)）的结合，实现了对占据轨道和虚轨道的平衡处理。虽然在 aug-cc-pVDZ 基组下，其平均误差（0.42 eV）略优于 CCSD(T)（0.54 eV），但标准差较大（0.63 eV 对比 0.30 eV）。
  • 成本： 基于 pCCD 的方法以极低的计算成本实现了这些结果。

意义与主张
作者声称，他们的修正 Koopmans 方法，特别是利用天然 pCCD 优化轨道时，提供了对“占据轨道和虚轨道能量的平衡处理”。这种平衡对于预测大型分子系统的电荷能隙至关重要，因为在这些系统中，传统的 HF 无法很好地描述虚轨道，而标准 DFT 则难以处理静态相关。

本文将这些方法定位为并非旨在取代高精度 EOM-pCCD 或 CCSD(T)，而是作为适用于有机电子领域新构建单元和掺杂策略高通量筛选的“经济型计算策略”。作者强调，这些模型为设计阶段提供了 HOMO-LUMO 能隙的“初步估计”，且具有较低的计算成本（$O^2V^2$），使其成为设计有机光伏材料的有力工具。这项工作并不声称解决了现有方法的所有缺陷，而是强调了基于双电子对（geminal）的方法在弥合建模 $\pi$ 共轭系统时在准确性与效率之间差距方面的潜力。