$G_0W_0$@HF and BSE methods in periodic systems from Hartree-Fock theory: gaussian orbital and density fitting approach

本文提出了一种基于高斯轨道和密度拟合的、适用于周期性体系的$G_0W_0$@HF 与 Bethe-Salpeter 方程框架，该框架通过采用不含等离激元极近似的确切 RPA 屏蔽以及针对虚态的混合收敛策略，修正了半导体和氧化物中 Hartree-Fock 方法对带隙和价带宽度的高估。

要点

想象一下，你正在试图理解一种固体材料（如钻石或一块硅）在受到光照或通电时的行为。为此，科学家需要计算材料内部电子的精确能级。将这些能级想象成电子居住的摩天大楼中的“楼层”。如果你确切知道楼层的位置，你就知道这座建筑如何运作。

几十年来，绘制这些楼层的标准方法是使用一种称为密度泛函理论（DFT）的方法。然而，DFT 就像使用一张略微模糊的地图；它能正确把握建筑的大致形状，但往往无法精确捕捉楼层的高度。为了获得更清晰的图像，科学家使用一种更先进的技术，称为GW（得名于方程中的符号 G 和 W）。这种方法就像从模糊的草图切换到高清 3D 模型，但它计算成本极高，并且通常需要使用一种特定的数学“网格”（称为平面波），这对于某些类型的材料来说难以处理。

新方法：不同的视角
本文由 Charles H. Patterson 撰写，介绍了一种构建这种高清 3D 模型的新方法。作者没有将标准的模糊地图（DFT）作为起点，而是从一个不同、非常锐利但过于僵硬的地图开始，称为哈特里 - 福克（HF）。

• 起点的缺陷：哈特里 - 福克方法就像是用一把过于严格的尺子绘制的地图。它预测楼层之间的距离过远（“带隙”过大），房间过宽（“带宽”过大）。如果仅使用这张地图，你的预测将是错误的。
• 解决方案：作者采用了一种巧妙的策略。他们从这张严格的哈特里 - 福克地图开始，然后应用一个“校正透镜”（GW 方法）来修正误差。论文表明，这个校正透镜非常擅长将过宽的房间缩小回其真实尺寸，从而生成一张与实验现实非常吻合的最终地图。

工具：高斯轨道与密度拟合
大多数 GW 计算使用“平面波”（像无限扁平 sheets 的网格）来描述电子。本文则使用高斯轨道。

• 类比：想象描述一座复杂的雕塑。平面波方法就像试图通过堆叠数百万个扁平的方形瓷砖来描述它。而高斯方法则像使用柔软、圆润的黏土块，这些黏土块可以完美地贴合雕塑的曲线。对于复杂的分子和晶体，这通常更高效。
• 密度拟合：为了让数学运算能够处理这些黏土块而不导致计算机崩溃，作者使用了一种称为密度拟合的技术。这就像一种“压缩算法”。该方法不是计算每一对黏土块之间的相互作用（这将耗时无穷），而是将它们分组为簇，并计算整个组的相互作用。这就像通过称量几个有代表性的人并相乘来估算人群的重量，而不是逐个称量每个人。

“无近似”技巧
在这些计算中，一个常见的捷径是“等离激元极点近似”。

• 类比：想象你试图预测鼓的声音。捷径方法说：“让我们假设鼓只发出一个特定的音符，忽略其余部分。”这很快，但会遗漏细微差别。
• 论文的声明：本文避免了这种捷径。它计算了鼓的完整、复杂的声音（电子相互作用的完整频率依赖性），而不假设它只是一个音符。这更准确，但需要为材料结构中的每个点求解一个巨大而复杂的谜题（Bethe-Salpeter 方程）。

他们发现了什么？
作者在四种材料上测试了这种新方法：钻石、硅、氧化镁（MgO）和二氧化钛（TiO2）。

1. 钻石和硅：标准的哈特里 - 福克方法预测“房间”（价带）过宽，大约宽了 25%。新方法对此进行了修正，将其缩小以完全匹配实验测量的结果。
2. 氧化物（MgO 和 TiO2）：该方法成功预测了能隙（楼层之间的距离）以及材料如何吸收光。虽然预测的能隙略大于实验中观察到的值（这是该领域的常见问题），但能级地图的整体形状非常准确。
3. 光吸收：在模拟这些材料如何吸收光（其“光学光谱”）时，该方法非常准确地再现了峰值的位置（吸收的颜色）。然而，对于氧化物，该方法预测的光吸收强度略强，就像麦克风可能捕捉到稍显过大的声音一样。

结论
本文证明，可以通过从严格的“哈特里 - 福克”模型开始并应用复杂的"GW"校正，同时使用灵活的“高斯”数学语言和智能的“压缩”技术（密度拟合），来构建固体中电子能量的高度精确、高清地图。它证明，你不需要标准的“平面波”网格也能获得出色的结果；事实上，这种替代方法可以修正起始方法的特定误差，从而产生与真实世界实验相符的结果。

技术摘要

技术摘要：基于 Hartree-Fock 理论的周期性体系 GoWo@HF 与 BSE 方法

问题陈述
从头算 GW 方法和 Bethe-Salpeter 方程（BSE）是计算晶体材料准粒子（QP）能级和光学性质的标准工具。然而，绝大多数现有实现依赖于平面波（PW）基组和密度泛函理论（DFT）哈密顿量（通常为 LDA 或 PBE）作为起点。尽管这些方法有效，但它们通常需要较大的基组以及特定的近似（如等离子体极化子近似，PPA）来高效处理介电函数的频率依赖性。此外，DFT 起点可能会在能带结构中引入特定的误差，这些误差难以与多体修正区分开来。有必要验证并将这些方法扩展至高斯轨道（GO）基组，后者是量子化学中的标准基组，且对于具有大晶胞或开放体积的体系具有优势，同时探索从 Hartree-Fock（HF）理论而非 DFT 出发作为起点的功效。

