Nonadiabatic Deactivation of 9H-Adenine: A Comprehensive Picture Based on Mixed Quantum−Classical Dynamics

Theoretical insight into photodeactivation mechanisms of adenine-uracil and adenine-thymine nucleobase pairs

In this work, several plausible intra- and intermolecular photoinduced processes of the Watson-Crick base pairs of adenine with uracil (A-U) or thymine (A-T) according to the results of spin component scaling variant of algebraic diagrammatic construction up to the second order [SCS-ADC(2)] calculations are discussed. Although widely explored, these systems lack complete characterization of possible intramolecular relaxation channels perturbed by intermolecular interactions. In particular, we address the still open debate on photodeactivation via purine-ring puckering at the C2 or C6-atom position of adenine. We also show that the presence of low-lying, long-lived 1nπ* states can be a significant factor in hindering relaxation via an electron-driven proton transfer process, as the population of these states can lead to an efficient intersystem crossing to a triplet manifold, the estimated rate of which is 1.6 × 1010 s-1 which exceeds the corresponding internal conversion to the ground state by an order of magnitude. Additionally, the SCS variant of the ADC(2) method is shown to provide a more balanced description of valence and charge-transfer excited states.

An Exhaustive Quantum-Classical Study of C6F6+ Using the Newly Formulated Parallel TDDVR Method

We recently implemented our parallelized quantum-classical dynamical approach, known as the Time-Dependent Discrete Variable Representation (TDDVR) method, which is applied to the spectroscopically important hexafluorobenzene (HFBz) radical cation, where several conical intersections exist in their seven lowest excited electronic states (S11B2u, S21E1g, S31B1u, S41E1u, and S51A2u) considering degeneracy among potential energy surfaces (PESs), to demonstrate their various dynamical aspects. This new parallel version shows almost linear scalability with increasing number of computing processors. To get photoelectron (PE) spectra, Mass-Analyzed Threshold Ionization (MATI) spectra, population dynamics, and many other dynamical observables, the first-principles dynamics is applied at the state-of-the-art level to the corresponding Hamiltonian, where the Jahn-Teller (JT) and pseudo-Jahn-Teller (PJT) type interactions are involved in those coupled seven electronic states. The quantum-classical method is used to thoroughly analyze the effects of these couplings on the nuclear dynamics of the involved electronic states, and the findings are compared with those observables obtained from experiments. Intrinsic dynamical properties are explained using the reduced densities of the wave packet (WP) in a coupled electronic manifold. The PE and MATI spectra of HFBz computed using TDDVR are found to be in good agreement with earlier experimental data and other theoretically simulated spectra.

Extreme ultraviolet time-resolved photoelectron spectroscopy of adenine, adenosine and adenosine monophosphate in a liquid flat jet

Time-resolved photoelectron spectroscopy using an extreme-ultraviolet (XUV) probe pulse was used to investigate the UV photoinduced dynamics of adenine (Ade), adenosine (Ado), and adenosine-5-monophosphate (AMP) in a liquid water jet. In contrast to previous studies using UV probe pulses, the XUV pulse at 21.7 eV can photoionize all excited states of a molecule, allowing for full relaxation pathways to be addressed after excitation at 4.66 eV. This work was carried out using a gas-dynamic flat liquid jet, resulting in considerably enhanced signal compared to a cylindrical jet. All three species decay on multiple time scales that are assigned based on their decay associated spectra; the fastest decay of ∼100 fs is assigned to ππ* decay to the ground state, while a smaller component with a lifetime of ∼500 fs is attributed to the nπ* state. An additional slower channel in Ade is assigned to the 7H Ade conformer, as seen previously. This work demonstrates the capability of XUV-TRPES to disentangle non-adiabatic dynamics in an aqueous solution in a state-specific manner and represents the first identification of the nπ* state in the relaxation dynamics of adenine and its derivatives.

