Influence of Selective Carbon 1s Excitation on Auger-Meitner Decay in the ESCA Molecule

Two-dimensional spectral mapping is used to visualize how resonant Auger-Meitner spectra are influenced by the site of the initial core-electron excitation and the symmetry of the core-excited state in the trifluoroethyl acetate molecule (ESCA). We observe a significant enhancement of electron yield for excitation of the COO 1s → π* and CF3 1s → σ* resonances unlike excitation at resonances involving the CH3 and CH2 sites. The CF3 1s → π* and CF3 1s → σ* resonance spectra are very different from each other, with the latter populating most valence states equally. Two complementary electronic structure calculations for the photoelectron cross section and Auger-Meitner intensity are shown to effectively reproduce the site- and state-selective nature of the resonant enhancement features. The site of the core-electron excitation and the respective final state hole locality increase the sensistivity of the photoelectron signal at specific functional group sites. This showcases resonant Auger-Meitner decay as a potentially powerful tool for selectively probing structural changes at specific functional group sites of polyatomic molecules.

Revealing the excited-state dynamics of cytidine and the role of excited-state proton transfer process

The high photostability of DNAs and RNAs is inextricably related to the photochemical and photophysical properties of their building blocks, nucleobases and nucleosides, which can dissipate the absorbed UV light energy in a harmless manner. The deactivation mechanism of the nucleosides, especially the decay pathways of cytidine (Cyd), has been a matter of intense debate. In the current study, we employ high-level electronic structure calculations combined with excited state non-adiabatic dynamic simulations to provide a clear picture of the excited state deactivation of Cyd in both gas phase and aqueous solution. In both environments, a barrierless decay path driven by the ring-puckering motion and a relaxation channel with a small energy barrier driven by the elongation motion of CO bond are assigned to <200 fs and sub-picosecond decay time component, respectively. The presence of ribose group has a subtle effect on the dynamic behavior of Cyd in gas phase as the ribose-to-base hydrogen/proton transfer process is energetically inaccessible with a sizable energy barrier of about 1.4 eV. However, this energy barrier is significantly reduced in water, especially when an explicit water molecule is present. Therefore, we argue that the long-lived decay channel found in aqueous solution could be assigned to the Cyd-water intermolecular hydrogen/proton transfer process. The present study postulates a novel scenario toward deep understanding the intrinsic photostability of DNAs and RNAs and provides solid evidence to disclose the long history debate of cytidine excited-state decay mechanism, especially for the assignment of experimentally observed time components.

Galván I, Rivalta I, Aquilante F, Garavelli M, Lindh R, Segarra-Martí J. Characterizing Conical Intersections in DNA/RNA Nucleobases with Multiconfigurational Wave Functions of Varying Active Space Size

We characterize the photochemically relevant conical intersections between the lowest-lying accessible electronic excited states of the different DNA/RNA nucleobases using Cholesky decomposition-based complete active space self-consistent field (CASSCF) algorithms. We benchmark two different basis set contractions and several active spaces for each nucleobase and conical intersection type, measuring for the first time how active space size affects conical intersection topographies in these systems and the potential implications these may have toward their description of photoinduced phenomena. Our results show that conical intersection topographies are highly sensitive to the electron correlation included in the model: by changing the amount (and type) of correlated orbitals, conical intersection topographies vastly change, and the changes observed do not follow any converging pattern toward the topographies obtained with the largest and most correlated active spaces. Comparison across systems shows analogous topographies for almost all intersections mediating population transfer to the dark 1nO/Nπ* states, while no similarities are observed for the "ethylene-like" conical intersection ascribed to mediate the ultrafast decay component to the ground state in all DNA/RNA nucleobases. Basis set size seems to have a minor effect, appearing to be relevant only for purine-based derivatives. We rule out structural changes as a key factor in classifying the different conical intersections, which display almost identical geometries across active space and basis set change, and we highlight instead the importance of correctly describing the electronic states involved at these crossing points. Our work shows that careful active space selection is essential to accurately describe conical intersection topographies and therefore to adequately account for their active role in molecular photochemistry.

