Excited-State Structure and Dynamics of Keto–Amino Cytosine: The 1ππ* State Is Nonplanar and Its Radiationless Decay Is Not Ultrafast

Analytic Gradients for Equation-of-Motion Coupled Cluster with Single, Double, and Perturbative Triple Excitations

Understanding the process of molecular photoexcitation is crucial in various fields, including drug development, materials science, photovoltaics, and more. The electronic vertical excitation energy is a critical property, for example in determining the singlet-triplet gap of chromophores. However, a full understanding of excited-state processes requires additional explorations of the excited-state potential energy surface and electronic properties, which is greatly aided by the availability of analytic energy gradients. Owing to its robust high accuracy over a wide range of chemical problems, equation-of-motion coupled cluster with single and double excitations (EOM-CCSD) is a powerful method for predicting excited-state properties, and the implementation of analytic gradients of many EOM-CCSD variants (excitation energies, ionization potentials, electron attachment energies, etc.) along with numerous successful applications highlights the flexibility of the method. In specific cases where a higher level of accuracy is needed or in more complex electronic structures, the inclusion of triple excitations becomes essential, for example, in the EOM-CCSD* approach of Saeh and Stanton. In this work, we derive and implement for the first time the analytic gradients of EOMEE-CCSD*, which also provides a template for analytic gradients of related excited-state methods with perturbative triple excitations. The capabilities of analytic EOMEE-CCSD* gradients are illustrated by several representative examples.

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.

Direct Observation of Excitation Wavelength-Dependent Ultrafast Intersystem Crossing in Cytosine Nucleoside Solution

A triplet excited state can lead to different DNA photolesions, especially in cytosine and its nucleoside/nucleotide as they are hotspots for DNA mutations. However, the triplet state generation mechanism is in controversy, and experimental evidence of ultrafast intersystem crossing (ISC) has not been registered in these molecules. In this work, ultrafast ISC is directly observed in 2'-deoxycytidine (dCyd) solution by using femtosecond transient absorption spectroscopy. Surprisingly, we demonstrate that ISC in dCyd is sensitive to the excitation wavelength, and a spin-vibronic ISC mechanism is proposed. This finding is the last piece of the dCyd excited-state deactivation mechanism puzzle and sets the base for further investigation of triplet state-involved photophysics and photochemistry in dCyd-containing DNA.

Ultrafast Dynamics of the Two Lowest Bright Excited States of Cytosine and 1-Methylcytosine: A Quantum Dynamical Study

The nonadiabatic quantum dynamics (QD) of cytosine and 1-methylcytosine in the gas phase is simulated for 250 fs after a photoexcitation to one of the first two bright states. The nuclear wavepacket is propagated on the coupled diabatic potential energy surfaces of the lowest seven excited states, including ππ*, nπ*, and Rydberg states along all the vibrational degrees of freedom. We focus in particular on the interplay between the bright and the dark nπ* states, not considering the decay to the ground electronic state. To run these simulations, we implemented an automatic general procedure to parametrize linear vibronic coupling (LVC) models with time-dependent density functional theory (DFT) computations and interfaced it with Gaussian package. The wavepacket was propagated with the multilayer version of the multiconfigurational time dependent Hartree method. Two different density functionals, PBE0 and CAM-B3LYP, which provide a different description of the relative stability of the lowest energy dark states, were used to parametrize the LVC Hamiltonian. Part of the photoexcited population on lowest HOMO-LUMO transition (π_Hπ_L^*) decays within less than 100 fs to a nπ* state which mainly involves a promotion of an electron from the oxygen lone pair to the LUMO (n_Oπ_L^*). The population of the second ππ* state decays almost completely, in <100 fs, not only to π_Hπ_L^* and to n_Oπ_L^* states but also to another nπ_L^* state involving the nitrogen lone pair. The efficiency of the adopted protocol allowed us to check the accuracy of the predictions by repeating the QD simulations with different LVC Hamiltonians parametrized either at the ground-state minimum or at stationary structures of different relevant excited states.

Locating Cytosine Conical Intersections by Laser Experiments and Ab Initio Calculations

