The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications

Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.

Excited state hydrogen or proton transfer pathways in microsolvated n-cyanoindole fluorescent probes

The sensitivity of the fluorescence properties of n-cyanoindole (n-CNI) fluorescent probes to the microenvironment makes them potential reporters of protein conformation and hydration. The fluorescence intensity of 5-CNI, 6-CNI, and 7-CNI is quenched when exposed to water solvent whereas substitution on position 4 of indoles dramatically increases it. A potential mechanism for this sensitivity to water may be similar to that found in indole. The fluorescence of indole is found to be quenched when interacting with water and ammonia solvent molecules via radiationless decay through an S1 (πσ*)/S0 conical intersection caused by excited state proton or hydrogen transfer to the solvent molecules. In this study we examine this fluorescence quenching mechanism along the N-H bond stretch of n-CNI probes using water cluster models and quantum mechanical calculations of the excited states. We find that n-CNI-(H2O)1-2 clusters form cyclic or non-cyclic structures via hydrogen bonds which lead to different photochemical reaction paths that can potentially quench the fluorescence by undergoing internal conversion from S1 to S0. However, the existence of a high energy barrier along the potential energy surface of the S1 state in most cases prevents this from occurring. We show that substitution on position 4 leads to the highest energy barrier that prevents the fluorophore from accessing these non-radiative channels, in agreement with its high intensity. We also find that the energy barrier in the S1 state of non-cyclic 5-CNI-(H2O)1-2 excited complexes decreases as the number of water molecules increases, which suggests great sensitivity of the fluorescence quenching on the aqueous environment.

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

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

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

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

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.

Photoinduced water-chromophore electron transfer causes formation of guanosine photodamage

UV-induced photolysis of aqueous guanine nucleosides produces 8-oxo-guanine and Fapy-guanine, which can induce various types of cellular malfunction. The mechanistic rationale underlying photodestructive processes of guanine nucleosides is still largely obscure. Here, we employ accurate quantum chemical calculations and demonstrate that an excited-state non-bonding interaction of guanosine and a water molecule facilitates the electron-driven proton transfer process from water to the chromophore fragment. This subsequently allows for the formation of a crucial intermediate, namely guanosine photohydrate. Further (photo)chemical reactions of this intermediate lead to the known products of guanine photodamage.

Excited state dynamics of 2'-deoxyisoguanosine and isoguanosine in aqueous solution

Photostability is thought to be an inherent property of nucleobases required to survive the extreme ultraviolet radiation conditions of the prebiotic era. Previous studies have shown that absorption of ultraviolet radiation by the canonical nucleosides results in ultrafast internal conversion to the ground state, demonstrating that these nucleosides efficiently dissipate the excess electronic energy to the environment. In recent years, studies on the photophysical and photochemical properties of nucleobase derivatives have revealed that chemical substitution influences the electronic relaxation pathways of purine and pyrimidine nucleobases. It has been suggested that amino or carbonyl substitution at the C6 position could increase the photostability of the purine derivatives more than the substitution at the C2 position. This investigation aims to elucidate the excited state dynamics of 2'-deoxyisoguanosine (dIsoGuo) and isoguanosine (IsoGuo) in aqueous solution at pH 7.4 and 1.4, which contain an amino group at the C6 position and a carbonyl group at the C2 position of the purine chromophore. The study of these derivatives is performed using absorption and emission spectroscopies, broadband transient absorption spectroscopy, and density functional and time-dependent density functional levels of theory. It is shown that the primary relaxation mechanism of dIsoGuo and IsoGuo involves nonradiative decay pathways, where the population decays from the S₁(ππ*) state through internal conversion to the ground state via two relaxation pathways with lifetimes of hundreds of femtoseconds and less than 2 ps, making these purine nucleosides photostable in aqueous solution.

Excited State Dynamics of Methylated Guanosine Derivatives Revealed by Femtosecond Time-resolved Spectroscopy

