From skin sensitizers to wastewater: the unknown photo-deactivation process of low-lying excited states of isothiazolinones. A non-adiabatic dynamics investigation

Isothiazolinones represent a class of heterocyclic compounds widely used in various applications, including as biocides in cosmetics, detergents, and paints, as well as in industrial wastewater treatment. Indeed, the presence of isothiazolinones in the environment and their associated potential health hazards have raised significant concerns. In this study, a non-adiabatic dynamics investigation was conducted using state-of-the-art methodologies to explore the photochemistry of isothiazolinones. A simplified model, isothiazol-3(2H)-one (ISO), was employed to represent this compound class. The study validated the model and demonstrated that ISO can return to its ground state through the cleavage of the S-N or S-C bonds, with no significant energy barrier observed. Non-adiabatic dynamics simulations provided insights into the time scales and detailed processes of isothiazolinone photodissociation. The preferred route for deactivation was found to be the cleavage of the S-N bond. This research enhances our understanding of the photodeactivation processes of isothiazolinones and their potential environmental impact.

Ribose and Deoxyribose Group Alter Excited-State Dynamics of 5-Azacytosine in Solution

5-Azacytosine (5-AC) is one of the best interesting noncanonical nucleobases due to its functionalization and structural imitation of natural bases. 5-AC can be used as the scaffold of two important chemotherapeutic medicines, 5-azacytidine and 2'-deoxy-5-azacytidine. Furthermore, increased sensitivity to UV leads to the photochemical effects of 5-AC also attracted attention. Yet, no study has been reported to explore the effect of glycosyl groups on the photophysical and photochemical properties of 5-AC, which can help to reveal the photostability of related actual clinic drugs. In this study, the excited-state dynamics of 5-azacytidine and 2'-deoxy-5-azacytidine are studied by femtosecond transient absorption and quantum-chemical calculations while revisiting that of 5-AC with a wider probe spectral range. It is shown that glycosyl substitution on the N1 position leads to ultrafast excited-state relaxation within several picoseconds in both nucleosides, which is distinct compared with the 17 ps lifetime seen in 5-AC. It is proposed that these changes are due to altering the energy level of the dark nπ* state. Moreover, our results suggest that it should be cautioned to simply replace sugar groups with methyl groups when doing a theoretical calculation study on nucleobases and their derivatives.

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.

Vertical Excitation Energies and Lifetimes of the Two Lowest Singlet Excited States of Cytosine, 5-Aza-cytosine, and the Triazine Family: Quantum Mechanics-Molecular Mechanics Studies

A swarm of semi-classical quantum mechanics/molecular mechanics molecular-dynamics simulations where OM2/MNDO is combined with the Gromacs program for consideration of explicit water is performed, solving the time-dependent Schrödinger equation in each step of the trajectories together with the Tully's fewest switches algorithm. Within this stochastic treatment, time dependent probabilities of the three lowest electronic states are determined. The fact that nucleobases are quickly deactivated is confirmed in the cytosine case where our best lifetime estimation is τ₁=0.82 ps for the model with 100 water molecules with the SPCE force field and a time step of 0.1 fs. Lifetimes of the remaining molecules are visibly longer: 5-azacytosine, 2,4-diamino-1,3,5-triazine (DT), and 2,4,6-triamino-1,3,5-triazine (TT) molecules have an S₁ → S₀ de-excitation time of slightly above 10 ps. The lifetimes of the triazine family increases with the increasing number of exocyclic amino groups, that is, s-triazine < 2-amino-1,3,5-triazine < DT < TT. This can be explained by a higher mobility of the carbon-bonded hydrogen atoms in comparison with heavier amino groups since their movement is slowed down due to a substantially higher mass than hydrogen atoms, which can easier reach the out-of-plane positions required in the conical intersection structures. Moreover, bulkier NH₂ ligands suffer due to greater friction caused by the surrounding water environment. These mechanical aspects caused a change in the explored lifetime dependences in comparison with our previous gas-phase study.

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.

Surface hopping dynamics reveal ultrafast triplet generation promoted by S1-T2-T1 spin-vibronic coupling in 2-mercaptobenzothiazole

We investigate the T₁ formation upon populating the optically "bright" S₂ in 2-mercaptobenzothiazole to interpret the underlying relaxation pathways associated with the experimental decay constants reported by D. Koyama and A. J. Orr-Ewing, Phys. Chem. Chem. Phys., 2016, 18, 26224-26235. Energetics, electronic populations and geometries of various stationary points of low-lying electronic states are computed using the semi-classical ab initio surface hopping dynamics simulations. Estimated decay constants of S₂-S₁ internal conversion (IC) and S₁-T₂ intersystem crossing (ISC) are in excellent agreement with the experiment. The observed ultrafast ISC is analyzed based on the S₁-T₂-T₁ spin-vibronic coupling mechanism. In contrast to the previous assignment of 6 ps to the T₂-T₁ IC, our findings enable us to attribute this decay constant to the combined events of T₂-T₁ IC followed by relaxation of vibrationally hot T₁.

