Time-dependent quantum wave-packet description of the 1pi sigma* photochemistry of phenol

The photoinduced hydrogen elimination reaction in phenol via the conical intersections of the dissociative 1pi sigma* state with the 1pi pi* state and the electronic ground state has been investigated by time-dependent quantum wave-packet calculations. A model including three intersecting electronic potential-energy surfaces (S0, 1pi sigma*, and 1pi pi*) and two nuclear degrees of freedom (OH stretching and OH torsion) has been constructed on the basis of accurate ab initio multireference electronic-structure data. The electronic population transfer processes at the conical intersections, the branching ratio between the two dissociation channels, and their dependence on the initial vibrational levels have been investigated by photoexciting phenol from different vibrational levels of its ground electronic state. The nonadiabatic transitions between the excited states and the ground state occur on a time scale of a few tens of femtoseconds if the 1pi pi*-1pi sigma* conical intersection is directly accessible, which requires the excitation of at least one quantum of the OH stretching mode in the 1pi pi* state. It is shown that the node structure, which is imposed on the nuclear wave packet by the initial preparation as well as by the transition through the first conical intersection (1pi pi*-1pi sigma*), has a profound effect on the nonadiabatic dynamics at the second conical intersection (1pi sigma*-S0). These findings suggest that laser control of the photodissociation of phenol via IR mode-specific excitation of vibrational levels in the electronic ground state should be possible.     

Ultrafast decay of electronically excited singlet cytosine via a pi,pi* to n(O),pi* state switch

Singlet fluorescence lifetimes of adenosine, cytidine, guanosine, and thymidine, determined by femtosecond pump-probe spectroscopy (Pecourt, J.-M. L.; Peon, J.; Kohler, B. J. Am. Chem. Soc. 2000, 122, 9348. Pecourt, J.-M. L.; Peon, J.; Kohler, B. J. Am. Chem. Soc. 2001, 123, 10370), show that the excited states produced by 263 nm light in these nucleosides decay in the subpicosecond range (290-720 fs). Ultrafast radiationless decay to the ground state greatly reduces the probability of photochemical damage. In this work we present a theoretical study of isolated cytosine, the chromophore of cytidine. The experimental lifetime of 720 fs indicates that there must be an ultrafast decay channel for this species. We have documented the possible decay channels and approximate energetics, using a valence-bond derived analysis to rationalize the structural details of the paths. The mechanism favored by our calculations and the experimental data involves (1) a two-mode decay coordinate composed of initial bond length inversion followed by internal vibrational energy redistribution (IVR) to populate a carbon pyramidalization mode, (2) a state switch between the pi,pi* and nO,pi* (excitation from oxygen lone pair) excited states, and (3) decay to the ground state through a conical intersention. A second decay path through the nN,pi* state (excitation from the nitrogen lone pair), with a higher barrier, involves out-of-plane bending of the amino substituent.     

Nonradiative decay mechanisms of the biologically relevant tautomer of guanine

We present the excited-state potential energy profiles of the biologically relevant 9H-keto-amino tautomer of guanine with respect to the radiationless decay via the out-of-plane deformation of the six-membered ring as well as the dissociation of NH bonds. The CASPT2//CASSCF method is employed for the reaction-path calculations. The reaction path for the out-of-plane deformation in the (1)pi pi* state leads in a barrierless way to a conical intersection with the electronic ground state. For the NH dissociation via the (1)pi sigma* state, the 9H-keto-amino tautomer is shown to have lower energy barriers than the 7H tautomers which we have studied recently. These two radiationless decay mechanisms explain the unexpected missing of the biologically relevant form in the resonant two-photon ionization spectrum of guanine in a supersonic jet. It is suggested that these ultrafast deactivation processes provide the biologically relevant tautomer of guanine with a high degree of photostability.     

Theoretical study on the singlet excited state of pterin and its deactivation pathway

The excited-state properties and related photophysical processes of the acidic and basic forms of pterin have been investigated by the density functional theory and ab initio methodologies. The solvent effects on the low-lying states have been estimated by the polarized continuum model and combined QM/MM calculations. Calculations reveal that the observed two strong absorptions arise from the strong pi --> pi* transitions to 1(pipi*L(a)) and 1(pipi*L(b)) in the acidic and basic forms of pterin. The first 1(pipi*L(a)) excited state is exclusively responsible for the experimental emission band. The vertical 1(n(N)pi*) state with a small oscillator strength, slightly higher in energy than the 1(pipi*L(a)) state, is less accessible by the direct electronic transition. The 1(n(N)pi*) state may be involved in the photophysical process of the excited pterin via the 1(pipi*L(a)/n(N)pi*) conical intersection. The radiationless decay of the excited PT to the ground state experiences a barrier of 13.8 kcal/mol for the acidic form to reach the (S(1)/S(0)) conical intersection. Such internal conversion can be enhanced with the increase in excitation energy, which will reduce the fluorescence intensity as observed experimentally.     

