Global uncertainty analysis for RRKM/master equation based kinetic predictions: A case study of ethanol decomposition

Abstract A precise understanding of the accuracy of reaction rate constants, whether determined experimentally or theoretically, is of considerable importance to kinetic modelers. While the uncertainties of experimentally measured rate constants are commonly provided, the “error bars” of computed (temperature- and pressure-dependent) rate constants are rarely evaluated rigorously. In this work, global uncertainty and sensitivity analysis is applied to the propagation of the uncertainties in the input parameters (e.g. barrier heights, frequencies and collisional energy transfer parameters et al.) to those in the rate constants computed by the RRKM/master equation method for the decomposition of ethanol. This case study provides a systematic exploration of the effect of temperature and pressure on the parametric uncertainties in RRKM/master equation calculations for a prototypical single-well multiple-channel dissociation. In the high pressure limit, the uncertainties in the theoretical predictions are controlled by the uncertainties in the input parameters involved in the transition state theory calculations, with the most important ones being those describing the energetics of the decomposition. At lower pressures, where fall-off is important, the uncertainties in the collisional energy transfer parameters play a significant role, particularly for the higher energy of the two channels. Remarkably, the competition between dissociation and collisional excitation leads to uncertainties of more than a factor of 100 in the predictions for the higher energy channel. These large uncertainties are related to the need for large-scale single-collision-induced transitions in energy in order to produce the higher energy products in the low pressure limit. The present study illustrates the value of detailed qualitative and quantitative studies of the uncertainties in theoretical kinetics predictions.

[1]  B. C. Garrett,et al.  Variational Transition State Theory , 1980 .

[2]  D. Golden,et al.  Master equation analysis of pressure-dependent atmospheric reactions. , 2003, Chemical reviews.

[3]  Yuyang Li,et al.  An experimental and kinetic modeling investigation on a rich premixed n-propylbenzene flame at low pressure , 2013 .

[4]  Frank M. Mourits,et al.  A critical evaluation of Lennard–Jones and Stockmayer potential parameters and of some correlation methods , 1977 .

[5]  C. Eckart The Penetration of a Potential Barrier by Electrons , 1930 .

[6]  Emilio Martínez-Núñez,et al.  Dinámica de reacciones unimoleculares en fase gas: Desviaciones del comportamiento estadístico , 2002 .

[7]  M. Lin,et al.  Thermal decomposition of ethanol. I. Ab Initio molecular orbital/Rice–Ramsperger–Kassel–Marcus prediction of rate constant and product branching ratios , 2002 .

[8]  Tiziano Faravelli,et al.  An experimental and kinetic modeling study of combustion of isomers of butanol , 2010 .

[9]  H. Rabitz,et al.  Practical Approaches To Construct RS-HDMR Component Functions , 2002 .

[10]  James A. Miller,et al.  Predictive a priori pressure-dependent kinetics , 2014, Science.

[11]  J. Barker,et al.  Master equation simulations of competing unimolecular and bimolecular reactions: application to OH production in the reaction of acetyl radical with O2. , 2007, Physical chemistry chemical physics : PCCP.

[12]  Upakarasamy Lourderaj,et al.  Theoretical and computational studies of non-RRKM unimolecular dynamics. , 2009, The journal of physical chemistry. A.

[13]  Alison S. Tomlin,et al.  GUI-HDMR - A software tool for global sensitivity analysis of complex models , 2009, Environ. Model. Softw..

[14]  S. Klippenstein,et al.  Rate constants for the thermal decomposition of ethanol and its bimolecular reactions with OH and D: reflected shock tube and theoretical studies. , 2010, The journal of physical chemistry. A.

[15]  James A. Miller,et al.  Lennard–Jones parameters for combustion and chemical kinetics modeling from full-dimensional intermolecular potentials , 2014 .

[16]  Kun Xu,et al.  Thermal decomposition of ethanol. 4. Ab initio chemical kinetics for reactions of H atoms with CH3CH2O and CH3CHOH radicals. , 2011, Journal of Physical Chemistry A.

[17]  H. Matsui,et al.  Shock tube study on the thermal decomposition of ethanol. , 2010, The journal of physical chemistry. A.

[18]  S. Klippenstein,et al.  Predictive theory for the combination kinetics of two alkyl radicals. , 2006, Physical chemistry chemical physics : PCCP.

[20]  Herschel Rabitz,et al.  Global uncertainty assessments by high dimensional model representations (HDMR) , 2002 .

[21]  M. Pilling,et al.  MESMER: an open-source master equation solver for multi-energy well reactions. , 2012, The journal of physical chemistry. A.