Disruptive Model That Explains for the Long-Lived Triplet States Observed for 2-Thiocytosine upon UVA Radiation

The possible role of radical species in the formation of the long-lived triplet states observed for 2-thiocytosine upon UV irradiation was theoretically investigated. It is predicted that the radical fragments arising from the homolytic rupture of the NH group of the thiobase can be yielded upon ultraviolet-A radiation. Recombination of the radicals through the most favorable singlet channel yields the lowest-lying tautomer of the 2-thiocytosine (the amino-thiol form) through a barrierless pathway. The rebounding of the radical fragments along the triplet channels that emerge from the attack of the hydrogen to the nitrogen atoms next to the C-S bond leads to stable structures for the amino-thion-N1H and amino-thion-N3H tautomers. These results allow for the rationalization of the near-unity triplet yields observed when this pure light-atom organic molecule is exposed to UV irradiation, without invoking intersystem crossings between the electronic states of different spin-multiplicities. A similar study for cytosine showed that the energy required to induce the homolytic breaking of the N-H bond of the nucleobase is not attainable under UVA radiation. This result is consistent with the experimental fact that no triplet states are observed when this molecule is exposed to that light.

The tautomer-specific excited state dynamics of 2,6-diaminopurine using resonance-enhanced multiphoton ionization and quantum chemical calculations

2,6-Diaminopurine (2,6-dAP) is an alternative nucleobase that potentially played a role in prebiotic chemistry. We studied its excited state dynamics in the gas phase by REMPI, IR-UV hole burning, and ps pump-probe spectroscopy and performed quantum chemical calculations at the SCS-ADC(2) level of theory to interpret the experimental results. We found the 9H tautomer to have a small barrier to ultrafast relaxation via puckering of its 6-membered ring. The 7H tautomer has a larger barrier to reach a conical intersection and also has a sizable triplet yield. These results are discussed relative to other purines, for which 9H tautomerization appears to be more photostable than 7H and homosubstituted purines appear to be less photostable than heterosubstituted or singly substituted purines.

Photodynamics of azaindoles in polar media: the influence of the environment

We have studied the relaxation dynamics of a family of azaindole (AI) structural isomers, 4-, 5-, 6- and 7-AI, by steady-state and time-resolved methods (fs-transient absorption and fluorescence up-conversion), in solvents of different polarity. The measurements in aprotic solvents show distinctive fluorescence yields and excited state lifetimes among the isomers, which are tuned by the polarity of the medium. Guided by simple TD-DFT calculations and based on the behavior observed in the isolated species, it has been possible to address the influence of the environment polarity on the relaxation route. According to the obtained picture, the energy of the nπ* state, which is strongly dependent on the position of the pyridinic nitrogen, controls the rate of the internal conversion channel that accounts for the distinctive photophysical behavior of the isomers. On the other hand, preliminary measurements in protic media (methanol) show a very different photodynamical behavior, in which the anomalous measured fluorescent patterns are very likely the result of reactive channels (proton transfer) triggered by the electronic excitation.

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.

Photophysics of tz Adenine and tz Guanine fluorescent nucleobases embedded into DNA and RNA

UV-VIS photoinduced events of tz A and tz G embedded into DNA and RNA are described by combining the Extended Multi-State Second-Order Perturbation Theory (XMS-CASPT2) and electrostatic embedding molecular mechanics methods (QM/MM). Our results point out that the S1 1 (ππ* La ) state is the bright state in both environments. After the photoexcitation to the S1 1 (ππ* La ) state, the electronic population evolves barrierless towards its minimum, from where the excess of energy can be dissipated by fluorescence. As the minimum energy crossing point structure between the ground and first bright states lies in a high-energy region, the direct internal conversion to the ground state is an unviable mechanism. Other spectroscopic properties (for instance, absorption and Stokes shifts) and comparisons with photochemical properties of canonical nucleobases are also provided.

