Analysis of Tertiary Butyl Radical + O2, Isobutene + HO2, Isobutene + OH, and Isobutene−OH Adducts + O2: A Detailed Tertiary Butyl Oxidation Mechanism

The multichannel i-propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol-1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2') lie above the reactants by (1.4, 2.5) kcal mol-1 and display large Born-Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol-1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol-1 of the corresponding lowest-energy minima.

Mechanisms and Kinetics Studies of Butylated Hydroxytoluene Degradation to Isobutene

2,6-Di-tert-butyl-hydroxytotulene (BHT) is a widely used antioxidant in various fields. In this study, we explored comprehensively the mechanisms and kinetics of BHT degradation to produce isobutene using the density functional theory method. Furthermore, the intrinsic chemical reactivity of BHT was investigated using the electrostatic potential, average local ionization energy, and Fukui function, and the most likely reaction site with OH radical was predicted. Two initiation pathways of BHT with OH radicals were reported. The OH addition pathways at the C2 site of BHT was found more likely to occur than the pathways of H abstracts from the t-butyl group due to the lower energy barrier. Rate constants of two initiation pathways were calculated by transition state theory, and they were promoted by the temperature rise. Mayer bond order and localized molecular orbitals analysis were conducted to reveal the variation of the chemical bonds in the reaction process. The tertiary butyl radical that had been generated in the OH-addition reaction was more likely to generate isobutene with the participation of oxygen. Overall, this research could help to reveal the transformation mechanism of isobutene produced by BHT degradation.

Infrared spectroscopic signature of a hydroperoxyalkyl radical (•QOOH)

Infrared (IR) action spectroscopy is utilized to characterize a prototypical carbon-centered hydroperoxyalkyl radical (•QOOH) transiently formed in the oxidation of volatile organic compounds. The •QOOH radical formed in isobutane oxidation, 2-hydroperoxy-2-methylprop-1-yl, •CH₂(CH₃)₂COOH, is generated in the laboratory by H-atom abstraction from tert-butyl hydroperoxide (TBHP). IR spectral features of jet-cooled and stabilized •QOOH radicals are observed from 2950 to 7050 cm^-1 at energies that lie below and above the transition state barrier leading to OH radical and cyclic ether products. The observed •QOOH features include overtone OH and CH stretch transitions, combination bands involving OH or CH stretch and a lower frequency mode, and fundamental OH and CH stretch transitions. Most features arise from a single vibrational transition with band contours well simulated at a rotational temperature of 10 K. In each case, the OH products resulting from unimolecular decay of vibrationally activated •QOOH are detected by UV laser-induced fluorescence. Assignments of observed •QOOH IR transitions are guided by anharmonic frequencies computed using second order vibrational perturbation theory, a 2 + 1 model that focuses on the coupling of the OH stretch with two low-frequency torsions, as well as recently predicted statistical •QOOH unimolecular decay rates that include heavy-atom tunneling. Most of the observed vibrational transitions of •QOOH are readily distinguished from those of the TBHP precursor. The distinctive IR transitions of •QOOH, including the strong fundamental OH stretch, provide a general means for detection of •QOOH under controlled laboratory and real-world conditions.

Accurate many-body expansion potential energy surface for SiH2 (11 A′) using a switching function formalism

An accurate many-body expansion potential energy surface for the ground state of SiH2 is reported. To warrant the correct behavior at the Si (1D) + H2 (X1Σ+g) dissociation channels...

Watching a hydroperoxyalkyl radical (•QOOH) dissociate

A prototypical hydroperoxyalkyl radical (•QOOH) intermediate, transiently formed in the oxidation of volatile organic compounds, was directly observed through its infrared fingerprint and energy-dependent unimolecular decay to hydroxyl radical and cyclic ether products. Direct time-domain measurements of •QOOH unimolecular dissociation rates over a wide range of energies were found to be in accord with those predicted theoretically using state-of-the-art electronic structure characterizations of the transition state barrier region. Unimolecular decay was enhanced by substantial heavy-atom tunneling involving O-O elongation and C-C-O angle contraction along the reaction pathway. Master equation modeling yielded a fully a priori prediction of the pressure-dependent thermal unimolecular dissociation rates for the •QOOH intermediate-again increased by heavy-atom tunneling-which are required for global models of atmospheric and combustion chemistry.

Thermochemical unification of molecular descriptors to predict radical hydrogen abstraction with low computational cost

