Reactivity–Structure-Based Rate Estimation Rules for Alkyl Radical H Atom Shift and Alkenyl Radical Cycloaddition Reactions

High-resolution infrared spectroscopy of jet cooled cyclobutyl in the α-CH stretch region: large-amplitude puckering dynamics in a 4-membered ring radical

Gas-phase cyclobutyl radical (c-C4H7) is generated at a rotational temperature of Trot = 26(1) K in a slit-jet discharge mixture of 70% Ne/30% He and 0.5-0.6% cyclobromobutane (c-C4H7Br). A fully rovibrationally resolved absorption spectrum of the α-CH stretch fundamental band between 3062.9 cm-1 to 3075.7 cm-1 is obtained and analyzed, yielding first precision structural and dynamical information for this novel radical species. The α-CH stretch band origin is determined to be 3068.7887(4) cm-1, which implies only a modest (≈0.8 cm-1) blue shift from rotationally unresolved infrared spectroscopic studies of cyclobutyl radicals in liquid He droplets [A. R. Brown, P. R. Franke and G. E. Douberly, J. Phys. Chem. A, 2017, 121, 7576-7587]. Of particular dynamical interest, a one-dimensional potential energy surface with respect to the ring puckering coordinate is computed at CCSD(T)/ANO2 level of theory and reveals a double minimum Cs puckered geometry, separated by an exceedingly shallow planar C2v transition state barrier (Ebarr ≈ 1 cm-1). Numerical solutions on this double minimum potential yield a zero-point energy for the ground state (Ezero-point ≈ 27 cm-1) greatly in excess of the interconversion barrier. This is indicative of highly delocalized, large amplitude motion of the four-membered ring structure, for which proper vibrationally averaging of the moment of inertia tensor reproduces the experimentally determined inertial defect remarkably well. Finally, intensity alternation in the experimental spectrum due to nuclear spin statistics upon exchange of three indistinguishable H atom pairs (IH = ½) matches Ka + Kc = even : odd = 36 : 28 predictions, implying that the unpaired electron in the radical center lies in an out-of-plane pπ orbital. Thus, high-resolution infrared spectroscopy provides first experimental confirmation of a shallow double minimum ring puckering potential with a highly delocalized ground state wave function peaked at a planar C2v transition state geometry consistent with a cyclobutyl π radical.

High-Pressure-Limit Rate Coefficients for HO2 Elimination Reactions of Hydroperoxyalkenylperoxy Radicals based on the Reaction Class Transition State Theory

Thermokinetic parameters and transport parameters are of great importance to the combustion model and the reaction rate rules are of great importance to construct the combustion reaction mechanism for hydrocarbon fuels. The HO₂ elimination reaction class for hydroperoxyalkenylperoxy radicals is one of the key reaction classes for olefin, for which the rate coefficients are lacking. Therefore, the rate coefficients and rate rules of the HO₂ elimination reaction class for hydroperoxyalkenylperoxy radicals are studied in this work. The reaction class transition state theory (RC-TST) is used to calculate the rate coefficients. In addition, the HO₂ elimination reaction class of hydroperoxyalkenylperoxy radicals is divided into four subclasses depending upon the type of H-C_β bond that is broken in the reactant molecules, and the rate rules are calculated by taking the average of rate coefficients from a representative set of reactions in a subclass. The calculated kinetics data would be valuable for the construction of the combustion reaction mechanism for olefin.

High-Pressure-Limit and Pressure-Dependent Rate Rules for Unimolecular Reactions Related to Hydroperoxy Alkyl Radicals in Normal-Alkyl Cyclohexane Combustion. 2. Cyclization Reaction Class

