Excited State Decay Mechanisms Research Articles

From indigo to isoindigo: Rationalizing the high efficiency of photoprotective molecular mechanisms

Indigo (IND) and its natural derivatives, including Tyrian purple and indirubin (INR), are examples of molecules with a history of millions of years. They also illustrate molecular evolution and longevity, linked to their extraordinary molecular stability associated with excited state decay mechanisms. In the excited state, these molecules efficiently undergo a radiationless deactivation process resulting in negligible fluorescence and triplet state formation (less than 0.1 % of the quantum loss in indigo for both processes). Here, we demonstrate that in isoindigo (ISO) this process is even more efficient than in IND, resulting in nearly 100 % of the excited state deactivation. In IND, the molecular mechanism behind this decay process involves intramolecular proton transfer (ESPT). In the case of INR, ESPT efficiently competes with the formation of a conformational species –the syn-rotamer– for deactivation to the ground state. With ISO, the ESPT process is absent. Instead, the deactivation mechanism involves a rotation of a few degrees around the central C–C bond, aRotational Isomerization (RI). This rotation leads to a sloped conical intersection (CI), making the radiationless deactivation process more efficient than the ESPT found with IND.

Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation.

Herein, we have used the "on-the-fly" ring-polymer surface-hopping simulation method with the centroid approximation (RPSH-CA), in combination with the multireference OM2/MRCI electronic structure calculations to study the photoinduced dynamics of a green fluorescent protein (GFP) chromophore analogue in the gas phase, i.e., o-HBI, at 50, 100, and 300 K with 1, 5, 10, and 15 beads (3600 1 ps trajectories). The electronic structure calculations identified five new minimum-energy conical intersection (MECI) structures, which, together with the previous one, play crucial roles in the excited-state decay dynamics of o-HBI. It is also found that the excited-state intramolecular proton transfer (ESIPT) occurs in an ultrafast manner and is completed within 20 fs in all the simulation conditions because there is no barrier associated with this ESIPT process in the S1 state. However, the other excited-state dynamical results are strongly related to the number of beads. At 50 and 100 K, the nuclear quantum effects (NQEs) are very important; therefore, the excited-state dynamical results change significantly with the bead number. For example, the S1 decay time deduced from time-dependent state populations becomes longer as the bead number increases. Nevertheless, an essentially convergent trend is observed when the bead number is close to 10. In contrast, at 300 K, the NQEs become weaker and the above dynamical results converge very quickly even with 1 bead. Most importantly, the NQEs seriously affect the excited-state decay mechanism of o-HBI. At 50 and 100 K, most trajectories decay to the S0 state via perpendicular keto MECIs, whereas, at 300 K, only twisted keto MECIs are responsible for the excited-state decay. The present work not only comprehensively explores the temperature-dependent photoinduced dynamics of o-HBI, but also demonstrates the importance and necessity of NQEs in nonadiabatic dynamics simulations, especially at relatively low temperatures.

Enhanced Visible Light Absorption in Heteroleptic Cuprous Phenanthrolines.

This work presents a series of Cu(I) heteroleptic 1,10-phenanthroline chromophores featuring enhanced UVA and visible-light-harvesting properties manifested through vectorial control of the copper-to-phenanthroline charge-transfer transitions. The molecules were prepared using the HETPHEN strategy, wherein a sterically congested 2,9-dimesityl-1,10-phenanthrolne (mesPhen) ligand was paired with a second phenanthroline ligand incorporating extended π-systems in their 4,7-positions. The combination of electrochemistry, static and time-resolved electronic spectroscopy, 77 K photoluminescence spectra, and time-dependent density functional theory calculations corroborated all of the experimental findings. The model chromophore, [Cu(mesPhen)(phen)]+ (1), lacking 4,7-substitutions preferentially reduces the mesPhen ligand in the lowest energy metal-to-ligand charge-transfer (MLCT) excited state. The remaining cuprous phenanthrolines (2-4) preferentially reduce their π-conjugated ligands in the low-lying MLCT excited state. The absorption cross sections of 2-4 were enhanced (εMLCTmax = 7430-9980 M-1 cm-1) and significantly broadened across the UVA and visible regions of the spectrum compared to 1 (εMLCTmax = 6494 M-1 cm-1). The excited-state decay mechanism mirrored those of long-lived homoleptic Cu(I) phenanthrolines, yielding three distinguishable time constants in ultrafast transient absorption experiments. These represent pseudo-Jahn-Teller distortion (τ1), singlet-triplet intersystem crossing (τ2), and the relaxed MLCT excited-state lifetime (τ3). Effective light-harvesting from Cu(I)-based chromophores can now be rationalized within the HETPHEN strategy while achieving directionality in their respective MLCT transitions, valuable for integration into more complex donor-acceptor architectures and longer-lived photosensitizers.

