2013/06/04 by E. Gerlach, Gerlach, Enrico, Siegfried Eggl +7
Mathematics · Physics and Astronomy · #Chaotic Dynamics (nlin.CD) #Computational Physics (physics.comp-ph) #FOS: Mathematics #FOS: Physical sciences #Nonlinear Photonic Systems #Numerical Analysis (math.NA) #Numerical methods for differential equations #Quantum chaos and dynamical systems
paper · pdf · doi:10.48550/arxiv.1306.0627
openalex publication_date 2013/06/04 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
Symplectic integration methods based on operator splitting are well established in many branches of science. For Hamiltonian systems which split in more than two parts, symplectic methods of higher order have been studied in detail only for a few special cases. In this work, we present and compare different ways to construct high order symplectic schemes for general Hamiltonian systems that can be split in three integrable parts. We use these techniques to numerically solve the equations of motion for a simple toy model, as well as the disordered discrete nonlinear Schrödinger equation. We thereby compare the efficiency of symplectic and non-symplectic integration methods. Our results show that the new symplectic schemes are superior to the other tested methods, with respect to both long term energy conservation and computational time requirements.