The hydroperoxy alkyl radicals are important intermediates in the low-temperature combustion for normal-alkyl cycloalkanes, and the cyclization reactions of hydroperoxy alkyl radicals to form cyclic ethers are responsible for a major fraction of the OH formation, which has the potential to promote ignition. In most of the previous modeling studies for normal-alkyl cycloalkane combustion, the kinetic data of the cyclization reactions in the detailed combustion mechanism were mainly taken from the analogous reactions in cyclohexane, methyl cyclohexane, and alkanes in published literature studies. In this work, the kinetics of the cyclization reaction class of hydroperoxy alkyl radicals in normal-alkyl cycloalkanes is studied, where the reaction class is divided into subclasses depending upon the ring size of the transition states, the types of the carbons on which the -OOH site is located and the types of the carbons on which the radical site is located, and the positions of the cyclization (on the alkyl side chain, on the cycle, or between the alkyl side chain and the cycle). Energy barriers and high-pressure-limit site rate constants and pressure-dependent rates for reactions in all subclasses are calculated, and rate rules for all subclasses are developed. The high-pressure-limit rate constants are determined from CBS-QB3 electronic structure calculations combined with canonical transition-state theory calculations, and pressure-dependent rate constants are calculated by using the Rice-Ramsberger-Kassel-Marcus/Master Equation theory at pressures varying from 0.01 to 100 atm. Comparisons of the rate constants for cyclization reactions of hydroperoxy alkyl cyclohexylperoxy radicals calculated in this work with the values of the corresponding reactions in some of the popular combustion mechanisms show that it is unreasonable to use the kinetic data of analogous reactions in alkanes, cyclohexanes, or smaller normal-alkyl cyclohexanes. Therefore, the accurate kinetic calculations and the construction of rate rules for normal-alkyl cycloalkanes are necessary and significant for the reliable modeling of the low-temperature combustion of normal-alkyl cyclohexanes.

High-Pressure-Limit and Pressure-Dependent Rate Rules for Unimolecular Reactions Related to Hydroperoxy Alkyl Radicals in Normal Alkyl Cyclohexane Combustion. 1. Concerted HO2 Elimination Reaction Class and β-Scission Reaction Class

The reactions of the concerted HO₂ elimination from alkyl peroxy radicals and the β-scission of the C-OOH bond from hydroperoxy alkyl radicals, which lead to the formation of olefins and HO₂ radicals, are two important reaction classes that compete with the second O₂ addition step of hydroperoxy alkyl radicals, which are responsible for the chain branching in the low-temperature oxidation of normal alkyl cycloalkanes. These two reaction classes are also believed to be responsible for the negative temperature coefficient behavior due to the formation of the relatively unreactive HO₂ radical, which has the potential to inhibit ignition of normal alkyl cycloalkanes. In this work, the kinetics of the above two reaction classes in normal alkyl cycloalkanes are studied, where reactions in the concerted elimination class are divided into subclasses depending upon the types of carbons from which the H atom is eliminated and the positions of the reaction center (on the alkyl side chain or on the cycle), and the reactions in the β-scission reaction class are divided into subclasses depending upon the types of the carbons on which the radical is located and the positions of the reaction center. Energy barriers by using quantum chemical methods at the CBS-QB3 level, high-pressure-limit rate constants by using canonical transition state theory, and pressure-dependent rate constants at pressures from 0.01 to 100 atm by using Rice-Ramsberger-Kassel-Marcus/Master Equation theory are calculated for a representative set of reactions from methyl cyclohexane to n-butyl cyclohexane in each subclass, from which high-pressure-limit rate rules and pressure-dependent rate rules for each subclass are derived from the average rate constants of reactions within each subclass. A comparison of the rate constants for the reactions in the two reaction classes calculated in this work is made with the rate constants of the same reactions from available mechanisms published in the literature, where most of the rate constants are approximately estimated from analogous reactions in alkanes or small alkyl cyclohexanes, and it is found that a large difference may exist between them, indicating that the present work, which provides more accurate kinetic parameters and reasonable rate rules for these reaction classes, can be helpful to construct higher-accuracy mechanism models for normal alkyl cyclohexane combustion.

An ab Initio/Transition State Theory Study of the Reactions of Ċ5H9 Species of Relevance to 1,3-Pentadiene, Part II: Pressure Dependent Rate Constants and Implications for Combustion Modeling