Cryogenic Ion Spectroscopy of Protonated and Sodiated Methyladenine Derivatives

Ultraviolet photodissociation (UVPD) spectra of protonated 9-methyladenine (H⁺9MA), protonated 7-methyl adenine (H⁺7MA), protonated 3-methyladenine (H⁺3MA), and sodiated 7-methyladenine (Na⁺7MA) near the origin bands of the S₀-S₁ transition were obtained using cryogenic ion spectroscopy. The UV-UV hole burning, infrared (IR) ion-dip, and IR-UV double resonance spectra showed that all the ions were present as single isomers in a cryogenic ion trap. The UVPD spectrum of H⁺9MA exhibited only a broad absorption band, whereas the spectra of H⁺7MA, H⁺3MA, and Na⁺7MA displayed moderately or well-resolved vibronic bands. Potential energy profiles were computed to understand the reason for the different bandwidths of the vibronic bands in the spectra. The broadening of the bands was correlated with the slopes between the Franck-Condon point and the conical intersection between the S₁ and S₀ states in the potential energy profiles, thus reflecting the deactivation rates in the S₁ state.

The impact of different geometrical restrictions on the nonadiabatic photoisomerization of biliverdin chromophores

The photoisomerization mechanism of the chromophore of bacterial biliverdin (BV) phytochromes is explored via nonadiabatic dynamics simulation by using the on-the-fly trajectory surface-hopping method at the semi-empirical OM2/MRCI level. Particularly, the current study focuses on the influence of geometrical constrains on the nonadiabatic photoisomerization dynamics of the BV chromophore. Here a rather simplified approach is employed in the nonadiabatic dynamics to capture the features of geometrical constrains, which adds mechanical restrictions to the specific moieties of the BV chromophore. This simplified method provides a rather quick approach to examine the influence of geometrical restrictions on photoisomerization. As expected, different constrains bring distinctive influences on the photoisomerization mechanism of the BV chromophore, giving either strong or minor modification of both involved reaction channels and excited-state lifetimes after the constrains are added in different ring moieties. These observations not only contribute to the primary understanding of the role of the spatial restriction caused by biological environments in photoinduced dynamics of the BV chromophore, but also provide useful ideas for the artificial regulation of the photoisomerization reaction channels of phytochrome proteins.

Rationalizing the Unexpected Sensitivity in Excited State Lifetimes of Adenine to Tautomerization by Nonadiabatic Molecular Dynamics

The remarkable photostability of canonical nucleobases makes them ideal building blocks for DNA and RNA. Even minor structural changes are expected to lead to drastic alteration of their subpicosecond excited state lifetimes. However, it is interesting to note that while the 9H- and 7H-amino tautomers of adenine possess drastically different lifetimes, 9H- and 7H-keto guanine possess similar excited state lifetimes. With an aim to explain this unexpected difference in sensitivity of lifetimes to tautomerization, we have investigated the excited state relaxation mechanism of UV-excited adenine and guanine tautomers using surface hopping based nonadiabatic molecular dynamics. We find that internal conversion in both guanine tautomers is almost barrierless while both adenine tautomers encounter significant barriers before they can deactivate. Moreover, the major deactivation channel (C2-puckering) in 9H-amino adenine is overall more efficient than the one (C6-puckering) in the 7H-amino form. We trace this difference to the frequent rotation of the amino group which disrupts its conjugation with the heterocyclic ring thereby reducing the strength of nonadiabatic coupling and, hence, delaying internal conversion.

Nonadiabatic Molecular Dynamics with Extended Density Functional Tight-Binding: Application to Nanocrystals and Periodic Solids

