[论文解读] 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

一句话总结

本文研究了在 Hénon-Heiles 型势能中，直线振子轨道的同频叉子分岔级联，从理论上证明了分岔能量几何级数的标度常数 δ 为 δ = e^{π/\sqrt{2}} ≈ 9.2206，且在两个方向上的空间标度常数相同，均为 α = β = √δ。自相似性与解析标度行为源于临界能量 e = 1 处的谐振鞍点结构及横向曲率。

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.

研究动机与目标

• 分析能量趋近于鞍点能量 e = 1 时，Hénon-Heiles 型势能中同频叉子分岔的级联行为。
• 确定在连续分岔中诞生的周期轨道自相似结构是否表现出相同的解析标度常数。
• 从势能参数（特别是横向曲率 ω⊥）出发，解析推导标度常数 δ。
• 将此标度行为与面积保测映射中的标准 Feigenbaum 场景进行对比，后者中 δ、α 和 β 不同且具有普适性。

提出的方法

• 通过能量 e 参数化 Hénon-Heiles 势能的缩放牛顿运动方程，使用数值积分与牛顿-拉夫森迭代法定位周期轨道。
• 分析周期轨道的稳定性矩阵，并追踪周期 TA 的轨迹 tr M(TA) 以识别分岔点。
• 利用 libration 周期 TA 在 e = 1 附近的渐近行为，通过对数发散关系 TA ∼ √(2 ln(64/ǫ))（其中 ǫ = 1 − e）推导其渐近形式。
• 结合 tr M(TA) 的渐近周期性 ∆T = 2π/ω⊥ 与能量依赖关系 TA ∼ d ln[c/(1 − e)]，推导出标度常数 δ = e^{2π/(ω⊥d)}。
• 通过分岔能量 en 的数值迭代验证解析结果，得到 δ ≈ 9.2203，标准差极低。
• 应用特殊函数（雅可比椭圆函数）描述分岔附近的微小振幅振荡，实现与数值结果的高精度匹配。

实验结果

研究问题

• RQ1 Hénon-Heiles 型势能中，控制分岔能量几何级数的标度常数 δ 的解析表达式是什么？
• RQ2 在连续分岔中生成的周期轨道自相似结构，在两个空间方向上是否表现出相同的标度常数 α 和 β？
• RQ3 该系统中的标度行为与二维面积保测映射中的标准 Feigenbaum 场景有何不同？
• RQ4 标度常数 δ 是否能从势能参数（特别是 ω⊥ 和 d）出发进行解析推导？
• RQ5 横向曲率 ω⊥ 与 Maslov 指数在决定鞍点附近渐近动力学中的作用是什么？

主要发现

• 分岔能量几何级数的标度常数 δ 被解析推导为 δ = e^{π/\sqrt{2}} ≈ 9.2206125，且获得高精度数值验证（δ ≈ 9.2203 ± 0.02）。
• 两个方向上的空间标度常数相同：α = β = √δ ≈ 3.0365，与标准 Feigenbaum 场景中 α ≠ β 的情况形成鲜明对比。
• librating 轨道的渐近周期 TA 随 e → 1 时对数发散，满足 TA ∼ √(2 ln(64/ǫ))（其中 ǫ = 1 − e），从而导出能量标度 δ = e^{π/\sqrt{2}}。
• 周期轨道尖端的自相似性仅在两代之间显现，这是由于分岔中诞生的轨道交替呈现稳定与不稳定特性，因此需通过 α² = δ 对连续代际进行缩放。
• 新不稳定轨道（e > 1 时）的 Lyapunov 指数 τ 满足在 e → 1 时 χ = 2πω∥/ω⊥，且 σ = χ/Tτ = ω∥。
• 使用雅可比椭函数 Ecm_p(at) 与 Esm_p(at) 能够精确描述分岔附近的微小振幅振荡，与数值结果高度一致。