Monitoring Non-Adiabatic Dynamics of the RNA Base Uracil by UV-Pump-IR-Probe Spectroscopy

Comparative Study of Uracil Excited-State Photophysics in Water and Acetonitrile via RMS-CASPT2-Driven Quantum-Classical Trajectories

We present a nonadiabatic molecular dynamics study of the ultrafast processes occurring in uracil upon UV light absorption, leading to electronic excitation and subsequent nonradiative decay. Previous studies have indicated that the mechanistic details of this process are drastically different depending on whether the process takes place in the gas phase, acetonitrile, or water. However, such results have been produced using quantum chemical methods that did not incorporate both static and dynamic electron correlation. In order to assess the previously proposed mechanisms, we simulate the photodynamics of uracil in the three environments mentioned above using quantum-classical trajectories and, for solvated uracil, hybrid quantum mechanics/molecular mechanics (QM/MM) models driven by the rotated multistate complete active space second-order perturbation (RMS-CASPT2) method. To do so, we exploit the gradient recently made available in OpenMolcas and compare the results to those obtained using the complete active space self-consistent field (CASSCF) method only accounting for static electron correlation. We show that RMS-CASPT2 produces, in general, a mechanistic picture different from the one obtained at the CASSCF level but confirms the hypothesis advanced on the basis of previous ROKS and TDDFT studies thus highlighting the importance of incorporating dynamic electron correlation in the investigation of ultrafast electronic deactivation processes.

Photophysics of uracil: an explicit time-dependent generating function-based method combining both nonadiabatic and spin-orbit coupling effects

We present a composite framework for calculating the rates of non-radiative deactivation processes, namely internal conversion (IC) and intersystem crossing (ISC), on an equal footing by explicitly computing the non-adiabatic coupling (NAC) and spin-orbit coupling (SOC) constants, respectively. The stationary-state approach uses a time-dependent generating function based on Fermi's golden rule. We validate the applicability of the framework by computing the rate of IC for azulene, obtaining comparable rates to experimental and previous theoretical results. Next, we investigate the photophysics associated with the complex photodynamics of the uracil molecule. Interestingly, our simulated rates corroborate experimental observations. Detailed analyses using Duschinsky rotation matrices, displacement vectors and NAC matrix elements are presented to interpret the findings alongside testing the suitability of the approach for such molecular systems. The suitability of the Fermi's golden rule based method is explained qualitatively in terms of single-mode potential energy surfaces.

Newton-X Platform: New Software Developments for Surface Hopping and Nuclear Ensembles

Newton-X is an open-source computational platform to perform nonadiabatic molecular dynamics based on surface hopping and spectrum simulations using the nuclear ensemble approach. Both are among the most common methodologies in computational chemistry for photophysical and photochemical investigations. This paper describes the main features of these methods and how they are implemented in Newton-X. It emphasizes the newest developments, including zero-point-energy leakage correction, dynamics on complex-valued potential energy surfaces, dynamics induced by incoherent light, dynamics based on machine-learning potentials, exciton dynamics of multiple chromophores, and supervised and unsupervised machine learning techniques. Newton-X is interfaced with several third-party quantum-chemistry programs, spanning a broad spectrum of electronic structure methods.

A molecular insight into the photophysics of barbituric acid, a candidate for canonical nucleobases' ancestor

This work investigates the photophysics of barbituric acid at different pH conditions using ab initio methods. Our calculations ascribe the most intense bands at ca. 260 nm at neutral pH and 210 nm at acidic pH conditions in the absorption spectra of this chromophore to the lowest lying ππ* transitions. Consistently with the ultrashort excited state lifetimes experimentally registered, the potential energy landscapes of both the neutral and deprotonated forms of barbituric acid combined with the interpretation of their transient absorption spectra suggest the deactivation of these systems along the singlet manifold. Compared to uracil, its closest natural nucleobase, barbituric acid presents a red shifted absorption spectrum, due to the lowering by more than 0.5 eV of the lowest-energy ππ* excited state, and a much more complex topography of the S₁ potential energy surface, with several energetically accessible local minima. This fact, however, does not affect the excited state lifetimes, which for barbituric acid were experimentally registered in the sub-ps time scale.

Time Resolved Photoelectron Spectroscopy as a Test of Electronic Structure and Nonadiabatic Dynamics

We compare different levels of theory for simulating excited state molecular dynamics and use time-resolved photoelectron spectroscopy measurements to benchmark the theory. We perform trajectory surface hopping simulations for uracil excited to the first bright state (ππ*) using three different levels of theory (CASSCF, MRCIS, and XMS-CASPT2) in order to understand the role of dynamical correlation in determining the excited state dynamics, with a focus on the coupling between different electronic states and internal conversion back to the ground state. These dynamics calculations are used to simulate the time-resolved photoelectron spectra. The comparison of the calculated and measured spectra allows us to draw conclusions regarding the relative insights and quantitative accuracy of the calculations at the three different levels of theory, demonstrating that detailed quantitative comparisons of time-resolved photoelectron spectra can be used to benchmark methodology.

Simulation of UV absorption spectra and relaxation dynamics of uracil and uracil-water clusters

