[論文レビュー] Stable and accurate schemes for Langevin dynamics with general kinetic energies

Gabriel Stoltz, Zofia Trstanova|arXiv (Cornell University)|Sep 9, 2016

Markov Chains and Monte Carlo Methods参考文献 3被引用数 5

ひとこと要約

本稿では、非標準的な運動エネルギー（例：適応的拘束型、非グローバルリプシッツ型）を有するランジュバンダイナミクスに対して、安定で弱2次精度の数値スキームを提示する。分割積分法とメトロポリス補正を組み合わせることで、正確性と安定性を確保する。この手法により、ボルツマン＝ギブズ分布の効率的サンプリングが可能となり、複雑系における準安定状態の低減が達成される。

ABSTRACT

We study Langevin dynamics with a kinetic energy different from the standard, quadratic one in order to accelerate the sampling of the Boltzmann-Gibbs distribution. We consider two cases: kinetic energies which are local perturbations of the standard kinetic energy around the origin, where they vanish (this corresponds to the so-called adaptively restrained Langevin dynamics); and more general non-globally Lipschitz energies. We develop numerical schemes which are stable and of weak order two, by considering splitting strategies where the discretizations of the fluctuation/dissipation are corrected by a Metropolis procedure. We use the newly developped schemes for two applications: optimizing the shape of the kinetic energy for the adaptively restrained Langevin dynamics, and reducing the metastability of some toy models with non-globally Lipschitz kinetic energies.

研究の動機と目的

• 非標準的な運動エネルギーを用いるランジュバンダイナミクスにおいて、ボルツマン＝ギブズ分布を効率的にサンプリングする課題に対処すること。
• グローバルリプシッツでない、あるいは標準的な二次形式の局所的摂動である運動エネルギーに対しても、安定かつ正確な数値スキームを構築すること。
• 適応的拘束ランジュエルンダイナミクスにおける運動エネルギーの形状を最適化し、サンプリング性能を向上させること。
• 非グローバルリプシッツ型運動エネルギーを有するモデルにおいて、改善された数値積分法により準安定状態を低減すること。

提案手法

• 動的挙動を確定的および確率的成分に分割する分割積分法を設計し、揺動／散逸項に補正を適用する。
• 各積分ステップの後にメトロポリス・ハスティングス補正を適用することで詳細平衡を満たし、弱2次精度を確保する。
• 分割と補正ステップを巧みにバランスさせることで、特異性や非リプシッツ性を示す運動エネルギーに対しても安定なスキームを構築する。
• 提案されたスキームを用いて、数値実験を通じて適応的拘束ランジュエルンダイナミクスにおける運動エネルギー関数形の最適化を実施する。
• 非グローバルリプシッツ型運動エネルギーを有するトイモデルにスキームを適用し、準安定状態の低減と高速な収束を実証する。

実験結果

リサーチクエスチョン

• RQ1 運動エネルギーがグローバルリプシッツでない場合、ランジュエルンダイナミクスの数値スキームをどのようにして安定かつ弱2次精度に保てるか？
• RQ2 非標準的な運動エネルギー設定下で、メトロポリス補正が分割積分法の安定性および精度に及ぼす影響は何か？
• RQ3 適応的拘束ランジュエルンダイナミクスにおける運動エネルギーの形状を最適化することで、ボルツマン＝ギブズ分布のサンプリング速度を向上させられるか？
• RQ4 提案されたスキームは、非グローバルリプシッツ型運動エネルギーを有するモデルにおいて、どの程度準安定状態を低減できるか？

主な発見

• 提案されたメトロポリス補正付き分割スキームは、弱2次精度を達成し、非グローバルリプシッツ型運動エネルギーに対しても安定である。
• 特異な運動エネルギー関数を有するトイモデルにおいて、準安定状態が著しく低減され、平衡状態への収束が高速化された。
• 適応的拘束ランジュエルンダイナミクスにおける運動エネルギーの形状最適化により、数値的に示されるようにサンプリング効率が向上した。
• この手法は、ターゲットのボルツマン＝ギブズ分布の不変測度を保持しており、長時間的な統計的挙動が正しく保たれている。