Ab initio studies on the photophysics of guanine tautomers: out-of-plane deformation and NH dissociation pathways to conical intersections

The radiationless decay mechanisms of the S1 excited states of the 7H-keto-amino, 7H-enol-amino, and 7H-keto-imino tautomers of guanine have been investigated with the CASPT2//CASSCF method. Out-of-plane deformation of the six-membered ring or the imino group as well as dissociation of NH bonds have been considered as photochemical pathways leading to conical intersections with the electronic ground state. It has been found that all three tautomers can reach S0-S1 conical intersections by out-of-plane deformation. However, only in the 7H-keto-amino tautomer the reaction path leading to the conical intersection is barrierless. This tautomer also has the lowest energy barrier for hydrogen detachment via the (1)pi sigma* state, whose potential energy surface intersects that of the (1)pi pi* state as well as that of the ground state. The other tautomers of guanine exhibit substantial energy barriers on their S1 potential energy surfaces with respect to both reaction mechanisms. These findings suggest that the 7H-keto-amino tautomer exhibits the shortest excited-state lifetime of the three tautomers due to particularly fast nonradiative deactivation processes through S0-S1 conical intersections. The computational results explain the remarkable observation that the energetically most stable 7H-keto-amino tautomer is missing in the resonant two-photon ionization spectrum of guanine in a supersonic jet. The results also explain that the energetically less stable 7H-enol-amino and 7H-keto-imino tautomers have longer excited-state lifetimes and are thus detectable by resonant two-photon ionization.     

CASCSF study on the photochemical transposition reactions of pyrazines

Several reaction pathways for the photochemical transformations of methyl-substituted pyrazine in its first excited state 1(pi --> pi*) have been determined using the CASSCF (six-orbital/six-electron active space) and MP2-CAS methods with the 6-311G(d) basis set. Our model investigations suggest that conical intersections play a crucial role in the photoisomerization of pyrazines. Moreover, the present theoretical findings indicate that all of the photoisomerizations of pyrazines adopt the same reaction path as follows: pyrazine --> Franck-Condon region --> conical intersection --> pyrimidine. That is, although an excited-state pyrazine molecule can initiate a phototransposition process easily, this process can be completed on the ground-state potential energy surface after passage through a conical intersection where a fast, radiationless decay is possible. The existence of these nonadiabatic reaction pathways is consistent with the available experimental observations of the photochemistry and photophysics of pyrazine and its methyl derivatives. In the present work, we propose a simple p-pi orbital topology model, which can be used as a diagnostic tool to predict the location of the conical intersections, as well as the geometries of the phototransposition products of various heterocycles.     

Solvent effect on the singlet excited-state dynamics of 5-fluorouracil in acetonitrile as compared with water

