A Plausible Mechanism of Uracil Photohydration Involves an Unusual Intermediate

Formation of Ground-State Intermediate during Electronic Relaxation of Pyrimidine Nucleobases

Ultrafast electronic relaxation of nucleobases from the 1ππ* state to the ground electronic state (S0) is crucial for the photostability of DNA and RNA. However, it has been suggested that electronic relaxation of pyrimidine nucleobases, nucleosides, and nucleotides in an aqueous environment generates an electronically excited intermediate state with a lifetime of tens to hundreds of picoseconds with a relatively high quantum yield (QY) of 0.2-0.5. The generation of such a long-lived excited state seems to be inconsistent with the photostability of these molecules. A recent extreme ultraviolet time-resolved photoelectron spectroscopy study by Miura et al. [J. Am. Chem. Soc. 2023, 145 (6), 3369] reinvestigated this problem and revealed that the accurately determined QY for long-lived excited states is much too low to allow an electronically excited reaction intermediate to exist. In the present study, we investigated the nature of the reaction intermediate using ultraviolet and infrared transient absorption spectroscopy, along with quantum chemical calculations to show that the intermediate is in the S0 state, and its infrared spectrum is compatible with a metastable twisted C═C species theoretically predicted by Park et al. [J. Phys. Chem. Lett. 2022, 13 (30), 7072].

A Computational Study Revealing the Unexplored Mechanisms of the Photo-isomerization of 4H-imidazole N-oxides and the Reverse Thermal Isomerization of their Photo-Products

The current computational study explores the photo-isomerization mechanisms of 4H-imidazole N-oxide systems leading to the formations of their experimentally reported photoproducts. Additionally, the reverse thermal isomerization reactions are also investigated. Quantum mechanical studies on 5-phenyl-4,4-dimethyl-4H-imidazole 3-oxides (IMO) and 1,3-dioxides(IMDO) with substitutions at 2-position reveal their excited state decay pathways forming the oxaziridine and trans,cis-dioxaziridine photoproducts, respectively. The vertically excited S2 state of 2-methyl-substituted IMO undergoes S2/S1 and S0/S1 conical intersections(CIs) through an upward CNO twist to form the respective oxaziridine. On the other hand, in the 2-phenyl-substituted IMO, a pathway through three consecutive CIs (S3/S2,S2/S1,S0/S1) leads to a downward twisted oxaziridine. The methyl-substituted IMDO is found to form both upward and downward twisted oxaziridines through low-lying S0/S1 CIs. The S0-S2 photo-excitation of these oxaziridines gives rise to their respective trans-dioxaziridines. The energy barriers of the oxaziridine→parent nitrone conversion processes are much lower in the studied IMO systems(27-33 kcal/mol) than their 2H-imidazole N-oxide analogues(~38 kcal/mol) which supports the experimentally reported faster thermal isomerizations of the former. These barrier heights in 4H-imidazole 3-oxides decrease from 2-methyl to 2-phenyl-substituted systems, and in the latter, it decreases with the increasing electron donating ability of the groups on phenyl (order of barrier height:-NO2>-OMe>-NMe2).

Expanding Horizons in Quantum Chemical Studies: The Versatile Power of MRSF-TDDFT

