Quantum chemical study of hydrogen abstraction reactions of the ethynyl radical with hydrogen compounds (C2H+HX)

We have theoretically investigated the hydrogen abstraction reactions of ethynyl radical with simple hydrogen compounds, C₂H+HX, using quantum chemical computations. Computations have been performed using the density functional theory with the recently proposed MPW1K functional and the 6-311++G(3df,2p) basis set. An analysis of the resulting energy barriers for hydrogen abstraction reactions has been carried out using the bond dissociation energy of the breaking X–H bond and DFT-based reactivity parameters to rationalize the reaction behavior.     

Ab-initio study of initial atmospheric oxidation reactions of C3 and C4 alkanes

Second-order, Møller–Plesset (MP2)-unrestricted Hartree–Fock calculations with full geometry optimization in the 6-31G(d, p) basis set were carried out to study the initial atmospheric oxidation reactions of alkanes. All structures in the initial hydrogen abstraction reaction by an OH radical and the subsequent addition of molecular oxygen to the alkyl radical were characterized for alkanes with three and four carbon atoms. The reaction paths for the formation of the peroxyl radicals were obtained and discussed in the light of similarities along series involving primary, secondary, and tertiary hydrogens. A 0.999 correlation was found between the height of our barriers for the OH abstraction of a primary hydrogen atom from alkanes containing one to four carbon atoms and the optimally estimated activation energies for this reaction recently presented. From the slope and the intersection at zero activation energy an equation was obtained that yields scaled values of the activation energies to account for the tunnel effect and for the error due to the basis set and the method employed. We present new results for the abstraction of the less favored primary hydrogens in propane, butane, and isobutane, which should be important at high temperatures. Negative net activation energies were obtained for the addition of molecular oxygen to all the alkyl radicals formed in the first reaction. The structure of the peroxyl radicals is discussed; and very good correlations are observed for similar compounds, regardless of the length of the carbon chain. A revision of some experimental values is suggested. Single point density functional calculations at the MP2 geometries were also performed with the B3LYP functional for comparison. The observed trends are exactly the same for the two methods. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 845–856, 1999     

Valence bond calculations of hydrogen transfer reactions: a general predictive pattern derived from theory

Hydrogen abstraction reactions of the type X(*) + H-H' --> X-H + H'(*) (X = F, Cl, Br, I) are studied by ab initio valence bond methods and the VB state correlation diagram (VBSCD) model. The reaction barriers and VB parameters of the VBSCD are computed by using the breathing orbital valence bond and valence bond configuration interaction methods. The combination of the VBSCD model and semiempirical VB theory leads to analytical expressions for the barriers and other VB quantities that match the ab initio VB calculations fairly well. The barriers are influenced by the endo- or exothermicity of the reaction, but the fundamental factor of the barrier is the average singlet-triplet gap of the bonds that are broken or formed in the reactions. Some further approximations lead to a simple formula that expresses the barrier for nonidentity and identity hydrogen abstraction reactions as a function of the bond strengths of reactants and products. The semiempirical expressions are shown to be useful not only for the model reactions that are studied in this work, but also for other nonidentity and identity hydrogen abstraction reactions that have been studied in previous articles.     

烷烃与氢过氧自由基氢提取反应类反应能垒与速率常数的精确计算

氢过氧自由基从烷烃分子中提取氢的反应是碳氢燃料中低温燃烧化学中非常重要的一类反应。本文用等键反应方法计算了这一类反应的动力学参数。所有反应物、过渡态、产物的几何结构均在HF/6-31+G(d)水平下优化得到。以反应中的过渡态反应中心的几何结构守恒为判据,该反应类可用等键反应处理。本文选取了乙烷和氢过氧自由基的氢提取反应为参考反应,其它反应作为目标反应,用等键反应方法对目标反应在HF/6-31+G(d)水平的近似能垒和反应速率常数进行了校正。为了验证方法的可靠性,选取C5以下的烷烃分子体系,对等键反应方法校正结果和高精度CCSD(T)/CBS直接计算结果进行了比较,最大绝对误差为5.58k J?mol~(-1),因此,采用等键反应方法只需用低水平HF从头算方法就可以再现高精度CCSD(T)/CBS计算结果,从而解决了该反应类中大分子体系的能垒的精确计算。本文的研究为碳氢化合物中低温燃烧模拟中重要的烷烃与氢过氧自由基氢提取反应提供了准确的动力学参数。     