Quantum mechanics/molecular mechanics studies on the excited-state decay mechanisms of cytidine aza-analogues: 5-azacytidine and 2'-deoxy-5-azacytidine in aqueous solution

The excited state properties and deactivation pathways of two DNA methylation inhibitors, i.e., 5-azacytidine (5ACyd) and 2'-deoxy-5-azacytidine (5AdCyd) in aqueous solution, are comprehensively explored with the QM(CASPT2//CASSCF)/MM protocol. We systematically map the feasible decay mechanisms based on the obtained excited-state decay paths involving all the identified minimum-energy structures, conical intersections, and crossing points driving the different internal conversion (IC) and intersystem crossing (ISC) routes in and between the 1ππ*, 1nπ*, 3ππ*, 3nπ*, and S0 states. Unlike the 1nπ* state below the 1ππ* state in 5ACyd, deoxyribose group substitution at the N1 position leads to the 1ππ* state becoming the S1 state in 5AdCyd. In 5ACyd and 5AdCyd, the initially populated 1ππ* state mainly deactivates to the S0 state through the direct 1ππ* → S0 IC or mediated by the 1nπ* state. The former nearly barrierless IC channel of 1ππ* → S0 occurs ultrafast via the nearby low-lying 1ππ*/S0 conical intersection. In the latter IC channel of 1ππ* → 1nπ* → S0, the initially photoexcited 1ππ* state first approaches the nearby S2/S1 conical section 1ππ*/1nπ* and then undergoes efficient IC to the 1nπ* state, followed by the further IC to the initial S0 state via the S1/S0 conical intersection 1nπ*/S0. The 1nπ*/S0 conical intersection is estimated to be located 6.0 and 4.9 kcal mol-1 above the 1nπ* state minimum in 5ACyd and 5AdCyd, respectively, at the QM(CASPT2)/MM level. In addition to the efficient singlet-mediated IC channels, the minor ISC routes would populate 1ππ* to T1(ππ*) through 1ππ* → T1 or 1ππ* → 1nπ* → T1. Relatively, the 1ππ* → 1nπ* → T1 route benefits from the spin-orbit coupling (SOC) of 1nπ*/3ππ* of 8.7 cm-1 in 5ACyd and 10.2 cm-1 in 5AdCyd, respectively. Subsequently, the T1 system will approach the nearby T1/S0 crossing point 3ππ*/S0 driving it back to the S0 state. Given the 3ππ*/S0 crossing point located above the T1 minimum and the small T1/S0 SOC, i.e., 8.4 kcal mol-1 and 2.1 cm-1 in 5ACyd and 6.8 kcal mol-1 and 1.9 cm-1 in 5AdCyd, respectively, the slow T1 → S0 would trap the system in the T1 state for a while. The present work could contribute to understanding the mechanistic photophysics and photochemistry of similar aza-nucleosides and their derivatives.

do Casal M, Pinheiro M, Toldo JM, Barbatti M. Modeling the heating and cooling of a chromophore after photoexcitation

The heating of a chromophore due to internal conversion and its cooling down due to energy dissipation to the solvent are crucial phenomena to characterize molecular photoprocesses. In this work, we simulated the ab initio nonadiabatic dynamics of cytosine, a prototypical chromophore undergoing ultrafast internal conversion, in three solvents—argon matrix, benzene, and water—spanning an extensive range of interactions. We implemented an analytical energy-transfer model to analyze these data and extract heating and cooling times. The model accounts for nonadiabatic effects, and excited- and ground-state energy transfer, and can analyze data from any dataset containing kinetic energy as a function of time. Cytosine heats up in the subpicosecond scale and cools down within 25, 4, and 1.3 ps in argon, benzene, and water, respectively. The time constants reveal that a significant fraction of the benzene and water heating occurs while cytosine is still electronically excited.

An analytical energy-transfer model is implemented to obtain a chromophore's heating and cooling times in a given solvent by using quantities available in nonadiabatic dynamics simulations.