[22]  James A. Miller,et al.  Theoretical unimolecular kinetics for CH4 + M ⇄ CH3 + H + M in eight baths, M = He, Ne, Ar, Kr, H2, N2, CO, and CH4. , 2011, The journal of physical chemistry. A.

[23]  D. Thompson Practical methods for calculating rates of unimolecular reactions , 1998 .

[24]  M. Olzmann,et al.  Pyrolysis of ethanol: A shock-tube/TOF-MS and modeling study , 2015 .

[25]  J. Bozzelli,et al.  Variational analysis of the phenyl + O2 and phenoxy + O reactions. , 2008, The journal of physical chemistry. A.

[26]  S. Klippenstein Variational optimizations in the Rice–Ramsperger–Kassel–Marcus theory calculations for unimolecular dissociations with no reverse barrier , 1992 .

[27]  Alison S. Tomlin,et al.  Global sensitivity analysis of a 3D street canyon model—Part I: The development of high dimensional model representations , 2008 .

[28]  H. Rabitz,et al.  Efficient input-output model representations , 1999 .

[29]  R. Reid,et al.  The Properties of Gases and Liquids , 1977 .

[30]  A. Mebel,et al.  An ab initio/RRKM study of the reaction mechanism and product branching ratios of the reactions of ethynyl radical with 1,2-butadiene , 2011 .

[31]  Herschel Rabitz,et al.  Efficient Implementation of High Dimensional Model Representations , 2001 .

[32]  R. Ghanem,et al.  Quantifying uncertainty in chemical systems modeling , 2004 .

[33]  W. Green,et al.  Reaction of phenyl radical with propylene as a possible source of indene and other polycyclic aromatic hydrocarbons: an ab initio/RRKM-ME study. , 2012, The journal of physical chemistry. A.

[34]  J. Troe Theory of thermal unimolecular reactions at low pressures. III. Superposition of weak and strong collisions , 1992 .

[35]  D. Dixon,et al.  Theoretical prediction of the heats of formation of C2H5O* radicals derived from ethanol and of the kinetics of beta-C-C scission in the ethoxy radical. , 2007, The journal of physical chemistry. A.

[36]  S. Nordholm,et al.  On the role of nonergodicity and slow IVR in unimolecular reaction rate theory—A review and a view , 2001 .

[37]  A. Jasper,et al.  Hydrogen-assisted isomerizations of fulvene to benzene and of larger cyclic aromatic hydrocarbons. , 2013 .

[38]  Jenghan Wang,et al.  Mechanism of Ethanol Reforming: Theoretical Foundations , 2009 .

[39]  James A. Miller,et al.  "Third-Body" Collision Efficiencies for Combustion Modeling: Hydrocarbons in Atomic and Diatomic Baths , 2015 .

[40]  H. Rabitz,et al.  High Dimensional Model Representations , 2001 .

[41]  H. Najm,et al.  Uncertainty quantification in reacting-flow simulations through non-intrusive spectral projection , 2003 .

[42]  Kun Xu,et al.  Kinetics for the reactions of phenyl with methanol and ethanol: Comparison of theory and experiment , 2013 .

[43]  T. Ziehn,et al.  A global sensitivity study of sulfur chemistry in a premixed methane flame model using HDMR , 2008 .

[44]  Sean C. Smith,et al.  Theory of Unimolecular and Recombination Reactions , 1990 .

[45]  Stephen J. Klippenstein,et al.  Uncertainty propagation in the derivation of phenomenological rate coefficients from theory: A case study of n-propyl radical oxidation , 2013 .

[46]  M. Pilling,et al.  Master equation models for chemical reactions of importance in combustion. , 2003, Annual review of physical chemistry.

[47]  M. Olzmann On the role of bimolecular reactions in chemical activation systems , 2002 .

[48]  H. Rabitz,et al.  General foundations of high‐dimensional model representations , 1999 .

[49]  Herschel Rabitz,et al.  Ratio control variate method for efficiently determining high‐dimensional model representations , 2006, J. Comput. Chem..

[50]  Jan M. L. Martin,et al.  Basis set convergence in second-row compounds. The importance of core polarization functions , 1998 .

[51]  Panos G Georgopoulos,et al.  Correlation method for variance reduction of Monte Carlo integration in RS‐HDMR , 2003, J. Comput. Chem..

[52]  Cosmin Safta,et al.  Uncertainty quantification of reaction mechanisms accounting for correlations introduced by rate rules and fitted Arrhenius parameters , 2013 .

[53]  Robert E. Wilson,et al.  Fundamentals of momentum, heat, and mass transfer , 1969 .

[54]  C. Togbé,et al.  Numerical and experimental study of ethanol combustion and oxidation in laminar premixed flames and in jet-stirred reactor , 2011 .