Excitation of Nucleobases from a Computational Perspective II: Dynamics

Time-resolved x-ray spectroscopy of nucleobases and their thionated analogs

The photoinduced relaxation dynamics of nucleobases and their thionated analogs have been investigated extensively over the past decades motivated by their crucial role in organisms and their application in medical and biochemical research and treatment. Most of these studies focused on the spectroscopy of valence electrons and fragmentation. The advent of ultrashort x-ray laser sources such as free-electron lasers, however, opens new opportunities for studying the ultrafast molecular relaxation dynamics utilizing the site- and element-selectivity of x-rays. In this review, we want to summarize ultrafast experiments on thymine and 2-thiouracil performed at free-electron lasers. We performed time-resolved x-ray absorption spectroscopy at the oxygen K-edge after UV excitation of thymine. In addition, we investigated the excited state dynamics of 2-tUra via x-ray photoelectron spectroscopy at sulfur. For these methods, we show a strong sensitivity to the electronic state or charge distribution, respectively. We also performed time-resolved Auger-Meitner spectroscopy, which shows spectral shifts associated with internuclear distances close to the probed site. We discuss the complementary aspects of time-resolved x-ray spectroscopy techniques compared to optical and UV spectroscopy for the investigation of ultrafast relaxation processes.

Parametrically Managed Activation Functions for Improved Global Potential Energy Surfaces for Six Coupled 5A' States and Fourteen Coupled 3A' States of O + O2

We report new potential energy surfaces for six coupled 5A' states and 14 coupled 3A' states of O3. The new surfaces are created by parametrically managed diabatization by deep neural network (PM-DDNN). The PM-DDNN method uses calculated adiabatic potential energy surfaces to discover and fit an underlying adiabatic-equivalent set of diabatic surfaces and their couplings and obtains the fit to the adiabatic surfaces by diagonalization of the diabatic potential energy matrix (DPEM). The procedure yields the adiabatic surfaces and their gradients, as well as the DPEM and its gradient. If desired one can also compute the nonadiabatic coupling due to the transformation. The present work improves on previous work by using a new coordinate to guide the decay of the neural network contribution to the many-body fit to the whole DPEM. The main objective was to obtain smoother potentials than the previous ones with better suitability for dynamics calculations, and this was achieved. Furthermore, we obtained suitably small deviations from the input reference data. For the six coupled 5A' surfaces, the 60,366 data below 10 eV are fit with a mean unsigned error (MUE) of 49 meV, and for the 14 coupled 3A' surfaces, the 76,733 data below 10 eV are fit with an MUE of 28 meV. The data below 5 eV fit even more accurately with MUEs of 37 meV (5A') and 20 meV (3A').

Ultrafast Intersystem Crossing Dynamics of 6-Selenoguanine in Water

Rationalizing the photochemistry of nucleobases where an oxygen is replaced by a heavier atom is essential for applications that exploit near-unity triplet quantum yields. Herein, we report on the ultrafast excited-state deactivation mechanism of 6-selenoguanine (6SeGua) in water by combining nonadiabatic trajectory surface-hopping dynamics with an electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. We find that the predominant relaxation mechanism after irradiation starts on the bright singlet S₂ state that converts internally to the dark S₁ state, from which the population is transferred to the triplet T₂ state via intersystem crossing and finally to the lowest T₁ state. This S₂ → S₁ → T₂ → T₁ deactivation pathway is similar to that observed for the lighter 6-thioguanine (6tGua) analogue, but counterintuitively, the T₁ lifetime of the heavier 6SeGua is shorter than that of 6tGua. This fact is explained by the smaller activation barrier to reach the T₁/S₀ crossing point and the larger spin-orbit couplings of 6SeGua compared to 6tGua. From the dynamical simulations, we also calculate transient absorption spectra (TAS), which provide two time constants (τ₁ = 131 fs and τ₂ = 191 fs) that are in excellent agreement with the experimentally reported value (τ_exp = 130 ± 50 fs) (Farrel et al. J. Am. Chem. Soc. 2018, 140, 11214). Intersystem crossing itself is calculated to occur with a time scale of 452 ± 38 fs, highlighting that the TAS is the result of a complex average of signals coming from different nonradiative processes and not intersystem crossing alone.

Excited state absorption of DNA bases in the gas phase and in chloroform solution: a comparative quantum mechanical study

