All-electron Dynamical Bethe-Salpeter Equation for Extended Systems with Atom-centered Orbital Basis Set

本文提出了一种用于扩展体系的全电子数值原子中心轨道实现动态 Bethe-Salpeter 方程的方法，该方法结合了动态屏蔽效应，并通过对萘等分子晶体的应用进行了验证。

要点

想象一下，你正试图预测一群人（电子）在体育场（晶体）中会对一声响亮的欢呼（光）做出何种反应。在量子化学的世界里，这被称为计算“激发态”。

长期以来，科学家们一直使用一种名为 Bethe-Salpeter 方程 (BSE) 的流行方法来解决这个问题。你可以把 BSE 想象成一本规则手册，规定了人群中的两个人——一名啦啦队员和一名喝倒彩的人（电子和由于电子离开而留下的“空穴”）——是如何相互作用的。

问题：“瞬时”与“实时”规则

标准的规则手册假设，当啦啦队员和喝倒彩的人发生相互作用时，它是瞬间发生的。这就像是在说：“如果我挥一下手，你在同一纳秒内就能看到。”这被称为静态近似 (static approximation)。

然而在现实中，存在一个微小的、极其短暂的延迟。人群的反应并不是瞬时的；存在一种涟漪效应。在物理学中，这被称为动力学屏蔽 (dynamical screening)。对于大多数材料，这种延迟非常小，以至于我们可以忽略它。但对于某些材料，比如有机晶体（想想萘，即樟脑丸的成分），这种延迟是非常巨大的。“涟漪”的影响至关重要。如果你忽略它，你对材料吸收光线的预测就会出错。

问题在于，计算这种“实时”延迟的成本极高。这就像是试图用慢动作拍摄体育场里每一个人对每一次欢呼的反应。这需要如此庞大的计算能力，以至于科学家通常无法对大型固体材料进行此类计算。

解决方案：一个更聪明的捷径

由周睿谊（音译）和金井洋介（音译）领导的研究小组开发了一种全新的、超高效的方法，无需使用城市规模的超级计算机即可计算这种“实时”延迟。

他们采用了一种巧妙的捷径方法，这种方法此前仅适用于一种特定的数学类型（使用“平面波”，就像平滑滚动的海浪），并将它翻译成了一种他们称之为数值原子中心轨道 (NAO) 的新语言。

以下是类比：

• 旧方法（平面波）： 想象一下，通过测量平坦网格上每一个点的水深来描述一座山的形状。这很精确，但需要测量数百万个点。
• 新方法 (NAO)： 想象一下，通过在地面上放置一些特定的、精细的雕塑（原子），并观察它们如何组合在一起，来描述同一座山。对于像分子这样复杂的形状，这种方法要高效得多。

作者们成功地教会了他们的“基于雕塑的”系统如何处理这种“实时延迟”（动力学屏蔽），使用的是一种称为有效介电函数 (Effective Dielectric Function) 的方法。他们并没有模拟每一秒的延迟，而是计算了一个捕捉了相互作用本质的单一“平均延迟”值。

“对称性”技巧

即便有了这个新捷径，计算晶体中每一个方向的延迟仍然太慢。因此，他们加入了第二个技巧：对称性映射 (Symmetry Mapping)。

想象一下雪花。它有六个完全相同的臂。如果你知道其中一个臂如何对热做出反应，你自动就知道其他五个臂是如何反应的，因为它们是完全一样的。你不需要测试全部六个臂。

作者们意识到，他们研究的晶体（萘）具有类似的对称性。他们没有计算晶体“地图”（布里渊区）中的每一个点，而是只计算了那些唯一的、不重复的部分（不可约布里渊区）。然后，他们利用数学方法将这些结果“镜像”到整个地图中。

这减少了大约 70% 的工作量，使计算变得足够快且具有实用性。

证明：樟脑晶体

为了证明他们的方法有效，他们以结晶萘作为测试对象。

