Efficient finite difference solutions to the time-dependent Schrödinger equation (Q676367)

For solving the initial value and time-dependent problem for the Schrödinger equation an expansion of the evolution operator in terms of unitary anti-Hermitian operators is used. The one-dimensional Laplacian is represented by the standard finite difference approximation of second order and the corresponding exponential is calculated exactly by means of Bessel functions. The proposed algorithms are generalized to two spatial dimensions.