[論文レビュー] Bifurcation cascades and self-similarity of periodic orbits with analytical scaling constants in Henon-Heiles type potentials

M. Brack|arXiv (Cornell University)|Jun 21, 2000

Quantum chaos and dynamical systems参考文献 34被引用数 38

ひとこと要約

本稿は、ヘノン＝ハイレース型ポテンシャルにおける直線的リブレーティング軌道の等時的ピッチfork分岐カスケードを調査し、臨界エネルギー e = 1 における調和的サドル構造と横方向の曲率に起因して、自己相似性と解析的スケーリングが生じることを示している。分岐エネルギーの等比数列のスケーリング定数 δ は、δ = e^{π/\sqrt{2}} ≈ 9.2206 として解析的に導出され、両方向における空間的スケーリング定数も α = β = √δ で同一である。

ABSTRACT

We investigate the isochronous bifurcations of the straight-line librating orbit in the Henon-Heiles and related potentials. With increasing scaled energy e, they form a cascade of pitchfork bifurcations that cumulate at the critical saddle-point energy e=1. The stable and unstable orbits created at these bifurcations appear in two sequences whose self-similar properties possess an analytical scaling behavior. Different from the standard Feigenbaum scenario in area preserving two-dimensional maps, here the scaling constants \alpha and \beta corresponding to the two spatial directions are identical and equal to the root of the scaling constant \delta that describes the geometric progression of bifurcation energies e_n in the limit n --> infinity. The value of \delta is given analytically in terms of the potential parameters.

研究の動機と目的

• ヘノン＝ハイレース型ポテンシャルにおける等時的ピッチfork分岐カスケードを、臨界エネルギー e = 1 に近づくに従って分析すること。
• 連続する分岐で生成された周期軌道の自己相似的構造が、同一のスケーリング定数 α と β を示すかどうかを特定すること。
• 特に横方向曲率 ω⊥ に依存するポテンシャルパラメータの観点から、スケーリング定数 δ を解析的に導出すること。
• 面積保存写像における標準的フェイゲンバック・スケーリングと比較し、δ、α、β がそれぞれ異なる普遍的定数である標準的フェイゲンバック状況とは対照的に、本系におけるスケーリング行動を明らかにすること。

提案手法

• エネルギー e を変数とするヘノン＝ハイレースポテンシャルのスケーリングされたニュートン運動方程式を数値積分とニュートン・ラプソン反復法を用いて解き、周期軌道を特定する。
• 周期軌道の安定性行列を解析し、周期 TA に関するトレース tr M(TA) を追跡することで分岐点を同定する。
• e = 1 に近い領域におけるリブレーティング周期 TA の漸近的挙動を、TA ∼ √(2 ln(64/ǫ))（ǫ = 1 − e）として導出する。
• tr M(TA) の漸近的周期性 ∆T = 2π/ω⊥ と、TA ∼ d ln[c/(1 − e)] におけるエネルギー依存性を組み合わせ、スケーリング定数 δ = e^{2π/(ω⊥d)} を導出する。
• 分岐エネルギー en の数値的反復により解析的結果を確認し、δ ≈ 9.2203（標準偏差が小さい）と得られた。
• 小振幅振動を記述するための特殊関数（ラマ函数）を用い、数値結果と高い精度で一致するように検証した。

実験結果

リサーチクエスチョン

• RQ1 ヘノン＝ハイレース型ポテンシャルにおける分岐エネルギーの等比数列のスケーリング定数 δ の解析的形は何か？
• RQ2 連続する分岐で生成された周期軌道の自己相似的構造は、両方向で同一のスケーリング定数 α と β を示すか？
• RQ3 本系におけるスケーリング行動は、2次元面積保存写像における標準的フェイゲンバック状況とどのように異なるか？
• RQ4 スケーリング定数 δ は、特に ω⊥ と d といったポテンシャルパラメータから解析的に導出可能か？
• RQ5 横方向曲率 ω⊥ とマスロフ指数は、サドル点付近の漸近的力学を決定づける役割を果たすか？

主な発見

• 分岐エネルギーの等比数列のスケーリング定数 δ は、解析的に δ = e^{π/\sqrt{2}} ≈ 9.2206125 として導出され、数値的確認（δ ≈ 9.2203 ± 0.02）がなされた。
• 両方向における空間的スケーリング定数は同一である：α = β = √δ ≈ 3.0365 であり、標準的フェイゲンバック状況における α ≠ β とは対照的である。
• リブレーティング軌道の漸近的周期 TA は、TA ∼ √(2 ln(64/ǫ))（ǫ = 1 − e）として対数発散する。この結果、エネルギースケーリング δ = e^{π/\sqrt{2}} が得られる。
• 周期軌道の先端における自己相似性は、分岐で生成された軌道が安定／不安定を交互に繰り返す性質のため、2世代にわたってのみ顕在化する。このため、連続する世代間のスケーリングには α² = δ が必要となる。
• e > 1 の新たな不安定軌道のリャプノフ指数 τ は、e → 1 の極限で χ = 2πω∥/ω⊥ に近づき、σ = χ/Tτ = ω∥ となる。
• 小振幅振動を記述するためのラマ函数 Ecm_p(at) と Esm_p(at) の使用により、分岐付近の振動を高精度で記述でき、数値結果と一致した。