1. 他们将这种新的“基于雕塑的方法”与旧的“海洋波方法”进行了对比。结果几乎完全一致（误差极小），证明了他们的翻译是成功的。
2. 随后，他们运行了完整的“实时”计算。他们发现，包含延迟（动力学屏蔽）会改变晶体吸收光的颜色。具体来说，它使光吸收的能量发生了约 0.12 电子伏特的偏移。

这为什么重要

这篇论文并不声称今天就能治愈疾病或制造出新的电池。相反，它为研究固体材料（如有机晶体）如何与光相互作用的科学家提供了一个更快速、更准确的工具。

通过使“实时”计算对于复杂的扩展系统成为可能，他们扫清了一个主要的障碍。现在，研究人员可以以前所未有的精度来研究具有强“电子-空穴”相互作用的材料（例如有机电子学中发现的材料），而无需等待数周让计算机完成数学运算。

简而言之： 他们将一个非常缓慢且复杂的计算过程，翻译成了一种更高效的语言，并添加了一个“镜像技巧”来加速，从而让科学家终于能够观察到固体晶体中电子微妙的实时相互作用。

技术摘要

技术摘要：针对扩展系统及原子中心轨道基组的全电子动力学 Bethe-Salpeter 方程

问题陈述
Bethe-Salpeter 方程 (BSE) 是多体微扰理论中用于计算电子激发（通过考虑粒子-空穴相互作用，即激子）的标准方法。然而，标准的 BSE 计算通常采用静态近似来处理屏蔽库仑相互作用核，忽略了其频率依赖性。虽然这种近似对于能量远小于等离子体频率的跃迁是有效的，但对于具有显著激子特征的系统（如分子、低维系统和有机固体，其中激子结合能较大），这种近似就会受到质疑。在这些情况下，动力学屏蔽效应可能会显著改变激发能和光谱形状。

由于需要密集的布里渊区 (BZ) 积分以实现收敛（特别是与 GW 计算结合进行准粒子能量计算时），求解扩展系统的全动力学 BSE 在计算上是极其昂贵的。现有的引入动力学效应的方法（例如在扩展空间中进行精确对角化或使用微扰理论）要么面临高阶计算复杂度的挑战，要么仅限于小尺寸系统。因此，需要一种计算高效且准确的方法，能够在全电子框架下处理扩展周期性系统中的动力学屏蔽，特别是在使用有利于化学应用的原子中心轨道 (NAO) 基组的情况下。

方法论
作者制定并实现了一种用于动力学 BSE 的有效介电函数方法，将最初为平面波 (PW) 基组开发的方案（Zhang 等人，Phys. Rev. B 107, 235205 (2023)）适配到了全电子数值原子中心轨道 (NAO) 框架中。

1. 理论框架：

   • 该方法依赖于 Shindo 近似，以解耦双频率极化函数的频率依赖性。
   • 作者并未求解处理频率依赖性 BSE 所需的复杂耦合特征值问题，而是利用最低激子态的激子结合能 ($E_b$) 取代了屏蔽库仑相互作用中的能量差项。
   • 这导致了一个由有效逆介电函数 ($\epsilon^{-1}_{\text{eff}}$) 推导出的有效频率无关屏蔽库仑算符 ($\hat{W}^{\text{eff}}$)。
   • 有效介电函数是通过等离子体极点模型构建的，其中频率积分通过解析方式完成。

2. NAO 基组中的实现：

   • 作者在其全电子 BSE@GW 代码中实现了该方法，并采用了带有原子中心辅助基函数 (ABF) 的解析度恒等 (RI) 技术。
   • 与通常采用对角近似（通过假设 $G=G'$ 来忽略局部场效应）的平面波表述不同，NAO 实现保留了辅助基组中完整的、密集的、非对角的介电矩阵。
   • 为了处理有效模型所需的介电响应非线性泛函变换，作者对静态逆介电矩阵进行了谱分解。该变换应用于特征值，随后在辅助基组中重建矩阵。这使得能够对局部场效应进行严格处理。

