Photodissociation dynamics of CH3C(O)SH in argon matrix: A QM/MM nonadiabatic dynamics simulation

Generalized Ab Initio Nonadiabatic Dynamics Simulation Methods from Molecular to Extended Systems

Nonadiabatic dynamics simulation has become a powerful tool to describe nonadiabatic effects involved in photophysical processes and photochemical reactions. In the past decade, our group has developed generalized trajectory-based ab initio surface-hopping (GTSH) dynamics simulation methods, which can be used to describe a series of nonadiabatic processes, such as internal conversion, intersystem crossing, excitation energy transfer and charge transfer of molecular systems, and photoinduced nonadiabatic carrier dynamics of extended systems with and without spin-orbit couplings. In this contribution, we will first give a brief introduction to our recently developed methods and related numerical implementations at different computational levels. Later, we will present some of our latest applications in realistic systems, which cover organic molecules, biological proteins, organometallic compounds, periodic organic and inorganic materials, etc. Final discussion is given to challenges and outlooks of ab initio nonadiabatic dynamics simulations.

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.

QM/MM nonadiabatic dynamics simulations on photoinduced Wolff rearrangements of 1,2,3-thiadiazole

The photoinduced rearrangement reaction mechanism of 1,2,3-thiadiazole remains experimentally elusive. Two possible mechanisms have been proposed to date. The first is a stepwise mechanism via a thiocarbene intermediate; the second is an excited-state concerted rearrangement mechanism. Herein we have adopted both the electronic structure calculations and nonadiabatic dynamics simulations to study the photoinduced rearrangement reactions of 1,2,3-thiadiazole in the S₂, S₁, and S₀ states in solution. On the basis of QM(CASPT2)/MM [quantum mechanics(complete active space self-consistent field second-order perturbation theory)/molecular mechanics] calculations, we have found that (1) the thiocarbene intermediate is not stable; thus, the stepwise mechanism should be unfavorable; (2) the excited-state decay from the S₂ via S₁ to S₀ state is ultrafast and completed within ca. 200 fs; therefore, both the S₂ and S₁ states should not have a long enough time for the excited-state rearrangements. Instead, we have computationally proposed a modified photoinduced rearrangement mechanism. Upon irradiation, the S₂ state is first populated (114.0 kcal/mol), followed by an ultrafast S₂ → S₁ → S₀ excited-state decay along the S-N bond fission, which eventually leads to a very "hot" intermediate with the S-N bond broken (18.3 kcal/mol). Then, thermal rearrangements to thioketene, thiirene, and ethynethiol occur in a concerted asynchronous way. This mechanistic scenario has been verified by full-dimensional trajectory-based nonadiabatic dynamics simulations at the QM(CASPT2)/MM level. Finally, our present computational work provides experimentally interesting mechanistic insights into the photoinduced rearrangement reactions of cyclic and acyclic diazo compounds.

Electronic structure calculations and nonadiabatic dynamics simulations of excited-state relaxation of Pigment Yellow 101

Pigment Yellow 101 (PY101) is widely used as a typical pigment due to its excellent excited-state properties. However, the origin of its photostability is still elusive. In this work, we have systematically investigated the photodynamics of PY101 by performing combined electronic structure calculations and trajectory-based nonadiabatic dynamics simulations. On the basis of the results, we have found that upon photoexcitation to the S₁ state, PY101 undergoes an essentially barrierless excited-state intramolecular single proton transfer generating an S₁ keto species. In the keto region, there is an energetically accessible S₁/S₀ conical intersection that funnels the system to the S₀ state quickly. In the S₀ state, the keto species either goes back to its trans-enol species through a ground-state reverse hydrogen transfer or arrives at the cis-keto region. In addition, we have found an additional excited-state decay channel for the S₁ enol species, which is directly linked to an S₁/S₀ conical intersection located in the enol region. This mechanism has also been confirmed by our dynamics simulations, in which about 54% of the trajectories decay to the S₀ state via the enol S₁/S₀ conical intersection; while the remaining ones employ the keto S₁/S₀ conical intersection. The gained mechanistic information helps us understand the photostability of the PY101 chromophore and its variants with the same molecular scaffold.

Photoinduced Curtius rearrangements of fluorocarbonyl azide, FC(O)N3: a QM/MM nonadiabatic dynamics simulation

Upon either photolysis or pyrolysis, carbonyl azide can eliminate molecular nitrogen along with the formation of isocyanate.

A theoretical study of the light-induced cross-linking reaction of 5-fluoro-4-thiouridine with thymine