Methylated DNA/RNA nucleobases are important epigenetic marks in living species and play an important role for targeted therapies. Moreover, they could bring significant changes to the photo-stability of nucleic acid, leading these sites become mutational hotspots for disease such as skin cancer. While a number of studies have demonstrated the relationship between excited state dynamics and the biological function of methylated cytosine in DNA, investigations aimed at unraveling the excited state dynamics of methylated guanosine in RNA have been largely overlooked. In this work, influence of methylation on the excited state dynamics of guanosine is studied by using femtosecond time-resolved spectroscopy. Our results suggest that the effect of methyl substitution on the photophysical properties of guanosine is position sensitive. N1-methylguanosine shows very similar excited state dynamics as that in guanosine, while almost one order of magnitude longer lifetime of the La state is observed in N2, N2-dimethylguanosine. Notably, N7-methylation can lead to a new minimum on the La state, which shows a two orders of magnitude longer excited state lifetime compared with guanosine. These findings not only help understanding excited state dynamics of methylated guanosines, but also lay the foundation for further studying DNA/RNA strands incorporated with these bases.

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

Quantum and semiclassical dynamical studies of nonadiabatic processes in solution: achievements and perspectives

We concisely review the main methodological approaches to model nonadiabatic dynamics in isotropic solutions and their applications. Three general classes of models are identified as the most used to include solvent effects in the simulations. The first model describes the solvent as a set of harmonic collective modes coupled to the solute degrees of freedom, and the second as a continuum, while the third explicitly includes solvent molecules in the calculations. The issues related to the use of these models in semiclassical and quantum dynamical simulations are discussed, as well as the main limitations and perspectives of each approach.

Excitation energy transport in DNA modelled by multi-chromophoric field-induced surface hopping

Revealing the extended excited state lifetime due to excitation energy transport in DNA by multi-chromophoric field-induced surface-hopping (McFISH).

MS-CASPT2 Studies on the Photophysics of Selenium-Substituted Guanine Nucleobase

The MS-CASPT2 method has been employed to optimize minimum-energy structures of 6-selenoguanine (6SeGua) and related two- and three-state intersection structures in and between the lowest five electronic states, i.e., S₂(¹ππ*), S₁(¹ nπ*), T₂(³ nπ*), T₁(³ππ*), and S₀. In combination with MS-CASPT2 calculated linearly interpolated internal coordinate paths, the photophysical mechanism of 6SeGua has been proposed. The initially populated S₂(¹ππ*) state decays to either S₁(¹ nπ*) or T₂(³ nπ*) states through a three-state S₂/S₁/T₂ intersection point. The large S₂/T₂ spin-orbit coupling of 435 cm^-1, according to the classical El-Sayed rule, benefits the S₂ → T₂ intersystem crossing process. The S₁(¹ nπ*) state that stems from the S₂ → S₁ internal conversion process at the S₂/S₁/T₂ intersection point can further jump to the T₂(³ nπ*) state through the S₁ → T₂ intersystem crossing process. This process does not comply with the El-Sayed rule, but it is still related to a comparatively large spin-orbit coupling of 39 cm^-1 and is expected to occur relatively fast. Finally, the T₂(³ nπ*) state, which is populated from the above S₂ → T₂ and S₁ → T₂ intersystem crossing processes, decays to the T₁(³ππ*) state via an internal conversion process. Because there is merely a small energy barrier of 0.11 eV separating the T₁(³ππ*) minimum and an energetically allowed two-state T₁/S₀ intersection point, the T₁(³ππ*) state still can decay to the S₀ state quickly, which is also enhanced by a large T₁/S₀ spin-orbit coupling of 252 cm^-1. Our proposed mechanism explains experimentally observed ultrafast intersystem crossing processes in 6SeGua and its 835-fold acceleration of the T₁ state decay to the S₀ state compared with 6tGua. Finally, we have found that the ground-state electronic structure of 6SeGua has more apparent multireference character.

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.

Ultrafast action chemistry in slow motion: atomistic description of the excitation and fluorescence processes in an archetypal fluorescent protein

