How Does the Central Atom Substitution Impact the Properties of a Criegee Intermediate? Insights from Multireference Calculations

Semiclassical Dynamics on Machine-Learned Coupled Multireference Potential Energy Surfaces: Application to the Photodissociation of the Simplest Criegee Intermediate

Determination of high-dimensional potential energy surfaces (PESs) and nonadiabatic couplings have always been quite challenging. To this end, machine learning (ML) models, trained with a finite set of ab initio data, allow accurate prediction of such properties. To express the PESs in terms of atomic contributions is the cornerstone of any ML based technique because it can be easily scaled to large systems. In this work, we have constructed high fidelity PESs and nonadiabatic coupling terms at the CASSCF level of ab initio data using a machine learning technique, namely, kernel-ridge regression. Additional MRCI-level calculations were carried out to assess the quality of the PESs. We use these machine-learned PESs and nonadiabatic couplings to simulate excited-state molecular dynamics based on Tully's fewest-switches surface hopping method (FSSH). FSSH is a semiclassical method in which nuclei move on the PESs due to the electrons according to the laws of classical mechanics. Nonadiabatic effects are taken into account in terms of transitions between PESs. We apply this scheme to study the O-O photodissociation of the simplest Criegee intermediate (CH₂OO). The FSSH trajectories were initiated on the lowest optically bright singlet excited state (S₂) and propagated along the three most important internal coordinates, namely, O-O and C-O bond distances and the COO bond angle. Some of the trajectories end up on energetically lower PESs as a result of radiationless transfer through conical intersections. All of the trajectories lead to the dissociation of the O-O bond due to the dissociative nature of the excited PESs through one of the two dissociative channels. The simulation reveals that there is about 88.4% probability of dissociation through the lower channel leading to the H₂CO (X¹A₁) and O (¹D) products, whereas there is only 11.6% probability of dissociation through the upper channel leading to H₂CO (a³A″) and O (³P) products.

Spectroscopy and photochemistry of ClSSO

Sulfur-chlorine cycles play a role in the atmosphere of Venus. It is thought that many sulfur-chlorine bearing molecules could be present in Venus's atmosphere and play an important role in its chemical processes. The goal of this work is to provide new insight into the electronic structure and spectroscopy of the [Cl, S, S, O] molecular system. Eight isomers could be formed, but only three were found to be thermodynamically stable relative to the first dissociation limit. We spectroscopically characterized the two lowest energy stable isomers, C1-ClSSO and trans-ClSSO, using the accurate CCSD(T)-F12/aug-cc-pVTZ method. The dipole moments of the two lowest energy stable isomers are predicted to be 1.90 and 1.60 debye, respectively. The C1-ClSSO isomer is suitable for laser induced fluorescence detection since the lowest excited electronic states absorb in the visible, ∼610 nm, and near UV region, 330 nm. We mapped the evolution of the low-lying excited electronic states along the ClS, SS, and SO bond lengths to find that the production of ClS, SO, or S₂O is plausible, whereas the production of ClS₂ is not allowed.

Hydrogen bond, ring tension and π-conjugation effects: methyl and vinyl substitutions dramatically change the photodynamics of Criegee intermediates

