Spectral Koopman Approach for Reconstructing State-space Geometry of Cislunar Restricted 3-Body Problem

技術要約：シスルナ制限三体問題の位相空間幾何学再構成のためのスペクトル・クープマン・アプローチ

問題提起
地球、月、太陽の重力相互作用によって支配されるシスルナ（月近傍）環境は、リブラシオン点、不変多様体、およびカオス的輸送を特徴とする高度に非線形な力学系である。従来のミッション設計および軌道解析は、円制限三体問題（CR3BP）と数値最適化手法（例：微分補正、マルチシューティング）に依存している。しかし、これらの既存手法には、以下の3つの持続的な限界が存在する。

1. 忠実度の不一致： 簡略化されたモデルと高忠実度な現実との間の乖離。
2. 計算コスト： ナビゲーションの更新が発生した際に、再計算が必要となるコスト。
3. スケーラビリティの問題： 長期ホライゾンの不確実性伝播におけるモンテカルロ法のスケーラビリティの低さ。

さらに、標準的な局所線形化手法は、共鳴や輸送チャネルのような強い非線形構造によって組織化されたグローバルな位相空間の幾何学を捉えることができない。一方、作用素論的手法（特にクープマン作用素）は、非線形系のグローバルな線形表現を提供するが、既存のデータ駆動型近似（例：拡張動的モード分解）は、慎重な基底選択を必要とし、固有のグローバルスペクトルを正確に回復することに苦慮する場合が多い。

手法
著者らは、アドホックな辞書選択に頼ることなく、平面CR3BPのグローバルな位相空間幾何学を再構成するための、パス積分定式化に基づくクープマン・スペクトルを用いた新しいアプローチを提案している。

1. クープマン作用素の枠組み： システムの動力学 $\dot{x} = f(x)$ は、観測量の進化を通じて解析される。クープマン作用素 $U_t$ は、非線形フローの線形な記述を提供する。主要な固有関数 $\phi_\lambda$ は、作用素の方程式 $\frac{\partial \phi_\lambda}{\partial x}f = \lambda \phi_\lambda$ を満たす。
2. パス積分定式化： 作用素を有限次元の部分空間で近似する代わりに、著者らは主要な固有関数を直接計算する。ベクトル場を線形成分 ($Ax$) と非線形成分 ($F_n(x)$) に分解することで、非線形部分の固有関数$h_\lambda(x)$ は、次のようなパス積分として表現される：
   $$h_\lambda(x) = \int_0^\infty e^{-\lambda \tau} w_\lambda^\top F_n(s_\tau(x)) d\tau$$
   ここで、$s_\tau(x)$ はシステムの軌跡を表す。
3. 数値実装： 固有関数は、有限の時間ホライゾン $T$ にわたる積分の評価によって近似される。不安定な固有関数（$\lambda_u > 0$ に関連するもの）は前方軌跡を用いて計算され、安定な固有関数（$\lambda_s < 0$ に関連するもの）は後方軌跡を用いて計算される。
4. 幾何学の再構成： 著者らは、これらの主要な固有関数の同時ゼロレベル曲線が、系の平衡点（ラグランジュ点）の不変多様体（安定、不安定、および中心多様体）に対応するという性質を利用している。

主な貢献

• グローバルな線形表現： 本論文は、クープマンの枠組みが、局所線形化の限界を克服し、CR3BPのグローバルな線形表現を生成できることを示している。
• 固有関数の直接計算： パス積分アプローチを活用することで、本手法はEDMDのような手法で必要とされる手動の辞書設計を回避し、時系列データから主要な固有関数を直接計算する。
• 不変多様体の特定： 本研究は、計算された固有関数のゼロレベル集合が位相空間の特性、具体的には5つのラグランジュ点（$L_1$ から $L_5$）に関連する安定、不安定、および中心多様体を符号化していることを確立している。
• 輸送チャネルの可視化： 本アプローチは、ある平衡点の不安定多様体が別の平衡点の安定多様体とどのように整列して輸送チャネルを形成するかを示すことで、シスルナ環境の「見えない高速道路」を可視化することに成功している。

結果
提案された手法は、平面4次元CR3BPモデル（地球・月系）に適用された。シミュレーション結果には以下が含まれる：

• ゼロ速度曲線（ZVCs）の再構成： 本手法は、ヤコビ定数によって定義されるエネルギー的に許容された領域と禁止領域を回復する。
• 平衡点における局所幾何学：
  • $L_1$ および $L_2$： 固有関数は輸送ゲートウェイを明確に描き出している。不安定多様体は「ネック」部を通じて外側へ延び（地球からの脱出や外部への進入を可能にする）、安定多様体は内側へと曲がっている（軌道を捕捉する）。
  • $L_3$： 外部領域へと漏れ出す脆弱な多様体を示しており、強い輸送回廊を形成していない。
  • $L_4$ および $L_5$： 中心固有関数が支配的であり、これらの方形点における線形安定性と一致する、トロポイド軌道やタドポール軌道に対応する閉じたループを生成する。
• グローバルな位相空間の分割： 固有関数のゼロレベル曲線は、位相空間をスペクトル的に分割する。エネルギーに基づくヤコビ面が運動の「場所」を定義するのに対し、クープマン曲線は、軌道がその領域内でどのように漸近的に組織化されているか（セパラトリクスや不変境界の特定）を明らかにする。
• 輸送メカニズム： 結果は、$L_1$ と $L_2$ の間の輸送チャネルを示しており、$L_1$ の不安定多様体と $L_2$ の安定多様体が不変チューブを形成し、地球領域と月領域の間の移動を促進している様子を描いている。

意義および主張
本論文は、スペクトル・クープマン・アプローチが、個々の軌道から観測量の進化へと焦点を移すことにより、シスルナ力学に対する強力なツールを提供すると主張している。主な意義は、以下の能力にある：

1. グローバルな幾何学の回復： ローカルな手法では捉えきれない非線形構造を捉えつつ、時系列データからCR3BPのグローバルな位相空間幾何学を再構成する。
2. ダウンストリーム・タスクの実現： 不確実性伝播、到達可能性解析、および制御設計のための基礎を提供する。クープマン作用素は線形であるため、著者らは、非線形なCR3BPに対する不確実性解析や制御戦略の設計が、従来の非線形手法と比較して簡素化され、スケーラビリティを高められると断じている。
3. データ駆動型の発見： 宇宙力学におけるデータ駆動型の発見に向けた新しい方向性を切り開く。作用素の線形性により、本質的に非線形なシステムに対する制御戦略の分析および設計が可能となる。

著者らは、現在の研究はCR3BPの幾何学の回復に焦点を当てているものの、この枠組みは、複雑な多体環境におけるリアルタイムの軌道最適化やロバスト制御への将来的な応用に対して有望であると結論付けている。