H-Abstraction reactions by OH, HO2, O, O2 and benzyl radical addition to O2 and their implications for kinetic modelling of toluene oxidation.

Alkylated aromatics constitute a significant fraction of the components commonly found in commercial fuels. Toluene is typically considered as a reference fuel. Together with n-heptane and iso-octane, it allows for realistic emulations of the behavior of real fuels by the means of surrogate mixture formulations. Moreover, it is a key precursor for the formation of poly-aromatic hydrocarbons, which are of relevance to understanding soot growth and oxidation mechanisms. In this study the POLIMI kinetic model is first updated based on the literature and on recent kinetic modelling studies of toluene pyrolysis and oxidation. Then, important reaction pathways are investigated by means of high-level theoretical methods, thereby advancing the present knowledge on toluene oxidation. H-Abstraction reactions by OH, HO2, O and O2, and the reactivity on the multi well benzyl-oxygen (C6H5CH2 + O2) potential energy surface (PES) were investigated using electronic structure calculations, transition state theory in its conventional, variational, and variable reaction coordinate forms (VRC-TST), and master equation calculations. Exploration of the effect on POLIMI model performance of literature rate constants and of the present calculations provides valuable guidelines for implementation of the new rate parameters in existing toluene kinetic models.

First-Principles Chemical Kinetic Modeling of Methyl trans-3-Hexenoate Epoxidation by HO2.

The design of innovative combustion processes relies on a comprehensive understanding of biodiesel oxidation kinetics. The present study aims at unraveling the reaction mechanism involved in the epoxidation of a realistic biodiesel surrogate, methyl trans-3-hexenoate, by hydroperoxy radicals using a bottom-up theoretical kinetics methodology. The obtained rate constants are in good agreement with experimental data for alkene epoxidation by HO2. The impact of temperature and pressure on epoxidation pathways involving H-bonded and non-H-bonded conformers was assessed. The obtained rate constant was finally implemented into a state-of-the-art detailed combustion mechanism, resulting in fairly good agreement with engine experiments.

Low Temperature Kinetics of the First Steps of Water Cluster Formation.

We present a combined experimental and theoretical low temperature kinetic study of water cluster formation. Water cluster growth takes place in low temperature (23-69 K) supersonic flows. The observed kinetics of formation of water clusters are reproduced with a kinetic model based on theoretical predictions for the first steps of clusterization. The temperature- and pressure-dependent association and dissociation rate coefficients are predicted with an ab initio transition state theory based master equation approach over a wide range of temperatures (20-100 K) and pressures (10^{-6}-10 bar).

Shock tube explorations of roaming radical mechanisms: the decompositions of isobutane and neopentane.

The thermal decompositions of isobutane and neopentane have been studied using both shock tube experiments and ab initio transition state theory based master equation calculations. Dissociation rate constants for these molecules have been measured at high temperatures (1260-1566 K) behind reflected shock waves using high-sensitivity H-ARAS detection. The two major dissociation channels at high temperature are iso-C(4)H(10) → CH(3) + i-C(3)H(7) (1a) and neo-C(5)H(12) → CH(3) + t-C(4)H(9) (2a). Ultrahigh-sensitivity ARAS detection of H-atoms produced from the rapid decomposition of the product radicals, i-C(3)H(7) in (1a) and t-C(4)H(9) in (2a), through i-C(3)H(7) + M → H + C(3)H(6) + M (3a) and t-C(4)H(9) + M → H + i-C(4)H(8) + M (4a) allowed measurements of both the total decomposition rate constants, k(total), and the branching to radical products, which were observed to be equivalent in both systems, k(1a)/k(total) and k(2a)/k(total) = 0.79 ± 0.05. Theoretical analyses indicate that in isobutane, the non-H-atom fraction has two contributions, the dominant fraction being due to the roaming radical mechanism leading to molecular products through iso-C(4)H(10) → CH(4) + C(3)H(6) (1b) with k(1b)/k(total) = 0.16, and a minor fraction that involves the isomerization of i-C(3)H(7) to n-C(3)H(7) that then subsequently forms methyl radicals, i-C(3)H(7) + M → n-C(3)H(7) + M → CH(3) + C(2)H(4) + M (3b). In contrast to isobutane, in neopentane, the contribution to the non-H-atom fraction is exclusively through the roaming radical mechanism that leads to neo-C(5)H(12) → CH(4) + i-C(4)H(8) (2b) with k(2b)/k(total) = 0.21. These quantitative measurements of larger contributions from the roaming mechanism for larger molecules are in agreement with the qualitative theoretical arguments that suggest long-range dispersion interactions (which become increasingly important for larger molecules) may enhance roaming.