In this work, we report a new methodology for nonadiabatic molecular dynamics calculations within the extended tight-binding (xTB) framework. We demonstrate the applicability of the developed approach to finite and periodic systems with thousands of atoms by modeling "hot" electron relaxation dynamics in silicon nanocrystals and electron-hole recombination in both a graphitic carbon nitride monolayer and a titanium-based metal-organic framework (MOF). This work reports the nonadiabatic dynamic simulations in the largest Si nanocrystals studied so far by the xTB framework, with diameters up to 3.5 nm. For silicon nanocrystals, we find a non-monotonic dependence of "hot" electron relaxation rates on the nanocrystal size, in agreement with available experimental reports. We rationalize this relationship by a combination of decreasing nonadiabatic couplings related to system size and the increase of available coherent transfer pathways in systems with higher densities of states. We emphasize the importance of proper treatment of coherences for obtaining such non-monotonic dependences. We characterize the electron-hole recombination dynamics in the graphitic carbon nitride monolayer and the Ti-containing MOF. We demonstrate the importance of spin-adaptation and proper sampling of surface hopping trajectories in modeling such processes. We also assess several trajectory surface hopping schemes and highlight their distinct qualitative behavior in modeling the excited-state dynamics in superexchange-like models depending on how they handle coherences between nearly parallel states.

Excited-state proton transfer in the 2-aminopyridine dimer: A surface hopping study

We present a theoretical investigation of the excited-state intermolecular proton transfer process in the 2-aminopyridine dimer. Prior experimental and theoretical studies on this doubly hydrogen bonded system have attributed an...

Excited State Deactivation Mechanisms of Protonated Adenine: a Theoretical study

Quantum chemical computational method as well as the adiabatic dynamics simulation have been employed to investigate the non-radiative relaxation mechanism of protonated 9H- and 7H-adenine (AH+). We have located three...

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.

Ultrafast Dynamics of the Isolated Adenosine-5'-triphosphate Dianion Probed by Time-Resolved Photoelectron Imaging

The excited state dynamics of the doubly deprotonated dianion of adenosine-5'-triphosphate, [ATP-H₂]^2-, has been spectroscopically explored by time-resolved photoelectron spectroscopy following excitation at 4.66 eV. Time-resolved photoelectron spectra show that two competing processes occur for the initially populated ¹ππ* state. The first is rapid electron emission by tunneling through a repulsive Coulomb barrier as the ¹ππ* state is a resonance. The second is nuclear motion on the ¹ππ* state surface leading to an intermediate that no longer tunnels and subsequently decays by internal conversion to the ground electronic state. The spectral signatures of the features are similar to those observed for other adenine-derivatives, suggesting that this nucleobase is quite insensitive to the nearby negative charges localized on the phosphates, except of course for the appearance of the additional electron tunneling channel, which is open in the dianion.

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.

Ultrafast nonradiative deactivation of photoexcited 8-oxo-hypoxanthine: a nonadiabatic molecular dynamics study

In the scientific endeavor to understand the chemical origins of life, the photochemistry of the smallest life building blocks, nucleobases, has been a constant object of focus and intense research. Here, we report the results of the first theoretical study on the photo-properties of an 8-oxo-hypoxanthine molecule, the chromophore of 8-oxo-inosine, which is relevant to the recently proposed, prebiotically plausible synthetic routes to the formation of purine- and pyrimidine-nucleotides. With ab initio and semi-empirical OM2/MRCI quantum-chemistry calculations, we predict a strong photostability of the 8-oxo-hypoxanthine system and see the origin of this effect in ultrafast nonradiative relaxation through puckering of the 6-membered heterocyclic ring.

An nπ* gated decay mediates excited-state lifetimes of isolated azaindoles

Aiming to serve as a guide to understand the relaxation mechanisms of more complex aza-aromatic compounds, such as purine bases, we have studied the non-radiative channels of a set of azaindole structural isomers: 4-, 5-, 6- and 7-azaindole (AI). The relaxation of the isolated molecules, after excitation at the low energy portion of their spectra, has been tracked by femtosecond time-resolved ionization, and the decay paths have been obtained with MS-CASPT2//TD-DFT calculations. Although the ultrashort measured lifetimes for 5- and 6-AI are in contrast to the long-living excited state found in 7-AI, the calculations describe a common relaxation pathway. Along it, the initially excited ππ* states decay to the ground state through a conical intersection accessed through an nπ* state that functions as a gate state. The work reveals that the position of the nitrogen atoms in the purine ring determines the barrier to access the gate state and therefore, the rate of the non-radiative relaxation.