Despite many studies, the mechanisms of nonradiative relaxation of uracil in the gas phase and in aqueous solution are still not fully resolved. Here we combine theoretical UV absorption spectroscopy with nonadiabatic dynamics simulations to identify the photophysical mechanisms that can give rise to experimentally observed decay time constants. We first compute and theoretically assign the electronic spectra of uracil using the second-order algebraic-diagrammatic-construction (ADC(2)) method. The obtained electronic states, their energy differences and state-specific solvation effects are the prerequisites for understanding the photodynamics. We then use nonadiabatic trajectory-surface-hopping dynamics simulations to investigate the photoinduced dynamics of uracil and uracil-water clusters. In contrast to previous studies, we found that a single mechanism - the ethylenic twist around the C[double bond, length as m-dash]C bond - is responsible for the ultrafast component of the nonradiative decay, both in the gas phase and in solution. Very good agreement with the experimentally determined ultrashort decay time constants is obtained.

Effect of dynamic correlation on the ultrafast relaxation of uracil in the gas phase

High level multi-reference non-adiabatic dynamics simulations reveal that uracil’s photoexcited S₂ state decays very quickly without any significant trapping.

Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going

Computational spectroscopy is becoming a mandatory tool for the interpretation of the complex, and often congested, spectral maps delivered by modern non-linear multi-pulse techniques. The fields of Electronic Structure Methods, Non-Adiabatic Molecular Dynamics, and Theoretical Spectroscopy represent the three pillars of the virtual ultrafast optical spectrometer, able to deliver transient spectra in silico from first principles. A successful simulation strategy requires a synergistic approach that balances between the three fields, each one having its very own challenges and bottlenecks. The aim of this Perspective is to demonstrate that, despite these challenges, an impressive agreement between theory and experiment is achievable now regarding the modeling of ultrafast photoinduced processes in complex molecular architectures. Beyond that, some key recent developments in the three fields are presented that we believe will have major impacts on spectroscopic simulations in the very near future. Potential directions of development, pending challenges, and rising opportunities are illustrated.

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.

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

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

Recent advances in multidimensional ultrafast spectroscopy

Multidimensional ultrafast spectroscopies are one of the premier tools to investigate condensed phase dynamics of biological, chemical and functional nanomaterial systems. As they reach maturity, the variety of frequency domains that can be explored has vastly increased, with experimental techniques capable of correlating excitation and emission frequencies from the terahertz through to the ultraviolet. Some of the most recent innovations also include extreme cross-peak spectroscopies that directly correlate the dynamics of electronic and vibrational states. This review article summarizes the key technological advances that have permitted these recent advances, and the insights gained from new multidimensional spectroscopic probes.

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

Probing the Interstate Coupling near a Conical Intersection by Optical Spectroscopy

Conical intersections are points where adiabatic potential energy surfaces cross. The interstate coupling between the potential energy surfaces plays a crucial role in many processes associated with conical intersections. Still no method exists to measure this coupling driving the chemical reactions between the potential energy surfaces involved. In this Letter, using a generic model for photoisomerization, we propose a novel experimental approach to estimate the coupling that mixes the electronic states near a conical intersection. The approach is based on analyzing the vibrational wavepacket of the reactant in the adiabatic ground and excited electronic states. The nuclear wavepacket dynamics are extracted from linear absorption and two-dimensional electronic spectroscopy. Comparing the frequencies of the coupling mode in the adiabatic ground and excited states from models with and without coupling between the potential energy surfaces suggests an experimental tool to determine the interstate coupling.

Monitoring Nonadiabatic Electron-Nuclear Dynamics in Molecules by Attosecond Streaking of Photoelectrons

Streaking of photoelectrons has long been used for the temporal characterization of attosecond extreme ultraviolet pulses. When the time-resolved photoelectrons originate from a coherent superposition of electronic states, they carry additional phase information, which can be retrieved by the streaking technique. In this contribution we extend the streaking formalism to include coupled electron and nuclear dynamics in molecules as well as initial coherences. We demonstrate how streaked photoelectrons offer a novel tool for monitoring nonadiabatic dynamics as it occurs in the vicinity of conical intersections and avoided crossings. Streaking can provide high time resolution direct signatures of electronic coherences, which affect many primary photochemical and biological events.

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.

Spectral lineshapes in nonlinear electronic spectroscopy

We outline a computational approach for nonlinear electronic spectra, which accounts for the electronic energy fluctuations due to nuclear degrees of freedom and explicitly incorporates the fluctuations of higher excited states, induced by the dynamics in the photoactive state(s). This approach is based on mixed quantum-classical dynamics simulations. Tedious averaging over multiple trajectories is avoided by employing the linearly displaced Brownian harmonic oscillator to model the correlation functions. The present strategy couples accurate computations of the high-lying excited state manifold with dynamics simulations. The application is made to the two-dimensional electronic spectra of pyrene, a polycyclic aromatic hydrocarbon characterized by an ultrafast (few tens of femtoseconds) decay from the bright S2 state to the dark S1 state. The spectra for waiting times t2 = 0 and t2 = 1 ps demonstrate the ability of this approach to model electronic state fluctuations and realistic lineshapes. Comparison with experimental spectra [Krebs et al., New Journal of Physics, 2013, 15, 085016] shows excellent agreement and allows us to unambiguously assign the excited state absorption features.

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.

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.

