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 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.

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.

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.

The Roaming Atom: Straying from the Reaction Path in Formaldehyde Decomposition

We present a combined experimental and theoretical investigation of formaldehyde (H2CO) dissociation to H2 and CO at energies just above the threshold for competing H elimination. High-resolution state-resolved imaging measurements of the CO velocity distributions reveal two dissociation pathways. The first proceeds through a well-established transition state to produce rotationally excited CO and vibrationally cold H2. The second dissociation pathway yields rotationally cold CO in conjunction with highly vibrationally excited H2. Quasi-classical trajectory calculations performed on a global potential energy surface for H2CO suggest that this second channel represents an intramolecular hydrogen abstraction mechanism: One hydrogen atom explores large regions of the potential energy surface before bonding with the second H atom, bypassing the saddle point entirely.

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.

Mapping the OH + CO-->HOCO reaction pathway through IR spectroscopy of the OH-CO reactant complex

A hydrogen-bonded complex composed of the OH and CO reactants has been identified along the OH + CO-->HOCO reaction pathway. IR action spectroscopy in the OH overtone region has been used to examine the vibrational modes of the linear OH-CO complex, including intermolecular bending modes that probe portions of the reaction path leading to HOCO. The spectroscopic measurements have accessed highly excited intermolecular levels, with energies up to 250 cm-1 above the zero-point level, which lie in close proximity to the transition state for reaction. The OH-CO binding energy, D0 < or = 430 cm-1, has also been established from the quantum state distribution of the OH fragments following vibrational predissociation of the OH-CO complex. Complementary electronic structure calculations have been performed to characterize the OH-CO and OH-OC complexes, the transition state for HOCO formation, and the direct reaction path that connects the experimentally observed OH-CO complex to the HOCO intermediate.

Theoretical Studies of the Energetics and Dynamics of Chemical Reactions

Computational studies of basic chemical processes not only provide numbers for comparison with experiment or for use in modeling complex chemical phenomena such as combustion, but also provide insight into the fundamental factors that govern molecular structure and change which cannot be obtained from experiment alone. We summarize the results of three case studies, on HCO, OH + H(2), and O + C(2)H(2), which illustrate the range of problems that can be addressed by using modern theoretical techniques. In all cases, the potential energy surfaces were characterized by using ab initio electronic structure methods. Collisions between molecules leading to reaction or energy transer were described with quantum dynamical methods (HCO), classical trajectory techniques (HCO and OH + H(2)), and statistical methods (HCO, OH + H(2), and O + C(2)H(2)). We can anticipate dramatic increases in the scope of this work as new generations of computers are introduced and as new chemistry software is developed to exploit these computers.