On the Intrinsic Population of the Lowest Triplet State of Thymine

Accurate Spin-Orbit Coupling by Relativistic Mixed-Reference Spin-Flip-TDDFT

Relativistic mixed-reference spin-flip (MRSF)-TDDFT is developed considering the spin-orbit coupling (SOC) within the mean-field approximation. The resulting SOC-MRSF faithfully reproduces the experiments with very high accuracy, which is also consistent with the values by four-component (4c) relativistic CASSCF and 4c-CASPT2 in the spin-orbit-energy splitting calculations of the C, Si, and Ge atoms. Even for the fifth-row element Sn, the SOC-MRSF yielded accurate splittings (∼ 3 % error). In the SOC calculations of the molecular 4-thiothymine with a third-row element, SOC-MRSF values are in excellent agreement with those of the SO-GMC-QDPT2 level, regardless of geometries and exchange-correlation functionals. The same SOC-MRSF predicted the anticipated chance of S₁ (nπ*) → T₁ (ππ*) intersystem crossing, even in thymine with only second-row elements. With its accuracy and practicality, thus, SOC-MRSF is a promising electronic structure protocol in challenging situations such as nonadiabatic molecular dynamics (NAMD) incorporating both internal conversions and intersystem crossings in large systems.

The Behavior of Triplet Thymine in a Model B-DNA Strand. Energetics and Spin Density Localization Revealed by ab initio Molecular Dynamics Simulations†

Among the naturally occurring nucleobases, thymine presents the lowest triplet state, hence it represents a hotspot for energy transfer and photosensitization. In turn, the population of the triplet state may lead to thymine dimerization and hence to the production of toxic DNA lesions and has been the subject of intensive theoretical and experimental investigations. Relying on QM/MM molecular dynamics simulations, we have sought to situate the energy of the lowest triplet state of thymine embedded in a B-DNA environment. The energy gap varies between 305 and 329 kJ mol^-1 when a single thymine is treated at the quantum chemistry level, depending on its position in the model double-stranded 16-bp oligonucleotide. The energy of triplet state decreases up to 300 kJ mol^-1 , due to polarization effects, when we consider coupled stacked nucleobases up to the inclusion of four nucleobases. Our value lies in good agreement with the energy inferred experimentally by Miranda and coworkers (270 kJ mol^-1 ), and our theoretical exploration opens the door to investigations toward other more complex and biologically relevant environments, such as thymines embedded in nucleosome core particles. Our investigations also provide a reference for further studies using semi-empirical approaches such as density functional-based tight-binding, allowing to further rationalize sequence effects.

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.

Stabilization of the Triplet Biradical Intermediate of 5-Methylcytosine Enhances Cyclobutane Pyrimidine Dimer (CPD) Formation in DNA

Cyclobutane pyrimidine dimer (CPD) is a photoproduct formed by two stacked pyrimidine bases through a cycloaddition reaction upon irradiation. Owing to its close association with skin cancer, the mechanism of CPD formation has been studied thoroughly. Among many aspects of CPD, its formation involving 5-methylcytosine (5mC) has been of special interest because the CPD yield is known to increase with C5-methylation of cytosine. In this work, high-level quantum mechanics/molecular mechanics (QM/MM) calculations are used to examine a previously experimentally detected pathway for CPD formation in hetero (thymine-cytosine and thymine-5mC) dipyrimidines, which is facilitated through intersystem crossing in thymine and formation of a triplet biradical intermediate. A DNA duplex model system containing a core sequence TmCG or TCG is used. The stabilization of a radical center in the biradical intermediate by the methyl group of 5mC can lead to increased CPD yield in TmCG compared with its non-methylated counterpart, TCG, thereby suggesting the existence of a new pathway of CPD formation enhanced by 5mC.

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.

Excited-State Dynamics in the RNA Nucleotide Uridine 5'-Monophosphate Investigated Using Femtosecond Broadband Transient Absorption Spectroscopy