Solvation Structures and Deactivation Pathways of Luminescent Isothiazole-Derived Nucleobases: tz, tz, and tz

The photophysical relaxation pathways of ^tzA, ^tzG, and ^tzI luminescent nucleobases were investigated with the MS-CASPT2 quantum-chemical method and double-ζ basis sets (cc-pVDZ) in gas and condensed phases (1,4-dioxane and water) with the sequential Monte Carlo/CASPT2 and free energy gradient (FEG) methods. Solvation shell structures, in the ground and excited states, were examined with the pairwise radial distribution function (G(r)) and solute-solvent hydrogen-bond networks. Site-specific hydrogen bonding analysis evidenced relevant changes between both electronic states. The three luminescent nucleobases share a common photophysical pattern, summarized as the lowest-lying ¹(ππ*) bright state that is populated directly after the absorption of radiation and evolves barrierless to the minimum energy structure, from where the excess of energy is released by fluorescence. From the ¹(ππ*)_min region, the conical intersection with the ground state ((ππ*/GS)_CI) is not accessible due to the presence of high energetic barriers. By combining the present results with those reported earlier by us for the pyrimidine fluorescent nucleobases, we present a comprehensive description of the photophysical properties of this important class of new fluorescent nucleosides.

Comparison of Spin-Flip TDDFT-Based Conical Intersection Approaches with XMS-CASPT2

Determining conical intersection geometries is of key importance to understanding the photochemical reactivity of molecules. While many small- to medium-sized molecules can be treated accurately using multireference approaches, larger molecules require a less computationally demanding approach. In this work, minimum energy crossing point conical intersection geometries for a series of molecules have been studied using spin-flip TDDFT (SF-TDDFT), within the Tamm-Dancoff Approximation, both with and without explicit calculation of nonadiabatic coupling terms, and compared with both XMS-CASPT2 and CASSCF calculated geometries. The less computationally demanding algorithms, which do not require explicit calculation of the nonadiabatic coupling terms, generally fare well with the XMS-CASPT2 reference structures, while the relative energetics are only reasonably replicated with the MECP structure as calculated with the BHHLYP functional and full nonadiabatic coupling terms. We also demonstrate that, occasionally, CASSCF structures deviate quantitatively from the XMS-CASPT2 structures, showing the importance of including dynamical correlation.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces

Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

metaFALCON: A Program Package for Automatic Sampling of Conical Intersection Seams Using Multistate Metadynamics

The multistate metadynamics for automatic exploration of conical intersection seams and systematic location of minimum energy crossing points in molecular systems and its implementation into the software package metaFALCON is presented. Based on a locally modified energy gap between two Born-Oppenheimer electronic states as a collective variable, multistate metadynamics trajectories are driven toward an intersection point starting from an arbitrary ground state geometry and are subsequently forced to explore the conical intersection seam landscape. For this purpose, an additional collective variable capable of distinguishing structures within the seam needs to be defined and an additional bias is introduced into the off-diagonal elements of an extended (multistate) electronic Hamiltonian. We demonstrate the performance of the algorithm on the examples of the 1,3-butadiene, benzene, and 9H-adenine molecules, where multiple minimum energy crossing points could be systematically located using the Wiener number or Cremer-Pople parameters as collective variables. Finally, with the example of 9H-adenine, we show that the multistate metadynamics potential can be used to obtain a global picture of a conical intersection seam. Our method can be straightforwardly connected with any ab initio or semiempirical electronic structure theory that provides energies and gradients of the respective electronic states and can serve for systematic elucidation of the role of conical intersections in the photophysics and photochemistry of complex molecular systems, thus complementing nonadiabatic dynamics simulations.

