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

The method of constructing the high-order stable operator-difference schemes for solving the time-dependent Schrödinger 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.