The temperature- and pressure-dependence of rate constants for several radicals and unsaturated hydrocarbons reactions (1,3-C₅H₈/1,4-C₅H₈/cyC₅H₈ + Ḣ, C₂H₄ + Ċ₃H₅-a, C₃H₆ + Ċ₂H₃) are analyzed in this paper. The abstraction reactions of these systems are also calculated and compared with available literature data. Ċ₅H₉ radicals can be produced via Ḣ atom addition reactions to the pentadiene isomers and cyclopentene, and also by H-atom abstraction reactions from 1- and 2-pentene and cyclopentane. Comprehensive Ċ₅H₉ potential energy surface (PES) analyses and high-pressure limiting rate constants for related reactions have been explored in part I of this work ( J. Phys. Chem. A 2019, 123 (22), 9019-9052). In this work, a chemical kinetic model is constructed based on the computed thermochemistry and high-pressure limiting rate constants from part I, to further understand the chemistry of different C₅H₈ molecules. The most important channels for these addition reactions are discussed in the present work based on reaction pathway analyses. The dominant reaction pathways for these five systems are combined together to generate a simplified Ċ₅H₉ PES including nine reactants, 25 transition states (TSs), and nine products. Spin-restricted single point energies are calculated for radicals and TSs on the simplified PES at the ROCCSD(T)/aug-cc-pVTZ level of theory with basis set corrections from MP2/aug-cc-pVXZ (where X = T and Q). Temperature- and pressure-dependent rate constants are calculated using RRKM theory with a Master Equation analysis, with restricted energies used for minima on the simplified Ċ₅H₉ PES and unrestricted energies for other species, over a temperature range of 300-2000 K and in the pressure range 0.01-100 atm. The rate constants calculated are in good agreement with those in the literature. The chemical kinetic model is updated with pressure-dependent rate constants and is used to simulate the species concentration profiles for Ḣ atom addition to cyclopentane and cyclopentene. Through detailed analyses and comparisons, this model can reproduce the experimental measurements of species qualitatively and quantitatively with reasonably good agreement.

Theoretical Study of Radical-Molecule Reactions with Negative Activation Energies in Combustion: Hydroxyl Radical Addition to Alkenes

Many of the radical-molecule reactions are nonelementary reactions with negative activation energies, which usually proceed through two steps. They exist extensively in the atmospheric chemistry and hydrocarbon fuel combustion, so they are extensively studied both theoretically and experimentally. At the same time, various models, such as a two transition state model, a steady-state model, an equilibrium-state model, and a direct elementary dynamics model are proposed to get the kinetic parameters for the overall reaction. In this paper, a conversion temperature T _C1 is defined as the temperature at which the standard molar Gibbs free energy change of the formation of the reaction complex is equal to zero, and it is found that when T ≫ T _C1, the direct elementary dynamics model with an inclusion of the tunneling correction of the second step reaction is applicable to calculate the overall reaction rate constants for this kind of reaction system. The reaction class of hydroxyl radical addition to alkenes is chosen as the objects of this study, five reactions are chosen as the representative for the reaction class, and their single-point energies are calculated using the method of CCSD(T)/CBS, and it is shown that the highest conversion temperature for the five reactions is 139.89 K, far below the usual initial low-temperature (550 K) oxidation chemistry of hydrocarbon fuels; therefore, the steady-state approximation method is applicable. All geometry optimizations are performed at the BH&HLYP/6-311+G(d,p) level, and the result shows that the geometric parameters in the reaction centers are conserved; hence, the isodesmic reaction method is applicable to this reaction class. To validate the accuracy of this scheme, a comparison of electronic energy difference at the BH&HLYP/6-311+G(d,p) level and the corrected electronic energy difference with the electronic energy difference at the CCSD(T)/CBS level is performed for the five representative reactions, and it is shown that the maximum absolute deviation of electronic energy difference can be reduced from 2.54 kcal·mol^-1 before correction to 0.58 kcal·mol^-1 after correction, indicating that the isodesmic reaction method is applicable for the accurate calculation of the kinetic parameters for large-size molecular systems with a negative activation energy reaction. The overall rate constants for 44 reactions of the reaction class of hydroxyl radical addition to alkenes are calculated using the transition-state theory in combination with the isodesmic correction scheme, and high-pressure limit rate rules for the reaction class are developed. In addition, the thermodynamic parameter is calculated and the results indicate that our dynamics model is applicable for our studied reaction class. A chemical kinetic modeling and sensitivity analysis using the calculated kinetic data is performed for the combustion of ethene, and the results indicate the studied reaction is important for the low-to-medium temperature combustion modeling of ethene.

