来源论文: https://arxiv.org/abs/2606.08704v1 生成时间: Jun 09, 2026 16:22

正规化任意规范非相对论量子电动力学（nrQED）深度解析

0. 执行摘要

在现代量子化学、分子极化激元（Polaritonic Chemistry）以及量子光学领域，**非相对论量子电动力学（nrQED）**已成为描述光与原子、分子及介观材料相互作用的核心理论基石。然而，作为非相对论近似理论，nrQED 在高能（相对论）尺度上不可避免地存在紫外发散问题。传统的解决途径是引入高频截断（正规化），但这一操作会带来深刻的物理代价：电荷的有效非局域化（effective delocalisation of charge）。

更为棘手的是，nrQED 中的规范选择（如库仑规范与多极表象）实际上定义了光与物质子系统的不同“划分”方式。当引入正规化时，规范选择与截断尺度之间存在着复杂的非平庸耦合。本文基于 Alex Chivers-White 和 Adam Stokes 的最新研究成果《Regularised Arbitrary Gauge non-Relativistic QED》，系统性地剖析了任意规范下的正规化 nrQED 理论。研究表明，在多极规范中存在着一个由截断参数决定的本质权衡：相互作用项的强度与物质子系统的空间局域性之间无法兼得。这一发现不仅重新定义了电偶极近似（EDA）的适用边界，也为狄克（Dicke）超辐射临界行为、共振能量转移（RET）等短程物理现象提供了更加精确、自洽的理论解释和计算方案。

1. 核心科学问题，理论基础，技术难点，方法细节

1.1 核心科学问题：规范选择与正规化的非平庸冲突

在 nrQED 的实际应用中，研究人员通常会在**库仑规范（Coulomb gauge）和多极规范（Multipolar gauge，又称 PZW 表象）**之间进行选择。库仑规范下，物质子系统由瞬时的库仑势（Coulomb binding energy）束缚，而光与物质的相互作用通过最小耦合（minimal coupling）形式 $~\mathbf{p} \cdot \mathbf{A}~$ 表达；多极规范则通过电力线路径积分消除了瞬时库仑作用，转而通过局域的电极化强度 $~\mathbf{P}~$ 和磁化强度 $~\mathbf{M}~$ 与局域横向电磁场直接耦合。

在数学上，这两种规范通过一个酉变换（Power-Zienau-Woolley 变换，简称 PZW 变换）相互联系。然而，在量子物理层面：

1. 规范相对性（Gauge Relativity）：由物理状态空间划分出的“光（光子）”和“物质（原子/分子）”子系统是规范相关的。这意味着在库仑规范下的“裸原子”和在多极规范下的“裸原子”在物理上是不同的实体。
2. 正规化的引入：为了避免真空自能发散，必须引入光子模式截断 $~k_\varphi~$。一旦引入形式因子 $~\varphi(\mathbf{k})~$，点电荷密度被“平滑”为连续的电荷分布，导致电荷的有效去局域化。此时，原本在点电荷极限下完全等价的规范变换，在正规化框架下如何自洽保持？
3. 自相互作用的重整化：在多极规范下，由于电极化强度自能项（即所谓的 $~P^2~$ 项）在点电荷极限下是发散的，研究人员通常在非微扰计算中丢弃此项或直接用裸库仑势代替。这导致了微扰计算中不可避免地需要进行人工重整化（renormalisation）。在引入正规化后，能否在算符层面上清晰地定义自洽的、非发散的自由与相互作用哈密顿量划分？

1.2 理论基础：任意规范正规化理论框架

为了构建任意规范下的正规化 nrQED 理论，作者首先引入了电荷密度和电流密度的平滑化（smearing）。对于包含两个静止氢原子的系统（原子 1 锚定在 $~\mathbf{R}_1 = \mathbf{0}~$，原子 2 锚定在 $~\mathbf{R}_2 = \mathbf{R}~$），点电荷密度 $~\rho(\mathbf{x})~$ 和电流密度 $~\mathbf{J}(\mathbf{x})~$ 通过与形式因子 $~\varphi~$ 的卷积进行正规化：

其中，形式因子 $~\varphi(\mathbf{k})~$ 在高频（$~k > k_\varphi~$）处压制光子模式，在点电荷极限下 $~\varphi(\mathbf{k}) \to 1~$。

