Quantum yields of singlet and triplet chemiexcitation of dimethyl 1,2-dioxetane: ab initio nonadiabatic molecular dynamic simulations

Ab initio molecular dynamics study of intersystem crossing dynamics for MH2 (M = Si, Ge, Sn, Pb) on spin-pure and spin-mixed potential energy surfaces

Recently, surface-hopping ab initio molecular dynamics (SH-AIMD) simulations have come to be used to discuss the mechanisms and dynamics of excited-state chemical reactions, including internal conversion and intersystem crossing. In dynamics simulations involving intersystem crossing, there are two potential energy surfaces (PESs) governing the motion of nuclei: PES in a spin-pure state and PES in a spin-mixed state. The former gives wrong results for molecular systems with large spin-orbit coupling (SOC), while the latter requires a potential gradient that includes a change in SOC at each point, making the computational cost very high. In this study, we systematically investigate the extent to which the magnitude of SOC affects the results of the spin-pure state-based dynamics simulations for the hydride MH2 (M = Si, Ge, Sn, Pb) by performing SH-AIMD simulations based on spin-pure and spin-mixed states. It is clearly shown that spin-mixed state PESs are indispensable for the dynamics simulation of intersystem crossing in systems containing elements Sn and Pb from the fifth period onward. Furthermore, in addition to the widely used Tully's fewest switches (TFS) algorithm, the Zhu-Nakamura (ZN) global switching algorithm, which is computationally less expensive, is applied to SH for comparison. The results from TFS- and ZN-SH-AIMD methods are in qualitative agreement, suggesting that the less expensive ZN-SH-AIMD can be successfully utilized to investigate the dynamics of photochemical reactions based on quantum chemical calculations.

Dissociation Time, Quantum Yield, and Dynamic Reaction Pathways in the Thermolysis of trans-3,4-Dimethyl-1,2-dioxetane

The thermolysis of trans-3,4-dimethyl-1,2-dioxetane is studied by trajectory surface hopping. The significant difference between long and short dissociation times is rationalized by frustrated dissociations and the time spent in triplet states. If the C-C bond breaks through an excited state channel, then the trajectory passes over a ridge of the potential energy surface of that state. The calculated triplet quantum yields match the experimental results. The dissociation half-times and quantum yields follow the same ascending order as per the product states, justifying the conjecture that the longer dissociation time leads to a higher quantum yield, proposed in the context of the methylation effect. The populations of the molecular Coulomb Hamiltonian and diagonal states reach equilibrium, but the triplet populations with different Sz components fluctuate indefinitely. Certain initial velocities, leading the trajectories to given product states, can be identified as the most characteristic features for sorting trajectories according to their product states.

Chemiexcitation: Mammalian Photochemistry in the Dark†

Light is one way to excite an electron in biology. Another is chemiexcitation, birthing a reaction product in an electronically excited state rather than exciting from the ground state. Chemiexcited molecules, as in bioluminescence, can release more energy than ATP. Excited states also allow bond rearrangements forbidden in ground states. Molecules with low-lying unoccupied orbitals, abundant in biology, are particularly susceptible. In mammals, chemiexcitation was discovered to transfer energy from excited melanin, neurotransmitters, or hormones to DNA, creating the lethal and carcinogenic cyclobutane pyrimidine dimer. That process was initiated by nitric oxide and superoxide, radicals triggered by ultraviolet light or inflammation. Several poorly understood chronic diseases share two properties: inflammation generates those radicals across the tissue, and cells that die are those containing melanin or neuromelanin. Chemiexcitation may therefore be a pathogenic event in noise- and drug-induced deafness, Parkinson's disease, and Alzheimer's; it may prevent macular degeneration early in life but turn pathogenic later. Beneficial evolutionary selection for excitable biomolecules may thus have conferred an Achilles heel. This review of recent findings on chemiexcitation in mammalian cells also describes the underlying physics, biochemistry, and potential pathogenesis, with the goal of making this interdisciplinary phenomenon accessible to researchers within each field.

Conical Intersection in Chemiluminescence of Cyclic Peroxides

