Dynamics in the quantum/classical limit based on selective use of the quantum potential

Multidimensional Tunneling Dynamics Employing Quantum-Trajectory Guided Adaptable Gaussian Bases

An efficient basis representation of time-dependent wavefunctions is essential for theoretical studies of high-dimensional molecular systems exhibiting large-amplitude motion. For fully coupled anharmonic systems, the complexity of a general wavefunction scales exponentially with the system size; therefore, for practical reasons, it is desirable to adapt the basis to the time-dependent wavefunction at hand. Often times on this quest for a minimal basis representation, time-dependent Gaussians are employed, in part because of their localization in both configuration and momentum spaces and also because of their direct connection to classical and semiclassical dynamics, guiding the evolution of the basis function parameters. In this work, the quantum-trajectory guided adaptable Gaussian (QTAG) bases method [ J. Chem. Theory Comput. 2020, 16, 18-34] is generalized to include correlated, i.e., non-factorizable, basis functions, and the performance of the QTAG dynamics is assessed on benchmark system/bath tunneling models of up to 20 dimensions. For the popular choice of initial conditions describing tunneling between the reactant/product wells, the minimal "semiclassical" description of the bath modes using essentially a single multidimensional basis function combined with the multi-Gaussian representation of the tunneling mode is shown to capture the dominant features of dynamics in a highly efficient manner.

Quantum Stochastic Trajectories: The Fokker-Planck-Bohm Equation Driven by the Reduced Density Matrix

The quantum molecular trajectory is the deterministic trajectory, arising from the Bohm theory, that describes the instantaneous positions of the nuclei of molecules by assuring the agreement with the predictions of quantum mechanics. Therefore, it provides the suitable framework for representing the geometry and the motions of molecules without neglecting their quantum nature. However, the quantum molecular trajectory is extremely demanding from the computational point of view, and this strongly limits its applications. To overcome such a drawback, we derive a stochastic representation of the quantum molecular trajectory, through projection operator techniques, for the degrees of freedom of an open quantum system. The resulting Fokker-Planck operator is parametrically dependent upon the reduced density matrix of the open system. Because of the pilot role played by the reduced density matrix, this stochastic approach is able to represent accurately the main features of the open system motions both at equilibrium and out of equilibrium with the environment. To verify this procedure, the predictions of the stochastic and deterministic representation are compared for a model system of six interacting harmonic oscillators, where one oscillator is taken as the open quantum system of interest. The undeniable advantage of the stochastic approach is that of providing a simplified and self-contained representation of the dynamics of the open system coordinates. Furthermore, it can be employed to study the out of equilibrium dynamics and the relaxation of quantum molecular motions during photoinduced processes, like photoinduced conformational changes and proton transfers.