Chemistry describes transformation of matter with reaction equations and corresponding rate constants. However, accurate rate constants are not always easy to get. Here we focus on radical oxidation reactions. Analysis of over 500 published rate constants of hydroxyl radicals led us to hypothesize that a modified linear free-energy relationship (LFER) could be used to predict rate constants speedily, reliably and accurately. LFERs correlate the Gibbs activation-energy with the Gibbs energy of reaction. We calculated the latter as the sum of one-electron transfer and, if appropriate, proton transfer. We parametrized specific transition state effects to orbital delocalizability and the polarity of the reactant. The calculation time for 500 reactions is less than 8 hours on a standard desktop-PC. Rate constants were also calculated for hydrogen and methyl radicals; these controls show that the predictions are applicable to a broader set of oxidizing radicals. An accuracy of 30-40% (standard deviation) with reference to reported experimental values was found suitable for the screening of complex chemical systems for possibly relevant reactions. In particular, potentially relevant reactions can be singled out and scrutinized in detail when prioritizing chemicals for environmental risk assessment.

A Systematic Theoretical Kinetics Analysis for the Waddington Mechanism in the Low-Temperature Oxidation of Butene and Butanol Isomers

The Waddington mechanism, or the Waddington-type reaction pathway, is crucial for low-temperature oxidation of both alkenes and alcohols. In this study, the Waddington mechanism in the oxidation chemistry of butene and butanol isomers was systematically investigated. Fundamental quantum chemical calculations were conducted for the rate constants and thermodynamic properties of the reactions and species in this mechanism. Calculations were performed using two different ab initio solvers: Gaussian 09 and Orca 4.0.0, and two different kinetic solvers: PAPR and MultiWell, comprehensively. Temperature- and pressure-dependent rate constants were performed based on the transition state theory, associated with the Rice Ramsperger Kassel Marcus and master equation theories. Temperature-dependent thermochemistry (enthalpies of formation, entropy, and heat capacity) of all major species was also conducted, based on the statistical thermodynamics. Of the two types of reaction, dissociation reactions were significantly faster than isomerization reactions, while the rate constants of both reactions converged toward higher temperatures. In comparison, between two ab initio solvers, the barrier height difference among all isomerization and dissociation reactions was about 2 and 0.5 kcal/mol, respectively, resulting in less than 50%, and a factor of 2–10 differences for the predicted rate coefficients of the two reaction types, respectively. Comparing the two kinetic solvers, the rate constants of the isomerization reactions showed less than a 32% difference, while the rate of one dissociation reaction (P1 ↔ WDT12) exhibited 1–2 orders of magnitude discrepancy. Compared with results from the literature, both reaction rate coefficients (R4 and R5 reaction systems) and species’ thermochemistry (all closed shell molecules and open shell radicals R4 and R5) showed good agreement with the corresponding values obtained from the literature. All calculated results can be directly used for the chemical kinetic model development of butene and butanol isomer oxidation.

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.

TOF MS Investigation of Nickel Oxide CVD

NiO layers were deposited by metal-organic chemical vapor deposition using bis-(ethylcyclopentadienyl) nickel (EtCp)₂Ni and oxygen or ozone. As a continuation of kinetic study of NiO MOCVD the gas-phase, transformations of (EtCp)₂Ni were studied in the temperature range of 380-830 K. Time of reactions corresponding to the residence time of the gas stream in hot zone of the reactor was about 0.1 s under conditions studied. The interaction of (EtCp)₂Ni with oxygen started at 450 K and its conversion rate reached the maximum at 700 K. The interaction of (EtCp)₂Ni with ozone started at 400 K and its conversion rate reached the maximum at 600 K. Transformations of the gas phase with the temperature in the reaction zone were studied, the model reaction schemes illustrating (EtCp)₂Ni transformations in the reaction systems containing oxygen and ozone have developed. In the reaction system (EtCp)₂Ni-O₂-Ar the main gas-phase products at 380-500 K were CO, CO₂, HCO, C₂H₅OH, CpCOOH, and CpO. Formation of the C₂H₂O, C₃H₄O, and C₅H₈O was found at 630-830 K. The same gas-phase species, (C₄H₃O)₂Ni and dialdehydes was formed in the reaction system (EtCp)₂Ni-O₃-O₂-Ar. Graphical Abstract ᅟ.

Low-temperature combustion chemistry of biofuels: pathways in the initial low-temperature (550 K-750 K) oxidation chemistry of isopentanol