ConspectusWhile traditional quantum chemical theories have long been central to research, they encounter limitations when applied to complex situations. Two of the most widely used quantum chemical approaches, Density Functional Theory (DFT) and Time-Dependent Density Functional Theory (TDDFT), perform well in cases with relatively weak electron correlation, such as the ground-state minima of closed-shell systems (Franck-Condon region). However, their applicability diminishes in more demanding scenarios. These limitations arise from the reliance of DFT on a single-determinantal framework and the inability of TDDFT to capture double and higher excited configurations in its response space.The recently developed Multi-Reference Spin-Flip Time-Dependent Density Functional Theory (MRSF-TDDFT) successfully overcomes these challenges, pushing the boundaries of DFT methods. MRSF-TDDFT is exceptionally versatile, making it suitable for various applications, including bond-breaking and bond-forming reactions, open-shell singlet systems such as diradicals, and a more accurate depiction of transition states. It also provides the correct topology for conical intersections (CoIns) and incorporates double excitations into the response space for a more precise description of excited states. With the help of its formal framework, core-hole relaxation for accurate X-ray absorption prediction can be also done readily. Notably, MRSF-TDDFT achieves an equal footing description of ground and excited states, with its dual-reference framework ensuring a balanced treatment of both dynamic and nondynamic electron correlations for high accuracy.In predictive tasks, such as calculating adiabatic singlet-triplet gaps, MRSF-TDDFT achieves accuracy comparable to that of far more computationally expensive coupled-cluster methods. The missing doubly excited state of H2 observed in TDDFT is accurately captured by MRSF-TDDFT, which also reproduces the correct asymptotic bond-breaking potential energy surface. Furthermore, the CoIns of butadiene, missed by both TDDFT and Complete-Active Space Self-Consistent Field (CASSCF) methods, are successfully recovered by MRSF-TDDFT, achieving results consistent with high-level theories, an important aspect for successful study of photochemical processes. Additionally, the common issue of CASSCF overestimating bright states (ionic states) due to the missing dynamic correlation is effectively resolved by MRSF-TDDFT.Despite its numerous advancements, MRSF-TDDFT retains the computational efficiency of conventional TDDFT, making it a practical tool for routine calculations. In addition, it has been demonstrated that the prediction accuracy of MRSF-TDDFT can be further enhanced through the development of tailor-made exchange-correlation functionals, paving the way for the creation of new, specialized functionals. Consequently, with its remarkable versatility, high accuracy, and computational practicality, this innovative method significantly expands scientists' ability to explore complex molecular behaviors and design advanced materials, including applications in photobiology, organic LEDs, photovoltaics, and spintronics, to name a few.

Superior Photostability of the Unnatural Base 6-Amino-5-nitropyridin-2-ol: A Case Study Using Ultrafast Broadband Fluorescence, Transient Absorption, and Theoretical Computation

6-Amino-5-nitropyridin-2-ol (Z), a nitroaromatic compound and a base for Hachimoji nucleic acids, holds significant potential in expanding the genetic alphabet, as well as in synthetic biology and biotechnology. Despite its promising applications, the spectral characterization and photoinduced properties of Z have remained largely unexplored until now. This study presents a comprehensive investigation into its excited state dynamics in various solvents, utilizing state-of-the-art ultrafast broadband time-resolved fluorescence and transient absorption spectroscopy, complemented by computational methods. The acquired results provide direct experimental evidence that, upon photoexcitation, Z emits prompt fluorescence from a nearly planar structure in its excited state, independent of solvent properties. This state deactivates nonradiatively within sub-picoseconds through internal conversion with a unitary yield, primarily mediated by the rotation of the nitro group. This unusually rapid deactivation pathway entirely excludes the involvement of long-lived nπ* states, triplet states, and photoproducts, which are commonly observed in most nitroaromatic compounds and natural DNA and RNA bases. Our findings underscore that Z, as an unnatural base, exhibits superior photostability compared to canonical natural bases. This provides valuable insights into the photodynamics of nitroaromatic compounds, which is beneficial for strategic substitution design in environmental and biological applications.

Temperature Controlled Decay and Pendulum Dynamics of Green Fluorescent Protein (GFP) Chromophore