The decay mechanism of S₀ → S₁ excited cytosine (Cyt) and the effect of substitution are studied combining jet-cooled spectroscopy (nanosecond resonant two-photon ionization (R2PI) and picosecond lifetime measurements) with CASPT2//CASSCF computations for eight derivatives. For Cyt and five derivatives substituted at N1, C5, and C6, rapid internal conversion sets in at 250-1200 cm^-1 above the 0₀⁰ bands. The break-off in the spectra correlates with the calculated barriers toward the "C5-C6 twist" conical intersection, which unambiguously establishes the decay mechanism at low S₁ state vibrational energies. The barriers increase with substituents that stabilize the charge shifts at C4, C5, and C6 following (¹ππ*) excitation. The R2PI spectra of the clamped derivatives 5,6-trimethyleneCyt (TMCyt) and 1-methyl-TMCyt (1M-TMCyt), which decay along an N3 out-of-plane coordinate, extend up to +3500 and +4500 cm^-1.

UV absorption spectra of DNA bases in the 350-190 nm range: assignment and state specific analysis of solvation effects

The theoretical assignment of electronic spectra of polyatomic molecules is a challenging problem that requires the specification of the character of a large number of electronic states. We propose a procedure for automatically determining the character of electronic transitions and apply it to the study of UV spectra of DNA bases in the gas phase and in the aqueous environment. The procedure is based on the computation of electronic wave function overlaps and accounts for an extensive sampling of nuclear geometries. Novelties of this work are the theoretical assignment of the electronic spectra of DNA bases up to 190 nm and a state specific analysis of solvation effects. By accounting for different effects contributing to the total solvent shift we obtained a good agreement between the computed and experimental spectra. Effects of vibrational averaging, temperature and solvent-induced structural changes shift excitation energies to lower values. Solvent-solute electrostatic interactions are state specific and strongly destabilize nRyd states, and to lesser extent nπ* and πRyd states. Altogether, this results in the stabilization of ππ* states and destabilization of nπ*, πRyd and nRyd states in solution.

UV-promoted radical formation, and near-IR-induced and spontaneous conformational isomerization in monomeric 9-methylguanine isolated in low-temperature Ar matrices

Three low-energy isomers of 9-methylguanine, the amino-oxo (AO) form and two amino-hydroxy (AH1 and AH2) conformers, were trapped from the gas phase into low-temperature argon matrices. The AH1 and AH2 isomers, differing in the orientation of the OH group, were found to transform into each other upon excitation with near-IR light. The population of the AO form of the compound was not changed upon any near-IR irradiation of the matrix samples. Using monochromatic near-IR light, generated by a frequency-tunable laser source, it was possible to selectively induce the AH1 → AH2 or AH2 → AH1 conversion. Photoreversibility of this conformational transformation was then demonstrated. Exposure of matrix-isolated monomers of 9-methylguanine to broadband near-IR light also led to conformational conversions within the amino-hydroxy tautomeric form; the final stage of this process was always the same photostationary state independent of the initial ratio of AH1 and AH2 populations. Spontaneous conformational conversion, transforming the higher-energy AH2 form into the lower-energy AH1 isomer, was observed for matrix-isolated monomers of 9-methylguanine kept in the dark. The mechanism of this process must rely on quantum tunneling of the light hydrogen atom. Irradiation of matrix-isolated 9-methylguanine with UV laser light at λ = 288 or 285 nm led to a substantial consumption of the two AH forms, while the amount of AO isomer remained unchanged. On the other hand, a decrease in the population of the AO isomer occurred upon excitations at shorter wavelengths, λ = 280 or 275 nm. The spectral changes observed after UV-irradiation suggest the generation (and stabilization in the matrix) of a radical species, resulting from the photocleavage of the O-H or N1-H bonds, in the AH or AO isomer, respectively.

Excited-state vibrations, lifetimes, and nonradiative dynamics of jet-cooled 1-ethylcytosine