Ab Initio/Transition-State Theory Study of the Reactions of Ċ5H9 Species of Relevance to 1,3-Pentadiene, Part I: Potential Energy Surfaces, Thermochemistry, and High-Pressure Limiting Rate Constants

In this study, the reactions of Ċ₅H₉ radicals are theoretically investigated, with a particular emphasis on hydrogen atom addition reactions to 1,3-pentadiene (C₅H₈) to form Ċ₅H₉ radicals, although the subsequent isomerization and decomposition reactions of the Ċ₅H₉ radicals are also of direct relevance to the radicals formed from the pyrolysis and oxidation of species including pentene and cyclopentane. Moreover, H-atom abstraction reactions by hydrogen atoms from 1,3-pentadiene are also investigated. The geometries and frequencies of 63 potential energy surface (PES) minima and 88 transition states are optimized at the ωB97XD/aug-cc-pVTZ level of theory. Spin-unrestricted open-shell single-point energies for all the species are calculated at the CCSD(T)/aug-cc-pVTZ level of theory with basis set corrections from MP2/aug-cc-pVXZ (where X = T and Q). A one-dimensional hindered rotor treatment is employed for torsional modes, with the M06-2X/6-311++G(D,P) method used to compute the potential energy as a function of the dihedral angle. The high-pressure limiting rate constants and the thermochemical properties for C₅ species are calculated using the Master Equation System Solver (MESS) with conventional transition-state theory and comparisons made with existing available literature data. A hydrogen atom can add to the terminal carbon atom of 1,3-pentadiene to form the 2,4-Ċ₅H₉ radical and/or the internal carbon atoms to form 2,5-Ċ₅H₉, 1,4-Ċ₅H₉, and 1,3-Ċ₅H₉ radicals. Among the four entrance channels for Ḣ atom addition reactions, the formation of 2,4-Ċ₅H₉ and 1,3-Ċ₅H₉ radicals is more exothermic in comparison to the other Ċ₅H₉ isomers (2,5-Ċ₅H₉, 1,4-Ċ₅H₉) because of the resonantly stabilized allylic structure. Consequently, the formation of the former is generally dominant in terms of barrier heights. Ḣ atom addition reactions to 1,3-pentadiene are compared to available C₃-C₅ alkenes and dienes, with external addition calculated to be kinetically favored over internal addition. However, the correlation between heats of formation and energy barriers for Ḣ atom addition to 1,2-dienes is different from that for 1,3- and 1,4-dienes. Hydrogen atom addition and abstraction rate constants are also compared for 1,3-pentadiene, with addition found to be dominant. The subsequent unimolecular reactions on the Ċ₅H₉ PES are found to be highly complex with reactions taking place on a multiple-well multiple-channel PES. For clarity, the reaction mechanism and kinetics of each Ċ₅H₉ radical are discussed individually in terms of the computed enthalpy of the reaction and activation, the transition-state structure/reaction class, and also in terms of the combustion species for which the reactions are of potential importance. The reactions on the Ċ₅H₉ PES are divided into three reaction classes (H-shift isomerization, cycloaddition, and β-scission reactions), and the reactivity-structure-based estimation rules for energy barriers are derived for these three reaction classes and compared to literature results for alkyl radicals.

Theoretical, Experimental, and Modeling Study of the Reaction of Hydrogen Atoms with 1- and 2-Pentene

Alkyl radicals are prominent in combustion chemistry as they are formed by hydrocarbon decomposition or from a radical attack on hydrocarbons. Accurate determinations of the thermochemistry and kinetics of their unimolecular isomerization and decomposition reactions and related addition reactions of alkenes are therefore important in simulating the combustion chemistry of virtually all hydrocarbon fuels. In this work, a comprehensive potential energy surface (PES) for Ḣ-atom addition to and abstraction from 1- and 2-pentene, and the subsequent C-C and C-H β-scission reactions, and H-atom transfer reactions has been considered. Thermochemical values for the species on the Ċ₅H₁₁ PES were calculated as a function of temperature (298-2000 K), with enthalpies of formation determined using a network of isodesmic reactions. High-pressure limiting and pressure-dependent rate constants were calculated using the Rice-Ramsperger-Kassel-Marcus theory coupled with a one-dimensional master equation. As a validation of our theoretical results, hydrogen atomic resonance absorption spectrometry experiments were performed on the Ḣ-atom addition and abstraction reactions of 1- and 2-pentene. By incorporating our calculations into a detailed chemical kinetic model (AramcoMech 3.0), excellent agreement with these experiments is observed. The theoretical results are further validated via a comprehensive series of simulations of literature data. Our a priori model is found to reproduce important absolute species concentrations and product ratios reported therein.