Solvent-Dependent Stabilization of a Charge Transfer State is the Key to Ultrafast Triplet State Formation in an Epigenetic DNA Nucleoside

2'-Deoxy-5-formylcytidine (5fdCyd), a naturally occurring nucleoside found in mammalian DNA and mitochondrial RNA, exhibits important epigenetic functionality in biological processes. Because it efficiently generates triplet excited states, it is an endogenous photosensitizer capable of damaging DNA, but the intersystem crossing (ISC) mechanism responsible for ultrafast triplet state generation is poorly understood. In this study, time-resolved mid-IR spectroscopy and quantum mechanical calculations reveal the distinct ultrafast ISC mechanisms of 5fdCyd in water versus acetonitrile. Our experiment indicates that in water, ISC to triplet states occurs within 1 ps after 285 nm excitation. PCM-TD-DFT computations suggest that this ultrafast ISC is mediated by a singlet state with significant cytosine-to-formyl charge-transfer (CT) character. In contrast, ISC in acetonitrile proceeds via a dark ¹ nπ* state with a lifetime of ∼3 ps. CT-induced ISC is not favored in acetonitrile because reaching the minimum of the gateway CT state is hampered by intramolecular hydrogen bonding, which enforces planarity between the aldehyde group and the aromatic group. Our study provides a comprehensive picture of the non-radiative decay of 5fdCyd in solution and new insights into the factors governing ISC in biomolecules. We propose that the intramolecular CT state observed here is a key to the excited-state dynamics of epigenetic nucleosides with modified exocyclic functional groups, paving the way to study their effects in DNA strands.

Enabling a Unified Description of Both Internal Conversion and Intersystem Crossing in Formaldehyde: A Global Coupled Quasi-Diabatic Hamiltonian for Its S0, S1, and T1 States

In our recent work, a diabatic Hamiltonian that couples the S₀ and S₁ states of formaldehyde was constructed using a robust fitting-and-diabatizing procedure with artificial neural networks, which is capable of representing adiabatic energies, energy gradients, and derivative couplings over a wide range of geometries including seams of conical intersection. In this work, based on the diabatization of S₀ and S₁, the spin-orbit couplings between singlet states (S₀, S₁) and triplet state T₁ are also determined in the same diabatic representation. The diabatized spin-orbit couplings are then fit with a symmetrized neural-network functional form. The ab initio spin-orbit couplings are well reproduced in large configuration space. Together with the neural-network-based potential energy surface for T₁, the full quasi-diabatic Hamiltonian for the S₀, S₁, and T₁ states is completed, enabling a unified description of both internal conversion and intersystem crossing in formaldehyde. The vibrational levels on the three adiabatic states are found to be in good agreement with known experimental band origins.

Behind the scenes of spin-forbidden decay pathways in transition metal complexes

The interpretation of the ultrafast photophysics of transition metal complexes following photo-absorption is quite involved as the heavy metal center leads to a complicated and entangled singlet-triplet manifold. This opens up multiple pathways for deactivation, often with competitive rates. As a result, intersystem crossing (ISC) and phosphorescence are commonly observed in transition metal complexes. A detailed understanding of such an excited-state structure and dynamics calls for state-of-the-art experimental and theoretical methodologies. In this review, we delve into the inability of non-relativistic quantum theory to describe spin-forbidden transitions, which can be overcome by taking into account spin-orbit coupling, whose importance grows with increasing atomic number. We present the quantum chemical theory of phosphorescence and ISC together with illustrative examples. Finally, a few applications are highlighted, bridging the gap between theoretical studies and experimental applications, such as photofunctional materials.

Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces

Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.

Photochemical Relaxation Pathways of 9 H-8-Azaguanine and 8 H-8-Azaguanine