We study the excited state absorption (ESA) properties of the four DNA bases (thymine, cytosine, adenine, and guanine) by different single reference quantum mechanical methods, namely, equation of motion coupled cluster singles and doubles (EOM-CCSD), singles, doubles and perturbative triples (EOM-CC3), and time-dependent density functional theory (TD-DFT), with the long-range corrected CAM-B3LYP functional. Preliminary results at the Tamm-Dancoff (TDA) CAM-B3LYP level using the maximum overlap method (MOM) are reported for thymine. In the gas phase, the three methods predict similar One Photon Absorption (OPA) spectra, which are consistent with the experimental results and with the most accurate computational studies available in the literature. The ESA spectra are then computed for the ππ* states (one for pyrimidine, two for purines) associated with the lowest-energy absorption band, and for the close-lying nπ* state. The EOM-CC3, EOM-CCSD and CAM-B3LYP methods provide similar ESA spectral patterns, which are also in qualitative agreement with literature RASPT2 results. Once validated in the gas phase, TD-CAM-B3LYP has been used to compute the ESA in chloroform, including solvent effects by the polarizable continuum model (PCM). The predicted OPA and ESA spectra in chloroform are very similar to those in the gas phase, most of the bands shifting by less than 0.1 eV, with a small increase of the intensities and a moderate destabilization of the nπ* state. Finally, ESA spectra have been computed from the minima of the lowest energy ππ* state, and found in line with the available experimental transient absorption spectra of the nucleosides in solution, providing further validation of our computational approach.

The Excited State Dynamics of a Mutagenic Cytidine Etheno Adduct Investigated by Combining Time-Resolved Spectroscopy and Quantum Mechanical Calculations

Joint femtosecond fluorescence upconversion experiments and theoretical calculations provide a hitherto unattained degree of characterization and understanding of the mutagenic etheno adduct 3,N4-etheno-2'-deoxycytidine (εdC) excited state relaxation. This endogenously formed lesion is attracting great interest because of its ubiquity in human tissues and its highly mutagenic properties. The εdC fluorescence is modified with respect to that of the canonical base dC, with a 3-fold increased lifetime and quantum yield at neutral pH. This behavior is amplified upon protonation of the etheno ring (εdCH⁺). Quantum mechanical calculations show that the lowest energy state ππ*1 is responsible for the fluorescence and that the main nonradiative decay pathway to the ground state goes through an ethene-like conical intersection, involving the out-of-plane motion of the C5 and C6 substituents. This conical intersection is lower in energy than the ππ* state (ππ*1) minimum, but a sizable energy barrier explains the increase of εdC and εdCH⁺ fluorescence lifetimes with respect to that of dC.

The Ultrafast Quantum Dynamics of Photoexcited Adenine-Thymine Basepair Investigated with a Fragment-based Diabatization and a Linear Vibronic Coupling Model

In this contribution we present a quantum dynamical study of the photoexcited hydrogen bonded base pair adenine-thymine (AT) in a Watson-Crick arrangement. To that end, we parametrize Linear Vibronic Coupling (LVC) models with Time-Dependent Density Functional Theory (TD-DFT) calculations, exploiting a fragment diabatization scheme (FrD) we have developed to define diabatic states on the basis of individual chromophores in a multichromophoric system. Wavepacket propagations were run with the multilayer extension of the Multiconfiguration Time-Dependent Hartree method. We considered excitations to the three lowest bright states, a ππ* state of thymine and two ππ* states (La and Lb) of adenine, and we found that on the 100 fs time scale the main decay pathways involve intramonomer population transfers toward nπ* states of the same nucleobase. In AT this transfer is less effective than in the isolated nucleobases, because hydrogen bonding destabilizes the nπ* states. The population transfer to the A → T charge transfer state is negligible, making the ultrafast (femtosecond) decay through the proton coupled electron transfer mechanism unlikely, in line with experimental results in apolar solvents. The excitation energy transfer is also very small. We carefully compare the predictions of LVC Hamiltonians obtained with different sets of diabatic states, defined so to match either local states of the two separated monomers or the base pair adiabatic states in the Franck-Condon region. To that end we also extend the flexibility of the FrD-LVC approach, introducing a new strategy to define fragments diabatic states that account for the effect of the rest of the multichromohoric system through a Molecular Mechanics potential.

Photophysical Deactivation Mechanisms of the Pyrimidine Analogue 1-Cyclohexyluracil

