[論文レビュー] Geometric Ergodicity of Two--dimensional Hamiltonian systems with a Lennard--Jones--like Repulsive Potential

Ben Cooke, David P. Herzog+3 more|arXiv (Cornell University)|Apr 19, 2011

Markov Chains and Monte Carlo Methods参考文献 19被引用数 4

ひとこと要約

本稿では、ストキャスティック平均化技術を用いて新たなリャプノフ関数を構築することで、レナード・ジョーンズ型の反発的ポテンシャルを有する2次元ランジュバン力学における幾何的エルゴドゥーシティを確立した。主な貢献は、特異性が存在する中でも指数的収束が幾何的レートで達成されることの厳密な証明であり、このような系における分子動力学シミュレーションの厳密な解析における長年の空白を解消した。

ABSTRACT

In this paper we establish the ergodicity of Langevin dynamics for simple two-particle system involving a Lennard-Jones type potential. To the best of our knowledge, this is the first such result for a system operating under this type of potential. Moreover we show that the dynamics are {\it geometrically} ergodic (have a spectral gap) and converge at a geometric rate. Methods from stochastic averaging are used to establish the existence of a Lyapunov function. The existence of a Lyapunov function in this setting seems resistant to more traditional approaches. This is a corrected version of the article.

研究の動機と目的

• 特異なレナード・ジョーンズ型ポテンシャルを有する2粒子系におけるランジュバン力学の幾何的エルゴドゥーシティを厳密に確立すること。
• 近接粒子相互作用に起因する特異性が存在する状況において収束速度を証明する課題に対処すること。
• 非特異なポテンシャルに限った既存のエルゴドゥーシティ結果を、特異的かつ非一様なポテンシャルに適したリャプノフ関数の構築によって拡張すること。
• エルゴドゥーシティがしばしば仮定されるが、証明されていない分子動力学シミュレーションにおける収束挙動の理解の基盤を提供すること。

提案手法

• 系の高エネルギー極限にストキャスティック平均化技術を適用し、リャプノフ関数を構築する。
• 補助的ポisson方程式手法を用いて、時間平均エネルギー散逸率を取り扱い、無限小の負のドリフトを持つ関数を導出する。
• 高エネルギー力学を決定的ハミルトニアン系の小さな摂動とみなす摂動的アプローチを採用する。
• 特異性付近および無限遠における挙動を制御するため、エネルギー領域（$\mathcal{S}_i^\pm$）ごとにリャプノフ関数を区分的構成する。
• [AKM12, HM15a, HM15b]からの漸近的解析および近似技術を用いて、無限遠におけるダイナミクスをモデル化する。
• リャプノフ関数の妥当性を検証し、$\mathcal{L}\Psi \leq -\gamma(\Lambda_* - \epsilon)H + \text{lower-order terms}$を示し、指数的減衰を保証する。

実験結果

リサーチクエスチョン

• RQ1 特異性が存在するにもかかわらず、レナード・ジョーンズ型ポテンシャルを有する2次元ランジュバン力学における幾何的エルゴドゥーシティを厳密に確立できるか？
• RQ2 非一様かつ特異的ポテンシャルを有する系において、指数的収束を保証するリャプノフ関数を構築することは可能か？
• RQ3 ストキャスティック平均化技術は、レナード・ジョーンズポテンシャルのようなスケール不変性のないポテンシャルに対応してどのように改造できるか？
• RQ4 特異的相互作用を有する系におけるリャプノフ関数の構築において、エネルギー散逸率の役割は何か？
• RQ5 現在の制限を考慮しても、この手法は多粒子系や二次的成長を示すポテンシャルへと拡張可能か？

主な発見

• 本稿では、レナード・ジョーンズ型ポテンシャルを有する2粒子ランジュバン力学における幾何的エルゴドゥーシティを証明し、全変動距離における指数的収束が確立された。
• 無限小生成作用素が $\mathcal{L}\Psi \leq -\gamma(\Lambda_* - \epsilon)H + \text{lower-order terms}$ を満たす、新規のリャプノフ関数が構築され、幾何的収束が保証された。
• 本手法はポテンシャルの非一様スケーリングに対処でき、従来の手法の主な障壁を克服した。
• ストキャスティック平均化と無限遠における漸近的解析に依拠することで、高エネルギー状態への逸脱を制御できた。
• この結果は、ある $\delta > 0$ に対して $\mathbf{E}_\mu(|f|^{2+\delta}) < \infty$ を満たす関数のエルゴドリック平均に対して中心極限定理の存在を示唆し、正確なサンプリングスキームを支持する。
• 本手法は、周期軌道が即時の平均化を妨げるため、二次的ポテンシャルへは直接拡張できないが、今後の道筋として周期平均エネルギーのアプローチが提案されている。