The S₁ excited-state lifetime of jet-cooled 1-ethylcytosine (1ECyt) is ∼1 ns, one of the longest lifetimes for cytosine derivatives to date. Here, we analyze its S₀ → S₁ vibronic spectrum using two-color resonant two-photon ionization and UV/UV holeburning spectroscopy. Compared to cytosine and 1-methylcytosine, the S₀ → S₁ spectrum of 1ECyt shows a progression in the out-of-plane "butterfly" mode ν₁ ^', identified by spin-component scaled-second-order coupled-cluster method ab initio calculations. We also report time-resolved S₁ state nonradiative dynamics at ∼20 ps resolution by the pump/delayed ionization technique. The S₁ lifetime increases with the number of ν₁ ^' quanta from τ = 930 ps at v₁ ^'=0 to 1030 ps at v₁ ^'=2, decreasing to 14 ps at 710 cm^-1 vibrational energy. We measured the rate constants for S₁ ⇝ S₀ internal conversion and S₁ ⇝ T₁ intersystem crossing (ISC): At the v' = 0 level, k_IC is 8 × 10⁸ s^-1 or three times smaller than 1-methylcytosine. The ISC rate constant from v' = 0 to the T₁(³ππ^*) state is k_ISC = 2.4 × 10⁸ s^-1, 10 times smaller than the ISC rate constants of cytosine, but similar to that of 1-methylcytosine. Based on the calculated S₁(¹ππ^*) state radiative lifetime τ_rad = 12 ns, the fluorescence quantum yield of 1ECyt is Φ_fl ∼ 7% and the intersystem crossing yield is Φ_ISC ∼ 20%. We measured the adiabatic ionization energy of 1-ethylcytosine via excitation of the S₁ state as 8.353 ± 0.008 eV, which is 0.38 eV lower than that of amino-keto cytosine. Measurement of the ionization energy of the long-lived T₁(ππ^*) state formed via ISC reveals that it lies 3.2-3.4 eV above the S₀ ground state.

Ultrafast Excited-State Dynamics of Hydrogen-Bonded Cytosine Microsolvated Clusters with Protic and Aprotic Polar Solvents

Microsolvation effects on the ultrafast excited-state deactivation dynamics of cytosine (Cy) were studied in hydrogen-bonded Cy clusters with protic and aprotic solvents using mass-resolved femtosecond pump-probe ionization spectroscopy. Two protic solvents, water (H₂O) and methanol (MeOH), and one aprotic solvent, tetrahydrofuran (THF), were investigated, and transients of Cy·(H₂O)_1-6, Cy·(MeOH)_1-3, and Cy·THF microsolvated clusters produced in supersonic expansions were measured. With the aid of electronic structure calculations, we assigned the observed dynamics to the low-energy isomers of various Cy clusters and discussed the microsolvation effect on the excited-state deactivation dynamics. With the protic solvents only the microsolvated clusters of Cy keto tautomer were observed. The observed decay time constants of Cy·(H₂O) _n are 0.5 ps for n = 1 and ∼0.2-0.25 ps for n = 2-6. For Cy·(MeOH) _n clusters, the decay time constant for n = 1 cluster is similar to that of the Cy monohydrate, but for n = 2 and 3 the decays are about a factor of 2 slower than the corresponding microhydrates. With the aprotic solvent, THF, hydrogen-bonded complexes of both keto and enol tautomers are present in the beam. The keto-Cy·THF shows a decay similar to that of the keto-Cy monomer, whereas the enol-Cy·THF exhibits a 2-fold slower decay than the enol-Cy monomer, suggesting an increase in the barrier to excited-state deactivation upon binding of one THF molecule to the enol form of Cy.

Probing the excited state relaxation dynamics of pyrimidine nucleosides in chloroform solution

