Secondary kinetics of methanol decomposition: theoretical rate coefficients for 3CH2 + OH, 3CH2 + 3CH2, and 3CH2 + CH3

Direct variable reaction coordinate transition state theory (VRC-TST) rate coefficients are reported for the (3)CH(2) + OH, (3)CH(2) + (3)CH(2), and (3)CH(2) + CH(3) barrierless association reactions. The predicted rate coefficient for the (3)CH(2) + OH reaction (approximately 1.2 x 10(-10) cm(3) molecule(-1) s(-1) for 300-2500 K) is 4-5 times larger than previous estimates, indicating that this reaction may be an important sink for OH in many combustion systems. The predicted rate coefficients for the (3)CH(2) + CH(3) and (3)CH(2) + (3)CH(2) reactions are found to be in good agreement with the range of available experimental measurements. Product branching in the self-reaction of methylene is discussed, and the C(2)H(2) + 2H and C(2)H(2) + H2 products are predicted in a ratio of 4:1. The effect of the present set of rate coefficients on modeling the secondary kinetics of methanol decomposition is briefly considered. Finally, the present set of rate coefficients, along with previous VRC-TST determinations of the rate coefficients for the self-reactions of CH(3) and OH and for the CH(3) + OH reaction, are used to test the geometric mean rule for the CH(3), (3)CH(2), and OH fragments. The geometric mean rule is found to predict the cross-combination rate coefficients for the (3)CH(2) + OH and (3)CH(2) + CH(3) reactions to better than 20%, with a larger (up to 50%) error for the CH(3) + OH reaction.     

Thermochemistry and accurate quantum reaction rate calculations for H2/HD/D2 + CH3

Accurate quantum-mechanical results for thermodynamic data, cumulative reaction probabilities (for J = 0), thermal rate constants, and kinetic isotope effects for the three isotopic reactions H2 + CH3 --> CH4 + H, HD + CH3 --> CH4 + D, and D2 + CH3 --> CH(3)D + D are presented. The calculations are performed using flux correlation functions and the multiconfigurational time-dependent Hartree (MCTDH) method to propagate wave packets employing a Shephard interpolated potential energy surface based on high-level ab initio calculations. The calculated exothermicity for the H2 + CH3 --> CH4 + H reaction agrees to within 0.2 kcal/mol with experimentally deduced values. For the H2 + CH3 --> CH4 + H and D2 + CH3 --> CH(3)D + D reactions, experimental rate constants from several groups are available. In comparing to these, we typically find agreement to within a factor of 2 or better. The kinetic isotope effect for the rate of the H2 + CH3 --> CH4 + H reaction compared to those for the HD + CH3 --> CH4 + D and D2 + CH3 --> CH(3)D + D reactions agree with experimental results to within 25% for all data points. Transition state theory is found to predict the kinetic isotope effect accurately when the mass of the transferred atom is unchanged. On the other hand, if the mass of the transferred atom differs between the isotopic reactions, transition state theory fails in the low-temperature regime (T < 400 K), due to the neglect of the tunneling effect.     

Mechanistic and kinetic study of the O + CH2OH reaction

The mechanism for the O + CH2OH reaction was investigated by various ab initio quantum chemistry methods. For the chemical activation mechanism, that is, the addition/elimination path, the couple-cluster methods including CCSD and CCSD(T) were employed with the cc-pVXZ (X = D, T, Q, 5) basis sets. For the abstraction channels, multireference methods including CASSCF, CASPT2, and MRCISD were used with the cc-pVDZ and cc-pVTZ basis sets. It has been shown that the production of H + HCOOH is the major channel in the chemical activation mechanism. The minor channels include HCO + H2O and OH + CH2O. The hydrogen abstraction by an O atom from the CH2OH radical produces either OH + CH2O or OH + HCOH. Moreover, the two abstraction reactions are essentially barrierless processes. The rate constants for the association of O with CH2OH have been calculated using the flexible transition state theory. A weak negative temperature dependence of the rate constants is found in the range 250-1000 K. Furthermore, it is estimated that the abstraction processes also play an important role in the O + CH2OH reaction. Additionally, the falloff behavior for the OCH2OH --> H + HCOOH reaction has been investigated. The present theoretical results are compared to the experimental measurements to understand the mechanism and kinetic behavior of the O + CH2OH reaction and the unimolecular reaction of the OCH2OH radical.     

