Reactions of Hydrocarbon Radicals and Biradicals

Predicting third-body collision efficiencies for water and other polyatomic baths

Low-pressure-limit microcanonical (collisional activation) and thermal rate constants are predicted using a combination of automated ab initio potential energy surface construction, classical trajectories, transition state theory, and a detailed kinetic...

Site-Specific Carbon-Carbon Bond Fission in Photoexcited Propyl Radicals Leads to Isomer-Selective Carbene and Radical Products

Although there have been many studies of C-H bond fission in the UV photochemistry of alkyl radicals, very little is known about the possible occurrence of C-C bond fission. Here, we report that upon excitation at 248 nm, gaseous 1-propyl radicals primarily undergo C-C bond fission, producing methylene (CH₂) and ethyl radicals (C₂H₅), rather than the more energetically favored methyl (CH₃) and ethylene (C₂H₄). In contrast, the exclusive C-C bond fission products from 2-propyl radicals were ethylidene (CHCH₃) plus methyl radicals (CH₃). The isomer-selective formation of high-energy carbene + radical products involves excited-state site-specific C-C bond fission at the radical carbon, with quantum yields comparable to those for C-H bond fission. Our observations suggest that a general feature of alkyl radical photochemistry is predissociation of the initially formed Rydberg states by high-lying valence states, yielding high-energy carbene plus alkyl radical products.

Low-Pressure Limit of Competitive Unimolecular Reactions

Barker and Ortiz found unusual falloff effects in the flux coefficients of the competitive unimolecular reactions of 2-methylhexyl radicals, and they concluded that this might have important effects on the rate constants of reactions with higher thresholds. To study this effect, we carried out master equation calculations of the same reaction system to learn whether this effect shows up in measurable rate constants, and the answer is yes. We also studied specially designed mechanisms to reveal that the various reactive pathways connecting the reagents can have a large effect on the rate constants, causing them to be quite different than if the reactions proceeded independently, and that reactions with significantly higher barriers may nevertheless have larger rate constants. This provides a new perspective for interpreting and predicting the kinetics of competitive unimolecular reactions.

Pyrolysis of Alkanes: A Computational Approach

In this study, we investigate the kinetic and thermodynamic aspects of thermal cracking reactions of short paraffin chains by density functional theory (DFT) methods. The thermal cracking reactions have been modeled for a series of shorter unbranched alkanes at 673 K by following a free-radical mechanism. Benchmark calculations have been carried out with different functionals (B3LYP, M06-2X, PBE0, BMK, B3PW91) and basis sets (6-31G(d,p), 6-311+G(d,p)) to determine the most suitable DFT method, and the results were compared to the available experimental data. Computations were also performed at the CBS-QB3 level to evaluate the accuracy of the DFT method. The thermodynamic and kinetic properties of the initiation, hydrogen atom transfer (HAT), and decomposition (β-scission) reactions are intensely discussed to better understand the trends in product distributions at high temperatures. Evans-Polanyi (EP) relations have been used to build a linear relationship between the enthalpy of reactions and their activation energies; this process may be useful for the determination of the kinetic parameters of longer paraffin chains as well. Finally, the preexponential factors of short-chain paraffin have been calculated and classified based on the identity of the radicalic products. The latter, together with the activation energies derived from the EP relations can be used safely for the prediction of the rate constants for long paraffin chains.

β-Bond Scission and the Yields of H and CH3 in the Decomposition of Isobutyl Radicals

The relative rates of C-C and C-H β-scission reactions of isobutyl radicals (2-methylprop-1-yl, C₄H₉) were investigated with shock tube experiments at temperatures of (950 to 1250) K and pressures of (200 to 400) kPa. We produced isobutyl radicals from the decomposition of dilute mixtures of isopentylbenzene and observed the stable decomposition products, propene and isobutene. These alkenes are characteristic of C-C and C-H bond scission, respectively. Propene was the main product, approximately 30 times more abundant than isobutene, indicating that C-C β-scission is the primary pathway. Uncertainty in the ratio of [isobutene]/[propene] from isobutyl decomposition is mainly due to a small amount of side chemistry, which we account for using a kinetics model based on JetSurF 2.0. Our data are well-described after adding chemistry specific to our system and adjusting some rate constants. We compare our data to other commonly used kinetics models: JetSurF 2.0, AramcoMech 2.0, and multiple models from Lawrence Livermore National Laboratory (LLNL). With the kinetics model, we have determined an upper limit of 3.0% on the branching fraction for C-H β-scission in the isobutyl radical for the temperatures and pressures of our experiments. While this agrees with previous high quality experimental results, many combustion kinetics models assume C-H branching values above this upper limit, possibly leading to large systematic inaccuracies in model predictions. Some kinetics models additionally assume contributions from 1,2-H shift reactions-which for isobutyl would produce the same products as C-H β-scission-and our upper limit includes possible involvement of such reactions. We suggest kinetics models should be updated to better reflect current experimental measurements.