Ultrafast transient electronic and vibrational absorption spectroscopy (TEAS and TVAS) of 2'-deoxy-cytidine (dC) and 2'-deoxy-thymidine (dT) dissolved in chloroform examines their excited-state dynamics and the recovery of ground electronic state molecules following absorption of ultraviolet light. The chloroform serves as a weakly interacting solvent, allowing comparisons to be drawn with prior experimental studies of the photodynamics of these nucleosides in the gas phase and in polar solvents such as water. The pyrimidine base nucleosides have some propensity to dimerize in aprotic solvents, but the monomer photochemistry can be resolved clearly and is the focus of this study. UV absorption at a wavelength of 260 nm excites a ¹ππ* ← S₀ transition, but prompt crossing of a significant fraction (50% in dC, 17% in dT) of the ¹ππ* population into a nearby ¹nπ* state is too fast for the experiments to resolve. The remaining flux on the ¹ππ* state leaves the vertical Franck-Condon region and encounters a conical intersection with the ground electronic state of ethylenic twist character. In dC, the ¹ππ* state decays to the ground state with a time constant of 1.1 ± 0.1 ps. The lifetime of the ¹nπ* state is much longer in the canonical forms of both molecules: recovery of the ground state population from these states occurs with time constants of 18.6 ± 1.1 ps in amino-oxo dC and ∼114 ps in dT, indicating potential energy barriers to the ¹nπ*/S₀ conical intersections. The small fraction of the imino-oxo tautomer of dC present in solution has a longer-lived ¹nπ* state with a lifetime for ground state recovery of 193 ± 55 ps. No evidence is found for photo-induced tautomerization of amino-oxo dC to the imino-oxo form, or for population of low lying triplet states of this nucleoside. In contrast, ∼8% of the UV-excited dT molecules access the long-lived T₁ (³ππ*) state through the ¹nπ* state. The primary influence of the solvent appears to be the degree to which it destabilizes the states of ¹nπ* character, with consequences for the lifetimes of these states as well as the triplet state yields.

Planarizing cytosine: The S1 state structure, vibrations, and nonradiative dynamics of jet-cooled 5,6-trimethylenecytosine

We measure the S₀ → S₁ spectrum and time-resolved S₁ state nonradiative dynamics of the "clamped" cytosine derivative 5,6-trimethylenecytosine (TMCyt) in a supersonic jet, using two-color resonant two-photon ionization (R2PI), UV/UV holeburning, and ns time-resolved pump/delayed ionization. The experiments are complemented with spin-component scaled second-order approximate coupled cluster (SCS-CC2), time-dependent density functional theory, and multi-state second-order perturbation-theory (MS-CASPT2) ab initio calculations. While the R2PI spectrum of cytosine breaks off ∼500 cm^-1 above its 0₀⁰ band, that of TMCyt extends up to +4400 cm^-1 higher, with over a hundred resolved vibronic bands. Thus, clamping the cytosine C⁵-C⁶ bond allows us to explore the S₁ state vibrations and S₀ → S₁ geometry changes in detail. The TMCyt S₁ state out-of-plane vibrations ν₁^', ν₃^', and ν₅^' lie below 420 cm^-1, and the in-plane ν₁₁^', ν₁₂^', and ν₂₃^' vibrational fundamentals appear at 450, 470, and 944 cm^-1. S₀  →  S₁ vibronic simulations based on SCS-CC2 calculations agree well with experiment if the calculated ν₁^', ν₃^', and ν₅^' frequencies are reduced by a factor of 2-3. MS-CASPT2 calculations predict that the ethylene-type S₁ ⇝ S₀ conical intersection (CI) increases from +366 cm^-1 in cytosine to >6000 cm^-1 in TMCyt, explaining the long lifetime and extended S₀ → S₁ spectrum. The lowest-energy S₁ ⇝ S₀ CI of TMCyt is the "amino out-of-plane" (OP_X) intersection, calculated at +4190 cm^-1. The experimental S₁ ⇝ S₀ internal conversion rate constant at the S₁(v^'=0) level is k_IC=0.98-2.2⋅10⁸ s^-1, which is ∼10 times smaller than in 1-methylcytosine and cytosine. The S₁(v^'=0) level relaxes into the T₁(³ππ*) state by intersystem crossing with k_ISC=0.41-1.6⋅10⁸ s^-1. The T₁ state energy is measured to lie 24 580±560 cm^-1 above the S₀ state. The S₁(v^'=0) lifetime is τ=2.9 ns, resulting in an estimated fluorescence quantum yield of Φ_fl=24%. Intense two-color R2PI spectra of the TMCyt amino-enol tautomers appear above 36 000 cm^-1. A sharp S₁ ionization threshold is observed for amino-keto TMCyt, yielding an adiabatic ionization energy of 8.114±0.002 eV.

A 'bottom up', ab initio computational approach to understanding fundamental photophysical processes in nitrogen containing heterocycles, DNA bases and base pairs

