Zero- and high-pressure mechanisms in the complex forming reactions of OH with methanol and formaldehyde at low temperatures

Analytical potential energy surface and dynamics for the OH + CH3OH reaction

Using as functional form a combination of valence bond and mechanic molecular terms a new full-dimensional potential energy surface was developed for the title reaction, named PES-2022, which was fitted to high-level ab initio calculations at the coupled-cluster singles, doubles, and perturbative triples-F12 explicitly correlated level on a representative number of points describing the reactive system. This surface simultaneously describes the two reaction channels, hydrogen abstraction from the methyl group [(R1) path] and from the alcohol group [(R2) path] of methanol to form water. PES-2022 is a smooth and continuous surface, which reasonably describes the topology of this reactive system from reactants to products, including the intermediate complexes present in the system. Based on PES-2022 an exhaustive dynamics study was performed using quasi-classical trajectory calculations under two different initial conditions: at a fixed room temperature, for direct comparison with the experimental evidence and at different collision energies, to analyze possible mechanisms of reaction. In the first case, the available energy was mostly deposited as water vibrational energy, with the vibrational population inverted in the stretching modes and not inverted in the bending modes, reproducing the experimental evidence. In the second case, the analysis of different dynamics magnitudes (excitation functions, product energy partitioning, and product scattering distributions), allows us to suggest different mechanisms for both (R1) and (R2) paths: a direct mechanism for the (R2) path vs an indirect one, related with "nearly trapped" trajectories in the intermediate complexes, for the (R1) path.

A comprehensive benchmark ab initio survey of the stationary points and products of the OH· + CH3OH system

Reactions between methanol and the hydroxyl radical are of significant interest for combustion-, atmospheric-, and astrochemistry. While the two primary product channels (the formation of H₂O with either CH₃O· or ·CH₂OH) have been the subject of numerous studies, the possibility of other products has seen little attention. Here, we present a comprehensive thermochemical survey of the stationary points and plausible products of the reaction, featuring 29 geometries optimized at the UCCSD(T)-F12b/aug-cc-pVTZ level, followed by accurate composite ab initio computations for all stationary points (including ·CH₂OH dissociation and isomerization) and five product channels, with a detailed evaluation of basis set convergence and efficiency. The computations reveal that the formation of methanediol and the hydroxymethoxy radical is thermodynamically favorable and the endothermicity of formaldehyde formation is low enough to be a plausible product channel. We also observe unexpectedly large energy deviations between the partially-spin-adapted ROHF-RCCSD(T) method and ROHF-UCCSD(T) as well as between UHF-UCCSDT(Q) and ROHF-UCCSDT(Q) results.

Probing the Exit Channel of the OH + CH3OH → H2O + CH3O Reaction by Photodetachment of CH3O-(H2O)

Transition state dynamics of bimolecular reactions can be probed by photodetachment of a precursor anion when the Franck-Condon region of the corresponding neutral potential energy surface is near a saddle point. In this study, photodetachment of anions at m/z = 49 enabled investigation of the exit channel of the OH + CH₃OH → H₂O + CH₃O reaction using photoelectron-photofragment coincidence spectroscopy. High-level coupled-cluster calculations of the stationary points on the anion surface show that the methoxide-water cluster CH₃O^-(H₂O) is the stable minimum on the anion surface. Photodetachment at a 3.20 eV photon energy leads to long-lived H₂O(CH₃O) complexes and H₂O + CH₃O products consistent with both direct dissociative photodetachment and resonance mediated processes on the neutral surface. The partitioning of total kinetic energy in the system indicates that water stretch and bend excitation is induced in dissociative photodetachment and evidence for long-lived complexes consistent with vibrational Feshbach resonances is reported.

Neural network potential energy surface for the low temperature ring polymer molecular dynamics of the H2CO + OH reaction

A new potential energy surface (PES) and dynamical study of the reactive process of H₂CO + OH toward the formation of HCO + H₂O and HCOOH + H are presented. In this work, a source of spurious long range interactions in symmetry adapted neural network (NN) schemes is identified, which prevents their direct application for low temperature dynamical studies. For this reason, a partition of the PES into a diabatic matrix plus a NN many-body term has been used, fitted with a novel artificial neural network scheme that prevents spurious asymptotic interactions. Quasi-classical trajectory (QCT) and ring polymer molecular dynamics (RPMD) studies have been carried on this PES to evaluate the rate constant temperature dependence for the different reactive processes, showing good agreement with the available experimental data. Of special interest is the analysis of the previously identified trapping mechanism in the RPMD study, which can be attributed to spurious resonances associated with excitations of the normal modes of the ring polymer.

Microcanonical Tunneling Rates from Density-of-States Instanton Theory

Semiclassical instanton theory is a form of quantum transition-state theory which can be applied to the computation of thermal reaction rates in complex molecular systems including quantum tunneling effects. There have been a number of attempts to extend the theory to treat microcanonical rates. However, the previous formulations are either computationally unfeasible for large systems due to an explicit sum over states or they involve extra approximations, which make them less reliable. We propose a robust and practical microcanonical formulation called density-of-states instanton theory, which avoids the sum over states altogether. In line with the semiclassical approximations inherent to the instanton approach, we employ the stationary-phase approximation to the inverse Laplace transform to obtain the densities of states. This can be evaluated using only post-processing of the data available from a small set of instanton calculations, such that our approach remains computationally efficient. We show that the new formulation predicts results that agree well with quantum scattering theory for an atom-diatom reaction and with experiments for a photoexcited unimolecular hydrogen transfer in a Criegee intermediate. When the thermal rate is evaluated from a Boltzmann average over our new microcanonical formalism, it can overcome some problems of conventional instanton theory. In particular, it predicts a smooth transition at the crossover temperature and is able to describe bimolecular reactions with pre-reactive complexes such as CH₃OH + OH.

Quantum Effects on the D + H3+ → H2D+ + H Deuteration Reaction and Isotopic Variants

The title reaction and its isotopic variants are studied using quasi-classical trajectory (QCT) (without taking into account corrections to account for the possible zero point energy breakdown) and ring polymer molecular dynamics (RPMD) methods with a full dimensional and accurate potential energy surface which presents an exchange barrier of approximately 0.144 eV. The QCT rate constant increases when the temperature decreases from 1500 to 10 K. On the contrary, the RPMD rate constant decreases with decreasing temperature, in semiquantitative agreement with recent experimental results. The present RPMD results are in between the thermal and translational experimental rate constants, extracted from the measured data to eliminate the initial vibrational excitation of H₃⁺, obtained in an arc discharge. The difference between the present RPMD results and experimental values is attributed to the possible existence of non thermal vibrational excitation of H₃⁺, not completely removed by the semiempirical model used for the analysis of the experimental results. Also, it is found that, below 200 K, the RPMD trajectories are trapped, forming long-lived collision complexes, with lifetimes longer than 1 ns. These collision complexes can fragment by either redissociating back to reactants or react to products, in the two cases tunneling through the centrifugal and reaction barriers, respectively. The contribution of the formation of the complex to the total deuteration rate should be calculated with more accurate quantum methods, as has been found recently for reactions of larger systems, and the present four atoms system is a good candidate to benchmark the adequacy of RPMD method at temperatures below 100 K.