Shock-tube study of the thermal decomposition of CH3CHO and CH3CHO + H reaction

The thermal decomposition of acetaldehyde, CH3CHO + M --> CH3 + HCO + M (eq 1), and the reaction CH3CHO + H --> products (eq 6) have been studied behind reflected shock waves with argon as the bath gas and using H-atom resonance absorption spectrometry as the detection technique. To suppress consecutive bimolecular reactions, the initial concentrations were kept low (approximately 10(13) cm(-3)). Reaction was investigated at temperatures ranging from 1250 to 1650 K at pressures between 1 and 5 bar. The rate coefficients were determined from the initial slope of the hydrogen profile via k1 = [CH3CHO]0(-1) x d[H]/dt, and the temperature dependences observed can be expressed by the following Arrhenius equations: k1(T, 1.4 bar) = 2.9 x 10(14) exp(-38 120 K/T) s(-1), k1(T, 2.9 bar) = 2.8 x 10(14) exp(-37 170 K/T) s(-1), and k1(T, 4.5 bar) = 1.1 x 10(14) exp(-35 150 K/T) s(-1). Reaction was studied with C2H5I as the H-atom precursor under pseudo-first-order conditions with respect to CH3CHO in the temperature range 1040-1240 K at a pressure of 1.4 bar. For the temperature dependence of the rate coefficient the following Arrhenius equation was obtained: k6(T) = 2.6 x 10(-10) exp(-3470 K/T) cm(3) s(-1). Combining our results with low-temperature data published by other authors, we recommend the following expression for the temperature range 300-2000 K: k6(T) = 6.6 x 10(-18) (T/K) (2.15) exp(-800 K/T) cm(3) s(-1). The uncertainties of the rate coefficients k1 and k6 were estimated to be +/-30%.     

Reflected shock tube studies of high-temperature rate constants for OH + CH4 --> CH3 + H2O and CH3 + NO2 --> CH3O + NO

The reflected shock tube technique with multipass absorption spectrometric detection of OH radicals at 308 nm has been used to study the reactions OH + CH(4) --> CH(3) + H(2)O and CH(3) + NO(2) --> CH(3)O + NO. Over the temperature range 840-2025 K, the rate constants for the first reaction can be represented by the Arrhenius expression k = (9.52 +/- 1.62) x 10(-11) exp[(-4134 +/- 222 K)/T] cm(3) molecule(-1) s(-1). Since this reaction is important in both combustion and atmospheric chemistry, there have been many prior investigations with a variety of techniques. The present results extend the temperature range by 500 K and have been combined with the most accurate earlier studies to derive an evaluation over the extended temperature range 195-2025 K. A three-parameter expression describes the rate behavior over this temperature range, k = (1.66 x 10(-18))T(2.182) exp[(-1231 K)/T] cm(3) molecule(-1) s(-1). Previous theoretical studies are discussed, and the present evaluation is compared to earlier theoretical estimates. Since CH(3) radicals are a product of the reaction and could cause secondary perturbations in rate constant determinations, the second reaction was studied by OH radical production from the fast reactions CH(3)O --> CH(2)O + H and H + NO(2) --> OH + NO. The measured rate constant is 2.26 x 10(-11) cm(3) molecule(-1) s(-1) and is not dependent on temperature from 233 to 1700 K within experimental error.     

An ab initio based global potential energy surface describing CH5+ --> CH3+ + H2

A full-dimensional, ab initio based potential energy surface (PES) for CH(5)(+), which can describe dissociation is reported. The PES is a precise fit to 36173 coupled-cluster [CCSD(T)] calculations of electronic energies done using an aug-cc-pVTZ basis. The fit uses a polynomial basis that is invariant with respect to permutation of the five H atoms, and thus describes all 120 equivalent minima. The rms fitting error is 78.1 cm(-1) for the entire data set of energies up to 30,000 cm(-1) and a normal-mode analysis of CH(5)(+) also verifies the accuracy of the fit. Two saddle points have been located on the surface as well and compared with previous theoretical work. The PES dissociates correctly to the fragments CH(3)(+) + H(2) and the equilibrium geometry and normal-mode analyses of these fragments are also presented. Diffusion Monte Carlo calculations are done for the zero-point energies of CH(5)(+) (and some isotopologs) as well as for the separated fragments of CH(5)(+), CH(3)(+) + H(2) and those of CH(4)D(+), CH(3)(+) + HD and CH(2)D(+) + H(2). Values of D(0) are reported for these dissociations. A molecular dynamics calculation of CH(4)D(+) dissociation at one total energy is also performed to both validate the applicability of the PES for dynamics studies as well as to test a simple classical statistical prediction of the branching ratio of the dissociation products.     

