[論文レビュー] Improving dynamical properties of metropolized discretizations of overdamped Langevin dynamics

Max Fathi, Gabriel Stoltz|arXiv (Cornell University)|May 19, 2015

Markov Chains and Monte Carlo Methods参考文献 29被引用数 17

ひとこと要約

この論文は、過減衰ランジュバン力学の数値解法において、強収束の順序を向上させたり、輸送係数推定のバイアスを低減させたりすることで、動的整合性を向上させる修正されたメトロポリス補正付き離散化スキームを提案する。標準のメトロポリス受容ルールをバーカー型ルールに置き換え、特化した提案分布を用いることで、不変測度を保持しつつ、拡散係数などの動的観測量の精度を顕著に向上させる。

ABSTRACT

The discretization of overdamped Langevin dynamics, through schemes such as the Euler-Maruyama method, can be corrected by some acceptance/rejection rule, based on a Metropolis-Hastings criterion for instance. In this case, the invariant measure sampled by the Markov chain is exactly the Boltzmann-Gibbs measure. However, rejections perturb the dynamical consistency of the resulting numerical method with the reference dynamics. We present in this work some modifications of the standard correction of discretizations of overdamped Langevin dynamics on compact spaces by a Metropolis-Hastings procedure, which allow us to either improve the strong order of the numerical method, or to decrease the bias in the estimation of transport coefficients characterizing the effective dynamical behavior of the dynamics. For the latter approach, we rely on modified numerical schemes together with a Barker rule for the acceptance/rejection criterion.

研究の動機と目的

• 過減衰ランジュバン力学の有限時間ステップ離散化が引き起こす輸送係数推定のバイアスを是正すること。
• 標準のメトロポリス＝ハスティングス補正が不変測度を保持するが、長時間ダイナミクスを歪める動的不整合を是正すること。
• 詳細つり合いを維持し、ギブス測度の正確なサンプリングを保ちつつ、動的観測量の収束性を向上させる数値スキームの開発。
• 受容ルール（特にバーカー則）が有効拡散係数などの輸送係数の精度に与える影響を調査すること。
• メトロポリス補正下でも強収束の順序を向上させるために、高次項を含む提案スキームを構築するフレームワークの提供。

提案手法

• 強収束の順序を向上させるために、ドリフトおよび拡散成分に高次項を組み込んだ修正された数値スキームを提案する。
• 輸送係数推定のバイアスを低減するため、標準のメトロポリス受容比をバーカー型受容ルールに置き換える。
• 時間ステップの高次項まで正確に、目的のスデの無限小生成作用素を再現するように設計された提案分布を構築する。
• 修正されたフォッカー＝プランクアプローチを用いて不変測度を導出し、新しい受容ルール下でも詳細つり合いが保たれることを保証する。
• 拡散行列の逆平方根およびその微分を含む新しい提案スキームクラスを導入し、動的整合性を向上させる。
• 摂動に対する自己拡散係数の線形応答を解析的に導出し、数値バイアスの正確なベンチマーク化を可能にする。

実験結果

リサーチクエスチョン

• RQ1 受容ルールの選択（メトロポリス対バーカー）が、拡散係数などの輸送係数推定のバイアスにどのように影響するか？
• RQ2 メトロポリス補正下でも、ドリフトおよび拡散成分に2次項を含む高次提案が、強収束の順序を向上させられるか？
• RQ3 拒否が、離散化された過減衰ランジュエルン過程の長時間ダイナミクスに与える影響は何か？
• RQ4 修正された提案が、数値的および真の有効拡散係数の乖離をどの程度低減できるか？
• RQ5 サンプリングの精度と連続ダイナミクスへの忠実度の両立を図る提案を体系的に構築する方法は何か？

主な発見

• 標準のメトロポリスルールの代わりにバーカー受容ルールを用いることで、自己拡散定数などの輸送係数推定におけるバイアスが顕著に低減する。
• ドリフトおよび拡散成分に2次項を含む修正提案は、メトロポリス補正下でも強収束の順序が向上する。
• 提案スキームは、正確な不変測度（ボルツマン＝ギブス分布）を保持しつつ、標準のMALA法よりも優れた動的整合性を達成する。
• 線形応答理論による自己拡散係数の解析的計算が、本研究で用いられた数値ベンチマークの正確性を確認する。
• 理論的解析により、新しいスキーム下で提案ノイズとドリフト補正の3次形式の期待値が消えることが示され、目的のスデと整合することが保証される。
• 数値実験により、新しい手法が、特に大きな時間ステップにおいて、標準のMALA法と比較して、有効輸送係数の参考値との一致が著しく良好であることが示された。