Conical Intersection Is Responsible for the Fluorescence Disappearance below 365 nm in Cyclopropanone

Electronic dynamics through conical intersections via non-Markovian stochastic Schrödinger equation with complex modes

Conical intersections (CIs) play a crucial role in photochemical reactions, offering an efficient channel for ultrafast non-adiabatic relaxation of excited states. This significantly influences the reaction pathways and the resulting products. In this work, we utilize the non-Markovian stochastic Schrödinger equation with complex modes method to explore the dynamics of electronic transitions through conical intersections (CIs) in pyrazine. The linear vibronic coupling model serves as the foundational framework, incorporating both intra-state and inter-state electron-vibrational interactions. The dynamics of the excited electronic transitions are analyzed across varying strengths of system-bath coupling and different bath relaxation times. The accuracy of this method is demonstrated by comparing its predictions with those from the hierarchical equations of motion method.

A Local Diabatisation Method for Two-State Adiabatic Conical Intersections

A methodology to locally characterize conical intersections (CIs) between two adiabatic electronic states for which no nonadiabatic coupling (NAC) vectors are available is presented. Based on the Hessian and gradient at the CI, the branching space coordinates are identified. The potential energy surface around the CI in the branching space is expressed in the diabatic representation, from which the NAC vectors can be calculated in a wave-function-free, energy-based approach. To demonstrate the universality of the developed methodology, the minimum-energy CI (MECI) between the first (S1) and second (S2) singlet excited states of formamide is investigated at the state-averaged complete active space self-consistent field (SA-CASSCF) and extended multistate complete active space second-order perturbation theory (XMS-CASPT2) levels of theory. In addition, the asymmetrical MECI between the ground state (S0) and S1 of cyclopropanone is evaluated using SA-CASSCF, as well as (ME)CIs between the S1 and S2 states of benzene using SA-CASSCF and time-dependent density functional theory (TDDFT). Finally, a CI between the S1 and S2 excited states of thiophene was analyzed using TDDFT.

What Controls the Quality of Photodynamical Simulations? Electronic Structure Versus Nonadiabatic Algorithm

The field of nonadiabatic dynamics has matured over the last decade with a range of algorithms and electronic structure methods available at the moment. While the community currently focuses more on developing and benchmarking new nonadiabatic dynamics algorithms, the underlying electronic structure controls the outcome of nonadiabatic simulations. Yet, the electronic-structure sensitivity analysis is typically neglected. In this work, we present a sensitivity analysis of the nonadiabatic dynamics of cyclopropanone to electronic structure methods and nonadiabatic dynamics algorithms. In particular, we compare wave function-based CASSCF, FOMO-CASCI, MS- and XMS-CASPT2, density-functional REKS, and semiempirical MRCI-OM3 electronic structure methods with the Landau-Zener surface hopping, fewest switches surface hopping, and ab initio multiple spawning with informed stochastic selection algorithms. The results clearly demonstrate that the electronic structure choice significantly influences the accuracy of nonadiabatic dynamics for cyclopropanone even when the potential energy surfaces exhibit qualitative and quantitative similarities. Thus, selecting the electronic structure solely on the basis of the mapping of potential energy surfaces can be misleading. Conversely, we observe no discernible differences in the performance of the nonadiabatic dynamics algorithms across the various methods. Based on the above results, we discuss the present-day practice in computational photodynamics.

The ΔSCF method for non-adiabatic dynamics of systems in the liquid phase

Computational studies of ultrafast photoinduced processes give valuable insights into the photochemical mechanisms of a broad range of compounds. In order to accurately reproduce, interpret, and predict experimental results, which are typically obtained in a condensed phase, it is indispensable to include the condensed phase environment in the computational model. However, most studies are still performed in vacuum due to the high computational cost of state-of-the-art non-adiabatic molecular dynamics (NAMD) simulations. The quantum mechanical/molecular mechanical (QM/MM) solvation method has been a popular model to perform photodynamics in the liquid phase. Nevertheless, the currently used QM/MM embedding techniques cannot sufficiently capture all solute-solvent interactions. In this Perspective, we will discuss the efficient ΔSCF electronic structure method and its applications with respect to the NAMD of solvated compounds, with a particular focus on explicit quantum mechanical solvation. As more research is required for this method to reach its full potential, some challenges and possible directions for future research are presented as well.

The photodissociation of solvated cyclopropanone and its hydrate explored via non-adiabatic molecular dynamics using ΔSCF