Quantum mechanics/molecular mechanics studies on the excited-state decay mechanisms of cytidine aza-analogues: 5-azacytidine and 2'-deoxy-5-azacytidine in aqueous solution.

The excited state properties and deactivation pathways of two DNA methylation inhibitors, i.e., 5-azacytidine (5ACyd) and 2'-deoxy-5-azacytidine (5AdCyd) in aqueous solution, are comprehensively explored with the QM(CASPT2//CASSCF)/MM protocol. We systematically map the feasible decay mechanisms based on the obtained excited-state decay paths involving all the identified minimum-energy structures, conical intersections, and crossing points driving the different internal conversion (IC) and intersystem crossing (ISC) routes in and between the 1ππ*, 1nπ*, 3ππ*, 3nπ*, and S0 states. Unlike the 1nπ* state below the 1ππ* state in 5ACyd, deoxyribose group substitution at the N1 position leads to the 1ππ* state becoming the S1 state in 5AdCyd. In 5ACyd and 5AdCyd, the initially populated 1ππ* state mainly deactivates to the S0 state through the direct 1ππ* → S0 IC or mediated by the 1nπ* state. The former nearly barrierless IC channel of 1ππ* → S0 occurs ultrafast via the nearby low-lying 1ππ*/S0 conical intersection. In the latter IC channel of 1ππ* → 1nπ* → S0, the initially photoexcited 1ππ* state first approaches the nearby S2/S1 conical section 1ππ*/1nπ* and then undergoes efficient IC to the 1nπ* state, followed by the further IC to the initial S0 state via the S1/S0 conical intersection 1nπ*/S0. The 1nπ*/S0 conical intersection is estimated to be located 6.0 and 4.9 kcal mol-1 above the 1nπ* state minimum in 5ACyd and 5AdCyd, respectively, at the QM(CASPT2)/MM level. In addition to the efficient singlet-mediated IC channels, the minor ISC routes would populate 1ππ* to T1(ππ*) through 1ππ* → T1 or 1ππ* → 1nπ* → T1. Relatively, the 1ππ* → 1nπ* → T1 route benefits from the spin-orbit coupling (SOC) of 1nπ*/3ππ* of 8.7 cm-1 in 5ACyd and 10.2 cm-1 in 5AdCyd, respectively. Subsequently, the T1 system will approach the nearby T1/S0 crossing point 3ππ*/S0 driving it back to the S0 state. Given the 3ππ*/S0 crossing point located above the T1 minimum and the small T1/S0 SOC, i.e., 8.4 kcal mol-1 and 2.1 cm-1 in 5ACyd and 6.8 kcal mol-1 and 1.9 cm-1 in 5AdCyd, respectively, the slow T1 → S0 would trap the system in the T1 state for a while. The present work could contribute to understanding the mechanistic photophysics and photochemistry of similar aza-nucleosides and their derivatives.

Intense Photoinduced Intervalence Charge Transfer in High-Valent Iron Mixed Phenolate/Carbene Complexes.

We report high‐valent iron complexes supported by N‐heterocyclic carbene (NHC)‐anchored, bis‐phenolate pincer ligands that undergo ligand‐to‐metal charge transfer (LMCT) upon photoexcitation. The resulting excited states – with a lifetime in the picosecond range – feature a ligand‐based, mixed‐valence system and intense intervalence charge transfer bands in the near‐infrared region. Upon oxidation of the complex, corresponding intervalence charge transfer absorptions are also observed in the ground state. We suggest that the spectroscopic hallmarks of such LMCT states provide useful tools to decipher excited‐state decay mechanisms in high‐valent NHC complexes. Our observations further indicate that NHC‐anchored, bis‐phenolate pincer ligands are not sufficiently strong donors to prevent the population of excited metal‐centered states in high‐valent iron complexes.

Ultrafast Excited-State Decay Mechanisms of 6-Thioguanine Followed by Sub-20 fs UV Transient Absorption Spectroscopy.

