Joint Approximate Diagonalization approach to Quasiparticle Self-Consistent $GW$ calculations

本文介绍了一种用于准粒子自洽 $GW$计算的联合近似对角化方法，该方法利用了全动力学自能以及由全格林函数导出的密度矩阵，在实现与标准$\mathrm{qs}GW$ 相当的准确度的同时，提供了与高层级 CCSD(T) 参考值更好的符合度。

要点

想象一下，你正试图为一个庞大且复杂的管弦乐团（一个原子或分子）进行调音，使其奏出完美的音符。在量子物理的世界里，这个“音符”就是将电子从系统中踢出的能量，被称为电离能（Ionization Potential）。

几十年来，科学家们一直使用一种称为 GW 的方法来预测这些音符。然而，这种标准方法就像是通过只听第一小提琴的声音，就假设其余的所有乐器都与它完美同步一样。这就是“单次射击”（single-shot）法：你做一个猜测，计算出音符，然后停止。如果你的初始猜测（“输入”）稍有偏差，最终的音符就会出错。

为了解决这个问题，科学家们开发了一种“自洽”的方法，称为 qsGW。你可以把它想象成一个反馈循环：你演奏一个音符，聆听结果，调整乐器的调音，再次演奏，并不断重复，直到声音趋于稳定。然而，标准的 qsGW 方法有一个捷径。为了让数学计算变得可行，它强行将复杂且变化的管弦乐声压缩成一个简单、静态且对称的形状。这就像是在说：“让我们假装这个管弦乐团只演奏一个完美且不变的和弦吧”，尽管在现实中，声音是动态且杂乱的。

新方法：“联合近似对角化”（JAD）

本文的作者 Ivan Duchemin 和 Xavier Blase 提出了一种调音的新方法。他们不再强迫声音变成简单的静态形状，而是使用了一种名为**联合近似对角化（Joint Approximate Diagonalization, JAD）**的技术。

这里有一个类比：
想象你有一张从奇怪角度拍摄的、模糊且混乱的人群照片。

• 旧方法（标准 qsGW）： 你试图强行让照片看起来像一个完美的、对称的网格。你抹去那些混乱的细节，以使其符合一个简单的规则。
• 新方法（JAD）： 你不是强迫照片发生改变，而是旋转摄像机（数学上的“基底”），直到这个混乱的人群尽可能完美地排列整齐。你没有抹去细节，你只是找到了一个最佳的角度，让每个人都能排成整齐的队列。

在这种新方法中，他们观察特定能量点上的“格林函数”（Green's function，类似于所有可能能量状态的地图）。他们旋转数学上的“摄像机”，直到这张地图看起来尽可能对角化（即变得笔直且清晰）。

核心区别：
这种新方法最重要的一点是，它并没有丢弃那些混乱、动态的细节。 它完整地保留了全套、复杂且随时间变化的“自能”（self-energy，即电子之间相互作用的方式）。它寻找的是观察这种复杂性的最佳角度，而不是将其简化为一个静态的、虚假的版本。

结果：调音管弦乐团

作者通过一个包含 100 种不同分子的“测试集”（GW100 集）测试了这种新方法。

1. 准确度： 尽管他们的新方法与现有的标准逻辑完全不同，但结果却出奇地相似。预测能量水平之间的差异微乎其微（大约只有一座大山中的一粒沙子那么大）。这表明两种方法都在寻找正确的“调音”，只是路径不同。
2. “中间地带”的改进： 他们还尝试了一个混合技巧。在标准方法中，他们计算“密度”（即电子在何处）时，仅仅是统计了占据席位的数量。但在完全自洽的方法中，他们会对整个“声波”进行时间积分。
   • 他们创建了一个新版本，称为 $\gamma$sGWJAD。这个版本通过对完整的、复杂的波进行积分来计算电子密度，而不是仅仅统计座位数。
   • 结果： 这种混合方法恰好落在标准方法与完全复杂方法之间。事实证明，它是所有方法中最准确的一个，甚至比其他方法更接近“金标准”参考计算（CCSD(T)）。

总结

• 问题： 计算电子能量的标准方法要么依赖于糟糕的初始猜测，要么将复杂的物理过程过度简化。
• 解决方案： 一种新的方法（JAD），它在不简化数据本身的前提下，找到了观察复杂数据的最佳“视角”。
• 结果： 它与目前的标准方法一样有效，但保留了更真实的物理特性。
• 额外收获： 通过将这种新方法与一种更彻底的电子计数方式相结合，他们创造了一种“金发姑娘”（意指恰到好处）方案，该方案比标准方法和完全复杂的方法都更准确，更接近真实的实验值。

简而言之，他们找到了一种通过旋转麦克风到完美位置来为量子管弦乐团调音的方法，而不是强迫音乐家们去演奏一首更简单的曲子。

技术摘要

技术摘要：用于准粒子自洽 GW 计算的联合近似对角化方法

问题陈述
本文探讨了从头算多体微扰理论中固有的挑战，特别是用于计算准粒子能量的 $GW$近似。虽然标准的非自洽$G_0W_0$ 方法效率很高，但其结果高度依赖于输入平均场轨道（例如来自 DFT 或 Hartree-Fock）的质量。为了缓解这一问题，研究者采用了自洽方案。