Shock tube and theoretical studies on the thermal decomposition of propane: evidence for a roaming radical channel.

The thermal decomposition of propane has been studied using both shock tube experiments and ab initio transition state theory-based master equation calculations. Dissociation rate constants for propane have been measured at high temperatures behind reflected shock waves using high-sensitivity H-ARAS detection and CH(3) optical absorption. The two major dissociation channels at high temperature are C(3)H(8) → CH(3) + C(2)H(5) (eq 1a) and C(3)H(8) → CH(4) + C(2)H(4) (eq 1b). Ultra high-sensitivity ARAS detection of H-atoms produced from the decomposition of the product, C(2)H(5), in (1a), allowed measurements of both the total decomposition rate constants, k(total), and the branching to radical products, k(1a)/k(total). Theoretical analyses indicate that the molecular products are formed exclusively through the roaming radical mechanism and that radical products are formed exclusively through channel 1a. The experiments were performed over the temperature range 1417-1819 K and gave a minor contribution of (10 ± 8%) due to roaming. A multipass CH(3) absorption diagnostic using a Zn resonance lamp was also developed and characterized in this work using the thermal decomposition of CH(3)I as a reference reaction. The measured rate constants for CH(3)I decomposition agreed with earlier determinations from this laboratory that were based on I-atom ARAS measurements. This CH(3) diagnostic was then used to detect radicals from channel 1a allowing lower temperature (1202-1543 K) measurements of k(1a) to be determined. Variable reaction coordinate-transition state theory was used to predict the high pressure limits for channel (1a) and other bond fission reactions in C(3)H(8). Conventional transition state theory calculations were also used to estimate rate constants for other tight transition state processes. These calculations predict a negligible contribution (<1%) from all other bond fission and tight transition state processes, indicating that the bond fission channel (1a) and the roaming channel (1b) are indeed the only active channels at the temperature and pressure ranges of the present experiments. The predicted reaction exo- and endothermicities are in excellent agreement with the current version of the Active Thermochemical Tables. Master equation calculations incorporating these transition state theory results yield predictions for the temperature and pressure dependence of the dissociation rate constants for channel 1a. The final theoretical results reliably reproduce the measured dissociation rate constants that are reported here and in the literature. The experimental data are well reproduced over the 500-2500 K and 1 × 10(-4) to 100 bar range (errors of ∼15% or less) by the following Troe parameters for Ar as the bath gas: k(∞) = 1.55 × 10(24)T(-2.034) exp(-45 490/T) s(-1), k(0) = 7.92 × 10(53)T(-16.67) exp(-50 380/T) cm(3) s(-1), and F(c) = 0.190 exp(-T/3091) + 0.810 exp(-T/128) + exp(-8829/T).

Direct observation of roaming radicals in the thermal decomposition of acetaldehyde.

The thermal dissociation of acetaldehyde has been studied with the reflected shock tube technique using H(D)-atom atomic resonance absorption spectrometry detection. The use of an unreversed light source yields extraordinarily sensitive H atom detection. As a result, we are able to measure both the total decomposition rate and the branching to radical versus molecular channels. This branching provides a direct measure of the contribution from the roaming radical mechanism since the contributions from the usual tight transition states are predicted by theory to be negligible. The experimental observations also provide a measure of the rate coefficient for H + CH(3)CHO. Another set of experiments employing C(2)H(5)I as an H-atom source provides additional data for this rate coefficient that extends to lower temperature. An evaluation of the available experimental results for H + CH(3)CHO can be expressed by a three-parameter Arrhenius expression as k = 7.66 x 10(-20)T(2.75) exp((-486 K)/T) cm(3) molecule(-1) s(-1) (298-1415 K). Analogous experiments employing C(2)D(5)I as a D-atom source allow for the study of the isotopically substituted reaction. The present experiments are the only direct measure for this reaction rate constant, and the results can be expressed by an Arrhenius expression as k = 5.20 x 10(-10) exp((-4430 K)/T) cm(3) molecule(-1) s(-1) (1151-1354 K). The H/D + CH(3)CHO reactions are also studied with ab initio transition-state theory, and the results are in remarkably good agreement with the current experimental data.