The photochemical reaction path approach, the MS-CASPT2 quantum-chemical method, and double-ζ basis sets (cc-pVDZ) were used to study 9 H-8-azaguanine and 8 H-8-azaguanine relaxation pathways. Several potential energy hypersurfaces were characterized by means of equilibrium geometries, surface crossings (conical intersections and singlet-triplet intersystem crossings), minimum energy paths, and linear interpolation in internal coordinates. The 9 H-8-azaguanine main photochemical event begins with the direct population of the ¹(ππ* L_a) state, which evolves toward a conical intersection with the ground state after surmounting a small energy barrier, explaining why it is nonfluorescent. For 8 H-8-azaguanine, two relaxation mechanisms are possible, depending on the excitation energy. If the S₁ ¹(ππ*) state is initially populated (lower-energy region), the system evolves barrierless to the S₁ ¹(ππ*)_min region, from where the excess energy is released. If the ¹(ππ* L_a) state is populated (higher-energy radiation), the system will encounter conical intersections with the S₂ ¹(n_Oπ*) and S₁ ¹(ππ*) states before evolving to the ¹(ππ* L_a)_min region, from where a conical intersection with the ground state is accessible, favoring radiationless deactivation to the ground state. However, because a fraction of the population can be transferred from ¹(ππ* L_a) to the S₁ ¹(ππ*) state, emission from the S₁ ¹(ππ*)_min region is also expected, although weaker than it would be if the S₁ ¹(ππ*) state were populated directly. Irrespective of the excitation energy, the emissive state is the same and a single fluorescence band is observed, with the strongest emission occurring upon excitation in the lower-energy region, as observed experimentally. Therefore, our computational study corroborates experimental results, attributing the emission of the neutral form of 8-azaguanine in solution to the presence of the minor 8 H-8-azaguanine tautomer, while the 9 H-8-azaguanine major tautomer is nonfluorescent.

Conical intersections and intersystem crossings explain product formation in photochemical reactions of aryl azides

Photochemistry of substituted aryl azides is governed by surface crossings. Internal conversion and intersystem crossing govern photodecomposition of 3-methoxyphenyl azide and 4-methoxyphenyl azide.

Highly efficient surface hopping dynamics using a linear vibronic coupling model

We report an implementation of the linear vibronic coupling (LVC) model within the surface hopping dynamics approach and present utilities for parameterizing this model in a blackbox fashion. This results in an extremely efficient method to obtain qualitative and even semi-quantitative information about the photodynamical behavior of a molecule, and provides a new route toward benchmarking the results of surface hopping computations. The merits and applicability of the method are demonstrated in a number of applications. First, the method is applied to the SO2 molecule showing that it is possible to compute its absorption spectrum beyond the Condon approximation, and that all the main features and timescales of previous on-the-fly dynamics simulations of intersystem crossing are reproduced while reducing the computational effort by three orders of magnitude. The dynamics results are benchmarked against exact wavepacket propagations on the same LVC potentials and against a variation of the electronic structure level. Four additional test cases are presented to exemplify the broader applicability of the model. The photodynamics of the isomeric adenine and 2-aminopurine molecules are studied and it is shown that the LVC model correctly predicts ultrafast decay in the former and an extended excited-state lifetime in the latter. Futhermore, the method correctly predicts ultrafast intersystem crossing in the modified nucleobase 2-thiocytosine and its absence in 5-azacytosine while it fails to describe the ultrafast internal conversion to the ground state in the latter.

Nonadiabatic dynamics: The SHARC approach

We review the Surface Hopping including ARbitrary Couplings (SHARC) approach for excited-state nonadiabatic dynamics simulations. As a generalization of the popular surface hopping method, SHARC allows simulating the full-dimensional dynamics of molecules including any type of coupling terms beyond nonadiabatic couplings. Examples of these arbitrary couplings include spin-orbit couplings or dipole moment-laser field couplings, such that SHARC can describe ultrafast internal conversion, intersystem crossing, and radiative processes. The key step of the SHARC approach consists of a diagonalization of the Hamiltonian including these couplings, such that the nuclear dynamics is carried out on potential energy surfaces including the effects of the couplings-this is critical in any applications considering, for example, transition metal complexes or strong laser fields. We also give an overview over the new SHARC2.0 dynamics software package, released under the GNU General Public License, which implements the SHARC approach and several analysis tools. The review closes with a brief survey of applications where SHARC was employed to study the nonadiabatic dynamics of a wide range of molecular systems. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and KineticsSoftware > Simulation MethodsSoftware > Quantum Chemistry.

Surface hopping dynamics including intersystem crossing using the algebraic diagrammatic construction method

We report an implementation for employing the algebraic diagrammatic construction to second order [ADC(2)] ab initio electronic structure level of theory in nonadiabatic dynamics simulations in the framework of the SHARC (surface hopping including arbitrary couplings) dynamics method. The implementation is intended to enable computationally efficient, reliable, and easy-to-use nonadiabatic dynamics simulations of intersystem crossing in organic molecules. The methodology is evaluated for the 2-thiouracil molecule. It is shown that ADC(2) yields reliable excited-state energies, wave functions, and spin-orbit coupling terms for this molecule. Dynamics simulations are compared to previously reported results using high-level multi-state complete active space perturbation theory, showing favorable agreement.

An MS-CASPT2 study of the photodecomposition of 4-methoxyphenyl azide: role of internal conversion and intersystem crossing

Photoexcitation of 4-methoxyphenyl azide at 266 nm yields triplet nitrene after decaying through an intersystem crossing (ISC1, 2¹A′/2³A′′) in a first step followed by internal conversion (CI2, 2³A′′/1³A′′).