The excited-state dynamics of 5-fluorouracil in acetonitrile has been investigated by femtosecond fluorescence upconversion spectroscopy in combination with quantum chemistry TD-DFT calculations ((PCM/TD-PBE0). Experimentally, it was found that when going from water to acetonitrile solution the fluorescence decay of 5FU becomes much faster. The calculations show that this is related to the opening of an additional decay channel in acetonitrile solution since the dark n/pi* excited state becomes near degenerate with the bright pi/pi* state, forming a conical intersection close to the Franck-Condon region. In both solvents, a S1-S0 conical intersection, governed by the out-of-plane motion of the fluorine atom, is active, allowing an ultrafast internal conversion to the ground state.     

The different photoisomerization efficiency of azobenzene in the lowest n pi* and pi pi* singlets: the role of a phantom state

Azobenzene E<==>Z photoisomerization, following excitation to the bright S(pi pi*) state, is investigated by means of ab initio CASSCF optimizations and perturbative CASPT2 corrections. Specifically, by elucidating the S(pi pi*) deactivation paths, we explain the mechanism responsible for azobenzene photoisomerization, the lower isomerization quantum yields observed for the S(pi pi*) excitation than for the S1(n pi*) excitation in the isolated molecule, and the recovery of the Kasha rule observed in sterically hindered azobenzenes. We find that a doubly excited state is a photoreaction intermediate that plays a very important role in the decay of the bright S(pi pi*). We show that this doubly excited state, which is immediately populated by molecules excited to S(pi pi*), drives the photoisomerization along the torsion path and also induces a fast internal conversion to the S1(n pi*) at a variety of geometries, thus shaping (all the most important features of) the S(pi pi*) decay pathway and photoreactivity. We reach this conclusion by determining the critical structures, the minimum energy paths originating on the bright S(pi pi*) state and on other relevant excited states including S1(n pi*), and by characterizing the conical intersection seams that are important in deciding the photochemical outcome. The model is consistent with the most recent time-resolved spectroscopic and photochemical data.     

Photochemistry of pyrrole: time-dependent quantum wave-packet description of the dynamics at the 1pi sigma*-S0 conical intersections

The photoinduced hydrogen-elimination reaction in pyrrole via the conical intersections of the two (1)pi sigma(*) excited states with the electronic ground states [(1)B(1)(pi sigma(*))-S(0) and (1)A(2)(pi sigma(*))-S(0)] have been investigated by time-dependent quantum wave-packet calculations. Model potential-energy surfaces of reduced dimensionality have been constructed on the basis of accurate multireference ab initio electronic-structure calculations. For the (1)B(1)-S(0) conical intersection, the model includes the NH stretching coordinate as the tuning mode and the hydrogen out-of-plane bending coordinate as the coupling mode. For the (1)A(2)-S(0) conical intersection, the NH stretching coordinate and the screwing coordinate of the ring hydrogens are taken into account. The latter is the dominant coupling mode of this conical intersection. The electronic population-transfer processes at the conical intersections, the branching ratio between the dissociation channels, and their dependence on the initial preparation of the system have been investigated for pyrrole and deuterated pyrrole. It is shown that the excitation of the NH stretching mode strongly enhances the reaction rate, while the excitation of the coupling mode influences the branching ratio of different dissociation channels. The results suggest that laser control of the photodissociation of pyrrole via mode-specific vibrational excitation should be possible. The calculations provide insight into the microscopic details of ultrafast internal-conversion processes in pyrrole via hydrogen-detachment processes, which are aborted at the (1)pi sigma(*)-S(0) conical intersections. These mechanisms are of relevance for the photostability of the building blocks of life (e.g., the DNA bases).     

The lowest singlet (n,pi*) and (pi,pi*) excited states of the hydrogen-bonded complex between water and pyrazine

The hydrogen bonding between water and pyrazine in its ground, lowest (n,pi*), and lowest (pi,pi*) states is investigated using density-functional theory (DFT), time-dependent density function theory (TD-DFT), coupled-cluster singles and doubles (CCSD) theory and equation-of-motion coupled cluster (EOM-CCSD) theory. For all states, the minimum-energy configuration is found to be an orthodox linear hydrogen-bonded species, with the bond strength increasing by 0.4 kcal mol-1 upon formation of the (pi,pi*) state and decreasing by 1.0 kcal mol-1 upon formation of the (n,pi*) state. The calculated solvent shifts for the complexes match experimental data and provide a basis for the understanding of the aqueous solvation of pyrazine, and the excited-state complexes are predicted to be only short-lived, explaining the failure of molecular beam experiments to observe them. Quite a different scenario for hydrogen bonding to the (n,pi*) excited state is found compared to those of H2O:pyridine and H2O:pyrimidine: for pyridine linear hydrogen bonds are unstable and hydrogen bonds to the electron-enriched pi cloud are strong, whereas for pyrimidine the excitation localizes on the nonbonded nitrogen leaving the hydrogen-bonding unaffected. For H2O:pyrazine, the (n,pi*) excitation remains largely delocalized, providing a distinct intermediary scenario.     

Photochemical reactivity of 2-vinylbiphenyl and 2-vinyl-1,3-terphenyl: the balance between nonadiabatic and adiabatic photocyclization

A mechanism for the photochemical conversion of 2-vinyl-1,3-terphenyl to 8,9a-dihydrophenanthrene (Lewis, F. D.; Zuo, X.; Gevorgyan, V.; Rubin, M. J. Am. Chem. Soc. 2002, 124, 13664-13665) is presented in this study, based on ab initio restricted active space self-consistent field calculations and a molecular mechanics-valence bond dynamics simulation of a model system: the syn isomer of 2-vinylbiphenyl. An extended crossing seam between the ground and first excited electronic states was found to be largely responsible for the efficient photocyclization of the photochemically active syn isomer. This mechanism is nonadiabatic in nature, with an excited-state reaction pathway approaching the crossing region during the initial stage of cyclization. Dynamics simulation shows that this seam is easily accessible by vibrational motion in the branching space, once a small barrier is passed on the S1 excited-state potential energy surface. Ultrafast radiationless decay to the ground state then follows, and the cyclization is completed on this surface. A second possible mechanism was identified, which involves complete adiabatic cyclization on the S1 surface, with decay to the ground state (at a different conical intersection) only taking place once the product is formed. Thus, there is a competition between these two mechanisms-nonadiabatic and adiabatic-governed by the dynamics of the system. A large quantum yield is predicted for the photocyclization of the syn isomer of 2-vinylbiphenyl and 2-vinyl-1,3-terphenyl, in agreement with experimental observations.     

Singlet excited-state dynamics of 5-fluorocytosine and Cytosine: an experimental and computational study

The photophysics of singlet excited 5-fluorocytosine (5FC) was studied in steady-state and time-resolved experiments and theoretically by quantum chemical calculations. Femtosecond transient absorption measurements show that replacement of the C5 hydrogen of cytosine by fluorine increases the excited-state lifetime by 2 orders of magnitude from 720 fs to 73 +/- 4 ps. Experimental evidence indicates that emission in both compounds originates from a single tautomeric form. The lifetime of 5FC is the same within experimental uncertainty in the solvents ethanol and dimethyl sulfoxide. The insensitivity of the S(1) lifetime to the protic nature of the solvent suggests that proton transfer is not the principal quenching mechanism for the excited state. Excited-state calculations were carried out for the amino-keto tautomer of 5FC, the dominant species in polar environments, in order to understand its longer excited-state lifetime. CASSCF and CAS-PT2 calculations of the excited states show that the minimum energy path connecting the minimum of the (1)pi,pi state with the conical intersection responsible for internal conversion has essentially the same energetics for cytosine and 5FC, suggesting that both bases decay nonradiatively by the same mechanism. The dramatic difference in lifetimes may be due to subtle changes along the decay coordinate. A possible reason may be differences in the intramolecular vibrational redistribution rate from the Franck-Condon active, in-plane modes to the out-of-plane modes that must be activated to reach the conical intersection region.     

Photophysical properties of 2,5-diphenyl-1,6,6a-trithiapentalene revealed by time-resolved spectroscopy

The fluorescence decay from S2(pi, pi*) state of 2,5-diphenyl-1,6,6a-trithiapentalene (DP-TTP) in cyclohexane, tetrahydrofuran and acetonitrile solutions of a quantum yield of approximately 0.02-0.04 were measured. The results indicate that, the dominant process of radiationless deactivation of the S2 state, is internal conversion to the S1 state. Upon laser pulse excitation (lambda(ex) = 532 nm) from the S1(pi, pi*) state, DP-TTP in deoxygenated benzonitrile, acetonitrile, ethanol and tetrahydrofuran solutions give rise to transient triplet triplet absorption (lambdaTmax = 700-720 nm). Kinetic data are presented for intrinsic triplet lifetimes, self-quenching and quenching by oxygen.     

The guanine tautomer puzzle: quantum chemical investigation of ground and excited states

Combined density functional and multireference configuration interaction methods have been employed to explore the ground and low-lying electronically excited states of the most important tautomeric and rotameric forms of guanine with the purpose of resolving the conflicting assignments of IR-UV bands found in the literature. The calculations predict sharp 1(pi-->pi*) origin transitions for the RN1 rotamer of the 7H-amino-hydroxy species and the RN7 rotamer of the 9H-amino-hydroxy species. The other 9H-amino-hydroxy rotamer, RN1, undergoes ultrafast nonradiative decay and is thus missing in the UV spectra. Because of its very small Franck-Condon factor and the presence of a conical intersection close by, it appears questionable, whether the 1(pi-->pi*) origin transition of 9H-amino-oxo-guanine can be observed experimentally. Vibrational overlap is more favorable for the 1(pi-->pi*) origin transition of the 7H- amino-oxo form, but also this tautomer is predicted to undergo ultrafast nonradiative decay of the 1(pi-->pi*) population. The good agreement of calculated IR frequencies of the amino-oxo species with recent IR spectra in He droplets and their mismatch with peaks observed in IR-UV spectra indicate that none of the bands stem from 7H- or 9H-amino-oxo guanine. Instead, our results suggest that these bands originate from 7H-imino-oxo guanine tautomers. In the excited-state dynamics of the biologically relevant 9H-amino-oxo tautomer, a diffuse charge transfer state is predicted to play a significant role.     

Energetics of cytosine singlet excited-state decay paths--a difficult case for CASSCF and CASPT2

Three deactivation paths for singlet excited cytosine are calculated at the CASPT2//CASSCF (complete active space second-order perturbation//complete active space self-consistent field) level of theory, using extended active spaces that allow for a reliable characterization of the paths and their energies. The lowest energy path, with a barrier of approximately 0.1 eV, corresponds to torsion of the C5-C6 bond, and the decay takes place at a conical intersection analogous to the one found for ethylene and its derivatives. There is a further path with a low energy barrier of approximately 0.2 eV associated with the (n(N),pi*) state which could also be populated with a low energy excitation. The path associated with a conical intersection between the ground and (n(O),pi*) states is significantly higher in energy (> 1 eV). The presence of minima on the potential energy surface for the (n,pi*) states that could contribute to the biexponential decay found in the gas phase was investigated, but could not be established unequivocally.     

A three-state model for the photophysics of adenine

An ab initio theoretical study at the CASPT2 level is reported on minimum energy reaction paths, state minima, transition states, reaction barriers, and conical intersections on the potential energy hypersurfaces of two tautomers of adenine: 9H- and 7H-adenine. The obtained results led to a complete interpretation of the photophysics of adenine and derivatives, both under jet-cooled conditions and in solution, within a three-state model. The ultrafast subpicosecond fluorescence decay measured in adenine is attributed to the low-lying conical intersection (gs/pipi* La)(CI), reached from the initially populated 1(pipi* La) state along a path which is found to be barrierless only in 9H-adenine, while for the 7H tautomer the presence of an intermediate plateau corresponding to an NH2-twisted conformation may explain the absence of ultrafast decay in 7-substituted compounds. A secondary picosecond decay is assigned to a path involving switches towards two other states, 1(pipi* Lb) and 1(npi*), ultimately leading to another conical intersection with the ground state, (gs/npi*), with a perpendicular disposition of the amino group. The topology of the hypersurfaces and the state properties explain the absence of secondary decay in 9-substituted adenines in water in terms of the higher position of the 1(npi*) state and also that the 1(pipi* Lb) state of 7H-adenine is responsible for the observed fluorescence in water. A detailed discussion comparing recent experimental and theoretical findings is given. As for other nucleobases, the predominant role of a pipi*-type state in the ultrafast deactivation of adenine is confirmed.     

A theoretical characterization of the photoisomerization channels of 1,2-cyclononadienes on both singlet and triplet potential-energy surfaces

The ground-, (1)(pipi*)-, and (3)(pipi*)-state potential-energy surfaces of 1,2-cyclononadiene and isomeric C(9)H(14) species, as well as 1-methyl-1,2-cyclononadiene and isomeric C(10)H(16) species were all mapped using CASSCF and the 6-31G(d) basis set. Theoretical results were found to be in good agreement with the available experimental observations for both 1,2-cyclononadiene and 1-methyl-1,2-cyclononadiene isomerization reactions under singlet and triplet direct or sensitized irradiation. Extremely efficient decay occurs from the first singlet excited state to the ground state through at least three different conical intersections (surface crossings). The first of these crossing points is accessed by a one-bond ring closure. From this conical intersection point (CI-A or CI-C), some possible subsequent ground-state reaction paths have been identified: 1) intramolecular C--H bond insertion to form the bicyclic photoproduct and 2) intramolecular C--H bond insertion to form tricyclic photoproducts. An excited state [1,3]-sigmatropic shift leads to the second conical intersection (CI-B or CI-E), which can give a three-bond cyclononyne species. Besides these, in the singlet photochemical reactions of 1-methyl-1,2-cyclononadiene, excited-state, one allenic C--H bond insertion leads to a third conical intersection (CI-D). Possible ground-state reaction pathways from this structure lead to the formation of a diene photoproduct or to transannular insertion photoproducts. Moreover, in the case of triplet 1,2-cyclononadiene and 1-methyl-1,2-cyclononadiene photoisomerization reactions, both chemical reactions will adopt a 1,3-biradical (T(1)/S(0)-1, T(1)/S(0)-2, and T(1)/S(0)-3), which may undergo intersystem crossings leading to the formation of tricyclic or bicyclic photoproducts. The results obtained allow a number of predictions to be made.     