Theoretical Studies on the Reaction Mechanism of 1-Chloroethane with Hydroxyl Radical

The reaction mechanism of 1-chloroethane with hydroxyl radical has been inves- tigated by using density functional theory (DFT) B3LYP/6-31G (d, p) method. All bond dissociation enthalpies were computed at the same theoretical level. It was found that hydrogen abstraction pathway is the most favorable. There are two hydrogen abstraction pathways with activation barriers of 0.630 and 4.988 kJ/mol, respectively, while chlorine abstraction pathway was not found. It was observed that activation energies have a more reasonable correlation with the reaction enthalpy changes (△Hr) than with bond dissociation enthalpies (BDE).     

The influence of hydrogen bonding interactions on the C-H bond activation step in class I ribonucleotide reductases

In order to model the C-H bond activation step in ribonucleotide reductases the hydrogen atom abstraction reaction from cis-tetrahydrofuran-2,3-diol (7) by methylthiyl (8) radical has been studied with theoretical methods. In order to identify an appropriate theoretical method for this system, the hydrogen transfer reaction between radical 8 and methanol (9) to give methanol radical (10) and methyl thiol (11) has been studied at several different levels of theory. While the reaction energy for this process is predicted equally well by the Becke3LYP and BH and HLYP hybrid functional methods, the reaction barrier is predicted to be significantly lower by the former. Compared to results obtained at CCSD(T)/cc-pVTZ level the BH and HLYP functional is better suited for the calculation of activation barriers for hydrogen abstraction reactions. This latter method was subsequently used to study the reaction of radical 8 with cis-tetrahydrofuran-2,3-diol 7 in the absence and in the presence of additional functional groups (acetate and acetamide) as models for the substrate reaction of class I ribonucleotide reductases (RNRs). The reaction barrier is lowest in those systems, in which acetate forms a double hydrogen bonded complex with the hydroxy groups of diol 7 (+8.2 kcal mol-1) and increases somewhat for side-on complexes between substrate 7 and acetate featuring only one hydrogen bond (+10.5 kcal mol-1). The barrier reduction of 6.5 kcal mol-1 obtained through complexation of diol 7 with acetate appears to be due to the formation of short strong hydrogen bonds in the transition. These effects can also be found in reactions of thiyl radical 8 with complexes of diol 7 with acetamide, but to a much smaller extent. The lowest reaction barrier is in this case calculated for the side-on complex (+11.2 kcal mol-1), while the bridging orientation between diol 7 and acetamide leads to a reaction barrier (+13.4) that is only slightly lower than that for the uncatalyzed process (+14.7 kcal mol-1). With respect to the structure of the active site of the RNR R1 subunit, only the side-on complexes appear to be relevant for the enzyme-catalyzed process. Under this condition the influence of the E441 side chan and thus the impact of the E441Q mutation in the initial C-H bond activation step will be rather small.     

Kinetics of the hydrogen abstraction reactions of 1,1‐ and 1,2‐difluoroethane with hydroxyl radical: an ab initio study