Influence of early stages of triglyceride pyrolysis on the formation of PAHs as coke precursors

Molecular beam (MB) time-of-flight mass spectrometry has been used to investigate thermal decomposition of triolein, to reveal the mechanisms of low temperature soot/coke formation characteristic for triglycerides (TGs). Mass detected pyrolysis products were observed at incremented temperatures using both VUV single photon ionization (general product detection) and REMPI based selective detection of aromatic products. To augment the simple mass characterizations, we have employed stoichiometric considerations; we have supplemented the analysis further by using the detailed information available from product analysis of batch reactor TG cracking. Both the VUV photoionization and batch reactor studies indicated that formation of C7-sized stable products is a marker of significant triolein decomposition that is coupled with PAH formation. A significant fraction of the C7 species observed likely formed as a result of a C-C bond scission at the allylic position to the ω-9 double bond of oleic acid. REMPI detection indicated a high specificity for PAH formation at three distinct molecular weight values, 276, 352 and 444 amu (the latter being a fullerene precursor). The stoichiometric analysis has shown that these PAHs likely arise from condensation reactions of either C7- or C8-sized fragments (three, four and five, respectively). The C8-sized intermediate would become essential whenever the PAH product of C7 fragment condensation contained an odd number of carbon atoms, resulting in a less stable aromatic structure with an incomplete double bond conjugation. MB experiments involving either addition or in situ generation of hydrogen resulted in an enhancement of lower molecular weight PAH formation, i.e., a decrease in the effective number of condensing fragments. In contrast, an increase in temperature yielded the opposite effect.

Chemical Kinetics of Hydrogen Atom Abstraction from Propargyl Sites by Hydrogen and Hydroxy Radicals

Hydrogen atom abstraction from propargyl C-H sites of alkynes plays a critical role in determining the reactivity of alkyne molecules and understanding the formation of soot precursors. This work reports a systematic theoretical study on the reaction mechanisms and rate constants for hydrogen abstraction reactions by hydrogen and hydroxy radicals from a series of alkyne molecules with different structural propargyl C-H atoms. Geometry optimizations and frequency calculations for all species are performed at M06-2X/cc-pVTZ level of theory and the hindered internal rotations are also treated at this level. The high-level W1BD and CCSD(T)/CBS theoretical calculations are used as a benchmark for a series of DFT calculations toward the selection of accurate DFT functionals for large reaction systems in this work. Based on the quantum chemistry calculations, rate constants are computed using the canonical transition state theory with tunneling correction and the treatment of internal rotations. The effects of the structure and reaction site on the energy barriers and rate constants are examined systematically. To the best of our knowledge, this work provides the first systematic study for one of the key initiation abstraction reactions for compounds containing propargyl hydrogen atoms.

Pressure-dependent rate rules for cycloaddition, intramolecular H-shift, and concerted elimination reactions of alkenyl peroxy radicals at low temperature

The reactions of cycloaddition, intramolecular H-shift and concerted elimination of alkenyl peroxy radicals are three kinds of important reactions in the low temperature combustion of alkenes. In this study, the cycloaddition reactions are divided into classes considering endo-cycloaddition, exo-cycloaddition and the size of the transition states; the intramolecular H-shift reactions are divided into classes depending upon the ring size of the transition states and the type of C-H bonds from which the hydrogen atom is transferred; the concerted elimination reactions are divided into classes according to the type of H-C_βC_αOO bond that is broken. All geometry optimizations are performed at the B3LYP/6-31G(2df,p) level. With the electronic structure calculations being performed using the composite Gaussian-4 (G4) method, high pressure limit rate constants and pressure-dependent rate constants at pressures varying from 0.01 to 100 atm are calculated by using canonical transition state theory and the Rice-Ramsberger-Kassel-Marcus/master equation method, respectively. All rate constants are given in the form of the modified Arrhenius expression. The high pressure limit rate rules and the pressure-dependent rate rules are derived by averaging the rate constants of a representative set of reactions in each class. The results show that the rate rules for these three classes of reactions have a large uncertainty and the impact of the pressure on the rate constants increases as temperature increases.

