[論文レビュー] DMRG Approach to Optimizing Two-Dimensional Tensor Networks

Katharine Hyatt, E. Miles Stoudenmire|arXiv (Cornell University)|Aug 23, 2019

Quantum many-body systems被引用数 18

ひとこと要約

この論文は、元のDMRG法の主要な特徴—正規形、効率的な固有値ソルバー、安定した大規模更新—を組み込むことで、2次元PEPSテンソルネットワークを最適化するDMRGに類似したアルゴリズムを提案する。この手法は2次元スピン模型において高速かつ高精度の収束を達成し、結合次元D=6で正確なQMC結果と同等の精度を示し、エネルギーの指数的収束を示している。

ABSTRACT

Tensor network algorithms have been remarkably successful solving a variety of problems in quantum many-body physics. However, algorithms to optimize two-dimensional tensor networks known as PEPS lack many of the aspects that make the seminal density matrix renormalization group (DMRG) algorithm so powerful for optimizing one-dimensional tensor networks known as matrix product states. We implement a framework for optimizing two-dimensional PEPS tensor networks which includes all of steps that make DMRG so successful for optimizing one-dimension tensor networks. We present results for several 2D spin models and discuss possible extensions and applications.

研究の動機と目的

• 1次元MPSに成功した強力なDMRGアルゴリズムを、2次元PEPSテンソルネットワークに拡張すること。
• 既存のPEPS最適化手法に欠けている正規形、効率的な固有値ソルバー、安定した大規模更新の問題を解決すること。
• テンソルネットワークを用いた2次元量子多体系の高精度でスケーラブルなシミュレーションを可能にすること。
• 将来の改良の基盤とすること、例えば自己適合的結合次元の増加や2サイト更新の導入。

提案手法

• 1サイトDMRGアプローチをPEPSに採用し、他のテンソルを固定したまま1つのテンソルだけを更新する。
• 最近のPEPS正規化の進展を活用して、正規直交性と数値的安定性を保証する正規形をネットワークに導入する。
• 固定されたテンソルが定義する正規直交基底にハミルトニアンを射影し、通常の固有値問題を形成する。
• 効率性と正確性を兼ね備えた反復的ソルバー（例：LanczosやDavidson）を用いて、射影された固有値問題を解く。
• 全体のネットワークの安定性を保ちながら、個々のテンソルに対して大規模で安定した更新を実行する。
• ITensorライブラリのGPUアクセceleratedバックエンドを活用し、特にテンソル結合とMPO-MPO乗算の処理でパフォーマンスを最適化する。

実験結果

リサーチクエスチョン

• RQ1 DMRGアルゴリズムの核となる利点—正規形、通常の固有値問題、安定した大規模更新—は、2次元PEPSに成功裏に適応可能か？
• RQ2 2次元スピン模型において、PEPS-DMRGの収束速度は正確な結果やQMC結果と比べてどうか？
• RQ3 有限な結合次元と正規化の精度の誤差が、PEPS-DMRGにおけるエネルギーとスピンスピン相関関数に与える影響は何か？
• RQ4 2サイトアルゴリズムを用いて、自己適合的結合次元の増加を拡張可能か？
• RQ5 有限サイズのPEPS-DMRGは、バルク性質の推定において無限大PEPS（iPEPS）と比べてどのように性能を発揮するか？

主な発見

• PEPS-DMRGアルゴリズムは、DMRGスイープ数に対してエネルギーが指数関数的に収束し、L=10系では約10スイープで最終エネルギーの10⁻²以内に到達する。
• L=10のヘイゼンベルグ模型において、結合次元D=6、環境次元χ=12を用いることで、エネルギー/サイトは正確なQMC結果と10⁻²以内に収束する。
• D=6のPEPSで計算されたスピンスピン相関関数は、正確なMPS-DMRG結果と比較して有限結合次元の差異を示すが、最新のiPEPS手法と同等の水準である。
• 最適化中に偶発的にエネルギーが上昇する現象は、正規化が不完全であることが原因であると特定され、正規化アルゴリズムの改善が安定性と正確性の向上に寄与する。
• 合計実行時間の46%がテンソル結合、14%がMPO-MPO乗算を占めており、アルゴリズム的最適化の余地が顕著に示されている。
• GPUアクセcelerated ITensorバックエンドは、シミュレーション速度を顕著に向上させ、より大きな系や高い結合次元の探索を効率的に行えるようにした。