In contrast to photophysics of thio-substituted nucleobases, their photoinduced cross-linking reactions with canonical nucleobases remain scarcely investigated computationally. In this work, we have adopted combined CASPT2/PCM//CASSCF and B3LYP-D3/PCM electronic structure methods to study this kind of photochemical reaction of 5-fluoro-4-thiouridine (truncated 5-fluoro-1-methyl-4-thiouracil used in calculations) and 1-methylthymine (referred to as thymine for clarity hereinafter). On the basis of CASPT2/PCM computed results, we have proposed two efficient excited-state relaxation pathways to populate the lowest T₁ state of the complex of 5-fluoro-1-methyl-4-thiouracil and thymine from its initially populated S₂(¹ππ*) state. In the first one, the S₂ system first hops to the S₁ state via an S₂/S₁ conical intersection, followed by a direct S₁ → T₁ intersystem crossing process enhanced by large S₁/T₁ spin-orbit coupling. In the second path, the resultant S₁ system first jumps to the T₂ state, from which an efficient T₂ → T₁ internal conversion occurs. The T₁ cross-linking reaction is overall divided into two phases. The first phase is a stepwise and nonadiabatic photocyclization reaction, which starts from the T₁ complex and ends up with an S₀ thietane intermediate. The second phase is a thermal reaction. The system first rearranges its four- and six-membered rings to form three new rings; then, an S₀ fluorine atom transfer occurs, followed by the formation of photoproducts. Finally, the present work paves the way for studying light-induced cross-linking reactions of thionucleobases with canonical bases in DNA and RNA.

Quantum Chemical Investigation on Photochemical Reactions of Nonanoic Acids at Air-Water Interface

Photoinduced chemical reactions of organic compounds at the marine boundary layer have recently attracted significant experimental attention because this kind of photoreactions has been proposed to have substantial impact on local new particle formation and their photoproducts could be a source of secondary organic aerosols. In this work, we have employed first-principles density functional theory method combined with cluster models to systematically explore photochemical reaction pathways of nonanoic acids (NAs) to form volatile saturated and unsaturated C₉ and C₈ aldehydes at air-water interfaces. On the basis of the results, we have found that the formation of C₉ aldehydes is not initiated by intermolecular Norrish type II reaction between two NAs but by intramolecular T₁ C-O bond fission of NA generating acyl and hydroxyl radicals. Subsequently, saturated C₉ aldehydes are formed through hydrogenation reaction of acyl radical by another intact NA. Following two dehydrogenation reactions, unsaturated C₉ aldehydes are generated. In parallel, the pathway to C₈ aldehydes is initiated by T₁ C-C bond fission of NA, which generates octyl and carboxyl radicals; then, an octanol is formed through recombination reaction of octyl with hydroxyl radical. In the following, two dehydrogenation reactions result into an enol intermediate from which saturated C₈ aldehydes are produced via NA-assisted intermolecular hydrogen transfer. Finally, two dehydrogenation reactions generate unsaturated C₈ aldehydes. In these reactions, water and NA molecules are found to play important roles. They significantly reduce relevant reaction barriers. Our work has also explored oxygenation reactions of NA with molecular oxygen and radical-radical dimerization reactions.

Photodynamics of oxybenzone sunscreen: Nonadiabatic dynamics simulations

Herein we have used combined static electronic structure calculations and "on-the-fly" global-switching trajectory surface-hopping dynamics simulations to explore the photochemical mechanism of oxybenzone sunscreen. We have first employed the multi-configurational CASSCF method to optimize minima, conical intersections, and minimum-energy reaction paths related to excited-state intramolecular proton transfer (ESIPT) and excited-state decays in the (1)ππ(∗), (1)nπ(∗), and S0 states (energies are refined at the higher MS-CASPT2 level). According to the mapped potential energy profiles, we have identified two ultrafast excited-state deactivation pathways for the initially populated (1)ππ(∗) system. The first is the diabatic ESIPT process along the (1)ππ(∗) potential energy profile. The generated (1)ππ(∗) keto species then decays to the S0 state via the keto (1)ππ(∗)/gs conical intersection. The second is internal conversion to the dark (1)nπ(∗) state near the (1)ππ(∗) /(1)nπ(∗) crossing point in the course of the diabatic (1)ππ(∗) ESIPT process. Our following dynamics simulations have shown that the ESIPT and (1)ππ(∗) → S0 internal conversion times are 104 and 286 fs, respectively. Finally, our present work demonstrates that in addition to the ESIPT process and the (1)ππ(∗) → S0 internal conversion in the keto region, the (1)ππ(∗) → (1)nπ(∗) internal conversion in the enol region plays as well an important role for the excited-state relaxation dynamics of oxybenzone.

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

We have employed computational methods to systematically explore the excited-state decay mechanism of 2,4-dithiothymine.

Intersystem crossing-branched excited-state intramolecular proton transfer for o-nitrophenol: An ab initio on-the-fly nonadiabatic molecular dynamic simulation

The 6SA-CASSCF(10, 10)/6-31G (d, p) quantum chemistry method has been applied to perform on-the-fly trajectory surface hopping simulation with global switching algorithm and to explore excited-state intramolecular proton transfer reactions for the o-nitrophenol molecule within low-lying electronic singlet states (S0 and S1) and triplet states (T1 and T2). The decisive photoisomerization mechanisms of o-nitrophenol upon S1 excitation are found by three intersystem crossings and one conical intersection between two triplet states, in which T1 state plays an essential role. The present simulation shows branch ratios and timescales of three key processes via T1 state, non-hydrogen transfer with ratio 48% and timescale 300 fs, the tunneling hydrogen transfer with ratios 36% and timescale 10 ps, and the direct hydrogen transfer with ratios 13% and timescale 40 fs. The present simulated timescales might be close to low limit of the recent experiment results.