Lyapunov and reversibility errors for Hamiltonian flows. (July 2018)

• Record Type: Journal Article
• Authors: Panichi, F.; Turchetti, G.
• Is Part Of: Chaos, solitons and fractals. Volume 112(2018)
• Page Start: 83
• Page End: 91
• Publication Date: 2018-07
• Subjects: Hamiltonian systems -- Stability analysis -- Reversibility -- Resonances; Chaotic behavior in systems -- Periodicals; Solitons -- Periodicals; Fractals -- Periodicals; Chaotic behavior in systems; Fractals; Solitons; Periodicals; 003.7
• DOI: 10.1016/j.chaos.2018.03.019
• Languages: English
• ISSNs: 0960-0779

Abstract

Abstract: We discuss the stability of a Hamiltonian system by comparing the standard Lyapunov error (LE) with the forward error (FE) due to a small random perturbation. We introduce also the reversibility error (RE) where the evolution is computed forward up to time t and backwards to t = 0 in presence of noise. This procedure has been investigated in the case of symplectic maps, but it turns out that the results are simpler in the case of a noisy flow, in the limit of zero noise amplitude. Indeed the stochastic processes defined by the displacement of the noisy orbit at time t for FE, or at time 0 for RE after the evolution up to time t, satisfy linear Langevin equations, are Gaussian processes, and the errors are just their root mean square deviations. All the errors are expressed in terms of the fundamental matrix L ( t ) of the tangent flow and can be evaluated numerically using a symplectic integrator. Letting eL ( t ) be the Lyapunov error and eR ( t ) be the reversibility error a very simple relation holds e R 2 ( t ) = ∫ 0 t e L 2 ( s ) d s . The integral relation is quite natural since the local errors due to a random perturbations accumulate during the evolution whereas for the Lyapunov case the error is introduced only at time zero and propagated. The plot of errors for initial conditions in a Poincaré section reflects the phase portrait, whereas in the action plane it allows to single out the resonance strips. We have applied the method to a 3D Hamiltonian model H = H 0 ( J ) + λ V ( Θ ), where analytic estimates can be obtained for the single resonances from perturbation theory. This allows to inspect the double resonance structure where the single resonance strips intersect. We have also considered the Hénon–Heiles Hamiltonian to show numerically the equivalence of the errors apart from a shift of 1/2 in the power law exponent in the case of regular orbits. The reversibility error method (REM), previously introduced as the error due to round off in the symplectic integration, appears to be comparable with RE also for the models considered here.