D-atom products in predissociation of CD2CD2OH from the 202-215 nm photodissociation of 2-bromoethanol.

Experimental observations of D fragments from the predissociation of rovibrationally excited partially deuterated 2-hydroxyethyl radicals, CD(2)CD(2)OH, are reported, and possible dissociation channels are analyzed by theory. The radicals are produced by photolysis of 2-bromoethanol at 202-215 nm, and some of them have sufficient internal energy to predissociate. D fragments are detected by 1 + 1' REMPI and their TOF distributions are determined. They can be associated with vinyl alcohol and/or acetaldehyde cofragments. From analysis of the maximum velocities and kinetic energies of the observed D fragments it is concluded that they originate from the decomposition of CD(2)CD(2)OH, but the experimental resolution is insufficient to distinguish between the two possible channels leading to D products. Theoretical analysis and RRKM calculations of microcanonical dissociation rates and branching ratios for the range of available excess energies (up to 5000-8000 cm(-1) above the OH + C(2)D(4) threshold) indicate that the D-producing channels are minor (about 1%) compared to the predominant OH + C(2)D(4) channel, and the branching ratio for D production is more favorable when the reactant radicals have low rotational energy. The vinyl alcohol channel is strongly favored over the acetaldehyde channel at all excess energies, except near the threshold of these channels.

Rate constants for the thermal decomposition of ethanol and its bimolecular reactions with OH and D: reflected shock tube and theoretical studies.

The thermal decomposition of ethanol and its reactions with OH and D have been studied with both shock tube experiments and ab initio transition state theory-based master equation calculations. Dissociation rate constants for ethanol have been measured at high T in reflected shock waves using OH optical absorption and high-sensitivity H-atom ARAS detection. The three dissociation processes that are dominant at high T are C2H5OH--> C2H4+H2O (A) -->CH3+CH2OH (B) -->C2H5+OH (C).The rate coefficient for reaction C was measured directly with high sensitivity at 308 nm using a multipass optical White cell. Meanwhile, H-atom ARAS measurements yield the overall rate coefficient and that for the sum of reactions B and C , since H-atoms are instantaneously formed from the decompositions of CH(2)OH and C(2)H(5) into CH(2)O + H and C(2)H(4) + H, respectively. By difference, rate constants for reaction 1 could be obtained. One potential complication is the scavenging of OH by unreacted ethanol in the OH experiments, and therefore, rate constants for OH+C2H5OH-->products (D)were measured using tert-butyl hydroperoxide (tBH) as the thermal source for OH. The present experiments can be represented by the Arrhenius expression k=(2.5+/-0.43) x 10(-11) exp(-911+/-191 K/T) cm3 molecule(-1) s(-1) over the T range 857-1297 K. For completeness, we have also measured the rate coefficient for the reaction of D atoms with ethanol D+C2H5OH-->products (E) whose H analogue is another key reaction in the combustion of ethanol. Over the T range 1054-1359 K, the rate constants from the present experiments can be represented by the Arrhenius expression, k=(3.98+/-0.76) x10(-10) exp(-4494+/-235 K/T) cm3 molecule(-1) s(-1). The high-pressure rate coefficients for reactions B and C were studied with variable reaction coordinate transition state theory employing directly determined CASPT2/cc-pvdz interaction energies. Reactions A , D , and E were studied with conventional transition state theory employing QCISD(T)/CBS energies. For the saddle point in reaction A , additional high-level corrections are evaluated. The predicted reaction exo- and endothermicities are in good agreement with the current Active Thermochemical Tables values. The transition state theory predictions for the microcanonical rate coefficients in ethanol decomposition are incorporated in master equation calculations to yield predictions for the temperature and pressure dependences of reactions A - C . With modest adjustments (<1 kcal/mol) to a few key barrier heights, the present experimental and adjusted theoretical results yield a consistent description of both the decomposition (1-3) and abstraction kinetics (4 and 5). The present results are compared with earlier experimental and theoretical work.

Shock tube and theory investigation of cyclohexane and 1-hexene decomposition.