We report quantum mechanical/molecular mechanical non-adiabatic molecular dynamics simulations on the electronically excited state of green fluorescent protein mutant S65T/H148D. We examine the driving force of the ultrafast (τ < 50 fs) excited-state proton transfer unleashed by absorption in the A band at 415 nm and propose an atomistic description of the two dynamical regimes experimentally observed [Stoner Ma et al., J. Am. Chem. Soc., 2008, 130, 1227]. These regimes are explained in terms of two sets of successive dynamical events: first the proton transfers quickly from the chromophore to the acceptor Asp148. Thereafter, on a slower time scale, there are geometrical changes in the cavity of the chromophore that involve the distance between the chromophore and Asp148, the planarity of the excited-state chromophore, and the distance between the chromophore and Tyr145. We find two different non-radiative relaxation channels that are operative for structures in the reactant region and that can explain the mismatch between the decay of the emission of A* and the rise of the emission of I*, as well as the temperature dependence of the non-radiative decay rate.

Photochemical relaxation pathways of S6-methylthioinosine and O6-methylguanosine in solution

S6-Methylthioinosine and O6-methylguanosine are byproducts resulting from the enzymatic reactions of sulfur-substituted prodrugs in cells and from the interaction of alkylating agents with cellular DNA, respectively. Their photochemistry has not been investigated, and it is currently unknown whether light absorption by these byproducts may pose any threat to the cell. In this contribution, their photoinduced processes upon absorption of UVB radiation are reported using broadband transient absorption spectroscopy. Plausible electronic relaxation mechanisms are proposed for both biological molecules, which are supported by steady-state absorption and emission measurements, and by singlet and triplet vertical excitation energies performed on a large subset of ground-state optimized conformational isomers in solution. The results are compared to the body of knowledge gathered in the scientific literature about the light-induced processes in the sulfur-substituted and canonical purine monomers. In particular, it is shown that S6-methylation decreases the rate to populate the lowest-energy triplet state and blueshifts the ground-state absorption spectrum compared to those for the sulfur-substituted prodrugs and for the 6-thioguanosine metabolite. Similarly, O6-methylation decreases the rate of internal conversion to the ground state observed in the guanine monomers by more than 10-fold in acetonitrile and 40-fold in aqueous solution, while it redshifts the ground-state absorption spectrum. Collectively, this investigation provides relevant new insights about the relationship between structural modifications of the purine chromophore and the electronic relaxation mechanisms in this important group of biological molecules.

Ultrafast unidirectional chiral rotation in the Z–E photoisomerization of two azoheteroarene photoswitches

Based on a large number of trajectories starting from the Z-isomer, for both azoheteroarenes, more than 99% of the trajectories decay through conical intersections with the same helicities as their initial geometries.

Excited-State Dynamics in O6-Methylguanosine: Impact of O6-Methylation on the Relaxation Mechanism of Guanine Monomers

Absorption of ultraviolet radiation by DNA bases results in ultrafast internal conversion to the ground state, which minimizes photodamage. However, exogenous and endogenous alkylating agents present in the cellular environment can methylate the nucleobases in DNA. In particular, methylation of guanosine at the O⁶ position in DNA leads to the formation of the O⁶-methylguanosine adduct, which may alter the photostability of DNA. This contribution demonstrates that O⁶-methylation of guanosine red shifts its ground-state absorption spectrum and slows down the rate of internal conversion to the ground state by ∼40-fold in aqueous solution. The 40-fold decrease in the rate of excited-state decay increases the probability of photodamage within cellular DNA. It is proposed that the longer decay lifetime corresponds to relaxation of the excited-state population in O⁶-methylguanosine along a C6-puckered reaction coordinate in the ¹ππ*(L_a) potential energy surface that runs parallel to an ultrafast internal conversion pathway along a C2-puckered coordinate.

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.

Ultraviolet Light Makes dGMP Floppy: Femtosecond Stimulated Raman Spectroscopy of 2'-Deoxyguanosine 5'-Monophosphate

