On Magnus integrators for time-dependent Schrödinger equations

Numerical methods based on the Magnus expansion are an efficient class of integrators for Schrödinger equations with time-dependent Hamiltonian. Though their derivation assumes an unreasonably small time step size as would be required for a standard explicit integrator, the methods perform well even for much larger step sizes. This favorable behavior is explained, and optimal-order error bounds are derived which require no or only mild restrictions of the step size. In contrast to standard integrators, the error does not depend on higher time derivatives of the solution, which is in general highly oscillatory.