The decomposition of cyclohexane (c-C(6)H(12)) was studied in a shock tube using the laser-schlieren technique over the temperature range 1300-2000 K and for 25-200 Torr in mixtures of 2%, 4%, 10%, and 20% cyclohexane in Kr. Vibrational relaxation of the cyclohexane was also examined in 10 experiments covering 1100-1600 K for pressures below 20 Torr, and relaxation was found to be too fast to allow resolution of incubation times. The dissociation of 1-hexene (1- C(6)H(12)), apparently the sole initial product of cyclohexane decomposition, was also studied over 1220-1700 K for 50 and 200 Torr using 2% and 3% 1-hexene in Kr. On heating, cyclohexane simply isomerizes to 1-hexene, and this then dissociates almost entirely by a more rapid C-C scission to allyl and n-propyl radicals. This two-step reaction results in an initial small density gradient from the slight endothermicity of the isomerization. The gradient then rises strongly as the product 1-hexene dissociates. For the lower temperatures, this behavior is fully resolved here. For the higher pressures, 1-hexene decomposition generates negative gradients (exothermic reaction) as the radicals formed begin to recombine. Cyclohexane also generates such gradients, but these are now much smaller because the radical pool is depleted by abstraction from the reactant. A complete mechanism for the 1-hexene decomposition and for that of cyclohexane involving 79 reactions and 30 species is used in the final modeling of the gradients. Rate constants and RRKM fit parameters for the initial reactions are provided for the entire range of conditions. The possibility of direct reaction to allyl and n-propyl radicals, without stabilization of the intermediate 1-hexene, is examined down to pressures as low as 25 Torr, without a clear resolution of the issue. High-pressure limit rate constants from RRKM extrapolation are k(infinity)(c-C(6)H(12) --> 1-C(6)H(12)) = (8.76 x 10(17)) exp((-91.94 kcal/mol)/RT) s(-1) (T = 1300-2000 K) and k(infinity)(1-C(6)H(12) --> (*)C(3)H(7) + (*)C(3)H(5)) = (1.46 x 10(16)) exp((-69.12 kcal/mol)/RT) s(-1) (T = 1200-1700 K). This high-pressure rate for cyclohexane is entirely consistent with the notion that the isomerization involves initial C-C fission to a diradical. These extrapolated high-pressure rates are in good agreement with much of the literature.

Thermal decomposition of NH2OH and subsequent reactions: ab initio transition state theory and reflected shock tube experiments.

Primary and secondary reactions involved in the thermal decomposition of NH2OH are studied with a combination of shock tube experiments and transition state theory based theoretical kinetics. This coupled theory and experiment study demonstrates the utility of NH2OH as a high temperature source of OH radicals. The reflected shock technique is employed in the determination of OH radical time profiles via multipass electronic absorption spectrometry. O-atoms are searched for with atomic resonance absorption spectrometry. The experiments provide a direct measurement of the rate coefficient, k1, for the thermal decomposition of NH2OH. Secondary rate measurements are obtained for the NH2 + OH (5a) and NH2OH + OH (6a) abstraction reactions. The experimental data are obtained for temperatures in the range from 1355 to 1889 K and are well represented by the respective rate expressions: log[k/(cm3 molecule(-1) s(-1))] = (-10.12 +/- 0.20) + (-6793 +/- 317 K/T) (k1); log[k/(cm3 molecule(-1) s(-1))] = (-10.00 +/- 0.06) + (-879 +/- 101 K/T) (k5a); log[k/(cm3 molecule(-1) s(-1))] = (-9.75 +/- 0.08) + (-1248 +/- 123 K/T) (k6a). Theoretical predictions are made for these rate coefficients as well for the reactions of NH2OH + NH2, NH2OH + NH, NH + OH, NH2 + NH2, NH2 + NH, and NH + NH, each of which could be of secondary importance in NH2OH thermal decomposition. The theoretical analyses employ a combination of ab initio transition state theory and master equation simulations. Comparisons between theory and experiment are made where possible. Modest adjustments of predicted barrier heights (i.e., by 2 kcal/mol or less) generally yield good agreement between theory and experiment. The rate coefficients obtained here should be of utility in modeling NOx in various combustion environments.

EStokTP: Electronic Structure to Temperature- and Pressure-Dependent Rate Constants-A Code for Automatically Predicting the Thermal Kinetics of Reactions.

