Ultrafast internal conversion of excited cytosine via the lowest pipi electronic singlet state

Computational evidence at the CASPT2 level supports that the lowest excited state pipi* contributes to the S1/S0 crossing responsible for the ultrafast decay of singlet excited cytosine. The computed radiative lifetime, 33 ns, is consistent with the experimentally derived value, 40 ns. The nOpi* state does not play a direct role in the rapid repopulation of the ground state; it is involved in a S2/S1 crossing. Alternative mechanisms through excited states pisigma* or nNpi* are not competitive in cytosine.     

Unified model for the ultrafast decay of pyrimidine nucleobases

Ultrafast decay processes detected after absorption of UV radiation in gas-phase pyrimidine nucleobases uracil, thymine, and cytosine are ascribed to the barrierless character of the pathway along the low-lying 1(pipi*) hypersurface connecting the Franck-Condon region with an out-of-plane distorted ethene-like conical intersection with the ground state. Longer lifetime decays and low quantum yield emission are on the other hand related to the presence of a 1(pipi*) state planar minimum on the S1 surface and the barriers to access other conical intersections. A unified model for the three systems is established on the basis of accurate multiconfigurational CASPT2 calculations, whereas the effect of the different levels of theory on the results is carefully analyzed.     

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.     

Quantum chemical investigation of the electronic spectra of the keto, enol, and keto-imine tautomers of cytosine

The low-lying excited singlet states of the keto, enol, and keto-imine tautomers of cytosine have been investigated employing a combined density functional/multireference configuration interaction (DFT/MRCI) method. Unconstrained geometry optimizations have yielded out-of-plain distorted structures of the pi --> pi and n --> pi excited states of all cytosine forms. For the keto tautomer, the DFT/MRCI adiabatic excitation energy of the pi --> pi state (4.06 eV including zero-point vibrational energy corrections) supports the resonant two-photon ionization (R2PI) spectrum (Nir et al. Phys. Chem. Chem. Phys. 2002, 5, 4780). On its S1 potential energy surface, a conical intersection between the 1pipi state and the electronic ground state has been identified. The barrier height of the reaction along a constrained minimum energy path amounts to merely 0.2 eV above the origin and explains the break-off of the R2PI spectrum. The 1pipi minimum of the enol tautomer is found at considerably higher excitation energies (4.50 eV). Because of significant geometry shifts with respect to the ground state, long vibrational progressions are expected, in accord with experimental observations. For the keto-imine tautomer, a crossing of the 1pipi potential energy surface with the ground-state surface has been found, too. Its n --> pi minimum (3.27 eV) is located well below the conical intersection between the pi --> pi and S0 states, but it will be difficult to observe because of its small transition moment. The identified conical intersections of the pi --> pi excited states of the keto cytosine tautomers are made responsible for the ultrafast decay to the electronic ground states and thus may explain their subpicoseconds lifetimes.     

On the origin of the ultrafast internal conversion of electronically excited pyrimidine bases

The ultrafast radiationless decay of photoexcited uracil and cytosine has been investigated by ab initio quantum chemical methods based on CIS and CR-EOM-CCSD(T) electronic energy calculations at optimized CIS geometries. The calculated potential energy profiles indicate that the S(1) --> S(0) internal conversion of the pyrimidine bases occurs through a barrierless state switch from the initially excited (1)pipi state to the out-of-plane deformed excited state of biradical character, which intersects the ground state at a lower energy. This three-state nonradiative decay mechanism predicts that replacement of the C5 hydrogen by fluorine introduces an energy barrier for the initial state switch, whereas replacement of the C6 hydrogen by fluorine does not. These predictions are borne out by the very different fluorescence yields of 5-fluorinated bases relative to the corresponding 6-fluorinated bases. It is concluded from these results that the origin of the ultrafast radiationless decay is the same for the two pyrimidine bases.     

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 study on deactivation pathways of excited 9H-guanine

The complete active space with second-order perturbation theory/complete active space self-consistent-field method was used to explore the nonradiative decay mechanism for excited 9H-guanine. On the 1pipi* (1L(a)) surface we determined a conical intersection (CI), labeled (S0pipi*)(CI), between the 1pipi* (1L(a)) excited state and the ground state, and a minimum, labeled (pipi*)min. For the 1pipi* (1L(a)) state, its probable deactivation path is to undergo a spontaneous relaxation to (pipi*)min first and then decay to the ground state through (S0pipi*)(CI), during which a small activation energy is required. On the 1n(N)pi* surface a CI between the 1n(N)pi* and 1pipi* (1L(a)) states was located, which suggests that the 1n(N)pi* excited state could transform to the 1pipi* (1L(a)) excited state first and then follow the deactivation path of the 1pipi* (1L(a)) state. This CI was also possibly involved in the nonradiative decay path of the second lowest 1pipi* (1L(b)) state. On the 1n(O)pi* surface a minimum was determined. The deactivation of the 1n(O)pi* state to the ground state was estimated to be energetically unfavorable. On the 1pisigma* surface, the dissociation of the N-H bond of the six-membered ring is difficult to occur due to a significant barrier.     

