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

技術要約：シュレディンガー作用素への適用を伴うスペクトル・ガラーキン法における保証付き下界

問題提起
スペクトル・ガラーキン法は、特に滑らかな問題に対して超代数的な収束を示すことから、固有値の高精度な近似において広く認識されている。しかし、スペクトル・フレームワーク内で、これらの固有値に対する厳密かつ計算可能な下界を求めることは、依然として困難な課題であった。古典的な包絡法（Kato、Weinstein–Temple、Lehmann–Goerischなど）は下界を提供できるが、通常、隣接する固有値に関する事前情報や、特定の問題構造を必要とする。対照的に、有限要素法（FEM）に基づく既存の保証付き下界フレームワークは、局所的なメッシュ幾何学や補間誤差定数（$C_h$）に依存しており、これらはスペクトル手法に対応する直接的な対応物を持たない。本論文は、シュレディンガー作用素（$-\Delta + V$）に適用される適合スペクトル・ガラーキン法のための、厳密かつ無条件な下界フレームワークを確立することで、この空白を埋めるものである。

手法およびフレームワーク
本論文は、離散的なガラーキン固有値 $\lambda_{k,N}$ を、射影誤差定数 $C_N$ を介して厳密な固有値 $\lambda_k$ に関連付ける、抽象的なヒルベルト空間フレームワーク（Liu [16, 17] によって元々開発されたもの）に基づいている。核となる不等式は以下の通りである：
$$\lambda_k \ge \frac{\lambda_{k,N}}{1 + C_N^2 \lambda_{k,N}}$$
本手法は、スペクトル空間における $C_N$ を決定するために、4つの異なるメカニズムを通じて進行する。

1. 厳密な固有基底空間： 作用素の厳密な固有関数によって張られる試行空間の場合、本論文は、メッシュ幾何学に依存しないパルセーバルの議論を通じて、閉形式の最適な射影定数 $C_N = \lambda_{M+1}^{-1/2}$（ここで $\lambda_{M+1}$ は最初に見落とされた固有値）を導出する。
2. 射影ギャップ・メカニズム（標準行列）： 非負のポテンシャル（$V \ge 0$）を持つシュレディンガー作用素において、試行空間（例：三角多項式）が固有基底ではない場合、標準的なリッツ射影 $P_N$ はラプラシアンのモード射影 $P_N^0$ と異なる。本論文は、補助的な離散解作用素を用いて、差 $e = P_N^0 u - P_N u$ を制御する「射影ギャップ」解析を導入する。これにより、ポテンシャルのフォームファクター $\eta_V$ に依存する定数 $C_N^{(V)}$ が得られる。有界なポテンシャルの場合、この因子は当初 $\|V\|_{L^\infty}$ に依存し、$X = \|V\|_{L^\infty}/\lambda_{1,N}$ およびラプラシアンのスペクトルギャップ $\nu^*$ を用いて $C_N^{(V)} = \sqrt{(1 + \min(X^2, X))/\nu^*}$ となる。
3. 複合離散化（$\|V\|_{L^\infty}$ 依存性の除去）： 大きなポテンシャル（非有界領域のドメイン切断で一般的）による定数の膨張を排除するために、本論文は複合離散化を提案する。これには以下が含まれる：
   • ポテンシャル $V$ の帯域制限された下近似 $W$ の構成。
   • 部分空間 $V_m$ への次数低減された重み付き射影 $\Pi_m^{\tilde{W}}$ の定義。
   • 高次空間の補間残差を重み付き射影が消去するという「スレイビング恒等式（exact aliasing property）」の利用。
     この構成により、ポテンシャルの大きさに依存しない、鋭い定数 $C_N = \nu^{*-1/2}$ が回復され、結果として下界の質がポテンシャルの大きさからデカップリングされる。
4. ドメイン切断： $\mathbb{R}^d$ 上の作用素に対しては、ノイマン境界条件を伴う切断されたボックス $\Omega$ 上のスペクトル下界を、外部ポテンシャルの最小値 $\sigma(\Omega)$ と比較する事後的な閉じ込めチェックと組み合わせることで、全空間に対する証明付きの両側境界を提供する。

