Multidimensional quantum dynamics with trajectories: a novel numerical implementation of Bohmian mechanics

A novel implementation of the de Broglie–Bohm mechanics is presented. The method employs the use of n-dimensional Delaunay tesselation for the purpose of computing the quantum potential term and is fully generalizable for the multidimensional case. We simulate the scattering of a Gaussian wavepacket from an Eckart barrier in two- and three-dimensions and compare our results against the dynamics obtained using a numerically exact propagation scheme.