Chemiluminescence (CL) utilizing chemiexcitation for energy transformation is one of the most highly sensitive and useful analytical techniques. The chemiexcitation is a chemical process of a ground-state reactant producing an excited-state product, in which a nonadiabatic event is facilitated by conical intersections (CIs), the specific molecular geometries where electronic states are degenerated. Cyclic peroxides, especially 1,2-dioxetane/dioxetanone derivatives, are the iconic chemiluminescent substances. In this Perspective, we concentrated on the CIs in the CL of cyclic peroxides. We first present a computational overview on the role of CIs between the ground (S₀) state and the lowest singlet excited (S₁) state in the thermolysis of cyclic peroxides. Subsequently, we discuss the role of the S₀/S₁ CI in the CL efficiency and point out misunderstandings in some theoretical studies on the singlet chemiexcitations of cyclic peroxides. Finally, we address the challenges and future prospects in theoretically calculating S₀/S₁ CIs and simulating the dynamics and chemiexcitation efficiency in the CL of cyclic peroxides.

A Look Inside the Black Box of Machine Learning Photodynamics Simulations

ConspectusPhotochemical reactions are of great importance in chemistry, biology, and materials science because they take advantage of a renewable energy source, mild reaction conditions, and high atom economy. Light absorption can excite molecules to a higher energy electronic state of the same spin multiplicity. The following nonadiabatic processes induce molecular transformations that afford exotic molecular architectures and high-energy-isomers that are inaccessible by thermal means. Computational simulations now complement time-resolved instrumentation to reveal ultrafast excited-state mechanistic information for photochemical reactions that is essential in disentangling elusive spectroscopic features, excited-state lifetimes, and excited-state mechanistic critical points. Nonadiabatic molecular dynamics (NAMD), powered by surface hopping techniques, is among the most widely applied techniques to model the photochemical reactions of medium-sized molecules. However, the computational efficiency is limited because of the requisite thousands of multiconfigurational quantum-chemical calculations multiplied by hundreds of trajectories. Machine learning (ML) has emerged as a revolutionary force in computational chemistry to predict the outcome of the resource-intensive multiconfigurational calculations on the fly. An ML potential trained with a substantial set of quantum-chemical calculations can predict the energies and forces with errors under chemical accuracy at a negligible cost. The integration of ML potentials in NAMD dramatically extends the maximum simulation time scale by ∼10 000-fold to the nanosecond regime.In this Account, we present a comprehensive demonstration of ML photodynamics simulations and summarize our most recent applications in resolving complex photochemical reactions. First, we address three fundamental components of ML techniques for photodynamics simulations: the quantum-chemical data set, the ML potential, and NAMD. Second, we describe best practices in building training data and our procedure toward training the ML photodynamics model with our recent literature contributions. We introduce a convenient training data generation scheme combining Wigner sampling and geometrical interpolation. It trains reliable and effective ML potentials suitable for subsequent active learning to detect undersampled data. We demonstrate how active learning automatically discovers new mechanistic pathways and reproduces experimental results. We point out that atomic permutation is an essential data augmentation approach to improve the learnability of distance-based molecular descriptors for highly symmetric molecules. Third, we demonstrate the utility of ML-photodynamics by showing the results of ML photodynamics simulations of (1) photo-torquoselective 4π disrotatory electrocyclic ring closing of norbornyl cyclohexadiene, which reveals a thermal conversion from experimentally unobserved intermediates to the reactant in 1 ns; (2) [2 + 2] photocycloaddition of substituted [3]-syn-ladderdienes in competition with 4π and 6π electrocyclic ring-opening reactions, uncovering substituent effects to explain the reported increased quantum yield of substituted cubane precursors; and (3) photochemical 4π disrotatory electrocyclic reactions of fluorobenzenes in nanoseconds with XMS-CASPT2-level training data. We expect this Account to broaden understanding of ML photodynamics and inspire future developments and applications to increasingly large molecules within complex environments on long time scales.

Trajectory surface hopping molecular dynamics simulations for retinal protonated Schiff-base photoisomerization