The hydrogen abstraction reactions of 1,1‐ and 1,2‐difluoroethane with the OH radical have been investigated by the ab initio molecular orbital theory. The geometries of the reactants, products, and transition states have been optimized at the (U)MP2=full level of theory in conjunction with 6‐311G(d,p) basis functions. Single‐point (U)MP2=full with larger basis set, such as 6‐311G(3d,2p), and QCISD(T)=full/6‐311G(d,p) calculations have also been carried out to observe the effects of basis sets utilized and higher order electron correlation. Three and four reaction channels have been identified for 1,1‐ and 1,2‐difluoroethane, respectively. In the case of 1,1‐difluoroethane, hydrogen abstraction from the α‐carbon has been found to be easier than that from the β‐carbon. The barriers of the four reaction channels for 1,2‐difluoroethane are close to each other. Weak hydrogen bonding interactions have been observed between hydroxyl hydrogen and a fluorine atom in the transition states. Rate constants for the reactions of 1,1‐ and 1,2‐difluoroethane with the OH radical have been calculated using the standard transition state theory and found to be in good agreement with the experimental results. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 1305–1318, 2000     

Mechanistic and kinetic study of the O + CH2OH reaction

The mechanism for the O + CH2OH reaction was investigated by various ab initio quantum chemistry methods. For the chemical activation mechanism, that is, the addition/elimination path, the couple-cluster methods including CCSD and CCSD(T) were employed with the cc-pVXZ (X = D, T, Q, 5) basis sets. For the abstraction channels, multireference methods including CASSCF, CASPT2, and MRCISD were used with the cc-pVDZ and cc-pVTZ basis sets. It has been shown that the production of H + HCOOH is the major channel in the chemical activation mechanism. The minor channels include HCO + H2O and OH + CH2O. The hydrogen abstraction by an O atom from the CH2OH radical produces either OH + CH2O or OH + HCOH. Moreover, the two abstraction reactions are essentially barrierless processes. The rate constants for the association of O with CH2OH have been calculated using the flexible transition state theory. A weak negative temperature dependence of the rate constants is found in the range 250-1000 K. Furthermore, it is estimated that the abstraction processes also play an important role in the O + CH2OH reaction. Additionally, the falloff behavior for the OCH2OH --> H + HCOOH reaction has been investigated. The present theoretical results are compared to the experimental measurements to understand the mechanism and kinetic behavior of the O + CH2OH reaction and the unimolecular reaction of the OCH2OH radical.     

A systematic computational study of the reactions of HO2 with RO2: the HO2 + C2H5O2 reaction

The reaction of HO2 with C2H5O2 has been studied using the density functional theory (B3LYP) and the coupled-cluster theory [CCSD(T)]. The reaction proceeds on the triplet potential energy surface via hydrogen abstraction to form ethyl hydroperoxide and oxygen. On the singlet potential energy surface, the addition-elimination mechanism is revealed. Variational transition state theory is used to calculate the temperature-dependent rate constants in the range 200-1000 K. At low temperatures (e.g., below 300 K), the reaction takes place predominantly on the triplet surface. The calculated low-temperature rate constants are in good agreement with the experimental data. As the temperature increases, the singlet reaction mechanism plays more and more important role, with the formation of OH radical predominantly. The isotope effect of the reaction (DO2 + C2D5O2 vs HO2 + C2H5O2) is negligible. In addition, the triplet abstraction energetic routes for the reactions of HO2 with 11 alkylperoxy radicals (CnHmO2) are studied. It is shown that the room-temperature rate constants have good linear correlation with the activation energies for the hydrogen abstraction.     

Thermally activated mechanisms of hydrogen abstraction by growth precursors during plasma deposition of silicon thin films

