Toward Consistent Predictions of Core/Valence Ionization Potentials and Valence Excitation Energies by MRSF-TDDFT

Optimizing exchange-correlation functionals for both core/valence ionization potentials (cIPs/vIPs) and valence excitation energies (VEEs) at the same time in the framework of MRSF-TDDFT is self-contradictory. To overcome the challenge, within the previous "adaptive exact exchange" or double-tuning strategy on Coulomb-attenuating XC functionals (CAM), a new XC functional specifically for cIPs and vIPs was first developed by enhancing exact exchange to both short- and long-range regions. The resulting DTCAM-XI functional achieved remarkably high accuracy in its predictions with errors of less than half eV. An additional concept of "valence attenuation", where the amount of exact exchange for the frontier orbital regions is selectively suppressed, was introduced to consistently predict both VEEs and IPs at the same time. The second functional, DTCAM-XIV, exhibits consistent overall prediction accuracy at ∼0.64 eV. By preferentially optimizing VEEs within the same "valence attenuation" concept, a third functional, DTCAM-VAEE, was obtained, which exhibits improved performance as compared to that of the previous DTCAM-VEE and DTCAM-AEE in the prediction of VEEs, making it an attractive alternative to BH&HLYP. As the combination of "adaptive exchange" and "valence attenuation" is operative, it would be exciting to explore its potential with a more tunable framework in the future.

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.

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

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

Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory

It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.218 and 0.327 eV, respectively]. On the other hand, DTCAM-AEE faithfully reproduced the qualitative aspects of conical intersections (CIs) of trans-butadiene and thymine and the nonadiabatic molecular dynamics (NAMD) simulations on thymine. The latter functional also remarkably exhibited the exact 1/R asymptotic behavior of the charge-transfer state of an ethylene-tetrafluoroethylene dimer and the accurate potential energy surfaces (PESs) along the two torsional angles of retinal protonated Schiff base model with six double bonds (rPSB6). Overall, DTCAM-AEE generally performs well, as its MAE (0.237) and IQR (0.41 eV) are much improved as compared to BH&HLYP. The current idea can also be applied to other XC functionals as well as other variants of linear response theories, opening a new way of developing XC functionals.

Significance of Energy Conservation in Coupled-Trajectory Approaches to Nonadiabatic Dynamics

Through approximating electron-nuclear correlation terms in the exact factorization approach, trajectory-based methods have been derived and successfully applied to the dynamics of a variety of light-induced molecular processes, capturing quantum (de)coherence effects rigorously. These terms account for the coupling among the trajectories, recovering the nonlocal nature of quantum nuclear dynamics that is completely overlooked in traditional independent-trajectory algorithms. Nevertheless, some of the approximations introduced in the derivation of some of these methods do not conserve the total energy. We analyze energy conservation in the coupled-trajectory mixed quantum-classical (CTMQC) algorithm and explore the performance of a modified algorithm, CTMQC-E, where some of the terms are redefined to restore energy conservation. A set of molecular models is used as a test, namely, 2-cis-penta-2,4-dienimium cation, bis(methylene) adamantyl radical cation, butatriene cation, uracil radical cation, and neutral pyrazine.

Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory: Multireference Advantages with the Practicality of Linear Response Theory

The density functional theory (DFT) and linear response (LR) time-dependent (TD)-DFT are of the utmost importance for routine computations. However, the single reference formulation of DFT suffers in the description of open-shell singlet systems such as diradicals and bond-breaking. LR-TDDFT, on the other hand, finds difficulties in the modeling of conical intersections, doubly excited states, and core-level excitations. In this Perspective, we demonstrate that many of these limitations can be overcome by recently developed mixed-reference (MR) spin-flip (SF)-TDDFT, providing an alternative yet accurate route for such challenging situations. Empowered by the practicality of the LR formalism, it is anticipated that MRSF-TDDFT can become one of the major workhorses for general routine tasks.

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

The exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron-ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron-nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

Ultrafast Photo-Ion Probing of the Relaxation Dynamics in 2-Thiouracil