Direct observation of an intramolecular charge transfer state in epigenetic nucleobase N6-methyladenine

N6-Methyladenine (6MeAde), the most abundant internal modification in mRNA, has proved to be an important epigenetic biomarker for gene regulation just like 5-methylcytosine in DNA. Recently, a unique UV-induced response of 6MeAde was reported, which makes it instructive and intriguing to reveal the excited state relaxation mechanism in this methylated adenine and its derivatives. In this work, we investigated 6MeAde and its ribose species N6-methyladenosine (6MeAdo) by using femtosecond time-resolved fluorescence up-conversion (FUC) and broadband transient absorption (TA) spectroscopy. Both 6MeAde and 6MeAdo exhibit a hundreds of femtoseconds lifetime, which originates from the efficient depletion of the ππ* (La) state. A several picoseconds lifetime is also observed and it should be attributed to the ππ* (Lb) state. Surprisingly, dual peak fluorescence emission is observed in 6MeAde and the long wavelength emission is ascribed to an intramolecular charge transfer (ICT) state. The lifetime of this ICT state is determined to be 107 ps. The kinetic isotope effect shows that the ICT state is closely associated with the solute-solvent H-bonding in aqueous solution. In 6MeAdo, the ICT state is apparently quenched and adenine-like excited state dynamics suggests that DNA/RNA containing such modification could still possess excellent photostability under UV irradiation. Our results contain an important insight for understanding excited state properties in epigenetic modified DNA/RNA.

pH Controlled Intersystem Crossing and Singlet Oxygen Generation of 8-Azaadenine in Aqueous Solution

Azabases are intriguing DNA and RNA analogues and have been used as effective antiviral and anticancer medicines. However, photosensitivity of these drugs has also been reported. Here, pH-controlled intersystem crossing (ISC) process of 9H 8-azaadenine (8-AA) in aqueous solution is reported. Broadband transient absorption measurements reveal that the hydrogen atom at N9 position can greatly affect ISC of 8-AA and ISC is more favorable when 8-AA is in its neutral form in aqueous solution. The initial excited ππ* (S₂ ) state evolves through ultrafast internal conversion (IC) (4.2 ps) to the lower-lying nπ* state (S₁ ), which further stands as a door way state for ISC with a time constant of 160 ps. The triplet state has a lifetime of 6.1 μs. On the other hand, deprotonation at N9 position promotes the IC from the ππ* (S₂ ) state to the ground state (S₀ ) and the lifetime of the S₂ state is determined to be 10 ps. The experimental results are further supported by time-dependent density functional theory (TDDFT) calculations. Singlet oxygen generation yield is measured to be 13.8 % for the neutral 8-AA while the deprotonated one exhibit much lower yield (<2 %), implying that this compound could be a potential pH-sensitized photodynamic therapy agent.

Mechanisms of photoreactivity in hydrogen-bonded adenine-H2O complexes