The ultrafast dynamics of 2'-deoxyguanosine 5'-monophosphate after excitation with ultraviolet light has been studied with femtosecond transient absorption (TA) and femtosecond stimulated Raman spectroscopy (FSRS). TA kinetics and transient anisotropy spectra reveal a rapid relaxation from the Franck-Condon region, producing an extremely red-shifted stimulated emission band at ∼440 nm that is formed after 200 fs and subsequent relaxation for 0.8-1.5 ps, consistent with prior studies. Viscosity dependence shows that the initial relaxation, before 0.5 ps, is the same in water or viscous glycerol/water mixtures, but after 0.5 ps the dynamics significantly slow down in a viscous solution. This indicates that large amplitude structural changes occur after 0.5 ps following photoexcitation. FSRS obtained with both 480 and 600 nm Raman pump pulses observe very broad Raman peaks at 509 and 1530 cm^-1, as well as a narrower peak at 1179 cm^-1. All of the Raman peaks decay with 0.7-1.3 ps time constants. The 1530 cm^-1 peak also shows an increasing inhomogeneous linewidth over the first 0.3 ps. Our TA and FSRS data are consistent with a structurally inhomogeneous population in the S₁ (L_a) state and, in particular, with previous theoretical models in which out-of-plane distortion at C2 and the amine move the molecule toward a conical intersection with the ground state. These FSRS data are the first to directly observe the structural inhomogeneity imparted upon the excited-state population by the broad, flat potential energy surface of the S₁ (L_a) state.

Linear photoabsorption spectra and vertical excitation energies of microsolvated DNA nucleobases in aqueous solution

Employing density functional theory (DFT) and time-dependent density functional theory (TDDFT) calculations in combination with the semiclassical nuclear ensemble method, we have simulated the photoabsorption spectra of the four canonical DNA nucleobases in aqueous solution. In order to model the effects of solvation, for each nucleobase, a number of solvating water molecules were explicitly included in the simulations, and additionally, the bulk solvent was represented by a continuous polarizable medium. We find that the effect of the solvation shell in general is significant, and its inclusion improves the realism of the spectral simulations. The involvement of lone electron pairs in the hydrogen bonding with the solvating water molecules has the effect of systematically increasing the energies of vertical excitation into the [Formula: see text]-type states. Apart from a systematic blue shift of around [Formula: see text][Formula: see text]eV observed in the absorption peaks, the calculated photoabsorption spectra reproduce the measured ones with good accuracy. The photoabsorption spectra are dominated by excited states with [Formula: see text] and partial [Formula: see text] character. No low-energy charge transfer states are observed with the use of the CAM-B3LYP and M06-2X functionals.

Restoring electronic coherence/decoherence for a trajectory-based nonadiabatic molecular dynamics

By utilizing the time-independent semiclassical phase integral, we obtained modified coupled time-dependent Schrödinger equations that restore coherences and induce decoherences within original simple trajectory-based nonadiabatic molecular dynamic algorithms. Nonadiabatic transition probabilities simulated from both Tully’s fewest switches and semiclassical Ehrenfest algorithms follow exact quantum electronic oscillations and amplitudes for three out of the four well-known model systems. Within the present theory, nonadiabatic transitions estimated from statistical ensemble of trajectories accurately follow those of the modified electronic wave functions. The present theory can be immediately applied to the molecular dynamic simulations of photochemical and photophysical processes involving electronic excited states.

Gas-Phase Femtosecond Particle Spectroscopy: A Bottom-Up Approach to Nucleotide Dynamics

We summarize how gas-phase ultrafast charged-particle spectroscopy has been used to provide an understanding of the photophysics of DNA building blocks. We focus on adenine and discuss how, following UV excitation, specific interactions determine the fates of its excited states. The dynamics can be probed using a systematic bottom-up approach that provides control over these interactions and that allows ever-larger complexes to be studied. Starting from a chromophore in adenine, the excited state decay mechanisms of adenine and chemically substituted or clustered adenine are considered and then extended to adenosine mono-, di-, and trinucleotides. We show that the gas-phase approach can offer exquisite insight into the dynamics observed in aqueous solution, but we also highlight stark differences. An outlook is provided that discusses some of the most promising developments in this bottom-up approach.

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.

Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Theory, Implementation, and Parameters