A priori rate predictions for gas phase reactions have undergone a gradual but dramatic transformation, with current predictions often rivaling the accuracy of the best available experimental data. The utility of such kinetic predictions would be greatly magnified if they could more readily be implemented for large numbers of systems. Here, we report the development of a new computational environment, namely, EStokTP, that reduces the human effort involved in the rate prediction for single channel reactions essentially to the specification of the methodology to be employed. The code can also be used to obtain all the necessary master equation building blocks for more complex reactions. In general, the prediction of rate constants involves two steps, with the first consisting of a set of electronic structure calculations and the second in the application of some form of kinetic solver, such as a transition state theory (TST)-based master equation solver. EStokTP provides a fully integrated treatment of both steps through calls to external codes to perform first the electronic structure and then the master equation calculations. It focuses on generating, extracting, and organizing the necessary structural properties from a sequence of calls to electronic structure codes, with robust automatic failure recovery options to limit human intervention. The code implements one or multidimensional hindered rotor treatments of internal torsional modes (with automated projection from the Hessian and with optional vibrationally adiabatic corrections), Eckart and multidimensional tunneling models (such as small curvature theory), and variational treatments (based on intrinsic reaction coordinate following). This focus on a robust implementation of high-level TST methods allows the code to be used in high accuracy studies of large sets of reactions, as illustrated here through sample studies of a few hundred reactions. At present, the following reaction types are implemented in EStokTP: abstraction, addition, isomerization, and beta-decomposition. Preliminary protocols for treating barrierless reactions and multiple-well and/or multiple-channel potential energy surfaces are also illustrated.

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.

Formation of NH3 and CH2NH in Titan's upper atmosphere.

The large abundance of NH3 in Titan's upper atmosphere is a consequence of coupled ion and neutral chemistry. The density of NH3 is inferred from the measured abundance of NH4+. NH3 is produced primarily through reaction of NH2 with H2CN, a process neglected in previous models. NH2 is produced by several reactions including electron recombination of CH2NH2+. The density of CH2NH2+ is closely linked to the density of CH2NH through proton exchange reactions and recombination. CH2NH is produced by reaction of N(2D) and NH with ambient hydrocarbons. Thus, production of NH3 is the result of a chain of reactions involving non-nitrile functional groups and the large density of NH3 implies large densities for these associated molecules. This suggests that amine and imine functional groups may be incorporated as well in other, more complex organic molecules.

A combined ab initio and photoionization mass spectrometric study of polyynes in fuel-rich flames.

Polyynic structures in fuel-rich low-pressure flames are observed using VUV photoionization molecular-beam mass spectrometry. High-level ab initio calculations of ionization energies for C2nH2 (n=1-5) and partially hydrogenated CnH4 (n=7-8) polyynes are compared with photoionization efficiency measurements in flames fuelled by allene, propyne, and cyclopentene. C2nH2 (n=1-5) intermediates are unambiguously identified, while HC[triple bond, length as m-dash]C-C[triple bond, length as m-dash]C-CH=C=CH2, HC[triple bond, length as m-dash]C-C[triple bond, length as m-dash]C-C[triple bond, length as m-dash]C-CH=CH2 (vinyltriacetylene) and HC[triple bond, length as m-dash]C-C[triple bond, length as m-dash]C-CH[double bond, length as m-dash]CH-C[triple bond, length as m-dash]CH are likely to contribute to the C7H4 and C8H4 signals. Mole fraction profiles as a function of distance from the burner are presented. C7H4 and C8H4 isomers are likely to be formed by reactions of C2H and C4H radicals but other plausible formation pathways are also discussed. Heats of formation and ionization energies of several combustion intermediates have been determined for the first time.

Thermal decomposition of CF3 and the reaction of CF2 + OH --> CF2O + H.

The reflected shock tube technique with multipass absorption spectrometric detection (at a total path length of approximately 1.75 m) of OH-radicals at 308 nm has been used to study the dissociation of CF3-radicals [CF3 + Kr --> CF2 + F + Kr (a)] between 1,803 and 2,204 K at three pressures between approximately 230 and 680 Torr. The OH-radical concentration buildup resulted from the fast reaction F + H2O --> OH + HF (b). Hence, OH is a marker for F-atoms. To extract rate constants for reaction (a), the [OH] profiles were modeled with a chemical mechanism. The initial rise in [OH] was mostly sensitive to reactions (a) and (b), but the long time values were additionally affected by CF2 + OH --> CF2O + H (c). Over the experimental temperature range, rate constants for (a) and (c) were determined from the mechanistic fits to be kCF3+Kr = 4.61 x 10-9 exp(-30,020 K/T) and kCF2+OH = (1.6 +/- 0.6) x 10-10, both in units of cm3 molecule-1 s-1. Reaction (a), its reverse recombination reaction reaction (-a), and reaction (c) are also studied theoretically. Reactions (c) and (-a) are studied with direct CASPT2 variable reaction coordinate transition state theory. A master equation analysis for reaction (a) incorporating the ab initio determined reactive flux for reaction (-a) suggests that this reaction is close to but not quite in the low-pressure limit for the pressures studied experimentally. In contrast, reaction (c) is predicted to be in the high-pressure limit due to the high exothermicity of the products. A comparison with past and present experimental results demonstrates good agreement between the theoretical predictions and the present data for both (a) and (c).