The photophysical relaxation mechanisms of 1-cyclohexyluracil, in vacuum and water, were investigated by employing the Multi-State CASPT2 (MS-CASPT2, Multi-State Complete Active-Space Second-Order Perturbation Theory) quantum chemical method and Dunning's cc-pVDZ basis sets. In both environments, our results suggest that the primary photophysical event is the population of the S11(ππ*) bright state. Afterwards, two likely deactivation pathways can take place, which is sustained by linear interpolation in internal coordinates defined via Z-Matrix scans connecting the most important characteristic points. The first one (Route 1) is the same relaxation mechanism observed for uracil, its canonical analogue, i.e., internal conversion to the ground state through an ethylenic-like conical intersection. The other route (Route 2) is the direct population transfer from the S11(ππ*) bright state to the T23(nπ*) triplet state via an intersystem crossing process involving the (S11(ππ*)/T23(nπ*))STCP singlet-triplet crossing point. As the spin-orbit coupling is not too large in either environment, we propose that most of the electronic population initially on the S11(ππ*) state returns to the ground following the same ultrafast deactivation mechanism observed in uracil (Route 1), while a smaller percentage goes to the triplet manifold. The presence of a minimum on the S11(ππ*) potential energy hypersurface in water can help to understand why experimentally it is noticed suppression of the triplet states population in polar protic solvent.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Photoprotecting Uracil by Coupling with Lossy Nanocavities

We analyze how the photorelaxation dynamics of a molecule can be controlled by modifying its electromagnetic environment using a nanocavity mode. In particular, we consider the photorelaxation of the RNA nucleobase uracil, which is the natural mechanism to prevent photodamage. In our theoretical work, we identify the operative conditions in which strong coupling with the cavity mode can open an efficient photoprotective channel, resulting in a relaxation dynamics twice as fast as the natural one. We rely on a state-of-the-art chemically detailed molecular model and a non-Hermitian Hamiltonian propagation approach to perform full-quantum simulations of the system dissipative dynamics. By focusing on the photon decay, our analysis unveils the active role played by cavity-induced dissipative processes in modifying chemical reaction rates, in the context of molecular polaritonics. Remarkably, we find that the photorelaxation efficiency is maximized when an optimal trade-off between light-matter coupling strength and photon decay rate is satisfied. This result is in contrast with the common intuition that increasing the quality factor of nanocavities and plasmonic devices improves their performance. Finally, we use a detailed model of a metal nanoparticle to show that the speedup of the uracil relaxation could be observed via coupling with a nanosphere pseudomode, without requiring the implementation of complex nanophotonic structures.

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.

First-principles characterization of the singlet excited state manifold in DNA/RNA nucleobases

TD-DFT characterization of the high-energy singlet excited state manifold of the canonical DNA/RNA nucleobasesin vacuumis assessed against RASPT2 reference computations for reliable simulations of linear and non-linear electronic spectra.

Approximate Analytical Gradients and Nonadiabatic Couplings for the State-Average Density Matrix Renormalization Group Self-Consistent-Field Method

We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wave function. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) ansatz, which employs a Lagrangian, and the corresponding Lagrange multipliers are obtained from a solution of the coupled-perturbed CASSCF (CP-CASSCF) equations. We introduce a definition of the matrix product state (MPS) Lagrange multipliers based on a single-site tensor in a mixed-canonical form of the MPS, such that a sweep procedure is avoided in the solution of the CP-CASSCF equations. We apply our implementation to the optimization of a conical intersection in 1,2-dioxetanone, where we are able to fully reproduce the SA-CASSCF result up to arbitrary accuracy.

Initial excited-state structural dynamics of 2′-deoxyadenosine

Purine nucleobases (adenosine and guanosine) are prone to spontaneous breaking of the nucleosidic bond to form abasic sites in both DNA and RNA. However, the purines also undergo photochemical reactions, including oxidation and cycloaddition, to form damage sites, albeit at lower photochemical quantum yields than the pyrimidines. In this study, we use ultraviolet resonance Raman spectroscopy to measure the initial excited-state structural dynamics in the nucleoside, 2′-deoxyadenosine. The resonance Raman-derived initial excited-state structural dynamics throughout the 260 nm La excited electronic state of adenine are found to be smaller in the nucleoside than in the previously reported 9-methyladenine nucleobase derivative, consistent with what is found for the pyrimidines thymine and uracil. Interestingly, resonance-enhanced vibrational modes in this electronic state also contain internal coordinates localized on the sugar, which may represent a different energy dissipation mechanism than in the pyrimidine nucleosides. The results will be discussed in terms of the initial excited-state photophysics and photochemistry of DNA and RNA.