Damage to RNA from ultraviolet radiation induces chemical modifications to the nucleobases. Unraveling the excited states involved in these reactions is essential; however, investigations aimed at understanding the electronic-energy relaxation pathways of the RNA nucleotide uridine 5'-monophosphate (UMP) have not received enough attention. In this Letter, the excited-state dynamics of UMP is investigated in aqueous solution. Excitation at 267 nm results in a trifurcation event that leads to the simultaneous population of the vibrationally excited ground state, a long-lived ¹nπ* state, and a receiver triplet state within 200 fs. The receiver state internally converts to the long-lived ³ππ* state in an ultrafast time scale. The results elucidate the electronic relaxation pathways and clarify earlier transient absorption experiments performed for uracil derivatives in solution. This mechanistic information is important because long-lived nπ* and ππ* excited states of both singlet and triplet multiplicities are thought to lead to the formation of harmful photoproducts.

UV-induced long-lived decays in solvated pyrimidine nucleosides resolved at the MS-CASPT2/MM level

The most relevant 'dark' electronic excited states in DNA/RNA pyrimidine nucleosides are mapped in water employing hybrid MS-CASPT2/MM optimisations with explicit solvation and including the sugar. Conical intersections (CIs) between initially accessed bright ¹ππ* and the lowest energy dark ¹nπ* excited states, involving the lone pair localised on the oxygen and/or nitrogen atoms are characterised. They are found in the vicinities of the Franck-Condon (FC) region and are shown to facilitate non-adiabatic population transfer. The excited state population of the ¹n_Oπ* state, localised in the carbonyl moiety on all pyrimidine nucleosides, is predicted to rapidly evolve to its minimum, displaying non-negligible potential energy barriers along its non-radiative decay, and accounting for the ps signal registered in pump-probe experiments as well as for an efficient population of the triplet state. Cytidine displays an additional ¹n_Nπ* state localised in the N3 atom and that leads to its excited state minimum displaying large potential energy barriers in the pathway connecting to the CI with the ground state. Sugar-to-base hydrogen/proton transfer processes are assessed in solution for the first time, displaying a sizable barrier along its decay and thus being competitive with other slow decay channels in the ps and ns timescales. A unified deactivation scheme for the long-lived channels of pyrimidine nucleosides is delivered, where the ¹n_Oπ* state is found to mediate the long-lived decay in the singlet manifold and act as the doorway for triplet population and thus accounting for the recorded phosphorescence and, more generally, for the transient/photoelectron spectral signals registered up to the ns timescale.

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.