Hydrogen abstraction by growth precursors is the dominant process responsible for reducing the hydrogen content of amorphous silicon thin films grown from SiH(4) discharges at low temperatures. Besides direct (Eley-Rideal) abstraction, gas-phase radicals may first adsorb on the growth surface and abstract hydrogen in a subsequent process, giving rise to thermally activated precursor-mediated (PM) and Langmuir-Hinshelwood (LH) abstraction mechanisms. Using results of first-principles density functional theory (DFT) calculations on the interaction of SiH(3) radicals with the hydrogen-terminated Si(001)-(2x1) surface, we show that precursor-mediated abstraction mechanisms can be described by a chemisorbed SiH(3) radical hopping between overcoordinated surface Si atoms while being weakly bonded to the surface before encountering a favorable site for hydrogen abstraction. The calculated energy barrier of 0.39 eV for the PM abstraction reaction is commensurate with the calculated barrier of 0.43-0.47 eV for diffusion of SiH(3) on the hydrogen-terminated Si(001)-(2x1) surface, which allows the radical to sample the entire surface for hydrogen atoms to abstract. In addition, using the same type of DFT analysis we have found that LH reaction pathways involve bond breaking between the silicon atoms of the chemisorbed SiH(3) radical and the film prior to hydrogen abstraction. The LH reaction pathways exhibit energy barriers of 0.76 eV or higher, confining the abstraction only to nearest-neighbor hydrogens. Furthermore, we have found that LH processes compete with radical desorption from the hydrogen-terminated Si(001)-(2x1) surface and may be suppressed by the dissociation of chemisorbed SiH(3) radicals into lower surface hydrides. Analysis of molecular-dynamics simulations of the growth process of plasma deposited silicon films have revealed that qualitatively similar pathways for thermally activated hydrogen abstraction also occur commonly on the amorphous silicon growth surface.     

氮原子与卤代甲烷反应的直接氢抽提过程研究

利用B3LYP理论研究了N(~4S)+CH_3X(X = H, F, Cl)反应体系的直接氢抽提过 程,分别得到了各反应物、产物和过渡态的优化构型和谐振频率。同时应用了6- 31G(d), 6-311+G(d,p)和6-311+ + G(2d,2p)基组,考察其大小对反应体系中各物 种构型及能量的影响。理论计算表明,随着基组的增加,反应势垒逐渐降低,反应 吸热减少。对比取代甲烷的情形,结果表明反应过程中卤素原子具有典型的诱导效 应,降低了抽提势垒。     

Efficient protocol for quantum Monte Carlo calculations of hydrogen abstraction barriers: Application to methanol

Accurate calculation of hydrogen abstraction reaction barriers is a challenging problem, often requiring high level quantum chemistry methods that scale poorly with system size. Quantum Monte Carlo (QMC) methods provide an alternative approach that exhibit much better scaling, but these methods are still computationally expensive. We describe approaches that can significantly reduce the cost of QMC calculations of barrier heights, using the hydrogen abstraction of methanol by a hydrogen atom as an illustrative example. By analysing the combined influence of trial wavefunctions and pseudopotential quadrature settings on the barrier heights, variance, and time‐step errors, we devise a simple protocol that minimizes the cost of the QMC calculations while retaining accuracy comparable to large‐basis coupled cluster theory. We demonstrate that this protocol is transferable to other hydrogen abstraction reactions.     

反应类等键反应方法及类反应热动力学参数计算

本文提出了反应类等键反应方法,将通常用于热力学性质计算的等键反应方法推广用于类反应中反应势垒和反应焓变的计算. 对碳氢燃料低温燃烧反应机理中的一类重要反应类：?烷基自由基过氧化氢裂解生成烯烃和HO2自由基的反应势垒和反应焓变进行了计算. 通过对该类16个反应中的5个代表反应分别在不同计算水平HF、DFT、MP2、CCSD(T)的比较计算发现,采用等键反应方法可在较低从头算级别计算得到类反应较高精度的反应势垒,提高了计算的效率和精度. 本文采用反应类等键反应方法在B3LYP/6-311G(d,p)计算水平对该类16个反应进行了反应势垒和反应焓变的计算,并建立了反应势垒和反应能的线性自由能关系：delta V=71.02+0.41?delta E.     

An assessment of theoretical procedures for predicting the thermochemistry and kinetics of hydrogen abstraction by methyl radical from benzene

