Gas-phase reactivity of CH3OH toward OH at interstellar temperatures (11.7-177.5 K): experimental and theoretical study

Computational study of the post-transition state dynamics for the OH + CH3OH reaction probed by photodetachment of the CH3O-(H2O) anion

Dissociative photodetachment dynamics simulations were conducted to study the CH3O-(H2O) → CH3O + H2O + e- reaction using classical molecular dynamics (MD) and ring-polymer molecular dynamics (RPMD) techniques on two newly formulated neutral potential energy surfaces (PES1 and PES2) by different research groups. While the dissociation dynamics exhibited similarities between classical MD and RPMD, there were noticeable differences in the fluctuation of probability densities for the internal modes due to nuclear quantum effects. Upon comparison of our findings with experimental data concerning the electron binding energy distribution and photofragment relative energy, it suggests that the potential energy landscapes of PES2 are reasonably precise. The time evolution of occupied vibrational states of the H2O photofragment is presented in this study.

Effect of temperature on the gas-phase reaction of CH3CN with OH radicals: experimental (T = 11.7-177.5 K) and computational (T = 10-400 K) kinetic study

Acetonitrile (CH3CN) is present in the interstellar medium (ISM) in a variety of environments. However, at the ultracold temperatures of the ISM, radical-molecule reactions are not widely investigated because of the experimental handicap of getting organic molecules in the gas phase by conventional techniques. The CRESU (French acronym for Reaction Kinetics in a Uniform Supersonic Flow) technique solves this problem. For this reason, we present in this work the kinetic study of the gas-phase reaction of CH3CN with one of the most ubiquitous radicals, the hydroxyl (OH) radical, as a function of temperature (11.7-177.5 K). The kinetic technique employed to investigate the CH3CN + OH reaction was the pulsed laser photolysis-laser induced fluorescence. The rate coefficient for this reaction k(T) has been observed to drastically increase from 177.5 K to 107.0 K (about 2 orders of magnitude), while the increase in k(T) from 107.0 K to 11.7 K was milder (around 4 times). The temperature dependent expressions for k(T) are provided in the two distinct T-ranges, excluding the upper limit obtained for k(177.5 K): In addition, the rate coefficients estimated by the canonical competitive unified statistical (CCUS) theory show a similar behaviour to the experimental results, when evaluated within the high-pressure limit. This is consistent with the experimentally observed independence of k(T) with total gas density at selected temperatures. Astrochemical networks, such as the KIDA database or UMIST, do not include the CH3CN + OH reaction as a potential depletion process for acetonitrile in the ISM because the current studies predict very low rate coefficients at IS temperatures. According to the model (T = 10 K), the impact of the titled reaction on the abundances of CH3CN appears to be negligible in dark molecular clouds of the ISM (∼1% of the total depletion reactions included in UMIST network). With respect to the potential formation of the CH2CN radical in those environments, even in the most favourable scenario, where this radical could be formed in a 100% yield from the CH3CN + OH reaction, this route would only contribute around 2% to the current assumed formation routes by the UMIST network.

A systematic mechanistic survey on the reactions between OH radical and CH3OH on ice

A systematic mechanistic survey was performed for the CH3OH + OH reaction on ice. ONIOM(ωB97X-D/Def2-TZVP:AMOEBA09) calculations suggested a range of binding energies for the CH2OH radical (0.29-0.69 eV) and CH3OH (0.15-0.72 eV) molecule on hexagonal water ice (Ih) and amorphous solid water (ASW). Computed average binding energies of CH2OH radical (0.49 eV) and CH3OH (0.41 eV) are relatively stronger compared to the CH3O radical binding energies (0.32 eV, Sameera et al., J. Phy. Chem. A, 2021, 125, 387-393). Thus, the CH3OH molecule, CH2OH and CH3O radicals can adsorb on ice, where the binding energies follow the order CH2OH > CH3OH > CH3O. The multi-component artificial force-induced reaction (MC-AFIR) method systematically determined the reaction mechanisms for the CH3OH + OH reaction on ice, where two reaction paths, giving rise to CH2OH and CH3O radicals, were confirmed. A range of reaction barriers, employing the ωB97X-D/Def2-TZVP level of theory, was found for each reaction (0.03-0.11 eV for CH2OH radical formation, and 0.03-0.44 eV for CH3O radical formation). Based on the lowest energy reaction paths, we suspect that both reactions operate on ice. The computed data in this study evidence that the nature of the binding site or the reaction site has a significant effect on the computed binding energies or reaction barriers. Thus, the outcomes of the present study will be very useful for the computational astrochemistry community to determine reliable binding energies and reaction barriers on ice.