Photoinduced formation mechanism of the thymine–thymine (6–4) adduct in DNA; a QM(CASPT2//CASSCF):MM(AMBER) study

Out-of-plane puckering of the C_4′ atom forming a new CC bond seems essential for oxetane formation.

Probing ultrafast ππ*/nπ* internal conversion in organic chromophores via K-edge resonant absorption

Many photoinduced processes including photosynthesis and human vision happen in organic molecules and involve coupled femtosecond dynamics of nuclei and electrons. Organic molecules with heteroatoms often possess an important excited-state relaxation channel from an optically allowed ππ* to a dark nπ* state. The ππ*/nπ* internal conversion is difficult to investigate, as most spectroscopic methods are not exclusively sensitive to changes in the excited-state electronic structure. Here, we report achieving the required sensitivity by exploiting the element and site specificity of near-edge soft X-ray absorption spectroscopy. As a hole forms in the n orbital during ππ*/nπ* internal conversion, the absorption spectrum at the heteroatom K-edge exhibits an additional resonance. We demonstrate the concept using the nucleobase thymine at the oxygen K-edge, and unambiguously show that ππ*/nπ* internal conversion takes place within (60 ± 30) fs. High-level-coupled cluster calculations confirm the method's impressive electronic structure sensitivity for excited-state investigations.Many photo-induced processes such as photosynthesis occur in organic molecules, but their femtosecond excited-state dynamics are difficult to track. Here, the authors exploit the element and site selectivity of soft X-ray absorption to sensitively follow the ultrafast ππ*/nπ* electronic relaxation of hetero-organic molecules.

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

Theory and Calculation of the Phosphorescence Phenomenon

Phosphorescence is a phenomenon of delayed luminescence that corresponds to the radiative decay of the molecular triplet state. As a general property of molecules, phosphorescence represents a cornerstone problem of chemical physics due to the spin prohibition of the underlying triplet-singlet emission and because its analysis embraces a deep knowledge of electronic molecular structure. Phosphorescence is the simplest physical process which provides an example of spin-forbidden transformation with a characteristic spin selectivity and magnetic field dependence, being the model also for more complicated chemical reactions and for spin catalysis applications. The bridging of the spin prohibition in phosphorescence is commonly analyzed by perturbation theory, which considers the intensity borrowing from spin-allowed electronic transitions. In this review, we highlight the basic theoretical principles and computational aspects for the estimation of various phosphorescence parameters, like intensity, radiative rate constant, lifetime, polarization, zero-field splitting, and spin sublevel population. Qualitative aspects of the phosphorescence phenomenon are discussed in terms of concepts like structure-activity relationships, donor-acceptor interactions, vibronic activity, and the role of spin-orbit coupling under charge-transfer perturbations. We illustrate the theory and principles of computational phosphorescence by highlighting studies of classical examples like molecular nitrogen and oxygen, benzene, naphthalene and their azaderivatives, porphyrins, as well as by reviewing current research on systems like electrophosphorescent transition metal complexes, nucleobases, and amino acids. We furthermore discuss modern studies of phosphorescence that cover topics of applied relevance, like the design of novel photofunctional materials for organic light-emitting diodes (OLEDs), photovoltaic cells, chemical sensors, and bioimaging.

How Does Thymine DNA Survive Ultrafast Dimerization Damage?

The photodimerization reaction between the two adjacent thymine bases within a single strand has been the subject of numerous studies due to its potential to induce DNA mutagenesis and possible tumorigenesis in human skin cells. It is well established that the cycloaddition photoreaction takes place on a picosecond time scale along barrierless or low barrier singlet/triplet pathways. However, the observed dimerization quantum yield in different thymine multimer is considerable lower than might be expected. A reasonable explanation is required to understand why thymine in DNA is able to survive ultrafast dimerization damage. In this work, accurate quantum calculations based on the combined CASPT2//CASSCF/AMBER method were conducted to map the excited state relaxation pathways of the thymine monomer in aqueous solution and of the thymine oligomer in DNA. A monomer-like decay pathway, induced by the twisting of the methyl group, is found to provide a bypass channel to ensure the photostability of thymine in single-stranded oligomers. This fast relaxation path is regulated by the conical intersection between the bright S_CT(¹ππ*) state with the intra-base charge transfer character and the ground state to remove the excess excitation energy, thereby achieving the ground-state recovery with high efficiency.

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

Cyclobutane Thymine Photodimerization Mechanism Revealed by Nonadiabatic Molecular Dynamics

The formation of cyclobutane thymine dimers is one of the most important DNA carcinogenic photolesions induced by ultraviolet irradiation. The long debated question whether thymine dimerization after direct light excitation involves singlet or triplet states is investigated here for the first time using nonadiabatic molecular dynamics simulations. We find that the precursor of this [2 + 2] cycloaddition reaction is the singlet doubly π²π*² excited state, which is spectroscopically rather dark. Excitation to the bright ¹ππ* or dark ¹nπ* excited states does not lead to thymine dimer formation. In all cases, intersystem crossing to the triplet states is not observed during the simulated time, indicating that ultrafast dimerization occurs in the singlet manifold. The dynamics simulations also show that dimerization takes place only when conformational control happens in the doubly excited state.

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.

Singlet-Triplet States Interaction Regions in DNA/RNA Nucleobase Hypersurfaces

The present study provides new insight into the intrinsic mechanisms for the population of the triplet manifold in DNA nucleobases by determining, at the multiconfigurational CASSCF/CASPT2 level, the singlet-triplet states crossing regions and the main decay paths for their lowest singlet and triplet states after near-UV irradiation. The studied singlet-triplet interacting regions are accessible along the minimum energy path of the initially populated singlet bright (1)ππ* state. In particular, all five natural DNA/RNA nucleobases have, at the end of the main minimum energy path and near a conical intersection of the ground and (1)ππ* states, a low-energy, easily accessible, singlet-triplet crossing region directly connecting the lowest singlet and triplet ππ* excited states. Adenine, thymine, and uracil display additional higher-energy crossing regions related to the presence of low-lying singlet and a triplet nπ* state. These funnels are absent in guanine and cytosine, which have the bright (1)ππ* state lower in energy and less accessible nπ* states. Knowledge of the location and accessibility of these regions, in which the singlet-triplet interaction is related to large spin-orbit coupling elements, may help to understand experimental evidence such as the wavelength dependence measured for the triplet formation quantum yield in nucleobases and the prevalence of adenine and thymine components in the phosphorescence spectra of DNA.

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.

Resolving the Benzophenone DNA-Photosensitization Mechanism at QM/MM Level

Benzophenone, the parent of the diarylketone family, is a versatile compound commonly used as a UV blocker. It may also trigger triplet-based DNA photosensitization. Therefore, benzophenone is involved in DNA photodamage induction. In the absence of experimentally resolved structure, the mechanism of DNA damage production remains elusive. Employing a hybrid quantum mechanics/molecular mechanics approach, here we address the spin transfer mechanism between this drug and proximal thymine, that is, the DNA nucleobase most prone to suffer triplet damages.

Non-radiative deactivation in phenol–pyridine complex: theoretical study

A comparative theoretical study of the nonradiative decay mechanisms to the ground state of phenol-pyridine complex under Cs symmetry (Pl and Pe conformers) and without symmetry constraints has been performed by ab initio calculations.

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.

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.

Physical quenching in competition with the formation of cyclobutane pyrimidine dimers in DNA photolesion

The potential energy profiles toward formation of cyclobutane pyrimidine dimers CPD and the physical quenching after UV excitation were explored for the dinucleotide thymine dinucleoside monophosphate (TpT) using density functional theory (ωB97XD) and the time-dependent density functional theory (TD-ωB97XD). The ωB97XD functional that includes empirical dispersion correction is shown to be an appropriate method to obtain rational results for the current large reaction system of TpT. Photophysical quenching is shown to be predominant over the photochemical CPD formation. Following the initial excitation to the (1)ππ* state, the underlying dark (1)nπ* state bifurcates the excited population to the prevailing IC to S0 and the small ISC to the long-lived triplet state T1 via T4 ((3)ππ*) state that has negligible energy gap with (1)nπ* state. Even for the reactive T1 state, two physical quenching pathways resulting in the conversion back to ground-state reactant via the T1/S0 crossing points are newly located, which are in strong competition with CPD formation. These results provide rationale for the recently observed nanosecond triplet decay rates in the single-stranded (dT)18 and inefficiency of deleterious CPD formation, which allow for a deeper understanding of DNA photostability.

Coupled-cluster and density functional theory studies of the electronic 0–0 transitions of the DNA bases

The 0–0 transitions of the electronic excitation spectra of the lowest tautomers of the four nucleotide (DNA) bases have been studied using linear-response approximate coupled-cluster singles and doubles (CC2) calculations.

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.

Photoinduced formation mechanism of the thymine-thymine (6-4) adduct

The photoinduced mechanism leading to the formation of the thymine-thymine (6-4) photolesion has been studied by using the CASPT2//CASSCF approach over a dinucleotide model in vacuo. Following light absorption, localization of the excitation on a single thymine leads to fast singlet-triplet crossing that populates the triplet (3)(nπ*) state of thymine. This state, displaying an elongated C(4)═O bond, triggers (6-4) dimer formation by reaction with the C(5)═C(6) double bond of the adjacent thymine, followed by a second intersystem crossing, which acts as a gate between the excited state of the reactant and the ground state of the photoproduct. The requirement of localized excitation on just one thymine, whose main decay channel (by radiationless repopulation of its ground state) is nonphotochemical, can rationalize the experimentally observed low quantum yield of formation for the thymine-thymine (6-4) adduct.

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.

Multiconfigurational Second-Order Perturbation Theory Restricted Active Space (RASPT2) Method for Electronic Excited States: A Benchmark Study

The recently developed second-order perturbation theory restricted active space (RASPT2) method has been benchmarked versus the well-established complete active space (CASPT2) approach. Vertical excitation energies for valence and Rydberg excited states of different groups of organic (polyenes, acenes, heterocycles, azabenzenes, nucleobases, and free base porphin) and inorganic (nickel atom and copper tetrachloride dianion) molecules have been computed at the RASPT2 and multistate (MS) RASPT2 levels using different reference spaces and compared with CASPT2, CCSD, and experimental data in order to set the accuracy of the approach, which extends the applicability of multiconfigurational perturbation theory to much larger and complex systems than previously. Relevant aspects in multiconfigurational excited state quantum chemistry such as the valence-Rydberg mixing problem in organic molecules or the double d-shell effect for first-row transition metals have also been addressed.