Guaranteed Lower Eigenvalue Bounds for Spectral Galerkin Methods with Application to Schrödinger Operators

技术摘要：谱伽辽金方法中保证的特征值下界及其在薛定谔算子中的应用

问题陈述
谱伽辽金方法因其在处理光滑问题时表现出的谱收敛（超代数收敛）特性，被广泛认为是在高精度近似特征值方面的卓越工具。然而，在谱框架内，建立一个严格且可计算的特征值下界仍然是一个难题。虽然经典的包络方法（如 Kato、Weinstein–Temple、Lehmann–Goerisch）可以提供下界，但它们通常需要关于相邻特征值的先验信息或特定的问题结构。相反，现有的基于有限元方法（FEM）的保证下界框架依赖于局部网格几何形状和插值误差常数（$C_h$），而这些常数在谱方法中没有直接的对应物。本文通过为应用于薛定谔算子（$-\Delta + V$）的保形谱伽辽金方法建立一个严格且无条件的下界框架，填补了这一空白。

方法论与框架
本文基于一个抽象的希尔伯特空间框架（最初由 Liu [16, 17] 开发），该框架将离散伽辽金特征值 $\lambda_{k,N}$ 与精确特征值 $\lambda_k$ 通过投影误差常数 $C_N$ 联系起来。核心不等式为：
$$\lambda_k \ge \frac{\lambda_{k,N}}{1 + C_N^2 \lambda_{k,N}}$$
该方法通过四种机制来确定谱空间的 $C_N$：

1. 精确特征基空间：对于由算子精确特征函数张成的试函数空间，本文推导出了闭合形式的最优投影常数 $C_N = \lambda_{M+1}^{-1/2}$，其中 $\lambda_{M+1}$ 是第一个被忽略的特征值。这是通过帕塞瓦尔（Parseval）论证建立的，且独立于网格几何形状。
2. 投影间隙机制（标准矩阵）：对于试函数空间（如三角多项式）不是特征基的非负势函数（$V \ge 0$）薛定谔算子，标准的 Ritz 投影 $P_N$ 与拉普拉斯模投影 $P_N^0$ 不同。本文引入了“投影间隙”分析，利用辅助离散解算子控制差值 $e = P_N^0 u - P_N u$。这产生了一个依赖于势函数形式因子 $\eta_V$ 的常数 $C_N^{(V)}$。对于有界势函数，该因子最初取决于 $\|V\|_{L^\infty}$，导致常数 $C_N^{(V)} = \sqrt{(1 + \min(X^2, X))/\nu^*}$，其中 $X = \|V\|_{L^\infty}/\lambda_{1,N}$，$\nu^*$ 是拉普拉斯算子的谱间隙。
3. 复合离散化（消除对 $\|V\|_{L^\infty}$ 的依赖）：为了消除由大势函数引起的常数膨胀（这在处理无界区域的域截断时很常见），本文提出了一种复合离散化方案，包括：
   • 构建势函数 $V$ 的带限欠近似 $W$。
   • 在子空间 $V_m$ 上定义降低阶数的加权投影 $\Pi_m^{\tilde{W}}$。
   • 利用一种“从属恒等式”（一种精确的混叠性质），其中加权投影消除了高阶空间的插值残差。
     该构造恢复了最优的、与势函数无关的常数 $C_N = \nu^{*-1/2}$，从而有效地将界限质量与势函数的大小解耦。
4. 域截断：对于 $\mathbb{R}^d$ 上的算子，本文结合了带有 Neumann 边界条件的截断框 $\Omega$ 上的谱下界，以及后验限制检查（通过将计算出的上界与外部势能极小值 $\sigma(\Omega)$ 进行比较），从而产生全空间的认证双侧界限。

主要贡献

• 最优谱常数：为特征基试函数空间推导出了闭合形式常数 $C_N = \lambda_{M+1}^{-1/2}$，这作为理论基准。
• 投影间隙分析：将辅助算子技术扩展到谱方法，为具有有界势函数的薛定谔算子的标准伽辽金矩阵提供了显式常数。
• 复合谱方法：开发了一种复合离散化方法，可以在不依赖于 $\|V\|_{L^\infty}$ 的情况下实现最优的谱间隙常数，克服了标准矩阵方法在处理束缚势时的局限性。
• 向奇异势函数的扩展：展示了辅助投影机制如何通过局部 Hardy 不等式扩展到奇异势函数（如库仑奇点），从而消除了对有界 $L^\infty$ 范数的要求。
• 双侧包络：为 $\mathbb{R}^d$ 上的薛定谔算子提供了完全可计算、经过认证的双侧界限，使用了 Neumann 域截断技术。

数值结果
本文在两个二维基准势函数上验证了理论：环形势 $V_1 = (|x|^2-1)^2$ 和各向异性双阱势 $V_2 = (x_1^2-1)^2 + x_2^2$。

• 效率：谱方法实现的认证界限在自由度（DOF）方面优于或等同于经过认证的 Crouzeix–Raviart 增强（CECR）有限元方法（来自 [18]），且使用的自由度减少了约两个数量级（例如，~2,400 DOF 对比 ~400,000 DOF）。
• 精度：对于环形势，当 $N=48$ 时，谱包络严格位于 FEM 包络内部。对于双阱势，在 $N=48$ 时，谱界限与 FEM 界限相当。
• 机制比较：在势函数适中的情形下，标准辅助投影法（定理 4.9）是高效且近乎最优的。然而，对于大束缚势（即 $X \gg 1$ 时），标准方法的间隙会显著膨胀（例如，72% 的间隙对比复合方法的 0.48%）。复合方法（定理 4.12）无论势函数大小如何，都能保持最优收敛率 $O(N^{-2})$。
• 尾部细化：本文展示了“尾部细化”形式因子（定理 4.5）——该因子可从标准矩阵中计算得出——能显著改善标准方法的界限（将 $V_1$ 的间隙从 72% 降至 7.7%），尽管对于极大的势函数，它仍逊于复合方法。

意义与主张
本文声称解决了直接从谱离散化中获得严格下界且无需先验谱信息的长期问题。通过将基于投影的框架扩展到谱方法，本文提供了：

1. 无条件界限：不需要了解相邻特征值的知识。
2. 最优常数：常数由拉普拉斯算子的谱间隙决定，而非取决于网格依赖的插值常数。
3. 高效率：能够利用上界的谱收敛特性，以显著更少的自由度实现认证精度。
4. 鲁棒性：复合方法确保了认证质量不会因势函数增大或由于域截断导致的束缚势而退化，这是涉及量子化学应用的至关重要要求。

作者指出，虽然上界呈谱收敛，但认证下界以代数速率 $O(N^{-2})$ 收敛。将指数精度转移到下界仍是一个开放性问题。对于奇异势函数（如库仑系统）的扩展被列为后续论文的主题。