The mechanisms of photoinduced reactions of adenine with water molecules in hydrogen-bonded adenine-water complexes were investigated with ab initio wave-function-based electronic-structure calculations. Two excited-state electron/proton transfer reaction mechanisms have been characterized: H-atom abstraction from water by photoexcited adenine as well as H-atom transfer from photoexcited adenine or the (adenine+H) radical to water. In the water-to-adenine H-atom transfer reaction, an electron from one of the p orbitals of the water molecule fills the hole in the n (π) orbital of the nπ* (ππ*) excited state of adenine, resulting in a charge-separated electronic state. The electronic charge separation is neutralized by the transfer of a proton from the water molecule to adenine, resulting in the (adenine+H)OH biradical in the electronic ground state. In the adenine-to-water H-atom transfer reaction, πσ* states localized at the acidic sites of adenine provide the mechanism for the photoejection of an electron from adenine, which is followed by proton transfer to the hydrogen-bonded water molecule, resulting in the (adenine-H)H3O biradical. The energy profiles of the photoreactions have been computed as relaxed scans with the ADC(2) electronic-structure method. These reactions, which involve the reactivity of adenine with hydrogen-bonded water molecules, compete with the well-established intrinsic excited-state deactivation mechanisms of adenine via ring-puckering or ring-opening conical intersections. By providing additional decay channels, the electron/proton exchange reactions with water can account for the observed significantly shortened excited-state lifetime of adenine in aqueous environments. These findings indicate that adenine possibly was not only a photostabilizer at the beginning of life, but also a primordial photocatalyst for water splitting.

Tracking the origin of photostability in purine nucleobases: the photophysics of 2-oxopurine

Molding purine PES through functionalization: whilst purine C₂-substitution maintains the features of the spectroscopic PES of the heterocycle, C₆-functionalization reshapes its topography leading to photostable systems.

Highly efficient surface hopping dynamics using a linear vibronic coupling model

We report an implementation of the linear vibronic coupling (LVC) model within the surface hopping dynamics approach and present utilities for parameterizing this model in a blackbox fashion. This results in an extremely efficient method to obtain qualitative and even semi-quantitative information about the photodynamical behavior of a molecule, and provides a new route toward benchmarking the results of surface hopping computations. The merits and applicability of the method are demonstrated in a number of applications. First, the method is applied to the SO2 molecule showing that it is possible to compute its absorption spectrum beyond the Condon approximation, and that all the main features and timescales of previous on-the-fly dynamics simulations of intersystem crossing are reproduced while reducing the computational effort by three orders of magnitude. The dynamics results are benchmarked against exact wavepacket propagations on the same LVC potentials and against a variation of the electronic structure level. Four additional test cases are presented to exemplify the broader applicability of the model. The photodynamics of the isomeric adenine and 2-aminopurine molecules are studied and it is shown that the LVC model correctly predicts ultrafast decay in the former and an extended excited-state lifetime in the latter. Futhermore, the method correctly predicts ultrafast intersystem crossing in the modified nucleobase 2-thiocytosine and its absence in 5-azacytosine while it fails to describe the ultrafast internal conversion to the ground state in the latter.

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.

Sub-50 fs excited state dynamics of 6-chloroguanine upon deep ultraviolet excitation

The photophysical properties of natural nucleobases and their respective nucleotides are ascribed to the sub-picosecond lifetime of their first singlet states in the UV-B region (260-350 nm). Electronic transitions of the ππ* type, which are stronger than those in the UV-B region, lie at the red edge of the UV-C range (100-260 nm) in all isolated nucleobases. The lowest energetic excited states in the UV-B region of nucleobases have been investigated using a plethora of experimental and theoretical methods in gas and solution phases. The sub-picosecond lifetime of these molecules is not a general attribute of all nucleobases but specific to the five primary nucleobases and a few xanthine and methylated derivatives. To determine the overall UV photostability, we aim to understand the effect of more energetic photons lying in the UV-C region on nucleobases. To determine the UV-C initiated photophysics of a nucleobase system, we chose a halogen substituted purine, 6-chloroguanine (6-ClG), that we had investigated previously using resonance Raman spectroscopy. We have performed quantitative measurements of the resonance Raman cross-section across the Bb absorption band (210-230 nm) and constructed the Raman excitation profiles. We modeled the excitation profiles using Lee and Heller's time-dependent theory of resonance Raman intensities to extract the initial excited state dynamics of 6-ClG within 30-50 fs after photoexcitation. We found that imidazole and pyrimidine rings of 6-ClG undergo expansion and contraction, respectively, following photoexcitation to the Bb state. The amount of distortions of the excited state structure from that of the ground state structure is reflected by the total internal reorganization energy that is determined at 112 cm(-1). The contribution of the inertial component of the solvent response towards the total reorganization energy was obtained at 1220 cm(-1). In addition, our simulation also yields an instantaneous response of the first solvation shell within an ultrafast timescale of less than 30 fs following photoexcitation.

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.

