Post-processed frozen-flow methods for the long time sampling of ergodic dynamics on Riemannian manifolds

技術要約：リーマン多様体上のエルゴード力学の長時間サンプリングのためのポストプロセス済みフローズンフロー法

問題提起
本論文は、リーマン多様体上のエルゴード的な確率微分方程式（SDE）の不変測度をサンプリングするための高次数値手法を設計するという課題に取り組んでいる。具体的には、著者らは以下の形式のSDEを検討している：
$$dX(t) = F(X(t))dt + \sqrt{2} \sum_{d=1}^D E_d(X(t)) \circ dW_d(t)$$
ここで、$M$は滑らかで完備な連結リーマン多様体であり、$E_d$は大域的な正規直交基底を形成し、$F$はベクトル場である（リーマン・ランジュバン動力学の場合、$F$はポテンシャル$V$から導出される可能性がある）。

既存のアプローチには、以下のような重大な制限がある：

1. 外在的手法（例：射影法）：多様体をより高次元のユークリッド空間に埋め込む。これらは便利ではあるが、しばしば厳しいタイムステップの制限や安定性の問題に直面する。
2. 既存の内在的手法：通常、強収束または弱1次収束に焦点を当てており、特定の多様体（リー群など）に限定されている。また、エルゴード測度に必要な長時間のサンプリング効率を優先していない。

目標は、幾何学的操作（測地線や平行移動）のみに依存し（埋め込みや局所座標を使用しない）、不変測度のための高次精度を実現する内在的な手法を開発することであり、有限時間の精度よりも長時間のサンプリング効率を優先するものである。

手法
著者らは、ポストプロセス済みフローズンフロー法に基づく新しいフレームワークを提案している。このアプローチは、2つの主要な概念を組み合わせている：

1. フローズンフロー法：Crouch-Grossmanおよび交換子フリーのリー群手法の確率的拡張である。これらの手法は、レヴィ＝チヴィタ接続ではなく、特定の曲率のない接続であるワイツェンベック接続に関連する基本的な幾何学的操作（測地線と平行移動）を用いて多様体上で進化する。
2. ポストプロセッシング（後処理）：軌道の最後に一度だけ適用される、恒等写像のランダムな摂動である。このテクニックは、ユークリッド空間におけるLeimkuhler-Matthews法に着想を得ており、ステップあたりの計算コストを大幅に増やすことなく、不変測度に対する高次の精度を実現する。

これらの手法の次数条件を導出するために、著者らはエキゾチックLie-Butcher級数に基づく代数的形式を利用している。彼らは以下を導入した：

• 平面エキゾチックフォレスト：ラプラス・ベルトラミ作用素を表す「エキゾチック」な装飾を取り入れた、標準的なButcher木のリーマン多様体への一般化。
• フォレスト上の部分積分（IBP）：ギブス測度に対するリーマン部分積分公式から導かれた新しい代数演算。これにより、数値作用素の随伴（$A_j^*$）を既約フォレストを用いて表現することが可能となる。
• 修正シャッフル代数：高次の不変測度の次数条件を特徴付けるために、修正シャッフル代数のキャラクタを定義する。

主な貢献

1. 内在的な高次特性化：本論文は、リーマン多様体上の不変測度の高次精度のための厳密な基準を提供している。これは、係数の条件が満たされれば、不変測度の誤差が$O(h^p)$で減衰することを証明することで、従来のユークリッド的な結果（例：[2, 24, 65]）をリーマン多様体の設定へと一般化したものである。
2. 新しい数値スキーム：著者らは3つの新しい手法を導出し、提示している：
   • 手法1：Leimkuhler-Matthews法をリーマン多様体に一般化したポストプロセス済みフローズンフロー法。1ステップあたり1回のベクトル場$F$の評価と3つの指数写像を必要とする。これは弱1次であるが、不変測度に対しては2次である。
   • 手法2および3：ポストプロセッシングを伴わない、代替スキーム（確率的フローズンフローHeunおよび2ステップ・フローズンフローRunge-Kutta）。これらも不変測度に対して2次である。
3. 代数的フレームワーク：本論文は、平面エキゾチックフォレスト上の部分積分演算を導入し、一般的なフォレストを既約なものへと写す簡約演算子（RED）を定義している。これにより、不変測度のための次数条件を体系的に導出することが可能となり、不変測度のための条件数が、弱誤差のための条件数（例：弱誤差で8条件に対し、不変測度では5条件）よりも大幅に少ないことを示している。
4. コスト効率：提案手法は、複数のベクトル場評価や複雑な座標変換を回避するため、既存の外在的射影法や他の内在的な高次スキームよりも計算コストが低いように設計されている。

結果
理論的な主張は、2つの古典的な多様体を用いた数値実験を通じて検証されている：

1. $SO(p)$（特殊直交群）：著者らは、二次および六次のポテンシャルを持つコンパクトなリー群を用いて手法をテストした。結果は、新しい手法（1, 2, 3）が不変測度に対して期待通りの2次収束を達成することを裏付けた。特に、手法1はSFF2法（参照用となる弱2次の内在的手法）と比較して計算コストが低く、二次形式の場合にモンテカルロ閾値に迅速に到達した。
2. $S^2$（球面）：Von-Mises Fisher動力学を用いて実験が行われた。誤差曲線は、やはり不変測度に対する期待通りの2次収束を示した。手法1はここでも優れた性能を示し、ほぼ瞬時にモンテカルロ閾値に到達した。

両方のケースにおいて、提案手法は、与えられた精度を達成するために必要なコストにおいて、既存の外在的および内在的なアプローチを上回った。

意義と主張
本論文は、リーマン設定へのLeimkuhler-Matthews法の直接的な一般化を提供し、最小限の計算オーバーヘッド（1回のベクトル場評価）で不変測度の2次精度を達成することを主張している。

著者らは、自身のアプローチが以下の点に重点を置いていることを強調している：

• 完全に内在的であり、埋め込みや局所座標を回避している。
• 有限時間の軌道精度ではなく、エルゴード的動力学に不可欠な長時間のサンプリング効率を優先している。
• 不変測度のための任意の高次のインテグレータを設計するための体系的な代数的メソドロジー（エキゾチックフォレストおよびpost-Lie代数による）を提供している。

論文は、現在の手法が（曲率を持たないがレヴィ＝チヴィタとは異なる）ワイツェンベック接続に依存している一方で、将来の研究では、これらのフローズンフロー法を（レヴィ＝チヴィタ接続のみに依存するように）拡張し、マルチプリカティブノイズを伴うSDEへと一般化することを目指していると述べて締めくくっている。また、開発された代数的形式は、粗いパス（rough paths）や繰り込み理論などの分野において、数値解析を超えた応用可能性を持っていることも強調されている。