Full rate constant matrix contraction method for obtaining branching ratio of unimolecular decomposition

The branching ratio of unimolecular decomposition can be evaluated by solving the rate equations. Recent advances in automated reaction path search methods have enabled efficient construction of the rate equations based on quantum chemical calculations. However, it is still difficult to solve the rate equations composed of hundreds or more elementary steps. This problem is especially serious when elementary steps that occur in highly different timescales coexist. In this article, we introduce an efficient approach to obtain the branching ratio from a given set of rate equations. It has been derived from a recently proposed rate constant matrix contraction (RCMC) method, and termed full‐RCMC (f‐RCMC). The f‐RCMC gives the branching ratio without solving the rate equations. Its performance was tested numerically for unimolecular decomposition of C3H5 and C4H5. Branching ratios obtained by the f‐RCMC precisely reproduced the values obtained by numerically solving the rate equations. It took about 95 h to solve the rate equations of C4H5 consisting of 234 elementary steps. In contrast, the f‐RCMC gave the branching ratio in less than 1 s. The f‐RCMC would thus be an efficient alternative of the conventional kinetic simulation approach. © 2016 Wiley Periodicals, Inc.

[1]  W. Yim,et al.  Application of ab initio molecular dynamics for a priori elucidation of the mechanism in unimolecular decomposition: the case of 5-nitro-2,4-dihydro-3H-1,2,4-triazol-3-one (NTO). , 2001, Journal of the American Chemical Society.

[2]  Julius Jellinek,et al.  Energy Landscapes: With Applications to Clusters, Biomolecules and Glasses , 2005 .

[3]  G. Whitten,et al.  Accurate and Facile Approximation for Vibrational Energy‐Level Sums , 1963 .

[4]  James A. Miller,et al.  From the Multiple-Well Master Equation to Phenomenological Rate Coefficients: Reactions on a C3H4 Potential Energy Surface , 2003 .

[5]  O. Isayev,et al.  Ab initio molecular dynamics study on the initial chemical events in nitramines: thermal decomposition of CL-20. , 2008, The journal of physical chemistry. B.

[6]  W. Tsang.  The stability of alkyl radicals , 1985 .

[7]  Michael A. Collins,et al.  Rigorous derivation of reaction path degeneracy in transition state theory , 1992 .

[8]  James A. Miller,et al.  Master equation methods in gas phase chemical kinetics. , 2006, The journal of physical chemistry. A.

[9]  S. Eser,et al.  Kinetics of Supercritical-Phase Thermal Decomposition of C10−C14 Normal Alkanes and Their Mixtures , 1997 .

[10]  Tetsuya Taketsugu,et al.  Artificial Force Induced Reaction (AFIR) Method for Exploring Quantum Chemical Potential Energy Surfaces. , 2016, Chemical record.

[11]  Stefan Grimme,et al.  Towards first principles calculation of electron impact mass spectra of molecules. , 2013, Angewandte Chemie.

[12]  J. Doye,et al.  The effect of the range of the potential on the structures of clusters , 1995 .

[13]  R. Kaiser,et al.  Chemistry of energetically activated cumulenes - from allene (H2CCCH2) to hexapentaene (H2CCCCCCH2). , 2008, Chemphyschem : a European journal of chemical physics and physical chemistry.

[14]  Robert J. Gelinas,et al.  Stiff systems of kinetic equations—A practitioner's view☆ , 1972 .

[15]  L. Drahos,et al.  MassKinetics: a theoretical model of mass spectra incorporating physical processes, reaction kinetics and mathematical descriptions. , 2001, Journal of mass spectrometry : JMS.

[16]  J. S. Francisco,et al.  Ab initio study of hydrogen migration across n-alkyl radicals. , 2011, The journal of physical chemistry. A.

[17]  K. Morokuma,et al.  Global investigation of potential energy surfaces for the pyrolysis of C(1)-C(3) hydrocarbons: toward the development of detailed kinetic models from first principles. , 2015, Physical chemistry chemical physics : PCCP.

