Fidelity and reversibility in the restricted three body problem

Abstract

We use the Reversibility Error Method and the Fidelity to analyze the global effects of a small perturbation in a non-integrable system. Both methods have already been proposed and used in the literature but the aim of this paper is to compare them in a physically significant example adding some considerations on the equivalence, observed in this case, between round-off and random perturbations.As a paradigmatic example we adopt the restricted planar circular three body problem. The cumulative effect of random perturbations or round-off leads to a divergence of the perturbed orbit from the reference one. Rather than computing the distance of the perturbed orbit from the reference one, after a given number n of iterations, a procedure we name the Forward Error Method (FEM), we measure the distance of the reversed orbit (n periods forward and backward) from the initial point. This approach, that we name Reversibility Error Method (REM), does not require the computation of the unperturbed map. The loss of memory of the perturbed map is quantified by the Fidelity decay rate whose computation requires a statistical average over an invariant region. Two distinct definitions of Fidelity are given. The asymptotic equivalence of REM and FEM is analytically proved for linear symplectic maps with random perturbations. For a given map, the REM plot provides a picture of the dynamic stability regions in the phase space, very easy to obtain for any kind of perturbation and very simple to implement numerically. The REM and FEM for linear symplectic maps are proved to be asymptotically equivalent. The global error growth follows a power law in the regions of integrable (or quasi integrable) motion and an exponential law in the regions of chaotic motion. We prove that the power law exponent is 3/2 for a generic anisochronous system, but drops down to 1/2 if the system is isochronous. Correspondingly the Fidelity F (t) exhibits an exponential decay and − ln F (t) grows just as the square of the FEM or REM error. The Reversibility Error and Fidelity can be used for a quantitative analysis of dynamical systems and are suited to investigate the transition regions from chaotic to regular motion even for Hamiltonian systems with many degrees of freedom such as the N -body problem.

The reversibility error method (REM): a new, dynamical fast indicator for planetary dynamics

“…To apply Eq. 2 numerically, we must guarantee that the map is invertible (Faranda et al 2012;Panichi et al 2016). For a numerical integrator affected by a round-off error of amplitude γ, we change Eq.…”

“…Unlike regular orbits, an ergodic motion is expected to result in large displacements of the initial condition x x x 0 after the forward and backward integration. Since SI are equivalent to symplectic maps, it makes it possible to determine and rigorously prove analytic properties of a numerical approach based on this idea developed in a series of papers (Turchetti et al 2010a,b;Faranda et al 2012;Panichi et al 2016).…”

The architecture and formation of the Kepler-30 planetary system

“…Here, we computed MEGNO with the symplectic, fourth-order integrator scheme SABA 4 (Laskar & Robutel 2001) and the symplectic tangent map (Mikkola & Innanen 1999;Goździewski et al 2008). The REM fast indicator (Faranda et al 2012;Panichi et al 2016) is the maximum Lyapunov exponent-like, CPU efficient fast indicator optimised to compute high-resolution dynamical maps for low-eccentric KE-PLER systems. The REM is based on the loose of time-reversibility propriety of chaotic orbits in conservative Hamiltonian systems.…”

Fast Indicators for Orbital Stability: A Survey on Lyapunov and Reversibility Errors

“…In the framework of the variational methods, we have proposed two indicators [10][11][12] the Lyapunov error (LE) and the reversibility error (RE) introducing also the modified reversibility error method (REM). The LE is due to a small displacement of the initial condition, the RE is due to an additive noise, and REM is due to roundoff.…”