Kinetics and product branching ratios of the reaction of (1)CH2 with H2 and D2

The reactions of singlet methylene (a(1)A1 (1)CH2) with hydrogen and deuterium have been studied by experimental and theoretical techniques. The rate coefficients for the removal of singlet methylene with H2 (k1) and D2 (k2) have been measured from 195 to 798 K and are essentially temperature-independent with values of k1 = (10.48 +/- 0.32) x 10(-11) cm(3) molecule(-1) s(-1) and k2 = (5.98 +/- 0.34) x 10(-11) cm(3) molecule(-1) s(-1), where the errors represent 2sigma, giving a ratio of k1/k2 = 1.75 +/- 0.11. In the reaction with H2, singlet methylene can be removed by reaction giving CH3 + H or deactivated to ground-state triplet methylene. Direct measurement of the H atom product showed that the fraction of relaxation decreased from 0.3 at 195 K to essentially zero at 398 K. For the reaction with deuterium, either H or D may be eliminated. Experimentally, the H:D ratio was determined to be 1.8 +/- 0.5 over the range 195-398 K. Theoretically, the reaction kinetics has been predicted with variable reaction coordinate transition state theory and with rigid-body trajectory simulations employing various high-level, ab initio-determined potential energy surfaces. The magnitudes of the calculated rate coefficients are in agreement with experiment, but the calculations show a significant negative temperature dependence that is not observed in the experimental results. The calculated and experimental H to D ratios from the reaction of singlet methylene with D2 are in good agreement, suggesting that the reaction proceeds entirely through the formation of a long-lived methane intermediate with a statistical distribution of energy.     

High-temperature rate constants for CH3OH + Kr --> products, OH + CH3OH --> products, OH + (CH3)(2)CO --> CH2COCH3 + H2O, and OH + CH3 --> CH) + H2O

The reflected shock tube technique with multipass absorption spectrometric detection of OH radicals at 308 nm (corresponding to a total path length of approximately 4.9 m) has been used to study the dissociation of methanol between 1591 and 2865 K. Rate constants for two product channels [CH3OH + Kr --> CH3 + OH + Kr (1) and CH3OH + Kr --> 1CH2 + H2O + Kr (2)] were determined. During the course of the study, it was necessary to determine several other rate constants that contributed to the profile fits. These include OH + CH3OH --> products, OH + (CH3)2CO --> CH2COCH3 + H2O, and OH + CH3 --> 1,3CH2 + H2O. The derived expressions, in units of cm(3) molecule(-1) s(-1), are k(1) = 9.33 x 10(-9) exp(-30857 K/T) for 1591-2287 K, k(2) = 3.27 x 10(-10) exp(-25946 K/T) for 1734-2287 K, kOH+CH3OH = 2.96 x 10-16T1.4434 exp(-57 K/T) for 210-1710 K, k(OH+(CH3)(2)CO) = (7.3 +/- 0.7) x 10(-12) for 1178-1299 K and k(OH+CH3) = (1.3 +/- 0.2) x 10(-11) for 1000-1200 K. With these values along with other well-established rate constants, a mechanism was used to obtain profile fits that agreed with experiment to within <+/-10%. The values obtained for reactions 1 and 2 are compared with earlier determinations and also with new theoretical calculations that are presented in the preceding article in this issue. These new calculations are in good agreement with the present data for both (1) and (2) and also for OH + CH3 --> products.     

Thermal decomposition of ethanol. III. A computational study of the kinetics and mechanism for the CH3+C2H5OH reaction

