Ozonolysis of 2-Methyl-2-pentenal: New Insights from Master Equation Modeling

Experimental and theoretical studies were carried out to investigate the ozonolysis of trans-2-methyl-2-pentenal. The experiments were conducted in atmospheric simulation chambers coupled to a Fourier transform infrared (FTIR) spectrometer and a gas chromatograph-mass spectrometer at room temperature and atmospheric pressure in the presence of an excess of cyclohexane in dry conditions (RH < 1%). The ozonolysis reaction was investigated theoretically from the results of accurate density functional (M06-2X) and ab initio [CCSD(T)] computations, employing the AVTZ basis set. The sequence of reaction steps was established, and the system of kinetics equations was modeled using MESMER. In the first step, a primary ozonide is formed, which then decomposes along two pathways. The principal ozonolysis products are propanal, methylglyoxal, ethylformate, and a secondary ozonide. An interesting competition between sequential reaction steps and well-skipping is found, which leads to an inversion of the expected methylglyoxal/propanal product ratio at temperatures below 210 K. The mechanism of the "hot ester" reaction channel of the Criegee intermediate was revisited. The computed ozonolysis rate constant and product branching ratio are in excellent agreement with the experimental data that are also reported in the present work.

Ab initio study of the O3-N2 complex: Potential energy surface and rovibrational states

The formation and destruction of O₃ within the Chapman cycle occurs as a result of inelastic collisions with a third body. Since N₂ is the most abundant atmospheric molecule, it can be considered as the most typical candidate when modeling energy-transfer dynamics. We report a new ab initio potential energy surface (PES) of the O₃-N₂ van der Waals complex. The interaction energies were calculated using the explicitly correlated single- and double-excitation coupled cluster method with a perturbative treatment of triple excitations [CCSD(T)-F12a] with the augmented correlation-consistent triple-zeta aug-cc-pVTZ basis set. The five-dimensional PES was analytically represented by an expansion in spherical harmonics up to eighth order inclusive. Along with the global minimum of the complex (D_e = 348.88 cm^-1), with N₂ being perpendicular to the O₃ plane, six stable configurations were found with a smaller binding energy. This PES was employed to calculate the bound states of the O₃-N₂ complex with both ortho- and para-N₂ for total angular momentum J = 0 and 1, as well as dipole transition probabilities. The nature of the bound states of the O₃-oN₂ and O₃-pN₂ species is discussed based on their rovibrational wave functions.

Inelastic scattering in isotopologues of O2-Ar: the effects of mass, symmetry, and density of states

The two species considered here, O2 (oxygen molecule) and Ar (argon-atom), are both abundant components of Earth's atmosphere and hence familiar collision partners in this medium. O2 is quite reactive and extensively involved in atmospheric chemistry, including Chapman's cycle of the formation and destruction of ozone; while Ar, like N2, typically plays the nevertheless crucial role of inert collider. Inert species can provide stabilization to metastable encounter-complexes through the energy transfer associated with inelastic collisions. The interplay of collision frequency and energy transfer efficiency, with state lifetimes and species concentrations, contributes to the rich and varied chemistry and dynamics found in diverse environments ranging from planetary atmospheres to the interstellar and circumstellar media. The nature and density of bound and resonance states, coupled electronic states, symmetry, and nuclear spin-statistics can all play a role. Here, we systematically investigate some of those factors by looking at the O2-Ar system, comparing rigorous quantum-scattering calculations for the 16O16O-40Ar, 18O16O-40Ar, and 18O18O-40Ar isotope combinations. A new accurate potential energy surface was constructed for this purpose holding the O2 bond distance at its vibrationally averaged distance.

Spectroscopy and Scattering Studies Using Interpolated Ab Initio Potentials

The Born-Oppenheimer potential energy surface (PES) has come a long way since its introduction in the 1920s, both conceptually and in predictive power for practical applications. Nevertheless, nearly 100 years later-despite astonishing advances in computational power-the state-of-the-art first-principles prediction of observables related to spectroscopy and scattering dynamics is surprisingly limited. For example, the water dimer, (H₂O)₂, with only six nuclei and 20 electrons, still presents a formidable challenge for full-dimensional variational calculations of bound states and is considered out of reach for rigorous scattering calculations. The extremely poor scaling of the most rigorous quantum methods is fundamental; however, recent progress in development of approximate methodologies has opened the door to fairly routine high-quality predictions, unthinkable 20 years ago. In this review, in relation to the workflow of spectroscopy and/or scattering studies, we summarize progress and challenges in the component areas of electronic structure calculations, PES fitting, and quantum dynamical calculations.

Computational study of the rovibrational spectrum of CO2-N2

The CO2-N2 complex is formed from two key components of Earth's atmosphere, and as such, has received some attention from both experimental and theoretical studies. On the theory side, a potential energy surface (PES) based on high level ab initio data was reported [Nasri et al., J. Chem. Phys., 2015, 142, 174301] and then used in more recently reported rovibrational calculations [Lara-Moreno et al., Phys. Chem. Chem. Phys., 2019, 21, 3550]. Accuracy of about 1 percent was achieved for calculated rotational transitions of the ground vibrational state of the complex, compared with previously reported microwave spectra. However, a very recent measurement of the geared bending mode frequency [Barclay et al., J. Chem. Phys., 2020, 153, 014303] recorded a value of 21.4 cm-1, which is wildly different from the corresponding calculated value of 45.9 cm-1. To provide some insight into this discrepancy, we have constructed a new more accurate PES, and used it to perform highly converged variational rovibrational calculations. Our new results yield a value of 21.1 cm-1 for that bending frequency, in close agreement with the experiment. We also obtain significantly improved predicted rotational transitions. Finally, we note that a very shallow well, previously reported as a distinct second isomer, is not found on our new PES, but rather a transition structure is seen in that location.

Rotational quenching of an interstellar gas thermometer: CH3CNHe collisions

Among all the molecular species found in the interstellar medium, molecules with threefold symmetry axes play a special role, as their rotational spectroscopy allows them to act as practical gas thermometers. Methyl-cyanide (CH3CN) is the second most abundant of those (after ammonia). We compute in this paper the collisional dynamics of methyl-cyanide in collision with helium, for both the A- and the E-symmetries of CH3CN. The potential energy surface is determined using the CCSD(T)-F12b formalism and fit with convenient analytic functions. We compute the rotationally inelastic cross sections for all levels up to 510 cm-1 of collision energy, employing at low energy exact Coupled Channels methods, and at higher energies, approximate Coupled States methods. For temperatures from 7 K up to 300 K, rates of quenching are computed and most are found to differ from those reported earlier (up to a factor of a thousand), calling for a possible reexamination of the temperatures assigned to low density gasses.

Benchmark ab initio and dynamical characterization of the stationary points of reactive atom + alkane and SN2 potential energy surfaces

We review composite ab initio and dynamical methods and their applications to characterize stationary points of atom/ion + molecule reactions.

Rotationally inelastic scattering of O3–Ar: state-to-state rates with the multiconfigurational time dependent Hartree method

The rates of state-changing collisions are compared for different isotopologues of ozone from quantum scattering calculations with the MCTDH method.