物理学における数値線形代数と固有値問題

物理演算子を離散化すると、物理学は行列に変換され、系のエネルギーとモードを見つけることは、大規模な線形システムを解き、固有値と固有ベクトルを計算するという数値問題になります。

Definition

物理学における数値線形代数とは、連続的な物理演算子が有限基底またはグリッド上で表現される際に生じる行列方程式および固有値問題を解くための一連のアルゴリズムを指します。

Scope

このトピックでは、物理学の中心となる行列計算について扱います。具体的には、直接法および反復法による線形システムの解法、ならびにQR法、ヤコビ法、ランチョス法、共役勾配法を用いた、大規模でしばしば疎なエルミート行列の固有値および固有ベクトルの計算です。疎性やエルミート性といった物理行列の構造に重点を置いています。

Core questions

• 離散化された物理学から生じる大規模な線形システムは、密な逆行列を形成せずにどのように解かれるのでしょうか？
• ハミルトニアン行列の固有値と固有ベクトルは、数値的にどのように計算されるのでしょうか？
• 大規模な疎行列に対して、直接因子分解よりも反復的なクリロフ法が好まれるのはなぜでしょうか？
• ランチョスアルゴリズムは、巨大な疎なエルミート行列のいくつかの極端な固有値をどのように抽出するのでしょうか？

Key theories

直接線形ソルバーと反復線形ソルバー

線形システムは、LU分解やコレスキー分解のような直接因子分解によって、丸め誤差を除いて厳密に解かれるか、または共役勾配法のような反復的なクリロフ法によって、疎性を利用し許容誤差まで収束させて解かれます。

固有値アルゴリズム

固有値と固有ベクトルは、密行列に対してQRアルゴリズムやヤコビ回転によって計算され、有限基底で表現された物理演算子の離散スペクトルを与えます。

ランチョス法とクリロフ部分空間法

ランチョスアルゴリズムは、クリロフ部分空間において大規模な疎なエルミート行列の小さな三重対角射影を構築し、完全な行列を記憶することなく、いくつかの極端な固有値と固有ベクトルを見つけることを可能にします。

Clinical relevance

これらのアルゴリズムは、量子力学におけるエネルギー準位と波動関数、振動の基準モード、固体中のバンド構造、および離散化された場の方程式の背後にある線形システムを計算するため、電子構造計算や凝縮系シミュレーションにおいて不可欠なものとなっています。

History

実用的な行列固有値計算は、1950年のランチョス反復法と1960年代初頭のQRアルゴリズムの登場により20世紀半ばに成熟しました。物理学における大規模な疎行列問題の増加により、クリロフ部分空間法が高次元ハミルトニアンのスペクトル解析における主要なツールとなりました。

Key figures

• Cornelius Lanczos
• Gene H. Golub
• James H. Wilkinson

Seminal works

• golub2013
• lanczos1950

Frequently asked questions

なぜ行列全体を対角化するのではなく、反復法を使用するのですか？

物理的なハミルトニアンは数十億の次元を持つことがありますが、疎であるため、それらを完全に記憶したり因子分解したりすることは不可能です。ランチョス法のような反復的なクリロフ法は、行列とベクトルの作用のみを必要とし、物理学で通常関心のあるいくつかの最低エネルギー固有状態を抽出することができます。

物理行列のエルミート性が数値的に重要なのはなぜですか？

エルミート行列は実数の固有値と直交する固有ベクトルを持つため、より安定で効率的な特殊なアルゴリズムを使用でき、計算されたエネルギーが実数であることを保証し、物理学と一致します。