Photodynamics Simulations of Thymine: Relaxation into the First Excited Singlet State

UV wavelength-dependent photoionization quantum yields for the dark 1nπ* state of aqueous thymidine

Despite the important role of the dark 1nπ* state in the photostability of thymidine in aqueous solution, no detailed ultraviolet (UV) wavelength-dependent investigation of the 1nπ* quantum yield (QY) in aqueous thymidine has been experimentally performed. Here, we investigate the wavelength-dependent photoemission spectra of aqueous thymidine from 266.7 to 240 nm using liquid-microjet photoelectron spectroscopy. Two observed ionization channels are assigned to resonant ionizations from 1ππ* to the cationic ground state D0 (π-1) and 1nπ* to the cationic excited state D1 (n-1). The weak 1nπ* → D1 ionization channel appears due to ultrafast 1ππ* → 1nπ* internal conversion within the pulse duration of ∼180 fs. The obtained 1nπ* quantum yields exhibit a strong wavelength dependence, ranging from 0 to 0.27 ± 0.01, suggesting a hitherto uncharacterized 1nπ* feature. The corresponding vertical ionization energies (VIEs) of D0 and D1 of aqueous thymidine are experimentally determined to be 8.47 ± 0.12 eV and 9.22 ± 0.29 eV, respectively. Our UV wavelength-dependent QYs might indicate that different structural critical points to connect the multidimensional 1ππ*/1nπ* conical intersection seam onto the multidimensional potential energy surface of the 1ππ* state might exist and determine the relaxation processes of aqueous thymidine upon UV excitation.

A Theoretical Exploration of the Photoinduced Breaking Mechanism of the Glycosidic Bond in Thymine Nucleotide

DNA glycosidic bond cleavage may induce cancer under the ultraviolet (UV) effect. Yet, the mechanism of glycosidic bond cleavage remains unclear and requires more detailed clarification. Herein, quantum chemical studies on its photoinduced mechanism are performed using a 5'-thymidine monophosphate (5'-dTMPH) model. In this study, four possible paths were examined to study the glycosidic bond cleavage. The results showed that, upon excitation, the electronic transition from the π bonding to π antibonding orbitals of the thymine ring leads to the damage of the thymine ring. Afterwards, the glycosidic bond is cleaved. At first, the doublet ground state (GS) path of glycosidic bond cleavage widely studied by other groups is caused by free electron generated by photoirradiation, with a kinetically feasible energy barrier of ~23 kcal/mol. Additionally, then, the other three paths were proposed that also might cause the glycosidic bond cleavage. The first one is the doublet excited state (ES) path, triggered by free electron along with UV excitation, which can result in a very-high-energy barrier ~49 kcal/mol that is kinetically unfavorable. The second one is the singlet ES path, induced by direct UV excitation, which assumes DNA is directly excited by UV light, which features a very low-energy barrier ~16 kcal/mol that is favored in kinetics. The third one is the triplet ES path, from the singlet state via intersystem crossing (ISC), which refers to a feasible ~27 kcal/mol energy barrier. This study emphasizes the pivotal role of the DNA glycosidic bond cleavage by our proposed direct UV excitation (especially singlet ES path) in addition to the authorized indirect free-electron-induced path, which should provide essential insights to future mechanistic comprehension and novel anti-cancer drug design.

Understanding the Photoprocesses in Biological Systems: Need for Accurate Multireference Treatment

Light-matter interaction is crucial to life itself and revolves around many of the central processes in biology. The need for understanding these photochemical and photophysical processes cannot be overemphasized. Interaction of light with biological systems starts with the absorption of light and subsequent phenomena that occur in the excited states of the system. However, excited states are typically difficult to understand within the mean field approximation of quantum chemical methods. Therefore, suitable multireference methods and methodologies have been developed to understand these phenomena. In this Perspective, we will describe a few methods and methodologies suitable for these descriptions and discuss some persisting difficulties.

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.

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.

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.

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.

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.

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

Observation of Ultrafast Intersystem Crossing in Thymine by Extreme Ultraviolet Time-Resolved Photoelectron Spectroscopy

We studied the photoinduced ultrafast relaxation dynamics of the nucleobase thymine using gas-phase time-resolved photoelectron spectroscopy. By employing extreme ultraviolet pulses from high harmonic generation for photoionization, we substantially extend our spectral observation window with respect to previous studies. This enables us to follow relaxation of the excited state population all the way to low-lying electronic states including the ground state. In thymine, we observe relaxation from the optically bright ¹ππ* state of thymine to a dark ¹nπ* state within 80 ± 30 fs. The ¹nπ* state relaxes further within 3.5 ± 0.3 ps to a low-lying electronic state. By comparison with quantum chemical simulations, we can unambiguously assign its spectroscopic signature to the ³ππ* state. Hence, our study draws a comprehensive picture of the relaxation mechanism of thymine including ultrafast intersystem crossing to the triplet manifold.

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.

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.

