On Magnus integrators for time-dependent Schrodinger equations

Numerical methods based on the Magnus expansion are an efficient class of integrators for Schrodinger 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 that 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.