Obtaining effective rate coefficients to describe the decomposition kinetics of the corannulene oxyradical at high temperatures

This work analyzes the effect of overlapping eigenvalues on the high-temperature kinetics of a large oxyradical based on master equation solutions.

Barrierless association of CF2 and dissociation of C2F4 by variational transition-state theory and system-specific quantum Rice-Ramsperger-Kassel theory

Bond dissociation is a fundamental chemical reaction, and the first principles modeling of the kinetics of dissociation reactions with a monotonically increasing potential energy along the dissociation coordinate presents a challenge not only for modern electronic structure methods but also for kinetics theory. In this work, we use multifaceted variable-reaction-coordinate variational transition-state theory (VRC-VTST) to compute the high-pressure limit dissociation rate constant of tetrafluoroethylene (C₂F₄), in which the potential energies are computed by direct dynamics with the M08-HX exchange correlation functional. To treat the pressure dependence of the unimolecular rate constants, we use the recently developed system-specific quantum Rice-Ramsperger-Kassel theory. The calculations are carried out by direct dynamics using an exchange correlation functional validated against calculations that go beyond coupled-cluster theory with single, double, and triple excitations. Our computed dissociation rate constants agree well with the recent experimental measurements.

Active Thermochemical Tables: Sequential Bond Dissociation Enthalpies of Methane, Ethane, and Methanol and the Related Thermochemistry

Active Thermochemical Tables (ATcT) thermochemistry for the sequential bond dissociations of methane, ethane, and methanol systems were obtained by analyzing and solving a very large thermochemical network (TN). Values for all possible C-H, C-C, C-O, and O-H bond dissociation enthalpies at 298.15 K (BDE298) and bond dissociation energies at 0 K (D0) are presented. The corresponding ATcT standard gas-phase enthalpies of formation of the resulting CHn, n = 4-0 species (methane, methyl, methylene, methylidyne, and carbon atom), C2Hn, n = 6-0 species (ethane, ethyl, ethylene, ethylidene, vinyl, ethylidyne, acetylene, vinylidene, ethynyl, and ethynylene), and COHn, n = 4-0 species (methanol, hydroxymethyl, methoxy, formaldehyde, hydroxymethylene, formyl, isoformyl, and carbon monoxide) are also presented. The ATcT thermochemistry of carbon dioxide, water, hydroxyl, and carbon, oxygen, and hydrogen atoms is also included, together with the sequential BDEs of CO2 and H2O. The provenances of the ATcT enthalpies of formation, which are quite distributed and involve a large number of relevant determinations, are analyzed by variance decomposition and discussed in terms of principal contributions. The underlying reasons for periodic appearances of remarkably low and/or unusually high BDEs, alternating along the dissociation sequences, are analyzed and quantitatively rationalized. The present ATcT results are the most accurate thermochemical values currently available for these species.

Theoretical investigation of the relaxation of the bending mode of CH₂(X̃) by collisions with helium

We have earlier determined the dependence on the bending angle of the interaction of the methylene radical (CH2) in its X̃³B₁ state with He [L. Ma, P. J. Dagdigian, and M. H. Alexander, J. Chem. Phys. 136, 224306 (2012)]. By integration over products of the bending vibrational wave function, in a quantum close-coupled treatment we have calculated cross sections for the ro-vibrational relaxation of CH 2(X̃). Specifically, we find that cross sections for a loss of one vibrational quantum (v(b) = 2 → 1 and 1 → 0) are roughly two orders of magnitude smaller, and those for a loss of two vibrational quanta (v(b) = 2 → 0) four orders of magnitude smaller, than those for pure rotational relaxation. In addition, no clear cut dependence on the energy gap is seen.