The availability of non-radiative decay mechanisms by which photoexcited molecules can revert to their ground electronic state, without experiencing potentially deleterious chemical transformation, is fundamental to molecular photostability. This Perspective Article combines results of new ab initio electronic structure calculations and prior experimental data in an effort to systematise trends in the non-radiative decay following UV excitation of selected families of heterocyclic molecules. We start with the prototypical uni- and bicyclic molecules phenol and indole, and explore the structural and photophysical consequences of incorporating progressively more nitrogen atoms within the respective ring structures en route to the DNA bases thymine, cytosine, adenine and guanine. For each of the latter, we identify low energy non-radiative decay pathways via conical intersections with the ground state potential energy surface accessed by out-of-plane ring deformations. This is followed by summary descriptions and illustrations of selected rival (electron driven H atom transfer) non-radiative excited state decay processes that demand consideration once the nucleobases are merely components in larger biomolecular systems like nucleosides, and both individual and stacked base-pairs.

How nature covers its bases

The response of nucleobases to UV radiation depends on structure in subtle ways, as revealed by gas-phase experiments.

Excited-State Dynamics of Isocytosine: A Hybrid Case of Canonical Nucleobase Photodynamics

We present resonant two-photon ionization (R2PI) spectra of isocytosine (isoC) and pump-probe results on two of its tautomers. IsoC is one of a handful of alternative bases that have been proposed in scenarios of prebiotic chemistry. It is structurally similar to both cytosine (C) and guanine (G). We compare the excited-state dynamics with the Watson-Crick (WC) C and G tautomeric forms. These results suggest that the excited-state dynamics of WC form of G may primarily depend on the heterocyclic substructure of the pyrimidine moiety, which is chemically identical to isoC. For WC isoC we find a single excited-state decay with a rate of ∼10¹⁰ s^-1, while the enol form has multiple decay rates, the fastest of which is 7 times slower than for WC isoC. The excited-state dynamics of isoC exhibits striking similarities with that of G, more so than with the photodynamics of C.

Simulation of Electronic Circular Dichroism of Nucleic Acids: From the Structure to the Spectrum

We present a quantum mechanical (QM) simulation of the electronic circular dichroism (ECD) of nucleic acids (NAs). The simulation combines classical molecular dynamics, to obtain the structure and its temperature-dependent fluctuations, with a QM excitonic model to determine the ECD. The excitonic model takes into account environmental effects through a polarizable embedding and uses a refined approach to calculate the electronic couplings in terms of full transition densities. Three NAs with either similar conformations but different base sequences or similar base sequences but different conformations have been investigated and the results were compared with experimental observations; a good agreement was seen in all cases. A detailed analysis of the nature of the ECD bands in terms of their excitonic composition was also carried out. Finally, a comparison between the QM and the DeVoe models clearly revealed the importance of including fluctuations of the excitonic parameters and of accurately determining the electronic couplings. This study demonstrates the feasibility of the ab initio simulation of the ECD spectra of NAs, that is, without the need of experimental structural or electronic data.

Life in the light: nucleic acid photoproperties as a legacy of chemical evolution

Photophysical investigations of the canonical nucleobases that make up DNA and RNA during the past 15 years have revealed that excited states formed by the absorption of UV radiation decay with subpicosecond lifetimes (i.e., <10(-12) s). Ultrashort lifetimes are a general property of absorbing sunscreen molecules, suggesting that the nucleobases are molecular survivors of a harsh UV environment. Encoding the genome using photostable building blocks is an elegant solution to the threat of photochemical damage. Ultrafast excited-state deactivation strongly supports the hypothesis that UV radiation played a major role in shaping molecular inventories on the early Earth before the emergence of life and the subsequent development of a protective ozone shield. Here, we review the general physical and chemical principles that underlie the photostability, or "UV hardiness", of modern nucleic acids and discuss the possible implications of these findings for prebiotic chemical evolution. In RNA and DNA strands, much longer-lived excited states are observed, which at first glance appear to increase the risk of photochemistry. It is proposed that the dramatically different photoproperties that emerge from assemblies of photostable building blocks may explain the transition from a world of molecular survival to a world in which energy-rich excited electronic states were eventually tamed for biological purposes such as energy transduction, signaling, and repair of the genetic machinery.

Non-radiative processes in protonated diazines, pyrimidine bases and an aromatic azine

