Abstract
We describe two different implementations of time-dependent density functional theory (TDDFT) for use in excited state molecular dynamics simulations. One is based on the linear response formulation (LR-TDDFT), whereas the other uses a time propagation scheme for the electronic wave functions (P-TDDFT). Photo-induced cis–trans isomerization of C=C, C=N and N=N double bonds is investigated in three model compounds, namely the 2,4-pentadiene-1-iminium cation (PSB), formaldimine and diimide. For formaldimine and diimide, the results obtained with both schemes are in agreement with experimental data and previously reported theoretical results. Molecular dynamics simulations yield new insights into the relaxation pathways in the excited state. For PSB, which is a model system for the retinal protonated Schiff base involved in the visual process, the forces computed from the LR-TDDFT S1 surface lead to an increased bond length alternation and, consequently, to single bond rotation. On the contrary, P-TDDFT dynamics lead to a decreased bond length alternation, in agreement with CASPT2 and restricted open-shell Kohn–Sham (ROKS) calculations.
Acknowledgements
We would like to thank J. Hutter for helpful discussions and implementation of the LR-TDDFT approach and M. Sprik for help with the P-TDDFT scheme. F. Buda is acknowledged for providing us with a preprint of his work on excited states. We thank M. Coutinho and L. Guidoni for discussions.