Thermal Activation of Methane by MgO${}^{+}$: Temperature Dependent Kinetics, Reactive Molecular Dynamics Simulations and Statistical Modeling

Abstract

The kinetics of MgO${}^{+}$ + CH${}_4$ was studied experimentally using the variable ion source, temperature adjustable selected ion flow tube (VISTA-SIFT) apparatus from 300 $-$ 600 K and computationally by running and analyzing reactive atomistic simulations. Rates and product branching fractions were determined as a function of temperature. The reaction proceeded with a rate of $k=5.9\pm 1.510^{-10}(T/300$ K${)}^{-0.5\pm 0.2}$ cm${}^3$ s${}^{-1}$. MgOH${}^{+}$ was the dominant product at all temperatures, but Mg${}^{+}$, the co-product of oxygen-atom transfer to form methanol, was observed with a product branching fraction of $0.08\pm 0.03(T/300$ K${)}^{-0.8\pm 0.7}$. Reactive molecular dynamics simulations using a reactive force field, as well as a neural network yield rate coefficients about one order of magnitude lower. This underestimation of the rates is traced back to the multireference character of the transition state [MgOCH${}_4$]${}^{+}$. Statistical modeling of the temperature-dependent kinetics provides further insight into the reactive potential surface. The rate limiting step was found to be consistent with a four-centered activation of the C-H bond, consistent with previous calculations. The product branching was modeled as a competition between dissociation of an insertion intermediate directly after the rate-limiting transition state, and traversing a transition state corresponding to a methyl migration leading to a Mg-CH${}_3$OH${}^{+}$ complex, though only if this transition state is stabilized significantly relative to the dissociated MgOH${}^{+}$ + CH${}_3$ product channel. An alternative non-statistical mechanism is discussed, whereby a post-transition state bifurcation in the potential surface could allow the reaction to proceed directly from the four-centered TS to the Mg-CH${}_3$OH${}^{+}$ complex thereby allowing a more robust competition between the product channels.

References in corpus (4)

• PhysNet: A Neural Network for Predicting Energies, Forces, Dipole Moments and Partial Charges
• Benchmark {\em ab initio} energy profiles for the gas-phase S${}_N$2 reactions Y${}^{-}$ + CH${}_3$X $\to$ CH${}_3$Y + X${}^{-}$ (X,Y = F,Cl,Br). Validation of hybrid DFT methods
• The $\mathrm{N}{(}^{4}S)+{\mathrm{O}}_2(X^3{\Sigma }_g^{-})\leftrightarrow \mathrm{O}{(}^{3}P)+\mathrm{NO}(X^2\Pi )$ Reaction: Thermal and Vibrational Relaxation Rates for the ${}^2$A${}^{\prime }$, ${}^4$A${}^{\prime }$ and ${}^2$A${}^{\prime \prime }$ States
• Reactive Atomistic Simulations of Diels-Alder Reactions: the Importance of Molecular Rotations

Cited by in corpus (4)

• Machine Learning Force Fields
• Conformational and state-specific effects in reactions of 2,3-dibromobutadiene with Coulomb-crystallized calcium ions
• Transfer Learned Potential Energy Surfaces: Accurate Anharmonic Vibrational Dynamics and Dissociation Energies for the Formic Acid Monomer and Dimer
• Reactive atomistic simulations of Diels-Alder-type reactions: Conformational and dynamic effects in the polar cycloaddition of 2,3-dibromobutadiene radical ions with maleic anhydride