方法论
作者提出了一个框架，用于使用高斯轨道基组和密度拟合（恒等式分解）方法对周期性体系进行$G_0W_0$（记为$GoW_0$）和 BSE 计算。该实现利用了Exciton代码，并具有以下关键技术特征：

• Hartree-Fock 起点： 与大多数固体 GW 应用不同，未微扰哈密顿量和单粒子轨道源自 HF 理论。作者承认 HF 通常会高估带隙和价带宽度，但认为随后的$GoW_0$自能修正可以将这些量重整化以匹配实验数据。
• 密度拟合： 为了管理周期性体系中四中心积分的计算成本，该方法采用基于库仑度规的密度拟合。波函数乘积被投影到辅助高斯基组上，从而降低了积分复杂度。库仑矩阵在该辅助基组中表示，其逆矩阵用于构建屏蔽相互作用。
• 无等离子体极化子近似的 RPA： 屏蔽库仑相互作用$W$通过关系式$W = v + v\Pi v$计算，其中$\Pi$是在随机相位近似（RPA）水平下处理的极化率。关键在于，该方法在布里渊区独特的有限波矢（$Q$）处求解 RPA 方程（对角化大型 RPA 哈密顿量）。这避免了使用等离子体极化子近似，从而在 RPA 和基组限制内对$W$的频率依赖性进行精确处理。
• 混合自能策略： 为了解决自能相对于虚拟态数目的收敛问题，采用了一种混合策略。自能的大部分使用$GoW_0$计算，仅使用有限的虚拟态集合（足以满足 RPA 屏蔽需求）。来自高能虚拟态的剩余贡献使用二阶微扰理论（$\Sigma^{(2)}$）进行评估。这种方法在计算成本与收敛性之间取得了平衡。
• BSE 与 TDA： Bethe-Salpeter 方程采用线性响应运动方程方法构建。采用了 Tamm-Dancoff 近似（TDA），忽略了共振项与反共振项之间的耦合，这被认为足以描述光吸收性质。
• 小$Q$极限： 对$Q \to 0$极限应用了特殊处理，以解决自能和屏蔽电子 - 空穴相互作用中的发散问题，利用了类似于立方固体中交换能计算的方法。

关键结果
该方法应用于元素半导体（金刚石、硅）和宽禁带氧化物（MgO、锐钛矿和金红石$TiO_2$）。

• 能带结构重整化：
  • 金刚石和硅： HF 理论显著高估了价带宽度（Si 和金刚石约为 25%）。$GoW_0@HF$修正成功地将这些宽度重整化，使其与实验光电子能谱数据高度一致（例如，Si 的价带宽度从 HF 的 17.00 eV 降低至$GoW_0@HF$的 13.03 eV，与实验值 12.5 eV 吻合）。
  • 带隙： 虽然 HF 高估了带隙，但$GoW_0@HF$得出的带隙通常大于$GoW_0@LDA$的结果。对于金刚石，间接带隙为 5.75 eV（实验值为 5.48 eV）；对于硅，为 1.38 eV（实验值为 1.17 eV）。对于 MgO，带隙为 9.13 eV，高于实验值 7.77 eV（即使考虑了零点重整化）。
• 光学性质（BSE）：
  • 通过 BSE-TDA 计算的介电函数在金刚石和硅的峰位上与实验结果吻合良好。
  • 对于氧化物（MgO 和$TiO_2$），峰位相当准确（与实验值的偏差在 0.1–0.4 eV 范围内），但该方法系统性地高估了振子强度（强度）。这种高估归因于使用了 HF 波函数，HF 波函数比 DFT 波函数更局域化，导致 BSE 核中的电子 - 空穴相互作用更强。
  • 对于$TiO_2$，计算光谱相对于实验观察到了向低能方向的小能量偏移（0.2–0.4 eV），作者将其归因于屏蔽电子 - 空穴吸引的小$Q$极限中的近似处理。

意义与主张
本文建立了一种稳健的密度拟合、高斯轨道方法，用于周期性体系中的$GoW_0$和 BSE 计算，且不依赖等离子体极化子近似。作者证明了：

1. HF 起点是可行的： 尽管 HF 在带隙和宽度方面存在已知缺陷，但$GoW_0$自能修正能够将 Si 和金刚石的能带结构重整化，得出与实验良好一致的结果，在带宽方面与$GoW_0@LDA$相当或更优。
2. 密度拟合的高效性： 该方法成功处理了周期性 GW/BSE 计算的计算需求，减少了所需的积分数量，使该方法适用于具有大晶胞的体系（如$TiO_2$），在这些体系中平面波方法可能效率较低或需要不同的收敛策略。
3. 无需等离子体极化子近似： 通过在有限$Q$矢量处对角化 RPA 哈密顿量，该方法捕捉了屏蔽相互作用的完整频率依赖性，为基于 PPA 的方法提供了更严格的替代方案。

作者得出结论，虽然$GoW_0@HF$方法得出的带隙通常略大于实验值（且大于$GoW_0@LDA$），但它为绝缘材料中的准粒子激发和光学光谱提供了一致且准确的描述，验证了高斯基组和 HF 起点在固体多体微扰理论中的应用。