The decay of cyclopropanone is a typical example of a photodecomposition process. Ethylene and carbon monoxide are formed following the excitation to the first singlet excited state through a symmetrical or asymmetrical pathway. The results obtained with non-adiabatic molecular dynamics (NAMD) using the delta self-consistent field (ΔSCF) method correspond well to previous experimental and multireference theoretical studies carried out in the gas phase. Moreover, this efficient methodology allows NAMD simulations of cyclopropanone in aqueous solution to be performed, which reveal analogue deactivation mechanisms, but a shorter lifetime and reduced photodissociation as compared to the gas-phase. The excited state dynamics of cyclopropanone hydrate, an enzyme inhibitor, in an aqueous environment are reported as well. Cyclopropanone hydrate strongly interacts with the surrounding solvent via the formation of hydrogen bonds. Excitation to the first singlet excited state shows an asymmetric pathway with cyclopropanone hydrate and propionic acid as the main photoproducts.

The lifetime and photodissociation of cyclopropanone are reduced in aqueous solution, while the excitation of solvated cyclopropanone hydrate yields a range of photoproducts.

A Simple Solution to Trivial Crossings: A Stochastic State Tracking Approach

We present a new state tracking algorithm based on a stochastic state reassignment that reflects the quantum mechanical interpretation of the state time-overlaps. We assess the new method with a range of model Hamiltonians and demonstrate that it yields the results generally consistent with the deterministic min-cost algorithm. However, the stochastic state tracking algorithm reduces magnitudes of the state population fluctuations as the quantum system evolves toward its equilibrium. The new algorithm facilitates the thermalization of quantum state populations and suppresses the population revivals and oscillations near the equilibrium in many-state systems. The new stochastic algorithm has a favorable computational scaling, is easy to implement due to its conceptual transparency, and treats various types of state identity changes (trivial or avoided crossings and any intermediate cases) on equal footing.

Nonadiabatic Dynamics Simulations on Early-Time Photochemistry of Spirobenzopyran

Photoinduced ring-opening, decay, and isomerization of spirobenzopyran have been explored by the OM2/MRCI nonadiabatic dynamics simulations based on Tully's fewest-switches surface hopping scheme. The efficient S₁ to S₀ internal conversion as observed in experiments is attributed to the existence of two efficient excited-state decay pathways. The first one is related to the C-N dissociation, and the second one is done to the C-O dissociation. The C-O dissociation pathway is dominant, and more than 90% trajectories decay to the S₀ state via the C-O bond-fission related S₁/S₀ conical intersections. Near these regions in the S₀ state, trajectories can either return to spirobenzopyran or proceed to various intermediates including merocyanine via a series of bond rotations. Our nonadiabatic dynamics simulations also demonstrate that the hydrogen-out-of-plane (HOOP) motion is important for efficient and ultrafast excited-state deactivation. On the other hand, we have also found that the replacement of methyl groups by hydrogen atoms in spirobenzopyran can artificially introduce different intramolecular hydrogen transfers leading to hydrogen-transferred intermediates. This finding is important for the community and demonstrates that such a kind of structural truncation, sometimes, could be problematic, leading to incorrect photodynamics. Our present work provides valuable insights into the photodynamics of spirobenzopyran, which could be helpful for the design of spiropyran-based photochromic materials.

Theoretical modelling of the dynamics of primary photoprocess of cyclopropanone

Photodecomposition of cyclopropanone is investigated by static quantum chemical calculations and non-adiabatic molecular dynamics (NAMD) simulations.

Unraveling Photodimerization of Cyclohexasilane from Molecular Dynamics Studies

Photoinduced reactions of a pair of cyclohexasilane (CHS) monomers are explored by time-dependent excited-state molecular dynamics (TDESMD) calculations. In TDESMD trajectories, one observes vivid reaction events including dimerization and fragmentation. A general reaction pathway is identified as (i) ring-opening formation of a dimer, (ii) rearrangement induced by bond breaking, and (iii) decomposition through the elimination of small fragments. The identified pathway supports the chemistry proposed for the fabrication of silicon-based materials using CHS as a precursor. In addition, we find dimers have smaller HOMO-LUMO gaps and exhibit a red shift and line-width broadening in the computed photoluminescence spectra compared with a pair of CHS monomers.

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.

Excited-State Proton-Transfer-Induced Trapping Enhances the Fluorescence Emission of a Locked GFP Chromophore

The chemical locking of the central single bond in core chromophores of green fluorescent proteins (GFPs) influences their excited-state behavior in a distinct manner. Experimentally, it significantly enhances the fluorescence quantum yield of GFP chromophores with an ortho-hydroxyl group, while it has almost no effect on the photophysics of GFP chromophores with a para-hydroxyl group. To unravel the underlying physical reasons for this different behavior, we report static electronic structure calculations and nonadiabatic dynamics simulations on excited-state intramolecular proton transfer, cis–trans isomerization, and excited-state deactivation in a locked ortho-substituted GFP model chromophore (o-LHBI). On the basis of our previous and present results, we find that the S₁ keto species is responsible for the fluorescence emission of the unlocked o-HBI and the locked o-LHBI species. Chemical locking does not change the parts of the S₁ and S₀ potential energy surfaces relevant to enol–keto tautomerization; hence, in both chromophores, there is an ultrafast excited-state intramolecular proton transfer that takes only 35 fs on average. However, the locking effectively hinders the S₁ keto species from approaching the keto S₁/S₀ conical intersections so that most of trajectories are trapped in the S₁ keto region for the entire 2 ps simulation time. Therefore, the fluorescence quantum yield of o-LHBI is enhanced compared with that of unlocked o-HBI, in which the S₁ excited-state decay is efficient and ultrafast. In the case of the para-substituted GFP model chromophores p-HBI and p-LHBI, chemical locking hardly affects their efficient excited-state deactivation via cis–trans isomerization; thus, the fluorescence quantum yields in these chromophores remain very low. The insights gained from the present work may help to guide the design of new GFP chromophores with improved fluorescence emission and brightness.