On-the-fly kinetics of hydrogen abstraction from polycyclic aromatic hydrocarbons by methyl/ethyl radicals

This work provides a rigorous procedure, within the framework of the Reaction Class Transition State Theory (RC-TST) and the Structure-Activity Relationship (SAR), for predicting reliable thermal rate constants on-the-fly for hydrogen abstraction reactions by methyl/ethyl radicals from Polycyclic Aromatic Hydrocarbons (PAHs) in a temperature range of 300-3000 K. All necessary RC-TST parameters were derived from ab initio calculations for a representative set of 36 reactions on which different error analyses and comparisons with available literature data were carried out. In addition to the good agreement between the RC-TST rate constants and the literature data, the detailed error analyses show that RC-TST/SAR, utilizing either the Linear Energy Relationship (LER) where only the reaction energy is needed or Barrier Height Grouping (BHG) where no additional data is needed, can predict the thermal rate constants for any reaction in the title reaction class with an average systematic error of less than 50% when compared to the explicit rate calculations. Therefore, the constructed RC-TST procedure can be confidently used to obtain reliable rate constants on the fly in an attempt to effectively construct detailed kinetic mechanisms for PAH-related fuels.

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

Ab initio derived group additivity model for intramolecular hydrogen abstraction reactions

A set of group additivity values for intramolecular hydrogen abstraction reactions of alkanes, alkenes and alkynes is reported. Calculating 448 reaction rate coefficients at the CBS-QB3 level of theory for 1-2 up to 1-7 hydrogen shift reactions allowed the estimation of ΔGAV° values for 270 groups. The influence of substituents on (1) the attacking radical, (2) the attacked carbon atom, and (3) the carbon chain between the attacking and attacked reactive atom has been systematically studied. Substituents have been varied between hydrogen atoms and sp3, sp2 and sp hybridized carbon atoms. It has been assumed that substituents further away from the reactive atoms or their connecting carbon chain have negligible influences on the kinetics. This group additivity model is applicable to a wide variety of reactions in the 300-1800 K temperature range. Correlations for tunneling coefficients have been generated which are complementary to the ΔGAV°'s to obtain accurate rate coefficients without the need for imaginary frequencies or electronic energies of activation. These correlations depend on the temperature and activation energy of the exothermic step. The group additivity model has been successfully applied to a test set of reactions also calculated at the CBS-QB3 level of theory. A mean absolute deviation of 1.18 to 1.71 has been achieved showing a good overall accuracy of the model.

Predicting pressure-dependent unimolecular rate constants using variational transition state theory with multidimensional tunneling combined with system-specific quantum RRK theory: a definitive test for fluoroform dissociation

Understanding the falloff in rate constants of gas-phase unimolecular reaction rate constants as the pressure is lowered is a fundamental problem in chemical kinetics, with practical importance for combustion, atmospheric chemistry, and essentially all gas-phase reaction mechanisms. In the present work, we use our recently developed system-specific quantum RRK theory, calibrated by canonical variational transition state theory with small-curvature tunneling, combined with the Lindemann-Hinshelwood mechanism, to model the dissociation reaction of fluoroform (CHF3), which provides a definitive test for falloff modeling. Our predicted pressure-dependent thermal rate constants are in excellent agreement with experimental values over a wide range of pressures and temperatures. The present validation of our methodology, which is able to include variational transition state effects, multidimensional tunneling based on the directly calculated potential energy surface along the tunneling path, and torsional and other vibrational anharmonicity, together with state-of-the-art reaction-path-based direct dynamics calculations, is important because the method is less empirical than models routinely used for generating full mechanisms, while also being simpler in key respects than full master equation treatments and the full reduced falloff curve and modified strong collision methods of Troe.

Nonmonotonic Temperature Dependence of the Pressure-Dependent Reaction Rate Constant and Kinetic Isotope Effect of Hydrogen Radical Reaction with Benzene Calculated by Variational Transition-State Theory