Spectral Fingerprint of Excited-State Energy Transfer in Dendrimers through Polarization-Sensitive Transient-Absorption Pump-Probe Signals: On-the-Fly Nonadiabatic Dynamics Simulations

The time-resolved polarization-sensitive transient-absorption (TA) pump-probe (PP) spectra are simulated using on-the-fly surface-hopping nonadiabatic dynamics and the doorway-window representation of nonlinear spectroscopy. A dendrimer model system composed of two linear phenylene ethynylene units (2-ring and 3-ring) is taken as an example. The ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) contributions as well as the total TA PP signals are obtained and carefully analyzed. It is shown that intramolecular excited-state energy transfer from the 2-ring unit to the 3-ring unit can be conveniently identified by employing pump and probe pulses with different polarizations. Our results demonstrate that time-resolved polarization-sensitive TA PP signals provide a powerful tool for the elucidation of excited-state energy-transfer pathways, notably in molecular systems possessing several optically bright nonadiabatically coupled electronic states with different orientations of transition dipole moments.

Fundamental photophysics of isomorphic and expanded fluorescent nucleoside analogues

Fluorescent nucleoside analogues (FNAs) are structurally diverse mimics of the natural essentially non-fluorescent nucleosides which have found numerous applications in probing the structure and dynamics of nucleic acids as well as their interactions with various biomolecules. In order to minimize disturbance in the labelled nucleic acid sequences, the FNA chromophoric groups should resemble the natural nucleobases in size and hydrogen-bonding patterns. Isomorphic and expanded FNAs are the two groups that best meet the criteria of non-perturbing fluorescent labels for DNA and RNA. Significant progress has been made over the past decades in understanding the fundamental photophysics that governs the spectroscopic and environmentally sensitive properties of these FNAs. Herein, we review recent advances in the spectroscopic and computational studies of selected isomorphic and expanded FNAs. We also show how this information can be used as a rational basis to design new FNAs, select appropriate sequences for optimal spectroscopic response and interpret fluorescence data in FNA applications.

Ab Initio Surface-Hopping Simulation of Femtosecond Transient-Absorption Pump-Probe Signals of Nonadiabatic Excited-State Dynamics Using the Doorway-Window Representation

An ab initio theoretical framework for the simulation of femtosecond time-resolved transient absorption (TA) pump-probe (PP) spectra with quasi-classical trajectories is presented. The simulations are based on the classical approximation to the doorway-window (DW) representation of third-order four-wave-mixing signals. The DW formula accounts for the finite duration and spectral shape of the pump and probe pulses. In the classical DW formalism, classical trajectories are stochastically sampled from a positive definite doorway distribution, and the signals are evaluated by averaging over a positive definite window distribution. Nonadiabatic excited-state dynamics is described by a stochastic surface-hopping algorithm. The method has been implemented for the pyrazine molecule with the second-order algebraic-diagrammatic construction (ADC(2)) ab initio electronic-structure method. The methodology is illustrated by ab initio simulations of the ground-state bleach, stimulated emission, and excited-state absorption contributions to the TA PP spectrum of gas-phase pyrazine.

Nonadiabatic Absorption Spectra and Ultrafast Dynamics of DNA and RNA Photoexcited Nucleobases

We have recently proposed a protocol for Quantum Dynamics (QD) calculations, which is based on a parameterisation of Linear Vibronic Coupling (LVC) Hamiltonians with Time Dependent (TD) Density Functional Theory (TD-DFT), and exploits the latest developments in multiconfigurational TD-Hartree methods for an effective wave packet propagation. In this contribution we explore the potentialities of this approach to compute nonadiabatic vibronic spectra and ultrafast dynamics, by applying it to the five nucleobases present in DNA and RNA. For all of them we computed the absorption spectra and the dynamics of ultrafast internal conversion (100 fs timescale), fully coupling the first 2-3 bright states and all the close by dark states, for a total of 6-9 states, and including all the normal coordinates. We adopted two different functionals, CAM-B3LYP and PBE0, and tested the effect of the basis set. Computed spectra are in good agreement with the available experimental data, remarkably improving over pure electronic computations, but also with respect to vibronic spectra obtained neglecting inter-state couplings. Our QD simulations indicate an effective population transfer from the lowest energy bright excited states to the close-lying dark excited states for uracil, thymine and adenine. Dynamics from higher-energy states show an ultrafast depopulation toward the more stable ones. The proposed protocol is sufficiently general and automatic to promise to become useful for widespread applications.