Understanding the primary steps following UV photoexcitation in sulphur-substituted DNA bases (thiobases) is fundamental for developing new phototherapeutic drugs. However, the investigation of the excited-state dynamics in sub-100 fs time scales has been elusive until now due to technical challenges. Here, we track the ultrafast decay mechanisms that lead to the electron trapping in the triplet manifold for 6-thioguanine in an aqueous solution, using broadband transient absorption spectroscopy with a sub-20 fs temporal resolution. We obtain experimental evidence of the fast internal conversion from the S2(ππ*) to the S1(nπ*) states, which takes place in about 80 fs and demonstrates that the S1(nπ*) state acts as a doorway to the triplet population in 522 fs. Our results are supported by MS-CASPT2 calculations, predicting a planar S2(ππ*) pseudo-minimum in agreement with the stimulated emission signal observed in the experiment.

Theoretical Studies on the Excited-State Decay Mechanism of Homomenthyl Salicylate in a Gas Phase and an Acetonitrile Solution

Here, we employ the CASPT2//CASSCF and QM(CASPT2//CASSCF)/MM approaches to explore the photochemical mechanism of homomenthyl salicylate (HMS) in vacuum and an acetonitrile solution. The results show that in both cases, the excited-state relaxation mainly involves a spectroscopically "bright" S1(1ππ*) state and the lower-lying T1 and T2 states. In the major relaxation pathway, the photoexcited S1 keto system first undergoes an essentially barrierless excited-state intramolecular proton transfer (ESIPT) to generate the S1 enol minimum, near which a favorable S1/S0 conical intersection decays the system to the S0 state followed by a reverse ground-state intramolecular proton transfer (GSIPT) to repopulate the initial S0 keto species. In the minor one, an S1/T2/T1 three-state intersection in the keto region makes the T1 state populated via direct and T2-mediated intersystem crossing (ISC) processes. In the T1 state, an ESIPT occurs, which is followed by ISC near a T1/S0 crossing point in the enol region to the S0 state and finally back to the S0 keto species. In addition, a T1/S0 crossing point near the T1 keto minimum can also help the system decay to the S0 keto species. However, small spin-orbit couplings between T1 and S0 at these T1/S0 crossing points make ISC to the S0 state very slow and make the system trapped in the T1 state for a while. The present work rationalizes not only the ultrafast excited-state decay dynamics of HMS but also its low quantum yield of phosphorescence at 77 K.

Tracking excited state decay mechanisms of pyrimidine nucleosides in real time

DNA owes its remarkable photostability to its building blocks—the nucleosides—that efficiently dissipate the energy acquired upon ultraviolet light absorption. The mechanism occurring on a sub-picosecond time scale has been a matter of intense debate. Here we combine sub-30-fs transient absorption spectroscopy experiments with broad spectral coverage and state-of-the-art mixed quantum-classical dynamics with spectral signal simulations to resolve the early steps of the deactivation mechanisms of uridine (Urd) and 5-methyluridine (5mUrd) in aqueous solution. We track the wave packet motion from the Franck-Condon region to the conical intersections (CIs) with the ground state and observe spectral signatures of excited-state vibrational modes. 5mUrd exhibits an order of magnitude longer lifetime with respect to Urd due to the solvent reorganization needed to facilitate bulky methyl group motions leading to the CI. This activates potentially lesion-inducing dynamics such as ring opening. Involvement of the 1nπ* state is found to be negligible.

Theoretical Determination of Rate Constants from Excited States: Application to Benzophenone.

A cost-effective method of theoretically predicting electronic-transition rate constants from the excited states of molecules is reported. This method is based on density functional theory calculations of electronic states and quantitative rate constant determination with the Fermi golden rule. The method is applied to the theoretical determination of the excited-state decay mechanism of photoexcited benzophenone, a representative molecule in photochemistry and biochemistry. Calculated rate constants for benzophenone are quantitatively consistent with experimental ones, which validates the reliability of our rate constant calculation. The calculated population kinetics indicate that S1 → T2 → T1 → S0 is the predominant decay pathway.

Photoprotective Properties of Eumelanin: Computational Insights into the Photophysics of a Catechol:Quinone Heterodimer Model System