The reaction between H and benzene is a prototype for reactions of radicals with aromatic hydrocarbons. Here we report calculations of the reaction rate constants and the branching ratios of the two channels of the reaction (H addition and H abstraction) over a wide temperature and pressure range. Our calculations, obtained with an accurate potential energy surface, are based on variational transition-state theory for the high-pressure limit of the addition reaction and for the abstraction reaction and on system-specific quantum Rice-Ramsperger-Kassel theory calibrated by variational transition-state theory for pressure effects on the addition reaction. The latter is a very convenient way to include variational effects, corner-cutting tunneling, and anharmonicity in falloff calculations. Our results are in very good agreement with the limited experimental data and show the importance of including pressure effects in the temperature interval where the mechanism changes from addition to abstraction. We found a negative temperature effect of the total reaction rate constants at 1 atm pressure in the temperature region where experimental data are missing and accurate theoretical data were previously missing as well. We also calculated the H + C₆H₆/C₆D₆ and D + C₆H₆/C₆D₆ kinetic isotope effects, and we compared our H + C₆H₆ results to previous theoretical data for H + toluene. We report a very novel nonmonotonic dependence of the kinetic isotope effect on temperature. A particularly striking effect is the prediction of a negative temperature dependence of the total rate constant over 300-500 K wide temperature ranges, depending on the pressure but generally in the range from 600 to 1700 K, which includes the temperature range of ignition in gasoline engines, which is important because aromatics are important components of common fuels.

Intermolecular oxidative decarbonylative [2 + 2 + 2] carbocyclization of N-(2-ethynylaryl)acrylamides with tertiary and secondary alkyl aldehydes involving C(sp3)-H functionalization

A new metal-free oxidative decarbonylative [2 + 2 + 2] carbocyclization of N-(2-ethynylaryl)acrylamides with tertiary and secondary alkyl aldehydes is described. This reaction enables the formation of three new C-C bonds in a single reaction by a sequence of oxidative decarbonylation, radical addition across C-C unsaturated bonds, C-H functionalization and annulation, and represents the first oxidative decarbonylative [2 + 2 + 2] carbocyclization approach using tertiary and secondary alkyl aldehydes as a two carbon unit for assembling six-membered carbocycle-fused polycycles.

Toward the Development of a Fundamentally Based Chemical Model for Cyclopentanone: High-Pressure-Limit Rate Constants for H Atom Abstraction and Fuel Radical Decomposition

Theoretical aspects of the development of a chemical kinetic model for the pyrolysis and combustion of a cyclic ketone, cyclopentanone, are considered. Calculated thermodynamic and kinetic data are presented for the first time for the principal species including 2- and 3-oxo-cyclopentyl radicals, which are in reasonable agreement with the literature. These radicals can be formed via H atom abstraction reactions by Ḣ and Ö atoms and ȮH, HȮ2, and ĊH3 radicals, the rate constants of which have been calculated. Abstraction from the β-hydrogen atom is the dominant process when ȮH is involved, but the reverse holds true for HȮ2 radicals. The subsequent β-scission of the radicals formed is also determined, and it is shown that recent tunable VUV photoionization mass spectrometry experiments can be interpreted in this light. The bulk of the calculations used the composite model chemistry G4, which was benchmarked in the simplest case with a coupled cluster treatment, CCSD(T), in the complete basis set limit.

Performance of First-Principles-Based Reaction Class Transition State Theory

Performance of the Reaction Class Transition State Theory (RC-TST) for prediction of rates constants of elementary reactions is examined using data from its previous applications to a number of different reaction classes. The RC-TST theory is taking advantage of the common structure denominator of all reactions in a given family combined with structure activity relationships to provide a rigorous theoretical framework to obtain rate expression of any reaction within a reaction class in a simple and cost-effective manner. This opens the possibility for integrating this methodology with an automated mechanism generator for "on-the-fly" generation of accurate kinetic models of complex reacting systems.

Ab initio study of the influence of resonance stabilization on intramolecular ring closure reactions of hydrocarbon radicals

The cyclization reactions of dieneyl radicals provide a low energy route to the formation of molecular weight growth products.

Reactions of allylic radicals that impact molecular weight growth kinetics

The reactions of allylic radicals have the potential to play a critical role in molecular weight growth (MWG) kinetics during hydrocarbon oxidation and/or pyrolysis.