来源论文: https://arxiv.org/abs/2606.08350v1 生成时间: Jun 09, 2026 06:05
全电子动力学Bethe-Salpeter方程（Dynamical BSE）：基于原子中心轨道基组的凝聚态激子物理新突破
0. 执行摘要
在现代凝聚态物理与量子化学领域，精确预测固体、分子晶体以及低维材料的电子激发态和光学性质是一项极具挑战性的任务。多体格林函数理论（Many-Body Perturbation Theory, MBPT）框架下的 $GW$ 近似与 Bethe-Salpeter 方程（BSE）方法已被公认为计算中性激子激发的“金标准”。然而，传统的 BSE 计算通常采用静态近似（Static Approximation），即忽略筛选库仑相互作用内核（Screened Coulomb Kernel, $W$）的频率依赖性。这种近似对于激子结合能较小的传统无机半导体（如单晶硅、砷化镓）非常有效，但在面对强激子效应系统（如有机分子晶体、低维纳米材料、核心能级激发）时，其局限性便暴露无遗。当激子结合能 $E_b$ 接近或达到系统的等离子体频率 $\omega_p$ 时，电子-空穴对之间的动态筛选效应（Dynamical Screening Effects）变得不可忽略。
为了攻克这一理论与计算瓶颈，Ruiyi Zhou 等人于 2026 年发表了题为 “All-electron Dynamical Bethe-Salpeter Equation for Extended Systems with Atom-centered Orbital Basis Set” 的重磅研究。该工作首次将基于平面波（Plane-Wave）的有效介电函数方法成功移植并推广到了全电子数值原子中心轨道（Numerical Atom-centered Orbital, NAO）框架中。通过引入谱分解（Spectral Decomposition）技术，巧妙地在非对角、密集的 NAO 辅助基组空间中处理了非线性介电响应变换，并结合对称性适应的不可约布里渊区（IBZ）映射技术，将高密度 $k$ 点网格下的计算开销削减了约 70%。该方法在全电子分子晶体（如萘晶体）的 $BSE@G_0W_0$ 计算中展现出了极高的数值精度与优异的计算效率，为研究强激子关联体系的动力学激发铺平了道路。
1. 核心科学问题、理论基础、技术难点与方法细节
1.1 核心科学问题：静态近似的失效与动态筛选的必要性
在很多体理论中，Bethe-Salpeter 方程描述了双粒子格林函数的演化，其核心是电子-空穴相互作用内核。在静态 BSE（Static BSE）中，假定相互作用是瞬时的：
$$W(\mathbf{r}, \mathbf{r}', \omega) \approx W(\mathbf{r}, \mathbf{r}', \omega=0)$$这一假设在物理上对应于“粒子间关联的建立时间为零”。然而，现实中电子周围的极化云（筛选电荷）建立需要有限的时间，该时间尺度与等离子体震荡周期（$\sim 1/\omega_p$）密切相关。当激子的运动特征时间（与激子结合能 $E_b$ 成反比）与等离子体震荡时间可比拟时，静态近似就会高估筛选作用，从而导致激子结合能偏小，光谱吸收峰位置出现偏差。
解决该问题的科学核心在于：如何在保证计算可行性的前提下，将频率相关的动力学筛选相互作用 $W(\omega)$ 引入到实际大体系的 BSE 计算中？
1.2 理论基础：Shindo 近似与有效介电函数方法
动力学 BSE（Dynamical BSE）本质上是一个耦合自洽特征值问题：
$$A(\omega)X_S(\omega) = E_S(\omega)X_S(\omega)$$$$E_S(\omega) = \hbar\omega$$对于大体系，在频域网格上逐点求解此方程（即精确对角化法）的计算成本是灾难性的。为了解决这一难题，该工作基于 Shindo 近似（Shindo’s Approximation）。Shindo 近似在弱动力学筛选极限下有效，通过引入单粒子格林函数，解耦了双频率极化传播子中的频率依赖性，从而将双频率 BSE 简化为单频率依赖的极化 BSE。
在此基础上，Zhang 等人（PRB 2023）提出了有效介电函数方法（Effective Dielectric Function Approach）。该方法通过将动力学 screened kernel 中的能级差项替换为最低激子态的激子结合能 $E_b$：
$$E_b = E_g^{QP} - E_{S=0}^{sta.}$$从而消除了相互作用对激发频率 $\omega$ 的显式依赖，推导出有效动力学筛选库仑相互作用：
$$\langle \mathbf{v}\mathbf{k}\mathbf{v}'\mathbf{k}' | \hat{W}^{\text{eff}} | \mathbf{c}\mathbf{k}\mathbf{c}'\mathbf{k}' \rangle$$在平面波表象下，利用等离子体极点模型（Plasmon-Pole Model），该相互作用可解析积分，化简为仅含静态反介电函数 $\epsilon^{-1}(\mathbf{q}+\mathbf{G}, \omega=0)$ 的有效介电响应矩阵：
$$\epsilon^{-1}_{\text{eff}}(\mathbf{q}+\mathbf{G}) = 1 - \frac{\hbar\omega_p \sqrt{1 - \epsilon^{-1}(\mathbf{q}+\mathbf{G}, \omega=0)}}{\hbar\omega_p \left[1 - \epsilon^{-1}(\mathbf{q}+\mathbf{G}, \omega=0)\right]^{-1/2} + E_b}$$其中 $\omega_p$ 为系统的等离子体频率。这一巧妙的物理重整化，避开了复杂的频域积分，使得动力学 BSE 的计算开销与静态 BSE 处于同一数量级。
1.3 技术难点：从平面波（PW）到数值原子轨道（NAO）的跨越
尽管有效介电函数法在平面波基组下非常成功，但在数值原子轨道（NAO）表象中，会遇到以下核心技术难点：
- 非对角矩阵与非线性变换的冲突：在平面波基组下，人们通常采用对角近似（Neglect of Local-Field Effects，即仅保留 $\mathbf{G}=\mathbf{G}'$ 的矩阵元），从而可以直接对每个对角元应用上述非线性标度公式。然而，在 NAO 基组中，由于原子轨道的局域性和非正交性，介电矩阵 $\epsilon^{-1}_{\mu\nu}(\mathbf{q}, \omega=0)$ 是高度非对角且密集的，无法直接进行元素级（element-wise）的非线性代数运算。若强行采用对角化近似，将彻底丢失局域场效应（Local-field effects）。
- 布里渊区（BZ）积分的超高计算开销：$GW$ 与 BSE 的计算成本极大程度上取决于 $k$ 点网格的密度。为了获得收敛的光学吸收光谱，通常需要极其致密的 $k$ 点网格，而在每个 $q$ 点上独立计算介电矩阵和进行矩阵对角化会带来难以承受的算力负担。
1.4 方法细节：NAO 谱分解与布里渊区对称性映射
针对上述难点，作者提出了两个极具创造性的解决方案：
1.4.1 NAO 空间中的介电矩阵谱分解
为了在保留非对角局部场效应的同时应用非线性有效介电变换，作者对辅助基组（Auxiliary Basis Functions, ABFs）表象下的静态反介电矩阵进行了谱分解（Spectral Decomposition）：
$$\epsilon^{-1}_{\mu\nu}(\mathbf{q}, \omega=0) = \sum_{\lambda} X_{\mu\lambda}(\mathbf{q}) \tilde{\epsilon}^{-1}_{\lambda}(\mathbf{q}, \omega=0) X^*_{ u\lambda}(\mathbf{q})$$其中，$\tilde{\epsilon}^{-1}_{\lambda}(\mathbf{q}, \omega=0)$ 为静态反介电矩阵在独立筛选通道（本征空间）中的本征值，$X$ 为相应的酉变换矩阵。由于在本征空间中各个筛选通道已经完全解耦，作者可以直接在本征值上应用非线性有效介电函数变换：
$$(\tilde{\epsilon}^{-1}_{\text{eff}})_{\lambda}(\mathbf{q}) = 1 - \frac{\hbar\omega_p \left[1 - \tilde{\epsilon}^{-1}_{\lambda}(\mathbf{q}, \omega=0)\right]}{\hbar\omega_p + E_b\left[1 - \tilde{\epsilon}^{-1}_{\lambda}(\mathbf{q}, \omega=0)\right]^{1/2}}$$变换完成后，再将有效介电矩阵投影回原始的辅助基组空间：
$$[\epsilon^{-1}_{\text{eff}}(\mathbf{q})]_{\alpha\beta} = \sum_{\lambda} X_{\alpha\lambda}(\mathbf{q}) (\tilde{\epsilon}^{-1}_{\text{eff}})_{\lambda}(\mathbf{q}) X^*_{eta\lambda}(\mathbf{q})$$这一数学处理既保证了严格的算符运算，又完整保留了 NAO 表象下的局部场效应。
1.4.2 对称性适应的不可约布里渊区（IBZ）映射
为了降低 $q$ 点网格计算开销，作者利用空间群对称性操作 $\hat{S} = \{R|\mathbf{f}\}$（旋转 $R$ 加分数平移 $\mathbf{f}$），将全布里渊区（BZ）中的任意 $q$ 点投影到不可约布里渊区（IBZ）的 $\mathbf{q}_{\text{IBZ}}$：
$$\mathbf{q} = \hat{S}\mathbf{q}_{\text{IBZ}} + \mathbf{G}$$在 $\hat{S}$ 操作下，Bloch 适应的辅助基组（ABF）变换关系为：
$$\hat{S} P^{\mathbf{q}_{\text{IBZ}}}_{\nu}(\mathbf{r}) = \sum_{\mu} \mathcal{D}_{\mu\nu}(\hat{S}, \mathbf{q}) P^{\mathbf{q}}_{\mu}(\mathbf{r})$$其中变换矩阵 $\mathcal{D}_{\mu\nu}$ 综合考虑了原子置换（Permutation）、几何旋转（Wigner 旋转矩阵 $R^l_{mm'}$）以及 Bloch 相位因子（Phase Factor）。
利用该映射，我们只需计算 IBZ 内各点的有效 screened kernel $W^{\text{eff}}(\mathbf{q}_{\text{IBZ}})$，其余全 BZ 中的矩阵元均可通过代数旋转快速重建：
$$W^{\text{eff}}_{\mu\nu}(\mathbf{q}) = \sum_{\alpha\beta} \mathcal{D}_{\mu\alpha}(\hat{S}, \mathbf{q}) W^{\text{eff}}_{\alpha\beta}(\mathbf{q}_{\text{IBZ}}) \mathcal{D}^*_{ u\beta}(\hat{S}, \mathbf{q})$$由于旋转操作的计算复杂度极低，该技术极大地释放了计算压力。
2. 关键 Benchmark 体系、计算所得数据与性能数据
为了验证全电子数值原子轨道（NAO）动力学 BSE 算法的精度和计算效率，论文选择了经典的单斜相萘分子晶体（Monoclinic Naphthalene, $\text{C}_{10}\text{H}_8$）作为 Benchmark 体系。该体系具有极强的各向异性和显著的激子效应（激子结合能接近 1 eV），是检验动力学筛选效应的绝佳平台。
2.1 体系参数设置
- 结构：单斜空间群 $P2_1/a$，胞内含 2 个萘分子（共 36 个原子，136 个电子），分子呈特征性的鲱骨式（Herringbone）紧密堆积（见论文 Figure 1）。
- 基组：采用 FHI-aims 的
intermediate级别数值原子轨道（NAO）基组。为了提高分辨率恒等（RI）近似的精度，在辅助基组中额外引入了 $4f$ 标性球面谐和函数。 - 网格：自洽 DFT 计算和基础网格采用 $\Gamma$ 点居中的 $5 \times 7 \times 5$ $k$ 点网格。
2.2 核心计算数据：BSE@DFT+$\Delta$ 验证
作为第一步，作者将基于 NAO 的新算法与原始的基于平面波投影全加势方法（PW-PAW，Zhang et al. 2023）进行了直接对比。两套计算均采用相同的“剪刀算符”（Scissor shift $\Delta = 1.55\text{ eV}$）来模拟准粒子能带隙（DFT+$\Delta$ Gap = 4.78 eV）。
表 1：NAO 动力学 BSE@DFT+$\Delta$ 与平面波（PW）结果的基准对比
| 方法 / 来源 | 准粒子带隙 $K_g$ (eV) | 光学带隙 $E_g^{\text{opt}}$ (eV) | 计算得到的激子结合能 $E_b$ (eV) | 输入激子结合能 (eV) |
|---|---|---|---|---|
| 文献 [38] (PW-PAW) | ||||
| DFT+$\Delta$ | 4.78 | NA | NA | NA |
| 静态 BSE | NA | 3.78 | 1.00 | NA |
| 动力学 BSE (输入 1.00) | NA | 3.68 | 1.10 | 1.00 |
| 动力学 BSE (输入 1.14) | NA | 3.67 | 1.11 | 1.14 |
| 本工作 (All-electron NAO) | ||||
| DFT+$\Delta$ | 4.78 | NA | NA | NA |
| 静态 BSE | NA | 3.81 | 0.97 | NA |
| 动力学 BSE | NA | 3.68 | 1.10 | 1.06 (自文献 [38]) |
数据解读：
- 在静态 BSE 下，NAO 算出的光学带隙（3.81 eV）与 PW 算出的（3.78 eV）非常接近，偏差仅为 30 meV。
- 在引入动力学筛选后，两者的光学带隙均红移至 3.68 eV，算得的激子结合能完全一致（均为 1.10 eV）。这一完美的一致性强有力地证明了基于 NAO 谱分解的有效介电响应表征在数值上的精确度。
2.3 核心计算数据：BSE@$G_0W_0$ 实战演示
进一步地，作者将动力学 BSE 方案与单激发 $G_0W_0$ 准粒子计算相结合。相比于常数剪刀差，BSE@$G_0W_0$ 考虑了状态特异（State-specific）且 $k$ 空间依赖的自能修正。
表 3：基于不同 $k_0$ 网格自能修正的 BSE@$G_0W_0$ 计算收敛性
| $G_0W_0$ $k_0$ 网格 | 物理计算步骤 | 准粒子带隙 (eV) | 光学带隙 (eV) | 计算激子结合能 $E_b$ (eV) | 输入 $E_b$ (eV) |
|---|---|---|---|---|---|
| $2 \times 3 \times 2$ | $G_0W_0$ | 5.37 | NA | NA | NA |
| 静态 BSE | NA | 4.39 | 0.98 | NA | |
| 动力学 BSE | NA | 4.27 | 1.10 | 0.98 | |
| $3 \times 4 \times 3$ | $G_0W_0$ | 5.37 | NA | NA | NA |
| 静态 BSE | NA | 4.40 | 0.97 | NA | |
| 动力学 BSE | NA | 4.27 | 1.10 | 0.97 |
数据解读：
- 准粒子带隙的一致性：即便使用较粗糙的 $2 \times 3 \times 2$ $\mathbf{k}_0$ 网格，其 $G_0W_0$ 准粒子带隙（5.37 eV）也已与更细致的 $3 \times 4 \times 3$ 网格完全一致，表明利用粗网格计算自能、并在细网格（$5 \times 7 \times 5$）上插值求解 BSE 的策略完全可行。
- 动力学红移效应：从静态 BSE 到动力学 BSE，光学带隙从 ~4.39 eV 降至 4.27 eV，产生了一个显著的 0.12 eV 的红移，激子结合能从 0.98 eV 增大到 1.10 eV（增幅达 12%）。这在物理上对应于：由于动态极化响应存在滞后性，激发态下的电荷屏蔽不够充分，使得电子-空穴间的有效吸引作用变强（即接近 Bare Coulomb 相互作用），从而增大了激子结合能。
2.4 计算效率与加速比：IBZ 映射分析
引入不可约布里渊区（IBZ）映射后，介电响应矩阵的对角化次数大幅度减少。
表 2：萘晶体中全 BZ 与不可约布里渊区（IBZ）的 $q$ 点数对比
| $q$ 点采样网格 | 全布里渊区（Full BZ）$q$ 点数 | 不可约布里渊区（IBZ）$q$ 点数 | 实际计算量占比 |
|---|---|---|---|
| $2 \times 3 \times 2$ | 12 | 8 | 66.7% |
| $3 \times 4 \times 3$ | 36 | 15 | 41.7% |
| $5 \times 7 \times 5$ | 175 | 52 | 29.7% |
效率评估： 在最终收敛的 $5 \times 7 \times 5$ 密集网格下，全布里渊区包含 175 个 $q$ 点，而 IBZ 仅包含 52 个点。通过 IBZ 对称性映射，直接避免了 70.3% 的介电自能和矩阵对角化计算！正如论文 Figure 3 所示，使用 IBZ 映射和全 BZ 计算所得的光学吸收光谱（介电常数虚部 $\epsilon_2(\omega)$）在整段能量范围内完全重合，精度未发生任何损失。
3. 代码实现细节、复现指南与开源工具链
本工作的算法核心已完全合入到著名全电子密度泛函计算程序 FHI-aims 中。对于量子化学或材料计算的科研人员，若想在延伸体系中开展动力学 BSE 模拟，可参考以下指南。
3.1 核心开源软件及 Repo 链接
- FHI-aims：全电子数值原子轨道计算平台。其内部集成了多体 perturbation 理论（包含自洽和单发射 $GW$、静态/动力学 BSE 表象）。
- 官方主页：https://fhi-aims.org/
- 代码托管与支持：可通过各学术研究机构获取许可。
- Gims：用于准备 FHI-aims 输入文件（
control.in和geometry.in）的图形化前端辅助工具。
3.2 关键数值方法细节：RI 与 LRI
在复现该工作时，须注意 FHI-aims 如何处理三中心积分。为了使 NAO 基组下的四中心库仑排斥积分转化为三中心积分，代码默认使用了**分辨率恒等（Resolution of Identity, RI）**技术，特别是其变体——局域分辨率恒等（Localized RI, LRI）。LRI 限制了辅助辅助基（ABFs）的展开范围，将计算扩展度从原先的 $O(N^4)$ 压制至 $O(N^3)$ 甚至更低，这对于求解静态反介电矩阵 $\epsilon^{-1}_{\mu\nu}$ 极为关键。
3.3 复现运行指南与参数配置说明
以下是复现类似于萘晶体 $BSE@G_0W_0$ 计算所需的核心输入配置。典型的 FHI-aims 计算需要两个关键文件：geometry.in（几何结构）和 control.in（控制参数）。
步骤 1：ground-state PBE 准备与 $G_0W_0$ 自能计算
在 control.in 中，激活自洽场（SCF）计算并设置准粒子能量修正参数：
# 基础物理体系控制
xc pbe
charge 0.0
spin none
use_density_matrix .true.
# 准粒子 $G_0W_0$ 设置
calculate_re_self_energy .true.
ananalytic_continuation pade
n_frequency_points 60
# k 点网格参数（对应于 coarse k0-grid）
k_grid 2 3 2
步骤 2：静态与动力学 BSE 谱求解
在得到准粒子本征值后，在 control.in 中加入激子表象配置：
# BSE 核心开关
calculate_bse .true.
bse_approximation tda # Tamm-Dancoff 近似
# 选择有源状态空间 (价带数与导带数)
bse_valence_states 5
bse_conduction_states 8
# 动力学控制参数
bse_dynamical .true. # 开启有效介电函数动力学修正
bse_plasma_frequency 17.9 # 萘晶体的 plasma 频率 (eV)
bse_input_binding_energy 0.98 # 输入激子结合能 Eb (可由静态 BSE 结果预先给定)
# 细化网格参数与 IBZ 映射
k_grid_bse 5 7 5 # 用于吸收光谱的 dense 细网格
use_symmetry_ibz .true. # 开启 IBZ 對称性映射
步骤 3：数据提取与后处理
运行结束后，FHI-aims 会生成光谱虚部文件 BSE_dielectric_tensor_yy.out（对应各偏振方向）。可以直接利用 Python 的 Matplotlib 对其进行可视化，观察动力学修正带来的红移峰值（如本文 Figure 4 所示偏振于 $b$ 轴的 $4.27 ext{ eV}$ 强激子吸收峰）。
4. 关键引用文献与方法局限性深度评论
4.1 关键引用文献
- 文献 [38] (Zhang, Leveillee, Schleife, Phys. Rev. B 2023, 107, 235205)：本工作最核心的灵感来源，首次提出了基于平面波的有效介电函数方案。
- 文献 [49] (Shindo, J. Phys. Soc. Jpn. 1970, 29, 287-296)：提出了浅能级激子的有效电荷和动态相互作用模型，是动力学 BSE 的理论基石（Shindo 近似）。
- 文献 [33] (Zhou, Yao, Blum, Ren, Kanai, J. Chem. Theory Comput. 2024, 21, 291-306)：发展了基于数值原子轨道（NAO）的周期性 $BSE@GW$ 静态算法框架，本工作是在其基础上的动力学升级。
- 文献 [37] (Bintrim, Berkelbach, J. Chem. Phys. 2022, 156, 044101)：提出了另一种精确计算动力学 BSE 的扩展本征空间方案（包含双激发），但其计算代价极高，不适合固体的超密 $k$ 点采样。
4.2 方法局限性与严苛评论
尽管本工作在技术上取得了令人瞩目的突破，但在物理假设、算法闭环和适用范围方面，依然存在不容忽视的局限性：
1. 对 Shindo 弱极化假设的强依赖
Shindo 近似本质上是一个微扰级的处理，其前提是极化效应（动态筛选）相对于静态筛选是一个“小扰动”。对于具有极高介电常数、或电荷密度剧烈波动的系统（例如发生金属-绝缘体转变附近的关联电子系统），Shindo 近似可能彻底失效，导致红移量严重偏离实验值。
2. Hybertsen-Louie 等离子体极点（Plasmon-Pole）模型的单一化假设
有效反介电公式（式 26 与 29）在推导过程中引入了经典等离子体极点近似。该模型假设介电函数谱主要被一个在特定频率 $\omega_p$ 处的单极点（单一等离子体共振模式）所主导。然而，在以下复杂材料中，等离子体激元色散往往极为复杂：
- 低维纳米结构（如 2D 莫尔超晶格）：等离子体激元呈现非平凡的声学支色散，无法用单一频点 $\omega_p$ 代替。
- 多带半导体：存在多个本征频率相差较大的带间跃迁极化支。在此类体系中，单一极点模型会导致动态极化过度简化，抹杀精细的极化耦合机制。
3. 激子结合能 $E_b$ 的非自洽性（Non-self-consistent Dependency）
目前的方法需要将静态 BSE 计算得到的激子结合能 $E_b$ 作为外界已知数输入到动力学 BSE 计算中（如表 3 所示，动力学计算需手动输入 $0.98 ext{ eV}$ 作为常数）。这种非自洽的迭代流程在美学上并不完美。在极端情况下，如果静态 BSE 对 $E_b$ 估计存在严重偏差，由于该值非自洽更新，可能会将误差直接遗传、甚至放大到动力学结果中。
4. 高阶激子与双激子物理的缺失
由 Bintrim 和 Berkelbach 发展的扩展空间方法表明，动力学筛选的微观物理源自激子与双激发态（Double Excitations）的耦合。本工作所采用的有效介电响应法虽然实现了激子态的红移，但无法直接给出任何双激子态的信息，也无法描述光吸收谱中伴随等离子体激元激发的卫星峰（Satellite Peaks）。对于想要探究深层次多激子产生（MEG）和多体关联动力学的学者，这一方案显得无能为力。
5. 补充探讨：多体格林函数理论与动态筛选物理本质
为了帮助读者更全面地理解这项工作，这里我们做进一步的理论与物理图像延伸。
5.1 从 DFT 到 $GW$+BSE：为什么单粒子图像不够？
在密度泛函理论（DFT）中，Kohn-Sham 轨道是虚拟的非相互作用粒子的波函数，其本征值没有严格的物理含义（除了最高占据轨道）。$GW$ 近似通过引入包含非局域、频率相关自能 $\Sigma(\omega)$ 的单粒子格林函数 $G$，给出了真实的**准粒子（Quasiparticle）**能量。这相当于回答了：“系统增加或减少一个电子时的激发态能量”。
然而，光学吸收（如阳光照射在太阳能电池板上产生电流）对应于电中性激发（Neutral Excitation）。在这个过程中，电子被激发到导带，但在价带留下了一个空穴。这两个带电荷的粒子由于库仑引力而相互吸引，形成了一个束缚态——激子（Exciton）。BSE 的作用正是描述这种电子-空穴的协同关联。如果直接用准粒子能带图（$GW$ 隙）去对比光学吸收谱，由于完全忽略了激子吸引作用，光谱会显著偏蓝。激子结合能 $E_b$ 就是连接准粒子间隙与光学间隙的桥梁。
5.2 动态筛选的生动物理图像
我们可以用一个通俗的物理图像来理解“静态筛选”与“动态筛选”的区别：
- 静态筛选（Static Screening）：想象一个带正电的电荷（空穴）和一个带负电的电荷（电子）在粘性介质（极化电子云）中移动。在静态近似下，我们假定介质的响应是无限迅速的。无论两电荷运动有多快，极化云始终完美地跟随着它们，将两者的相互吸引力牢牢屏蔽，使其保持较弱的相互吸引。
- 动态筛选（Dynamical Screening）：如果这两个电荷由于极强的吸引力而旋转、靠近得非常快（大激子结合能，高频运动），极化介质的电荷来不及瞬间反应（有极化惯性）。在这段短暂的时间窗口内，极化云落在了电子和空穴的后面。由于缺乏足够的极化电荷来削弱两者的库仑引力，电子和空穴感受到了极接近“裸库仑力”的强大吸引。这种由于响应滞后导致局部筛选不充分、从而强化引力的现象，正是激发态红移（带隙减小）的物理根源。
5.3 为什么全电子 NAO 在这一领域潜力巨大？
在以往，材料物理学家更偏爱平面波（Plane-Wave）基组，因为平面波天然具有正交性且非常容易通过快速傅里叶变换（FFT）求解三维泊松方程。然而，平面波基组的尺寸取决于空间的大小（晶胞体积），对于包含大量真空层的低维材料、或者是原子极其密集的分子晶体，平面波基组会包含大量的冗余，计算开销巨大。此外，平面波无法直接描述原子核附近的芯电子，必须使用伪势（Pseudopotentials）。
相比之下，**数值原子中心轨道（NAO）**具有以下无可比拟的优势：
- 全电子特性（All-electron）：不使用伪势，可以对核心层电子（Core electrons）和价电子（Valence electrons）进行同等精度的计算。这使得 NAO 动力学 BSE 能够无缝推广到**核心能级激发（X 射线吸收光谱，XAS）**的研究。在核心激发中，激子结合能可能高达数十至上百电子伏特，动态筛选效应占据主导地位，本工作为其指明了方向。
- 空间局域性与高可伸缩性：NAO 基组仅在有原子的地方有非零值。结合 LRI 技术，其构建的介电矩阵高度稀疏（在大尺寸下），极其适合未来在多节点超级计算机上进行大规模并行。对于含数百原子的超大分子晶体，NAO 将是唯一可行的计算路线。
5.4 总结与学术前沿展望
全电子数值原子轨道（NAO）动力学 Bethe-Salpeter 方程的成功公式化，标志着凝聚态激发态计算领域迈出了坚实的一步。通过谱分解法和 IBZ 映射，该研究彻底攻克了“非对角高维介电变换”与“高密度布里渊区积分”两大技术天堑。
展望未来，该方法的演进方向可能包括：
- 自洽迭代计算 $E_b$：开发自动寻找最低激发能并将其反馈给有效介电算符的自洽算法，彻底告别参数手动输入。
- 超越等离子体极点模型：引入双极点或非解析实频介电扫描，以高精度还原二维超晶格和拓扑材料中精细的等离子体共振特性。
- 光伏材料靶向设计：将该方法应用于非富勒烯有机太阳能电池、超卤化物闪烁体等强激子效应工业产品的能带对齐和激子寿命预测，加速新型清洁能源材料的实用化开发。