Directing excited state dynamics via chemical substitution: A systematic study of π-donors and π-acceptors at a carbon-carbon double bond

Functional group substituents are a ubiquitous tool in ground-state organic chemistry often employed to fine-tune chemical properties and obtain desired chemical reaction outcomes. Their effect on photoexcited electronic states, however, remains poorly understood. To help build an intuition for these effects, we have studied ethylene, substituted with electron acceptor (cyano) and/or electron donor (methoxy) substituents, both theoretically and experimentally: using ab initio quantum molecular dynamics and time-resolved photoelectron spectroscopy. Our results show the consistent trend that photo-induced ethylenic dynamics is primarily localized to the carbon with the greater electron density. For doubly substituted ethylenes, the trend is additive when both substituents are located on opposite carbons, whereas the methoxy group (in concert with steric effects) dominates when both substituents are located on a single carbon atom. These results point to the development of rules for structure-dynamics correlations; in this case, a novel mechanistic ultrafast photochemistry for conjugated carbon chains employing long-established chemical concepts.

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.

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.

Photochemical pathways in nucleobases measured with an X-ray FEL

The conversion of light energy into other molecular energetic degrees of freedom is often dominated by ultrafast, non-adiabatic processes. Femtosecond spectroscopy with optical pulses has helped in shaping our understanding of crucial processes in molecular energy-conversion. The advent of new, ultrashort and bright X-ray free electron laser sources opens the possibility to use X-ray-typical element and site sensitivity for ultrafast molecular research. We present two types of spectroscopy, ultrafast Auger and ultrafast X-ray absorption spectroscopy, and discuss their sensitivity to molecular processes. While Auger spectroscopy is able to monitor bond distance changes in the vicinity of an X-ray created core hole, near-edge absorption spectroscopy can deliver high-fidelity information on non-adiabatic transitions involving lone-pair orbitals. We demonstrate these features on the example of the UV-excited nucleobase thymine, investigated at the oxygen K-edge. We find a C-O bond elongation in the Auger data in addition to ππ*/ nπ* non-adiabatic transition in X-ray near-edge absorption. We compare the results from both methods and draw a conclusive scenario of non-adiabatic molecular relaxation after UV excitation. This article is part of the theme issue 'Measurement of ultrafast electronic and structural dynamics with X-rays'.

Two-dimensional electronic-vibrational spectroscopic study of conical intersection dynamics: an experimental and electronic structure study

The relaxation from the lowest singlet excited state of the triphenylmethane dyes, crystal violet and malachite green, is studied via two-dimensional electronic-vibrational (2DEV) spectroscopy.

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.

Photophysical properties of pyrrolocytosine, a cytosine fluorescent base analogue

The photophysical behavior of pyrrolocytosine (PC), a fluorescent base analogue of cytosine, has been investigated using theoretical approaches. The similarities between the PC and cytosine structures allow PC to maintain the pseudo-Watson-Crick base-pairing arrangement with guanine. Cytosine, similar to the other natural nucleobases, is practically non-fluorescent, because of ultrafast radiationless decay occurring through conical intersections. PC displays a much higher fluorescence quantum yield than cytosine, making it an effective fluorescent marker to study the structure, function, and dynamics of DNA/RNA complexes. Similar to 2-aminopurine, a constitutional isomer of adenine that base-pairs with thymine, PC's fluorescence is quenched when it is incorporated into a dinucleotide or a trinucleotide. In this work we examine the photophysical properties of isolated PC, microhydrated PC, as well as, complexes where PC is either base-stacked or hydrogen-bonded with guanine. Our results indicate that hydration affects the radiationless decay pathways in PC by destabilizing conical intersections. The calculations of dimers and trimers show that the radiative decay is affected by π stacking, while the presence of charge transfer states between PC and guanine may contribute to radiationless decay.

