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.     

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.     

A three-state model for the photophysics of guanine

The nonadiabatic photochemistry of the guanine molecule (2-amino-6-oxopurine) and some of its tautomers has been studied by means of the high-level theoretical ab initio quantum chemistry methods CASSCF and CASPT2. Accurate computations, based by the first time on minimum energy reaction paths, states minima, transition states, reaction barriers, and conical intersections on the potential energy hypersurfaces of the molecules lead to interpret the photochemistry of guanine and derivatives within a three-state model. As in the other purine DNA nucleobase, adenine, the ultrafast subpicosecond fluorescence decay measured in guanine is attributed to the barrierless character of the path leading from the initially populated 1(pi pi* L(a)) spectroscopic state of the molecule toward the low-lying methanamine-like conical intersection (gs/pi pi* L(a))CI. On the contrary, other tautomers are shown to have a reaction energy barrier along the main relaxation profile. A second, slower decay is attributed to a path involving switches toward two other states, 1(pi pi* L(b)) and, in particular, 1(n(O) pi*), ultimately leading to conical intersections with the ground state. A common framework for the ultrafast relaxation of the natural nucleobases is obtained in which the predominant role of a pi pi*-type state is confirmed.     

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.     

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).     

Ultrafast nonradiative decay by hypoxanthine and several methylxanthines in aqueous and acetonitrile solution

Excited state lifetimes of hypoxanthine and four methylxanthine compounds (paraxanthine, theophylline, theobromine, and caffeine) were studied by femtosecond transient absorption spectroscopy in aqueous and acetonitrile solution. Transient absorption signals recorded at visible and UV probe wavelengths reveal that internal conversion from the photoexcited state to the electronic ground state takes place in water on the hundreds of femtoseconds time scale. Excited-state relaxation occurs approximately threefold more slowly in acetonitrile solution than in water. Structural considerations suggest that the deactivating conical intersection for the methylxanthines differs from that responsible for nonradiative decay by hypoxanthine, adenine, and guanine. All compounds studied have ultrashort excited state lifetimes similar to those of adenine and guanine, suggesting that these xanthine derivatives could have been photostable building blocks in prebiotic environments exposed to intense UV radiation.     

Theoretical study of the photoinduced transfer among the ground state and two metastable states in [Fe(CN)5NO]2-

Ab initio calculations were performed to investigate photoinduced transfers among the ground state (GS) and two metastable states (MS1 and MS2) of [Fe(CN)5NO]2-. We obtained the global potential energy surface of the electronic ground state by a scheme of multireference singly and doubly excited configuration interaction followed by a Davidson-type quadruple correction (MRSDCI+Q). The ground state surface has three local minima corresponding to GS, MS1, and MS2. The character of bond between Fe and the nitrosyl group are discussed. We carried out calculations of the lower five electronic excited states by MRSDCI+Q. The main configurations of these lower five excited states were represented by the dFe-->pi*NO transition accompanied by considerable back-donation. The potential energy surfaces of the six states, including the ground state, were obtained by state averaged complete active space self-consistent field calculations. The surfaces have several conical intersections and avoided crossings in the reaction pathway. The photoinduced transfers among GS, MS1, and MS2 are caused by the nonadiabatic effect near these crossings.     

Near-UV resonant two-photon ionization spectroscopy of gas phase guanine: evidence for the observation of three rare tautomers

In light of a recently published study on the IR spectroscopy of guanine in He droplets (Choi, M. Y.; Miller, R. E. J. Am. Chem. Soc. 2006, 128, 7320), the present letter proposes a new interpretation of the resonant two-photon ionization (R2PI) experiments on gas phase guanine, which is supported by quantum chemistry calculations. Whereas He droplet experiments detect the most stable forms, only one of these forms is observed (very marginally) in the R2PI spectrum, which is actually dominated by three less stable "rare" tautomers, whose stabilities lie in the 3-7 kcal/mol range. The absence of the most stable forms in the R2PI spectrum suggests that a tautomer-dependent ultrafast relaxation process takes place in the excited state of these stable tautomers. The present reinterpretation modifies qualitatively the picture of the excited state of guanine tautomers and should contribute to the understanding of the deactivation mechanisms taking place in the excited state of DNA bases.     

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.     

Theoretical analysis of photoinduced H-atom elimination in thiophenol

