Comparative study of variational chaos indicators and ODEs' numerical integrators

Abstract

The reader can find in the literature a lot of different techniques to study the dynamics of a given system and also, many suitable numerical integrators to compute them. Notwithstanding the recent work of Maffione et al. (2011a) for mappings, a detailed comparison among the widespread indicators of chaos in a general system is still lacking. Such a comparison could lead to select the most efficient algorithms given a certain dynamical problem. Furthermore, in order to choose the appropriate numerical integrators to compute them, more comparative studies among numerical integrators are also needed. This work deals with both problems. We first extend the work of Maffione et al. (2011) for mappings to the 2D Hénon & Heiles (1964) potential, and compare several variational indicators of chaos: the Lyapunov Indicator (LI); the Mean Exponential Growth Factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI) and its generalized version, the Generalized Alignment Index (GALI); the Fast Lyapunov Indicator (FLI) and its variant, the Orthogonal Fast Lyapunov Indicator (OFLI); the Spectral Distance (D) and the Dynamical Spectras of Stretching Numbers (SSNs). We also include in the record the Relative Lyapunov Indicator (RLI), which is not a variational indicator as the others. Then, we test a numerical technique to integrate Ordinary Differential Equations (ODEs) based on the Taylor method implemented by Jorba & Zou (2005) (called taylor), and we compare its performance with other two well-known efficient integrators: the Prince & Dormand (1981) implementation of a Runge-Kutta of order 7-8 (DOPRI8) and a Bulirsch-Stöer implementation. These tests are run under two very different systems from the complexity of their equations point of view: a triaxial galactic potential model and a perturbed 3D quartic oscillator.

35 pages, 16 figures

References in corpus (6)

• Detecting chaos, determining the dimensions of tori and predicting slow diffusion in Fermi--Pasta--Ulam lattices by the Generalized Alignment Index method
• A comparison of different indicators of chaos based on the deviation vectors. Application to symplectic mappings
• Regular and chaotic orbits in barred galaxies - I. Applying the SALI/GALI method to explore their distribution in several models
• Long term stability and dynamical environment of the PSR 1257+12 planetary system
• Are retrograde resonances possible in multi-planet systems?
• Testing a Fast Dynamical Indicator: The MEGNO

Cited by in corpus (14)

• Chaotic behavior of three interacting vortices in a confined Bose-Einstein condensate
• Chaos detection tools: application to a self-consistent triaxial model
• On the relevance of chaos for halo stars in the Solar Neighbourhood
• Orbital dynamics in the post-Newtonian planar circular Sun-Jupiter system
• Orbit classification in the planar circular Pluto-Charon system
• On the relevance of chaos for halo stars in the solar neighbourhood II
• Orbital and escape dynamics in barred galaxies - I. The 2D system
• Finding regions of bounded motion in binary asteroid environment using Lagrangian descriptors
• Comparing the escape dynamics in tidally limited star cluster models
• Polynomial Stochastic Dynamical Indicators
• Propagation of rays in 2D and 3D waveguides: a stability analysis with Lyapunov and Reversibility fast indicators
• Escape dynamics in a binary system of interacting galaxies
• Fugitive stars in active galaxies
• Smart: A program to automatically compute accelerations and variational equations