Solvent effects on the excited state characteristics of adenine–thymine base pairs

A systematic analysis of the excited state characteristics of the DNA base pair adenine–thymine in stacked and Watson–Crick hydrogen bonded configurations has been carried out in this study.

Ultrafast structural dynamics of photoexcited adenine

Ultraviolet Resonance Raman (UVRR) spectroscopy derives distinct electronic properties of adenine in the L_a (260 nm) and B_b (210 nm) excited states.

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.

Vibrationally resolved optical spectra and ultrafast electronic relaxation dynamics of diamantane

We present theoretical simulations of the vibrationally resolved photoabsorption and photoemission spectra of diamantane combined with nonadiabatic dynamics simulations in order to identify the state responsible for the measured photoluminescence of diamantane and to determine the mechanism and the time-scales of the electronic state relaxation. Diamantane is a prototype representative of the diamondoid class of hydrocarbons which have recently gained significant interest due to their unique electronic properties. This molecule is characterised by an almost dark first excited state, which therefore cannot be directly excited. Moreover, the calculated vertical transition from the geometrically relaxed first excited state to the ground state also bears no intensity. However, recent experiments suggest that the observed photoluminescence originates from the lowest excited state. We have performed spectral simulations in the frame of the Herzberg-Teller approximation for vibronic transitions, which goes beyond the Franck-Condon approximation of constant transition dipole moments and takes into account their linear dependence on the geometrical deformations. In this way, the available experimental spectrum could be fully reproduced, resolving the issue about the origin of the photoluminescence. Moreover, the photoemission from the first excited state also implies that ultrafast nonradiative processes have to take place after the initial excitation of the bright electronic states. We have determined the mechanism and time-scales of these relaxation processes by performing nonadiabatic dynamics simulations in the manifold of s- and p-type Rydberg excited states. The simulations demonstrate that the lowest excited electronic state of diamantane gains significant population from higher-lying states already after several hundreds of femtoseconds. Thus, our dynamics simulations combined with spectra calculated using the Herzberg-Teller approximation allow us to fully explain the observed photoabsorption and photoemission properties of diamantane.

Multiple Decay Mechanisms and 2D-UV Spectroscopic Fingerprints of Singlet Excited Solvated Adenine-Uracil Monophosphate

The decay channels of singlet excited adenine uracil monophosphate (ApU) in water are studied with CASPT2//CASSCF:MM potential energy calculations and simulation of the 2D-UV spectroscopic fingerprints with the aim of elucidating the role of the different electronic states of the stacked conformer in the excited state dynamics. The adenine (1) La state can decay without a barrier to a conical intersection with the ground state. In contrast, the adenine (1) Lb and uracil S(U) states have minima that are separated from the intersections by sizeable barriers. Depending on the backbone conformation, the CT state can undergo inter-base hydrogen transfer and decay to the ground state through a conical intersection, or it can yield a long-lived minimum stabilized by a hydrogen bond between the two ribose rings. This suggests that the (1) Lb , S(U) and CT states of the stacked conformer may all contribute to the experimental lifetimes of 18 and 240 ps. We have also simulated the time evolution of the 2D-UV spectra and provide the specific fingerprint of each species in a recommended probe window between 25 000 and 38 000 cm(-1) in which decongested, clearly distinguishable spectra can be obtained. This is expected to allow the mechanistic scenarios to be discerned in the near future with the help of the corresponding experiments. Our results reveal the complexity of the photophysics of the relatively small ApU system, and the potential of 2D-UV spectroscopy to disentangle the photophysics of multichromophoric systems.