Reflected shock tube and theoretical studies of high-temperature rate constants for OH+CF3H<-->CF3+H2O and CF3+OH-->products.

The reflected shock tube technique with multipass absorption spectrometric detection of OH radicals at 308 nm, using either 36 or 60 optical passes corresponding to total path lengths of 3.25 or 5.25 m, respectively, has been used to study the bimolecular reactions, OH+CF3H-->CF3+H2O (1) and CF3+H2O-->OH+CF3H (-1), between 995 and 1663 K. During the course of the study, estimates of rate constants for CF3+OH-->products (2) could also be determined. Experiments on reaction -1 were transformed through equilibrium constants to k1, giving the Arrhenius expression k1=(9.7+/-2.1)x10(-12) exp(-4398+/-275K/T) cm3 molecule(-1) s(-1). Over the temperature range, 1318-1663 K, the results for reaction 2 were constant at k2=(1.5+/-0.4)x10(-11) cm3 molecule(-1) s(-1). Reactions 1 and -1 were also studied with variational transition state theory (VTST) employing QCISD(T) properties for the transition state. These a priori VTST predictions were in good agreement with the present experimental results but were too low at the lower temperatures of earlier experiments, suggesting that either the barrier height was overestimated by about 1.3 kcal/mol or that the effect of tunneling was greatly underestimated. The present experimental results have been combined with the most accurate earlier studies to derive an evaluation over the extended temperature range of 252-1663 K. The three parameter expression k1=2.08x10(-17) T1.5513 exp(-1848 K/T) cm3 molecule(-1) s(-1) describes the rate behavior over this temperature range. Alternatively, the expression k1,th=1.78x10(-23) T3.406 exp(-837 K/T) cm3 molecule(-1) s(-1) obtained from empirically adjusted VTST calculations over the 250-2250 K range agrees with the experimental evaluation to within a factor of 1.6. Reaction 2 was also studied with direct CASPT2 variable reaction coordinate transition state theory. The resulting predictions for the capture rate are found to be in good agreement with the mean of the experimental results and can be represented by the expression k2,th=2.42x10(-11) T-0.0650 exp(134 K/T) cm3 molecule(-1) s(-1) over the 200-2500 K temperature range. The products of this reaction are predicted to be CF2O+HF.

Initial steps of aromatic ring formation in a laminar premixed fuel-rich cyclopentene flame.

A fuel-rich, nonsooting, premixed laminar cyclopentene flame (phi = 2.0) at 37.6 Torr (50 mbar) is investigated by flame-sampling photoionization molecular-beam mass spectrometry utilizing vacuum-ultraviolet synchrotron radiation. Mole fractions as a function of distance from the burner are measured for 49 intermediates with ion masses ranging from 2 (H2) to 106 (C8H10), providing a broad database for flame modeling studies. The isomeric composition is resolved for most species, and the identification of several C4Hx, C7H6, and C7H8 isomers is discussed in detail. The presence of C5H5CCH/C5H4CCH2 and cycloheptatriene is revealed by comparisons between flame-sampled photoionization efficiency data and theoretical simulations, based on calculated ionization energies and Franck-Condon factors. This insight suggests a new potential molecular- weight growth mechanism that is characterized by C5-C7 ring enlargement reactions.

A theoretical analysis of the reaction between propargyl and molecular oxygen.

