A Coordinate System for Dynamical Instabilities in Hierarchical Systems in REBOUND

技術的サマリー：REBOUND における階層系内の動力学的不安定性のための座標系

問題提起
$N$ 体問題は、特にギガ年スケールでの惑星系のシミュレーションにおいて、古典力学における中心的な課題であり続けています。ウィズダム・ホーマン（WH）ハイブリッド積分器は、ハミルトニアンを支配的なケプラー項と摂動的な相互作用項に分割することで長期的な積分のゴールドスタンダードとなっていますが、その有効性は座標系の選択に大きく依存します。

• ヤコビ座標は、2 体問題を厳密に解くことで高い精度を提供しますが、柔軟性に欠けます。動力学的不安定性により軌道の順序が変化した場合には機能しません。
• **民主的日心座標（DHC）**は、軌道交差を許容し、ハイブリッド積分器（例：MERCURY、SyMBA、TRACE）の標準となっていますが、単一の支配的な中心質量を持たない系、すなわち惑星系が連星の一方の恒星を周回する広連星系のような環境では性能が劣ります。

広連星系において、DHC は効率的な積分に必要な階層性を維持できず、一方ヤコビ座標は、これらの環境で一般的に見られる接近遭遇や軌道的不安定性を処理できません。惑星、宿主星、伴星の任意の 2 体間の接近遭遇を支援しつつ、WH 手法の速度の利点を維持する座標系が必要です。

手法
著者らは、J. E. Chambers ら（2002 年）によって導出された座標系、通称「広連星座標（WB）」を、オープンソースパッケージREBOUND内のTRACEハイブリッド積分器に実装しました。

1. 座標変換:
   WB 系は、中心恒星の重心運動、伴星の相対運動、および宿主星に対する惑星の相対運動を分離する新しい正準座標（$\vec{Q}$）と運動量（$\vec{P}$）を定義します。これにより、ハミルトニアンを以下の 3 つの成分に分割できます。

   • $H_{Jump}$：中心恒星の重心運動。
   • $H_{Interaction}$：惑星 - 惑星間および惑星 - 伴星間の相互作用。
   • $H_{Kepler}$：宿主星周りを周回する惑星および伴星のケプラー運動。
     伴星が存在しない場合、これらの座標は標準的な DHC に帰着します。

2. 積分スキーム:
   著者らは、TRACE の時間可逆的なスイッチング機構を WB 座標に適応させました。積分器は、$H_{Interaction}$、$H_{Jump}$、$H_{Kepler}$ を順次進めるサブステップ（リープフロッグ方式）で動作します。

   • スイッチング論理: 接近遭遇（惑星 - 惑星、惑星 - 恒星、または惑星 - 伴星）が発生すると、積分器は動的に $H_{Interaction}$ から項を $H_{Kepler}$ へ移動させ、ケプラー項が支配的であることを保証します。
   • 遭遇処理:
     • 惑星が宿主星または伴星に臨界半径（デフォルト：伴星に対して 0.1 Hill 半径）以内まで接近した場合、系は慣性座標系に戻り、運動方程式全体を解くために従来の積分器（Bulirsch-Stoer または IAS15）を使用します。
     • 惑星同士が接近する場合、特定の相互作用項を $H_{Kepler}$ へ移動させ、Bulirsch-Stoer 積分器を用いて解きます。一方、残りの系は WH ソルバーで継続します。

3. ベンチマーク:
   WB 座標における TRACE の性能は、以下のものに対してベンチマークされました。

   • IAS15: 高精度な適応時間ステップ積分器（真値として機能）。
   • DHC における TRACE: 比較のための標準実装。
   • ヤコビ座標における WHFast: 安定領域での比較に使用。

主要な結果

• パラメータ空間分析: 伴星の質量（$1 M_J$ から $3 M_\odot$）と半長径（50–1000 AU）を変化させた 10,000 回のシミュレーションの調査により、WB 座標はより大きく遠方の伴星に対して DHC を大幅に上回ることが明らかになりました。著者らは、WB 座標が必要となる時期を判断するためのべき乗則のガイドライン（式 11）を提供します。
  $$\log_{10}\left(\frac{M}{M_\odot}\right) \approx 0.9 \log_{10}\left(\frac{a}{1 \text{ AU}}\right) - 4.67$$
  WB 座標は、検討されたパラメータ空間全体において、ヤコビ座標における WHFast と競争力のある性能を発揮します。

• 惑星 - 惑星散乱: 3 つの木星質量惑星と 0.5 $M_\odot$ の伴星を含むシミュレーションにおいて、DHC における TRACE は IAS15 と比較して統計的に誤った結果（半長径と軌道傾斜角の過大評価、離心率の過小評価）を生成しました。一方、WB 座標における TRACE は IAS15 と統計的に類似した結果を回復しつつ、9 倍の高速化を達成しました。

• 恒星のフライバイ: 海王星のような惑星に対する恒星のフライバイをシミュレーションしたテストにおいて、DHC における TRACE はエネルギー、離心率、軌道傾斜角の変化に関する解析的予測を再現できませんでした。WB 座標における TRACE は、解析的予測と IAS15 の結果の両方に一致し、IAS15 に対して2 倍の高速化を提供しました。

• ZLK 振動: von Zeipel-Lidov-Kozai（ZLK）サイクルにおいて、固定時間ステップ積分器は、高離心率相で到達する極端な離心率に対してしばしば困難をきたします。

  • DHC では、TRACE が IAS15 の精度に一致する有効なパラメータ空間は極めて限定的でした。
  • WB 座標では、TRACE は連星の軌道傾斜角が $70^\circ$ まで（$e_{max} \sim 0.9$）の範囲で IAS15 の精度に一致しました。
  • 内側に地球型惑星が存在する特定のテストケース（小さな大域的時間ステップを強制）において、WB 座標における TRACE は精度を維持しつつ IAS15 よりも20 倍高速であり、一方 DHC は正しい振動周期を再現できませんでした。

意義とガイドライン
本論文は、TRACE への WB 座標の実装が、従来のハイブリッド手法が失敗する特定の領域、すなわち動力学的不安定性と接近遭遇に陥りやすい広連星系に対処していると主張しています。著者らは、連星系における混沌とした系の大規模な集合に対して、この手法が精度と計算の実用性のバランスを取る唯一の実用的な積分アプローチである可能性を強調しています。

著者らは、REBOUND における積分器の選択に関する以下のガイドラインを提供しています。

1. IAS15: 計算時間が制約とならない場合、または正確な軌跡が必要な場合に使用します。
2. WHFast: 計算時間が懸念される場合、かつ接近遭遇が予想されない場合に使用します。
3. TRACE: 計算時間が懸念される場合、接近遭遇が予想される場合、かつ定性的に正しい軌跡で十分である場合に使用します。
   • 座標の選択: 太陽型恒星を周回する系については、式 11 を用いて DHC と WB の間で判断します。真の恒星伴星に対しては WB が強く推奨されます。

この研究は、WB 座標における TRACE が、高精度な適応積分器に伴う prohibit 的な計算コストなしに、連星系における動力学的不安定性、惑星形成、および恒星フライバイの研究に十分な精度を提供することを示しています。