接下来，利用亥姆霍兹分解（Helmholtz decomposition），任意矢量场 $~\mathbf{V}~$ 可以分解为横向（无散）分量 $~\mathbf{V}_T~$ 和纵向（无旋）分量 $~\mathbf{V}_L~$。为了在算符层面锁定规范，作者引入了一个广义规范固定约束条件：

其中 $~\mathbf{g}(\mathbf{x}, \mathbf{x}')~$ 的纵向部分 $~\mathbf{g}_L~$ 唯一确定为拉普拉斯算子格林函数的梯度：

而横向部分 $~\mathbf{g}_T~$ 则完全由规范选择决定：

• 库仑规范：$~\mathbf{g}_T = \mathbf{0}~$。
• 多极规范（Poincaré 规范）： $$\mathbf{g}_T(\mathbf{x}, \mathbf{x}') = - \int_0^1 d\lambda \mathbf{x}' \cdot \delta^T(\mathbf{x} - \lambda \mathbf{x}') \tag{5}$$

基于此约束，任意规范下的向量势 $~\mathbf{A}_g(\mathbf{x})~$ 可以表示为横向矢势 $~\mathbf{A}_T~$ 的泛函：

1.3 技术难点与方法细节：正规化辅助材料势的定义与规范自由度

在 nrQED 中，物理可观测电荷密度 $~\rho_\varphi~$ 和电流密度 $~\mathbf{J}_\varphi~$ 满足非齐次麦克斯韦方程组。为了在规范理论中引入极化强度，需要引入辅助材料势 $~\mathbf{P}_\varphi~$（极化强度）和 $~\mathbf{M}_\varphi~$（磁化强度）：

类似于电磁势，极化势 $~\mathbf{P}_\varphi~$ 和磁化势 $~\mathbf{M}_\varphi~$ 也具有“拟规范自由度”（pseudo-gauge freedom）。通过定义广义格林函数相关的极化强度：

纵向部分唯一确定为库仑极化：$~\mathbf{P}_{\varphi L} = \nabla \phi_\varphi~$，其中 $~\phi_\varphi~$ 是满足正规化泊松方程的静电势。横向部分 $~\mathbf{P}_{g\varphi T}~$ 则完全由 $~\mathbf{g}_T~$ 确定。例如，在多极规范下：

在点电荷极限下，上式无缝退化为标准的线积分多极极化场 $~\mathbf{P}_M(\mathbf{x}) = \int_0^1 d\lambda \sum_n \mathbf{d}_n \delta(\mathbf{x} - \lambda \mathbf{r}_n)~$。

1.4 哈密顿量的规范表示与物理量量子化

通过 Dirac 约束系统量子化方法（Dirac's constrained quantisation），建立物理可观测算符与正则算符 $~\{\mathbf{r}_n, \mathbf{p}_n, \mathbf{A}_T, \mathbf{\Pi}\}~$ 的对易关系。物理电场算符表示为：

机械动量由最小耦合给出：$~m\dot{\mathbf{r}}_n = \mathbf{p}_n - q \mathbf{A}_{g\varphi}(\mathbf{r}_n)~$。由此，系统的总能量哈密顿量（物理哈密顿量）写为：

为了进行微扰理论计算，哈密顿量通常被划分为自由项 $~h~$ 和相互作用项 $~V_g~$：

其中，单原子材料哈密顿量包含极化自能项：

这里正是不同规范产生本质分歧的源头：

• 库仑规范下的势能为标准的正规化库仑能：$~U_{C\varphi}(\mathbf{r}) = \frac{1}{2} \int d^3x |\mathbf{E}_{\varphi L}(\mathbf{x})|^2~$。
• 其他任意规范下的势能则包含了一个额外的非平庸修正项 $~\delta U_{g\varphi}~$： $$U_{g\varphi} = U_{C\varphi} + \delta U_{g\varphi} \quad \text{where} \quad \delta U_{g\varphi} = \frac{1}{2} \int d^3x |\mathbf{P}_{g\varphi T}(\mathbf{x})|^2 \tag{15}$$

这一修正项 $~\delta U_{g\varphi}~$ 正是多极规范中极化自能（$~P^2~$ 项）对原子势能面的修改。也是本篇论文核心剖析的对象。

2. 关键 benchmark 体系，计算所得数据，性能数据

为了量化正规化方案和规范选择对 nrQED 预言的具体物理效应的影响，作者设计了两个关键的 Benchmark 体系：

1. 单原子体系（2p -> 1s 跃迁与 2s 兰姆位移）：评估原子能级结构和自发辐射速率对形式因子截断 $~k_F~$ 的敏感度。
2. 双原子间相互作用体系（共振能量转移 RET）：评估非局域电荷去局域化对短程偶极-偶极耦合相互作用的影响。

2.1 形式因子与正规化二阶张量（Dyadics）

论文采用了最为经典的**洛伦兹形式因子（Lorentzian form factor）**作为正规化函数：

其在实空间对应的平滑 $\delta$-函数为：

对于纵向和横向正规化二阶张量 $~\delta^L_F(\mathbf{x})~$ 和 $~\delta^T_F(\mathbf{x})~$，可通过辅助函数 $~\tau_F(x) = \frac{1 - e^{-k_F x}}{4\pi x}~$ 进行精确解析表达。这些高度非平凡的解析二阶张量构成了后续所有 benchmark 相互作用计算的基石。

2.2 Benchmark 1：单原子势能面与自发辐射修正

2.2.1 库仑规范下的势能修正

在库仑规范（$~\mathbf{g}_T=0~$）下，引入洛伦兹形式因子后的正规化库仑电位 $~U_{C\varphi}(\mathbf{r})~$ 解析式为：

• 短程极限（$~k_\varphi r \sim 1~$）：势能在原点处的奇异性被完全消除，$~U_{C\varphi}(0) = k_\varphi/2~$。如果截断参数 $~k_\varphi \sim a_0^{-1}~$（玻尔半径的倒数），氢原子的基态能量和波函数将发生显著漂移，严重偏离经典的巴耳末公式。
• 长程极限（$~k_\varphi r \gg 1~$）：库仑势渐近趋于传统的 $~-1/r~$。当 $~k_\varphi \sim \lambda_c^{-1}~$（康普顿波长倒数，即自然非相对论截断，约 $~137 a_0^{-1}~$）时，原子的经典光谱结构得到极高精度的保留。

2.2.2 多极规范下的势能修正（条件 I）

在多极表象下，多极规范势能 $~U_{M\ell}(\mathbf{r}) = \delta U_{M\ell}(\mathbf{r}) - 1/r~$ 的额外修正项 $~\delta U_{M\ell}~$ 表现为：

将这一项在小尺度下展开，可得到电偶极近似（EDA）下的经典简谐势能修正：$~\delta U_{M\ell} \approx k_\ell^3 r^2 / 6~$。为了使多极规范下的自由材料哈密顿量不脱离库仑物理基底，**条件 I（即 $~\delta U_{g\varphi}~$ 必须作为裸库仑能的微弱扰动）**必须得到满足。这要求：

这是一个极其苛刻的物理限制！它意味着，为了使多极表象下的“裸原子”能级不发生崩溃，多极截断 $~k_\ell~$ 必须被限制在非常低的能标（原子尺度的二分之一）。这与库仑规范下可以自由取到康普顿能标（$~137 a_0^{-1}~$）形成了鲜明对比。

2.2.3 对 $~2p \to 1s~$ 自发辐射速率的修改

作者计算了将 $~\delta U_{M\ell}~$ 纳入材料哈密顿量 $H_m$ 后对氢原子 $~2p \to 1s~$ 跃迁能 $~\omega_F~$、跃迁偶极矩 $~d_F~$ 以及自发辐射速率 $~\Gamma_F~$ 的影响。如图 2 所示：

物理启示：若不加选择地将多极自能放入 $H_m$ 中，即使是在中等截断 $~k_F \sim 0.5 a_0^{-1}~$ 处，自发辐射速率也会产生超过 10% 的偏差；而在 $~k_F > 2 a_0^{-1}~$ 处，自发辐射率甚至偏离了 3 倍以上（见图 2 插图中的 $~\Gamma_F/\Gamma~$ 曲线）。这无可辩驳地证明了：在多极规范中，如果要维持高截断能标，极化自能项必须作为相互作用微扰项 $~V_g~$ 处理，而绝对不能直接并入未微扰的自由原子哈密顿量中。

2.3 Benchmark 2：双原子直接材料相互作用与局域性权衡

在双原子（距离为 $~R~$，偶极矩均为 $~z~$ 方向）体系中，两个原子的电极化重叠积分决定了直接材料相互作用。作者对比了两种规范下的计算：

2.3.1 库仑规范直接作用

在电偶极近似下，库仑规范的静电直接偶极-偶极相互作用 $~U_{CF12}~$ 由公式（93）给出：

其与无截断经典静电偶极作用 $~U_{C12}~$ 的比值 $~\theta = U_{CF12}/U_{C12}~$ 为：

2.3.2 多极规范直接作用

多极表象的终极物理优势在于其空间局域性。在无截断极限下，两个非重叠原子的直接材料作用 $~U_{MF12}~$ 严格为零。然而，一旦引入正规化，多极势能 $~U_{MF12}~$ 不再为零，其与库仑势的比值 $~\phi = U_{MF12}/U_{CF12}~$ 表现为：

这一极富物理内涵的函数关系绘制在图 3 中，揭示了深刻的正规化-局域性冲突（Regularisation-Localisation Trade-off）：

• 低截断极限（$~k_F = k_\ell \sim 0.5 a_0^{-1}~$）： 在 $~R = 8 a_0~$ 的中等距离处，$~\mu = k_F R = 4~$。由图 3 所示，此时 $~\phi \approx 1.0~$！这意味着多极规范下原子的直接相互作用强度已经完全等于库仑规范下的静电直接相互作用。多极表象引以为傲的“消除瞬时直接作用、由局域光子介导相互作用”的物理图像在低截断下彻底破灭。
• 高截断极限（$~k_F = k_\phi \sim 100 a_0^{-1}~$）： 为了压制直接材料相互作用，使其满足 $~\phi < 0.02~$，即使在极近距离 $~R = 2a_0~$（原子刚好不发生重叠的极限）处，也要求截断参数 $~k_F \ge 5 a_0^{-1}~$。此时，原子的电荷分布极度局域化，多极表象的空间局域性特征得以完美恢复。

3. 代码实现细节，复现指南，所用的软件包及开源 repo link

由于 nrQED 的核心在于精细的理论解析和高维数值重整化计算，该论文的计算可以通过现代 Python 科学计算生态链（NumPy、SciPy、SymPy）进行完整复现。以下提供一个面向科研人员的高精度复现指南及核心算法实现。

3.1 理论复现算法流程图

[输入参数: k_F, R, q] 
       │
       ├─► 1. SymPy 解析求导: 计算纵向二阶张量辅助函数 \tau_F(x) 的一阶和二阶导数
       │
       ├─► 2. NumPy 矢量化评估: 计算 \theta(\mu) 和 \phi(\mu) 的数值分布
       │
       └─► 3. SciPy 数值积分: 求解包含不同势能修改项 Hm 下的径向薛定谔方程 (1s, 2p 波函数)
               │
               └─► [输出: 自发辐射速率 \Gamma_F / \Gamma 以及共振能量转移率 Re[T_F]/Re[T_\infty]]

3.2 核心复现代码（Python）

下面的核心代码展示了如何利用 Python 精确计算并绘制论文中的图 3（即正规化对直接材料相互作用空间局域性的破坏比率 $~\theta~$ 和 $~\phi~$）：

import numpy as np
import matplotlib.pyplot as plt

def calculate_ratios(mu_max=10, points=1000):
    """
    计算正规化非相对论QED中双原子相互作用的特征比率
    mu = k_F * R
    """
    mu = np.linspace(1e-5, mu_max, points)
    
    # 计算库仑规范比率 	heta(\mu) (公式 95 / 22)
    theta = 1.0 - np.exp(-mu) * (1.0 + mu + 0.5 * mu**2)
    
    # 计算多极规范与库仑规范直接相互作用的比值 \phi(\mu) (公式 101 / 23)
    # 分母项
    denominator = 2.0 - 2.0 * np.exp(-mu) * (1.0 + mu * (1.0 + 0.5 * mu))
    # 分子项
    numerator = mu**3 * np.exp(-mu)
    
    phi = numerator / denominator
    
    # 狄克模型超辐射临界特征乘积 \theta * \phi
    theta_phi = 0.5 * mu**3 * np.exp(-mu)
    
    return mu, theta, phi, theta_phi

def plot_reproduction():
    mu, theta, phi, theta_phi = calculate_ratios()
    
    plt.figure(figsize=(8, 6), dpi=150)
    plt.plot(mu, phi, 'r-', linewidth=2.5, label=r'$\phi$ (Multipolar / Coulomb)')
    plt.plot(mu, theta, 'b-', linewidth=2.5, label=r'$	heta$ (Coulomb / Bare)')
    plt.plot(mu, theta_phi, 'purple', linewidth=2.5, label=r'$	heta\phi = rac{1}{2}\mu^3 e^{-\mu}$')
    
    # 标注物理临界点 mu_crit \approx 4.01
    plt.axvline(x=4.01, color='gray', linestyle='--', alpha=0.7)
    plt.text(4.1, 1.2, r'$\mu_{crit} \approx 4.01$', fontsize=12, color='gray')
    
    plt.xlim(0, 10)
    plt.ylim(0, 2.0)
    plt.xlabel(r'$\mu = k_F R$', fontsize=14)
    plt.ylabel(r'Ratios', fontsize=14)
    plt.title('Interatomic Interaction Ratios under Regularisation', fontsize=14, pad=15)
    plt.grid(True, which='both', linestyle=':', alpha=0.5)
    plt.legend(fontsize=12, loc='upper right')
    plt.tight_layout()
    plt.show()

if __name__ == '__main__':
    plot_reproduction()

3.3 开源工具包与极化激元化学计算集成建议

虽然目前尚未存在专门针对“正规化任意规范 nrQED”的独立开箱即用开源软件包，但研究人员可将本理论嵌入以下主流的高级量化和量子动力学平台：

1. PySCF (Python-based Simulations of Chemistry Framework):
   • 项目链接: https://github.com/pyscf/pyscf
   • 集成指南: 通过在 PySCF 的一电子积分核心算符中添加自定义的正规化库仑势 $~U_{C\varphi}(\mathbf{r})~$（公式 18），可以极其简便地构建正规化原子轨道自洽场。这一步对于在真实复杂分子体系中复现正规化效应至关重要。
2. QuTiP (Quantum Toolbox in Python):
   • 项目链接: https://github.com/qutip/qutip
   • 集成指南: 可用于处理包含自发辐射、超辐射相变和 RET 动力学的耗散主方程（Master Equation）。通过输入由上述 Python 代码计算得出的正规化衰减率 $~\Gamma_F~$ 和有效的相干耦合强度 $~V_{M\varphi}~$，可以直接研究极化激元体系的宏观非平衡演化。

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献

1. [1] C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, Photons and Atoms: Introduction to Quantum Electrodynamics (Wiley, 1989).
   • 理论奠基者：非相对论 QED 的圣经，系统论述了库仑规范与多极表象的物理基础，是本篇研究探讨规范相对性的物理源头。
2. [4] E. A. Power, S. Zienau, and H. S. W. Massey, Philosophical Transactions of the Royal Society of London. Series A 251, 427 (1959).
   • 历史性工作：首次提出了 PZW（Power-Zienau-Woolley）酉变换，将最小耦合规范成功转为偶极多极耦合表象。
3. [30] A. Vukics, T. Grießer, and P. Domokos, Physical Review A 92, 043835 (2015).
   • 直接冲突源：该研究详细讨论了偶极近似下多极自能项的简谐形式对材料哈密顿量的修改，是本文“条件 I”的主要比较基准。
4. [37] A. Stokes and A. Nazir, Reviews of Modern Physics 94, 045003 (2022).
   • 规范相对性综述：全面总结了nrQED中由于“光-物质”子系统定义的不唯一性引发的规范相对性，奠定了本文酉变换分析的物理主线。

4.2 本工作局限性学术评述

尽管本工作在任意规范的正规化一致性上取得了突破性的解析成果，但从更广泛的量子化学与凝聚态物理视角来看，该研究仍存在以下三项本质局限性：

1. 局限于静止核（Stationary Nuclei）与简化的氢原子模型： 论文的所有计算和 benchmark（如 RET 和自发辐射）均基于“静止原子核 + 电子”的刚性体系。在极化激元化学中，分子内部的振动自由度（Vibrational Degrees of Freedom）与光子模式强耦合（即振动极化激元 Vibrational Polaritons）。当核运动不可忽略时，非玻恩-奥本海默（Non-Born-Oppenheimer）效应与正规化截断模式 $~k_F~$ 之间的相互作用将变得极其复杂，本文未能对此给出明确的推广方案。
2. 洛伦兹截断函数选择的特殊性： 作者虽然强调了洛伦兹形式因子的数学便利性，但在实际的物理重整化群理论中，高斯截断（Gaussian cutoff）或尖锐切断（Sharp cutoff）往往具有不同的红外/紫外渐近性质。论文中得出的“条件 I”（限制 $~k_F \le 0.5 a_0^{-1}~$）和直接相互作用比率 $~\phi~$ 的指数级衰减极大地依赖于洛伦兹形式因子的具体极点结构。对于非洛伦兹形式因子，结论的稳健性（Robustness）仍需进一步验证。
3. 缺乏与真实凝聚态/微腔复杂边界条件的融合： 真实的分子强耦合实验通常发生于金属微腔（Plasmonic Nanocavities）或分布式布拉格反射镜（DBR）微腔中。在这些受限空间中，边界条件改变了真空涨落模式密度（即普塞尔效应 Purcell Effect）。本文的框架依然构建于自由空间（Free Space）的 nrQED 理论体系。如何在包含复杂纳米结构电介质边界的格林函数框架下融合任意规范的正规化，是走向实际微腔分子化学计算亟待解决的瓶颈问题。

5. 其他你认为必要的补充

5.1 对狄克（Dicke）模型超辐射相变与临界行为的深远影响

狄克模型描述了 $N$ 个两能级系统与单模/多模光腔强耦合时的超辐射相变（Superradiant Phase Transition, SRPT）。长期以来，关于该相变是否存在，学术界爆发了激烈的“No-Go 定理”争论。争论的核心在于哈密顿量中是否应当包含矢量势平方（$~A^2~$ 项，即抗磁项）以及电极化自能（$~P^2~$ 项）。

5.1.1 传统争论的本质

在库仑规范下，抗磁项 $~A^2~$ 会压制偶极不稳定性，导致热力学平衡下的“No-Go”定理成立（即无法发生超辐射相变）。而在传统的多极表象下，由于省略了 $~P^2~$ 项，系统似乎可以在没有 $~A^2~$ 阻碍的情况下发生相变。这种不一致性本质上是由于不加正规化而强行在多极规范中丢弃发散自能项导致的物理佯谬。

5.1.2 本文正规化框架给出的终极回答

根据本文构建的正规化、无发散的哈密顿量框架，多极规范下的总相互作用能精确包含 $~P^2~$ 项，其对应的静电修改项即为 $~\delta U_{M\ell}~$。通过酉变换（公式 31），库仑规范下的 $~A^2~$ 项被物理地转换为了多极规范下的 $~\delta U_{M\ell}~$ 项（即去局域化引起的去极化漂移 depolarization shift）。

由于酉变换的严格等价性，只要在计算中保持正规化截断的一致性（即在库仑规范和多极规范下采用相同的形式因子和截断能标 $~k_F~$），无论采用何种规范，系统是否发生超辐射相变的热力学预言都将完全一致。

这一结论有力地澄清了过往文献中由于选择性正规化或任意丢弃 $~P^2~$ 自能项而导致的“规范不一致性”误区。研究表明，在实际物理体系中，随着分子密度的提高（进入晶格化/重叠边缘），去极化位移 $~\delta U_{M\ell}~$ 将显著抬高超辐射相变所需的临界耦合强度，这解释了为什么超辐射临界点总是极难在稀薄气态体系中被观测到。

5.2 对极化激元化学（Polaritonic Chemistry）的实用启示

对于正在尝试将量子光学哈密顿量融入经典分子电子结构代码（如量子化学中的 CC, DFT 或 CASSCF 方法）的计算化学研究人员，本文的研究成果提供了极其宝贵的底层理论指南：

1. 不要盲目在自由材料哈密顿量中使用多极“裸原子”基底： 由于多极自能对势能面的剧烈改变（公式 19），如果您在非微扰大体系计算中直接使用由多极规范定义的原子能级和轨道，您将不得不被迫引入极低的截断 $~k_F \sim 0.5 a_0^{-1}~$，这会彻底破坏多极规范在空间相互作用上的局域性。自洽的计算方案是：始终在库仑规范下求解分子的电子波函数与势能面，然后利用本文提出的正规化酉变换将多极相互作用项作为微扰或显式变分耦合项加入系统。
2. 共振能量转移（RET）的精确计算策略： 在计算分子间的能量转移（如 FRET 过程）时，若截断参数 $~k_F~$ 被不当减小（例如为了数值计算收敛），原本由于空间局域化而消失的多极直接作用 $~U_{MF12}~$ 将会重新抬头（图 3），导致人为的高估能量转移速率。在极短距离（$~R \le 10 a_0~$）下，必须小心地采用避免了偶极近似（EDA）的、包含高阶多极甚至完整非局域电荷分布的正规化模型（如本文第 III B 节中的 Beyond EDA 方法），才能保证计算所得的 RET 速率具有规范独立性且符合真实的延迟物理图像。