Convergence analysis of high-order commutator-free quasi-Magnus exponential integrators for nonautonomous linear Schrodinger equations

Abstract

[EN] This work is devoted to the derivation of a convergence result for high-order commutator-free quasi-Magnus (CFQM) exponential integrators applied to nonautonomous linear Schrodinger equations; a detailed stability and local error analysis is provided for the relevant special case where the Hamilton operator comprises the Laplacian and a regular space-time-dependent potential. In the context of nonautonomous linear ordinary differential equations, CFQM exponential integrators are composed of exponentials involving linear combinations of certain values of the associated time-dependent matrix; this approach extends to nonautonomous linear evolution equations given by unbounded operators. An inherent advantage of CFQM exponential integrators over other time integration methods such as Runge-Kutta methods or Magnus integrators is that structural properties of the underlying operator family are well preserved; this characteristic is confirmed by a theoretical analysis ensuring unconditional stability in the underlying Hilbert space and the full order of convergence under low regularity requirements on the initial state. Due to the fact that convenient tools for products of matrix exponentials such as the Baker-Campbell-Hausdorff formula involve infinite series and thus cannot be applied in connection with unbounded operators, a certain complexity in the investigation of higher-order CFQM exponential integrators for Schrodinger equations is related to an appropriate treatment of compositions of evolution operators; an effective concept for the derivation of a local error expansion relies on suitable linearisations of the evolution equations for the exact and numerical solutions, representations by the variation-ofconstants formula and Taylor series expansions of parts of the integrands, where the arising iterated commutators determine the regularity requirements on the problem data.