Intramolecular hydrogen bonding plays a crucial role in the photophysics and photochemistry of the GFP chromophore

In commonly studied GFP chromophore analogues such as 4-(4-hydroxybenzylidene)-1,2-dimethyl-1H-imidazol-5(4H)-one (PHBDI), the dominant photoinduced processes are cis-trans isomerization and subsequent S(1) → S(0) decay via a conical intersection characterized by a highly twisted double bond. The recently synthesized 2-hydroxy-substituted isomer (OHBDI) shows an entirely different photochemical behavior experimentally, since it mainly undergoes ultrafast intramolecular excited-state proton transfer, followed by S(1) → S(0) decay and ground-state reverse hydrogen transfer. We have chosen 4-(2-hydroxybenzylidene)-1H-imidazol-5(4H)-one (OHBI) to model the gas-phase photodynamics of such 2-hydroxy-substituted chromophores. We first use various electronic structure methods (DFT, TDDFT, CC2, DFT/MRCI, OM2/MRCI) to explore the S(0) and S(1) potential energy surfaces of OHBI and to locate the relevant minima, transition state, and minimum-energy conical intersection. These static calculations suggest the following decay mechanism: upon photoexcitation to the S(1) state, an ultrafast adiabatic charge-transfer induced excited-state intramolecular proton transfer (ESIPT) occurs that leads to the S(1) minimum-energy structure. Nearby, there is a S(1)/S(0) minimum-energy conical intersection that allows for an efficient nonadiabatic S(1) → S(0) internal conversion, which is followed by a fast ground-state reverse hydrogen transfer (GSHT). This mechanism is verified by semiempirical OM2/MRCI surface-hopping dynamics simulations, in which the successive ESIPT-GSTH processes are observed, but without cis-trans isomerization (which is a minor path experimentally with less than 5% yield). These gas-phase simulations of OHBI give an estimated first-order decay time of 476 fs for the S(1) state, which is larger but of the same order as the experimental values measured for OHBDI in solution: 270 fs in CH(3)CN and 230 fs in CH(2)Cl(2). The differences between the photoinduced processes of the 2- and 4-hydroxy-substituted chromophores are attributed to the presence or absence of intramolecular hydrogen bonding between the two rings.

Conical intersections in solution: formulation, algorithm, and implementation with combined quantum mechanics/molecular mechanics method

The significance of conical intersections in photophysics, photochemistry, and photodissociation of polyatomic molecules in gas phase has been demonstrated by numerous experimental and theoretical studies. Optimization of conical intersections of small- and medium-size molecules in gas phase has currently become a routine optimization process, as it has been implemented in many electronic structure packages. However, optimization of conical intersections of small- and medium-size molecules in solution or macromolecules remains inefficient, even poorly defined, due to large number of degrees of freedom and costly evaluations of gradient difference and nonadiabatic coupling vectors. In this work, based on the sequential quantum mechanics and molecular mechanics (QM/MM) and QM/MM-minimum free energy path methods, we have designed two conical intersection optimization methods for small- and medium-size molecules in solution or macromolecules. The first one is sequential QM conical intersection optimization and MM minimization for potential energy surfaces; the second one is sequential QM conical intersection optimization and MM sampling for potential of mean force surfaces, i.e., free energy surfaces. In such methods, the region where electronic structures change remarkably is placed into the QM subsystem, while the rest of the system is placed into the MM subsystem; thus, dimensionalities of gradient difference and nonadiabatic coupling vectors are decreased due to the relatively small QM subsystem. Furthermore, in comparison with the concurrent optimization scheme, sequential QM conical intersection optimization and MM minimization or sampling reduce the number of evaluations of gradient difference and nonadiabatic coupling vectors because these vectors need to be calculated only when the QM subsystem moves, independent of the MM minimization or sampling. Taken together, costly evaluations of gradient difference and nonadiabatic coupling vectors in solution or macromolecules can be reduced significantly. Test optimizations of conical intersections of cyclopropanone and acetaldehyde in aqueous solution have been carried out successfully.