Ultrafast deactivation of an excited cytosine-guanine base pair in DNA

Multiconfigurational ab initio calculations and QM/MM molecular dynamics simulations of a photoexcited cytosine-guanine base pair in both gas phase and embedded in the DNA provide detailed structural and dynamical insights into the ultrafast radiationless deactivation mechanism. Photon absorption promotes transfer of a proton from the guanine to the cytosine. This proton transfer is followed by an efficient radiationless decay of the excited state via an extended conical intersection seam. The optimization of the conical intersection revealed that it has an unusual topology, in that there is only one degeneracy-lifting coordinate. This is the central mechanistic feature for the decay both in vacuo and in the DNA. Radiationless decay occurs along an extended hyperline nearly parallel to the proton-transfer coordinate, indicating the proton transfer itself is not directly responsible for the deactivation. The seam is displaced from the minimum energy proton-transfer path along a skeletal deformation of the bases. Decay can thus occur anywhere along the single proton-transfer coordinate, accounting for the remarkably short excited-state lifetime of the Watson-Crick base pair. In vacuo, decay occurs after a complete proton transfer, whereas in DNA, decay can also occur much earlier. The origin of this effect lies in the temporal electrostatic stabilization of dipole in the charge-transfer state in DNA.     

Radiationless decay mechanism of cytosine: an ab initio study with comparisons to the fluorescent analogue 5-methyl-2-pyrimidinone