In this work, we investigate the relaxation processes of 2-thiouracil after UV photoexcitation to the S₂ state through the use of ultrafast, single-colour, pump-probe UV/UV spectroscopy. We place focus on investigating the appearance and subsequent decay signals of ionized fragments. We complement this with VUV-induced dissociative photoionisation studies collected at a synchrotron, allowing us to better understand and assign the ionisation channels involved in the appearance of the fragments. We find that all fragments appear when single photons with energy > 11 eV are used in the VUV experiments and hence appear through 3+ photon-order processes when 266 nm light is used. We also observe three major decays for the fragment ions: a sub-autocorrelation decay (i.e., sub-370 fs), a secondary ultrafast decay on the order of 300-400 fs, and a long decay on the order of 220 to 400 ps (all fragment dependent). These decays agree well with the previously established S₂ → S₁ → Triplet → Ground decay process. Results from the VUV study also suggest that some of the fragments may be created by dynamics occurring in the excited cationic state.

Photochemistry of Thymine in Solution and DNA Revealed by an Electrostatic Embedding QM/MM Combined with Mixed-Reference Spin-Flip TDDFT

The photochemistry of nucleobases, important for their role as building blocks of DNA, is largely affected by the electrostatic environment in which they are soaked. For example, despite the numerous studies of thymine in solution and DNA, there is still a debate on the photochemical deactivation pathways after UV absorption. Many theoretical models are oversimplified due to the lack of computationally accurate and efficient electronic structure methodologies that capture excited state electron correlation effects when nucleobases are embedded in large electrostatic media. Here, we combine mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) with electrostatic embedding QM/MM using electrostatic potential fittingfitted (ESPF) atomic charges, as a strategy to accurately and efficiently describe the electronic structure of chromophores polarized by an electrostatic medium. In particular, we develop analytic expressions for the energy and gradient of MRSF/MM based on the ESPF coupling using atom-centered grids and total charge conservation. We apply this methodology to the study of solvation effects on thymine photochemistry in water and thymine dimers in DNA. In the former, the combination of trajectory surface hopping (TSH) nonadiabatic molecular dynamics (NAMD) with MRSF/MM remarkably revealed accelerated deactivation decay pathways, which is consistent with the experimental decay time of ∼400 fs. The enhanced hopping rate can be explained by the preferential stabilization of corresponding conical interactions due to their increased dipole moments. Structurally, it is a consequence of characteristic methyl puckered geometries near the conical intersection region. For the thymine dimer in B-DNA, we found new photochemical pathways through conical intersections that could explain the formation of cyclobutadiene dimers and 6-4 photoproducts.

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 (L_a and L_b) 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 Photo-ion Probing of the Ring-Opening Process in Trans-Stilbene Oxide

The ultrafast photo‐induced ring opening of the oxirane derivative trans‐stilbene oxide has been studied through the use of ultrafast UV/UV pump‐probe spectroscopy by using photo‐ion detection. Single‐ and multiphoton probe paths and final states were identified through comparisons between UV power studies and synchrotron‐based vacuum ultraviolet (VUV) single‐photon ionization studies. Three major time‐dependent features of the parent ion (sub‐450 fs decay, (1.5±0.2) ps, and >100 ps) were observed. These decays are discussed in conjunction with the primary ring‐opening mechanism of stilbene oxide, which occurs through C−C dissociation in the oxirane ring. The appearance of fragments relating to the masses of dehydrogenated diphenylmethane (167 amu) and dehydrogenated methylbenzene (90 amu) were also investigated. The appearance of the 167 amu fragment could suggest an alternative ultrafast ring‐opening pathway via the dissociation of one of the C−O bonds within the oxirane ring.

Impact of the Dynamic Electron Correlation on the Unusually Long Excited-State Lifetime of Thymine

