Roaming radicals in the thermal decomposition of dimethyl ether: Experiment and theory

Authors

R. Sivaramakrishnan, Argonne National Laboratory, Argonne, IL
J. Michael, Argonne National Laboratory, Argonne, IL
A. Wagner, Argonne National Laboratory, Argonne, IL
R. Dawes, Missouri University of Science and Technology, Rolla, MO
A. Jasper, Sandia National Laboratories, Livermore, CA
L. Harding, Argonne National Laboratory, Argonne, IL
Y. Georgievskii, Argonne National Laboratory, Argonne, IL
S. Klippenstein, Argonne National Laboratory, Argonne, IL

Date of this Version

2011

Citation

Combustion and Flame 158 (2011) 618–632;

doi:10.1016/j.combustflame.2010.12.017

Abstract

The thermal dissociation of dimethyl ether has been studied with a combination of reflected shock tube experiments and ab initio dynamics simulations coupled with transition state theory based master equation calculations. The experiments use the extraordinary sensitivity provided by H-atom ARAS detection with an unreversed light source to measure both the total decomposition rate and the branching to radical products versus molecular products, with the molecular products arising predominantly through roaming according to the theoretical analysis. The experimental observations also provide a measure of the rate coefficient for H + CH3OCH3. An evaluation of the available experimental results for H + CH₃OCH₃ can be expressed by a three parameter Arrhenius expression as,

k = 6.54 x 10^-24T^4.13 exp(-896/T)cm³ molecule^-1s^-1(273-1465 K)

The potential energy surface is explored with high level ab initio electronic structure theory. The dynamics of roaming versus radical formation is studied with a reduced dimensional trajectory approach. The requisite potential energy surface is obtained from an interpolative moving least squares fit to wide-ranging ab initio data for the long-range interactions between methyl and methoxy. The predicted roaming and radical micro-canonical fluxes are incorporated in a master equation treatment of the temperature and pressure dependence of the dissociation process. The tight (i.e., non-roaming) transition states leading to a variety of additional molecular fragments are also included in the master equation analysis, but are predicted to have a negligible contribution to product formation. The final theoretical results reliably reproduce the measured dissociation rate to radical products reported here and are well reproduced over the 500–2000 K temperature range and the 0.01–300 bar pressure range by the following modified Arrhenius parameters for the Troe falloff format:

k_1,∞(T)= 2.33 x 10¹⁹T^-0.661exp(-42345/T)s^-1

k_1,0(T) = 2.86 x 10³⁵T^-11.4exp(-46953/T)cm³ molecule^-1s^-1

F_cent(T)= exp(-T/880)

The experimentally observed branching ratio of 0.19 ± 0.07 provides a direct measure of the contribution from the roaming radical mechanism. The theoretical analysis predicts a much smaller roaming contribution of 0.02.