We describe the Reversibility Error Method (REM) and its applications to planetary dynamics. REM is based on the time-reversibility analysis of the phase-space trajectories of conservative Hamiltonian systems. The round-off errors break the time reversibility and the displacement from the initial condition, occurring when we integrate it forward and backward for the same time interval, is related to the dynamical character of the trajectory. If the motion is chaotic, in the sense of non-zero maximal Characteristic Lyapunov Exponent (mLCE), then REM increases exponentially with time, as exp \lambda t, while when the motion is regular (quasi-periodic) then REM increases as a power law in time, as t^{\alpha}, where {\alpha} and {\lambda} are real coefficients. We compare the REM with a variant of mLCE, the Mean Exponential Growth factor of Nearby Orbits (MEGNO). The test set includes the restricted three body problem and five resonant planetary systems: HD 37124, Kepler-60, Kepler-36, Kepler-29 and Kepler-26. We found a very good agreement between the outcomes of these algorithms. Moreover, the numerical implementation of REM is astonishing simple, and is based on solid theoretical background. The REM requires only a symplectic and time-reversible (symmetric) integrator of the equations of motion. This method is also CPU efficient. It may be particularly useful for the dynamical analysis of multiple planetary systems in the KEPLER sample, characterized by low-eccentricity orbits and relatively weak mutual interactions. As an interesting side-result, we found a possible stable chaos occurrence in the Kepler-29 planetary system.

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.