Nonadiabatic quantum dynamics without potential energy surfaces

We present an ab initio algorithm for quantum dynamics simulations that reformulates the traditional “curse of dimensionality” that plagues all state-of-the-art techniques for solving the time-dependent Schrödinger equation. Using a stochastic wave-function ansatz that is based on a set of interacting single-particle conditional wave functions, we show that the difficulty of the problem becomes dominated by the number of trajectories needed to describe the process, rather than simply the number of degrees of freedom involved. This highly parallelizable technique achieves quantitative accuracy for situations in which mean-field theory drastically fails to capture qualitative aspects of the dynamics, such as quantum decoherence or the reduced nuclear probability density, using orders of magnitude fewer trajectories than a mean-field simulation. We illustrate the performance of this method for two fundamental nonequilibrium processes: a photoexcited proton-coupled electron transfer problem, and nonequilibrium dynamics in a cavity bound electron-photon system in the ultrastrong-coupling regime