Magnus-factorized method for numerical solving the time-dependent Schrödinger equation

Abstract The method of constructing the high-order stable operator-difference schemes for solving the time-dependent Schrodinger equation are proposed as a generalization of the Crank–Nicolson scheme. The schemes are constructed in terms of the Magnus expansion for the evolution operator. The case of the non-homogeneous equation is also considered.