How nature covers its bases

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

Ultrafast Nuclear Dynamics of Photoexcited Guanosine-5'-Monophosphate in Three Singlet States

We report measurement of resonance Raman (RR) spectra of guanosine-5'-monophosphate (GMP), a DNA nucleotide at excitation wavelengths throughout its ππ* absorption band (B_b) in the 210-230 nm range. From these data, we constructed wavelength-dependent Raman intensity excitation profiles (REPs) for all observed modes. These profiles and the absorption spectrum were then modeled using self-consistent simulations based on the time-dependent wave packet propagation formalism. We inferred the initial structural dynamics of GMP immediately after photoexcitation in terms of dimensionless displacements. The simulations also provide linewidth-broadening parameters that in turn report on the time scale of dynamics. We compared deduced structural changes in the purine ring upon photoabsorption into the B_b state with those deduced for the two lowest lying ππ* (L_a and L_b at 280 and 248 nm, respectively) excited states of GMP. We find that excitation to the L_b state lengthens C₆-N₁ and C₂═N₃ bonds, which lie along the formation coordinate of various oxidative adducts but B_b excitation does not. We also find that photoabsorption by the B_b state weakens the C₈-N₉ bond and thus might assist imidazole ring opening via cleavage of the same bond. Electronic excitation to different ππ* states of the guanine chromophore results in contrasting structural changes; although absorption by the L_a and L_b states induces expansion of pyrimidine and contraction of imidazole rings, excitation results in overall shrinkage of both the rings. Computed absolute changes in internal coordinates imply that photoexcitation to any of the three singlet states of GMP does not lead directly to the formation of a cation radical of guanine.

Nonadiabatic dynamics simulation of keto isocytosine: a comparison of dynamical performance of different electronic-structure methods

Different reaction channels are obtained in the nonadiabatic dynamics simulations of isocytosine at CASSCF and ADC(2) levels.

Challenges in Simulating Light-Induced Processes in DNA

In this contribution, we give a perspective on the main challenges in performing theoretical simulations of photoinduced phenomena within DNA and its molecular building blocks. We distinguish the different tasks that should be involved in the simulation of a complete DNA strand subject to UV irradiation: (i) stationary quantum chemical computations; (ii) the explicit description of the initial excitation of DNA with light; (iii) modeling the nonadiabatic excited state dynamics; (iv) simulation of the detected experimental observable; and (v) the subsequent analysis of the respective results. We succinctly describe the methods that are currently employed in each of these steps. While for each of them, there are different approaches with different degrees of accuracy, no feasible method exists to tackle all problems at once. Depending on the technique or combination of several ones, it can be problematic to describe the stacking of nucleobases, bond breaking and formation, quantum interferences and tunneling or even simply to characterize the involved wavefunctions. It is therefore argued that more method development and/or the combination of different techniques are urgently required. It is essential also to exercise these new developments in further studies on DNA and subsystems thereof, ideally comprising simulations of all of the different components that occur in the corresponding experiments.

Photoelectron spectra of 2-thiouracil, 4-thiouracil, and 2,4-dithiouracil

Ground- and excited-state UV photoelectron spectra of thiouracils (2-thiouracil, 4-thiouracil, and 2,4-dithiouracil) have been simulated using multireference configuration interaction calculations and Dyson norms as a measure for the photoionization intensity. Except for a constant shift, the calculated spectrum of 2-thiouracil agrees very well with experiment, while no experimental spectra are available for the two other compounds. For all three molecules, the photoelectron spectra show distinct bands due to ionization of the sulphur and oxygen lone pairs and the pyrimidine π system. The excited-state photoelectron spectra of 2-thiouracil show bands at much lower energies than in the ground state spectrum, allowing to monitor the excited-state population in time-resolved UV photoelectron spectroscopy experiments. However, the results also reveal that single-photon ionization probe schemes alone will not allow monitoring all photodynamic processes existing in 2-thiouracil. Especially, due to overlapping bands of singlet and triplet states the clear observation of intersystem crossing will be hampered.

