密度行列繰り込み群法 – Wikipedia

出典は列挙するだけでなく、脚注などを用いてどの記述の情報源であるかを明記してください。記事の信頼性向上にご協力をお願いいたします。（2016年10月）

密度行列繰り込み群法（みつどぎょうれつくりこみぐんほう 英: density matrix renormalization group; DMRG）は、量子多体系における低エネルギー物理を高精度に計算するために考案された数値変分法である。1992年に Steven R. White により開発された。

DMRG の背景にある考え方[編集]

量子多体系の物理に関して主に問題となるのは、ヒルベルト空間が指数関数的に大きくなることである。例えば、長さ L のスピン 1/2（英語版） チェインは、2^L の自由度をもつ。DMRG法は反復的な変分法であり、問題の量子状態についてもっとも重要な自由度にのみ有効自由度を絞り込むことができる。問題とされるのは基底状態であることが多い。

この手法では、ウォームアップサイクル後に系を（同じサイズとは限らない）二つのブロックと、その間に位置するの二つのサイトに分ける。ウォームアップ中に、各ブロックを「代表する」一連の状態を選定する。左ブロック + 二つのサイト + 右ブロックを合わせてスーパーブロックと呼ばれる。スーパーブロックは全系よりも自由度が低減しており、基底状態の候補が見付けやすい。その代償として精度は低下するが、下記の反復法により向上させることができる。

見付かった基底状態の候補を、名前の通り密度行列を用いて各ブロックに対応する部分空間上に射影する。これにより、各ブロックの「関連する状態」が更新される。

ここで、片方のブロックを大きくし、もう片方を小さくして同じ手続きを繰り返す。大きくしたブロックが最大サイズに到達したら、かわりにもう片方を大きくする。最初の（等しいサイズの）状況に立ち戻ったとき、「掃引」が完了したという。1 次元格子ならば通常、数回の掃引で 10¹⁰ 分の 1 の精度を得るのに十分である。

DMRG法は Steven White と Reinhard Noack により、1 次元箱内のスピン 0 粒子のスペクトルを求めるというトイモデルに対して始めて適用された。このモデルはケネス・ウィルソンにより、何らかの新しいくりこみ群の方法をテストするために考案された。このような単純な問題でも、正しく解けない方法ばかりだったのである。DMRG法は従来のくりこみ群の方法にあった問題点を、系を一つのブロックと一つのサイトに分けるのではなく二つのブロックを二つのサイトで繋ぐように分け、さらに各ステップの最後に最も重要で保存するべき状態を密度行列を用いて識別することにより克服している。 このトイモデルを解くことに成功したのち、DRMG法はハイゼンベルクモデル（英語版）にも適用され、成功している。

実装上の技術的詳細[編集]

DMRGアルゴリズムの実用的実装は長大な作業である。次のような計算上のトリックが主に用いられている。

• スーパーブロックの基底状態はランチョス法による行列対角化により求める。他にも、特に非エルミート行列を扱う場合はアーノルディ法（英語版）を使うこともある。
• ランチョス法はランダムなシードから開始することが多いが、DMRG法ではあるステップで得られた解を適切に変換してシードとする法が良い場合がある。
• 対称性のある系の場合、例えばハイゼンベルクモデルにおける総スピンなど、保存される量子数がある場合がある。これによりヒルベルト空間を分割し、そのそれぞれについて基底状態を求めるのが便利である。
• 例として、ハイゼンベルクモデルに対するDMRG法（英語版）が挙げられる。

DMRG法は、横磁場イジングモデルやハイゼンベルクモデル（英語版）など、およびハバードモデルなどのフェルミオン系、近藤効果などの欠陥のある問題、ボソン系、量子ワイヤー（英語版）に接続された量子ドットの物理など、スピンチェインの低エネルギー物性を得るための応用が成功している。樹状グラフを扱えるよう拡張されたものもあり、デンドリマーの研究に応用されている。片方の次元がもう片方よりも非常に大きいような二次元系も精度よく扱えるため、ラダーの研究にも有用であることが知られている。

二次元系の平衡状態についての統計物理学（英語版）^{[要リンク修正]}的研究向けや、一次元系の非平衡（英語版）現象の解析向けの拡張も存在する。

量子化学分野においては強相関系を扱うための応用もされている。

行列積仮設[編集]

DRMG法が一次元系で成功した背景には、これが行列積状態空間上における変分法であるという事実がある。行列積状態とは、次の形式で表わされる状態である。

${displaystyle sum _{s_{1}cdots s_{N}}operatorname {Tr} (A^{s_{1}}cdots A^{s_{N}})|s_{1}cdots s_{N}rangle }$

ここで、s₁ … s_N は例えばスピンチェイン上のスピンの z 成分であり、 A^s_i は任意の m 次元行列である。m → ∞ の極限において、この表現は厳密となる。このことは S. Rommer と S. Ostlund により理論化された。

DMRG の拡張[編集]

2004年、行列積状態の実時間発展向けに時間発展ブロックデシメーション法（英語版）が実装された。このアイデアは量子コンピュータの古典シミュレーションに基いている。続いて、DRMG形式の実時間発展を計算する新手法が考案された。これについては A. Feiguin と S.R. White による論文を参照のこと。

近年、行列積状態の定義を拡張することにより、二次元および三次元へと拡張する提案がなされている。これについては F. Verstraete と I. Cirac による論文を参照のこと。

参考文献[編集]

関連ソフトウェア[編集]

• Powder with Power: Fortran で書かれた時間依存DMRG法用フリープログラム
• The ALPS Project: C++ で書かれた時間非依存DMRGおよび量子モンテカルロ用フリープログラム
• DMRG++: C++ で書かれたフリーなDMRG実装
• ITensor (Intelligent Tensor) ライブラリ: C++ で書かれたテンソル積状態および行列積状態に基くDMRG計算用のライブラリ
• Snake DMRG program: C++ で書かれたDMRG, tDMRG および有限温度DMRG用プログラム
• CheMPS2: C++ で書かれた非経験的分子軌道法向けスピン適合DRMG用オープンソースプログラム
• Block: C++ で書かれた量子化学およびモデルハミルトニアン向けのDMRGオープンソースフレームワーク。SU(2) および一般的な非アーベル対称性をサポート。

関連項目[編集]