Probing electronic and vibrational dynamics in molecules by time-resolved photoelectron, Auger-electron, and X-ray photon scattering spectroscopy

We present a unified description for time-resolved electron and photon scattering spectroscopies from molecules prepared in nonstationary states. Signals are expressed in terms of superoperator Green's functions and a systematic procedure for treating various degrees of freedom consistently at different levels of theory is developed. The standard Fermi Golden Rule expressions for photoelectron spectra, which are limited to broad, slowly-varying signals, are obtained as a limiting case of our more general theory that applies to broader parameter regimes.

Ultrafast intersystem crossing dynamics in uracil unravelled by ab initio molecular dynamics

Ab initio molecular dynamics simulations have been performed in order to investigate the relaxation dynamics of uracil after UV excitation in gas phase. Intersystem crossing (ISC) has been included for the first time into time-dependent simulations of uracil, allowing the system to relax in the singlet as well as in the triplet states. The results show a qualitatively different picture than similar simulations that include singlet states only. The inclusion of ISC effectively quenches the relaxation to the singlet ground state and instead privileges transitions from the low-lying nπ* state (S1) to a ππ* triplet state (T2) followed by rapid internal conversion to the lowest triplet state.

Femtosecond stimulated Raman spectroscopy of the cyclobutane thymine dimer repair mechanism: a computational study

Cyclobutane thymine dimer, one of the major lesions in DNA formed by exposure to UV sunlight, is repaired in a photoreactivation process, which is essential to maintain life. The molecular mechanism of the central step, i.e., intradimer C-C bond splitting, still remains an open question. In a simulation study, we demonstrate how the time evolution of characteristic marker bands (C═O and C═C/C-C stretch vibrations) of cyclobutane thymine dimer and thymine dinucleotide radical anion, thymidylyl(3'→5')thymidine, can be directly probed with femtosecond stimulated Raman spectroscopy (FSRS). We construct a DFT(M05-2X) potential energy surface with two minor barriers for the intradimer C₅-C₅' splitting and a main barrier for the C₆-C₆' splitting, and identify the appearance of two C₅═C₆ stretch vibrations due to the C₆-C₆' splitting as a spectroscopic signature of the underlying bond splitting mechanism. The sequential mechanism shows only absorptive features in the simulated FSRS signals, whereas the fast concerted mechanism shows characteristic dispersive line shapes.

Time-resolved broadband Raman spectroscopies: a unified six-wave-mixing representation

Excited-state vibrational dynamics in molecules can be studied by an electronically off-resonant Raman process induced by a probe pulse with variable delay with respect to an actinic pulse. We establish the connection between several variants of the technique that involve either spontaneous or stimulated Raman detection and different pulse configurations. By using loop diagrams in the frequency domain, we show that all signals can be described as six wave mixing which depend on the same four point molecular correlation functions involving two transition dipoles and two polarizabilities and accompanied by a different gating. Simulations for the stochastic two-state-jump model illustrate the origin of the absorptive and dispersive features observed experimentally.

Probing the Conical Intersection Dynamics of the RNA Base Uracil by UV-Pump Stimulated-Raman-Probe Signals; Ab Initio Simulations

Nonadiabatic electron and nuclear dynamics of photoexcited molecules involving conical intersections is of fundamental importance in many reactions such as the self-protection mechanism of DNA and RNA bases against UV irradiation. Nonlinear vibrational spectroscopy can provide an ultrafast sensitive probe for these processes. We employ a simulation protocol that combines nonadiabatic on-the-fly molecular dynamics with a mode-tracking algorithm for the simulation of femtosecond stimulated Raman spectroscopy (SRS) signals of the high frequency C–H- and N–H-stretch vibrations of the photoexcited RNA base uracil. The simulations rely on a microscopically derived expression that takes into account the path integral of the excited state evolution and the pulse shapes. Analysis of the joint time/frequency resolution of the technique reveals a matter chirp contribution that limits the inherent temporal resolution. Characteristic signatures of relaxation dynamics mediated in the vicinity of conical intersection are predicted. The C–H and N–H spectator modes provide high sensitivity to their local environment and act as local probes with submolecular and high temporal resolution.

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.

Broadband infrared and Raman probes of excited-state vibrational molecular dynamics: simulation protocols based on loop diagrams

Vibrational motions in electronically excited states can be observed either by time and frequency resolved infrared absorption or by off resonant stimulated Raman techniques. Multipoint correlation function expressions are derived for both signals. Three representations which suggest different simulation protocols for the signals are developed. These are based on the forward and the backward propagation of the wavefunction, sum over state expansion using an effective vibrational Hamiltonian or a semiclassical treatment of a bath. We show that the effective temporal (Δt) and spectral (Δω) resolution of the techniques is not controlled solely by experimental knobs but also depends on the system dynamics being probed. The Fourier uncertainty ΔωΔt > 1 is never violated.