[論文レビュー] Renormalization algorithms for Quantum-Many Body Systems in two and higher dimensions

Frank Verstraete, J. I. Cirac|arXiv (Cornell University)|Jul 2, 2004

Quantum many-body systems被引用数 503

ひとこと要約

本稿では、2次元以上への行列積状態（MPS）の自然な拡張として、投影もつれペア状態（PEPS）を導入し、量子多体系の効率的バリエーション的シミュレーションを可能にする。著者らは、相関関数の計算に適したスケーラブルなアルゴリズムを開発し、基底状態および虚時間発展演算に応用した。小さな結合次元Dでも高い精度が得られ、2次元ヘイゼンベルグ模型やフラストレーテッドスピン系において実証された。

ABSTRACT

We describe quantum many--body systems in terms of projected entangled--pair states, which naturally extend matrix product states to two and more dimensions. We present an algorithm to determine correlation functions in an efficient way. We use this result to build powerful numerical simulation techniques to describe the ground state, finite temperature, and evolution of spin systems in two and higher dimensions.

研究の動機と目的

• 2次元以上における量子多体系のシミュレーションのためのスケーラブルな数値的手法を開発すること。従来のDMRGやモンテカルロ法では限界がある。
• 補助系のもつれペアを用いて行列積状態（MPS）を高次元に一般化し、投影もつれペア状態（PEPS）の概念を導くこと。
• バリエーション的手法や時間発展演算シミュレーションに不可欠な、PEPSにおける相関関数の効率的計算を可能にすること。
• PEPSに基づくアルゴリズムが、スピン系の基底状態の特定および実時間・虚時間発展演算（特にフラストレーションを持つモデルを含む）に正確に応用可能であることを示すこと。
• さまざまな格子構造、有限温度状態、散乱を伴うオープン量子系に適用可能なフレームワークを確立すること。

提案手法

• 各物理スピンが結合次元Dの補助系のもつれペアを射影することで形成されるPEPSをMPSの一般化として提案する。
• 局所テンソル（射影子）を用いて量子状態をテンソルネットワークとして表現し、隣接するサイト間でもつれを生じさせる。補助自由度から物理スピンへ写像する。
• テンソルネットワークの縮約を用いてPEPSにおける相関関数を効率的に計算するアルゴリズムを開発し、バリエーション最適化を可能にする。
• Trotter分解を用いた時間発展演算アルゴリズムを実装し、時間ステップを部分ステップに分割。各部分ステップでは1つの相互作用のみが有効になるようにする。
• 各部分ステップ後に最適近似により結合次元を低減し、行単位で順次射影子を最適化することで、計算コストを低く維持する。
• 計算複雑性を扱うために多項式時間の縮約技術を用い、系のサイズNと結合次元Dに対してスケーラブルであることを保証する。

実験結果

リサーチクエスチョン

• RQ1 行列積状態（MPS）は2次元以上に一般化可能であり、強い相関を持つ系のスケーラブルなシミュレーションを可能にするか？
• RQ2 PEPSフレームワーク内での相関関数計算のための効率的アルゴリズムを開発可能か？ これはバリエーション手法を支援する。
• RQ3 PEPSに基づくアルゴリズムは、2次元スピン系（特にフラストレーションを持つモデルを含む）の基底状態および虚時間発展演算を正確にシミュレート可能か？
• RQ4 2次元ヘイゼンベルグ模型のような系において、結合次元Dの増加に伴いPEPS近似の精度はどのように変化するか？
• RQ5 PEPS形式は有限温度状態、実時間発展演算、散乱を伴う系へ拡張可能か？

主な発見

• PEPS形式は、DMRGに類似したバリエーション的手法を2次元以上へ自然に拡張し、強い相関を持つ系の正確なシミュレーションを可能にした。
• 4×4ヘイゼンベルグ反強磁性体において、D=3のとき変分エネルギーが正確な基底状態エネルギーと0.004%以内に収束しており、Dの増加に伴う収束の速さが明確に示された。
• D=2でさえも、正確な結果と優れた一致を示し、エネルギー差はD=1の35%からD=2の2%にまで減少した。
• アルゴリズムは基底状態に迅速に収束し、フラストレーションを持つヘイゼンベルグ模型のエネルギー差はDの増加に伴い精度が向上した。
• 10×10格子上でのシミュレーションでは、D=3がD=2に比べて基底状態エネルギーの推定において顕著な改善を示した。
• 反強磁性的相互作用およびフラストレーションを持つ相互作用の両方に対して、この手法は頑健であり、一般性と広範な応用可能性を示した。