The excited state lifetimes of DNA bases are often very short due to very efficient non-radiative processes assigned to the ππ*-nπ* coupling. A set of protonated aromatic diazine molecules (pyridazine, pyrimidine and pyrazine C4H5N2(+)) and protonated pyrimidine DNA bases (cytosine, uracil and thymine), as well as the protonated pyridine (C5H6N(+)), have been investigated. For all these molecules except one tautomer of protonated uracil (enol-keto), electronic spectroscopy exhibits vibrational line broadening. Excited state geometry optimization at the CC2 level has been conducted to find out whether the excited state lifetimes measured from line broadening can be correlated to the calculated ordering of the ππ* and nπ* states and the ππ*-nπ* energy gap. The short lifetimes, observed when one nitrogen atom of the ring is not protonated, can be rationalized by relaxation of the ππ* state to the nπ* state or directly to the electronic ground state through ring puckering.

Intersystem crossing rates of S1 state keto-amino cytosine at low excess energy

The amino-keto tautomer of supersonic jet-cooled cytosine undergoes intersystem crossing (ISC) from the v = 0 and low-lying vibronic levels of its S1((1)ππ(∗)) state. We investigate these ISC rates experimentally and theoretically as a function of S1 state vibrational excess energy Eexc. The S1 vibronic levels are pumped with a ∼5 ns UV laser, the S1 and triplet state ion signals are separated by prompt or delayed ionization with a second UV laser pulse. After correcting the raw ISC yields for the relative S1 and T1 ionization cross sections, we obtain energy dependent ISC quantum yields QISC (corr)=1%-5%. These are combined with previously measured vibronic state-specific decay rates, giving ISC rates kISC = 0.4-1.5 ⋅ 10(9) s(-1), the corresponding S1⇝S0 internal conversion (IC) rates are 30-100 times larger. Theoretical ISC rates are computed using SCS-CC2 methods, which predict rapid ISC from the S1; v = 0 state with kISC = 3 ⋅ 10(9) s(-1) to the T1((3)ππ(∗)) triplet state. The surprisingly high rate of this El Sayed-forbidden transition is caused by a substantial admixture of (1)nOπ(∗) character into the S1((1)ππ(∗)) wave function at its non-planar minimum geometry. The combination of experiment and theory implies that (1) below Eexc = 550 cm(-1) in the S1 state, S1⇝S0 internal conversion dominates the nonradiative decay with kIC ≥ 2 ⋅ 10(10) s(-1), (2) the calculated S1⇝T1 ((1)ππ(∗)⇝(3)ππ(∗)) ISC rate is in good agreement with experiment, (3) being El-Sayed forbidden, the S1⇝T1 ISC is moderately fast (kISC = 3 ⋅ 10(9) s(-1)), and not ultrafast, as claimed by other calculations, and (4) at Eexc ∼ 550 cm(-1) the IC rate increases by ∼50 times, probably by accessing the lowest conical intersection (the C5-twist CI) and thereby effectively switching off the ISC decay channels.

Remarkable effects of solvent and substitution on the photo-dynamics of cytosine: a femtosecond broadband time-resolved fluorescence and transient absorption study

Cytosine (Cyt) among all the nucleic acid bases features the most complex and least understood nonradiative deactivation, a process that is crucially important for its photostability. Herein, the excited state dynamics of Cyt and a series of its N1- and C5-derivatives, including the full set of Cyt nucleosides and nucleotides in DNA and RNA and the nucleosides of 5-methyl cytosine, 5-methylcytidine and 2'-deoxy-5-methylcytidine, have been investigated in water and in methanol employing femtosecond broadband time-resolved fluorescence coupled with fs transient absorption spectroscopy. The results reveal remarkable state-specific effects of the substitution and solvent in tuning distinctively the timescales and pathways of the nonradiative decays. For Cyt and the N1-derivatives, the nonradiative deactivations occur in a common two-state process through three channels, two from the light-absorbing ππ* state with respectively the sub-picosecond (∼0.2 ps) and the picosecond (∼1.5 ps) time constant, and the third is due to an optically dark nπ* state with the lifetime ranging from several to hundreds of picoseconds depending on solvents and substitutions. Compared to Cyt, the presence of the ribose or deoxyribose moiety at the N1 position of N1-derivatives facilitates the formation of the nπ* at the sub-picosecond timescale and at the same time increases its lifetime by ∼4-6 times in both water and methanol. In sharp contrast, the existence of the methyl group at the C5 position of the C5-derivatives eliminates completely the sub-picosecond ππ* channel and the channel due to the nπ*, but on the other hand slows down the decay of the ππ* state which after relaxation exhibits a single time constant of ∼4.1 to ∼7.6 ps depending on solvents. Varying the solvent from water to methanol accelerates only slightly the decay of the ππ* state in all the compounds; while for Cyt and its N1-derivatives, this change of solvent also retards strongly the nπ* channel, prolongs its lifetime from such as ∼7.7 ps in water to ∼52 ps in methanol for Cyt and from ∼30 ps in water to ∼186 ps in methanol for deoxycytidine. The spectral signatures we obtained for the ππ* and the nπ* states allow unambiguous evidence for clarifying uncertainties in the excited states of Cyt and the derivatives. The results provide a unifying experimental characterization at an improved level of detail about the photophysics of Cyt and its analogues under biologically relevant conditions and may help in understanding the photostability as well as photo-damages of the bases and related DNAs.