Semiempirical orthogonalization-corrected methods (OM1, OM2, and OM3) go beyond the standard MNDO model by explicitly including additional interactions into the Fock matrix in an approximate manner (Pauli repulsion, penetration effects, and core–valence interactions), which yields systematic improvements both for ground-state and excited-state properties. In this Article, we describe the underlying theoretical formalism of the OMx methods and their implementation in full detail, and we report all relevant OMx parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. For a standard set of mostly organic molecules commonly used in semiempirical method development, the OMx results are found to be superior to those from standard MNDO-type methods. Parametrized Grimme-type dispersion corrections can be added to OM2 and OM3 energies to provide a realistic treatment of noncovalent interaction energies, as demonstrated for the complexes in the S22 and S66×8 test sets.

Excited-state deactivation in 8-oxo-deoxyguanosine: comparison between anionic and neutral forms

Ab initio explorations of excited-state potential-energy surfaces show that a radiationless deactivation mechanism via intramolecular excited-state proton transfer is available in neutral 8-oxo-deoxyguanosine, whereas it is not available in the anionic form.

Generalized trajectory surface-hopping method for internal conversion and intersystem crossing

Trajectory-based fewest-switches surface-hopping (FSSH) dynamics simulations have become a popular and reliable theoretical tool to simulate nonadiabatic photophysical and photochemical processes. Most available FSSH methods model internal conversion. We present a generalized trajectory surface-hopping (GTSH) method for simulating both internal conversion and intersystem crossing processes on an equal footing. We consider hops between adiabatic eigenstates of the non-relativistic electronic Hamiltonian (pure spin states), which is appropriate for sufficiently small spin-orbit coupling. This choice allows us to make maximum use of existing electronic structure programs and to minimize the changes to available implementations of the traditional FSSH method. The GTSH method is formulated within the quantum mechanics (QM)/molecular mechanics framework, but can of course also be applied at the pure QM level. The algorithm implemented in the GTSH code is specified step by step. As an initial GTSH application, we report simulations of the nonadiabatic processes in the lowest four electronic states (S0, S1, T1, and T2) of acrolein both in vacuo and in acetonitrile solution, in which the acrolein molecule is treated at the ab initio complete-active-space self-consistent-field level. These dynamics simulations provide detailed mechanistic insight by identifying and characterizing two nonadiabatic routes to the lowest triplet state, namely, direct S1 → T1 hopping as major pathway and sequential S1 → T2 → T1 hopping as minor pathway, with the T2 state acting as a relay state. They illustrate the potential of the GTSH approach to explore photoinduced processes in complex systems, in which intersystem crossing plays an important role.

Constructing polyatomic potential energy surfaces by interpolating diabatic Hamiltonian matrices with demonstration on green fluorescent protein chromophore

Simulating molecular dynamics directly on quantum chemically obtained potential energy surfaces is generally time consuming. The cost becomes overwhelming especially when excited state dynamics is aimed with multiple electronic states. The interpolated potential has been suggested as a remedy for the cost issue in various simulation settings ranging from fast gas phase reactions of small molecules to relatively slow condensed phase dynamics with complex surrounding. Here, we present a scheme for interpolating multiple electronic surfaces of a relatively large molecule, with an intention of applying it to studying nonadiabatic behaviors. The scheme starts with adiabatic potential information and its diabatic transformation, both of which can be readily obtained, in principle, with quantum chemical calculations. The adiabatic energies and their derivatives on each interpolation center are combined with the derivative coupling vectors to generate the corresponding diabatic Hamiltonian and its derivatives, and they are subsequently adopted in producing a globally defined diabatic Hamiltonian function. As a demonstration, we employ the scheme to build an interpolated Hamiltonian of a relatively large chromophore, para-hydroxybenzylidene imidazolinone, in reference to its all-atom analytical surface model. We show that the interpolation is indeed reliable enough to reproduce important features of the reference surface model, such as its adiabatic energies and derivative couplings. In addition, nonadiabatic surface hopping simulations with interpolation yield population transfer dynamics that is well in accord with the result generated with the reference analytic surface. With these, we conclude by suggesting that the interpolation of diabatic Hamiltonians will be applicable for studying nonadiabatic behaviors of sizeable molecules.