The ultrafast radiationless decay mechanism of photoexcited cytosine has been theoretically supported by exploring the important potential energy surfaces using multireference configuration-interaction ab initio methods for the gas-phase keto-tautomer free base. At vertical excitation, the bright state is S1 (pipi*) at 5.14 eV, with S2 (nNpi*) and S3 (nOpi*) being dark states at 5.29 and 5.93 eV, respectively. Minimum energy paths connect the Franck-Condon region to a shallow minimum on the pipi* surface at 4.31 eV. Two different energetically accessible conical intersections with the ground state surface are shown to be connected to this minimum. One pathway involves N3 distorting out of plane in a sofa conformation, and the other pathway involves a dihedral twist about the C5-C6 bond. Each of these pathways from the minimum contains a low barrier of 0.14 eV, easily accessed by low vibronic levels. The path involving the N3 sofa distortion leads to a conical intersection with the ground state at 4.27 eV. The other pathway leads to an intersection with the ground state at 3.98 eV, lower than the minimum by about 0.3 eV. Comparisons with our previously reported study of the fluorescent cytosine analogue 5-methyl-2-pyrimidinone (5M2P) reveal remarkably similar conformational distortions throughout the decay pathways of both bases. The different photophysical behavior between the two molecules is attributed to energetic differences. Vertical excitation in cytosine occurs at a much higher energy initially, creating more vibrational energy than 5M2P in the Franck-Condon region, and the minimum S1 energy for 5M2P is too low to access an intersection with the ground state, causing population trapping and fluorescence. Calculations of vertical excitation energies of 5-amino-2-pyrimidinone and 2-pyrimidinone reveal that the higher excitation energy of cytosine is likely due to the presence of the amino group at the 4-position.     

The decay from the dark npi* excited state in uracil: an integrated CASPT2/CASSCF and PCM/TD-DFT study in the gas phase and in water

By integrating the results of MS-CASPT2/CASSCF and TD-PBE0 calculations, we propose a mechanism for the decay of the excited dark state in pyrimidine, fully consistent with all the available experimental results. An effective conical intersection (CI-npi) exists between the spectroscopic pi/pi* excited state (Spi) and a dark n/pi* state (Sn), and a fraction of the population decays to the minimum of Sn (Sn-min). The conical intersection between Sn and the ground-state is not involved in the decay mechanism, because of its high energy gap with respect to Sn-min. On the other hand, especially in hydrogen bonding solvents, the energy gap between Sn-min and CI-npi is rather small. After thermalization in Sn-min, the system can thus recross CI-npi and then quickly proceed on the Spi barrierless path toward the conical intersection with the ground state.