https://doi.org/10.1140/epjst/e2016-02657-0
Regular Article
On the symplectic integration of the discrete nonlinear Schrödinger equation with disorder
1 Lohrmann Observatory, Technical University Dresden, 01062 Dresden, Germany
2 Department of Mathematics and Applied Mathematics, University of Cape Town, Rondebosch 7701, South Africa
a e-mail: enrico.gerlach@tu-dresden.de
Received: 28 March 2016
Revised: 26 July 2016
Published online: 30 September 2016
We present several methods, which utilize symplectic integration techniques based on two and three part operator splitting, for numerically solving the equations of motion of the disordered, discrete nonlinear Schrödinger (DDNLS) equation, and compare their efficiency. Our results suggest that the most suitable methods for the very long time integration of this one-dimensional Hamiltonian lattice model with many degrees of freedom (of the order of a few hundreds) are the ones based on three part splits of the system's Hamiltonian. Two part split techniques can be preferred for relatively small lattices having up to N ≈ 70 sites. An advantage of the latter methods is the better conservation of the system's second integral, i.e. the wave packet's norm.
© EDP Sciences, Springer-Verlag, 2016