A theoretical investigation of P-hydroxyphenacyl caged phototrigger compounds: how water induces the photodeprotection and subsequent rearrangement reactions

Complete active-space self-consistent field (CASSCF) calculations with a (14,11) active space and density functional theory calculations followed by Car-Parrinello molecular dynamic simulations are reported for the p-hydroxyphenacyl acetate, diethyl phosphate, and diphenyl phosphate phototrigger compounds. These calculations considered the explicit hydrogen bonding of water molecules to the phototrigger compound and help reveal the role of water in promoting the photodeprotection and subsequent rearrangement reactions for the p-hydroxyphenacyl caged phototrigger compounds experimentally observed in the presence of appreciable amounts of water but not observed in neat nonproton solvents like acetonitrile. The 267 nm excitation of the phototrigger compounds leads to an instantaneous population of the S3(1pipi*) state Franck-Condon region, which is followed by an internal conversion deactivation route to the S1(1npi*) state via a 1pipi*/1npi* vibronic coupling. The shorter lifetime of the S1(1npi*) state (approximately 1 ps) starting from the FC geometry is terminated by a fast intersystem crossing at a 3pipi*/3npi* intersection with a structure of mixed pipi*/npi* excitation in the triplet state. The deprotection reaction is triggered by a proton (or hydrogen atom) transfer assisted by water bridges and emanates from this pipi*/npi* triplet state intersection. With the departure of the leaving group, the reaction evolves into a water-mediated post-deprotection phase where the spin inversion of pQM (X, 3A) leads to a spiroketone in the ground state by a cyclization process that is followed by an attack of water to produce a 1,1'-di-hydroxyl-spiroketone. Finally, the H atom of the hydroxyl in 1,1'-di-hydroxyl-spiroketon transfers back to the p-O atom aided by water molecules to generate the p-hydroxyphenyl-acetic acid final rearrangement product.     

Theoretical studies of the photochemical dynamics of acetylacetone: isomerzation, dissociation, and dehydration reactions

The potential energy surfaces of the C-O cleavage, rotational isomerization, keto-enolic tautomerization, and dehydration reactions of acetylacetone in the lowest triplet and ground states have been determined using the complete active space self-consistent field and density functional theory methods. The main photochemical mechanism obtained indicates that the acetylacetone molecule in the S(2)((1)pipi*) state can relax to the T(1)((3)pipi*) state via the S(2)-S(1) vibronic interaction and an S(1)/T(1)/T(2) intersection. The C-O fission pathway is the predominant dissociation process in the T(1)((3)pipi) state. Rotational isomerization reactions proceed difficultly in the ground state but very easily in the T(1)((3)pipi*) state. Keto-enolic tautomerization takes place with little probability for acetylacetone in the gas phase.     

Photostability versus photodegradation in the excited-state intramolecular proton transfer of nitro enamines: competing reaction paths and conical intersections

The phototautomerization mechanism of a model nitro enamine (NEA) chromophore (incorporated in the structure of a highly photolabile pesticide, tetrahydro-2-(nitromethylene)-2H-1,3-thiazine) has been studied using complete active space self-consistent field reaction path computations. The optically accessible 1pipi* excited state of NEA involves separation of charge and correlates diabatically with the ground state of the tautomerized acinitro imine (ANI) form. For optimum photostabilization, the 1pipi* state of NEA should be S1: in this case, the tautomer would be efficiently formed via a diabatic intramolecular proton-transfer pathway passing through an S1/S0 conical intersection, followed by a facile thermal back proton-transfer reaction. However, in NEA itself the lowest excited states correspond to nitro group 1npi* states, and there are additional surface crossings that provide a mechanism for populating the 1npi* manifold. The above results indicate that the high photolability observed for the pesticide [Kleier, D.; Holden, I.; Casida, J. E.; Ruzo, L. O. J. Agric. Food Chem. 1985, 33, 998-1000] has to be ascribed to photochemistry originating on the 1npi* manifold of states, populated indirectly from the 1pipi* state.     

Direct observation of electronic relaxation dynamics in adenine via time-resolved photoelectron spectroscopy

