An efficient algorithm to compute the exponential of skew-Hermitian matrices for the time integration of the Schr\"odinger equation

DOI10.1016/J.MATCOM.2021.12.002arXiv2103.10132WikidataQ114149929 ScholiaQ114149929MaRDI QIDQ6363142

Fernando Casas, Philipp Bader, Sergio Blanes, Muaz Seydaoğlu

Publication date: 18 March 2021

Abstract: We present a practical algorithm to approximate the exponential of skew-Hermitian matrices up to round-off error based on an efficient computation of Chebyshev polynomials of matrices and the corresponding error analysis. It is based on Chebyshev polynomials of degrees 2, 4, 8, 12 and 18 which are computed with only 1, 2, 3, 4 and 5 matrix-matrix products, respectively. For problems of the form $exp(-iA)$, with $A$ a real and symmetric matrix, an improved version is presented that computes the sine and cosine of $A$ with a reduced computational cost. The theoretical analysis, supported by numerical experiments, indicates that the new methods are more efficient than schemes based on rational Pad'e approximants and Taylor polynomials for all tolerances and time interval lengths. The new procedure is particularly recommended to be used in conjunction with exponential integrators for the numerical time integration of the Schr"odinger equation.

Mathematics Subject Classification ID

Numerical analysis (65-XX) Quantum theory (81-XX)

This page was built for publication: An efficient algorithm to compute the exponential of skew-Hermitian matrices for the time integration of the Schr\"odinger equation