Global switching trajectory surface hopping molecular dynamics simulations are performed using accurate on-the-fly (TD)CAM-B3LYP/6-31G potential energy surfaces to study retinal protonated Schiff-base photoisomerization up to S₁ excitation. The simulations detected two-layer conical intersection networks: one is at an energy as high as 8 eV and the other is in the energy range around 3-4 eV. Six conical intersections within the low-layer energy region that correspond to active conical intersections under experimental conditions are found via the use of pairwise isomers, within which nonadiabatic molecular dynamics simulations are performed. Eight isomer products are populated with simulated sampling trajectories from which the simulated quantum yield in the gas phase is estimated to be 0.11 (0.08) moving from the all-trans isomer to the 11-cis (11-cis to all-trans) isomer in comparison with an experimental value of 0.09 (0.2) in the solution phase. Each conical intersection is related to one specific twist angle accompanying a related CC double bond motion during photoisomerization. Nonplanar distortion of the entire dynamic process has a significant role in the formation of the relevant photoisomerization products. The present simulation indicates that all hopping points show well-behaved potential energy surface topology, as calculated via the conventional TDDFT method, at conical intersections between S₁ and S₀ states. Therefore, the present nonadiabatic dynamics simulations with the TDDFT method are very encouraging for simulating various large systems related to retinal Schiff-base photoisomerization in the real world.

H2O-Bridged Proton-Transfer Channel in Emitter Species Formation in Obelin Bioluminescence

Bioluminescence of a number of marine organisms is conditioned by Ca²⁺-regulated photoprotein (CaRP) with coelenterazine as the reaction substrate. The reaction product, coelenteramide, at the first singlet excited state (S₁) is the emitter of CaRP. The S₁-state coelenteramide is produced via the decomposition of coelenterazine dioxetanone. Experiments suggested that the neutral S₁-coelenteramide is the primary emitter species. This supposition contradicts with theoretical calculations showing that the anionic S₁-coelenteramide is a primary product of the decomposition of coelenterazine dioxetanone. In this study, applying molecular dynamic (MD) simulations and the hybrid quantum mechanics/molecular mechanics (QM/MM) method, we investigated a proton-transfer (PT) process taking place in CaRP obelin from Obelia longissima for emitter formation. Our calculations demonstrate a concerted PT process with a water molecule as a bridge between anionic S₁-coelenteramide and the nearest histidine residue. The low activation barrier as well as the strong hydrogen-bond network between the proton donor and the proton acceptor suggests a fast PT process comparable with that of the lifetime of excited anionic S₁-coelenteramide. The existence of the PT process eliminates the discrepancy between experimental and theoretical studies. The fast PT process at emitter formation can also take place in other CaRPs.

Photoisomerization-mechanism-associated excited-state hydrogen transfer in 2'-hydroxychalcone revealed by on-the-fly trajectory surface-hopping molecular dynamics simulation

By performing global-switching on-the-fly trajectory surface-hopping molecular dynamics simulation at the OM2/MRCI (14,15) quantum level, we probed the S3(ππ*) photoisomerization mechanisms associated with excited-state intramolecular hydrogen transfer for 2'-hydroxychalcone (2HC) within the interwoven conical intersection networks from four singlet electronic states (S3, S2, S1, and S0). The simulated quantum yields of 0.03 for cis-to-trans and zero for trans-to-cis photoisomerization were due to almost all the conical intersections being localized either in the cis-2HC or in trans-2HC region, and there was little chance for sampling trajectories to reach the rotation conical intersection (S1/S0) in between cis-2HC and trans-2HC that is key for reactive isomerization. The potential energy well on the S1 state in the trans-2HC region prevents trajectories from trans-to-cis photoisomerization, while the fact there is no well on S1 state in cis-2HC region opens a few chances for trajectories to reach the rotation conical intersections. The present simulation found that excited-state intramolecular hydrogen transfers in 2HC have a negative impact for reactive isomerization, and that hydrogen transfers take place on the S1 state, while back-transfer on the S0 state prevents sampling trajectories reaching rotational conical intersections. It was realized that it could be possible to enhance the quantum yield of 2HC photoisomerization by suppressing the hydrogen transfer (such as by changing an electron-donating substitution or adjusting the substitution position to decrease the acidity of the phenol group). From a perspective view of the potential energy surfaces, the theoretical design of such 2HC derivatives could enhance (control) the quantum yield by shifting the conical intersections away from the cis- and trans-region.