The atmospheric relevance of primary alcohols and imidogen reactions

Organic alcohols as very volatile compounds play a crucial role in the air quality of the atmosphere. So, the removal processes of such compounds are an important atmospheric challenge. The main goal of this research is to discover the atmospheric relevance of degradation paths of linear alcohols by imidogen with the aid of simulation by quantum mechanical (QM) methods. To this end, we combine broad mechanistic and kinetic results to get more accurate information and to have a deeper insight into the behavior of the designed reactions. Thus, the main and necessary reaction pathways are explored by well-behaved QM methods for complete elucidation of the studying gaseous reactions. Moreover, the potential energy surfaces as  a main factor are computed for easier judging of the most probable pathways in the simulated reactions. Our attempt to find the occurrence of the considered reactions in the atmospheric conditions is completed by precisely evaluating the rate constants of all elementary reactions. All of the computed bimolecular rate constants have a positive dependency on both temperature and pressure. The kinetic results show that H-abstraction from the α carbon is dominant relative to the other sites. Finally, by the results of this study, we conclude that at moderate temperatures and pressures primary alcohols can degrade with imidogen, so they can get atmospheric relevance.

Experimental and theoretical study of the low-temperature kinetics of the reaction of CN with CH2O and implications for interstellar environments

Rate coefficients for the reaction of CN with CH₂O were measured for the first time below room temperature in the range 32-103 K using a pulsed Laval nozzle apparatus together with the Pulsed Laser Photolysis-Laser-Induced Fluorescence technique. The rate coefficients exhibited a strong negative temperature dependence, reaching (4.62 ± 0.84) × 10^-11 cm³ molecule^-1 s^-1 at 32 K, and no pressure dependence was observed at 70 K. The potential energy surface (PES) of the CN + CH₂O reaction was calculated at the CCSD(T)/aug-cc-pVTZ//M06-2X/aug-cc-pVTZ level of theory, with the lowest energy channel to reaction characterized by the formation of a weakly-bound van der Waals complex, bound by 13.3 kJ mol^-1, prior to two transition states with energies of -0.62 and 3.97 kJ mol^-1, leading to the products HCN + HCO or HNC + HCO, respectively. For the formation of formyl cyanide, HCOCN, a large activation barrier of 32.9 kJ mol^-1 was calculated. Reaction rate theory calculations were performed with the MESMER (Master Equation Solver for Multi Energy well Reactions) package on this PES to calculate rate coefficients. While this ab initio description provided good agreement with the low-temperature rate coefficients, it was not capable of describing the high-temperature experimental rate coefficients from the literature. However, increasing the energies and imaginary frequencies of both transition states allowed MESMER simulations of the rate coefficients to be in good agreement with data spanning 32-769 K. The mechanism for the reaction is the formation of a weakly-bound complex followed by quantum mechanical tunnelling through the small barrier to form HCN + HCO products. MESMER calculations showed that channel generating HNC is not important. MESMER simulated the rate coefficients from 4-1000 K which were used to recommend best-fit modified Arrhenius expressions for use in astrochemical modelling. The UMIST Rate12 (UDfa) model yielded no significant changes in the abundances of HCN, HNC, and HCO for a variety of environments upon inclusion of rate coefficients reported here. The main implication from this study is that the title reaction is not a primary formation route to the interstellar molecule formyl cyanide, HCOCN, as currently implemented in the KIDA astrochemical model.

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.

Fragmentation of interstellar methanol by collisions with He˙+: an experimental and computational study

Methanol is a key species in astrochemistry as its presence and reactivity provides a primary route to the synthesis of more complex interstellar organic molecules (iCOMs) that may eventually be incorporated in newly formed planetary systems. In the interstellar medium, methanol is formed by hydrogenation of CO ices on grains, and its fate upon collisions with interstellar ions should be accounted for to correctly model iCOM abundances in objects at various stages of stellar evolution. The absolute cross sections (CSs) and branching ratios (BRs) for the collisions of He˙⁺ ions with CH₃OH are measured, as a function of the collision energy, using a Guided Ion Beam Mass Spectrometer (GIB-MS). Insights into the dissociative electron (charge) exchange mechanism have been obtained by computing the entrance and exit multidimensional Potential Energy Surfaces (PESs) and by modelling the non-adiabatic transitions using an improved Landau-Zener-Stückelberg approach. Notably, the dynamical treatment reproducing the experimental findings includes a strong orientation effect of the system formed by the small He˙⁺ ion and the highly polar CH₃OH molecule, in the electric field gradient associated to the strongly anisotropic intermolecular interaction. This is a stereodynamical effect that plays a fundamental role in collision events occurring under a variety of conditions, with kinetic energy confined within intervals ranging from the sub-thermal to the hyper-thermal regime.