The photoinduced hydrogen elimination reaction in thiophenol via the conical intersections of the dissociative (1)πσ? excited state with the bound (1)ππ? excited state and the electronic ground state has been investigated with ab initio electronic-structure calculations and time-dependent quantum wave-packet calculations. A screening of the coupling constants of the symmetry-allowed coupling modes at the (1)ππ?-(1)πσ? and (1)πσ?-S(0) conical intersection shows that the SH torsional mode is by far the most important coupling mode at both conical intersections. A model including three intersecting potential-energy surfaces (S(0), (1)ππ?, (1)πσ?) and two nuclear degrees of freedom (SH stretch and SH torsion) has been constructed on the basis of ab initio complete-active-space self-consistent field and multireference second-order perturbation theory calculations. The nonadiabatic quantum wave-packet dynamics initiated by optical excitation of the (1)ππ? and (1)πσ? states has been explored for this three-state two-coordinate model. The photodissociation dynamics is characterized in terms of snapshots of time-dependent wave packets, time-dependent electronic population probabilities, and the branching ratio of the (2)σ/(2)π electronic states of the thiophenoxyl radical. The dependence of the timescale of the photodissociation process and the branching ratio on the initial excitation of the SH stretching and SH torsional vibrations has been analyzed. It is shown that the node structure, which is imposed on the nuclear wave packets by the initial vibrational preparation as well as by the transitions through the conical intersections, has a profound effect on the photodissociation dynamics. The effect of additional weak coupling modes of CC twist (ν(16a)) and ring-distortion (ν(16b)) character has been investigated with three-dimensional and four-dimensional time-dependent wave-packet calculations, and has been found to be minor.     

Excited State Photophysics of Curcumin and its Modulation in Alkaline Non-Aqueous Medium

Curcumin, a well-known medicinal pigment, has seen limited applications in biology despite having great potential as a therapeutic drug. Deprotonation is one of the possible ways to enhance solubility of curcumin in polar solvent. Here, we have explored the effect of deprotonation on the ultrafast dynamics of this biomolecule with the help of the time-resolved fluorescence spectroscopic measurements using the femtosecond fluorescence upconversion technique. The excited state photophysics of fully deprotonated curcumin significantly differs from that of neutral curcumin. We have observed that the completely deprotonated curcumin not only has higher quantum yield, but also higher excited state lifetime and slower solvation dynamics in comparison to neutral curcumin. We propose solvation dynamics and intramolecular charge transfer as the excited state processes associated with the radiative decay of the completely deprotonated molecule, while ruling out the possibility of excited state proton exchange or proton transfer. Our results are well supported by time-dependent density-functional theory calculations. Lastly, we have also demonstrated the possibility of modulating the ultrafast dynamics of fully deprotonated curcumin using non-aqueous alkaline binary solvent mixtures. We believe our results will provide significant physical insight towards unveiling the excited state dynamics of this molecule.     

Distinguishing between relaxation pathways by combining dissociative ionization pump probe spectroscopy and ab initio calculations: a case study of cytosine

We present a general method for tracking molecular relaxation along different pathways from an excited state down to the ground state. We follow the excited state dynamics of cytosine pumped near the S(0)-S(1) resonance using ultrafast laser pulses in the deep ultraviolet and probed with strong field near infrared pulses which ionize and dissociate the molecules. The fragment ions are detected via time of flight mass spectroscopy as a function of pump probe delay and probe pulse intensity. Our measurements reveal that different molecular fragments show different timescales, indicating that there are multiple relaxation pathways down to the ground state. We interpret our measurements with the help of ab initio electronic structure calculations of both the neutral molecule and the molecular cation for different conformations en route to relaxation back down to the ground state. Our measurements and calculations show passage through two seams of conical intersections between ground and excited states and demonstrate the ability of dissociative ionization pump probe measurements in conjunction with ab initio electronic structure calculations to track molecular relaxation through multiple pathways.     

Solvent effect on conical intersections in excited-state 9H-adenine: radiationless decay mechanism in polar solvent

The ground- and excited-state free energy minima and the conical intersections among these states of 9H-adenine in aqueous and acetonitrile solutions are studied theoretically to elucidate the mechanism of radiationless decay. We employ the recently proposed linear-response free energy (LRFE) to locate the energy minima and conical intersections in solution. The LRFE is calculated by using the reference interaction site model self-consistent field method. The geometry optimizations are carried out at the complete active space self-consistent field level, and the dynamic electron correlation energies are estimated by the multireference M?ller-Plesset method. We find that the conical intersection between the (1)L(a) and (1)L(b) states in aqueous solution occurs at a wide area of the free energy surface, indicating a strong vibronic coupling between them. On the other hand, the (1)npi(*) state is largely blue-shifted at planar geometries in solution, which suggests that the nonadiabatic transition to this state is suppressed. The importance of the (1)pisigma(*) channel is also examined in both the gas phase and solution. Based on the free energy characteristics obtained by the calculations, we intend to explain the experimental observations that the excited state of 9H-adenine decays monoexponentially with shorter lifetimes in polar solvents than that in the gas phase.     

Ultrafast dynamics in thiophene investigated by femtosecond pump probe photoelectron spectroscopy and theory

