Efficient Splitting Methods Based on Modified Potentials: Numerical Integration of Linear Parabolic Problems and Imaginary Time Propagation of the Schrodinger Equation

Abstract

We present a new family of fourth-order splitting methods with positive coefficients especially tailored for the time integration of linear parabolic problems and, in particular, for the time dependent Schrödinger equation, both in real and imaginary time. They are based on the use of a double commutator and a modified processor, and are more efficient than other widely used schemes found in the literature. Moreover, for certain potentials, they achieve order six. Several examples in one, two and three dimensions clearly illustrate the computational advantages of the new schemes.