An efficient algorithm based on splitting for the time integration of the Schrödinger equation

From MaRDI portal
Publication:2374876

DOI10.1016/J.JCP.2015.09.047zbMATH Open1349.65393arXiv1502.06401OpenAlexW1882722114MaRDI QIDQ2374876FDOQ2374876


Authors: Sergio Blanes, Fernando Casas, Ander Murua Edit this on Wikidata


Publication date: 5 December 2016

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: We present a practical algorithm based on symplectic splitting methods to integrate numerically in time the Schr"odinger equation. When discretized in space, the Schr"odinger equation can be recast as a classical Hamiltonian system corresponding to a generalized high-dimensional separable harmonic oscillator. The particular structure of this system combined with previously obtained stability and error analyses allows us to construct a set of highly efficient symplectic integrators with sharp error bounds and optimized for different tolerances and time integration intervals. They can be considered, in this setting, as polynomial approximations to the matrix exponential in a similar way as methods based on Chebyshev and Taylor polynomials. The theoretical analysis, supported by numerical experiments, indicates that the new methods are more efficient than schemes based on Chebyshev polynomials for all tolerances and time intervals. The algorithm we present incorporates the new splitting methods and automatically selects the most efficient scheme given a tolerance, a time integration interval and an estimate on the spectral radius of the Hamiltonian.


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




Recommendations




Cites Work


Cited In (20)

Uses Software





This page was built for publication: An efficient algorithm based on splitting for the time integration of the Schrödinger equation

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