Insights into the deactivation of 5-bromouracil after ultraviolet excitation

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

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.

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

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

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

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.

Monitoring Ultrafast Chemical Dynamics by Time-Domain X-ray Photo- and Auger-Electron Spectroscopy

The directed flow of charge and energy is at the heart of all chemical processes. Extraordinary efforts are underway to monitor and understand the concerted motion of electrons and nuclei with ever increasing spatial and temporal sensitivity. The element specificity, chemical sensitivity, and temporal resolution of ultrafast X-ray spectroscopy techniques hold great promise to provide new insight into the fundamental interactions underlying chemical dynamics in systems ranging from isolated molecules to application-like devices. Here, we focus on the potential of ultrafast X-ray spectroscopy techniques based on the detection of photo- and Auger electrons to provide new fundamental insight into photochemical processes of systems with various degrees of complexity. Isolated nucleobases provide an excellent testing ground for our most fundamental understanding of intramolecular coupling between electrons and nuclei beyond the traditionally applied Born-Oppenheimer approximation. Ultrafast electronic relaxation dynamics enabled by the breakdown of this approximation is the major component of the nucleobase photoprotection mechanisms. Transient X-ray induced Auger electron spectroscopy on photoexcited thymine molecules provides atomic-site specific details of the extremely efficient coupling that converts potentially bond changing ultraviolet photon energy into benign heat. In particular, the time-dependent spectral shift of a specific Auger band is sensitive to the length of a single bond within the molecule. The X-ray induced Auger transients show evidence for an electronic transition out of the initially excited state within only ∼200 fs in contrast to theoretically predicted picosecond population trapping behind a reaction barrier. Photoinduced charge transfer dynamics between transition metal complexes and semiconductor nanostructures are of central importance for many emerging energy and climate relevant technologies. Numerous demonstrations of photovoltaic and photocatalytic activity have been performed based on the combination of strong light absorption in dye molecules with charge separation and transport in adjacent semiconductor nanostructures. However, a fundamental understanding of the enabling and limiting dynamics on critical atomic length- and time scales is often still lacking. Femtosecond time-resolved X-ray photoelectron spectroscopy is employed to gain a better understanding of a short-lived intermediate that may be linked to the unexpectedly limited performance of ZnO based dye-sensitized solar cells by delaying the generation of free charge carriers. The transient spectra strongly suggest that photoexcited dye molecules attached to ZnO nanocrystals inject their charges into the substrate within less than 1 ps but the electrons are then temporarily trapped at the surface of the semiconductor in direct vicinity of the injecting molecules. The experiments are extended to monitor the electronic response of the semiconductor substrate to the collective injection from a monolayer of dye molecules and the subsequent electron-ion recombination dynamics. The results indicate some qualitative similarities but quantitative differences between the recombination dynamics at molecule-semiconductor interfaces and previously studied bulk-surface electron-hole recombination dynamics in photoexcited semiconductors.

Excited-state relaxation of hydrated thymine and thymidine measured by liquid-jet photoelectron spectroscopy: experiment and simulation

Time-resolved photoelectron spectroscopy is performed on thymine and thymidine in aqueous solution to study the excited-state relaxation dynamics of these molecules. We find two contributions with sub-ps lifetimes in line with recent excited-state QM/MM molecular dynamics simulations (J. Chem. Phys. 2013, 139, 214304). The temporal evolution of ionization energies for the excited ππ* state along the QM/MM molecular dynamics trajectories were calculated and are compatible with experimental results, where the two contributions correspond to the relaxation paths in the ππ* state involving different conical intersections with the ground state. Theoretical calculations also show that ionization from the nπ* state is possible at the given photon energies, but we have not found any experimental indication for signal from the nπ* state. In contrast to currently accepted relaxation mechanisms, we suggest that the nπ* state is not involved in the relaxation process of thymine in aqueous solution.

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 intersystem crossing dynamics in uracil unravelled by ab initio molecular dynamics

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

Intersystem Crossing Enables 4-Thiothymidine to Act as a Photosensitizer in Photodynamic Therapy: An Ab Initio QM/MM Study

Motivated by its potential use as a photosensitizer in photodynamic therapy, we report the first ab initio quantum mechanics/molecular mechanics (QM/MM) study of 4-thiothymidine in aqueous solution. The core chromophore 4-thiothymine was described using the multiconfigurational CASSCF and CASPT2 QM methods, while the ribose and the solvent water molecules were treated at the MM level (CHARMM and TIP3P, respectively). The minima of the five lowest electronic states (S0, S1, S2, T1, and T2) and six minimum-energy intersections were fully optimized at the QM(CASSCF)/MM level, and their energies were further refined by single-point QM(CASPT2)/MM and CASPT2 calculations. The relevant spin-orbit couplings were also computed. We find that (1) there are three efficient photophysical pathways that account for the experimentally observed ultrafast formation of the lowest triplet state with a quantum yield of nearly unity, (2) the striking qualitative differences in the photophysical behavior of 4-thiothymine and thymine originate from the different electronic structure of their S1 states, and (3) environmental effects play an important role. The present QM/MM calculations provide mechanistic insight that may guide the design of improved photosensitizers for photodynamic therapy.