The mechanism for the CH3+C2H5OH reaction has been investigated by the modified Gaussian-2 method based on the geometric parameters of the stationary points optimized at the B3LYP/6-311+G(d,p) level of theory. Five transition states have been identified for the production of CH4+CH3CHOH (TS1), CH4+CH3CH2O (TS2), CH4+CH2CH2OH (TS3), CH3OH+CH3CH2 (TS4), and CH3CH2OCH3+H (TS5) with the corresponding barriers 12.0, 13.2, 16.0, 44.7, and 49.9 kcal/mol, respectively. The predicted rate constants and branching ratios for the three lower-energy H-abstraction reactions were calculated using the conventional and variational transition state theory with quantum-mechanical tunneling corrections for the temperature range 300-3000 K. The predicted total rate constant, kt=8.36 x 10(-76) T(20.00) exp(5258/T) cm3 mol(-1) s(-1) (300-600 K) and 6.10 x 10(-25) T(4.10)exp(-4058/T) cm3 mol(-1) s(-1) (600-3000 K), agrees closely with existing experimental data in the temperature range 403-523 K. Similarly, the predicted rate constants for CH3+CH3CD2OH and CD3+C2H5OD are also in reasonable agreement with available low temperature kinetic data.     

Ab initio potential energy surface of CH and reaction dynamics of H + CH+

This work presents a new ground state potential energy surface (PES) for CH. The potential is tested using quasi classical trajectory (QCT) and quantum reactive scattering methods for the H + CH(+) reaction. Cross sections and rate coefficients for all reaction channels up to 300 K are calculated. The abstraction rate coefficients follow the expected slightly decreasing behaviour above 90 K, but have a positive gradient with lower temperatures. The inelastic collision and exchange reaction rate constants are increasing monotonically with temperature. The rate coefficients of the exchange reaction differ significantly between QCT and quantum reactive scattering, due to intrinsic shortcomings of the QCT final state distributions.     

Theoretical study of the complex-forming CH + H2 --> CH2 + H reaction

The complex-forming CH + H2 --> CH2 + H reaction is studied employing a recently developed global potential energy function. The reaction probability in the total angular momentum J = 0 limit is estimated with a four-atom quantum wave packet method and compared with classical trajectory and statistical theory results. The formation of complexes from different reactant internal states is also determined with wave packet calculations. While there is no barrier to reaction along the minimum energy path, we find that there are angular constraints to complex formation. Trajectory-based estimates of the low-pressure rate constants are made and compared with experimental results. We find that zero-point energy violation in the trajectories is a particularly severe problem for this reaction.     

Rate constants and hydrogen isotope substitution effects in the CH3 + HCl and CH3 + Cl2 reactions

The kinetics of the CH3 + Cl2 (k2a) and CD3 + Cl2 (k2b) reactions were studied over the temperature range 188-500 K using laser photolysis-photoionization mass spectrometry. The rate constants of these reactions are independent of the bath gas pressure within the experimental range, 0.6-5.1 Torr (He). The rate constants were fitted by the modified Arrhenius expression, k2a = 1.7 x 10(-13)(T/300 K)(2.52)exp(5520 J mol(-1)/RT) and k2b = 2.9 x 10(-13)(T/300 K)(1.84)exp(4770 J mol(-1)/RT) cm(3) molecule(-1) s(-1). The results for reaction 2a are in good agreement with the previous determinations performed at and above ambient temperature. Rate constants of the CH3 + Cl2 and CD3 + Cl2 reactions obtained in this work exhibit minima at about 270-300 K. The rate constants have positive temperature dependences above the minima, and negative below. Deuterium substitution increases the rate constant, in particular at low temperatures, where the effect reaches ca. 45% at 188 K. These observations are quantitatively rationalized in terms of stationary points on a potential energy surface based on QCISD/6-311G(d,p) geometries and frequencies, combined with CCSD(T) energies extrapolated to the complete basis set limit. 1D tunneling as well as the possibility of the negative energies of the transition state are incorporated into a transition state theory analysis, an approach which also accounts for prior experiments on the CH3 + HCl system and its various deuterated isotopic substitutions [Eskola, A. J.; Seetula, J. A.; Timonen, R. S. Chem. Phys. 2006, 331, 26].     

The heats of formation of diazene, hydrazine, N2H3+, N2H5+, N2H, and N2H3 and the Methyl Derivatives CH3NNH, CH3NNCH3, and CH3HNNHCH3

