[論文レビュー] Least Action Principle for the Real-Time Density Matrix Renormalization Group

Kouji Ueda, Chenglong Jin+3 more|arXiv (Cornell University)|Dec 19, 2006

Quantum many-body systems被引用数 3

ひとこと要約

本稿では、時間発展演算子の累積的切断誤差に起因する時間反転対称性の破れを回復し、誤差の蓄積を低減するために、リアルタイム密度行列密度行列群（DMRG）における最小作用原理を導入する。時間発展演算子の誤差の二乗に基づく離散的作用関数を最小化することで、反復的かつ並列的な更新を用いて、すべての時間スライスにおいて同時に行列積状態（MPS）表現を最適化し、対称的かつ高精度な時間発展演算を実現する。

ABSTRACT

A kind of least action principle is introduced for the discrete time evolution of one-dimensional quantum lattice models. Based on this principle, we obtain an optimal condition for the matrix product states on succeeding time slices generated by the real-time density matrix renormalization group method. This optimization can also be applied to classical simulations of quantum circuits. We discuss the time reversal symmetry in the fully optimized MPS.

研究の動機と目的

• 時間発展演算における累積的切断誤差によって生じるリアルタイムDMRGにおける時間非対称性を解消する。
• 時間区間全体にわたる変分原理を定式化することで、数値的時間発展演算における時間反転対称性を回復する。
• 時間発展演算MPSにおける繰り返しの正規化ステップに起因する数値誤差の蓄積を最小化する。
• 時間発展演算MPSの全時間スライスを同時に最適化するグローバル変分関数を用いることで、より正確で対称的な時間発展演算を可能にする。
• 適応的時間ステップ、シンプレクティック積分法、非局所的または非可逆的演算子を含めるために、変分的再定式化を用いた拡張可能なフレームワークを提供する。

提案手法

• シュレーディンガー方程式からの逸脱を定量化するため、離散的作用関数 $ I = \int \left| \left( \frac{\partial}{\partial t} - \frac{H}{i\hbar} \right) |\Psi(t)\rangle \right|^2 dt $ を定義する。
• 時間を $ N $ 個の区間に離散化し、ステップサイズを $ \Delta t $ とする。短時間発展演算を $ \mathcal{T} = \exp(\Delta t H / i\hbar) $ と近似する。
• 各時間スライスにおける時間発展演算状態を行列積状態（MPS）で表現し、局所テンソルを反復的に最適化する。
• 全時間スライスにわたる作用関数の最小化を変分最適化により実行し、時間発展演算全体を一つの変分問題として扱う。
• 前向きおよび後向きの時間方向に並列スイープを実施し、各時間点におけるMPSテンソルの対称的最適化を可能にする。
• 伝搬行列技術および正規化スキーム（例：$ B[s] \Lambda^B B[s] $）を用いて、結合次元 $ m $ を動的に調整し、精度を向上させる。

実験結果

リサーチクエスチョン

• RQ1 リアルタイムDMRGシミュレーションにおいて、切断誤差が時間反転対称性を破る場合、どのようにして時間反転対称性を回復できるか？
• RQ2 繰り返しの正規化ステップに起因する累積誤差を最小化するための変分原理として、どのようなものが適しているか？
• RQ3 全時間スライスにわたるグローバル最適化は、リアルタイムDMRG時間発展演算の正確性と対称性を向上させられるか？
• RQ4 作用関数 $ I = \int \mathcal{L} dt $、ここで $ \mathcal{L} = \left| \left( \frac{\partial}{\partial t} - \frac{H}{i\hbar} \right) |\Psi(t)\rangle \right|^2 $ が時間発展演算の安定化に果たす役割は何か？
• RQ5 本手法は、非局所的または非可逆的演算子を扱うために、あるいは適応的時間ステップやシンプレクティック積分法を組み込むために、どのように拡張可能か？

主な発見

• 最小作用原理により、全時間スライスにわたるグローバル作用関数の最小化が実現され、リアルタイムDMRGにおける時間反転対称性が回復され、前向きと後向きの時間発展演算における非対称性が低減される。
• グローバル変分最適化により、切断誤差の蓄積を回避することで、数値誤差の指数的増大を防ぐ。
• 作用関数に基づく変分法により、前向きおよび後向きの時間スイープを並列で計算可能となり、計算効率と対称性が向上する。
• 拡張された正規化波動関数 $ B[s_\ell] \Lambda^B B[s_{\ell+1}] $ を用いることで、結合次元 $ m $ の動的調整が可能となり、発展演算中の精度が向上する。
• マルチターゲティングスキームを用いることで、長距離相互作用を有する系に対しても適用可能であり、可逆的伝搬行列を有する古典的統計的モデルや量子回路に対しても適用可能である。
• 本手法は、シンプレクティック積分法、適応的時間ステップ、その他の高度な数値技法を統合可能な変分フレームワークを提供する。