Revealing Deactivation Pathways Hidden in Time-Resolved Photoelectron Spectra

Time-resolved photoelectron spectroscopy is commonly employed with the intention to monitor electronic excited-state dynamics occurring in a neutral molecule. With the help of theory, we show that when excited-state processes occur on similar time scales the different relaxation pathways are completely obscured in the total photoionization signal recorded in the experiment. Using non-adiabatic molecular dynamics and Dyson norms, we calculate the photoionization signal of cytosine and disentangle the transient contributions originating from the different deactivation pathways of its tautomers. In the simulations, the total signal from the relevant keto and enol tautomers can be decomposed into contributions either from the neutral electronic state populations or from the distinct mechanistic pathways across the multiple potential surfaces. The lifetimes corresponding to these contributions cannot be extracted from the experiment, thereby illustrating that new experimental setups are necessary to unravel the intricate non-adiabatic pathways occurring in polyatomic molecules after irradiation by light.

The origin of efficient triplet state population in sulfur-substituted nucleobases

Elucidating the photophysical mechanisms in sulfur-substituted nucleobases (thiobases) is essential for designing prospective drugs for photo- and chemotherapeutic applications. Although it has long been established that the phototherapeutic activity of thiobases is intimately linked to efficient intersystem crossing into reactive triplet states, the molecular factors underlying this efficiency are poorly understood. Herein we combine femtosecond transient absorption experiments with quantum chemistry and nonadiabatic dynamics simulations to investigate 2-thiocytosine as a necessary step to unravel the electronic and structural elements that lead to ultrafast and near-unity triplet-state population in thiobases in general. We show that different parts of the potential energy surfaces are stabilized to different extents via thionation, quenching the intrinsic photostability of canonical DNA and RNA nucleobases. These findings satisfactorily explain why thiobases exhibit the fastest intersystem crossing lifetimes measured to date among bio-organic molecules and have near-unity triplet yields, whereas the triplet yields of canonical nucleobases are nearly zero.

Investigation of the mechanisms of photo-induced formation of cyclobutane dimers of cytosine and 2,4-diaminopyrimidine

The mechanisms of the formation of cyclobutane dimers (CBD) of cytosine and 2,4-diaminopyrimidine were studied at the CC2 theoretical level and cc-pVDZ basis functions. Four orientations of the two monomers are explored: cys-syn, cis-anti, trans-syn, and trans-anti. The research revealed that in all cases the cyclobutane structures are formed along the (1)ππ* excited-state reaction paths of the stacked aggregates. We localized the S1/S0 conical intersections mediating those transformations. The results obtained agree well with the previously reported investigations on the cis-syn cyclodimer formations of other pyrimidines.

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.

Singlet-Triplet States Interaction Regions in DNA/RNA Nucleobase Hypersurfaces

The present study provides new insight into the intrinsic mechanisms for the population of the triplet manifold in DNA nucleobases by determining, at the multiconfigurational CASSCF/CASPT2 level, the singlet-triplet states crossing regions and the main decay paths for their lowest singlet and triplet states after near-UV irradiation. The studied singlet-triplet interacting regions are accessible along the minimum energy path of the initially populated singlet bright (1)ππ* state. In particular, all five natural DNA/RNA nucleobases have, at the end of the main minimum energy path and near a conical intersection of the ground and (1)ππ* states, a low-energy, easily accessible, singlet-triplet crossing region directly connecting the lowest singlet and triplet ππ* excited states. Adenine, thymine, and uracil display additional higher-energy crossing regions related to the presence of low-lying singlet and a triplet nπ* state. These funnels are absent in guanine and cytosine, which have the bright (1)ππ* state lower in energy and less accessible nπ* states. Knowledge of the location and accessibility of these regions, in which the singlet-triplet interaction is related to large spin-orbit coupling elements, may help to understand experimental evidence such as the wavelength dependence measured for the triplet formation quantum yield in nucleobases and the prevalence of adenine and thymine components in the phosphorescence spectra of DNA.