[論文レビュー] A Monte Carlo Time-Dependent Variational Principle

F. W. G. Transchel, Ashley Milsted+1 more|arXiv (Cornell University)|Nov 20, 2014

Advanced Thermodynamics and Statistical Mechanics被引用数 4

ひとこと要約

本稿では、 Lindblad マスター方程式に従う散乱量子ダイナミクスをシミュレートするための、時刻依存変分原理（TDVP）のモンテカルロ拡張を提案する。Lindblad生成子から導かれる Fokker-Planck 方程式のサンプリングを通じて、純粋状態の変分クラス（例：行列積状態（MPS））を用いて、混合状態のダイナミクスを、並列性が著しく高く、より大きな系へスケーラブルに効率的にシミュレート可能にする。

ABSTRACT

We generalize the Time-Dependent Variational Principle (TDVP) to dissipative systems using Monte Carlo methods, allowing the application of existing variational classes for pure states, such as Matrix Product States (MPS), to the simulation of Lindblad master equation dynamics. The key step is to use sampling to approximately solve the Fokker-Planck equation derived from the Lindblad generators. An important computational advantage of this method, compared to other variational approaches to mixed state dynamics, is that it is "embarrassingly parallel".

研究の動機と目的

• Lindblad マスター方程式で記述される散乱系に、時刻依存変分原理（TDVP）を拡張すること。
• 混合状態に対して一意な情報距離測度が存在しないという問題を、純粋状態波動関数のモンテカルロサンプリングを用いることで克服すること。
• MPS などの既存の変分クラスを用いて、オープンな量子多体系の効率的かつ明示的に並列化可能なシミュレーションを可能にすること。
• 制御された散乱を伴う可解性・非可解性を持つスピンチェーンに対して、この手法をテストし、解析解と比較すること。
• この手法が、散乱系における量子もつれと相関ダイナミクスを捉えるために、結合次元がどの程度必要となるかを評価すること。

提案手法

• Lindblad マスター方程式から Fokker-Planck 方程式を導出し、波動関数パラメータの確率的進化を記述する。
• Fokker-Planck 方程式の解を近似するためにモンテカルロサンプリングを用い、変分パラメータの確率的進化を可能にする。
• サンプルされた波動関数に TDVP 形式を適用し、変分パラメータの時間発展を記述する流れ方程式を導出する。
• 密度行列を純粋状態射影子の重み付き積分として表現する：$\rho_t = \int p_t(\bar{a},a) |\Psi(a)\rangle\langle\Psi(a)| \, dad\bar{a}$、ここで $|\Psi(a)\rangle$ は MPS などの変分クラスである。
• 数値的統合のために、時刻に伴うブロック縮約（TEBD）に類似したスキームを用いて、行列積状態（MPS）で手法を実装する。
• 各プロセッサに独立してサンプリング軌道を分散させることで、計算の並列性を確保し、明示的に並列化可能な実行を実現する。

実験結果

リサーチクエスチョン

• RQ1 純粋状態の変分アンザッツのモンテカルロサンプリングを用いることで、TDVP 形式を散乱ダイナミクスに一般化できるか？
• RQ2 この手法の精度は、モンテカルロサンプル数および MPS の結合次元にどのように依存するか？
• RQ3 エッジ散乱を伴う可解系（例：XXZ ハイゼンベルグチェーン）の非平衡ダイナミクスを、この手法が正確に捉えることができるか？
• RQ4 空間的に非一様な散乱を伴う非可解系において、干渉やドメイン形成を捉える能力はいかがなものか？
• RQ5 散乱系における長距離相関およびもつれを捉えるために、MPS の結合次元 $D$ はどの程度増大する必要があるか？

主な発見

• 高いモンテカルロサンプル数（例：48,000）を用いる場合、エッジ散乱を伴う XXZ ハイゼンベルグチェーンに対して、解析解と良好な一致を示し、$t=80$ から収束が観察された。
• より大きな系（$N_{\text{sites}}=16$）においては、結合次元 $D$ を増加させることで相関クラスタリングの捉えが向上し、$D=32$ は $D=2$ よりも顕著に良好な一致を示した。
• 両端に異なるポンプ（上・下）を施した二重均一散乱系では、中心付近に干渉パターンを捉え、特に弱い相互作用および大きな $\lambda$ の場合に顕著な非自明なもつれと相関ダイナミクスを示した。
• 中心への情報伝播が遅いにもかかわらず、臨界的・収束が遅い状態（例：エッジ駆動の XXZ チェーン）においても、解析結果と定性的に一致した。
• 干渉やもつれを捉えるために、結合次元 $D$ が $2^{N/2}$ に達する必要はなく、小さな $D$ でも関連するもつれ構造を効率的に捉えることができた。
• この手法は数値的に安定でスケーラブルであり、サンプル数と結合次元を増加させることで精度が体系的に向上する傾向が示された。これにより、正確な対角化を超える系のシミュレーションに応用可能であることが裏付けられた。