Computational efficiency of numerical integration methods for the tangent dynamics of many-body Hamiltonian systems in one and two spatial dimensions

Abstract

We investigate the computational performance of various numerical methods for the integration of the equations of motion and the variational equations for some typical classical many-body models of condensed matter physics: the Fermi-Pasta-Ulam-Tsingou (FPUT) chain and the one- and two-dimensional disordered, discrete nonlinear Schrödinger equations (DDNLS). In our analysis we consider methods based on Taylor series expansion, Runge-Kutta discretization and symplectic transformations. The latter have the ability to exactly preserve the symplectic structure of Hamiltonian systems, which results in keeping bounded the error of the system's computed total energy. We perform extensive numerical simulations for several initial conditions of the studied models and compare the numerical efficiency of the used integrators by testing their ability to accurately reproduce characteristics of the systems' dynamics and quantify their chaoticity through the computation of the maximum Lyapunov exponent. We also report the expressions of the implemented symplectic schemes and provide the explicit forms of the used differential operators. Among the tested numerical schemes the symplectic integrators $ABA864$ and $SRK{N}_{14}^a$ exhibit the best performance, respectively for moderate and high accuracy levels in the case of the FPUT chain, while for the DDNLS models $s9\mathcal{A}\mathcal{BC}6$ and $s11\mathcal{A}\mathcal{BC}6$ (moderate accuracy), along with $s17\mathcal{A}\mathcal{BC}8$ and $s19\mathcal{A}\mathcal{BC}8$ (high accuracy) proved to be the most efficient schemes.

arXiv commentAccepted for publication in Mathematics in Engineering

References in corpus (32)

• Anderson localization of a non-interacting Bose-Einstein condensate
• Direct observation of Anderson localization of matter-waves in a controlled disorder
• Anderson localization and nonlinearity in one dimensional disordered photonic lattices
• Destruction of Anderson localization by a weak nonlinearity
• Observation of breathers in Josephson ladders
• The Lyapunov Characteristic Exponents and their computation
• Universal spreading of wavepackets in disordered nonlinear systems
• High order symplectic integrators for perturbed Hamiltonian systems
• Effect of optical disorder and single defects on the expansion of a Bose-Einstein condensate in a one-dimensional waveguide
• Statistical mechanics of a nonlinear discrete system
• Delocalization of wave packets in disordered nonlinear chains
• q-Breathers and the Fermi-Pasta-Ulam Problem
• Intermittent many-body dynamics at equilibrium
• Delocalization induced by nonlinearity in systems with disorder
• New families of symplectic splitting methods for numerical integration in dynamical astronomy
• On the universality of anomalous one-dimensional heat conductivity
• Numerical integration of variational equations
• Spreading of wave packets in disordered systems with tunable nonlinearity
• High precision Symplectic Integrators for the Solar System
• Dynamical Thermalization of Disordered Nonlinear Lattices
• Weakly nonergodic dynamics in the Gross--Pitaevskii lattice
• Nonequilibrium chaos of disordered nonlinear waves
• Nonlinear waves in disordered chains: probing the limits of chaos and spreading
• Chaotic Dynamics of N-degree of Freedom Hamiltonian Systems
• The Fermi-Pasta-Ulam problem: periodic orbits, normal forms and resonance overlap criteria
• Dynamical glass in weakly non-integrable Klein-Gordon chains
• Energy transport in one-dimensional disordered granular solids
• Heterogeneity and chaos in the Peyrard-Bishop-Dauxois DNA model
• High order three part split symplectic integrators: Efficient techniques for the long time simulation of the disordered discrete nonlinear Schroedinger equation
• Computational efficiency of symplectic integration schemes: Application to multidimensional disordered Klein-Gordon lattices
• Dedicated symplectic integrators for rotation motions
• High order explicit symplectic integrators for the Discrete Non Linear Schrödinger equation