Gas-Phase Cytosine and Cytosine-N1-Derivatives Have 0.1-1 ns Lifetimes Near the S1 State Minimum

Ultraviolet radiative damage to DNA is inefficient because of the ultrafast S1 ⇝ S0 internal conversion of its nucleobases. Using picosecond pump-ionization delay measurements, we find that the S1((1)ππ*) state vibrationless lifetime of gas-phase keto-amino cytosine (Cyt) is τ = 730 ps or ∼ 700 times longer than that measured by femtosecond pump-probe ionization at higher vibrational excess energy, Eexc. N1-Alkylation increases the S1 lifetime up to τ = 1030 ps for N1-ethyl-Cyt but decreases it to 100 ps for N1-isopropyl-Cyt. Increasing the vibrational energy to Eexc = 300-550 cm(-1) decreases the lifetimes to 20-30 ps. The nonradiative dynamics of S1 cytosine is not solely a property of the amino-pyrimidinone chromophore but is strongly influenced by the N1-substituent. Correlated excited-state calculations predict that the gap between the S2((1)nOπ*) and S1((1)ππ*) states decreases along the series of N1-derivatives, thereby influencing the S1 state lifetime.

Quantum Mechanical Studies on the Photophysics and the Photochemistry of Nucleic Acids and Nucleobases

The photophysics and photochemistry of DNA is of great importance due to the potential damage of the genetic code by UV light. Quantum mechanical studies have played a key role in interpretating the results of modern time-resolved pump-probe spectroscopy, and in elucidating the main photoactivated reactive paths. This review provides a concise, complete picture of the computational studies carried out, approximately, in the past decade. We start with an overview of the photophysics of the nucleobases in the gas phase and in solution. We discuss the proposed mechanisms for ultrafast decay to the ground state, that involve conical intersections, consider the role of triplet states, and analyze how the solvent modulates the photophysics. Then we move to larger systems, from dinucleotides to single- and double-stranded oligonucleotides. We focus on the possible role of charge transfer and delocalized or excitonic states in the photophysics of these systems and discuss the main photochemical paths. We finish with an outlook on the current challenges in the field and future directions of research.

Time-resolved photoelectron imaging of the isolated deprotonated nucleotides

Time-resolved photoelectron spectroscopy of deprotonated nucleotides provides new insights into their relaxation dynamics.

Excited-state tautomerization of gas-phase cytosine

In order to investigate experimentally observed phototautomerization of gas-phase cytosine, several excited-state tautomerization mechanisms were characterized at the EOM-CCSD and TDDFT levels. All pathways that took place exclusively on the S1 surface were found to have significant barriers that were much higher than the barriers involved in radiationless decay of cytosine tautomers through conical intersections back to the ground state; tautomerization in this fashion cannot compete with radiationless relaxation. However, an alternative possibility is that the conical intersections that facilitate radiationless decay could also facilitate tautomerization. Barrierless pathways indicate that it is energetically possible that bifurcation at the conical intersections can lead to a subset of the population reaching different tautomers. This could be an explanation for the observed tautomerization of keto cytosine after exposure to low-energy UV light.