Intramolecular vibrational redistribution in the non-radiative excited state decay of uracil in the gas phase: an ab initio molecular dynamics study

We report a study of intramolecular vibrational distribution (IVR) occurring in the electronic ground state of uracil (S-0) in the gas phase, following photoexcitation in the lowest energy bright excited state (S pi) and decay through the ethylene-like S pi/S-0 Conical Intersection (CI-0 pi). To this aim we have performed 20 independent ab initio molecular dynamics simulations starting from CI-0 pi (ten of them with 1 eV kinetic energy randomly distributed over the different molecular degrees of freedom) and 10 starting from the ground state minimum (Franck-Condon, FC, point), with the excess kinetic energy equal to the energy gap between CI-0 pi and the FC point. The simulations, exploiting PBE0/6-31G(d) calculations, were performed over an overall period of 10 ps. A thorough statistical analysis of the variation of the geometrical parameters of uracil during the simulation time and of the distribution of the kinetic energy among the different vibrational degrees of freedom provides a consistent picture of the IVR process. In the first 0-200 fs the structural dynamics involve mainly the recovery of the average planarity. In the 200-600 fs time range, a substantial activation of CO and NH degrees of freedom is observed. After 500-600 fs most of the geometrical parameters reach average values similar to those found after 10 ps, though the system cannot be considered to be in equilibrium yet.