Excited-State Dynamics of Thienoguanosine, an Isomorphic Highly Fluorescent Analogue of Guanosine

Thienoguanosine (^th G) is an isomorphic analogue of guanosine with promising potentialities as fluorescent DNA label. As a free probe in protic solvents, ^th G exists in two tautomeric forms, identified as the H1, being the only one observed in nonprotic solvents, and H3 keto-amino tautomers. We herein investigate the photophysics of ^th G in solvents of different polarity, from water to dioxane, by combining time-resolved fluorescence with PCM/TD-DFT and CASSCF calculations. Fluorescence lifetimes of 14.5-20.5 and 7-13 ns were observed for the H1 and H3 tautomers, respectively, in the tested solvents. In methanol and ethanol, an additional fluorescent decay lifetime (≈3 ns) at the blue emission side (λ≈430 nm) as well as a 0.5 ns component with negative amplitude at the red edge of the spectrum, typical of an excited-state reaction, were observed. Our computational analysis explains the solvent effects observed on the tautomeric equilibrium. The main radiative and nonradiative deactivation routes have been mapped by PCM/TD-DFT calculations in solution and CASSCF in the gas phase. The most easily accessible conical intersection, involving an out-of plane motion of the sulfur atom in the five-membered ring of ^th G, is separated by a sizeable energy barrier (≥0.4 eV) from the minimum of the spectroscopic state, which explains the large experimental fluorescence quantum yield.

Photochemical Relaxation Pathways of 9 H-8-Azaguanine and 8 H-8-Azaguanine

The photochemical reaction path approach, the MS-CASPT2 quantum-chemical method, and double-ζ basis sets (cc-pVDZ) were used to study 9 H-8-azaguanine and 8 H-8-azaguanine relaxation pathways. Several potential energy hypersurfaces were characterized by means of equilibrium geometries, surface crossings (conical intersections and singlet-triplet intersystem crossings), minimum energy paths, and linear interpolation in internal coordinates. The 9 H-8-azaguanine main photochemical event begins with the direct population of the ¹(ππ* L_a) state, which evolves toward a conical intersection with the ground state after surmounting a small energy barrier, explaining why it is nonfluorescent. For 8 H-8-azaguanine, two relaxation mechanisms are possible, depending on the excitation energy. If the S₁ ¹(ππ*) state is initially populated (lower-energy region), the system evolves barrierless to the S₁ ¹(ππ*)_min region, from where the excess energy is released. If the ¹(ππ* L_a) state is populated (higher-energy radiation), the system will encounter conical intersections with the S₂ ¹(n_Oπ*) and S₁ ¹(ππ*) states before evolving to the ¹(ππ* L_a)_min region, from where a conical intersection with the ground state is accessible, favoring radiationless deactivation to the ground state. However, because a fraction of the population can be transferred from ¹(ππ* L_a) to the S₁ ¹(ππ*) state, emission from the S₁ ¹(ππ*)_min region is also expected, although weaker than it would be if the S₁ ¹(ππ*) state were populated directly. Irrespective of the excitation energy, the emissive state is the same and a single fluorescence band is observed, with the strongest emission occurring upon excitation in the lower-energy region, as observed experimentally. Therefore, our computational study corroborates experimental results, attributing the emission of the neutral form of 8-azaguanine in solution to the presence of the minor 8 H-8-azaguanine tautomer, while the 9 H-8-azaguanine major tautomer is nonfluorescent.

Nonadiabatic dynamics: The SHARC approach

We review the Surface Hopping including ARbitrary Couplings (SHARC) approach for excited-state nonadiabatic dynamics simulations. As a generalization of the popular surface hopping method, SHARC allows simulating the full-dimensional dynamics of molecules including any type of coupling terms beyond nonadiabatic couplings. Examples of these arbitrary couplings include spin-orbit couplings or dipole moment-laser field couplings, such that SHARC can describe ultrafast internal conversion, intersystem crossing, and radiative processes. The key step of the SHARC approach consists of a diagonalization of the Hamiltonian including these couplings, such that the nuclear dynamics is carried out on potential energy surfaces including the effects of the couplings-this is critical in any applications considering, for example, transition metal complexes or strong laser fields. We also give an overview over the new SHARC2.0 dynamics software package, released under the GNU General Public License, which implements the SHARC approach and several analysis tools. The review closes with a brief survey of applications where SHARC was employed to study the nonadiabatic dynamics of a wide range of molecular systems. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and KineticsSoftware > Simulation MethodsSoftware > Quantum Chemistry.