Reaction of OH radicals with CH3NH2 in the gas phase: experimental (11.7-177.5 K) and computed rate coefficients (10-1000 K)

Nitrogen-bearing molecules, like methylamine (CH₃NH₂), can be the building blocks of amino acids in the interstellar medium (ISM). At the ultralow temperatures of the ISM, it is important to know its gas-phase reactivity towards interstellar radicals and the products formed. In this work, the kinetics of the OH + CH₃NH₂ reaction was experimentally and theoretically investigated at low- and high-pressure limits (LPL and HPL) between 10 and 1000 K. Moreover, the CH₂NH₂ and CH₃NH yields were computed in the same temperature range for both pressure regimes. A pulsed CRESU (French acronym for Reaction Kinetics in a Uniform Supersonic Flow) apparatus was employed to determine the rate coefficient, k(T), in the 11.7-177.5 K range. A drastic increase of k(T) when the temperature is lowered was observed in agreement with theoretical calculations, evaluated by the competitive canonical unified statistical (CCUS) theory, below 300 K in the LPL regime. The same trend was observed in the HPL regime below 350 K, but the theoretical k(T) values were higher than the experimental ones. Above 200 K, the calculated rate coefficients are improved with respect to previous computational studies and are in excellent agreement with the experimental literature data. In the LPL, the formation of CH₃NH becomes largely dominant below ca. 100 K. Conversely, in the HPL regime, CH₂NH₂ is the only product below 100 K, whereas CH₃NH becomes dominant at 298 K with a branching ratio similar to the one found in the LPL regime (≈70%). At T > 300 K, both reaction channels are competitive independently of the pressure regime.

Interstellar hide and go seek: C3H4O. There and back (again)

The molecular species C₃H₄O represents a striking example of an astrochemical conundrum. With more than 60 structural isomers theoretically possible, to date only acrolein (CH₂CHCHO) has been identified in the Sgr B2(N) region of the interstellar medium (ISM). The topography of the singlet potential energy surface is complicated, with three low-lying minima predicted to be almost isoenergetic: cis and trans-acrolein, and methylketene (CH₃CHCO). Our CCSD(T)/cc-pVTZ calculations confirm that methylketene is energetically lower than cis-acrolein, lying only 1.9 kJ mol^-1 above the trans-isomer, which is the global minimum. In this respect, methylketene is a promising candidate for interstellar observation. Unfortunately, however, despite several searches its astronomical detection has been unsuccessful. To this end, the key question is whether in fact methylketene exists as a discrete chemical entity in the ISM at all? In this paper, we present a detailed examination of the C₃H₄O potential energy surface, with specific focus on formation pathways. CCSD(T)/cc-pVTZ calculations enable a more elaborate interpretation of reaction mechanisms than was published hitherto. Our results show that gauche-propargyl alcohol and syn and anti-allenol emerge as interesting new targets for observational astronomers in TMC-1: given the recent discovery of the propargyl radical in this region, barrierless product channels involving OH˙ lend support to their candidacy as possible interstellar species. Finally, this work provides accurate spectral data of these three potential molecules, to be used for searches in interstellar space.

Low-temperature reaction dynamics of paramagnetic species in the gas phase

Radicals are abundant in a range of important gas-phase environments. They are prevalent in the atmosphere, in interstellar space, and in combustion processes. As such, understanding how radicals react is essential for the development of accurate models of the complex chemistry occurring in these gas-phase environments. By controlling the properties of the colliding reactants, we can also gain insights into how radical reactions occur on a fundamental level. Recent years have seen remarkable advances in the breadth of experimental methods successfully applied to the study of reaction dynamics involving paramagnetic species-from improvements to the well-known crossed molecular beams approach to newer techniques involving magnetically guided and decelerated beams. Coupled with ever-improving theoretical methods, quantum features are being observed and interesting insights into reaction dynamics are being uncovered in an increasingly diverse range of systems. In this highlight article, we explore some of the exciting recent developments in the study of chemical dynamics involving paramagnetic species. We focus on low-energy reactive collisions involving neutral radical species, where the reaction parameters are controlled. We conclude by identifying some of the limitations of current methods and exploring possible new directions for the field.

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.

