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

DOI10.3934/MINE.2019.3.447MaRDI QIDQ2305083zbMATH OpenOpenAlexWikidataFDO

Authors Carlo Danieli, Bertin Many Manda, Thudiyangal Mithun, Ch. Skokos

Publication date 10 March 2020

Published in Mathematics in Engineering (Search for Journal in Brave)

Copyright license Creative Commons Attribution 4.0 International

Full work available at URL https://arxiv.org/abs/1812.01870

zbMATH Keywords

optimizationLyapunov exponentordinary differential equationscomputational efficiencysymplectic integratorsvariational equationsclassical many-body systems

Mathematics Subject Classification ID

Numerical optimization and variational techniques (65K10) NLS equations (nonlinear Schrödinger equations) (35Q55) Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations (65L06) Numerical methods for Hamiltonian systems including symplectic integrators (65P10) Many-body theory; quantum Hall effect (81V70) Disordered systems (random Ising models, random Schrödinger operators, etc.) in equilibrium statistical mechanics (82B44) Characteristic and Lyapunov exponents of ordinary differential equations (34D08) Computational molecular dynamics in statistical mechanics (82M37)

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"odinger 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 $s9mathcalABC6$ and $s11mathcalABC6$ (moderate accuracy), along with $s17mathcalABC8$ and $s19mathcalABC8$ (high accuracy) proved to be the most efficient schemes.

Recommendations

• Numerical integration of variational equations for Hamiltonian systems with long range interactions
• Efficient integration of the variational equations of multidimensional Hamiltonian systems: application to the Fermi-Pasta-Ulam lattice
• Publication:4941673
• Methods for the numerical integration of Hamiltonian systems
• scientific article; zbMATH DE number 177016
• A high-precision numerical method for the dissipative dynamics of the Hamiltonian systems
• Effective computation of the dynamics around a two-dimensional torus of a Hamiltonian system
• Efficient numerical approximation of the effective Hamiltonian
• Closed-form modified Hamiltonians for integrable numerical integration schemes
• Energy-preserving numerical schemes of high accuracy for one-dimensional Hamiltonian systems

Cites work

• scientific article; zbMATH DE number 3856471 (Why is no real title available?)
• scientific article; zbMATH DE number 5353839 (Why is no real title available?)
• scientific article; zbMATH DE number 1031368 (Why is no real title available?)
• scientific article; zbMATH DE number 1745051 (Why is no real title available?)
• scientific article; zbMATH DE number 3241199 (Why is no real title available?)
• scientific article; zbMATH DE number 3253340 (Why is no real title available?)
• A concise introduction to geometric numerical integration
• Algorithm 924, TIDES, a Taylor series integrator for differential equations
• Analyzing chaos in higher order disordered quartic-sextic Klein-Gordon lattices using \(q\)-statistics
• CHAOTIC DYNAMICS OF N-DEGREE OF FREEDOM HAMILTONIAN SYSTEMS
• Canonical Runge-Kutta methods
• Chaotic behavior in nonlinear Hamiltonian systems and equilibrium statistical mechanics.
• Chopping time of the FPU \({\alpha }\)-model
• Comparing the efficiency of numerical techniques for the integration of variational equations
• Complex statistics and diffusion in nonlinear disordered particle chains
• Composition constants for raising the orders of unconventional schemes for ordinary differential equations
• Composition methods in the presence of small parameters
• Dedicated symplectic integrators for rotation motions
• Derivation of symmetric composition constants for symmetric integrators
• Efficient integration of the variational equations of multidimensional Hamiltonian systems: application to the Fermi-Pasta-Ulam lattice
• Fourth-order symplectic integration
• Geometric integration for particle accelerators
• Geometric integrators for ODEs
• High order three part split symplectic integrators: Efficient techniques for the long time simulation of the disordered discrete nonlinear Schrödinger equation
• High-order symplectic integrators for perturbed Hamiltonian systems
• Interaction of `solitons' in a collisionless plasma and the recurrence of initial states
• Lyapunov characteristic exponents for smooth dynamical systems and for Hamiltonian systems; a method for computing all of them. I: Theory
• New families of symplectic splitting methods for numerical integration in dynamical astronomy
• On the numerical integration of FPU-like systems
• Performance of the Taylor series method for ODEs/DAEs
• Practical symplectic partitioned Runge-Kutta and Runge-Kutta-Nyström methods
• Recent progress in the theory and application of symplectic integrators
• Runge-Kutta schemes for Hamiltonian systems
• Simulating Hamiltonian Dynamics
• Solving Ordinary Differential Equations I
• Splitting and composition methods in the numerical integration of differential equations
• Splitting methods
• The Fermi-Pasta-Ulam problem: periodic orbits, normal forms and resonance overlap criteria
• Wave interactions in localizing media -- a coin with many faces

Cited in

(11)

• Dynamical chaos in the integrable Toda chain induced by time discretization
• Numerical integration of variational equations for Hamiltonian systems with long range interactions
• Effect of discrete breathers on the specific heat of a nonlinear chain
• Exact numerical methods for a many-body Wannier-Stark system
• Identifying localized and spreading chaos in nonlinear disordered lattices by the generalized alignment index (GALI) method
• Efficient integration of the variational equations of multidimensional Hamiltonian systems: application to the Fermi-Pasta-Ulam lattice
• Hamiltonian lattice dynamics
• Frequency Map Analysis of Spatiotemporal Chaos in the Nonlinear Disordered Klein–Gordon Lattice
• Wave-packet spreading in disordered soft architected structures
• Effective computation of the dynamics around a two-dimensional torus of a Hamiltonian system
• Fast numerical computation of quasi-periodic equilibrium states in 1D statistical mechanics, including twist maps

Describes a project that uses

Uses Software

• TIDES
• Matlab

This page was built for publication: Computational efficiency of numerical integration methods for the tangent dynamics of many-body Hamiltonian systems in one and two spatial dimensions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2305083)