Excited State Decay Pathways of 2'-Deoxy-5-methylcytidine and Deoxycytidine Revisited in Solution: A Comprehensive Kinetic Study by Femtosecond Transient Absorption

Methylated cytosine is proved to have an important role as an epigenetic signal in gene regulation and is often referred to "the fifth base of DNA". A comprehensive understanding of the electronic excited state relaxation in cytosine and its methylated derivatives is crucial for revealing UV-induced photodamage to the biological genome. Because of the existence of multiple closely lying "bright" and "dark" excited states, the decay pathways in these DNA nucleosides are the most complex and the least understood so far. In this study, femtosecond transient absorption with different excitation wavelengths (240-296 nm) was used to study the relaxation of excited electronic states of 5-methylcytosine (5mC) and 2'-deoxy-5-methylcytidine (5mdCyd) in phosphate buffered aqueous solution and in acetonitrile solution. Two distinct nonradiative decay channels were directly observed. The first one is a several picosecond internal conversion channel that involves two bright ππ* states (ππ*₂ and ππ*₁) when ππ*₂ state is initially populated. The second channel contains the lower energy ππ*₁ state and a so far experimental unidentified long-lived state which exhibits a several nanosecond lifetime. The long-lived state can only be accessed by the initially excited ππ*₁ state. Inspired by this new discovery in 5mC and 5mdCyd, we revisited the decay of excited state of 2'-deoxycytidine (dCyd), revealing very similar decay pathways. Additionally, a well-known dark n_Oπ* state (carbonyl lone pair) with ∼30 ps lifetime is present in both decay channels in dCyd. With our detailed experimental results, we successfully reconcile the long history debate of cytosine excited state relaxation mechanism by pointing out that the reason for the complex dynamics under traditional 266 nm excitation is mixed signals from the above-mentioned two distinct decay pathways. Our findings lead to a dramatically different and new picture of electronic energy relaxation in 5mdCyd/dCyd and could help to understand photostability as well as UV-induced photodamage of these nucleotides and related DNAs.

Decay Pathways of Thymine Revisited

The decay of electronically excited states of thymine (Thy) and thymidine 5'-monophosphate (TMP) was studied by time-resolved UV/vis and IR spectroscopy. In addition to the well-established ultrafast internal conversion to the ground state, a so far unidentified UV-induced species is observed. In D₂O, this species decays with a time constant of 300 ps for thymine and of 1 ns for TMP. The species coexists with the lowest triplet state and is formed with a comparably high quantum yield of about 10% independent of the solvent. The experimentally determined spectral signatures are discussed in the light of quantum chemical calculations of the singlet and triplet excited states of thymine.

Internal conversion and intersystem crossing pathways in UV excited, isolated uracils and their implications in prebiotic chemistry

The photodynamic properties of molecules determine their ability to survive in harsh radiation environments. As such, the photostability of heterocyclic aromatic compounds to electromagnetic radiation is expected to have been one of the selection pressures influencing the prebiotic chemistry on early Earth. In the present study, the gas-phase photodynamics of uracil, 5-methyluracil (thymine) and 2-thiouracil-three heterocyclic compounds thought to be present during this era-are assessed in the context of their recently proposed intersystem crossing pathways that compete with internal conversion to the ground state. Specifically, time-resolved photoelectron spectroscopy measurements evidence femtosecond to picosecond timescales for relaxation of the singlet (1)ππ* and (1)nπ* states as well as for intersystem crossing to the triplet manifold. Trapping in the excited triplet state and intersystem crossing back to the ground state are investigated as potential factors contributing to the susceptibility of these molecules to ultraviolet photodamage.

Solvatochromic Effects on the Absorption Spectrum of 2-Thiocytosine

The solvatochromic effects of six different solvents on the UV absorption spectrum of 2-thiocytosine have been studied by a combination of experimental and theoretical techniques. The steady-state absorption spectra show significant shifts of the absorption bands, where in more polar solvents the first absorption maximum shifts to higher transition energies and the second maximum to lower energies. The observed solvatochromic shifts have been rationalized using three popular solvatochromic scales and with high-level multireference quantum chemistry calculations including implicit and explicit solvent effects. It has been found that the dipole moments of the excited states account for some general shifts in the excitation energies, whereas the explicit solvent interactions explain the differences in the spectra recorded in the different solvents.