• **全自洽 $scGW$：** 自洽地更新格林函数（$G$）和屏蔽库仑势（$W$）。虽然在理论上是严谨的，但其得到的电离能（IPs）往往与高水平参考方法（如 CCSD(T)）存在显著偏差，通常会低估小分子的电离能。
• 准粒子自洽 $qsGW$：** 一种受限方案，它更新准粒子能量和轨道，但依赖于一个**静态、对称化且与能量无关的自能（$\Sigma$）假设（ansatz）。虽然这使得结果与实验相比比 $G_0W_0$ 更为一致，但也引入了一种舍弃了自能完整动力学特性的近似。此外，$qsGW$和$scGW$ 的结果通常相差数个 0.1 eV，$qsGW$ 往往相对于 CCSD(T) 高估电离能。

作者确定了对一种既能保留完整动力学自能，又无需采用静态对称化假设，同时仍能提供稳定且准确的一体轨道方案的需求。

方法论
作者提出了一种基于**联合近似对角化（JAD）**的一体 GW 格林函数的新方案，称为 $qsGW^{JAD}$。

1. 核心概念： 该方法不再通过对称化的自能构建有效静态哈密顿量（如标准 $qsGW$ 中那样），而是寻求对输入的 Kohn-Sham 分子轨道（MOs）进行一个幺正旋转（$U$），以最小化在一组输入准粒子能量处评估的格林函数矩阵集 $\{G(\varepsilon^{QP}_n)\}$ 的非对角元素。

   • 目标函数为：$\text{argmin}_U \sum_n || \text{off-diag}(U^\dagger G(\varepsilon^{QP}_n) U) ||$。
   • 这种方法保留了完整的动力学自能（$\Sigma_{GW}^{XC}(\omega)$），未做任何近似。其近似仅在于：没有任何单一的旋转基组能够严格地在整个能量谱范围内同时对所有 $G(\varepsilon^{QP}_n)$ 矩阵进行对角化。

2. 自洽循环：

   • 优化的轨道用于更新密度矩阵和非相互作用易受性（$\chi_0$）。
   • 构建新的动力学自能。
   • 从由更新后的格林函数导出的谱函数 $A_n(\omega)$ 的主极点中提取准粒子能量。
   • 重复 JAD 过程以在下一次迭代中更新轨道。

3. 中间方案 ($\gamma sGW^{JAD}$)：

   • 作者引入了一种变体，其中密度矩阵并非简单地通过求和占据准粒子轨道来构建（如标准 $qsGW$中那样），而是通过沿虚轴积分完整的格林函数来构建（如$scGW$ 中那样）。
   • 这捕捉了准粒子峰和非相干背景，从而创建了一个介于 $qsGW$和$scGW$ 之间的方案。

主要贡献

• 新颖的形式体系： 引入了应用于 $qsGW$ 的 JAD 方法，该方法绕过了需要静态对称化自能假设的需求，同时保持了准粒子自洽框架。
• 动力学保留： 该方法允许使用完整的频率相关自能，解决了标准 $qsGW$ 的一个关键局限。
• 混合密度矩阵： 证明了通过对完整格林函数进行虚轴积分来构建密度矩阵，可以得到一种中间方案（$\gamma sGW^{JAD}$），该方案连接了 $qsGW$与$scGW$。

结果
这些方法在 GW100 分子测试集（为排除缺乏全电子基组的分子，集合缩减至 93 个分子）上进行了基准测试，并与参考 CCSD(T) 数据进行了对比。

• $qsGW^{JAD}$ 对比标准 $qsGW$：

  • 新方案产生的结果与传统的 $qsGW$ 方法非常接近。
  • $qsGW^{JAD}$ 与标准 $qsGW$ 之间的平均绝对误差（MAE）为 65 meV（平均符号误差为 3 meV）。
  • 两者相对于 CCSD(T) 的偏差相似（MAE $\approx$ 0.21–0.22 eV，MSE $\approx$ -0.15 eV），表明 JAD 方法成功地在不使用静态近似的情况下复制了标准静态假设方法的准确性。

• $\gamma sGW^{JAD}$ 对比 CCSD(T)：

  • 通过积分更新密度矩阵的中间方案显示出与 CCSD(T) 参考值更好的符合度。
  • MAE 降低至 156 meV，MSE 降低至 62 meV。
  • 这使得结果处于标准 $qsGW$（高估 IP）和 $scGW$（低估 IP）之间，更接近 CCSD(T) 基准。

• 谱函数分析：

  • 在自能在输入基组中具有强非对角性的情况下（例如 CO 的 LUMO），JAD 旋转成功产生了一个由单个准粒子峰主导的谱函数，这与强制对角化（标准 $qsGW$ 的结果）相似，但它是通过轨道优化而非自能近似实现的。

意义与主张
论文声称 $qsGW^{JAD}$ 方法为标准的准粒子自洽方案提供了一种稳健的替代方案。其主要意义在于，它在实现与广泛使用的 $qsGW$ 方法相当的准确性的同时，保留了自能完整的动力学特性，从而避免了静态对称化假设所固有的特定近似。

此外，作者证明了密度矩阵的构建方式至关重要。通过采用由完整格林函数导出的密度矩阵（$\gamma sGW^{JAD}$ 变体），该方法可以进行调节，使其产生的结果比标准 $qsGW$或$scGW$ 都更接近 CCSD(T) 参考值。作者认为，这种中间方案为提高分子系统准粒子计算的准确性提供了一条有前景的路径，且不会产生全自洽 $scGW$ 的全部计算成本或准确性惩罚。