[論文レビュー] Boltzmann-Gibbs Preserving Stochastic Variational Integrator

Nawaf Bou‐Rabee, Houman Owhadi|arXiv (Cornell University)|Dec 26, 2007

Model Reduction and Neural Networks被引用数 1

ひとこと要約

本稿では、変分積分法とオーランシュタイン＝ウーレンボルグ過程を組み合わせたLie-Trotter分割法を、慣性のあるランジュヴィン方程式に適用し、ボルツマン＝ギブズ不変測度を保存する手法を提案する。理論的に、離散的不変測度が変分積分法のエネルギー精度の範囲内で真の不変測度を近似することを証明しており、エネルギー誤差がゼロのときには正確なサンプリングが達成される。

ABSTRACT

This paper presents a Lie-Trotter splitting for inertial Langevin equations (Geometric Langevin Algorithm) and analyzes its long-time statistical properties. The splitting is defined as a composition of a variational integrator with an Ornstein-Uhlenbeck flow. Assuming the exact solution and the splitting are geometrically ergodic, the paper proves the discrete invariant measure of the splitting approximates the invariant measure of inertial Langevin to within the accuracy of the variational integrator in representing the Hamiltonian. In particular, if the variational integrator admits no energy error, then the method samples the invariant measure of inertial Langevin without error. Numerical validation is provided using explicit variational integrators with first, second, and fourth order accuracy.

研究の動機と目的

• 長期間にわたり統計的精度を保つ構造保存型数値積分法を、慣性のあるランジュヴィン方程式に対して開発すること。
• 変分積分法とオーランシュタイン＝ウーレンボルグ過程の組み合わせからなる分割積分法の不変測度を分析すること。
• 離散的不変測度が、確率的ダイナミクスの真の不変測度と一致する条件を確立すること。
• 1次、2次、4次精度の明示的変分積分法を用いて、手法の妥当性を検証すること。

提案手法

• 本手法は、Lie-Trotter分割を用い、慣性のあるランジュヴィンダイナミクスを2つの部分に分解する。1つは変分積分法によって制御されるハミルトニアン流れであり、もう1つはオーランシュタイン＝ウーレンボルグ過程としてモデル化された確率的減衰項である。
• 変分積分法は、元のハミルトニアン系のシンプレクティック構造とエネルギー保存性を保持するように構築される。
• オーランシュタイン＝ウーレンボルグ流れは、慣性のあるランジュヴィン方程式における摩擦項とノイズ項をモデル化し、正しい統計的挙動を保証する。
• これらの2つの流れの合成により、離散的確率過程が得られ、その不変測度が真のボルツマン＝ギブズ測度への収束について分析される。
• 解析では、正確な解と分割スキームの両方について幾何的エルゴード性を仮定し、不変測度の収束を確立する。
• 数値的妥当性は、1次、2次、4次精度の明示的変分積分法を用いて確認され、不変測度のサンプリングにおける精度が示された。

実験結果

リサーチクエスチョン

• RQ1 変分積分法とオーランシュタイン＝ウーレンボルグ過程の組み合わせによる分割積分法は、慣性のあるランジュヴィンダイナミクスの不変測度を保存できるか？
• RQ2 分割法の離散的不変測度は、真の不変測度をどの程度正確に近似するか？
• RQ3 どのような条件下で、この手法が不変測度を正確にサンプリングするか？
• RQ4 変分積分法の精度の次数は、サンプリングの統計的精度にどのように影響するか？

主な発見

• 分割積分法の離散的不変測度は、変分積分法のエネルギー誤差の範囲内で、慣性のあるランジュヴィン方程式の真の不変測度を近似する。
• 変分積分法にエネルギー誤差が存在しない場合、この手法は不変測度を正確にサンプリングし、統計的誤差がゼロとなる。
• 正確な解と同様の条件下で、幾何的エルゴード性が維持され、長期的な安定性と収束性が保証される。
• 数値実験により、1次、2次、4次精度の変分積分法が、それぞれ対応する精度で不変測度のサンプリングを達成することが確認された。