The substitution effect in chemistry is a concept that is probably too common to mention, while for a molecule with an elusive electronic structure, substitution can introduce an unusual effect that dramatically tunes the chemical process. To reveal the substitution effects on the photodynamics of Criegee Intermediates (CIs), we carried out the multireference CASSCF trajectory surface-hopping (TSH) molecular dynamics and CASPT2 electronic-structure calculations for a methyl-substituted CI (MCI) and a vinyl-substituted CI (VCI). The results show that for different substituents, the hydrogen bond, ring tension and π-conjugation not only alter the relative stabilities of the conformers/configurations, but also dramatically change the photo-induced channel of CIs. For an anti-MCI, the dominant channel starting from the S₁ state is the ring-closure process leading to dioxirane, while in the syn configuration, the intramolecular CHO hydrogen bond hinders the rotation around the C-O bond and thus leads to a high yield of in-plane O-O dissociation towards acetaldehyde (X¹A') and the O(¹D) atom. In a VCI with an unsaturated substituent, the π-conjugation greatly strengthens the O-O bond and therefore no O-O dissociation is observed in all configurations. In addition, the CHO hydrogen bond in the syn_(CO)-VCI further stabilizes the S₁-state intermediates and makes them less reactive; in contrast, isomerization to dioxirane becomes the globally dominant channel in the anti_(CO)-VCI. The dramatic substitution effects by saturated and unsaturated substituents on CIs found here will deepen the understanding of Criegee-intermediate chemistry.

Photochemistry from low-lying states of HOSO. J Chem Phys 2020;152:134302. [PMID: 32268736 DOI: 10.1063/5.0001867]

Using configuration interaction ab initio methods, the evolution of the lowest electronic states of singlet and triplet spin multiplicities of HOSO⁺ along the stretching and bending coordinates of is investigated. Equilibrium geometries, rotational constants, and harmonic vibrational frequencies of the lowest electronic states are calculated, i.e., X¹A', 1¹A″, 1³A', and 1³A″. The global minimum of the 1¹A″ state is located below the first dissociation limit and its calculated lifetime is predicted to be 0.40 µs, making it suitable for detection by laser-induced fluorescence. According to the potential energy surfaces, HOSO⁺ should produce SO₂ ⁺ and H after ultraviolet photon absorption to the 2¹A' state. This work opens the door to investigate the branching ratio and the production rates of SO₂ ⁺, SO⁺, and OH from HOSO⁺. These insights can help understand the SO₂ cycle in the earth's atmosphere and its effect on cooling our planet.

Gauging stability and reactivity of carbonyl O-oxide Criegee intermediates

In this study, we evaluated the effect of substitution on the stability and reactivity of carbonyl O-oxide Criegee intermediates (CIs). In this regard, we computed a set of more than 50 carbonyl oxides at the CBS-QB3 level of theory and assessed their stability by means of an isodesmic reaction equation defining a carbonyl oxide stabilization energy (COSE). Almost all substituents are stabilizing and amino groups in particular leading to COSE values of almost 60 kcal mol^-1. As opposed to π-donors, substituents with a strong σ-electron pull destabilize the C[double bond, length as m-dash]O-O group. Furthermore, we studied how the intrinsic stabilization of the Criegee intermediate is reflected in its C[double bond, length as m-dash]O and O-O bond lengths as well as the partial charges on the individual atoms of the carbonyl oxide moiety. As a potential measure for reactivity, we determined the adiabatic singlet-triplet energy gap of all carbonyl oxides. Amino substituted CIs exhibit high-lying triplet states and have relatively large barriers towards addition of water or the OH radical. However, the ΔE_S-T cannot serve as a rigorous measure for carbonyl oxide reactivity.

An Extended Computational Study of Criegee Intermediate-Alcohol Reactions

High-level ab initio calculations (DF-LCCSD(T)-F12a//B3LYP/aug-cc-pVTZ) are performed on a range of stabilized Criegee intermediate (sCI)-alcohol reactions, computing reaction coordinate energies, leading to the formation of α-alkoxyalkyl hydroperoxides (AAAHs). These potential energy surfaces are used to model bimolecular reaction kinetics over a range of temperatures. The calculations performed in this work reproduce the complicated temperature-dependent reaction rates of CH₂OO and (CH₃)₂COO with methanol, which have previously been experimentally determined. This methodology is then extended to compute reaction rates of 22 different Criegee intermediates with methanol, including several intermediates derived from isoprene ozonolysis. In some cases, sCI-alcohol reaction rates approach those of sCI-(H₂O)₂. This suggests that in regions with elevated alcohol concentrations, such as urban Brazil, these reactions may generate significant quantities of AAAHs and may begin to compete with sCI reactions with other trace tropospheric pollutants such as SO₂. This work also demonstrates the ability of alcohols to catalyze the 1,4-H transfer unimolecular decomposition of α-methyl substituted sCIs.

Photochemistry of the Simplest Criegee Intermediate, CH2OO: Photoisomerization Channel toward Dioxirane Revealed by CASPT2 Calculations and Trajectory Surface-Hopping Dynamics

The photochemistry of Criegee intermediates plays a significant role in atmospheric chemistry, but it is relatively less explored compared with their thermal reactions. Using multireference CASPT2 electronic structure calculations and CASSCF trajectory surface-hopping molecular dynamics, we have revealed a dark-state-involved A¹A → X¹A photoisomerization channel of the simple Criegee intermediate (CH₂OO) that leads to a cyclic dioxirane. The excited molecules on the A¹A state, which can have either originated from the B¹A state via B¹A → A¹A internal conversion or formed by state-selective electronic excitation, is driven by the out-of-plane motion toward a perpendicular A/X¹A minimal-energy crossing point (MECI) then radiationless decay to the ground state with an average time constant of ∼138 fs, finally forming dioxirane at ∼254 fs. The dynamics starting from the A¹A state show that the quantum yield of photoisomerization from the simple Criegee intermediate to dioxirane is 38%. The finding of the A¹A → X¹A photoisomerization channel is expected to broaden the reactivity profile and deepen the understanding of the photochemistry of Criegee intermediates.