[論文レビュー] On Temple--Kato like inequalities and applications

Luka Grubišić|ArXiv.org|Nov 16, 2005

Spectral Theory in Mathematical Physics参考文献 33被引用数 5

ひとこと要約

本稿では、シュール補行列や二次形式などの行列理論的道具を用いて、ヒルベルト空間内に定義された非有界な正定値作用素の離散固有値について、スケーリングに頑健で、最適な上界と下界を確立する。主な貢献は、テンプル–カトー不等式の二次形式版を提案することであり、これにより、特異的摂動を受けるシュレーディンガー作用素や、収束性と最適性が保証された適応型有限要素法に適用可能な、高性能な残差型固有値推定が可能になる。

ABSTRACT

We give both lower and upper estimates for eigenvalues of unbounded positive definite operators in an arbitrary Hilbert space. We show scaling robust relative eigenvalue estimates for these operators in analogy to such estimates of current interest in Numerical Linear Algebra. Only simple matrix theoretic tools like Schur complements have been used. As prototypes for the strength of our method we discuss a singularly perturbed Schroedinger operator and study convergence estimates for finite element approximations. The estimates can be viewed as a natural quadratic form version of the celebrated Temple--Kato inequality.

研究の動機と目的

• ヒルベルト空間内に定義された非有界な正定値作用素に対して、スペクトルギャップや特異的摂動が存在する状況でも、頑健で最適な固有値推定を構築すること。
• 作用素の構造に制限を課さない形で、行列に基づく相対的摂動理論を無限次元設定に拡張するための二次形式を用いたアプローチ。
• メッシュのスケーリングに依存しない、明示的かつ計算可能な誤差境界を、有限要素近似に対して提供すること。
• 適応型有限要素法における固有値の重複度と残差に基づく誤差指標を自然に扱えるフレームワークを確立すること。
• 剛性／特異的摂動作用素や大結合定数極限の文脈において、推定の最適性を実証すること。

提案手法

• カトーやリード＝シンモンの二次形式理論を用いて、非有界な正定値作用素の正則化フレームワークを定式化する。
• 特にシュール補行列とブロック作用素分解を用いた行列解析的手法を適用し、有限次元相対摂動理論の結果に類似した残差方程式を導出する。
• 直交射影を用いたブロック作用素系を構築し、ガラーキン型離散化により作用素の逆行列を近似する。
• 残差の $‘H^{-1}$ ノルムにおける誤差推定を導出し、これは形状正則性にのみ依存する定数を除いて、実際の固有値誤差と等価であることが示される。
• ユニタリ不変ノルムと $‘H_{2,d}^{-1}$ 内の残差ノルムを用いて、完全な線形系を解かずに効率的かつ低コストな誤差推定を可能にする。
• 既知の結果 [2, 9, 23] を用いたガラーキン直交性条件と組み合わせ、トレース型およびノルムに基づく誤差境界を導出する。

実験結果

リサーチクエスチョン

• RQ1 スケーリングに頑健な固有値推定を、非有界な正定値作用素のヒルベルト空間内に、基本的な行列ツールのみを用いて導出可能か？
• RQ2 テンプル–カトー不等式を、最適性を保ちつつスペクトルギャップを扱える二次形式設定に一般化できるか？
• RQ3 得られた固有値誤差推定は、特異的摂動シュレーディンガー作用素および有限要素近似に対して、どの程度最適か？
• RQ4 完全な線形系を解かずに、$H^{-1}$-類似計量における残差ノルムを効率的に近似可能か？
• RQ5 本手法は、適応型有限要素法の文脈において、固有値の重複度とクラスタ解像度をどのように処理するか？

主な発見

• 本稿では、シュール補行列技術を用いて、ヒルベルト空間内に定義された正定値作用素の離散固有値について、スケーリングに頑健で最適な上界と下界を導出する。
• 剛性／特異的摂動作用素の文脈において、最低次の固有値に対して推定が最適であることが示され、モデルシュレーディンガー作用素を用いた数値的検証が行われた。
• $‘H_{2,d}^{-1}$ ノルムにおける残差に基づく誤差推定は、メッシュの形状正則性にのみ依存する定数を除いて、真の固有値誤差と等価であることが証明された。
• 誤差境界 $c_* \frac{\|r_{i,d}\|^2_{\mathbf{H}_{2,d}^{-1}}}{\mu_{i,d}} \leq \frac{\|\mathbf{H}u_{i,d}-\mu_{i,d}u_{i,d}\|^2_{\mathbf{H}^{-1}}}{\|\mathbf{H}u_{i,d}\|^2_{\mathbf{H}^{-1}}} \leq C_* \frac{\|r_{i,d}\|^2_{\mathbf{H}_{2,d}^{-1}}}{\mu_{i,d}}$ が、メッシュの形状正則性にのみ依存する定数 $c_*$, $C_*$ を用いて成り立つ。
• 補空間内の補助ベクトルを用いることで、線形系を解くよりも安価な、$ \mathbf{H}^{-1}$-ノルムの残差の効率的かつ低コストな推定が可能になる。
• このフレームワークは固有値の重複度を自然に扱い、同一の残差に基づくアプローチにより、最適な固有ベクトル誤差推定が可能になる。