Melanins are skin-centered molecular structures that block harmful UV radiation from the sun and help protect chromosomal DNA from UV damage. Understanding the photodynamics of the chromophores that make up eumelanin is therefore paramount. This manuscript presents a multi-reference computational study of the mechanisms responsible for the experimentally observed photostability of a melanin-relevant model heterodimer comprising a catechol (C)–benzoquinone (Q) pair. The present results validate a recently proposed photoinduced intermolecular transfer of an H atom from an OH moiety of C to a carbonyl-oxygen atom of the Q. Photoexcitation of the ground state C:Q heterodimer (which has a π-stacked “sandwich” structure) results in population of a locally excited ππ* state (on Q), which develops increasing charge-transfer (biradical) character as it evolves to a “hinged” minimum energy geometry and drives proton transfer (i.e., net H atom transfer) from C to Q. The study provides further insights into excited state decay mechanisms that could contribute to the photostability afforded by the bulk polymeric structure of eumelanin.

A theoretical study on the excited-state deactivation paths for the A-5FU dimer.

The photostability of DNA plays a key role in the normal function of organisms. A-5FU is a base pair derivative of the A-T dimer where the methyl group is replaced by a F atom. Here, accurate static TDDFT calculations and non-adiabatic dynamic simulations are used to systematically investigate the excited-state decay paths of the A-5FU dimer related to the proton transfer and the out-of-plane twisting deformation motion of A and 5FU in the 1ππ* and 1nπ* states. CC2 is used to check the accuracy of the current TDDFT calculations. Our results show that the deformation of the C[double bond, length as m-dash]C or C[double bond, length as m-dash]N double bond in A and 5FU provides an efficient pathway for the depopulation of the lowest excited states, which can compete with the excited-state proton transfer paths in the dimer. This finding indicates that monomer-like decay paths could be important for the photostability of weakly hydrogen-bonded DNA base pairs and provide a new insight into the excited-state decay paths in base pairs and their analogues.

Ultrafast excited state dynamics of biliverdin dimethyl ester coordinate with zinc ions

As one of the biological endogenous pigments, biliverdin (BV) and its dimethyl ester (BVE) have extremely weak fluorescence in solution with quantum yield less than 0.01%. However, the situation reverses with the addition of zinc ions. The strength for fluorescence of BVE-Zn2+ complex is greatly enhanced and fluorescence quantum yield can increase to 5%. Herein, we studied ultrafast excited state dynamics of BVE-Zn2+ complex in ethanol, n-propanol, and DMSO solutions in order to reveal the mechanism of fluorescence quantum yield enhancement. The results show that BVE can form a stable coordination complex with zinc with 1:1 stoichiometry in solution. BVE is structurally and energetically more stable in the complex. Using picosecond time-resolve fluorescence and femtosecond transient absorption spectroscopy, we show that smaller non-radiative rate constant of BVE-Zn2+ complex in DMSO is the key to increasing its fluorescence quantum yield and the excited state decay mechanism is also revealed. These results provide valuable information about the fluorescence property change after BVE binding to metal ions and may provide a guidance for the study of phytochromes or other fluorescence proteins in which BV/BVE acts as chromophores.

Potential‐Energy Surface and Dynamics Simulation of THBDBA: An Annulated Tetraphenylethene Derivative Combining Aggregation‐Induced Emission and Switch Behavior

AbstractElectronic structure calculations and nonadiabatic dynamics simulations were carried out on the excited state decay mechanism of a typical aggregation‐induced emission (AIE)‐active molecule, 10,10′,11,11′‐tetrahydro‐5,5′‐bidibenzo[a,d][7]annulenylidene (THBDBA), which was experimentally found to be almost non‐emissive in THF solution but highly luminescent as an aggregate. At the MS‐CASPT2//CASSCF level, the decay path from S1 leads to a conical intersection (CI) seam with the ground state (S0) via a cyclization process of two neighboring phenyl groups. Nonadiabatic dynamics simulations of 367 surface hopping trajectories show that all the structures of hopping points resemble those of the S1/S0 CI seam with cyclized geometry. The final structures of trajectories hopping to S0 belong to either Franck‐Condon (FC) like structures (66.8 %) or cyclized conformers (33.2 %), which agrees well with the experimental result of 23 % cyclization product (Y. Cai et al. Chem. Sci., 2011, 2, 2029–2034). Our results also indicate that long‐time irradiation leads to a photostationary state where the open and cyclized forms interconvert photochemically. The preference for photocyclization over double bond isomerization is due to the rigidity of the dihydrofulvene rings. From the point of view of AIE, the calculations suggest that THBDBA follows our restricted access to a CI (RACI) model because the photocyclization requires large displacements of the phenyl rings. These displacements must be blocked in the aggregate phase, inhibiting the decay and leading to luminescence. Our study contributes to understand the AIE and switch behavior of THBDBA and will be useful to support the design of new molecules with these features.

