Regarding the validity of the time-dependent Kohn-Sham approach for electron-nuclear dynamics via trajectory surface hopping.

The implementation of fewest-switches surface-hopping (FSSH) within time-dependent Kohn-Sham (TDKS) theory [Phys. Rev. Lett. 95, 163001 (2005)] has allowed us to study successfully excited state dynamics involving many electronic states in a variety of molecular and nanoscale systems, including chromophore-semiconductor interfaces, semiconductor and metallic quantum dots, carbon nanotubes and graphene nanoribbons, etc. At the same time, a concern has been raised that the KS orbital basis used in the calculation provides only approximate potential energy surfaces [J. Chem. Phys. 125, 014110 (2006)]. While this approximation does exist in our method, we show here that FSSH-TDKS is a viable option for computationally efficient calculations in large systems with straightforward excited state dynamics. We demonstrate that the potential energy surfaces and nonadiabatic transition probabilities obtained within the TDKS and linear response (LR) time-dependent density functional theories (TDDFT) agree semiquantitatively for three different systems, including an organic chromophore ligating a transition metal, a quantum dot, and a small molecule. Further, in the latter case the FSSH-TDKS procedure generates results that are in line with FSSH implemented within LR-TDDFT. The FSSH-TDKS approach is successful for several reasons. First, single-particle KS excitations often give a good representation of LR excitations. In this regard, DFT compares favorably with the Hartree-Fock theory, for which LR excitations are typically combinations of multiple single-particle excitations. Second, the majority of the FSSH-TDKS applications have been performed with large systems involving simple excitations types. Excitation of a single electron in such systems creates a relatively small perturbation to the total electron density summed over all electrons, and it has a small effect on the nuclear dynamics compared, for instance, with thermal nuclear fluctuations. In such cases an additional, classical-path approximation can be made. Third, typical observables measured in time-resolved experiments involve averaging over many initial conditions. Such averaging tends to cancel out random errors that may be encountered in individual simulated trajectories. Finally, if the flow of energy between electronic and nuclear subsystems is insignificant, the ad hoc FSSH procedure is not required, and a straightforward mean-field, Ehrenfest approach is sufficient. Then, the KS representation provides rigorously a convenient and efficient basis for numerically solving the TDDFT equations of motion.