The branched C(5) alcohol isopentanol (3-methylbutan-1-ol) has shown promise as a potential biofuel both because of new advanced biochemical routes for its production and because of its combustion characteristics, in particular as a fuel for homogeneous-charge compression ignition (HCCI) or related strategies. In the present work, the fundamental autoignition chemistry of isopentanol is investigated by using the technique of pulsed-photolytic Cl-initiated oxidation and by analyzing the reacting mixture by time-resolved tunable synchrotron photoionization mass spectrometry in low-pressure (8 Torr) experiments in the 550-750 K temperature range. The mass-spectrometric experiments reveal a rich chemistry for the initial steps of isopentanol oxidation and give new insight into the low-temperature oxidation mechanism of medium-chain alcohols. Formation of isopentanal (3-methylbutanal) and unsaturated alcohols (including enols) associated with HO(2) production was observed. Cyclic ether channels are not observed, although such channels dominate OH formation in alkane oxidation. Rather, products are observed that correspond to formation of OH viaβ-C-C bond fission pathways of QOOH species derived from β- and γ-hydroxyisopentylperoxy (RO(2)) radicals. In these pathways, internal hydrogen abstraction in the RO(2)⇄ QOOH isomerization reaction takes place from either the -OH group or the C-H bond in α-position to the -OH group. These pathways should be broadly characteristic for longer-chain alcohol oxidation. Isomer-resolved branching ratios are deduced, showing evolution of the main products from 550 to 750 K, which can be qualitatively explained by the dominance of RO(2) chemistry at lower temperature and hydroxyisopentyl decomposition at higher temperature.

Innate C-H trifluoromethylation of heterocycles

Direct methods for the trifluoromethylation of heteroaromatic systems are in extremely high demand in nearly every sector of chemical industry. Here we report the discovery of a general procedure using a benchtop stable trifluoromethyl radical source that functions broadly on a variety of electron deficient and rich heteroaromatic systems and demonstrates high functional group tolerance. This C-H trifluoromethylation protocol is operationally simple (avoids gaseous CF(3)I), scalable, proceeds at ambient temperature, can be used directly on unprotected molecules, and is demonstrated to proceed at the innately reactive positions of the substrate. The unique and orthogonal reactivity of the trifluoromethyl radical relative to aryl radicals has also been investigated on both a complex natural product and a pharmaceutical agent. Finally, preliminary data suggest that the regioselectivity of C-H trifluoromethylation can be fine-tuned simply by judicious solvent choice.

Enthalpies of formation and bond dissociation energies of lower alkyl hydroperoxides and related hydroperoxy and alkoxy radicals

The enthalpies of formation and bond dissociation energies, D(ROO-H), D(RO-OH), D(RO-O), D(R-O 2) and D(R-OOH) of alkyl hydroperoxides, ROOH, alkyl peroxy, RO, and alkoxide radicals, RO, have been computed at CBS-QB3 and APNO levels of theory via isodesmic and atomization procedures for R = methyl, ethyl, n-propyl and isopropyl and n-butyl, tert-butyl, isobutyl and sec-butyl. We show that D(ROO-H) approximately 357, D(RO-OH) approximately 190 and D(RO-O) approximately 263 kJ mol (-1) for all R, whereas both D(R-OO) and D(R-OOH) strengthen with increasing methyl substitution at the alpha-carbon but remain constant with increasing carbon chain length. We recommend a new set of group additivity contributions for the estimation of enthalpies of formation and bond energies.

Contra-thermodynamic Behavior in Intermolecular Hydrogen Transfer of Alkylperoxy Radicals

Quantum chemical investigation of bimolecular hydrogen transfer involving alkylperoxy radicals, a key reaction family in the free-radical oxidation of hydrocarbons, was performed to establish structure-reactivity relationships. Eight different reactions were investigated featuring four different alkane substrates (methane, ethane, propane and isobutane) and two different alkylperoxy radicals (methylperoxy and iso-propylperoxy). Including forward and reverse pairs, sixteen different activation energies and enthalpies of reaction were used to formulate structure-reactivity relationships to describe this chemistry. We observed that the enthalpy of formation of loosely bound intermediate states has a strong inverse correlation with the overall heat of reaction and that this results in unique contra-thermodynamic behavior such that more exothermic reactions have higher activation barriers. A new structure-reactivity relationship was proposed that fits the calculated data extremely well: E(A)=E(o)+alphaDeltaH(rxn) where alpha=-0.10 for DeltaH(rxn)<0, and alpha=1.10 for DeltaH(rxn)>0 and E(o)=3.05 kcal mol(-1).

Thermochemical Properties, ΔfH°(298), S°(298), and Cp°(T), for n-Butyl and n-Pentyl Hydroperoxides and the Alkyl and Peroxy Radicals, Transition States, and Kinetics for Intramolecular Hydrogen Shift Reactions of the Peroxy Radicals