The reaction enthalpy (298 K), barrier (0 K), and activation energy and preexponential factor (600-800 K) have been examined computationally for the abstraction of hydrogen from benzene by the methyl radical, to assess their sensitivity to the applied level of theory. The computational methods considered include high-level composite procedures, including W1, G3-RAD, G3(MP2)-RAD, and CBS-QB3, as well as conventional ab initio and density functional theory (DFT) methods, with the latter two classes employing the 6-31G(d), 6-31+G(d,p) and/or 6-311+G(3df,2p) basis sets, and including ZPVE/thermal corrections obtained from 6-31G(d) or 6-31+G(d,p) calculations. Virtually all the theoretical procedures except UMP2 are found to give geometries that are suitable for subsequent calculation of the reaction enthalpy and barrier. For the reaction enthalpy, W1, G3-RAD, and URCCSD(T) give best agreement with experiment, while the large-basis-set DFT procedures slightly underestimate the endothermicity. The reaction barrier is slightly more sensitive to the choice of basis set and/or correlation level, with URCCSD(T) and the low-cost BMK method providing values in close agreement with the benchmark G3-RAD value. Inspection of the theoretically calculated rate parameters reveals a minor dependence on the level of theory for the preexponential factor. There is more sensitivity for the activation energy, with a reasonable agreement with experiment being obtained for the G3 methods and the hybrid functionals BMK, BB1K, and MPW1K, especially in combination with the 6-311+G(3df,2p) basis set. Overall, the high-level G3-RAD composite procedure, URCCSD(T), and the cost-effective DFT methods BMK, BB1K, and MPW1K give the best results among the methods assessed for calculating the thermochemistry and kinetics of hydrogen abstraction by the methyl radical from benzene.     

Skewered reactions in free‐radical crosslinking (co)polymerizations of liquid polybutadiene as internal olefinic multivinyl crosslinker

As an extension of our continuing studies concerned with the mechanistic discussion of network formation in the free‐radical crosslinking (co)polymerization of multivinyl monomers, this work refers to the skewered reactions in the crosslinking (co)polymerizations of liquid polybutadiene rubber (LBR) as an internal olefinic multivinyl monomer or crosslinker, especially focused on the competitive occurrence of both addition or skewered reaction to internal carbon–carbon (CC) double bonds and abstraction reaction of allylic hydrogens in LBR by growing polymer radical. Thus, LBR is regarded as an internal olefinic multiallyl monomer‐linked allyl groups (? CH?CH? CH₂? ) with methylene units (? CH₂? ). First, gelation in the polymerization of LBR was explored in detail, especially at elevated temperatures. The occurrence of intermolecular crosslinking was easier in the order LBR > LBR containing 20 mol % of 1,2‐structural units > liquid polyisoprene rubber. Then, we pursued the polymerization of LBR using dicumyl peroxide (DCPO) as typical organic peroxide used at elevated temperatures. The primary cumyloxy radical generated by the thermal decomposition of DCPO may add to CC double bond or abstract allylic hydrogen or undergo β‐scission to generate a secondary methyl radical. The initiation by the cumyloxy radical was omitted. The ratio of allylic hydrogen abstraction to β‐scission reaction was estimated; thus, only 39% of cumyloxy radical was used for the allylic hydrogen abstraction reaction. The addition of methyl radical to CC double bond was clearly observed. Finally, we pursued the intermolecular and intramolecular skewered reactions in free‐radical crosslinking LBR/vinyl pivalate copolymerizations. © 2011 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem, 2012     

Exploring the H-abstraction reactions of CHCl·-/CCl2·- with CX3H（X = F,Cl,Br and I） using the density functional theory method

Gas-phase hydrogen abstraction reactions have been compared using the popular density functional theory(DFT) functional BHandHLYP/aug-cc-pVTZ/RECP level of theory,on the basis of the model reaction CHCl·-/CCl2·-+ CX3H(X = F,Cl,Br and I).Our theoretical findings suggest the efficiency of the H-abstraction reactions induced by either CHCl·-or CCl2·-increases as the substrate is changed from CF3H to CI3H,and that CHCl·-has a higher activity in hydrogen abstraction than CCl2·-for a given substrate.The entropy effect at 298 K does not significantly change the trend in reactivity of the various reactions,which is in general controlled by the heights of activation energies △E≠.Therefore,we have explored the origin of the energy barriers △E≠ of the reactions using the activation strain model of chemical reactivity.     

