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

技術的概要：有限要素法による特異ポテンシャルを伴うシュレーディンガー作用素の明示的な両側固有値 bound

問題の定義
本論文は、非有界領域 $\mathbb{R}^N$（$N=2, 3$）上のシュレーディンガー作用素 $H = -\Delta + V$ の固有値 $\lambda_k$ に対する、明示的かつ計算可能な両側 bound の計算に取り組む。主な課題は、水素原子をモデル化する $V(x) = -Z/|x|$ および $H_2^+$ 分子イオンをモデル化する $V(x) = -Z_1/|x-a_1| - Z_2/|x-a_2|$ という、特異な引力性クーロン型ポテンシャルの扱いにある。

既存の数値計算手法は、通常、適合ガレルキン（レイリー・リッツ）法を通じて上限 bound を提供する。一方、厳密な下限 bound を得ることは、特にポテンシャルが下方に無界である非有界領域および特異ポテンシャルにおいて、著しく困難である。古典的な下限 bound 手法（Lehmann–Maehly–Goerisch、Kato–Temple など）は、非適合法や特定の条件を必要とし、ポテンシャル $V$ が負で特異である場合には機能しない。さらに、下限 bound に対する標準的な有限要素アプローチは、しばしばポテンシャルの負の部分の $L^\infty$ ノルム（$\gamma_h = \|c_h^-\|_{L^\infty}$）によるシフトに依存する。クーロン特異点の場合、$\gamma_h$ はメッシュが細かくなるにつれて発散（$O(h^{-1})$）し、下限 bound の収束を阻害する。

手法
提案されたアルゴリズムは、領域の切断と、複合強化 Crouzeix–Raviart（CECR）有限要素法の拡張を組み合わせる。このアプローチは、以下の 3 つの主要な構成要素からなる。

1. 領域の切断：非有界領域 $\mathbb{R}^N$ を有界領域 $D(R)$ に切断する。

   • 上限 bound：標準的な適合 $P_1$ ガレルキン要素を用い、ディリクレ境界条件（$\partial D(R)$ 上で $u=0$）を通じて得られる。
   • 下限 bound：ノイマン境界条件（$\partial u/\partial n = 0$）を通じて得られる。重要な補題により、外部ポテンシャルの最小値 $\sigma_{ext}(D(R))$ が $k$ 番目の固有値を超えている場合、$D(R)$ 上のノイマン固有値が完全空間の固有値に対する有効な下限 bound となることが示される。この条件は、計算された上限 bound を用いて事後検証される。

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)$、メッシュの細分化に依存）する。
   • ずらされた双線形形式の強制性を保証するため、メッシュサイズ $h$ に依存しない固定シフト $\sigma > 0$ をポテンシャルに適用する。これにより、発散するメッシュ依存シフトの必要性を回避する。

3. 摂動解析：

   • 連続的なポテンシャル $V$ は、その要素ごとの定数平均 $c_h$ によって近似される。
   • 摂動補題により、$V$ を $c_h$ で置き換えることで生じる誤差が bound される。摂動比 $\epsilon_h$ は、特異点近傍における特定のメッシュ細分化の仮定の下で、ゼロに収束することが示される（細分化メッシュ上で $O(h^2)$）。

主要な貢献
本論文は、以下の具体的な貢献を主張する。

• 最初の明示的収束下限 bound：非有界領域上の非有界特異ポテンシャルを伴うシュレーディンガー作用素に対して、明示的かつ収束する両側固有値 bound を提供する数値アルゴリズムを初めて提示する。
• ペア空間抽象定理：CECR 解析を 2 成分空間に持ち上げるペア空間フレームワークを開発し、符号不定ポテンシャルに対する下限 bound 定理の自己完結的な max-min 証明を可能にする。
• 要素ごとの反応定数（$\Gamma_h$）：発散する $L^\infty$ シフトに代わる $\Gamma_h$ を導入する。これにより、CECR のギャップが細分化メッシュ上で $O(h^2)$ の速度で閉じ、計算不可能な形式ノルム定数の必要性を回避する。
• 固定シフト両側 bound：ポテンシャルの負の部分の $L^\infty$ 制御を必要としない、完全に計算可能な囲み $L_k \le \lambda_k \le U_k$ を構築する。
• $H_2^+$ への適用：2 中心クーロンポテンシャルに理論を特化し、2 次元および 3 次元の両方において許容されるシフトと摂動パラメータの明示的な式を提供する。

数値結果
本アルゴリズムは、4 つのテストケースで検証された。

1. 2 次元水素原子（$V = -1/|x|$）：本手法は、正確な基底状態 $\lambda_1 = -1$ を見事に囲い込んだ。認証された下限 bound 区間はメッシュの細分化に伴い著しく狭まり、最適シフト bound は最細レベルで幅 $0.002$ を達成した。
2. 2 次元 $H_2^+$ イオン：2 中心ポテンシャルの結果は収束を示し、基底状態の認証された区間はメッシュが細かくなるにつれて $[-4.43, -1.317]$ から $[-2.36, -1.317]$ に狭まった。
3. 3 次元水素原子：球状領域において、本手法は正確な $\lambda_1 = -0.25$ を囲い込んだ。認証されたノイマン囲いの幅は $O(h^2)$ で収束した。
4. 3 次元 $H_2^+$ イオン：箱状領域上の分子イオンに対して、本手法は基底状態の認証された区間 $[-1.045, -0.544]$（幅 $0.501$）を提供し、これは参照値$\approx -0.5513$ を含んでいる。

すべてのケースにおいて、CECR 下限 bound ギャップ $(\mu_{k,h}^\sigma - \sigma) - L_{k}^{\mu,\sigma}$ は、細分化メッシュ上で予測された $O(h^2)$ 速度で閉じることが観察された。本論文は、粗い四面体メッシュ上の要素ごとの定数の減衰が遅いため、3 次元のギャップは 2 次元よりも広いものの、収束挙動は理論と一致していると指摘している。

意義と主張
本論文は、この仕事が非有界特異ポテンシャルに関連する固有値に対して明示的かつ収束する下限 bound を提供できる最初の有限要素フレームワークであると主張する。その意義は、引力性特異ポテンシャルによって引き起こされる「強制性のギャップ」を克服し、それまで標準的な下限 bound 手法が発散していた問題を解決した点にある。

著者らは、本手法が明示的であることを強調している。すべての定数（補間定数、摂動 bound、シフトパラメータ）は、メッシュとポテンシャルデータから閉じた形式で導出される。現在の実装は倍精度算術を使用しているが、著者らは、完全に検証された（IEEE-754 安全、区間算術）囲みへの変換は、将来の作業に委ねられた別個の検証ステップであると明確に述べている。また、本論文は、現在の bound が射影ベースであるものの、区間をさらに狭めることができる Lehmann–Goerisch bound のような鋭い後処理手法に対する厳密な出発点（$\rho$）として機能すると指摘している。

本作業は、多電子問題を直接解決するものではないが、磁気ハミルトニアンおよび反対称性の制約を伴う多電子系への拡張を、自然な将来の方向性として特定している。