We present femtosecond time-resolved photoelectron spectra of adenine in a molecular beam, recorded at pump wavelengths of 250, 267, and 277 nm. This leads to initial excitation of the bright S2(pipi*). Close to the band origin (277 nm), the lifetime is several picoseconds. Higher vibronic levels (267 and 250 nm excitation) show much shorter lifetimes of t < 50 fs, and we observe strong coupling between S2(pipi*) and S1(npi*). Rapid internal conversion (t < 50 fs) populates the lower lying S1(npi*) state which has a lifetime of 750 fs. At 267 nm, we found evidence for an additional channel which is consistent with the dissociative S3(pisigma*) state, previously proposed as an ultrafast relaxation pathway from S2(pipi*).     

Theoretical study toward understanding ultrafast internal conversion of excited 9H-adenine

The CASPT2/CASSCF method with the 6-311G basis set and an active space up to (14, 11) was used to explore the ultrafast internal conversion mechanism for excited 9H-adenine. Three minima, two transition states, and seven conical intersections were obtained to build up the two deactivation pathways for the internal conversion mechanism. Special efforts were made to explore the excited-state potential energy surfaces near the Franck-Condon region and determine the various barriers in the processes of deactivation. The barrier required from the 1pipi (1La) state to deactivate nonradiatively is found to be lower than that required from the 1pipi (1Lb) state. On 250 nm excitation, the 1pipi (1La) state is populated, and the transition from 1pipi (1La) to the lowest 1npi state involves very low barriers, which may account for the observed short (<50 fs) lifetime of the 1pipi excited state. The deactivation of the lowest 1npi state is required to overcome a barrier of 3.15 kcal/mol, which should be responsible for the 750 fs lifetime of the npi excited state. On 267 nm excitation, the vibrationally active 1pipi (1Lb) state is populated. Excitation at 277 nm prepares the 1pipi (1Lb) state without much excessive vibrational energy, which may be responsible for the observed >2 ps lifetime.     

Conical intersections in thymine

The mechanisms which are responsible for the radiationless deactivation of the npi* and pipi* excited singlet states of thymine have been investigated with multireference ab initio methods (the complete-active-space self-consistent-field (CASSCF) method and second-order perturbation theory with respect to the CASSCF reference (CASPT2)) as well as with the CC2 (approximated singles and doubles coupled-cluster) method. The vertical excitation energies, the equilibrium geometries of the 1npi*and 1pipi* states, as well as their adiabatic excitation energies have been determined. Three conical intersections of the S1 and S0 energy surfaces have been located. The energy profiles of the excited states and the ground state have been calculated with the CASSCF method along straight-line reaction paths leading from the ground-state equilibrium geometry to the conical intersections. All three conical intersections are characterized by strongly out-of-plane distorted geometries. The lowest-energy conical intersection (CI1) arises from a crossing of the lowest 1pipi* state with the electronic ground state. It is found to be accessible in a barrierless manner from the minimum of the 1pipi* state, providing a direct and fast pathway for the quenching of the population of the lowest optically allowed excited states of thymine. This result explains the complete diffuseness of the absorption spectrum of thymine in supersonic jets. The lowest vibronic levels of the optically nearly dark 1npi* state are predicted to lie below CI1, explaining the experimental observation of a long-lived population of dark excited states in gas-phase thymine.     

The fluorescence mechanism of 5-methyl-2-pyrimidinone: an ab initio study of a fluorescent pyrimidine analog

The photophysically important potential energy surfaces of the fluorescent pyrimidine analog 5-methyl-2-pyrimidinone have been explored using multireference configuration-interaction ab initio methods at three levels of dynamical correlation, all of which support a fluorescence mechanism. At vertical excitation S1 (dark, n(N)pi*) and S2 (bright, pipi*) are almost degenerate at 4.4 eV, with S3 (dark, n(O)pi*) at 5.1 eV. The excited system can follow the S1-S2 seam of conical intersections, accessible from the Franck-Condon region, to its minimum and then evolve from this conical intersection on the S1 (pipi*) surface to a global minimum. At lower levels of correlation, the S1 surface shows two minima separated by a barrier of up to 0.18 eV. The secondary minimum found at the lower levels of correlation becomes the global minimum with higher correlation. The S1 population at this minimum can be trapped from accessing the lowest energy S0-S1 (pipi*/gs) conical intersection by an energy gap at least 0.3-0.4 eV higher than the S1 minimum. The calculated emission energy from this minimum is 2.80 eV. Gradient pathways connecting important S1 geometries are presented, as well as other excited state conical intersections.     

Ultrafast nonradiative dynamics in electronically excited hexafluorobenzene by femtosecond time-resolved mass spectrometry