Computational reference data for the photochemistry of cyclobutane pyrimidine dimers

The cis-syn cyclobutane pyrimidine dimer is one of the major classes of carcinogenic UV-induced DNA photoproducts. In this work, diverse high-level quantum-chemical methods were used to determine the spectroscopic properties of neutral (singlet and triplet) and charged (cation and anion) species of thymine dimers. Maps of potential energy, charge distribution, electron affinity, and ionization potential of the thymidine dimers were computed along the two dimerization coordinates for neutral and charged species, as well as for the singlet excited state. This set of data aims at providing consistent results computed with the same methods as for photodamage and repair. Based on these results, several different photo-, heat-, and charge-induced mechanisms of dimerization and repair are characterized and discussed. Additionally, a new stable dimer with methylmethylidene-hexahydropyrimidine structure was found in the S0 state.

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.

Photoinduced processes in nucleic acids

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

State-specific heavy-atom effect on intersystem crossing processes in 2-thiothymine: a potential photodynamic therapy photosensitizer

Thiothymidine has a potential application as a photosensitizer in cancer photodynamic therapy (PDT). As the chromophore of thiothymidine, 2-thiothymine exhibits ultrahigh quantum yield of intersystem crossing to the lowest triplet state T(1) (ca. 100%), which contrasts with the excited-state behavior of the natural thymine that dissipates excess electronic energy via ultrafast internal conversion to the ground state. In this work, we employed high-level complete-active space self-consistent field and its second-order perturbation methods to explore the photophysical mechanism of a 2-thiothymine model. We have optimized the minimum energy structures in the low-lying seven electronic states, as well as ten intersection points. On the basis of the computed potential energy profiles and spin-orbit couplings, we proposed three competitive, efficient nonadiabatic pathways to the lowest triplet state T(1) from the initially populated singlet state S(2). The suggested mechanistic scenario explains well the recent experimental phenomena. The origin responsible for the distinct photophysical behaviors between thymine and 2-thiothymine is ascribed to the heavy-atom effect, which is significantly enhanced in the latter. Additionally, this heavy-atom effect is found to be state-specific, which could in principle be used to tune the photophysics of 2-thiothymine. The present high-level electronic structure calculations also contribute to understand the working mechanism of thiothymidine in PDT.

Comparison of vertical and adiabatic harmonic approaches for the calculation of the vibrational structure of electronic spectra

The calculation of the vibrational structure associated to electronic spectra in large molecules requires a Taylor expansion of the initial and final state potential energy surface (PES) around some reference nuclear structure. Vertical (V) and adiabatic (A) approaches expand the final state PES around the initial-state (V) or final-state (A) equilibrium structure. Simplest models only take into account displacements of initial- and final-state minima, intermediate ones also allow for difference in frequencies and more accurate models introduce the Dushinsky effect through the computation of the Hessians of both the initial and final state. In this contribution we summarize and compare the mathematical expressions of the complete hierarchy of V and A harmonic models and we implement them in a numerical code, presenting a detailed comparison of their performance on a number of prototypical systems. We also address non-Condon effects through linear expansions of the transition dipole as a function of nuclear coordinates (Herzberg-Teller effect) and compare the results of expansions around initial and final state equilibrium geometries. By a throughout analysis of our results we highlight a number of general trends in the relative performance of the models that can provide hints for their proper choice. Moreover we show that A and V models including final state PES Hessian outperform the simpler ones and that discrepancies in their predictions are diagnostic for failure of harmonic approximation and/or of Born-Oppenheimer approximation (existence of remarkable geometry-dependent mixing of electronic states).

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.

ONIOM approach for non-adiabatic on-the-fly molecular dynamics demonstrated for the backbone controlled Dewar valence isomerization

Non-adiabatic on-the-fly molecular dynamics (NA-O-MD) simulations require the electronic wavefunction, energy gradients, and derivative coupling vectors in every timestep. Thus, they are commonly restricted to the excited state dynamics of molecules with up to ≈20 atoms. We discuss an approximation that combines the ONIOM(QM:QM) method with NA-O-MD simulations to allow calculations for larger molecules. As a proof of principle we present the excited state dynamics of a (6-4)-lesion containing dinucleotide (63 atoms), and especially the importance to include the confinement effects of the DNA backbone. The method is able to include electron correlation on a high level of theory and offers an attractive alternative to QM:MM approaches for moderate sized systems with unknown force fields.