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

一句话总结

本文提出了一种全局辛积分器，通过辛格式同时积分哈密顿系统中的轨道与偏差向量，显著提升了混沌检测的准确性。该方法在使用大时间步长（最大可达非辛方法的四倍）时仍能保持能量守恒，并正确识别规则轨道与混沌轨道，同时减少CPU计算时间，并避免在Hénon-Heiles模型和限制性三体问题中出现混沌检测的误报。

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等混沌检测指标的准确性。
• 在Hénon-Heiles模型和限制性三体问题等混沌系统的长期模拟中，降低计算成本与能量损失。
• 证明偏差向量的辛积分优于非辛方法，在保持动力学保真度和正确识别混沌行为方面表现更优。

提出的方法

• 该方法采用对称辛积分器（SABA_n）同时积分主轨道与两个偏差向量，以保持系统的辛结构。
• 采用哈密顿量分裂形式 H(x) = A(p) + B(q)，从而可应用如SABA_10等辛积分器实现高精度的长期积分。
• 全局辛积分器即使在大时间步长下也能保持微小的能量损失与稳定的数值行为，这对长期模拟至关重要。
• 偏差向量通过切向映射 d𝐯/dt = M(φ(t))𝐯 进行演化，其中 M 为哈密顿向量场的雅可比矩阵，并使用与轨道相同的辛格式进行积分。
• 通过SALI混沌指标对方法进行验证，该指标通过追踪偏差向量的演化来区分规则运动与混沌运动。
• 将该方法与Bulirsch-Stoer等非辛格式进行对比，后者虽使用更小的时间步长，但表现出显著的能量漂移，且无法正确识别混沌行为。

实验结果

研究问题

• RQ1 同时对轨道与偏差向量进行辛积分是否能提升哈密顿系统中混沌检测的准确性？
• RQ2 在不损害准确性或能量守恒的前提下，混沌检测中可使用多大的时间步长？
• RQ3 在长时间积分中，辛积分器与非辛方法在区分规则与混沌轨道方面的性能差异如何？
• RQ4 全局辛积分器是否能正确识别Hénon-Heiles模型与限制性三体问题中分离子的影响？

主要发现

• 该全局辛积分器即使在时间步长为非辛方法四倍的情况下，也能正确识别Hénon-Heiles模型中几乎所有轨道的性质。
• 非辛积分器如Bulirsch-Stoer即使使用更小的时间步长，仍产生显著的能量损失，且无法正确区分规则与混沌轨道。
• 该辛方法在长达100,000年的积分中能量守恒精度达到10^-12，确保了长期模拟的可靠性。
• 与自适应非辛格式如Bulirsch-Stoer相比，该方法将CPU时间减少至最多三分之一，同时保持精度。
• 在限制性三体问题中，该辛积分器能清晰识别分离子的影响，而传统非辛方法则无法检测到混沌行为。
• 该全局辛积分器通过支持使用大而稳定的长时间步长，实现了对限制性三体问题中长期共振的准确检测。