[18]  P. Savage.  Mechanisms and Kinetics Models for Hydrocarbon Pyrolysis , 2000 .

[19]  R. Kaiser,et al.  Theoretical study of the C6H3 potential energy surface and rate constants and product branching ratios of the C2H(2Sigma+) + C4H2(1Sigma(g)+) and C4H(2Sigma+) + C2H2(1Sigma(g)+) reactions. , 2008, The Journal of chemical physics.

[20]  R. Walsh,et al.  Thermal unimolecular reactions of hydrocarbons , 1969 .

[21]  Johanna L. Miller,et al.  Dissociation channels of the 1-buten-2-yl radical and its photolytic precursor 2-bromo-1-butene. , 2005, The journal of physical chemistry. A.

[22]  I. R. Slagle,et al.  Unimolecular Decomposition of n-C4H9 and iso-C4H9 Radicals , 1996 .

[23]  M. Karplus,et al.  The topology of multidimensional potential energy surfaces: Theory and application to peptide structure and kinetics , 1997 .

[24]  Z. Belohlav,et al.  Experimental study of hydrocarbon structure effects on the composition of its pyrolysis products , 2010 .

[25]  A. Mebel,et al.  Ab initio molecular orbital/Rice–Ramsperger–Kassel–Marcus theory study of multichannel rate constants for the unimolecular decomposition of benzene and the H+C6H5 reaction over the ground electronic state , 2001 .

[26]  K. Morokuma,et al.  Intrinsic reaction coordinate: Calculation, bifurcation, and automated search , 2015 .

[27]  R. Kaiser,et al.  Understanding the chemical dynamics of the reactions of dicarbon with 1-butyne, 2-butyne, and 1,2-butadiene--toward the formation of resonantly stabilized free radicals. , 2014, Physical chemistry chemical physics : PCCP.

[28]  L. H. Andersen,et al.  Branching ratios for the dissociative recombination of hydrocarbon ions. I: The cases of C4H9+ and C4H5+ , 2003 .

[29]  D. Neumark,et al.  An ab initio/RRKM study of product branching ratios in the photodissociation of buta-1,2- and -1,3-dienes and but-2-yne at 193 nm. , 2003, Chemistry.

[30]  Tamiki Komatsuzaki,et al.  Kinetic Analysis for the Multistep Profiles of Organic Reactions: Significance of the Conformational Entropy on the Rate Constants of the Claisen Rearrangement. , 2015, The journal of physical chemistry. A.

[31]  W. Goddard,et al.  The Mechanism for Unimolecular Decomposition of RDX (1,3,5-Trinitro-1,3,5-triazine), an ab Initio Study , 2000 .

[32]  H. Schlegel,et al.  Ab initio classical trajectory study of H2CO+H2 + CO dissociation , 1994 .

[33]  M. Koshi,et al.  Direct investigations on the thermal unimolecular isomerization reaction of 1-pentyl radicals , 2002 .

[34]  H. Bernhard Schlegel,et al.  Ab initio molecular dynamics studies of the photodissociation of formaldehyde, H2CO→H2+CO: Direct classical trajectory calculations by MP2 and density functional theory , 2000 .

[35]  Donald G Truhlar,et al.  Modeling the kinetics of bimolecular reactions. , 2006, Chemical reviews.

[36]  K. Morokuma,et al.  Ab Initio Study of the Mechanism for the Thermal Decomposition of the Phenoxy Radical , 1996 .

[37]  Michael A. Collins,et al.  Molecular potential-energy surfaces for chemical reaction dynamics , 2002 .

[38]  Tamás Turányi,et al.  On the error of the quasi-steady-state approximation , 1993 .

[39]  Yu Zhuang,et al.  Direct chemical dynamics simulations: coupling of classical and quasiclassical trajectories with electronic structure theory , 2013 .