The benzene+OH potential energy surface: intermediates and transition states

The potential energy surface for the interaction between benzene and hydroxyl radical is studied in detail using quantum mechanical methods, with a particular focus on the hydrogen abstraction pathway. Geometric parameters are optimized using a variety of density functional methods as well as perturbation theory. Energies are refined using coupled cluster singles and doubles with perturbative triples [CCSD(T)] extrapolated to the complete basis set limit. At our most reliable level of theory, complexation energies are found to be (with zero-point corrected energies in parentheses) 3.7 (2.8) kcal/mol for the benzene-hydroxyl radical complex and 2.9 (-1.7) kcal/mol for the phenyl radical-water complex. The barrier to H abstraction lies 6.5 (4.2) kcal/mol above the infinitely separated benzene and hydroxyl radical monomers.     

Reactivity of alkanes on zeolites: a computational study of propane conversion reactions

In this work, quantum chemical methods were used to study propane conversion reactions on zeolites; these reactions included protolytic cracking, primary hydrogen exchange, secondary hydrogen exchange, and dehydrogenation reactions. The reactants, products, and transition-state structures were optimized at the B3LYP/6-31G level and the energies were calculated with CBS-QB3, a complete basis set composite energy method. The computed activation barriers were 62.1 and 62.6 kcal/mol for protolytic cracking through two different transition states, 30.4 kcal/mol for primary hydrogen exchange, 29.8 kcal/mol for secondary hydrogen exchange, and 76.7 kcal/mol for dehydrogenation reactions. The effects of basis set for the geometry optimization and zeolite acidity on the reaction barriers were also investigated. Adding extra polarization and diffuse functions for the geometry optimization did not affect the activation barriers obtained with the composite energy method. The largest difference in calculated activation barriers is within 1 kcal/mol. Reaction activation barriers do change as zeolite acidity changes, however. Linear relationships were found between activation barriers and zeolite deprotonation energies. Analytical expressions for each reaction were proposed so that accurate activation barriers can be obtained when using different zeolites as catalysts, as long as the deprotonation energies are first acquired.     

A DFT‐Based Investigation of Hydrogen Abstraction Reactions from Methylated Polycyclic Aromatic Hydrocarbons

The growth of polycyclic aromatic hydrocarbons (PAHs) is in many areas of combustion and pyrolysis of hydrocarbons an inconvenient side effect that warrants an extensive investigation of the underlying reaction mechanism, which is known to be a cascade of radical reactions. Herein, the focus lies on one of the key reaction classes within the coke formation process: hydrogen abstraction reactions induced by a methyl radical from methylated benzenoid species. It has been shown previously that hydrogen abstractions determine the global PAH formation rate. In particular, the influence of the polyaromatic environment on the thermodynamic and kinetic properties is the subject of a thorough exploration. Reaction enthalpies at 298 K, reaction barriers at 0 K, rate constants, and kinetic parameters (within the temperature interval 700–1100 K) are calculated by using B3LYP/6‐31+G(d,p) geometries and BMK/6‐311+G(3df,2p) single‐point energies. This level of theory has been validated with available experimental data for the abstraction at toluene. The enhanced stability of the product benzylic radicals and its influence on the reaction enthalpies is highlighted. Corrections for tunneling effects and hindered (or free) rotations of the methyl group are taken into account. The largest spreading in thermochemical and kinetic data is observed in the series of linear acenes, and a normal reactivity–enthalpy relationship is obtained. The abstraction of a methyl hydrogen atom at one of the center rings of large methylated acenes is largely preferred. Geometrical and electronic aspects lie at the basis of this striking feature. Comparison with hydrogen abstractions leading to arylic radicals is also made.