The excited-state dynamics of the GFP chromophore, HBDI- (anionic p-hydroxybenzylidene-2,3-dimethylimidazolinone), were investigated through a combination of theoretical nonadiabatic molecular dynamics (NAMD) simulations and femtosecond transient absorption spectroscopy (fs-TA). The NAMD simulations revealed that the primary dynamics in excited states involve the formation of a P-twisted intermediate (S1min,P), which undergoes pendulum-like oscillations with respect to ϕ = 90°. This motion serves as a reservoir for the excited-state population and the primary source of fluorescence. Rather than a direct channel from the major S1min,P, a coordinated pathway of S1min,P → S1min → S1min,I → S0 is responsible for the decay to the ground state, emphasizing the importance of planar intermediate (S1min) formation. The experimental fs-TA spectra confirmed these dynamics, revealing three distinct time scales (340-470 fs, 1.4 ps, and 8.3 ps), corresponding to the formation of S1min,P and its decay governed by the coordinated pathway. At low temperatures, the coordinated decay pathway is suppressed, leading to prolonged fluorescence lifetimes, consistent with low-temperature experimental results. This study presents a new model for the excited-state dynamics of GFP chromophore, suggesting that pendulum motion and the coordinated decay pathway play a crucial role in regulating fluorescence intensity.

OpenQP: A Quantum Chemical Platform Featuring MRSF-TDDFT with an Emphasis on Open-Source Ecosystem

The OpenQP (Open Quantum Platform) is a new open-source quantum chemistry library developed to tackle sustainability and interoperability challenges in the field of computational chemistry. OpenQP provides various popular quantum chemical theories as autonomous modules such as energy and gradient calculations of HF, DFT, TDDFT, SF-TDDFT, and MRSF-TDDFT, thereby allowing easy interconnection with third-party software. A scientifically notable feature is the innovative mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) and its customized exchange-correlation functionals such as the DTCAM series of VAEE, XI, XIV, AEE, and VEE, which significantly expand the applicability scope of DFT and TDDFT. OpenQP also supports parallel execution and is optimized with BLAS and LAPACK for high performance. Future enhancements such as extended Koopman's theorem (EKT)-MRSF-TDDFT and spin-orbit coupling (SOC)-MRSF-TDDFT will further expand OpenQP's capabilities. Additionally, a Python wrapper PyOQP is provided that performs tasks such as geometry optimization, conical intersection searches, and nonadiabatic coupling calculations, among others, by prototyping the modules of the OpenQP library in combination with third-party libraries. Overall, OpenQP aligns with modern trends in high-performance scientific software development by offering flexible prototyping and operation while retaining the performance benefits of compiled languages like Fortran and C. They enhance the sustainability and interoperability of quantum chemical software, making OpenQP a crucial platform for accelerating the development of advanced quantum theories like MRSF-TDDFT.

UMRSF-TDDFT: Unrestricted Mixed-Reference Spin-Flip-TDDFT

An unrestricted version of Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory (UMRSF-TDDFT) was developed based on unrestricted Kohn-Sham orbitals (UKS) with a new molecular orbital (MO) reordering scheme. Additionally, a simple yet accurate method for estimating ⟨S2⟩ expectation values was devised. UMRSF-TDDFT was benchmarked against cases where DFT, TDDFT, and SF-TDDFT traditionally fail to provide accurate descriptions. In an application to the ground and excited states of a Be atom, UMRSF-TDDFT successfully recovers the degenerate states, with its energies slightly reduced compared to its RO counterpart, due to the additional variational flexibility of UKS. A clear difference between UMRSF and U-SF-TDDFT is evident in the bond breaking of the hydrogen fluoride system, as the latter misses an important configuration. In the case of the Jahn-Teller distortion of trimethylenemethane (TMM), the relative singlet energy compared to the triplet is lower by 0.1 and 0.2 eV for UMRSF and U-SF-TDDFT, respectively, than that of MRSF-TDDFT. The reduction in UMRSF energy is attributed to spatial orbital relaxations, whereas the reduction in U-SF-TDDFT energy results from spin contamination. Overall, the additional orbital relaxations afforded by unrestricted Kohn-Sham (UKS) orbitals in UMRSF-TDDFT lead to lower total system energies compared to their restricted open-shell counterparts. This enhancement adds a practical and accurate quantum chemical theory to the existing RO variant for addressing challenging systems where traditional quantum theories suffer.

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.

