QM/MM Studies on the Photophysical Mechanism of a Truncated Octocrylene Model

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.

Quantum mechanics/molecular mechanics studies on mechanistic photophysics of cytosine aza-analogues: 2,4-diamino-1,3,5-triazine and 2-amino-1,3,5-triazine in aqueous solution

The excited-state properties and photophysics of cytosine aza-analogues, i.e., 2,4-diamino-1,3,5-triazine (2,4-DT) and 2-amino-1,3,5-triazine (2-AT) in solution have been systematically explored using the QM(MS-CASPT2//CASSCF)/MM approach. The excited-state nonradiative relaxation mechanisms for the initially photoexcited S₁(ππ*) state decay back to the S₀ state are proposed in terms of the present computed minima, surface crossings (conical intersections and singlet-triplet crossings), and excited-state decay paths in the S₁, S₂, T₁, T₂, and S₀ states. Upon photoexcitation to the bright S₁(ππ*) state, 2,4-DT quickly relaxes to its S₁ minimum and then overcomes a small energy barrier of 5.1 kcal mol^-1 to approach a S₁/S₀ conical intersection, where the S₁ system hops to the S₀ state through S₁ → S₀ internal conversion (IC). In addition, at the S₁ minimum, the system could partially undergo intersystem crossing (ISC) to the T₁ state, followed by further ISC to the S₀ state via the T₁/S₀ crossing point. In the T₁ state, an energy barrier of 7.9 kcal mol^-1 will trap 2,4-DT for a while. In parallel, for 2-AT, the system first relaxes to the S₁ minimum and then S₁ → S₀ IC or S₁ → T₁ → S₀ ISCs take place to the S₀ state by surmounting a large barrier of 15.3 kcal mol^-1 or 11.9 kcal mol^-1, respectively, which heavily suppress electronic transition to the S₀ state. Different from 2,4-DT, upon photoexcitation in the Franck-Condon region, 2-AT can quickly evolve in an essentially barrierless manner to nearby S₂/S₁ conical intersection, where the S₂ and T₁ states can be populated. Once it hops to the S₂ state, the system will overcome a relatively small barrier (6.6 kcal mol^-1vs. 15.3 kcal mol^-1) through IC to the S₀ state. Similarly, an energy barrier of 11.9 kcal mol^-1 heavily suppresses the T₁ state transformation to the S₀ state. The present work manifests that the amination/deamination of the triazine rings can affect some degree of different vertical and adiabatic excitation energies and nonradiative decay pathways in solution. It not only rationalizes excited-state decay dynamics of 2,4-DT and 2-AT in aqueous solution but could also provide insights into the understanding of the photophysics of aza-nucleobases.

A theoretical insight into excited-state decay and proton transfer of p-nitrophenylphenol in the gas phase and methanol solution

The excited-state decay (ESD) and proton transfer (EPT) of p-nitrophenylphenol (NO₂-Bp-OH), especially in the triplet states, were not characterized with high-level theoretical methods to date. Herein, the MS-CASPT2//CASSCF and QM(MS-CASPT2//CASSCF)/MM methods were employed to gain an atomic-level understanding of the ESD and EPT of NO₂-Bp-OH in the gas phase and its hydrogen-bonded complex in methanol. Our calculation results revealed that the S₁ and S₂ states of NO₂-Bp-OH are of ¹ππ* and ¹nπ* characters at the Franck-Condon (FC) point, which correspond to the ICT-EPT and intramolecular charge-transfer (ICT) states in spectroscopic experiments. The former state has a charge-transfer property that could facilitate the EPT reaction, while the latter one might be unfavorable for EPT. The vertical excitation energies of these states are almost degenerate at the FC region and the electronic configurations of ¹ππ* and ¹nπ* will exchange from the S₁ FC region to the S₁ minimum, which means that the ¹nπ* state will participate in ESD once NO₂-Bp-OH departs from the S₁ FC region. Besides, we found that three triplets lie below the first bright state and will play very important roles in intersystem crossing processes. In terms of several pivotal surface crossings and relevant linearly interpolated internal coordinate (LIIC) paths, three feasible but competing ESD channels that could effectively lead the system to the ground state or the lowest triplet state were put forward. Once arrived at the T₁ state, the system has enough time and internal energy to undergo the EPT reaction. The methanol solvent has a certain effect on the relative energies and spin-orbit couplings, but does not qualitatively change the ESD processes of NO₂-Bp-OH. By contrast, the solvent effects will remarkably stabilize the proton-transferred product by the hydrogen bond networks and assist to form the triplet anion. Our present work would pave the road to properly understand the mechanistic photochemistry of similar hydroxyaromatic compounds.

Quantum mechanics/molecular mechanics studies on the mechanistic photophysics of sunscreen oxybenzone in methanol solution

Herein, we have employed the QM(CASPT2//CASSCF)/MM method to explore the photophysical and photochemical mechanism of oxybenzone (OB) in methanol solution. Based on the optimized minima, conical intersections and crossing points, and minimum-energy reaction paths related to excited-state intramolecular proton transfer (ESIPT) and excited-state decay paths in the ¹ππ*, ¹nπ*, ³ππ*, ³nπ*, and S₀ states, we have identified several feasible excited-state relaxation pathways for the initially populated S₂(¹ππ*) state to decay to the initial enol isomer' S₀ state. The major one is the singlet-mediated and stretch-torsion coupled ESIPT pathway, in which the system first undergoes an essentially barrierless ¹ππ* ESIPT process to generate the ¹ππ* keto species, and finally realizes its ground state recovery through the subsequent carbonyl stretch-torsion facilitating S₁ → S₀ internal conversion (IC) and the reverse ground-state intramolecular proton transfer (GSIPT) process. The minor ones are related to intersystem crossing (ISC) processes. At the S₂(¹ππ*) minimum, an S₂(¹ππ*)/S₁(¹nπ*)/T₂(³nπ*) three-state intersection region helps the S₂ system branch into the T₁ state through a S₂ → S₁ → T₁ or S₂ → T₂ → T₁ process. Once it has reached the T₁ state, the system may relax to the S₀ state via direct ISC or via subsequent nearly barrierless ³ππ* ESIPT to yield the T₁ keto tautomer and ISC. The resultant S₀ keto species significantly undergoes reverse GSIPT and only a small fraction yields the trans-keto form that relaxes back more slowly. However, due to small spin-orbit couplings at T₁/S₀ crossing points, the ISC to S₀ state occurs very slowly. The present work rationalizes not only the ultrafast excited-state decay dynamics of OB but also its phosphorescence emission at low temperature.

Substituent and Solvent Effects on the Photoisomerization of Cinnamate Derivatives: An XMS-CASPT2 Study

Cinnamate derivatives show a variety of photo-induced reactions. Among them is trans-cis photoisomerization, which may involve the nonradiative decay (NRD) process. The extended multistate complete active space second-order perturbation (XMS-CASPT2) method was used in this study as a suitable theory for treating multireference electronic nature, which was frequently manifested in the photoisomerization process. The minimum energy paths of the trans-cis photoisomerization process of cinnamate derivatives were determined, and the activation energies were estimated using the resulting intrinsic reaction coordinate (IRC) paths. Natural orbital analysis revealed that the transition state's (TS) electronic structure is zwitterionic-like, elucidating the solvent and substituent effect on the energy barrier of photoisomerization paths. Furthermore, it was found that the charge on the pyramidalized carbon atom at the TS structure was strongly correlated with the activation energy barrier for the cinnamate derivatives. Thus, it seemingly provided a physical picture of the photoisomerization of cinnamates and was a good descriptor potentially applicable to molecular design for controlling the rate constant of the photoisomerization reaction.

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" S₁(¹ππ*) state and the lower-lying T₁ and T₂ states. In the major relaxation pathway, the photoexcited S₁ keto system first undergoes an essentially barrierless excited-state intramolecular proton transfer (ESIPT) to generate the S₁ enol minimum, near which a favorable S₁/S₀ conical intersection decays the system to the S₀ state followed by a reverse ground-state intramolecular proton transfer (GSIPT) to repopulate the initial S₀ keto species. In the minor one, an S₁/T₂/T₁ three-state intersection in the keto region makes the T₁ state populated via direct and T₂-mediated intersystem crossing (ISC) processes. In the T₁ state, an ESIPT occurs, which is followed by ISC near a T₁/S₀ crossing point in the enol region to the S₀ state and finally back to the S₀ keto species. In addition, a T₁/S₀ crossing point near the T₁ keto minimum can also help the system decay to the S₀ keto species. However, small spin-orbit couplings between T₁ and S₀ at these T₁/S₀ crossing points make ISC to the S₀ state very slow and make the system trapped in the T₁ 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.

Mechanistic photophysics and photochemistry of unnatural bases and sunscreen molecules: insights from electronic structure calculations

Photophysics and photochemistry are basic subjects in the study of light-matter interactions and are ubiquitous in diverse fields such as biology, energy, materials, and environment. A full understanding of mechanistic photophysics and photochemistry underpins many recent advances and applications. This contribution first provides a short discussion on the theoretical calculation methods we have used in relevant studies, then we introduce our latest progress on the mechanistic photophysics and photochemistry of two classes of molecular systems, namely unnatural bases and sunscreens. For unnatural bases, we disclose the intrinsic driving forces for the ultrafast population to reactive triplet states, impacts of the position and degree of chalcogen substitutions, and the effects of complex environments. For sunscreen molecules, we reveal the photoprotection mechanisms that dissipate excess photon energy to the surroundings by ultrafast internal conversion to the ground state. Finally, relevant theoretical challenges and outlooks are discussed.

Quantum Mechanics/Molecular Mechanics Studies on the Photophysical Mechanism of Methyl Salicylate

Methyl salicylate (MS) as a subunit of larger salicylates found in commercial sunscreens has been shown to exhibit keto-enol tautomerization and dual fluorescence emission via excited-state intramolecular proton transfer (ESIPT) after the absorption of ultraviolet (UV) radiation. However, its excited-state relaxation mechanism is unclear. Herein, we have employed the quantum mechanics(CASPT2//CASSCF)/molecular mechanics method to explore the ESIPT and excited-state relaxation mechanism of MS in the lowest three electronic states, that is, S₀, S₁, and T₁ states, in a methanol solution. Based on the optimized geometric and electronic structures, conical intersections and crossing points, and minimum-energy paths combined with the computed linearly interpolated Cartesian coordinate paths, the photophysical mechanism of MS has been proposed. The S₁ state is a spectroscopically bright ¹ππ* state in the Franck-Condon region. From the initially populated S₁ state, there exist three nonradiative relaxation paths to repopulate the S₀ state. In the first one, the S₁ system (i.e., ketoB form) first undergoes an ESIPT path to generate an S₁ tautomer (i.e., enol form) that exhibits a large Stokes shift in experiments. The generated S₁ enol tautomer further evolves toward the nearby S₁/S₀ conical intersection and then hops to the S₀ state, followed by the backward ground-state intramolecular proton transfer (GSIPT) to the initial ketoB form S₀ state. In the second one, the S₁ system first hops through the S₁ → T₁ intersystem crossing (ISC) to the T₁ state, which then further decays to the S₀ state via T₁ → S₀ ISC at the T₁/S₀ crossing point. In the third path, the T₁ system that stems from the S₁ → T₁ ISC process via the S₁/T₁ crossing point first takes place a T₁ ESIPT to generate a T₁ enol tautomer, which can further decay to the S₀ state via T₁-to-S₀ ISC. Finally, the GSIPT occurs to back the system to the initial ketoB form S₀ state. Our present work could contribute to understanding the photophysics of MS and its derivatives.

Theoretical studies on the photochemistry of 2-nitrofluorene in the gas phase and acetonitrile solution

The photophysical and photochemical mechanisms of 2-nitrofluorene (2-NF) in the gas phase and acetonitrile solution have been studied theoretically. Upon ∼330 nm irradiation to the first bright state (¹ππ*), the 2-NF system can decay to triplet excited states via rapid intersystem crossing (ISC) processes through different surface crossing points or to the ground state via an ultrafast internal conversion (IC) process through the S₁/S₀ conical intersection. The ¹nπ* dark state will serve as a bridge when the system leaves the Franck-Condon (FC) region and approaches to the S₁ minimum. The molecule maintains a planar geometry during the excited-state relaxation processes. The differences on excitation properties such as electronic configurations and spin-orbit coupling (SOC) interactions between those in the gas phase and acetonitrile solution cannot be neglected, indicating possible changes on the efficiency of the related ISC processes for the 2-NF system in solution. Once arrived at the T₁ state, it would further decay to the S₀ state or photodegrade into the Ar-O˙ and NO˙ free radicals. During the intramolecular rearrangement process, the twisting of the nitro group out of the aromatic-ring plane is regarded as a critical structural variation for the photodegradation of the 2-NF system. The free radicals finally form through oxaziridine-type intermediate and transition state structures. The present work provides important mechanistic insights to the photochemistry of nitro-substituted polyaromatic compounds.

Ultrafast Nonadiabatic Photoisomerization Dynamics Mechanism for the UV Photoprotection of Stilbenoids in Grape Skin

Natural UV photoprotection plays a vital role in physiological protection. It has been reported that UVC radiation can make resveratrol (RSV) and piceatannol (PIC) accumulate in grape skin. In this work, we demonstrated that RSV and PIC could significantly absorb UVA and UVB, and confirmed their satisfactory photostability. Furthermore, we clarified the UV photoprotection mechanism of typical stilbenoids of RSV and PIC for the first time by using combined femtosecond transient absorption (FTA) spectroscopy and time-dependent density functional theory (TD-DFT) calculations. RSV and PIC can be photoexcited to the excited state after UVA and UVB absorption. Subsequently, the photoisomerized RSV and PIC quickly relax to the ground state via nonadiabatic transition from the S₁ state at a conical intersection (CI) position between potential energy surfaces (PESs) of S₁ and S₀ states. This ultrafast trans-cis photoisomerization will take place within a few tens of picoseconds. As a result, the UV energy absorbed by RSV and PIC could be dissipated by an ultrafast nonadiabatic photoisomerization process.