Insights into the deactivation of 5-bromouracil after ultraviolet excitation

5-Bromouracil is a nucleobase analogue that can replace thymine in DNA strands and acts as a strong radiosensitizer, with potential applications in molecular biology and cancer therapy. Here, the deactivation of 5-bromouracil after ultraviolet irradiation is investigated in the singlet and triplet manifold by accurate quantum chemistry calculations and non-adiabatic dynamics simulations. It is found that, after irradiation to the bright ππ* state, three main relaxation pathways are, in principle, possible: relaxation back to the ground state, intersystem crossing (ISC) and C-Br photodissociation. Based on accurate MS-CASPT2 optimizations, we propose that ground-state relaxation should be the predominant deactivation pathway in the gas phase. We then employ different electronic structure methods to assess their suitability to carry out excited-state dynamics simulations. MRCIS (multi-reference configuration interaction including single excitations) was used in surface hopping simulations to compute the ultrafast ISC dynamics, which mostly involves the 1nOπ* and 3ππ* states.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

Assessment of the Potential Energy Hypersurfaces in Thymine within Multiconfigurational Theory: CASSCF vs. CASPT2

The present study provides new insights into the topography of the potential energy hypersurfaces (PEHs) of the thymine nucleobase in order to rationalize its main ultrafast photochemical decay paths by employing two methodologies based on the complete active space self-consistent field (CASSCF) and the complete active space second-order perturbation theory (CASPT2) methods: (i) CASSCF optimized structures and energies corrected with the CASPT2 method at the CASSCF geometries and (ii) CASPT2 optimized geometries and energies. A direct comparison between these strategies is drawn, yielding qualitatively similar results within a static framework. A number of analyses are performed to assess the accuracy of these different computational strategies under study based on a variety of numerical thresholds and optimization methods. Several basis sets and active spaces have also been calibrated to understand to what extent they can influence the resulting geometries and subsequent interpretation of the photochemical decay channels. The study shows small discrepancies between CASSCF and CASPT2 PEHs, displaying a shallow planar or twisted ¹(ππ*) minimum, respectively, and thus featuring a qualitatively similar scenario for supporting the ultrafast bi-exponential deactivation registered in thymine upon UV-light exposure. A deeper knowledge of the PEHs at different levels of theory provides useful insight into its correct characterization and subsequent interpretation of the experimental observations. The discrepancies displayed by the different methods studied here are then discussed and framed within their potential consequences in on-the-fly non-adiabatic molecular dynamics simulations, where qualitatively diverse outcomes are expected.

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.

From non-covalent binding to irreversible DNA lesions: nile blue and nile red as photosensitizing agents

We report a molecular modeling study, coupled with spectroscopy experiments, on the behavior of two well known organic dyes, nile blue and nile red, when interacting with B-DNA. In particular, we evidence the presence of two competitive binding modes, for both drugs. However their subsequent photophysical behavior is different and only nile blue is able to induce DNA photosensitization via an electron transfer mechanism. Most notably, even in the case of nile blue, its sensitization capabilities strongly depend on the environment resulting in a single active binding mode: the minor groove. Fluorescence spectroscopy confirms the presence of competitive interaction modes for both sensitizers, while the sensitization via electron transfer, is possible only in the case of nile blue.

Simulation of Multi-Dimensional Signals in the Optical Domain: Quantum-Classical Feedback in Nonlinear Exciton Propagation

We present an algorithm for the simulation of nonlinear 2D spectra of molecular systems in the UV-vis spectral region from atomistic molecular dynamics trajectories subject to nonadiabatic relaxation. We combine the nonlinear exciton propagation (NEP) protocol, that relies on a quasiparticle approach with the surface hopping methodology to account for quantum-classical feedback during the dynamics. Phenomena, such as dynamic Stokes shift due to nuclear relaxation, spectral diffusion, and population transfer among electronic states, are thus naturally included and benchmarked on a model of two electronic states coupled to a harmonic coordinate and a classical heatbath. The capabilities of the algorithm are further demonstrated for the bichromophore diphenylmethane that is described in a fully microscopic fashion including all 69 classical nuclear degrees of freedom. We demonstrate that simulated 2D signals are especially sensitive to the applied theoretical approximations (i.e., choice of active space in the CASSCF method) where population dynamics appears comparable.

Intersystem Crossing Pathways in the Noncanonical Nucleobase 2-Thiouracil: A Time-Dependent Picture

