コンテンツにスキップ

スパース対角化

局所的な相互作用を持つ系の量子力学におけるハミルトニアン行列はスパース(疎)であり、要素の大部分がゼロです。サイズ N×NN \times N の密行列を保存・操作するには O(N2)O(N^2) のメモリが必要であり、対角化には O(N3)O(N^3) の計算時間がかかります。NN が大きい場合、これは実行不可能になります。Lanczos 法はこのスパース性を活用するように設計されており、行列を明示的に保存・操作する必要がなく、行列ベクトル積だけで済みます。

Lanczos 法は、大規模なスパース行列について、少数の端の固有値とそれに対応する固有ベクトルを求めるための強力なアルゴリズムです。特に量子格子系に有用です。量子格子系は多くの場合、密行列の手法では扱えないほど大きな高次元の行列で記述されるためです。Lanczos 法はこれらの行列のスパース性を利用して、求めたい固有値と固有ベクトルを効率よく計算します。

ALPS におけるスパース対角化の実装についても説明します。