Computational modeling of photoexcitation in DNA single and double strands

The photoexcitation of DNA strands triggers extremely complex photoinduced processes, which cannot be understood solely on the basis of the behavior of the nucleobase building blocks. Decisive factors in DNA oligomers and polymers include collective electronic effects, excitonic coupling, hydrogen-bonding interactions, local steric hindrance, charge transfer, and environmental and solvent effects. This chapter surveys recent theoretical and computational efforts to model real-world excited-state DNA strands using a variety of established and emerging theoretical methods. One central issue is the role of localized vs delocalized excitations and the extent to which they determine the nature and the temporal evolution of the initial photoexcitation in DNA strands.

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.

Excited-state dynamics of guanosine in aqueous solution revealed by time-resolved photoelectron spectroscopy: experiment and theory

Time-resolved photoelectron spectroscopy is performed on aqueous guanosine solution to study its excited-state relaxation dynamics.

Solvent effects on the ultrafast nonradiative deactivation mechanisms of thymine in aqueous solution: excited-state QM/MM molecular dynamics simulations

On-the-fly excited-state quantum mechanics/molecular mechanics molecular dynamics (QM/MM-MD) simulations of thymine in aqueous solution are performed to investigate the role of solvent water molecules on the nonradiative deactivation process. The complete active space second-order perturbation theory (CASPT2) method is employed for a thymine molecule as the QM part in order to provide a reliable description of the excited-state potential energies. It is found that, in addition to the previously reported deactivation pathway involving the twisting of the C-C double bond in the pyrimidine ring, another efficient deactivation pathway leading to conical intersections that accompanies the out-of-plane displacement of the carbonyl group is observed in aqueous solution. Decay through this pathway is not observed in the gas phase simulations, and our analysis indicates that the hydrogen bonds with solvent water molecules play a key role in stabilizing the potential energies of thymine in this additional decay pathway.

The solvent effect and identification of a weakly emissive state in nonradiative dynamics of guanine nucleosides and nucleotides--a combined femtosecond broadband time-resolved fluorescence and transient absorption study

A combined method of femtosecond broadband time-resolved fluorescence (fs-TRF) and transient absorption (fs-TA) was employed to investigate the excited state dynamics of 2'-deoxyguanosine (dG) and 2'-deoxyguanosine 5'-monophosphate (dGMP). Comparative fs-TRF and fs-TA measurements were conducted on dG and dGMP in neutral water, deuterated water, and methanol with excitation wavelengths of 245, 267 and 285 nm. Very similar results were observed with dG and dGMP. The data provide compelling evidence for the co-existence of two nonradiative pathways. One is the generally recognized Laππ* mediated channel, the other involves an unprecedented weakly emissive state which plays a significant role in the overall deactivation processes. The Laππ* channel features biphasic dynamics with time constants (τ1/τ2) of ~0.2/0.8 ps in water and ~0.25/1.0 ps in methanol. The biphasic decay arises due to a partial transfer with τ1 of the Laππ* population to the newly identified state followed by conversion in τ2 of the remaining Laππ* molecules into the electronic ground state. The channel mediated by the weakly emissive species shows solvent-dependent dynamics with time constants (τ3) of ~2.0 ps in water, ~2.3 ps in deuterated water, and ~4.1 ps in methanol. The species features absorption at UV wavelengths (~300-400 nm) and exhibits deeply red-shifted fluorescence (λmax ~ 520 nm) with polarization direction varied markedly from that of the Laππ* but close to the Lbππ*. This species acts as an effective quenching state to the radiative decay of the brightly emissive Laππ* and Lbππ*. It sets in promptly (<~50 fs) after the photoexcitation and is further populated through nonadiabatic coupling with the Laππ*. The overall involvement of this state is facilitated with excitation at high energy and is favoured in methanol over water. According to the spectral character and the solvent effect in particular the kinetic isotope effect, the species is tentatively associated to the πσ* state with charge transfer (CT) character which is considered to be preferentially stabilized by hydrogen-bonding between the guanine amino and surrounding solvent molecules. The result of this study leads to a dramatically different picture of guanine deactivation. It demonstrates a crucial role of the solvent in shaping the nonradiative dynamics of guanine nucleosides and nucleotides. The data presented are important for understanding the detailed photophysics of not only the monomeric guanine but also DNA assemblies that contain guanine in base pairs or have a guanine tetrad as the structural motif.

