Prototropically Controlled Dynamics of Cytosine Photodecay

The effect of the existence of several prototropic tautomers of cytosine on its UV/vis spectra and the excited state decay dynamics is studied by spectral and nonadiabatic molecular dynamics (NAMD) simulations in connection with the mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. Simulated UV/vis spectra provide a strong indication that the H3N keto-amino cytosine tautomer (the least anticipated species) may be present under experimental conditions. The NAMD simulations yield a wide range of excited state decay constants for various tautomers of cytosine, ranging from ∼1.3 ps for the biologically relevant H1N keto-amino tautomer to ∼0.1 ps for the keto-imino tautomer. The slowness of the H1N decay dynamics follows from the presence of a barrier on the excited state energy surface separating the Franck-Condon structure from the major decay funnel, the conical intersection seam. It is suggested that the experimentally observed photodecay dynamics may result from a combination of the decay processes of various tautomers (H3N in particular) present simultaneously under the experimental conditions.