The fast nonradiative decay dynamics of the lowest two excited pipi(*) electronic states (S(2) and S(3)) of hexafluorobenzene have been investigated by using femtosecond time-resolved time-of-flight mass spectrometry. The molecules were excited at wavelengths between 265 nm > or = lambda(pump) > or = 217 nm and probed by four- and three-photon ionization at lambda(probe)=775 nm. The observed temporal profiles exhibit two exponential decay times (tau(1)=0.54-0.1 ps and tau(2)=493-4.67 ps, depending on the excitation wavelength) and a superimposed coherent oscillation with vibrational frequency nu(osc)=97 cm(-1) and damping time tau(D) that is two to three times longer than the respective tau(1). The first decay component (tau(1)) is assigned to rapid radiationless transfer from the excited optically bright pipi(*) electronic state (S(2) or S(3), respectively) through a conical intersection (CI) to the lower-lying optically dark pisigma(*) state (S(1)) of the molecule; the second component (tau(2)) is attributed to the subsequent slower relaxation from the S(1) state back to the electronic ground state (S(0)). tau(2) dramatically decreases with increasing vibronic excitation energy up to the CI connecting the pisigma(*) with the S(0) state. The coherent oscillation is identified as nuclear motion along the out-of-plane vibration nu(16a) (notation as for benzene), which has e(2u) symmetry and acts as coupling mode between the pipi(*) and pisigma(*) states.     

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.     

Environmental photochemistry of nitro-PAHs: direct observation of ultrafast intersystem crossing in 1-nitropyrene

Femtosecond broadband transient absorption experiments of 1-nitropyrene, a nitro-polycyclic aromatic hydrocarbon of environmental concern are presented in cyclohexane and hexane solutions. The transient absorption spectra show the presence of three species that are assigned to the Franck-Condon excited lowest singlet (S1) state, the structurally relaxed S1 state, and the lowest excited triplet state. The spectral changes at early times are interpreted in terms of conformational dynamics; primarily due to an ultrafast rotation of the nitro group in the S1 state. This excited state relaxation is followed by intersystem crossing with a time constant of 7 ps. CIS/6-31G(d,p) calculations predict planarization of the nitro-aromatic torsional angle as the major nuclear relaxation coordinate, from 32.8 degrees at the HF/6-31G(d,p) level of theory in the ground state (27.46 degrees at B3LYP/6-31++G(d,p)) to 0.07 degrees in the S1 state. Vertical excitation energies at the TDDFT/6-31++G(d,p) and TDDFT/IEFPCM/6-31++G(d,p) levels of theory predict a small energy gap (<0.12 eV) between the S1(pipi*) state and the third excited triplet state T3(npi*) in the gas phase and in cyclohexane, respectively. The small energy gap suggests a large spin-orbit coupling between the S1(pipi*) and T3(npi*) states, which explains the ultrafast intersystem crossing of 1-nitropyrene in nonpolar solvents.     

Mechanism and dynamics of azobenzene photoisomerization

The excited-state dynamics of trans-azobenzene were investigated by femtosecond time-resolved photoelectron spectroscopy and ab initio molecular dynamics. Two near-degenerate pipi* excited states, S2 and S3,4, were identified in a region hitherto associated with only one excited state. These results help to explain contradictory reports about the photoisomerization mechanism and the wavelength dependence of the quantum yield. A new model for the isomerization mechanism is proposed.     

Properties of a new fluorescent cytosine analogue, pyrrolocytosine

Pyrrolocytosine is a novel, environment sensitive, fluorescent base that can be used in place of cytosine as a fluorescent marker in nucleic acids. In this work the results of a detailed computational investigation into the hybridization and photochemical properties of the base are reported. The interaction energy of the base pair formed between pyrrolocytosine and guanine, calculated at the MP2/6-31G(d,0.25)//HF/6-31G(d,p) level, was found to be -27.2 kcal mol(-1), comparing very favorably with the value calculated for the cytosine and guanine base pair, -25.8 kcal mol(-1). The wavelengths for the vertical transitions of pyrrolocytosine and cytosine were determined using both the configuration interaction technique, with singly excited configurations (CIS) and time-dependent density functional theory using the B3LYP functional (TDB3LYP). It was found that the spacing between the first pipi state and the first npi state was significantly larger in the case of pyrrolocytosine than cytosine, providing a rationale for the higher fluorescence quantum yield of the former. Hydrogen bonding of pyrrolocytosine to guanine did not affect the predicted fluorescence properties of pyrrolocytosine whereas stacking guanine above pyrrolocytosine, in a manner appropriate to B-form DNA, significantly reduced the predicted fluorescence. Calculations on the two base systems using the TDB3LYP method produced low-lying charge-transfer states which are not predicted when the CIS method is used and are not thought to be physically meaningful.