Molecular Simulation for Atmospheric Reaction Exploration and Discovery: Non-Equilibrium Dynamics, Roaming and Glycolaldehyde Formation Following Photo-Induced Decomposition of syn-Acetaldehyde Oxide

Meenu Upadhyay,Kai Töpfer,Markus Meuwly

2023-07-06

Abstract:The decomposition and chemical dynamics for vibrationally excited syn-CH$_3$CHOO is followed based on statistically significant numbers of molecular dynamics simulations. Using a neural network-based reactive potential energy surface, transfer learned to the CASPT2 level of theory, the final total kinetic energy release and rotational state distributions of the OH fragment are in quantitative agreement with experiment. In particular the widths of these distributions are sensitive to the experimentally unknown strength of the O--O bond strength, for which values $D_e \in [22,25]$ kcal/mol are found. Due to the non-equilibrium nature of the process considered, the energy-dependent rates do not depend appreciably on the O--O scission energy. Roaming dynamics of the OH-photoproduct leads to formation of glycolaldehyde on the picosecond time scale with subsequent decomposition into CH$_2$OH+HCO. Atomistic simulations with global reactive machine-learned energy functions provide a viable route to quantitatively explore the chemistry and reaction dynamics for atmospheric reactions.

Chemical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to explore and discover non - equilibrium dynamics, roaming mechanisms and the formation of glyceraldehyde in atmospheric reactions through molecular dynamics simulations. Specifically, the study focused on the decomposition and chemical kinetic processes of vibrationally excited syn - CH3CHOO. By using a neural - network - based reactive potential energy surface and combining transfer learning at the CASPT2 theoretical level, the authors aimed to:

1. **Quantify experimental results**: Verify the consistency between the final total kinetic energy release (TKER) of the OH fragment and the rotational state distribution obtained from the simulation and the experimental results, especially that the widths of these distributions are very sensitive to the O - O bond strength which is unknown in the experiment. The value of the O - O bond strength found in the study is between 22 and 25 kcal/mol.

2. **Explore new reaction paths**: Explore the formation of glyceraldehyde on the picosecond time scale and its subsequent decomposition into CH2OH + HCO by simulating the roaming dynamics of the OH photoproduct. This reveals that there are other possible reaction paths in addition to the known CH2CHO + OH product.

3. **Provide new insights into atmospheric chemical reactions**: Use a machine - learning energy function of global reactivity for atomic - level simulations, providing a feasible way to quantitatively explore the chemistry and reaction kinetics of atmospheric reactions.

Through detailed molecular dynamics simulations, the paper not only verifies experimental observations but also discovers new reaction mechanisms, such as the roaming behavior of the OH fragment leading to the formation of glyceraldehyde, thus providing important theoretical support for understanding the complex chemical processes in the atmosphere.