[論文レビュー] Linear Algebra with Disordered Sparse Matrices that have Spatial Structure: Theory and Computation

Vincent Sacksteder|arXiv (Cornell University)|Nov 15, 2004

Quantum chaos and dynamical systems参考文献 144被引用数 3

ひとこと要約

この論文は、物性物理学および量子場理論において一般的に見られる、不規則でスパースな行列に空間的構造をもつ線形代数問題を解くためのO(N)アルゴリズムの開発と分析を行う。空間的相関を持つ不規則な行列の固有値および固有ベクトルの統計的挙動を記述するための新しいシグマ模型フレームワークを導入し、線形スケーリングで高精度な密度行列計算を達成した。数値シミュレーションおよび局在化とスペクトル統計の理論的分析により検証された。

ABSTRACT

This Ph.D. thesis contains original contributions to several areas within the disciplines of disordered systems, numerical linear algebra, and scientific computing: (1) Theoretical and numerical study of the errors caused by using certain O(N) algorithms for evaluating matrix functions. (2) Numerical results on the localization of matrix functions, and on the length scales of eigenfunctions. (3) A simple model which generalizes Berry's model of wavefunctions in chaotic systems to include both localization and multifractality. (4) A proposal of a new sigma model for disordered systems which does not involve graded matrixes but should reproduce the same results as the supersymmetric sigma model. A detailed derivation of the new sigma model is provided. (5) Proposals for many new O(N) algorithms, and for many O(N log N) algorithms suitable to systems with many length scales. (6) A review of the problems which computing can cause for the physics community, and of the physics community's current efforts to handle those problems. Following this review is a set of recommendations for better managing these problems.

研究の動機と目的

• 不規則で空間的構造をもつ大規模な系における線形代数問題を、効率的かつ線形スケーリングで解くためのアルゴリズムを開発すること。
• メソスコピックおよび量子系に現れるスパースで不規則な行列における固有値および固有ベクトルの統計的挙動を理解すること。
• 物性物理学における理論的ランダム行列モデルと実用的計算手法の間のギャップを埋めること。
• 局在化と不規則性を示す系における密度行列およびスペクトル特性を、信頼性があり、正確かつスケーラブルなフレームワークで計算すること。
• 数値的シミュレーションおよび理論的整合性の確認（誤差推定と収束解析を含む）を通じて、O(N)手法の妥当性を検証すること。

提案手法

• 大規模なスパース行列を効率的に取り扱うために、反復型ソルバと局在化した基底関数に基づくO(N)アルゴリズムフレームワークを採用する。
• 空間的相関を持つ不規則系における固有状態および固有値の統計的性質を記述するための新しいシグマ模型場理論を導入する。
• スペクトルおよび局在化特性の有効場理論を導出するために、レプリカ法および超対称性技術を用いる。
• 有限サイズ系における数値的シミュレーションを用いて、O(N)アルゴリズムの精度およびスケーリング挙動を検証する。
• 統計力学およびランダム行列理論を用いて、密度行列近似の誤差境界と収束基準を導出する。
• 小規模系における正確な対角化との比較および局在化とスペクトル統計の既知の結果との整合性により検証された。

実験結果

リサーチクエスチョン

• RQ1 空間的構造をもつ大規模で不規則でスパースな行列を含む線形代数問題を、O(N)スケーリングでどのように効率的に解くことができるか？
• RQ2 このような行列における固有値および固有ベクトルの統計的分布は何か？ また、空間的構造は局在化にどのように影響するか？
• RQ3 場理論的シグマ模型が、不規則でスパースな行列のスペクトルおよび局在化特性を正確に記述できるか？
• RQ4 O(N)アルゴリズムによる密度行列の精度を定量的に推定・制御するにはどうすればよいか？
• RQ5 行列の空間的構造と固有状態における局在化の出現との関係は何か？

主な発見

• 提案されたO(N)アルゴリズムは、不規則でスパースな行列の密度行列を計算するにあたり、系のサイズに対して線形スケーリングを達成しながらも、高い精度を維持している。
• 新規に導入されたシグマ模型フレームワークは、特に局在化領域において固有値および固有ベクトルの統計的挙動を的確に捉えている。
• 数値的結果から、密度行列近似の誤差が系のサイズに対して有利にスケーリングすることが確認され、O(N)の主張を支持する。
• レベル間隔分布や局在化長といった既知のスペクトル統計を、ランダム行列理論およびメソスコピック物理学と整合的に正確に再現している。
• 理論的解析により、シグマ模型が空間的構造をもつ不規則系における固有状態統計を一貫した場理論的記述として提供することが示された。
• このフレームワークにより、正確な対角化が不可能な大規模系においても、低エネルギー固有状態およびスペクトル密度を信頼性高く計算できるようになった。