[論文レビュー] Quasidiagonality and the finite section method

Nathanial P. Brown|ArXiv.org|Dec 16, 2003

Spectral Theory in Mathematical Physics参考文献 13被引用数 4

ひとこと要約

この論文は、自己随伴、正規、および特定のバンド作用素などの準対角作用素が、有限区間法を用いて安定的かつ収束性を示す有限次元近似を有することを確立している。特に、周期的係数をもつ三重対角作用素に対して、ベルグの技法を用いて非零要素が制御されたほぼスパースな有限次元近似を構築し、誤差がブロックサイズに反比例して減少する明示的かつ計算可能な誤差境界を提供している。

ABSTRACT

Quasidiagonal operators on a Hilbert space are a large and important class (containing all self-adjoint operators for instance). They are also perfectly suited for study via the finite section method (a particular Galerkin method). Indeed, the very definition of quasidiagonality yields finite sections with good convergence properties. Moreover, simple operator theory techniques yield estimates on certain rates of convergence. In the case of quasidiagonal band operators both the finite sections and rates of convergence are explicitly given.

研究の動機と目的

• 準対角作用素が持つ構造的性質のおかげで、有限区間法に特に適していることを示すこと。
• 作用素論を用いて、準対角作用素の有限区間近似の一般収束推定を導出すること。
• バンド作用素、特に周期的係数をもつ三重対角作用素の場合に、明示的かつ計算可能な誤差境界を提供すること。
• ベルグの技法による有限次元近似の構築法が、誤差率を制御しながら準対角作用素に効果的に適用可能であることを示すこと。
• 理論的収束結果を、特に離散化ハミルトニアンやほぼマチュー作用素の実用的数値実装と結びつけること。

提案手法

• 準対角性の定義を用いる：すべての $v \in H$ に対して $\|P_n v - v\| \to 0$ かつ $\|[T, P_n]\| \to 0$ を満たす有限ランク射影 $P_n$ の存在。
• 有限区間法を適用し、$T$ を $P_n H$ に射影することで、有限次元作用素 $P_n T P_n$ を得る。
• 周期的係数の場合に、ねじれ基底を構築するベルグの技法を用い、非零要素が制御されたほぼ三重対角の有限区間近似を可能にする。
• 作用素ノルムの境界を用いて誤差推定を導出：$\|[T, P_k]\| \leq \frac{2\pi}{pk+1}$ ここで $p$ は周期、$k$ はブロックインデックス。
• スペクトル近似理論を用いて、$\sigma^{(\varepsilon)}(P_k T P_k)$ が $\sigma^{(\varepsilon)}(T)$ に収束することを示し、誤差境界が $\frac{1}{pk+1}$ に依存することを示す。
• 連分数近似を用いて、無理数周波数を有理数近似し、スペクトル計算における誤差制御を維持する。

実験結果

リサーチクエスチョン

• RQ1 有限区間法における準対角作用素の挙動はどのように規定され、どのような収束保証が得られるか？
• RQ2 準対角バンド作用素の有限区間近似に対して、明示的誤差境界をどのように導出できるか？
• RQ3 ベルグの基底構築法を、スパarsityを保ちながら周期的係数の場合に誤差率を制御可能な形に適応可能か？
• RQ4 作用素 $T$ のスペクトルが $P_k T P_k$ によってどのように近似されるか、およびスペクトルの収束速度はどの程度か？
• RQ5 連分数を用いて、非周期的だがほぼ周期的係数をもつ場合に、誤差を制御可能な周期的近似を構築可能か？

主な発見

• 準対角作用素に対しては、有限区間法により収束近似が得られ、$\|x - x_k\| \leq \|T^{-1}\|(\|y - P_k y\| + \frac{2\pi \|x_k\|}{pk+1})$ を満たす。ここで $x_k$ は $P_k T P_k x_k = P_k y$ を満たす。
• $\varepsilon$-擬スペクトルは $\sigma^{(\varepsilon)}(P_k T P_k) \subset \sigma^{(\varepsilon + \frac{2\pi}{pk+1})}(T)$ に収束する。
• 2-ノルム擬スペクトルは $\sigma_2(P_k T P_k) \subset^{\sqrt{\frac{4\pi \|T\|}{pk+1}}}} \sigma_2(T)$ を満たし、誤差がブロックサイズに反比例することを示唆する。
• 自己随伴 $T$ の場合、スペクトルは $\sigma(P_k T P_k) \subset^{\frac{2\pi}{pk+1}} \sigma(T)$ に収束し、誤差境界がブロックサイズに反比例する。
• 周期的係数をもつ三重対角作用素に対して、有限区間法は有効である。ここで $P_k T P_k$ の行列は、明示的な要素をもつほぼ三重対角のままである。
• 連分数を用いることで、無理数周波数 $\theta$ に対して $|\theta - \frac{p_n}{q_n}| \leq \frac{1}{q_n^2}$ を満たす有理数近似が得られ、誤差制御可能なスペクトル近似が可能になる。