[論文レビュー] Hamiltonian Monte Carlo on Symmetric and Homogeneous Spaces via Symplectic Reduction

Alessandro Barp, A.D. Kennedy+1 more|arXiv (Cornell University)|Mar 7, 2019

Advanced Algebra and Geometry参考文献 47被引用数 5

ひとこと要約

この論文は、対称的かつ自然に再帰的な同次空間におけるハミルトニアン・モンテカルロ（HMC）の一般枠組みを提示する。シンプレクティック還元を用いて、グローバル座標と定数計量を持つ行列リー群に動的を上げることで、座標特異性や制約を回避する。この手法により、球面、スタイリル・マニフォールド、グラスマン多様体、対称正定値行列などの曲がった多様体上でも、幾何学的に一貫した効率的なサンプリングが可能になる。

ABSTRACT

The Hamiltonian Monte Carlo method generates samples by introducing a mechanical system that explores the target density. For distributions on manifolds it is not always simple to perform the mechanics as a result of the lack of global coordinates, the constraints of the manifold, and the requirement to compute the geodesic flow. In this paper we explain how to construct the Hamiltonian system on naturally reductive homogeneous spaces using symplectic reduction, which lifts the HMC scheme to a matrix Lie group with global coordinates and constant metric. This provides a general framework that is applicable to many manifolds that arise in applications, such as hyperspheres, hyperbolic spaces, symmetric positive-definite matrices, Grassmannian, and Stiefel manifolds.

研究の動機と目的

• 局所座標による特異性や数値的不安定性を引き起こすため、曲がった多様体上でのハミルトニアン・モンテカルロの実行に直面する課題に対処すること。
• 多様体上でのHMCにおいて、制約力やラグランジュ乗数の必要性を排除するために、グローバル座標を持つリー群に動的を埋め込むこと。
• 球面、双曲空間、スタイリル・マニフォールドおよびグラスマン多様体、対称正定値行列を含む広範なクラスの多様体に適用可能な、一般化された幾何学的に一貫したHMC枠組みを提供すること。
• 行列群の演算を用いて、持ち上げた位相空間上で対称なシンプレクティック積分子を構築することで、効率的なサンプリングを実現すること。
• 群不変性と運動量写像の保存を通じて、多様体制約を自然に維持する動的を確保し、数値的ドリフトを回避すること。

提案手法

• シンプレクティック還元を用いて、$\nG$ という行列リー群上に、目的の多様体 $\nM = \nG/\nK$ に射影されるハミルトニアン系を構築し、幾何的構造を保存する。
• $\nK$-作用に関連する運動量写像 $J$ を用いて、制約 $J = 0$ を強制し、明示的な制約なしに、軌道が同次空間上に留まるように保証する。
• HMCの動的をグローバル座標と行列乗算により統合が容易なリー群 $\nG$ に持ち上げることで、局所パrameter化に起因する特異性を回避する。
• シンプレクティック構造と群不変性を保存する基本ステップの合成として、対称なシンプレクティック積分子を構築し、長期的な安定性を確保する。
• $\nG$ 上のモーラー＝カルタン形式を用いて、リー代数値1形式として動的を表現し、行列指数関数による測地線フローの効率的計算を可能にする。
• ネーターの定理を応用して対称性を活用する：ハミルトニアンの $\nK$-不変性から、保存量（一般化運動量）が得られ、位相空間の次元削減とサンプリングの簡略化に用いられる。

実験結果

リサーチクエスチョン

• RQ1 球面やグラスマン多様体のような曲がった非ユークリッド多様体上でのHMCを、局所座標に依存せずに一般化する方法は何か？
• RQ2 シンプレクティック還元を用いて、対称的かつ自然に再帰的な同次空間上に、グローバルに有効なHMC枠組みを構築できるか？
• RQ3 多様体制約とシンプレクティック構造を保存しながら、グローバル座標を持つリー群に動的を持ち上げる方法は何か？
• RQ4 運動量写像は、ポテンシャル障壁やラグランジュ乗数を用いずに、多様体制約をどのように強制するか？
• RQ5 行列群の演算を用いて、曲がった多様体上で対称的かつシンプレクティックな積分子をどのように設計するか？これにより、数値的安定性と正確性が保証されるか？

主な発見

• 提案手法により、$S^n$、$\text{St}(k,n)$、$\text{Gr}(k,n)$、$\text{Sym}^+(n)$ などの対称的かつ自然に再帰的な同次空間上でのHMCサンプリングが、座標特異性や制約力なしに可能になる。
• 動的を行列リー群に持ち上げることで、グローバル座標表現と定数計量を達成し、数値積分の簡素化と座標境界での問題の回避が実現される。
• シンプレクティック還元の使用により、ハミルトニアン系が多様体の幾何学的・位相的構造を尊重し、運動量写像による自然な制約により、軌道が同次空間上に留まる。
• シンプレクティック構造と群不変性を保存する、効率的な対称シンプレクティック積分子がサポートされ、安定した長期的サンプリング行動が得られる。
• 2次元球面上での実験では、軌道長6,371 km、ステップサイズ63.71 kmの条件下で、制約違反なしに安定したサンプリングが達成された。
• このアプローチは、すべての自然に再帰的な同次空間に一般化可能であり、非ユークリッドなサンプリングを伴う広範な統計的・物理的問題に適用可能な統一枠組みを提供する。