The temperature- and pressure-dependent kinetics of the reaction between propargyl and molecular oxygen have been studied with a combination of electronic structure theory, transition state theory, and the time-dependent master equation. The stationary points on the potential energy surface were located with B3LYP density functional theory. Approximate QCISD(T,Full)/6-311++G(3df,2pd) energies were obtained at these stationary points. At low temperatures the reaction is dominated by addition to the CH2 side of the propargyl radical followed by stabilization. However, addition to the CH side, which is followed by one of various possible internal rearrangements, becomes the dominant process at higher temperatures. These internal rearrangements involve a splitting of the O2 bond via the formation of 3-, 4- or 5-membered rings, with the apparent products being CH2CO + HCO. Rearrangement via the 3-membered ring is found to dominate the kinetics. Rearrangement from the CH2 addition product, via a 4-membered ring, would yield H2CO + HCCO, but the barrier to this rearrangement is too high to be kinetically significant. Other possible products require H transfers and, as a result, appear to be kinetically irrelevant. Modest variations in the energetics of a few key stationary points (most notably the entrance barrier heights) yield kinetic results that are in good agreement with the experimental results of Slagle and Gutman (I. R. Slagle and D. Gutman, Proc. Combust. Inst., 1986, 21, 875) and of Atkinson and Hudgens (D. B. Atkinson and J. W. Hudgens, J. Phys. Chem. A, 1999, 103, 4242).

A direct transition state theory based analysis of the branching in NH2 + NO.

A combination of high-level quantum-chemical simulations and sophisticated transition state theory analyses is employed in a study of the temperature dependence of the N2H + OH-->HNNOH recombination reaction. The implications for the branching between N2H + OH and N2 + H2O in the NH2 + NO reaction are also explored. The transition state partition function for the N2H + OH recombination reaction is evaluated with a direct implementation of variable reaction coordinate (VRC) transition state theory (TST). The orientation dependent interaction energies are directly determined at the CAS + 1 + 2/cc-pvdz level. Corrections for basis set limitations are obtained via calculations along the cis and trans minimum energy paths employing an approximately aug-pvtz basis set. The calculated rate constant for the N2H + OH-->HNNOH recombination is found to decrease significantly with increasing temperature, in agreement with the predictions of our earlier theoretical study. Conventional transition state theory analyses, employing new coupled cluster estimates for the vibrational frequencies and energies at the saddlepoints along the NH2 + NO reaction pathway, are coupled with the VRC-TST analyses for the N2H + OH channels to provide estimates for the branching in the NH2 + NO reaction. Modest variations in the exothermicity of the reaction (1-2 kcal mol-1), and in a few of the saddlepoint energies (2-4 kcal mol-1), yield TST based predictions for the branching fraction that are in satisfactory agreement with related experimental results. The unmodified results are in reasonable agreement for higher temperatures, but predict too low a branching ratio near room temperature, as well as too steep an initial rise.

Infrared frequency-modulation probing of product formation in alkyl + O2 reactions. Part IV. Reactions of propyl and butyl radicals with O2.

The time-resolved production of HO2 in the Cl-initiated oxidation of iso- and n-butane is measured using continuous-wave (CW) infrared frequency modulation spectroscopy between 298 and 693 K. The yield of HO2 is determined relative to the Cl2/CH3OH/O2 system. As in studies of smaller alkanes, the branching fraction to HO2 + alkene in butyl + O2 displays a dramatic rise with increasing temperature between about 550 and 700 K (the "transition region") which is accompanied by a qualitative change in the time behavior of the HO2 production. At low temperatures the HO2 is formed promptly; a second, slower production of HO2 is responsible for the bulk of the increased yield in the transition temperature region. In contrast to reactions of smaller alkyl radicals with O2, the total HO2 yield in the butyl radical reactions appears to remain significantly below 1 up to 700 K, implying a significant role for OH-producing channels. The slower HO2 production in butane oxidation displays an apparent activation energy similar to that measured for smaller alkyl + O2 reactions, suggesting that the energetics of the HO2 elimination transition state are similar for a broad range of R + O2 systems. A combination of QCISD(T) based characterizations of the propyl and butyl + O2 potential energy surfaces and master equation based characterization of the propyl + O2 kinetics provide the framework for explanation of the experimentally observed HO2 production in Cl-initiated propane and butane oxidation. These calculations suggest that the HO2 elimination channel is similar in all reaction systems, and that hydroperoxyalkyl (QOOH) species produced by internal H-atom abstraction in RO2 can provide a path to OH formation. However, the QOOH formed by the energetically favorable 1,5 isomerization (via a six-membered ring transition state) generally experiences significant barriers (relative to the radical + O2 reactants) to the production of an oxetane + OH. In contrast, the barriers to forming OH + an oxirane or an oxolane, via 1,4 or 1,6 isomerizations, respectively, are generally below reactants.