The deactivation mechanism after ultraviolet irradiation of 2-thiouracil has been investigated using nonadiabatic dynamics simulations at the MS-CASPT2 level of theory. It is found that after excitation the S2 quickly relaxes to S1, and from there intersystem crossing takes place to both T2 and T1 with a time constant of 400 fs and a triplet yield above 80%, in very good agreement with recent femtosecond experiments in solution. Both indirect S1 → T2 → T1 and direct S1 → T1 pathways contribute to intersystem crossing, with the former being predominant. The results contribute to the understanding of how some noncanonical nucleobases respond to harmful ultraviolet light, which could be relevant for prospective photochemotherapeutic applications.

Excimers and Exciplexes in Photoinitiated Processes of Oligonucleotides

A lot has been learned about the physical and chemical transformations that originate from the absorption of light by DNA, and computational chemistry has played a critical role in revealing the mechanisms of how these transformations occur. Nucleic acids consist of chromophores interacting via π stacking and hydrogen bonding. The fate of these systems after they absorb light is determined by the interplay and competition between pathways involving one chromophore or interacting chromophores. This Perspective highlights the role of π stacking in photophysical and photochemical processes in oligonucleotides and reveals the importance of excimers and exciplexes. Special types of excimers/exciplexes, characterized as bonded excimers/exciplexes, are also found to be important.

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.

Direct Observation of Triplet-State Population Dynamics in the RNA Uracil Derivative 1-Cyclohexyluracil

Investigation of the excited-state dynamics in nucleic acid monomers is an area of active research due to the crucial role these early events play in DNA and RNA photodamage. The dynamics and rate at which the triplet state is populated are key mechanistic pathways yet to be fully elucidated. Direct spectroscopic evidence is presented in this contribution for intersystem crossing dynamics in a uracil derivative, 1-cyclohexyluracil. It is shown that intersystem crossing to the triplet manifold occurs in one picosecond or less in acetonitrile solution-at least an order of magnitude faster than previously estimated experimentally. Broadband transient absorption measurements also reveal the primary electronic relaxation pathways of the uracil chromophore, including the absorption spectra of the (1)ππ*, (1)nπ*, and (3)ππ* states and the rates of vibrational cooling in the ground and (3)ππ* states. The experimental results are supported by density functional calculations.

Intramolecular vibrational redistribution in the non-radiative excited state decay of uracil in the gas phase: an ab initio molecular dynamics study

We report a study of intramolecular vibrational distribution (IVR) occurring in the electronic ground state of uracil (S0) in the gas phase, following photoexcitation in the lowest energy bright excited state (Sπ) and decay through the ethylene-like Sπ/S0 Conical Intersection (CI-0π). To this aim we have performed 20 independent ab initio molecular dynamics simulations starting from CI-0π (ten of them with 1 eV kinetic energy randomly distributed over the different molecular degrees of freedom) and 10 starting from the ground state minimum (Franck-Condon, FC, point), with the excess kinetic energy equal to the energy gap between CI-0π and the FC point. The simulations, exploiting PBE0/6-31G(d) calculations, were performed over an overall period of 10 ps. A thorough statistical analysis of the variation of the geometrical parameters of uracil during the simulation time and of the distribution of the kinetic energy among the different vibrational degrees of freedom provides a consistent picture of the IVR process. In the first 0-200 fs the structural dynamics involve mainly the recovery of the average planarity. In the 200-600 fs time range, a substantial activation of CO and NH degrees of freedom is observed. After 500-600 fs most of the geometrical parameters reach average values similar to those found after 10 ps, though the system cannot be considered to be in equilibrium yet.

Deciphering the photochemical mechanisms describing the UV-induced processes occurring in solvated guanine monophosphate