The impact of water vapor on the OH reactivity toward CH3CHO at ultra-low temperatures (21.7-135.0 K): Experiments and theory

The role of water vapor (H2O) and its hydrogen-bonded complexes in the gas-phase reactivity of organic compounds with hydroxyl (OH) radicals has been the subject of many recent studies. Contradictory effects have been reported at temperatures between 200 and 400 K. For the OH + acetaldehyde reaction, a slight catalytic effect of H2O was previously reported at temperatures between 60 and 118 K. In this work, we used Laval nozzle expansions to reinvestigate the impact of H2O on the OH-reactivity with acetaldehyde between 21.7 and 135.0 K. The results of this comprehensive study demonstrate that water, instead, slows down the reaction by factors of ∼3 (21.7 K) and ∼2 (36.2-89.5 K), and almost no effect of added H2O was observed at 135.0 K.

Simplified Protocol for the Calculation of Multiconformer Transition State Theory Rate Constants Applied to Tropospheric OH-Initiated Oxidation Reactions

Chemical kinetics plays a fundamental role in the understanding and modeling of tropospheric chemical processes, one of the most important being the atmospheric degradation of volatile organic compounds. These potentially harmful molecules are emitted into the troposphere by natural and anthropogenic sources and are chemically removed by undergoing oxidation processes, most frequently initiated by reaction with OH radicals, the atmosphere's "detergent". Obtaining the respective rate constants is therefore of critical importance, with calculations based on transition state theory (TST) often being the preferred choice. However, for molecules with rich conformational variety, a single-conformer method such as lowest-conformer TST is unsuitable while state-of-the-art TST-based methodologies easily become unmanageable. In this Feature Article, the author reviews his own cost-effective protocol for the calculation of bimolecular rate constants of OH-initiated reactions in the high-pressure limit based on multiconformer transition state theory. The protocol, which is easily extendable to other oxidation reactions involving saturated organic molecules, is based on a variety of freeware and open-source software and tested against a series of oxidation reactions of hydrofluoropolyethers, computationally very challenging molecules with potential environmental relevance. The main features, advantages and disadvantages of the protocol are presented, along with an assessment of its predictive utility based on a comparison with experimental rate constants.

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.

A new instrument for kinetics and branching ratio studies of gas phase collisional processes at very low temperatures

A new instrument dedicated to the kinetic study of low-temperature gas phase neutral-neutral reactions, including clustering processes, is presented. It combines a supersonic flow reactor with vacuum ultra-violet synchrotron photoionization time-of-flight mass spectrometry. A photoion-photoelectron coincidence detection scheme has been adopted to optimize the particle counting efficiency. The characteristics of the instrument are detailed along with its capabilities illustrated through a few results obtained at low temperatures (<100 K) including a photoionization spectrum of n-butane, the detection of formic acid dimer formation, and the observation of diacetylene molecules formed by the reaction between the C₂H radical and C₂H₂.

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.

Gas phase reaction kinetics of complex organic molecules at temperatures of the interstellar medium: The OH + CH3OH case

Recent experimental and theoretical works concerning gas-phase radical-neutral reactions involving Complex Organic Molecules are reviewed in the context of cold interstellar objects with a special emphasis on the OH + CH3OH reaction and its potential impact on the formation of CH3O.

Gas-phase kinetics of CH3CHO with OH radicals between 11.7 and 177.5 K

Gas-phase reactions in the interstellar medium (ISM) are a source of molecules in this environment. The knowledge of the rate coefficient for neutral-neutral reactions as a function of temperature, k(T), is essential to improve astrochemical models. In this work, we have experimentally measured k(T) for the reaction between the OH radical and acetaldehyde, both present in many sources of the ISM. Laser techniques coupled to a CRESU system were used to perform the kinetic measurements. The obtained modified Arrhenius equation is k(T = 11.7-177.5 K) = (1.2 ± 0.2) × 10-11 (T/300 K)-(1.8±0.1) exp-{(28.7 ± 2.5)/T} cm3 molecule-1 s-1. The k(T) value of the title reaction has been measured for the first time below 60 K. No pressure dependence of k(T) was observed at ca. 21, 50, 64 and 106 K. Finally, a pure gas-phase model indicates that the title reaction could become the main CH3CO formation pathway in dark molecular clouds, assuming that CH3CO is the main reaction product at 10 K.

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.

