“…Extending adiabatic ab initio molecular dynamics to the propagation of classical nuclei in a single non-interacting excited state is straightforward in the framework of wavefunction-based methods such as Hartree-Fock, 78-88 generalized valence bond (GVB), [89][90][91][92][93] multi-configuration self-consistent field (MCSCF), 79 complete active space SCF (CASSCF), 82,94,95 or full configuration interaction (FCI) 96 approaches. However, these methods are computationally rather demanding -even given present-day methods, algorithms, and hardware.…”