Optimal conditions for the numerical calculation of the largest Lyapunov exponent for systems of ordinary differential equations

Abstract

A general indicator of the presence of chaos in a dynamical system is the largest Lyapunov exponent. This quantity provides a measure of the mean exponential rate of divergence of nearby orbits. In this paper, we show that the so-called two-particle method introduced by Benettin et al. could lead to spurious estimations of the largest Lyapunov exponent. As a comparator method, the maximum Lyapunov exponent is computed from the solution of the variational equations of the system. We show that the incorrect estimation of the largest Lyapunov exponent is based on the setting of the renormalization time and the initial distance between trajectories. Unlike previously published works, we here present three criteria that could help to determine correctly these parameters so that the maximum Lyapunov exponent is close to the expected value. The results have been tested with four well known dynamical systems: Ueda, Duffing, Rössler and Lorenz.

12 pages, 8 figures. Accepted in the International Journal of Modern Physics C

Cited by in corpus (5)

• Pseudo-Newtonian planar circular restricted 3-body problem
• On the conservation of the Jacobi integral in the post-Newtonian circular restricted three-body problem
• Dynamical analysis of bounded and unbounded orbits in a generalized Hénon-Heiles system
• Signatures of the interplay between chaos and local criticality on the dynamics of scrambling in many-body systems
• Orbital dynamics in realistic galaxy models: NGC 3726, NGC 3877 and NGC 4010