Beyond Born-Oppenheimer Time-Dependent Density Functional Theory

We formulate a time-dependent density functional theory for the coupled dynamics of electrons and nuclei that goes beyond the Born-Oppenheimer (BO) approximation. We prove that the time-dependent marginal nuclear probability density $\mathrm{\mid }\chi (\text{bdu}R,t){\mathrm{\mid }}^2$, the conditional electronic density $n_{\text{bdu}R}(\text{br},t)$, and the current density ${\mathit{J}}_{\text{bdu}R}(\text{br},t)$ are sufficient to uniquely determine the full time-evolving electron-nuclear wave function, and thus the dynamics of all observables. Moreover, we propose a time-dependent Kohn-Sham scheme which reproduces the exact conditional electronic density and current density and the exact N-body nuclear density. The remaining task is to look for functional approximations for the Kohn-Sham exchange-correlation scalar and vector potentials. Using a model driven proton transfer system, we numerically demonstrate that the adiabatic extension of a beyond-BO ground state functional captures the dominant nonadiabatic effects in the regime of slow driving.