Experimental and theoretical investigation on the OH + CH3C(O)CH3 reaction at interstellar temperatures (T=11.7-64.4 K)

The rate coefficient, k(T), for the gas-phase reaction between OH radicals and acetone CH3C(O)CH3, has been measured using the pulsed CRESU (French acronym for Reaction Kinetics in a Uniform Supersonic Flow) technique (T = 11.7-64.4 K). The temperature dependence of k(T = 10-300 K) has also been computed using a RRKM-Master equation analysis after partial revision of the potential energy surface. In agreement with previous studies we found that the reaction proceeds via initial formation of two pre-reactive complexes both leading to H2O + CH3C(O)CH2 by H-abstraction tunneling. The experimental k(T) was found to increase as temperature was lowered. The measured values have been found to be several orders of magnitude higher than k(300 K). This trend is reproduced by calculations, with a special good agreement with experiments below 25 K. The effect of total gas density on k(T) has been explored. Experimentally, no pressure dependence of k(20 K) and k(64 K) was observed, while k(50 K) at the largest gas density 4.47×1017 cm-3 is twice higher than the average values found at lower densities. The computed k(T) is also reported for 103 cm-3 of He (representative of the interstellar medium). The predicted rate coefficients at 10 K surround the experimental value which appears to be very close to the low pressure regime prevailing in the interstellar medium. For gas-phase model chemistry of interstellar molecular clouds, we suggest using the calculated value of 1.8×10-10 cm3 molecule-1 s-1 at 10 K and the reaction products are water and CH3C(O)CH2 radicals.

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

A recent Ring Polymer Molecular Dynamics study of the reactions of OH with methanol and formaldehyde, at zero pressure and below 100 K, has shown the formation of long lived complexes, with long lifetimes, longer than 100 ns for the lower temperatures studied, 20-100 K (del Mazo-Sevillano et al., 2019). These long lifetimes support the existence of multi collision events with the He buffer-gas atoms under experimental conditions, as suggested by several transition state theory studies of these reactions. In this work we study these secondary collisions, as a dynamical approach to study pressure effects on these reactions. For this purpose, the potential energy surfaces of He with H2CO, OH, H2O and HCO are calculated at highly accurate ab initio level. The stability of some of the complexes is studied using Path Integral Molecular dynamics techniques, determining that OH-H2CO complexes can be formed up to 100 K or higher temperatures, while the weaker He-H2CO complexes dissociate at approximately 50 K. The predicted IR intensity spectra shows new features which could help the identification of the OH-H2CO complex. Finally, the He-H2CO + OH and OH-H2CO + He collisions are studied using quassi-classical trajectories, finding that the cross section to produce HCO + H2O products increases with decreasing collision energy, and that it is ten times higher in the He-H2CO + OH case.

Quantum Roaming in the Complex-Forming Mechanism of the Reactions of OH with Formaldehyde and Methanol at Low Temperature and Zero Pressure: A Ring Polymer Molecular Dynamics Approach

The quantum dynamics of the title reactions are studied using the ring polymer molecular dynamics (RPMD) method from 20 to 1200 K using recently proposed full dimensional potential energy surfaces which include long-range dipole-dipole interactions. A V-shaped dependence of the reaction rate constants is found with a minimum at 200-300 K, in rather good agreement with the current experimental data. For temperatures above 300 K the reaction proceeds following a direct H-abstraction mechanism. However, below 100 K the reaction proceeds via organic-molecule···OH collision complexes, with very long lifetimes, longer than 10^-7 s, associated with quantum roaming arising from the inclusion of quantum effects by the use of RPMD. The long lifetimes of these complexes are comparable to the time scale of the tunnelling to form reaction products. These complexes are formed at zero pressure because of quantum effects and not only at high pressure as suggested by transition state theory (TST) calculations for OH + methanol and other OH reactions. The zero-pressure rate constants reproduce quite well measured ones below 200 K, and this agreement opens the question of how important the pressure effects on the reaction rate constants are, as implied in TST-like formalisms. The zero-pressure mechanism is applicable only to very low gas density environments, such as the interstellar medium, which are not repeatable by experiments.