[論文レビュー] Detecting Order and Chaos by the Linear Dependence Index (LDI) Method

Chris G. Antonopoulos, Tassos Bountis|ArXiv.org|Nov 2, 2007

Quantum chaos and dynamical systems参考文献 18被引用数 13

ひとこと要約

本稿では、多次元ハミルトニアン系におけるカオスの検出のための、一般化されたアライメント指数（GALI）の計算的に効率的な代替手法として、線形従属インデックス（LDI）法を導入する。特異値分解（SVD）を用いて発展ベクトルに適用することで、LDIはGALIと同一の数値的結果を達成するが、特に高次元系においてはCPU時間を著しく削減する。

ABSTRACT

We introduce a new methodology for a fast and reliable discrimination between ordered and chaotic orbits in multidimensional Hamiltonian systems which we call the Linear Dependence Index (LDI). The new method is based on the recently introduced theory of the Generalized Alignment Indices (GALI). LDI takes advantage of the linear dependence (or independence) of deviation vectors of a given orbit, using the method of Singular Value Decomposition at every time step of the integration. We show that the LDI produces estimates which numerically coincide with those of the GALI method for the same number of $m$ deviation vectors, while its main advantage is that it requires considerable less CPU time than GALI especially in Hamiltonian systems of many degrees of freedom

研究の動機と目的

• 多次元ハミルトニアン系における秩序的軌道とカオス的軌道を識別する、より高速で効率的な手法の開発。
• 各時間ステップで多数の行列式を計算する必要があるGALI法の計算的ボトル neck を克服すること。
• 発展ベクトルの線形従属をカオスの兆候として利用し、SVDを用いてリアルタイム計算を可能にすること。
• さまざまな自由度とエネルギー領域において、LDI法とGALI法の妥当性を検証すること。
• GALI法が時間的・計算的コストが膨大になり、実用的でなくなる高次元系において、LDI法の優位性を示すこと。

提案手法

• LDI法は、位相空間における$m$個の発展ベクトルからなる$2N \times m$行列の特異値分解（SVD）から得られる特異値の積を計算する。
• これらの特異値の時間発展を追跡することで、発展ベクトルが線形独立のままである（秩序的運動）か、線形従属に収束する（カオス的運動）かを判定する。
• カオス的軌道では最大特異値が支配的となり、特異値の積は指数関数的に減少し、GALIインデックスと類似した振る舞いを示す。
• 順序$m$のLDIインデックスは、$m$個の特異値の積として定義され、数値的に対応するGALI$_m$インデックスと一致する。
• 数値積分の各時間ステップで本手法を適用することで、動的挙動のリアルタイム検出が可能になる。
• 理論的裏付けは、発展ベクトルの線形従属が最大リャプノフ方向に一致することに起因し、これはカオスの特徴的兆候である。

実験結果

リサーチクエスチョン

• RQ1 LDI法は、GALI法と同等の精度でハミルトニアン系におけるカオスを信頼性高く検出できるか？
• RQ2 自由度の増加に伴い、LDI法の計算コストはGALI法と比べてどの程度異なるか？
• RQ3 LDI法はGALI法の物理的洞察を保持しつつ、より迅速な数値的評価を可能にするか？
• RQ4 LDI法は低次元系および高次元系の両方において、秩序的軌道とカオス的軌道を正確に区別できるか？
• RQ5 自由度$N \gg 10$の系において、LDI法とGALI法のCPU時間のスケーリング特性はいかなるものか？

主な発見

• 同じ発展ベクトル数$m$と基準軌道に対して、LDI法はGALI法と同一の数値的結果を生成する。
• 5自由度のFPU格子系では、LDI 5とGALI 5が$t = 80$までの計算にそれぞれ約1.5秒を要し、低次元系では同等の性能を示した。
• 15自由度の系では、GALI 8の計算に約186秒を要したが、LDI 8は同じ時間間隔でわずか1秒で完了し、186倍の高速化が達成された。
• LDI 8とGALI 8の時間発展は、同じ指数関数的減少則$\propto e^{-0.385t}$に従い、両者の動的挙動が同一であることを確認した。
• LDI法は、$E = 2$における不安定なSPO1軌道付近のFPUモデルでカオス的運動を的確に同定し、リャプノフ指数が正の値に収束することを示した。
• GALI法が各時間ステップで数百万個の行列式を計算する必要があるため、高次元系では計算的に非現実的になるが、LDI法はそのような状況において特に優位である。