[論文レビュー] Continuum limit and preconditioned Langevin sampling of the path integral molecular dynamics

Jianfeng Lu, Yulong Lu+1 more|arXiv (Cornell University)|Nov 27, 2018

Quantum, superfluid, helium dynamics参考文献 44被引用数 5

ひとこと要約

本稿では、N→∞におけるリングポリマーの連続極限を活用して、経路積分分子動力学（PIMD）における効率的サンプリングを実現するための前処理付きランジュバンダイナミクスを提案する。二次的ポテンシャルにおいて次元に依存しないレートで指数的収束を確立し、非アディアバティックな多準位量子系とも整合するフレームワークを提供する。

ABSTRACT

We investigate the continuum limit that the number of beads goes to infinity in the ring polymer representation of thermal averages. Studying the continuum limit of the trajectory sampling equation sheds light on possible preconditioning techniques for sampling ring polymer configurations with large number of beads. We propose two preconditioned Langevin sampling dynamics, which are shown to have improved stability and sampling accuracy. We present a careful mode analysis of the preconditioned dynamics and show their connections to the normal mode, the staging coordinate and the Matsubara mode representation for ring polymers. In the case where the potential is quadratic, we show that the continuum limit of the preconditioned mass modified Langevin dynamics converges to its equilibrium exponentially fast, which suggests that the finite-dimensional counterpart has a dimension-independent convergence rate. In addition, the preconditioning techniques can be naturally applied to the multi-level quantum systems in the nonadiabatic regime, which are compatible with various numerical approaches.

研究の動機と目的

• ビーズ数Nの増加に伴い硬いばねポテンシャルが生じるため、経路積分分子動力学（PIMD）におけるサンプリングの非効率性という課題に対処すること。
• 高ビーズ数系のための前処理戦略を導くために、N→∞におけるリングポリマー動力学の連続極限を調査すること。
• PIMDシミュレーションにおける数値的安定性とサンプリング精度を向上させる前処理付きランジュバンダイナミクスを開発すること。
• 二次的ポテンシャルの場合における前処理付きダイナミクスの理論的収束保証を確立すること。
• 断熱的でない多準位量子系へとこのフレームワークを拡張し、断熱的でない基底におけるリングポリマー表現を用いること。

提案手法

• 配置空間におけるブラウン運動ループとしてリングポリマーの連続極限を形式化し、エネルギーをLℝ^d（ループ空間）上で定義する。
• E(𝒬) = ∫₀^β [½|∂τ𝒬(τ)|² + V(𝒬(τ))] dτ を用いて、ループ空間上に無限次元のギブス測度πを定義する。
• 質量の変更に基づくものと、前処理行列L^αを用いるものの、2種類の前処理付きランジュバンダイナミクスを提案する。
• モード分解を用いて、それらのダイナミクスを正規モード、ステージ座標、マツバラモードと関連付ける。
• 前処理付きダイナミクスのフォッカー・プランク方程式を導出し、不変測度が目的のギブス測度に対応することを検証する。
• 有効ハミルトニアンと観測量の重み関数を用いて、二準位量子系に対して断熱的でない基底における形式を拡張する。

実験結果

リサーチクエスチョン

• RQ1 ビーズ数N→∞におけるリングポリマー軌道の連続極限は、どのように振る舞うか？
• RQ2 高ビーズ数PIMDシミュレーションにおけるサンプリングの安定化と加速に適した前処理戦略は何か？
• RQ3 前処理付きランジュバンダイナミクスは、二次的ポテンシャルの場合に次元に依存しない収束レートを達成できるか？
• RQ4 提案されたダイナミクスは、正規モード、ステージ座標、マツバラモードといった既存の表現とどのように関係しているか？
• RQ5 前処理フレームワークは、非アディアバティックな多準位量子系へと拡張可能か？

主な発見

• 前処理付き質量変更ランジュバンダイナミクスは、二次的ポテンシャルの場合に、ビーズ数に依存しない収束レートで指数的に平衡状態に収束することが示された。
• 軌道のサンプリング方程式の連続極限が、有限NにおけるPIMDシミュレーションのための効果的な前処理戦略を導く構造的洞察を明らかにした。
• 提案された前処理付きダイナミクスは、さまざまな数値的手法と整合しており、非アディアバティック系へ自然に拡張可能であることが示された。
• モード解析により、前処理付きダイナミクスとマツバラモード表現との間の関連が明らかとなり、改善されたスペクトル性質が示唆された。
• 二準位量子系において、断熱的でない基底におけるリングポリマー表現により、提案されたダイナミクスを用いて一貫した熱平均のサンプリングが可能である。
• 前処理付きダイナミクスの不変測度が、目的のギブス測度と一致しており、サンプリングスキームの正しさが確認された。