The photophysics and photochemistry of water-solvated guanine monophosphate (GMP) are here characterized by means of a multireference quantum-chemical/molecular mechanics theoretical approach (CASPT2//CASSCF/AMBER) in order to elucidate the main photo-processes occurring upon UV-light irradiation. The effect of the solvent and of the phosphate group on the energetics and structural features of this system are evaluated for the first time employing high-level ab initio methods and thoroughly compared to those in vacuo previously reported in the literature and to the experimental evidence to assess to which extent they influence the photoinduced mechanisms. Solvated electronic excitation energies of solvated GMP at the Franck-Condon (FC) region show a red shift for the ππ(*) La and Lb states, whereas the energy of the oxygen lone-pair nπ(*) state is blue-shifted. The main photoinduced decay route is promoted through a ring-puckering motion along the bright lowest-lying La state toward a conical intersection (CI) with the ground state, involving a very shallow stationary point along the minimum energy pathway in contrast to the barrierless profile found in gas-phase, the point being placed at the end of the minimum energy path (MEP) thus endorsing its ultrafast deactivation in accordance with time-resolved transient and photoelectron spectroscopy experiments. The role of the nπ(*) state in the solvated system is severely diminished as the crossings with the initially populated La state and also with the Lb state are placed too high energetically to partake prominently in the deactivation photo-process. The proposed mechanism present in solvated and in vacuo DNA/RNA chromophores validates the intrinsic photostability mechanism through CI-mediated non-radiative processes accompanying the bright excited-state population toward the ground state and subsequent relaxation back to the FC region.

Electronic and structural elements that regulate the excited-state dynamics in purine nucleobase derivatives

The excited-state dynamics of the purine free base and 9-methylpurine are investigated using experimental and theoretical methods. Femtosecond broadband transient absorption experiments reveal that excitation of these purine derivatives in aqueous solution at 266 nm results primarily in ultrafast conversion of the S2(ππ*) state to the vibrationally excited (1)nπ* state. Following vibrational and conformational relaxation, the (1)nπ* state acts as a doorway state in the efficient population of the triplet manifold with an intersystem crossing lifetime of hundreds of picoseconds. Experiments show an almost 2-fold increase in the intersystem crossing rate on going from polar aprotic to nonpolar solvents, suggesting that a solvent-dependent energy barrier must be surmounted to access the singlet-to-triplet crossing region. Ab initio static and surface-hopping dynamics simulations lend strong support to the proposed relaxation mechanism. Collectively, the experimental and computational results demonstrate that the accessibility of the nπ* states and the topology of the potential energy surfaces in the vicinity of conical intersections are key elements in controlling the excited-state dynamics of the purine derivatives. From a structural perspective, it is shown that the purine chromophore is not responsible for the ultrafast internal conversion in the adenine and guanine monomers. Instead, C6 functionalization plays an important role in regulating the rates of radiative and nonradiative relaxation. C6 functionalization inhibits access to the (1)nπ* state while simultaneously facilitating access to the (1)ππ*(La)/S0 conical intersection, such that population of the (1)nπ* state cannot compete with the relaxation pathways to the ground state involving ring puckering at the C2 position.

Ultrafast Excited-State Deactivation of 8-Hydroxy-2'-deoxyguanosine Studied by Femtosecond Fluorescence Spectroscopy and Quantum-Chemical Calculations

The fluorescence properties of the 8-hydroxy-2'-deoxyguanosine (8-oxo-dG) in aqueous solution at pH 6.5 are studied by steady-state spectroscopy and femtosecond fluorescence up-conversion and compared with those of 2'-deoxyguanine (dG) and 2'-deoxyguanine monophosphate (dGMP). The steady-state fluorescence spectrum of 8-oxo-dG is composed of a broad band that peaks at 356 nm and extends over the entire visible spectral region, and its fluorescence quantum yield is twice that of dG/dGMP. After excitation at 267 nm, the initial fluorescence anisotropy at all wavelengths is lower than 0.1, giving evidence of an ultrafast internal conversion (<100 fs) between the two lowest excited ππ* states (Lb and La). The fluorescence decays of 8-oxo-dG are biexponential with an average lifetime of 0.7 ± 0.1 ps, which is about two times longer than that of dGMP. In contrast with dGMP, only a moderate dynamical shift (∼1400 vs 10,000 cm(-1)) of the fluorescence spectra of 8-oxo-dG is observed on the time scale of a few picoseconds without modification of the spectral shape. PCM/TD-DFT calculations, employing either the PBE0 or the M052X functionals, provide absorption spectra in good agreement with the experimental one and show that the deactivation path is similar to that proposed for dGMP, with a fast motion toward an energy plateau, where the purine ring keeps an almost planar geometry, followed by decay to S0, via out-of-the plane motion of amino substituent.

Photoinduced processes in nucleic acids

Photoinduced processes in nucleic acids are phenomena of fundamental interest in diverse fields, from prebiotic studies, through medical research on carcinogenesis, to the development of bioorganic photodevices. In this contribution we survey many aspects of the research across the boundaries. Starting from a historical background, where the main milestones are identified, we review the main findings of the physical-chemical research of photoinduced processes on several types of nucleic-acid fragments, from monomers to duplexes. We also discuss a number of different issues which are still under debate.