A hybrid of a time-of-flight mass spectrometer and a time-of-flight "magnetic-bottle type" photoelectron (PE) spectrometer is used for fs pump-probe investigations of the excited state dynamics of thiophene. A resonant two-photon ionization spectrum of the onset of the excited states has been recorded with a tunable UV laser of 190 fs pulse width. With the pump laser set to the first intense transition we find by UV probe ionization first a small time shift of the maxima in the PE spectrum and then a fast decay to a low constant intensity level. The fitted time constants are 80+/-10 fs, and 25+/-10 fs, respectively. Theoretical calculations show that upon geometry relaxation the electronic state order changes and conical intersections between excited states exist. We use the vertical state order S1, S2, S3 to define the terms S1, S2, and S3 for the characterization of the electron configuration of these states. On the basis of our theoretical result we discuss the electronic state order in the UV spectra and identify in the photoelectron spectrum the origin of the first cation excited state D1. The fast excited state dynamics agrees best with a vibrational dynamics in the photo-excited S1 (80+/-10 fs) and an ultrafast decay via a conical intersection, presumably a ring opening to the S3 state (25+/-10 fs). The subsequently observed weak constant signal is taken as an indication, that in the gas phase the ring-closure to S0 is slower than 50 ps. An ultrafast equilibrium between S1 and S2 before ring opening is not supported by our data.     

Conical intersections involving the dissociative 1pisigma* state in 9H-adenine: a quantum chemical ab initio study

The conical intersections of the dissociative 1pisigma* excited state with the lowest 1pipi* excited state and the electronic ground state of 9H-adenine have been investigated with multireference electronic structure calculations. Adiabatic and quasidiabatic potential energy surfaces and coupling elements were calculated as a function of the NH stretch coordinate of the azine group and the out-of-plane angle of the hydrogen atom, employing MultiReference Configuration-Interaction (MRCI) as well as Complete-Active-Space Self-Consistent-Field (CASSCF) methods. Characteristic properties of the 1pipi*-1pisigma* and 1pisigma*-S0 conical intersections, such as the diabatic-to-adiabatic mixing angle, the geometric phase of the adiabatic electronic wavefunctions, the derivative coupling, as well as adiabatic and diabatic transition dipole moment surfaces were investigated in detail. These data are a prerequisite for future quantum wavepacket simulations of the photodissociation and internal-conversion dynamics of adenine.     

Relaxation pathways of excited N-(triphenylmethyl)salicylidenimine in solutions

Excited state relaxation of N-(triphenylmethyl)-salicylidenimine (MS1) in protic and aprotic solvents has been investigated by means of absorption pump-probe spectroscopy with femtosecond time resolution and fluorescence spectroscopy with picosecond time resolution. Short-lived excited states and long-lived photoproducts have been identified from the differential absorption spectra. Excited states and photoproducts were different under excitation of enol-closed and cis-keto tautomers. As a result, the commonly accepted excited state relaxation model of aromatic anils, which assumes an ultrafast transformation of excited enol-closed tautomers into cis-keto tautomers, has been modified. Performed quantum chemical calculations suggest that hydrogen-bonded ethanol molecules facilitate formation of cis-keto tautomers and are responsible for their different relaxation pathways in comparison with relaxation of excited enol-closed tauromers. Fluorescence decay on a nanosecond time scale was attributed to aggregated MS1 molecules.     

Photodeactivation paths in norbornadiene

The first high level ab initio quantum‐chemical calculations of potential energy surfaces (PESs) for low‐lying singlet excited states of norbornadiene in the gas phase are presented. The optimization of the stationary points (minima and conical intersections) and the recalculation of the energies were performed using the multireference configuration interaction with singles (MR‐CIS) and the multiconfigurational second‐order perturbation (CASPT2) methods, respectively. It was shown that the crossing between valence V₂ and Rydberg R₁ states close to the Franck–Condon (FC) point permits an easy population switch between these states. Also, a new deactivation path in which the doubly excited state with (π₃)² configuration (DE) has a prominent role in photodeactivation from the R₁ state due to the R₁/DE and the DE/V₁ conical intersections very close to the R₁ and DE minima, respectively, was proposed. Subsequent deactivation from the V₁ to the ground state goes through an Olivucci–Robb‐type conical intersection that adopts a rhombic distorted geometry. The deactivation path has negligible barriers, thereby making ultrafast radiationless decay to the ground state possible. © 2013 Wiley Periodicals, Inc.     

Excited electronic states and photophysics of uracil–water complexes

The electronically excited singlet states of complexes of uracil with one water molecule have been studied theoretically using ab initio multireference configuration interaction methods. In agreement with previous theoretical and experimental results, four cyclic isomers of uracil forming hydrogen bonds with the water molecule have been located with energies within 0.2 eV from the lowest energy isomer. Focus has been given on the mechanism for radiationless decay to the ground state after initial UV absorption and on the effect of complexation with water on previously reported radiationless decay pathways. Features on the excited state potential energy surfaces, such as minima, transition states and conical intersections, have been located for all isomers and compared with those of free uracil. The hydrogen-bonded water molecule changes the relative energies of these features and may lead to different excited state dynamics and lifetimes, in agreement with experimental observations.     

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.     

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.