The heats of formation of N(2)H, diazene (cis- and trans-N(2)H(2)), N(2)H(3), and hydrazine (N(2)H(4)), as well as their protonated species (diazenium, N(2)H(3)(+), and hydrazinium, N(2)H(5)(+)), have been calculated by using high level electronic structure theory. Energies were calculated by using coupled cluster theory with a perturbative treatment of the triple excitations (CCSD(T)) and employing augmented correlation consistent basis sets (aug-cc-pVnZ) up to quintuple-zeta, to perform a complete basis set extrapolation for the energy. Geometries were optimized at the CCSD(T) level with the aug-cc-pVDZ and aug-cc-pVTZ basis sets. Core-valence and scalar relativistic corrections were included, as well as scaled zero point energies. We find the following heats of formation (kcal/mol) at 0 (298) K: DeltaH(f)(N(2)H) = 60.8 (60.1); DeltaH(f)(cis-N(2)H(2)) = 54.9 (53.2); DeltaH(f)(trans-N(2)H(2)) = 49.9 (48.1) versus >/=48.8 +/- 0.5 (exptl, 0 K); DeltaH(f)(N(2)H(4)) = 26.6 (23.1) versus 22.8 +/- 0.2 (exptl, 298 K); DeltaH(f)(N(2)H(3)) = 56.2 (53.6); DeltaH(f)(N(2)H(3)(+)) = 231.6 (228.9); and DeltaH(f)(N(2)H(5)(+)) = 187.1 (182.7). In addition, we calculated the heats of formation of CH(3)NH(2), CH(3)NNH, and CH(3)HNNHCH(3) by using isodesmic reactions and at the G3(MP2) level. The calculated results for the hydrogenation reaction RNNR + H(2) --> RHNNHR show that substitution of an organic substituent for H improved the energetics, suggesting that these types of compounds may be possible to use in a chemical hydrogen storage system.     

Rate coefficients and equilibrium constant for the CH2CHO + O2 reaction system

The kinetics of the CH2CHO + O2 reaction was experimentally studied in two quasi-static reactors and a discharge flow-reactor at temperatures ranging from 298 to 660 K and pressures between 1 mbar and 46 bar with helium as the bath gas. The CH2CHO radicals were produced by the laser-flash photolysis of ethyl vinyl ether at 193 nm and by the reaction F + CH3CHO, respectively. Laser-induced fluorescence excited at 337 or 347.4 nm was used to monitor the CH2CHO concentration. The reaction proceeded via reversible complex formation with subsequent isomerization and fast decomposition: CH2CHO + O2 <= => O2CH2CHO --> HO2CH2CO --> products. The rate coefficients for the first and second steps were determined (k1, k-1, k2) and analyzed by a master equation with specific rate coefficients from the Rice-Ramsperger-Kassel-Marcus (RRKM) theory. Molecular and transition-state parameters were obtained from quantum chemical calculations. A third-law analysis led to the following thermodynamic parameters for the first step: Delta(R)S degrees 300K(1) = -144 J K(-1) mol(-1) (1 bar) and Delta(R)H degrees 300K(1) = (-101 +/- 4) kJ mol(-1). From the falloff analysis, the following temperature dependencies for the low- and high-pressure limiting rate coefficients were obtained: k1(0) = 5.14 x 10(-14) exp(210 K/T) cm(-3) s(-1); k1(infinity) = 1.7 x 10(-12) exp(-520 K/T) cm(-3) s(-1); and k2(infinity) = 1.3 x 10(12) exp[-(82 +/- 4) kJ mol(-1)/RT] s(-1). Readily applicable analytical representations for the pressure and temperature dependence of k1 were derived to be used in kinetic modeling.     

