Explicit Two-Sided Eigenvalue Bounds for Schrödinger Operators with Singular Potentials via Finite Element Method

技术摘要：通过有限元方法获得奇异势 Schrödinger 算子的显式双侧特征值界

问题陈述
本文解决了在无界区域 $\mathbb{R}^N$（$N=2, 3$）上计算 Schrödinger 算子 $H = -\Delta + V$ 的特征值 $\lambda_k$ 的显式、可计算双侧界的问题。主要挑战在于处理库仑型奇异吸引势，具体为 $V(x) = -Z/|x|$（模拟氢原子）和双中心势 $V(x) = -Z_1/|x-a_1| - Z_2/|x-a_2|$（模拟 $H_2^+$ 分子离子）。

现有的数值方法通常通过协调 Galerkin（Rayleigh–Ritz）方法提供上界。获得严格的下界要困难得多，特别是在无界区域和势函数下方无界的奇异势情况下。经典的下界技术（例如 Lehmann–Maehly–Goerisch、Kato–Temple）通常需要非协调方法或特定条件，而当势函数 $V$ 为负且奇异时，这些条件往往失效。此外，用于下界的标准有限元方法通常依赖于将势函数按其负部分的 $L^\infty$ 范数进行平移（$\gamma_h = \|c_h^-\|_{L^\infty}$）。对于库仑奇点，随着网格细化，$\gamma_h$ 会发散（$O(h^{-1})$），导致下界无法收敛。

方法论
所提出的算法结合了区域截断与复合增强 Crouzeix–Raviart（CECR）有限元方法的扩展。该方法由三个主要部分组成：

1. 区域截断：将无界区域 $\mathbb{R}^N$ 截断为有界区域 $D(R)$。

   • 上界：通过 Dirichlet 边界条件（在 $\partial D(R)$ 上 $u=0$）使用标准协调 $P_1$ Galerkin 单元获得。
   • 下界：通过 Neumann 边界条件（$\partial u/\partial n = 0$）获得。一个关键引理表明，如果外部势最小值 $\sigma_{ext}(D(R))$ 超过第 $k$ 个特征值，则 $D(R)$ 上的 Neumann 特征值即为全空间特征值的有效下界。该条件利用计算出的上界进行后验验证。

2. 适用于变号势的 CECR 框架：

   • 该方法采用对空间框架，将函数提升为对 $(u, \Pi_{0,h}u)$，从而将动能（梯度）与反应能（势）分离。
   • 为了在不产生发散平移的情况下处理奇点，作者引入了单元反应常数 $\Gamma_h := \max_K c_{h|K}^- \cdot h_K^2 / \pi^2$。与全局 $L^\infty$ 平移 $\gamma_h$ 不同，即使对于库仑奇点，随着网格细化，$\Gamma_h$ 也会消失（$O(h^2)$ 或 $O(h)$，取决于分级情况）。
   • 对势函数施加一个固定平移 $\sigma > 0$（独立于网格尺寸 $h$），以确保平移后的双线性形式的强制性。这避免了需要依赖于网格且会发散的平移。

3. 扰动分析：

   • 连续势 $V$ 由其分段常数单元平均值 $c_h$ 近似。
   • 一个扰动引理界定了用 $c_h$ 替换 $V$ 所引入的误差。在奇点附近特定的网格分级假设下，扰动比 $\epsilon_h$ 被证明收敛于零（在分级网格上为 $O(h^2)$）。

主要贡献
本文声称具有以下具体贡献：

• 首个显式收敛下界：提出了首个数值算法，为无界区域上具有无界奇异势的 Schrödinger 算子提供显式、收敛的双侧特征值界。
• 对空间抽象定理：开发了一种对空间框架，将 CECR 分析提升到双分量空间，允许对变号势的下界定理进行自包含的极大极小证明。
• 单元反应常数（$\Gamma_h$）：引入 $\Gamma_h$ 以替代发散的 $L^\infty$ 平移。这使得 CECR 间隙在分级网格上以 $O(h^2)$ 的速率闭合，从而避开了对不可计算的范数常数的需求。
• 固定平移双侧界：构建了一个完全可计算的包围区间 $L_k \le \lambda_k \le U_k$，无需对势的负部分进行 $L^\infty$ 控制。
• 应用于 $H_2^+$：将理论专门化于双中心库仑势，提供了 2D 和 3D 中允许平移和扰动参数的显式公式。

数值结果
该算法在四个测试案例中得到了验证：

1. 2D 氢原子（$V = -1/|x|$）：该方法成功地将精确基态 $\lambda_1 = -1$ 括定在内。随着网格细化，认证的下界区间显著变窄，最优平移界在最精细级别达到了 $0.002$ 的宽度。
2. 2D $H_2^+$ 离子：双中心势的结果显示了收敛性，基态的认证区间随着网格细化从 $[-4.43, -1.317]$ 缩小到 $[-2.36, -1.317]$。
3. 3D 氢原子：在球域上，该方法将精确的 $\lambda_1 = -0.25$ 括定在内。认证的 Neumann 包围宽度以 $O(h^2)$ 收敛。
4. 3D $H_2^+$ 离子：对于盒域上的分子离子，该方法为基态提供了认证区间 $[-1.045, -0.544]$（宽度 $0.501$），该区间包含参考值$\approx -0.5513$。

在所有情况下，CECR 下界间隙 $(\mu_{k,h}^\sigma - \sigma) - L_{k}^{\mu,\sigma}$ 在分级网格上均以预测的 $O(h^2)$ 速率闭合。论文指出，虽然由于粗四面体网格上单元常数的衰减较慢，3D 间隙比 2D 更宽，但收敛行为与理论一致。

意义与声明
本文断言，这项工作构成了首个能够提供与无界奇异势相关的特征值的显式且收敛下界的有限元框架。其意义在于克服了由吸引奇异势引起的“强制性间隙”，此前这导致标准下界技术发散。

作者强调，该方法是显式的：所有常数（插值常数、扰动界、平移参数）均从网格和势数据中以闭式推导得出。虽然当前的实现使用双精度算术，但作者明确指出，转换为完全验证的（IEEE-754 安全、区间算术）包围是留给未来工作的单独验证步骤。论文还指出，虽然当前的界是基于投影的，但它们为更锐利的后处理技术（如 Lehmann–Goerisch 界）提供了严格的起点（$\rho$），这些技术可以进一步收紧区间。

这项工作并未声称直接解决多电子问题，但指出将方法扩展到磁 Hamilton 量和具有反对称约束的多电子系统是自然的未来方向。