Jacobi Hamiltonian Integrators: construction and applications

技術要約：ヤコビ・ハミルトニアン積分器：構成と応用

問題提起
シンプレクティック多様体上で定義される古典的なハミルトニアン系は、幾何学的に十分に理解されており、その数値積分におけるシンプレクティック積分器は成熟した分野である。しかし、散逸、時間依存性、または特異性を伴う多くの物理系は、より一般的な幾何学的構造、具体的にはヤコビ多様体（Jacobi manifolds）によってモデル化される方が適切である。ヤコビ幾何学はポアソン幾何学とコンタクト幾何学を統合するものであり、保存的および散逸的力学の両方の枠組みを提供する。ポアソン系やコンタクト系のための幾何学的積分器は存在するものの、一般的なヤコビ・ハミルトニアン系の特定の幾何学的構造を保存する積分器を構築するための体系的な枠組みは欠けていた。標準的な積分器は、ヤコビ流の定性的な幾何学的特性を保持できないことが多く、その結果、長期的な挙動の悪化やエネルギー散逸プロファイルの不正確さを招く。

手法
著者らは、ヤコビ多様体と斉次ポアソン多様体の幾何学的関係を利用することで、**ヤコビ・ハミルトニアン積分器（JHI）**を構築するための体系的な枠組みを提案している。その手法は以下のステップで進行する：

1. ポアソニゼーション（Poissonization）：ハミルトニアン $H$ を持つヤコビ多様体 $(J, \Lambda, E)$ を、斉次ポアソン多様体 $(P, \Pi, h_z)$ へと持ち上げる。これは、補助変数 $t$（膨張作用を表す）を導入し、斉次ポアソン構造 $\Pi = \frac{1}{t}\Lambda + \frac{\partial}{\partial t} \wedge E$ および 1-斉次ハミルトニアン $\hat{H} = tH$ を定義することによって達成される。ヤコビ力学は、$P$ 上の斉次ポアソン力学の射影として回収される。
2. 斉次シンプレクティック・バイリアリゼーション（Homogeneous Symplectic Bi-realization）：ポアソニゼーションされた系は、斉次シンプレクティック・バイリアリゼーションを用いて、さらにシンプレクティックな設定へと持ち上げられる。これには、余接束を多様体に接続し、かつポアソン構造と斉次性を保存する写像 $\alpha, \beta: U \subset T^*P \to P$ の構築が含まれる。ハミルトニアン流は、このシンプレクティック・グロイド内のラグランジュ・ビセクションとして表現される。
3. マグヌス展開による反復構成：積分器は、これらのラグランジュ・ビセクションの生成関数を近似することによって構築される。時間依存のハミルトニアンの場合、著者らはマグヌス展開を用いて関数の級数 $S_i$ を生成する。フローはハミルトニアンベクトル場の単一の指数関数として近似され、これにより基礎となるリー代数の構造が保証される。
4. 再帰的アルゴリズム：明示的な再帰的アルゴリズムにより、任意の次数 $k$ の係数 $S_i$ を生成する。この手法は、写像 $\alpha$ を通じて中間点 $y_n$ を解き、計算された $S_i$ 関数の勾配を用いて状態 $x_{n+1}$ を更新するプロセスを含む。

主な貢献

• 体系的な枠組み：本論文は、ヤコビ・ハミルトニアン系を斉次ポアソン系に変換する一般的な手順を確立し、既存のポアソン積分技術をより広いヤコビの文脈へ適用することを可能にした。
• 任意次数のアルゴリズム：マグヌス展開と厳密な斉次ラグランジュ・ビセクションに基づき、任意の次数 $k$ のJHIを構築できる明示的な再帰的アルゴリズムを導入した。
• 構造保存：得られた積分器は、斉次ポアソン構造、ポアソニゼーションされたハミルトニアンのレベルセット、およびヤコビ・カシミール関数を、当該手法の次数まで保存することが証明されている。
• 誤差の逆解析（Backward Error Analysis）：誤差の逆解析を通じて、JHIが修正された時間依存ヤコビ・ハミルトニアン系の厳密なフローと一致することを示す。これにより、積分器が元のハミルトニアン $H$ の（ヤコビ系では保存されないが、持ち上げられた斉次系では保存される）正しい散逸挙動を含む、長期的な定性的忠実性を継承することが保証される。

結果および数値実験
提案された手法は、低次元の例示的問題から古典的なモデルに至るまで、多様な例でテストされた：

• コンタクト・ハミルトニアン系：3次元コンタクト系において、第1次JHIは標準的な2次ルンゲ・クッタ法（RK-2）と比較して座標の「爆発」を遅延させ、形式的な次数は低いものの優れた定性的挙動を示した。
• 厳密解：二次ハミルトニアンを持つ特定の2次元ヤコビ問題において、構築されたJHIは厳密解を再現し、高次の補正項はすべて消失した。
• 収束率：2次元、3次元、および4次元のヤコビ問題に関する数値実験により、理論的な収束次数が確認された。特に、近似的なシンプレクティック・バイリアリゼーションを使用した場合（4次元の場合）でも、手法は2次の収束を達成した。
• 古典的モデル：
  • 減衰調和振動子：JHIは、軌道とハミルトニアン誤差の特性を保持する点で、セミインプリシット・シンプレクティック・オイラー法よりも優れた性能を示した。
  • ロトカ・ヴォルテラ系：JHIは修正された系の閉軌道構造を正常に保持したが、RK-2は軌道の閉包を維持できなかった。
  • 剛体回転：線形双ベクトルと二次ベクトル場を用いた3次元剛体回転モデルに適用された。近似的なバイリアリゼーションの使用により、不安定な構成において厳密解からのわずかな偏差が生じたものの、JHIは閉軌道を維持し、ハミルトニアン誤差を有界に保ち、長期的な幾何学的整合性においてRK-2を上回った。

意義と主張
本論文は、ヤコビ・ハミルトニアン積分器が、幾何学的積分技術の自然な拡張であることを主張している。主な意義は、標準的な手法が失敗しやすい散逸的および時間依存的なシステムの固有の幾何学的構造を保存する積分器を構築できる点にある。著者らは、JHIが単に解を近似するのではなく、修正されたハミルトニアン系の厳密なフローを提供することで、エネルギー誤差の有界性と正しい長期挙動を保証することを強調している。本研究は、数値実験を通じて、厳密なシンプレクティック・バイリアリゼーションを解析的に得ることが困難な場合であっても、JHIが強い幾何学的整合性と精度を備えていることを実証している。この枠組みは、適応的時間刻み幅や、より大規模な高次元非保存系への応用に関する将来の研究への基礎として提示されている。