[論文レビュー] Symplectic integration of deviation vectors and chaos determination. Application to the Hénon-Heiles model and to the restricted three-body problem

Anne-Sophie Libert, Charles Hubaux+1 more|arXiv (Cornell University)|May 31, 2010

Quantum chaos and dynamical systems参考文献 22被引用数 5

ひとこと要約

本稿では、ハミルトニアン系における軌道と乖離ベクトルを、シンプレクティックスキームを用いて同時に統合するグローバルシンプレクティック積分法を提案する。この手法により、カオス検出の精度が顕著に向上する。本手法はエネルギー保存を維持し、大きな時間刻み（非シンプレクティック手法の4倍まで）を用いても、正則軌道およびカオス的軌道を正しく特定可能であり、CPU時間を削減し、ヘノン＝ハイレス系および制限3体問題におけるカオス検出の誤検出を回避する。

ABSTRACT

In this work we propose a new numerical approach to distinguish between regular and chaotic orbits in Hamiltonian systems, based on the simultaneous integration of both the orbit and the deviation vectors using a symplectic scheme, hereby called {\em global symplectic integrator}. In particular, the proposed method allows us to recover the correct orbits character with very large integration time steps, small energy losses and short CPU times. To illustrate the numerical performances of the global symplectic integrator we will apply it to two well-known and widely studied problems: the Hénon-Heiles model and the restricted three-body problem.

研究の動機と目的

• 長時間にわたる積分において、ハミルトニアン系の正則軌道とカオス的軌道を信頼性高く区別する課題に対処すること。
• SALIなどのカオス検出指標の精度を向上させるために、シンプレクティックスキームを用いて軌道と乖離ベクトルを同時に統合すること。
• ヘノン＝ハイレスモデルや制限3体問題のようなカオス的系の長期的シミュレーションにおいて、計算コストとエネルギー損失を低減すること。
• シンプレクティック積分を用いた乖離ベクトルの統合が、動的忠実性を保ち、カオス的振る舞いを正しく同定する点で、非シンプレクティック手法を上回ることを示すこと。

提案手法

• 本手法は、対称なシンプレクティック積分法（SABA_n）を用い、主軌道と2つの乖離ベクトルを同時に統合することで、系のシンプレクティック構造を保存する。
• H(x) = A(p) + B(q) というハミルトニアンの分割を用い、SABA_10 などのシンプレクティック積分法を高精度な長期的統合に活用できる。
• グローバルシンプレクティック積分法は、大きな時間刻みでも小さなエネルギー損失と安定した数値的挙動を保証し、長期的シミュレーションに不可欠である。
• 乖離ベクトルは、d𝐯/dt = M(φ(t))𝐯 という接線写像に従い、M はハミルトニアンベクトル場のヤコビ行列である。このベクトルは、軌道と同一のシンプレクティックスキームで統合される。
• SALIカオス指標を用いて手法を検証し、乖離ベクトルの進化を追跡することで、正則運動とカオス的運動を区別する。
• Bulirsch-Stoer などの非シンプレクティックスキームと比較し、エネルギーの著しいドリフトを示し、小さな時間刻みでもカオスを正しく同定できないことが判明した。

実験結果

リサーチクエスチョン

• RQ1 軌道と乖離ベクトルの両方をシンプレクティックに積分することで、ハミルトニアン系におけるカオス検出の精度が向上するか？
• RQ2 精度やエネルギー保存を損なわずに、どの程度大きな時間刻みをカオス検出に使用できるか？
• RQ3 長時間にわたる積分において、シンプレクティック積分法と非シンプレクティック積分法は、正則軌道とカオス的軌道をどの程度正確に区別できるか？
• RQ4 グローバルシンプレクティック積分法は、ヘノン＝ハイレス系や制限3体問題のような分離子の影響を正しく同定できるか？

主な発見

• ヘノン＝ハイレスモデルにおいて、非シンプレクティック手法に比べて4倍の時間刻みを用いても、ほぼすべての軌道の性質を正しく同定できる。
• 非シンプレクティック積分法（Bulirsch-Stoer など）は、小さな時間刻みでも著しいエネルギー損失を示し、正則軌道とカオス的軌道を区別できない。
• シンプレクティック法は、100,000年間の積分においてもエネルギー保存を10^-12のオーダーで維持し、信頼性の高い長期的シミュレーションを可能にする。
• 適応的非シンプレクティックスキーム（Bulirsch-Stoer など）と比較して、CPU時間を最大3倍まで削減しながらも、精度を保持する。
• 制限3体問題における分離子の影響は、シンプレクティック積分法を用いることで明確に同定可能であるのに対し、非シンプレクティック手法ではカオス的振るまいを検出できない。
• グローバルシンプレクティック積分法により、大きな安定した時間刻みを用いて長時間積分が可能となり、制限3体問題における長期的共振の正確な検出が可能になった。