The photochemical trans → cis and thermal cis → trans isomerization pathways of azobenzo-13-crown ether: A computational study on a strained cyclic azobenzene system

The isomerization of azobenzo-13-crown ether can be expected to be hindered due to the polyoxyethylene linkage connecting the 2,2'-positions of azobenzene. The mixed reference spin-flip time-dependent density functional theory results reveal that the planar and rotational minima of the first photo-excited singlet state (S1) of the trans-isomer pass through a barrier (2.5-5.0 kcal/mol) as it goes toward the torsional conical intersection (S0/S1) geometry ( Collapse

Key Words Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Dilawar Singh Sisodiya Department of Chemistry, Birla Institute of Technology and Science (BITS), Pilani, K. K. Birla Goa Campus, Zuarinagar, India
Anjan Chattopadhyay Department of Chemistry, Birla Institute of Technology and Science (BITS), Pilani, K. K. Birla Goa Campus, Zuarinagar, India

Collapse

Comparative Study of Uracil Excited-State Photophysics in Water and Acetonitrile via RMS-CASPT2-Driven Quantum-Classical Trajectories

We present a nonadiabatic molecular dynamics study of the ultrafast processes occurring in uracil upon UV light absorption, leading to electronic excitation and subsequent nonradiative decay. Previous studies have indicated that the mechanistic details of this process are drastically different depending on whether the process takes place in the gas phase, acetonitrile, or water. However, such results have been produced using quantum chemical methods that did not incorporate both static and dynamic electron correlation. In order to assess the previously proposed mechanisms, we simulate the photodynamics of uracil in the three environments mentioned above using quantum-classical trajectories and, for solvated uracil, hybrid quantum mechanics/molecular mechanics (QM/MM) models driven by the rotated multistate complete active space second-order perturbation (RMS-CASPT2) method. To do so, we exploit the gradient recently made available in OpenMolcas and compare the results to those obtained using the complete active space self-consistent field (CASSCF) method only accounting for static electron correlation. We show that RMS-CASPT2 produces, in general, a mechanistic picture different from the one obtained at the CASSCF level but confirms the hypothesis advanced on the basis of previous ROKS and TDDFT studies thus highlighting the importance of incorporating dynamic electron correlation in the investigation of ultrafast electronic deactivation 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.

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.

High-performance strategies for the recent MRSF-TDDFT in GAMESS

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. This paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.

Ultrafast Electronic Relaxation in 6-Methyluracil and 5-Fluorouracil in Isolated and Aqueous Conditions: Substituent and Solvent Effects

We report ultrafast extreme ultraviolet photoelectron spectroscopy of 6-methyluracil (6mUra) and 5-fluorouracil (5FUra) in the gas phase and 6mUra and 5-fluorouridine in an aqueous environment. In the gas phase, internal conversion (IC) occurs from ¹ππ* to ¹nπ* states in tens of femtoseconds, followed by intersystem crossing to the ³ππ* state in several picoseconds. In an aqueous solution, 6mUra undergoes IC almost exclusively to the ground state (S₀) in about 100 fs, which is essentially the same process as that for unsubstituted uracil, but much faster than that for thymine (5-methyluracil). The different dynamics for C5 and C6 methylation suggest that IC from ¹ππ* to S₀ is facilitated by out-of-plane (OOP) motion of the C5 substituent. The slow IC for C5-substituted molecules in an aqueous environment is ascribed to the solvent reorganization that is required for this OOP motion to occur. The slow rate for 5FUrd may arise in part from an increased barrier height due to C5 fluorination.

Formation of Long-Lived Dark States during Electronic Relaxation of Pyrimidine Nucleobases Studied Using Extreme Ultraviolet Time-Resolved Photoelectron Spectroscopy