A QM/MM MD insight into photodynamics of hypoxanthine: distinct nonadiabatic decay behaviors between keto-N7H and keto-N9H tautomers in aqueous solution

QM/MM MD simulations reveal different H-bonding networks around DNA base analogues, and the π-electron H-bond in the solvated keto-N7H may facilitate its S₁ → S₀ nonadiabatic decay.

Time-resolved photoelectron imaging of the isolated deprotonated nucleotides

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

Ultrafast photo-initiated molecular quantum dynamics in the DNA dinucleotide d(ApG) revealed by broadband transient absorption spectroscopy

The ultrafast photo-initiated quantum dynamics of the adenine-guanine dinucleotide d(ApG) in aqueous solution (pH 7) has been studied by femtosecond time-resolved spectroscopy after excitation at lambda = 260 nm. The results reveal a hierarchy of processes on time scales from tau < 100 fs to tau > 100 ps. Characteristic spectro-temporal signatures are observed indicating the transformation of the molecules in the electronic relaxation from the photo-excited state to a long-lived exciplex. In particular, broadband UV/VIS excited-state absorption (ESA) measurements detected a distinctive absorption by the excited dinucleotide around lambda = 335 nm, approximately 0.5 eV to the blue compared to the maximum of the broad and unstructured ESA spectrum after excitation of an equimolar mixture of the mononucleotides dAMP and dGMP. A similar feature has been identified as signature of the excimer in the dynamics of the adenine dinucleotide d(ApA). The lifetime of the d(ApG) exciplex was found to be tau = 124 +/- 4 ps both from the ESA decay time and from the ground-state recovery time, far longer than the sub-picosecond lifetimes of excited dAMP or dGMP. Fluorescence-time profiles measured by the up-conversion technique indicate that the exciplex state is reached around approximately 6 ps after excitation. Very weak residual fluorescence at longer times red-shifted to the emission from the photo-excited state shows that the exciplex is almost optically dark, but still has enough oscillator strength to give rise to the dual fluorescence of the dinucleotide in the static fluorescence spectrum.

Ultrafast electronic deactivation dynamics of the inosine dimer--a model case for H-bonded purine bases

The structural properties and ultrafast electronic deactivation dynamics of the inosine dimer in CHCl3 have been investigated by two-dimensional (1)H NMR and static FTIR spectroscopy and by femtosecond time-resolved transient absorption spectroscopy, respectively. The (1)H NMR and IR spectra show the formation of a well-defined, symmetric dimer with an association equilibrium constant of KI·I = 690 ± 100 M(-1). The excited-state dynamics after photoexcitation at λpump = 260 nm monitored by ultrafast absorption spectroscopy show great similarity with those of the monomer inosine in an aqueous solution and are governed by a decay time of τ = 90 ± 10 fs, which is one of the shortest electronic lifetimes of all nucleobases and nucleobase dimers studied so far. On the basis of these observations, the inosine dimer is expected to follow a similar relaxation pathway as the monomer, involving an out-of-plane deformation of the six-membered ring. The importance of the C(2) position for the electronic deactivation of hypoxanthine and guanine is discussed. The obtained well-determined structure and straightforward dynamics qualify the inosine dimer as an excellent reference case for more complicated systems such as the G·G dimer and the G·C and A·T Watson-Crick pairs.