Non-radiative relaxation of the photoexcited thymine in the gas phase shows an unusually long excited-state lifetime, and, over the years, a number of models, i.e., S₁-trapping, S₂-trapping, and S₁&S₂-trapping, have been put forward to explain its mechanism. Here, we investigate this mechanism using non-adiabatic molecular dynamics (NAMD) simulations in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. We show that the previously predicted S₂-trapping model was due to an artifact caused by an insufficient account of the dynamic electron correlation. The current work supports the S₁-trapping mechanism with two lifetimes, τ₁ = 30 ± 1 fs and τ₂ = 6.1 ± 0.035 ps, quantitatively consistent with the recent time-resolved experiments. Upon excitation to the S₂ (ππ*) state, thymine undergoes an ultrafast (ca. 30 fs) S₂→S₁ internal conversion and resides around the minimum on the S₁ (n_Oπ*) surface, slowly decaying to the ground state (ca. 6.1 ps). While the S₂→S₁ internal conversion is mediated by fast bond length alternation distortion, the subsequent S₁→S₀ occurs through several conical intersections, involving a slow puckering motion.

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.

ESR, STESR, DFT, and MD Study of the Dynamical Structure of Cucurbituril[7]-Spin Probe Guest-Host Complexes

We study the molecular dynamics and structures of the guest-host complexes of cucurbituril, CB[7], with spin probes through the conventional electron spin resonance (ESR), saturation transfer ESR (STESR), density functional theory (DFT), and molecular dynamics (MD) computations. Protonated TEMPOamine (I), a derivative of TEMPO having a positive charge and an octyl group on the quaternary nitrogen atom (II), and the neutral spin-labeled indole (III) are used as guests. To eliminate the overall complex rotation, the solutions of complexes in a solid CB[7] matrix were prepared. Resultantly, for all of the spin probes, the combined study of the conventional ESR and STESR spectra indicates the librational character of the rotational motion within the CB[7] cavity as opposed to the diffusional rotation over the whole solid angle. The kinetic accessibilities of the reporter NO groups to the paramagnetic complexes in aqueous solutions, determined by Heisenberg exchange broadening of the ESR spectra, together with the environment polarities from the hyperfine interaction values, as well as DFT computation results and MD simulations, were used to estimate the spin probe location relative to CB[7]. Utilizing the concept of the aqueous clusters surrounding the spin probes and CB[7] molecules and MD simulations has allowed the application of DFT to estimate the aqueous environment effects on the complexation energy and spatial structure of the guest-host complexes.

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

Photochemical pathways in nucleobases measured with an X-ray FEL

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

Multistate hybrid time-dependent density functional theory with surface hopping accurately captures ultrafast thymine photodeactivation

We report an efficient analytical implementation of first-order nonadiabatic derivative couplings between arbitrary Born–Oppenheimer states in the hybrid time-dependent density functional theory (TDDFT) framework using atom-centered basis functions.

Conical-intersection dynamics and ground-state chemistry probed by extreme-ultraviolet time-resolved photoelectron spectroscopy

Time-resolved photoelectron spectroscopy (TRPES) is a useful approach to elucidate the coupled electronic-nuclear quantum dynamics underlying chemical processes, but has remained limited by the use of low photon energies. Here, we demonstrate the general advantages of XUV-TRPES through an application to NO₂, one of the simplest species displaying the complexity of a non-adiabatic photochemical process. The high photon energy enables ionization from the entire geometrical configuration space, giving access to the true dynamics of the system. Specifically, the technique reveals dynamics through a conical intersection, large-amplitude motion and photodissociation in the electronic ground state. XUV-TRPES simultaneously projects the excited-state wave packet onto many final states, offering a multi-dimensional view of the coupled electronic and nuclear dynamics. Our interpretations are supported by ab initio wavepacket calculations on new global potential-energy surfaces. The presented results contribute to establish XUV-TRPES as a powerful technique providing a complete picture of ultrafast chemical dynamics from photoexcitation to the final products.

Probing the excited state relaxation dynamics of pyrimidine nucleosides in chloroform solution