Ultrafast electronic relaxation of nucleobases from ¹ππ* states to the ground state (S₀) is considered essential for the photostability of DNA. However, transient absorption spectroscopy (TAS) has indicated that some nucleobases in aqueous solutions create long-lived ¹nπ*/³ππ* dark states from the ¹ππ* states with a high quantum yield of 0.4-0.5. We investigated electronic relaxation in pyrimidine nucleobases in both aqueous solutions and the gas phase using extreme ultraviolet (EUV) time-resolved photoelectron spectroscopy. Femtosecond EUV probe pulses cause ionization from all electronic states involved in the relaxation process, providing a clear overview of the electronic dynamics. The ¹nπ* quantum yields for aqueous cytidine and uracil (Ura) derivatives were found to be considerably lower (<0.07) than previous estimates reported by TAS. On the other hand, aqueous thymine (Thy) and thymidine exhibited a longer ¹ππ* lifetime and a higher quantum yield (0.12-0.22) for the ¹nπ* state. A similar trend was found for isolated Thy and Ura in the gas phase: the ¹ππ* lifetimes are 39 and 17 fs and the quantum yield for ¹nπ* are 1.0 and 0.45 for Thy and Ura, respectively. The result indicates that single methylation to the C₅ position hinders the out-of-plane deformation that drives the system to the conical intersection region between ¹ππ* and S₀, providing a large impact on the photophysics/photochemistry of a pyrimidine nucleobase. The significant reduction of ¹nπ* yield in aqueous solution is ascribed to the destabilization of the ¹nπ* state induced by hydrogen bonding.

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.

Prototropically Controlled Dynamics of Cytosine Photodecay

The effect of the existence of several prototropic tautomers of cytosine on its UV/vis spectra and the excited state decay dynamics is studied by spectral and nonadiabatic molecular dynamics (NAMD) simulations in connection with the mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. Simulated UV/vis spectra provide a strong indication that the H3N keto-amino cytosine tautomer (the least anticipated species) may be present under experimental conditions. The NAMD simulations yield a wide range of excited state decay constants for various tautomers of cytosine, ranging from ∼1.3 ps for the biologically relevant H1N keto-amino tautomer to ∼0.1 ps for the keto-imino tautomer. The slowness of the H1N decay dynamics follows from the presence of a barrier on the excited state energy surface separating the Franck-Condon structure from the major decay funnel, the conical intersection seam. It is suggested that the experimentally observed photodecay dynamics may result from a combination of the decay processes of various tautomers (H3N in particular) present simultaneously under the experimental conditions.

Ultrafast Excited State Aromatization in Dihydroazulene

Excited-state aromatization dynamics in the photochemical ring opening of dihydroazulene (DHA) is investigated by nonadiabatic molecular dynamics simulations in connection with the mixed-reference spin-flip (MRSF)-TDDFT method. It is found that, in the main reaction channel, the ring opening occurs in the excited state in a sequence of steps with increasing aromaticity. The first stage lasting ca. 200 fs produces an 8π semiaromatic S₁ minimum (S_1, min) through an ultrafast damped bond length alternation (BLA) movement synchronized with a partial planarization of the cycloheptatriene ring. An additional ca. 200 fs are required to gain the vibrational energy needed to overcome a ring-opening transition state characterized by an enhanced Baird aromaticity. Unlike other BLA motions of ππ* state, it was shown that their damping is a characteristic feature of aromatic bond-equalization process. In addition, some minor channels of the reaction have also been discovered, where noticeably higher barriers of the S₁ non/antiaromatic transition structures must be surmounted. These anti-Baird channels led to reformation of DHA or other closed-ring products. The observed competition between the Baird and anti-Baird channels suggests that the quantum yield of photochemical products can be controllable by tipping their balance. Hence, here we suggest including the concept of anti-Baird, which would expand the applicability of Baird rule to much broader situations.

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.