主な貢献

• 最適なスペクトル定数： 固有基底試行空間に対する閉形式の定数 $C_N = \lambda_{M+1}^{-1/2}$ の導出であり、これは理論的なベンチマークとなる。
• 射影ギャップ解析： 補助作用素の手法をスペクトル法へと拡張し、有界ポテンシャルを持つシュレディンガー作用素の標準的なガラーキン行列に対する明示的な定数を提供する。
• 複合スペクトル法： 大きなポテンシャル（閉じ込めポテンシャル）に対して品質が低下しない、ポテンシャルの $\|V\|_{L^\infty}$ に依存しない最適なスペクトルギャップ定数を実現する複合離散化の開発。
• 特異ポテンシャルへの拡張： ローカライズされたハーディ不等量を介して、補助射影メカニズムがどのように特異なポテンシャル（例：クーロン特異性）へ拡張されるかのプレビューを提供し、$L^\infty$ ノルムの必要性を排除する。
• 両側包絡： ノイマン・ドメイン切断を用いて、$\mathbb{R}^d$ 上のシュレディンガー作用素に対する完全に計算可能で証明付きの両側境界を提供する。

数値結果
本論文は、2次元における2つのベンチマークポテンシャル、環状ポテンシャル $V_1 = (|x|^2-1)^2$ と異方性ダブルウェル $V_2 = (x_1^2-1)^2 + x_2^2$ について理論を検証している。

• 効率性： スペクトル法は、[18] の証明付き Crouzeix–Raviart 拡張（CECR）有限要素法と比較して、約2桁少ない自由度（例：約400,000 DOFs に対し 約2,400 DOFs）で、同等またはそれを上回る証明付きの精度を達成する。
• 精度： 環状ポテンシャルにおいて、$N=48$ でのスペクトル包絡は、FEMの包絡の厳密に内側に位置する。ダブルウェルにおいては、スペクトル下界は $N=48$ でFEMの下界と同等である。
• メカニズムの比較： 中程度のポテンシャルの領域では、標準的な補助射影法（定理4.9）は効率的でほぼ最適である。しかし、大きな閉じ込めポテンシャル（$X \gg 1$）の場合、標準的な手法のギャップは著しく膨張する（例：複合メソッドの0.48%に対し、標準メソッドは72%のギャップ）。複合メソッド（定理4.12）は、ポテンシャルの大きさに関わらず、最適な収束率 $O(N^{-2})$ を維持する。
• テイル精緻化： 標準的な行列から計算可能な「テイル精緻化（tail-refined）」フォームファクター（定理4.5）は、標準的な手法の精度を大幅に改善することを示す（$V_1$ においてギャップを72%から7.7%に減少させる）が、非常に大きなポテンシャルに対しては複合メソッドには及ばない。

意義および主張
本論文は、事前情報を必要とせずに、スペクトル離散化から直接、厳密な下界を得るという長年の問題を解決したと主張している。スペクトル法に射影ベースのフレームワークを拡張することで、以下を実現している：

1. 無条件の境界： 隣接する固有値を必要としない境界。
2. 最適な定数： メッシュ依存の補間定数ではなく、ラプラシアンのスペクトルギャップによって決定される定数。
3. 高い効率性： 上界のスペクトル収束を利用することで、低次のFEMと比較して、大幅に少ない自由度で証明付きの精度を達成できる。
4. 堅牢性： 複合メソッドは、ドメイン切断を伴う量子化学への応用において極めて重要となる、大きなポテンシャルや閉じ込めポテンシャルに対しても、証明の質が低下しないことを保証する。

著者らは、上界はスペクトル収束する一方で、証明付きの下界は代数的な収束率 $O(N^{-2})$ で閉じることを指摘している。指数的な精度を下界へと転移させることは、未解決の課題として残されている。特異なポテンシャル（例：クーロン系）への拡張は、随伴論文の主題として認められている。