This article reviews several numerical methods for the time-dependent Schrödinger Equation (TDSE). We consider both the most commonly used approach—short-time propagation, which solves the TDSE by assuming that the Hamiltonian is time-independent over sufficiently small (time) intervals—as well as a number of higher-order alternatives. Our goal is to dispel the notion that the latter are too computationally demanding for practical use. To that end, we cover methods whose numerical building blocks are shared by short-time propagators or can be handled by standard libraries. Moreover, we make the case that these methods are best positioned to take advantage of parallel computing environments. One of the alternatives considered is a “double DVR” solver, which applies an expansion in a product basis of functions in space and time to obtain a solution (over all space and at multiple time points simultaneously) with a single linear system solve. To our knowledge, and despite its simplicity, this approach has not previously been applied to the TDSE.
Schneider et al. (Mon,) studied this question.