New insights into the dissociation dynamics of methylated anilines.

Aniline, an important model system for biological chromophores, undergoes ultrafast H-atom loss upon absorption of an ultraviolet photon. By varying the number and position of methyl substituents on both the aromatic ring and amine functional group, we explore the ultrafast production of photofragments as a function of molecular structure. Both N-methyl aniline and 3,5-dimethyl aniline show altered H-atom loss behaviour compared to aniline, while no evidence for CH3 loss was found from either N-methyl aniline or N,N-dimethyl aniline. With the addition of time-resolved photoelectron spectroscopy, the photofragment appearance times are matched to excited state relaxation pathways. Evidence for a sequential excited state relaxation mechanism, potentially involving a valence-to-Rydberg decay mechanism, will be presented. Such a global, bottom-up approach to molecular photochemistry is crucial to understanding the dissociative pathways and excited state decay mechanisms of biomolecule photoprotection in nature.

UV-Vis, Fluorescence, and Resonance Raman Spectroscopic and Density Functional Theoretical Studies on 3-Amino-1,2,4-triazole: Microsolvation and Solvent-Dependent Nonadiabatic Excited State Decay in Solution.

The microsolvation and photophysics of 3-amino-1,2,4-triazole (3AT) after excitation to the light-absorbing S2(nπ*) state were studied by using resonance Raman spectroscopy and single component artificial force-induced reaction (SC-AFIR) in a global reaction route mapping (GRRM) strategy. The vibrational spectra were assigned on the basis of experimental data and density functional theory (DFT) calculations. The resonance Raman spectra of 3AT were measured to probe the excited state structural dynamics in the Franck-Condon region. The conformations of 3AT(CH3CN)1, 3AT(CH3OH)2, and 3AT(H2O)2 clusters were determined by combining vibrational spectrum experiments and B3LYP/6-311++G(d,p) computations. DFT calculations were carried out to obtain the minimal excitation energies of the lower-lying singlet excited states, and the curve-crossing points. It was revealed that the short-time structural dynamics of 3AT were dominated by the N-N stretching coordinates. An excited state decay mechanism is proposed: 3AT is initially excited to the S2(nπ*) state, then the conical intersection (CI) of the S2(nπ*)/S1(ππ*) potential energy surfaces is crossed, and 3AT then decays to the lower solvent-dependent excited state S1(ππ*). It subsequently returns to the S0 state, accompanied by a large Stokes fluorescence shift, which was interpreted as the stabilized S1(ππ*) excited state bonding to several water molecules via intermolecular hydrogen bonding.

Ultrafast Excited-State Dynamics of Cytosine Aza-Derivative and Analogues.

Excited state dynamics of 5-azacytosine (5-AC), 2,4-diamino-1,3,5-triazine (2,4-DT), and 2-amino-1,3,5-triazine (2-AT) were comprehensively investigated by steady state absorption, fluorescence, and femtosecond transient absorption measurements. Time-dependent density functional theory (TDDFT) calculations were performed to help assign the absorption bands and understand the excited state decay mechanisms. The experimental results of excited singlet state dynamics for 5-AC, 2,4-DT, and 2-AT with femtosecond time resolution were reported for the first time. Two distinct decay pathways, with ∼1 ps and tens of picosecond lifetimes, were observed in 5-AC. Only one decay pathway with 17 ps lifetime was observed in 2,4-DT while an emissive state was found in 2-AT. TDDFT calculations suggest that 5-AC has a dark nπ* (S1) state below the first allowed ππ* (S2) state, which leads to the ultrafast decay of the ππ* state. In 2,4-DT, there is no dark nπ* state below the bright ππ* (S1) state and the 17 ps lifetime is assigned to the relaxation from the ππ* (S1) state to ground state. Two dark nπ* states (S1 and S2) were found in 2-AT, which exhibits much more complex excited state dynamics compared with the other two. Photoluminescence in 2-AT has been confirmed to be fluorescence emission from its bright ππ* (S3) state. Our results strongly suggest that electronic structures are very sensitive to the substitution on the triazine ring and that the photophysical properties of nucleic acid analogues depend highly on their molecular structures.

Excited-State Decay Paths in Tetraphenylethene Derivatives.