Delta(DeltaS) and Delta(DeltaS) for the competing bond cleavage reactions in (CH3CN)(ROH)H+ [R = CH3, C2H5, C3H7, (CH3)2CH

Microcanonical variational transition-state theory was used to determine the entropies of activation for hydrogen-bond cleavage reactions leading to CH(3)CN + ROH(2)(+) in a series of acetonitrile-alcohol proton-bound pairs (CH(3)CN)(ROH)H(+) (where R = CH(3), CH(3)CH(2), CH(3)CH(2)CH(2), and (CH(3))(2)CH). In each case, the dissociation potential surface was modelled at the MP2/6-31 + G(d) level of theory. The dissociating configurations having the minimum sums-of-states were identified in each case and the resulting entropies of activation were calculated. Combined with previous work on the competing reaction leading to CH(3)CNH(+) + ROH, the results permitted the determination of the Delta(DeltaS) in each proton-bound pair. For the (CH(3)CN)(CH(3)OH)H(+) and (CH(3)CN)(CH(3)CH(2)OH)H(+) proton-bound pairs, the entropies of activation for the two dissociating channels are essentially the same [i.e., Delta(DeltaS) = 0], while Delta(DeltaS) for the propanol-containing pairs ranged between 40 and 45 J K(-1) mol(-1). The latter non-zero values are due to a combination of the location of the dividing surface in each dissociation and the rapidity with the frequencies of the vanishing vibrational modes go to zero as they are converted to product translations and rotations during the dissociation.     

Ab initio studies of alkyl radical reactions: Combination and disproportionation reactions of CH3 with C2H5, and the decomposition of chemically activated C3H8

This paper reports the first quantitative ab initio prediction of the disproportionation/combination ratio of alkyl+alkyl reactions using CH3+C2H5 as an example. The reaction has been investigated by the modified Gaussian-2 method with variational transition state or Rice-Ramsperger-Kassel-Marcus calculations for several channels producing (1) CH4+CH2CH2, (2) C3H8, (3) CH4CH3CH, (4) H2+CH3CHCH2, (5) H2+CH3CCH3, and (6) C2H6+CH2 by H-abstraction and association/decomposition mechanisms through singlet and triplet potential energy paths. Significantly, the disproportionation reaction (1) producing CH4+C2H4 was found to occur primarily by the lowest energy path via a loose hydrogen-bonding singlet molecular complex, H3CHC2H4, with a 3.5 kcal/mol binding energy and a small decomposition barrier (1.9 kcal/mol), instead of a direct H-abstraction process. Bimolecular reaction rate constants for the formation of the above products have been calculated in the temperature range 300-3000 K. At 1 atm, formation of C3H8 is dominant below 1200 K. Over 1200 K, the disproportionation reaction becomes competitive. The sum of products (3)-(6) accounts for less than 0.3% below 1500 K and it reaches around 1%-4% above 2000 K. The predicted rate constant for the disproportionation reaction with multiple reflections above the complex well, k1=5.04 x T(0.41) exp(429/T) at 200-600 K and k1=1.96 x 10(-20) T(2.45) exp(1470/T) cm3 molecule(-1) s(-1) at 600-3000 K, agrees closely with experimental values. Similarly, the predicted high-pressure rate constants for the combination reaction forming C3H8 and its reverse dissociation reaction in the temperature range 300-3000 K, k2(infinity)=2.41 x 10(-10) T(-0.34) exp(259/T) cm3 molecule(-1) s(-1) and k(-2)(infinity)=8.89 x 10(22) T(-1.67)exp(-46 037/T) s(-1), respectively, are also in good agreement with available experimental data.     

Kinetics and mechanism of the C6H5 + CH3CHO reaction: experimental measurement and theoretical prediction of the reactivity toward four molecular sites.

The kinetics and mechanism of the reaction of C6H5 with CH3CHO have been investigated experimentally and theoretically. The total rate constant for the reaction has been measured by means of the cavity ring-down spectrometry (CRDS) in the temperature range 299-501 K at pressures covering 20-75 Torr. The overall bimolecular rate constant can be represented by the expression k = (2.8 +/- 0.2) x 10(11) exp[-(700 +/- 30)/T] cm3 mol-1 s-1, which is slightly faster than for the analogous C6H5 + CH2O reaction determined with the same method in the same temperature range. The reaction mechanism for the C6H5 + CH3CHO reaction was also explored with quantum-chemical calculations at various hybrid density functional theories (DFTs) and using ab initio high-level composite methods. The theories predict that the reaction may occur by two hydrogen-abstraction and two addition channels with the aldehydic hydrogen-abstraction reaction being dominant. The rate constant calculated by the transition state theory for the aldehydic hydrogen-abstraction reaction is in good agreement with the experimental result after a very small adjustment of the predicted reaction barrier (+0.3 kcal mol-1). Contributions from other product channels are negligible under our experimental conditions. For combustion applications, we have calculated the rate constants for key product channels in the temperature range of 298-2500 K under atmospheric-pressure conditions; they can be represented by the following expressions in units of cm 3mol-1 s-1: k1,cho = 8.8 x 10(3)T2.6 exp(-90/T), k2,ch3 = 6.0 x 10(1)T3.3 exp(-950/T), k3a(C6H5COCH3 + H) = 4.2 x 10(5)T0.6 exp(-410/T) and k3b(C6H5CHO + CH3) = 6.6 x 10(9)T-0.5 exp(-310/T).