3. 计算优化 (IBZ 映射)：

   • 为了解决屏蔽库仑相互作用 $W(q)$ 所需的密集 BZ 采样瓶颈，作者采用了一种对称性适配策略。
   • 他们仅在不可约布里渊区 (IBZ) 内的 $q$ 点处计算静态介电函数及生成的有效屏蔽相互作用。
   • 通过晶体学对称操作（旋转和相位移动）将 IBZ 矩阵映射回完整的 BZ 矩阵。这避免了对对称等效点的重复计算。

核心贡献

• 公式化： 本文提出了首个用于动力学 BSE 计算的全电子 NAO 基组有效介电函数方法，将该方法从原有的平面波实现扩展到了更广的领域。
• 算法开发： 作者开发了一种特定算法，用于处理辅助基组中的非对角介电矩阵，从而能够包含通常在对角近似中被忽略的局部场效应。
• 对称性适配： 通过集成 IBZ 到全 BZ 的映射来构建屏蔽库仑相互作用，显著降低了扩展系统中实现收敛所需的密集 $q$ 点采样的计算成本。
• 实现： 该方法已实现在 FHI-aims 代码中，并基于其现有的全电子 BSE@GW 基础设施。

结果
该方法在具有强激子效应（结合能 $\sim$ 1 eV）的有机晶体——结晶萘上进行了验证和演示。

1. 与平面波结果的验证：

   • 作者将他们的 NAO 基动力学 BSE 结果（使用“剪切位移”DFT+$\Delta$ 方法）与 Zhang 等人的参考平面波结果进行了对比。
   • 计算的光学吸收光谱显示出极佳的一致性，峰值位置和强度匹配误差在 30 meV 以内。
   • 计算出的激子结合能与参考值之差小于 30 meV，验证了 NAO 实现的准确性。

2. IBZ 映射的效率：

   • 对于萘（单斜晶系，$C_{2h}$ 对称性），对于一个 $5\times7\times5$ 的网格，IBZ 映射将所需的 $q$ 点数量从 175 个（全 BZ）减少到了 52 个（IBZ）。
   • 测试确认，通过 IBZ 映射获得的吸收光谱与通过全 BZ 采样获得的结果完全一致，同时显著加速了屏蔽库仑相互作用的计算。

3. Dynamical BSE@G$_0$W$_0$ 应用：

   • 作者进行了单次 G$_0$W$_0$ 计算，随后进行了静态和动力学 BSE 计算。
   • 与静态 BSE 结果相比，引入动力学屏蔽导致激子峰发生了约 0.12 eV 的红移（从 4.39 eV 变为 4.27 eV）。
   • 这种红移对应于激子结合能的增加，表明在考虑动力学效应时，电子-空穴相互作用得到了增强。
   • 结果与在小分子研究中观察到的微扰研究趋势一致，证实了该方法捕捉扩展系统中动力学效应的能力。

意义与主张
本文声称，这项工作成功地弥合了有效介电函数方法的计算效率与扩展系统全电子 NAO 基组准确性之间的差距。通过利用已建立的平面波结果验证实现过程，并展示其在挑战性有机晶体上的应用，作者建立了一个研究分子晶体及其他扩展材料中动力学屏蔽效应的稳健框架。

作者谦虚地指出，虽然目前的实现是一个重要的进步，但未来的工作将侧重于进一步提高密集 BZ 采样的可扩展性和效率。他们提议引入先进的数值技术，如迭代特征值求解器、粗粒化 k 点网格插值以及 Wannier 插值，以将 BSE@GW 方法的应用范围扩展到日益复杂的材料中。本文并未声称解决了所有的计算挑战，而是提出了一个经过验证的、处理在最相关系统中处理动力学效应的实用方法。