The photophysical properties of tetraphenylethene (TPE) compounds may differ widely depending on the substitution pattern, for example, with regard to the fluorescence quantum yield ϕf and the propensity to exhibit aggregation-induced emission (AIE). We report combined electronic structure calculations and nonadiabatic dynamics simulations to study the excited-state decay mechanisms of two TPE derivatives with four methyl substituents, either in the meta position (TPE-4mM, ϕf = 0.1%) or in the ortho position (TPE-4oM, ϕf = 64.3%). In both cases, two excited-state decay pathways may be relevant, namely, photoisomerization around the central ethylenic double bond and photocyclization involving two adjacent phenyl rings. In TPE-4mM, the barrierless S1 cyclization is favored; it is responsible for the ultralow fluorescence quantum yield observed experimentally. In TPE-4oM, both the S1 photocyclization and photoisomerization paths are blocked by non-negligible barriers, and fluorescence is thus feasible. Nonadiabatic dynamics simulations with more than 1000 surface hopping trajectories show ultrafast cyclization upon photoexcitation of TPE-4mM, whereas TPE-4oM remains unreactive during the 1 ps simulations. We discuss the chances for spectroscopic detection of the postulated cyclic photoproduct of TPE-4mM and the relevance of our findings for the AIE process.

The Position of the N Atom Plays a Significant Role for Excited-State Decay of Heterocycles.

We have employed combined electronic structure calculations and nonadiabatic dynamics simulations to study the S1 radiationless deactivation mechanism of pyrazole. In terms of MS-CASPT2 computed results, we propose that the 1πσ* state-driven nonadiabatic N-N dissociation is a major relaxation path; the ring-puckering deformation path as well as the 1πσ* state-driven N-H dissociation are less favorable. This excited-state decay mechanism is supported by MS-CASPT2 nonadiabatic dynamics simulations. The present study demonstrates that pyrazole has a different excited-state radiationless deactivation mechanism compared with its structural isomer imidazole, in which the 1πσ* state-driven nonadiabatic N-H dissociation plays a more important role. However, such a channel is suppressed in pyrazole; instead, the 1πσ* state-driven nonadiabatic N-N dissociation is dominant.

Sequential Proton-Coupled Electron Transfer Mediates Excited-State Deactivation of a Eumelanin Building Block.

Skin photoprotection is commonly believed to rely on the photochemistry of 5,6-dihydroxyindole (DHI)- and 5,6-dihydroxyindole-2-carboxylic acid (DHICA)-based eumelanin building blocks. Attempts to elucidate the underlying excited-state relaxation mechanisms have been partly unsuccessful due to the marked instability to oxidation. We report a study of the excited-state deactivation of DHI using steady-state and time-resolved fluorescence accompanied by high-level quantum-chemistry calculations including solvent effects. Spectroscopic data show that deactivation of the lowest excited state of DHI in aqueous buffer proceeds on the 100 ps time scale and is 20 times faster than in methanol. Quantum-chemical calculations reveal that the excited-state decay mechanism is a sequential proton-coupled electron transfer, which involves the initial formation of a solvated electron from DHI, followed by the transfer of a proton to the solvent. This unexpected finding would prompt a revision of current notions about eumelanin photophysics and photobiology.

The excited-state decay mechanism of 2,4-dithiothymine in the gas phase, microsolvated surroundings, and aqueous solution.

The photophysics of thiothymines has been extensively studied computationally in the past few years due to their significant potential as photosensitizers in photodynamic therapy. However, the corresponding computational studies of the photophysical mechanism of 2,4-dithiothymine are scarce. Herein we have employed the CASPT2//CASSCF and QM(CASPT2//CASSCF)/MM methods to systematically explore the excited-state decay mechanism of 2,4-dithiothymine in isolated, microsolvated, and aqueous surroundings. First, we have optimized minima and conical intersections in and between the lowest six excited singlet and triplet states i.e., , , , , and ; then, based on computed excited-state decay paths and spin-orbit couplings, we have proposed several nonadiabatic pathways that efficiently populate the lowest triplet state to explain the experimentally observed ultrahigh triplet-state quantum yield. Moreover, we have found that the excited-state decay mechanism in microsolvated and aqueous environments is more complicated than that in the gas phase. The solute-solvent interaction has significant effects on the excited-state potential energy surfaces of 2,4-dithiothymine and eventually on its excited-state decay mechanism. Finally, the present computational efforts contribute important mechanistic knowledge to the understanding of the photophysics of thiothymine-based photosensitizers.