Ultrafast transient electronic and vibrational absorption spectroscopy (TEAS and TVAS) of 2'-deoxy-cytidine (dC) and 2'-deoxy-thymidine (dT) dissolved in chloroform examines their excited-state dynamics and the recovery of ground electronic state molecules following absorption of ultraviolet light. The chloroform serves as a weakly interacting solvent, allowing comparisons to be drawn with prior experimental studies of the photodynamics of these nucleosides in the gas phase and in polar solvents such as water. The pyrimidine base nucleosides have some propensity to dimerize in aprotic solvents, but the monomer photochemistry can be resolved clearly and is the focus of this study. UV absorption at a wavelength of 260 nm excites a ¹ππ* ← S₀ transition, but prompt crossing of a significant fraction (50% in dC, 17% in dT) of the ¹ππ* population into a nearby ¹nπ* state is too fast for the experiments to resolve. The remaining flux on the ¹ππ* state leaves the vertical Franck-Condon region and encounters a conical intersection with the ground electronic state of ethylenic twist character. In dC, the ¹ππ* state decays to the ground state with a time constant of 1.1 ± 0.1 ps. The lifetime of the ¹nπ* state is much longer in the canonical forms of both molecules: recovery of the ground state population from these states occurs with time constants of 18.6 ± 1.1 ps in amino-oxo dC and ∼114 ps in dT, indicating potential energy barriers to the ¹nπ*/S₀ conical intersections. The small fraction of the imino-oxo tautomer of dC present in solution has a longer-lived ¹nπ* state with a lifetime for ground state recovery of 193 ± 55 ps. No evidence is found for photo-induced tautomerization of amino-oxo dC to the imino-oxo form, or for population of low lying triplet states of this nucleoside. In contrast, ∼8% of the UV-excited dT molecules access the long-lived T₁ (³ππ*) state through the ¹nπ* state. The primary influence of the solvent appears to be the degree to which it destabilizes the states of ¹nπ* character, with consequences for the lifetimes of these states as well as the triplet state yields.

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

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

New Insights into the State Trapping of UV-Excited Thymine

After UV excitation, gas phase thymine returns to a ground state in 5 to 7 ps, showing multiple time constants. There is no consensus on the assignment of these processes, with a dispute between models claiming that thymine is trapped either in the first (S₁) or in the second (S₂) excited states. In the present study, a nonadiabatic dynamics simulation of thymine is performed on the basis of ADC(2) surfaces, to understand the role of dynamic electron correlation on the deactivation pathways. The results show that trapping in S₂ is strongly reduced in comparison to previous simulations considering only non-dynamic electron correlation on CASSCF surfaces. The reason for the difference is traced back to the energetic cost for formation of a CO π bond in S₂.

Multiconfigurational Second-Order Perturbation Theory with Frozen Natural Orbitals Extended to the Treatment of Photochemical Problems

A new flavor of the frozen natural orbital complete active space second-order perturbation theory method (FNO-CASPT2, Aquilante et al., J. Chem. Phys. 131, 034113) is proposed herein. In this new implementation, the virtual space in Cholesky decomposition-based CASPT2 computations (CD-CASPT2) is truncated by excluding those orbitals that contribute the least toward preserving a predefined value of the trace of an approximate density matrix, as that represents a measure of the amount of dynamic correlation retained in the model. In this way, the amount of correlation included is practically constant at all nuclear arrangements, thus allowing for the computation of smooth electronic states surfaces and energy gradients-essential requirements for theoretical studies in photochemistry. The method has been benchmarked for a series of relevant biochromophores for which large speed-ups have been recorded while retaining the accuracy achieved in the corresponding CD-CASPT2 calculations. Both vertical excitation energies and gradient calculations have been carried out to establish general guidelines as to how much correlation needs to be retained in the calculation for the results to be consistent with the CD-CASPT2 findings. Our results feature errors within a tenth of an eV for the most difficult cases and have been validated to be used for gradient computations where an up to 3-fold speed-up is observed depending on the size of the system and the basis set employed.

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.

Characterizing the dark state in thymine and uracil by double resonant spectroscopy and quantum computation

