We report a new theoretical approach to solve adiabatic quantum molecular dynamics halfway between wave function and trajectory-based methods. The evolution of a N-body nuclear wave function moving on a 3N-dimensional Born− Oppenheimer potential-energy hyper-surface is rewritten in terms of single-nuclei wave functions evolving nonunitarily on a 3-dimensional potential-energy surface that depends parametrically on the configuration of an ensemble of generally defined trajectories. The scheme is exact and, together with the use of trajectory-based statistical techniques, can be exploited to circumvent the calculation and storage of many-body quantities (e.g., wave function and potential-energy surface) whose size scales exponentially with the number of nuclear degrees of freedom. As a proof of concept, we present numerical simulations of a 2-dimensional model porphine where switching from concerted to sequential double proton transfer (and back) is induced quantum mechanically. O n the basis of the Born−Huang expansion of the molecular wave function, 1 an exact description of adiabatic molecular dynamics requires the propagation of a nuclear wavepacket on the ground-state Born−Oppenheimer potential-energy surface (gs-BOPES). This propagation scheme is, somehow, computationally doubly prohibitive. Besides the computational burden associated with the propagation of the (many-body) nuclear wave function, the calculation of the gs-BOPES constitutes, per se, a time-independent problem that grows exponentially with the number of electrons and nuclei. In this respect, two main classes of computational methods have emerged depending on whether the knowledge of the gs-BOPES is required in the full configuration space, that is, fullquantum methods, 2,3 or only at certain reduced number of points, namely trajectory-based or direct methods.4,5 While methods for computing the energy of any configuration of nuclei have become quicker and more accurate, full-quantum dynamics calculations still become rapidly unfeasible for large molecules. Alternatively, direct dynamics notably reduce the computational cost of the simulations by avoiding partially, sometimes completely, the calculation of the full gs-BOPES (this can be done, for instance, by the use of reaction-path Hamiltonians 6,7 ). Nuclear quantum effects, however, can be hardly included systematically in this second class of methods. Up to date, only quantum-trajectory methods have the particularity of being able to describe all nuclear quantum effects (just as full-quantum methods) and being on-the-fly simultaneously.8−11 Unfortunately, these methods have serious problems in dealing with the so-called quantum potential, which gathers, by definition, all quantum information on the system. The mathematical structure of the quantum potential depends on the inverse of the quantum probability density, and thus, its manipulation entails serious instability problems.
12−14We report here an exact theoretical approach to solve adiabatic quantum molecular dynamics based o...