Alkyl radicals in atmospheric and combustion environments undergo a rapid association with molecular oxygen (3O2) to form an alkyl peroxy radical (ROO*). One important reaction of these peroxy radicals is the intramolecular H-shift (intramolecular abstraction) to form a hydroperoxide alkyl radical (R'*COOH), where the hydroperoxide alkyl radical may undergo chemical activation reaction with O2 and result in chain branching at moderate to low temperatures. The thermochemistry and trends in kinetic parameters for the hydrogen shift reactions from each carbon (4-8-member-ring TST's) in n-butyl and n-pentyl peroxy radicals (CCCCOO* and CCCCCOO*) are analyzed using density functional and ab initio calculation methods. Thermochemical properties, DeltafH degrees (298 K), C-H bond energies, S degrees (298 K), and Cp degrees (T) of saturated linear C4 and C5 aliphatic peroxides (ROOH), as well as the corresponding hydroperoxide alkyl radicals (R'*COOH), are determined. DeltafH degrees (298 K) are obtained from isodesmic reactions and the total energies of the CBS-QB3 and B3LYP computational methods. Contributions to the entropy and the heat capacity from translation, vibration, and external rotation are calculated using the rigid-rotor-harmonic-oscillator approximation based on the CBS-QB3 frequencies and structures. The results indicate that pre-exponential factors, A(T), decrease with the increase of the ring size (4-8-member-ring TS, H-atom included). The DeltaH for 4-, 5-, 6-, and 7-member rings in n-butyl (and n-pentyl) peroxy are 40.8 (40.8), 31.4 (31.5), 20.5 (20.0), 22.6-p (19.4) kcal mol(-1), respectively. The DeltaH for the 8-member ring in n-pentylperoxy is 23.8-p kcal mol(-1), All abstractions are from secondary (-CH2-) groups except those marked (-p), which are from primary sites. Enthalpy and barrier values from the B3LYP/6-311++G(2d,p) and BHandHLYP/6-311G(d,p) methods are compared with CBS-QB3 results. The B3LYP results show good agreement with the higher level CBS-QB3 calculation method; the BHandH barriers for the intramolecular peroxy H-shifts are not acceptable.

Thermochemical and Kinetic Analysis on the Reactions of O2 with Products from OH Addition to Isobutene, 2-Hydroxy-1,1-dimethylethyl, and 2-Hydroxy-2-methylpropyl Radicals: HO2 Formation from Oxidation of Neopentane, Part II

Unimolecular dissociation of a neopentyl radical to isobutene and methyl radical is competitive with the neopentyl association with O2 ((3)Sigma(g)-) in thermal oxidative systems. Furthermore, both isobutene and the OH radical are important primary products from the reactions of neopentyl with O2. Consequently, the reactions of O2 with the 2-hydroxy-1,1-dimethylethyl and 2-hydroxy-2-methylpropyl radicals resulting from the OH addition to isobutene are important to understanding the oxidation of neopentane and other branched hydrocarbons. Reactions that correspond to the association of radical adducts with O2((3)Sigma(g)-) involve chemically activated peroxy intermediates, which can isomerize and react to form one of several products before stabilization. The above reaction systems were analyzed with ab initio and density functional calculations to evaluate the thermochemistry, reaction paths, and kinetics that are important in neopentyl radical oxidation. The stationary points of potential energy surfaces were analyzed based on the enthalpies calculated at the CBS-Q level. The entropies, S(degrees)298, and heat capacities, C(p)(T), (0 <or= T/K <or= 1500), from vibration, translation, and external rotation contributions were calculated using statistical mechanics based on the vibrational frequencies and structures obtained from the density functional study. The hindered internal rotor contributions to S(degrees)298 and C(p)(T) were calculated by solving the Schrödinger equation with free rotor wave functions, and the partition coefficients were treated by direct integration over energy levels of the internal rotation potentials. Enthalpies of formation (DeltaH(f)(degrees)298) were determined using isodesmic reaction analysis. The DeltaH(f)(degrees)298 values of (CH3)2C*CH(2)OH, (CH3)2C(OO*)CH(2)OH, (CH3)2C(OH)C*H2, and (CH3)2C(OH)CH(2)OO* radicals were determined to be -23.3, -62.2, -24.2, and -61.8 kcal mol(-1), respectively. Elementary rate constants were calculated from canonical transition state theory, and pressure-dependent rate constants for multichannel reaction systems were calculated as functions of pressure and temperature using multifrequency quantum Rice-Ramsperger-Kassel (QRRK) analysis for k(E) and a master equation for pressure falloff. Kinetic parameters for intermediate and product formation channels of the above reaction systems are presented as functions of temperature and pressure.