Double resonant spectroscopy characterizes both grounds state and dark excited state of uracil and thymine.

Ultrafast non-radiative decay of gas-phase nucleosides

De-excitation of DNA nucleosides on picosecond timescales was measured and found to be twice as fast as the equivalent nucleobases.

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.

Protonation effect on the electronic properties of 2-pyridone monomer, dimer and its water clusters: a theoretical study

The CC2 (second order approximate coupled cluster method) has been applied to investigate protonation effect on electronic transition energies of 2-pyridone (2PY), 2-pyridone dimer, and micro-solvated 2-pyridone (0-2 water molecules). The PE profiles of protonated 2-pyridone (2PYH(+)) as well as monohydrated 2PYH(+) at the different electronic states have been investigated. The (1)πσ∗ state in protonated species (2PYH(+)) is a barrier free and dissociative state along the O-H stretching coordinate. In this reaction coordinate, the lowest lying (1)πσ∗ predissociates the bound S1((1)ππ∗) state, connecting the latter to a conical intersection with the S0 state. These conical intersections lead the (1)ππ∗ state to proceed as predissociative state and finally direct the excited system to the ground state. Furthermore, in presence of water molecule, the (1)πσ∗ state still remains dissociative but the conical intersection between (1)πσ∗ and ground state disappears. In addition, according to the CC2 calculation results, it has been predicted that protonation significantly blue shifts the S1-S0 electronic transition of monomer, dimer, and microhydrated 2-pyridone.

Time-resolved photoelectron imaging of the isolated deprotonated nucleotides

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

Photophysics and photochemistry of thymine deoxy-dinucleotide in water: a PCM/TD-DFT quantum mechanical study

We here report a fully quantum mechanical study of the main photochemical and photophysical decay routes in aqueous solution of thymine deoxy-dinucleotide (TpT(-) and TpTNa) and of its analogue locked in C3-endo puckering, characterizing five different representative backbone conformers and discussing the chemical physical effects modulating the yield of the different photoproducts. Our approach is based on time-dependent DFT calculations, using the last generation M052X functional, whereas solvent effects are included by means of the polarizable continuum model. Especially when at least one of the sugars adopts C3-endo puckering, a barrierless path on the bright ππ* excitons leads to the S(1)/S(0) crossing region corresponding to the formation of cyclobutane pyrimidine dimer. Charge transfer excited states involving the transfer of an electron from the 5' Thy toward the 3' Thy are involved in the formation of the oxetane intermediate in the path leading to 6-4 pyrimidine pyrimidinone adducts. A non-negligible energy barrier is associated with this latter pathway, which is possible only when one of the two nucleotides adopts C2-endo puckering. Monomer-like decay pathways, involving ππ* or nπ* excited states localized on a single base, are shown to be operative also for loosely stacked bases.

Specific photodynamics in thymine clusters: the role of hydrogen bonding

A photoionization detected IR study of thymine and 1-methylthymine monohydrates and of their homodimers was carried out to shed some light on the structure of the thymine clusters whose complex photodynamics has recently been the subject of great interest. Under supersonic jet conditions, thymine forms doubly H-bonded cyclic clusters with water or another base preferentially via its N1-H group and the adjacent carbonyl group. This hydrate is of no biological relevance since the N1-H group is the sugar binding site in thymidine. On the other hand, 1-methylthymine forms the donor H-bonds only via the N3-H group. Hence, properties of the N1-H and the N3-H bound clusters of thymine can be studied using thymine and 1-methylthymine molecules, respectively. No biologically relevant conformations of the dimers and hydrates of thymine, contrary to those of 1-methylthymine, are observed under supersonic jet conditions. Thymine homodimer, which extensively fragments upon UV ionization by formation of a protonated monomer, exhibits two N1-H···O═C2 hydrogen bonds. The photodynamics of hydrated thymines is found to be extremely sensitive to the hydration site: ranging from an ultrafast relaxation in less than 100 fs up to formation of a dark state with the lifetime on the microsecond time scale.