Ab initio thermochemistry and kinetics for carbon-centered radical addition and beta-scission reactions

A quantitative comparison of ab initio calculated rate coefficients using five computational methods and five different approaches of treating hindered internal rotation and tunneling with experimental values of rate coefficients for nine carbon-centered radical additions/beta scissions at 300, 600, and 1000 K is performed. The high-accuracy compound methods, CBS-QB3 and G3B3, and the density functionals, MPW1PW91, BB1K, and BMK, have been evaluated using the following approaches: (i) the harmonic oscillator approximation; (ii) the hindered internal rotor approximation for the internal rotation about the forming/breaking bond in the transition state and product; and the hindered internal rotation approximation combined with (iii) Wigner, (iv) Skodje and Truhlar, and (v) Eckart zero-curvature tunneling corrections. The density functional theory (DFT) based values for beta-scission rate coefficients deviate significantly from the experimental ones at 300 K, and the DFT methods do not accurately predict the equilibrium coefficient. The hindered rotor approximation offers a significant improvement in the agreement with experimental rate coefficients as compared to the harmonic oscillator treatment, especially at higher temperatures. Tunneling correction factors are smaller than 1.40 at 300 K and 1.03 at 1000 K. For both the CBS-QB3 method, including the hindered rotor treatment but excluding tunneling corrections, and the G3B3 method, including hindered rotor and Eckart tunneling corrections, a mean factor of deviation with experimentally observed values of 3 is found.     

Criegee中间体RCHOO(R=H,CH_3)与NCO的反应机理

采用CCSD(T)/aug-cc-p VTZ//B3LYP/6-311+G(2df,2p)方法对Criegee中间体RCHOO(R=H,CH_3)与NCO反应的机理进行了研究,利用经典过渡态理论(TST)并结合Eckart校正模型计算了标题反应在298~500 K范围内优势通道的速率常数.结果表明,上述反应包含亲核加成、氧化和抽氢3类机理,其中每类又包括NCO中N和O分别进攻的两种形式.亲核加成反应中O端进攻为优势通道,氧化和抽氢反应则是N端进攻为优势通道;甲基取代使CH_3CHOO反应活性高于CH2OO;anti-CH_3CHOO的加成及氧化反应活性高于syn-CH_3CHOO,而抽氢反应则是syn-CH_3CHOO的活性高于anti-CH_3CHOO.anti-构象对总速率常数的贡献大于syn-构象,且总速率常数具有显著的负温度效应.     

DFT and ab initio direct dynamics studies on the hydrogen abstraction reactions of chlorine atoms with CH(4-n)F(n) (n = 1-3)

A direct dynamics study is carried out for the hydrogen abstraction reactions Cl + CH(4-n)F(n) (n = 1-3) in the temperature range of 200-1,000 K. The minimum energy paths (MEPs) of these reactions are calculated at the BH&H-LYP/6-311G(d,p) level, and the energies along the MEPs are further refined at the QCISD(T)/6-311+G(2df,2p) and QCISD(T)/6-311+G(d,p) (single-point) level. The rate constants obtained by using the improved canonical variational transition state theory incorporating small-curvature tunneling correction (ICVT/SCT) are in good agreement with the available experimental results. It is shown that the vibrational adiabatic potential energy curves for these reactions have two barriers, a situation similar to the analogous reactions CH(3)X+Cl (X=Cl, Br). The theoretical results show that for the title reactions the variational effect should not be neglected over the whole considered temperature range, while the small-curvature tunneling effect is only important in the lower temperature range. The effects of fluorine substitution on the rate of this kind of reactions are also examined.     

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.     

Ab initio and kinetic calculations for the reactions of Cl with SiH(n)Cl(4-n) (n=1,2,3,4)

The direct hydrogen abstraction reactions of Cl atom with SiH(n)Cl(4-n) (n=1,2,3,4) have been studied systematically using ab initio molecular orbital theory. Geometries have been optimized at the MP2 level with 6-311+G(d) basis set, QCISD(T)/6-311+G(d,p) has been used in the final single point energy calculation. The kinetic calculations of these reactions have been explored using the canonical variational transition (CVT) state theory method with small-curvature tunneling (SCT) effect correction over the temperature range of 200-2000 K. The CVT/SCT rate constants exhibit typical non-Arrhenius behavior and three-parameter rate-temperature formulas have been fitted for the reactions of Cl with SiH4, SiH3Cl, SiH2Cl2, and SiHCl3, respectively (in unit of cm(3) molecule(-1) s(-1)). The calculated CVT/SCT rate constants are in agreement with the available experimental values.     

Reaction-path dynamics and theoretical rate constants for the CH(n)F(4-n) + O3 --> HOOO + CH(n-1)F(4-n) (n = 2,3) reactions

We present a theoretical study of the hydrogen abstraction reactions from CH(3)F and CH(2)F(2) by an ozone molecule. The geometries and harmonic vibrational frequencies of all stationary points are calculated at the MPW1K, BHandHLYP, and MPWB1K levels of theory. The energies of all of the stationary points were refined by using both higher-level (denoted as HL) energy calculations and QCISD(T)/6-311++G(2df,2pd) calculations based on the optimized geometries at the MPW1K/6-31+G(d,p) level of theory. The minimum energy paths (MEPs) were obtained by the MPW1K/6-31+G(d,p) level of theory. Energetic information of the points along the MEPs is further refined by the HL method. The rate constants were evaluated on the basis of the MEPs from the HL level of theory in the temperature range 200-2500 K by using the conventional transition-state theory (TST), the canonical variational transition-state theory (CVT), the microcanonical variational transition-state theory (microVT), the CVT coupled with the small-curvature tunneling (SCT) correction (CVT/SCT), and the microVT coupled with the Eckart tunneling correction (microVT/Eckart) based on the ab initio calculations. A general agreement was found among the TST, CVT, and microVT theories. The fitted three-parameter Arrhenius expressions of the calculated forward CVT/SCT and microVT/Eckart rate constants of the ozonolysis of fluoromethane are k(CVT/SCT)(T) = 2.76 x 10(-34)T(5.81)e((-13975/)(T)) and k(microVT/Eckart)(T) = 1.15 x 10(-34)T(5.97)e((-14530.7/)(T)), respectively. The fitted three-parameter Arrhenius expressions of the calculated forward CVT/SCT and microVT/Eckart rate constants of the ozonolysis of difluoromethane are k(CVT/SCT)(T) = 2.29 x 10(-36)T(6.42)e((-15451.6/)(T)) and k(microVT/Eckart)(T) = 1.31 x 10(-36)T(6.45)e((-15465.8/)(T)), respectively.     

Tunneling in hydrogen-transfer isomerization of n-alkyl radicals

The role of quantum tunneling in hydrogen shift in linear heptyl radicals is explored using multidimensional, small-curvature tunneling method for the transmission coefficients and a potential energy surface computed at the CBS-QB3 level of theory. Several one-dimensional approximations (Wigner, Skodje and Truhlar, and Eckart methods) were compared to the multidimensional results. The Eckart method was found to be sufficiently accurate in comparison to the small-curvature tunneling results for a wide range of temperature, but this agreement is in fact fortuitous and caused by error cancellations. High-pressure limit rate constants were calculated using the transition state theory with treatment of hindered rotations and Eckart transmission coefficients for all hydrogen-transfer isomerizations in n-pentyl to n-octyl radicals. Rate constants are found in good agreement with experimental kinetic data available for n-pentyl and n-hexyl radicals. In the case of n-heptyl and n-octyl, our calculated rates agree well with limited experimentally derived data. Several conclusions made in the experimental studies of Tsang et al. (Tsang, W.; McGivern, W. S.; Manion, J. A. Proc. Combust. Inst. 2009, 32, 131-138) are confirmed theoretically: older low-temperature experimental data, characterized by small pre-exponential factors and activation energies, can be reconciled with high-temperature data by taking into account tunneling; at low temperatures, transmission coefficients are substantially larger for H-atom transfers through a five-membered ring transition state than those with six-membered rings; channels with transition ring structures involving greater than 8 atoms can be neglected because of entropic effects that inhibit such transitions. The set of computational kinetic rates were used to derive a general rate rule that explicitly accounts for tunneling. The rate rule is shown to reproduce closely the theoretical rate constants.     

First principles study of the reaction of formic and acetic acids with hydroxyl radicals

The oxidation of formic and acetic acids with hydroxyl radicals was studied as a model for the oxidation of larger carboxylic acids using first principles calculations. For formic acid, the CBS-QB3 activation barriers of 14.1 and 12.4 kJ/mol for the acid and for the formyl channel, respectively, are within 3 kJ/mol of benchmark W1U values. Tunneling significantly enhances the rate coefficient for the acid channel and is responsible for the dominance of the acid channel at 298 K. At 298 K, tunneling correction factors of 339 and 2.0 were calculated for the acid and the formyl channel using the small-curvature tunneling method and the CBS-QB3 potential energy surface. The Wigner, Eckart, and zero-curvature tunneling methods severely underestimate the importance of tunneling for the acid channel. The resulting reaction rate coefficient of 0.98 x 10(5) m(3)/(mol x s) at 298 K is within a factor 2-3 of experimental values. For acetic acid, an activation barrier of 11.0 kJ/mol and a tunneling correction factor of 199 were calculated for the acid channel. Two mechanisms compete for hydrogen abstraction at the methyl group, with activation barriers of 11.9 and 12.5 kJ/mol and tunneling correction factors of 9.1 and 4.1 at 298 K. The resulting rate coefficient of 1.2 x 10(5) m(3)/(mol x s) at 298 K and branching ratio of 94% compare well with experimental data.     

Theoretical studies on dynamics and thermochemistry of the reactions CF(3)CHCl(2)+Cl-->CF(3)CCl(2)+HCl and CF(3)CHFCl+Cl-->CF(3)CFCl+HCl

A dual-level direct dynamics study has been carried out for the two hydrogen abstraction reactions CF(3)CHCl(2)+Cl and CF(3)CHFCl+Cl. The geometries and frequencies of the stationary points are optimized at the BHLYP/6-311G(d,p), B3LYP/6-311G(d,p), and MP2/6-31G(d) levels, respectively, with single-point calculations for energy at the BHLYP/6-311++G(3df,2p), G3(MP2), and QCISD(T)/6-311G(d,p) levels. The enthalpies of formation for the species CF(3)CHCl(2), CF(3)CHFCl, CF(3)CCl(2), and CF(3)CFCl are evaluated at higher levels. With the information of the potential energy surface at BHLYP/6-311++G(3df,2p)//6-311G(d,p) level, we employ canonical variational transition-state theory with small-curvature tunneling correction to calculate the rate constants. The agreement between theoretical and experimental rate constants is good in the measured temperature range 276-382 K. The effect of fluorine substitution on reactivity of the C-H bond is discussed.     

DFT and ab initio dual-level direct dynamics studies on the reactions of fluorine atom with HOCl and HOBr

The multichannel reactions (1) HOCl + F --> products and (2) HOBr + F --> products have been investigated using the dual-level direct dynamics method. The minimum energy paths (MEPs) are calculated at both the MPW1K/6-311G(d,p) and QCISD/6-311G(d,p) levels, then the single-point energies are further corrected at the QCISD(T)/6-311++G(3df,3pd) level of theory. There are hydrogen-bonded complexes with the energies less than those of the reactants or products located at the entrance or exit channel of both hydrogen abstraction reactions; while for the halogen abstraction channels only one complex exists at the reactant side in the bromine abstraction channel. The rate constants are evaluated by the improved canonical variational transition-state theory (ICVT). The agreement of the rate constants with available experimental values for two reactions at room temperature is good. Theoretical results indicate that for the reaction HOCl + F, hydrogen abstraction channel leading to the formation of HF + ClO will predominate the reaction over the whole temperature range, and the reaction of HOBr + F may proceed mainly through the bromine abstraction channel at the lower temperature while the contribution of hydrogen abstraction channel will become significant as the temperature increases. Because of lack of the kinetic data of these reactions, the present theoretical results are expected to be useful and reasonable to estimate the dynamical properties of these reactions over a wide temperature range where no experimental value is available.     

Kinetics of α hydrogen abstractions from thiols, sulfides and thiocarbonyl compounds

Hydrogen abstraction reactions involving organosulfur compounds play an important role in many industrial, biological and atmospheric processes. Despite their chemical relevance, little is known about their kinetics. In this work a group additivity model is developed that allows predicting the Arrhenius parameters for abstraction reactions of α hydrogen atoms from thiols, alkyl sulfides, alkyl disulfides and thiocarbonyl compounds by carbon-centered radicals at temperatures ranging from 300 to 1500 K. Rate coefficients for 102 hydrogen abstractions were obtained using conventional transition state theory within the high-pressure limit. Electronic barriers were calculated using the CBS-QB3 method and the rate coefficients were corrected for tunneling and hindered rotation about the transitional bond. Group additivity values for 46 groups are determined. To account for resonance and hyperconjugative stabilization in the transition state, 8 resonance corrections were fitted to a set of 32 reactions. The developed group additivity scheme was validated using a test set containing an additional 30 reactions. The group additivity scheme succeeds in reproducing the rate coefficients on average within a factor of 2.4 at 300 K and 1.4 at 1000 K. Mean absolute deviations of the Arrhenius parameters amount to, respectively, 2.5 kJ mol(-1) for E(a) and 0.13 for log A, both at 300 and 1000 K. This work hence illustrates that the recently developed group additivity methods for Arrhenius parameters extrapolate successfully to hetero-element containing compounds.     

A study of the hydrogen abstraction reactions of C2H radical with CH3CN, C2H5CN, and C3H7CN by dual-level generalized transition state theory

The hydrogen abstraction reactions C2H + CH3CN --> products (R1), C2H + CH3CH2CN --> products (R2), and C2H + CH3CH2CH2CN --> products (R3) have been investigated by dual-level generalized transition state theory. Optimized geometries and frequencies of all the stationary points and extra points along the minimum-energy path (MEP) are performed at the BH&H-LYP and MP2 methods with the 6-311G(d, p) basis set, and the energy profiles are further refined at the MC-QCISD level of theory. The rate constants are evaluated using canonical variational transition state theory (CVT) with a small-curvature tunneling correction (SCT) over a wide temperature range 104-2000 K. The calculated CVT/SCT rate constants are in good agreement with the available experimental values. Our calculations show that for reaction R2, the alpha-hydrogen abstraction channel and beta-hydrogen abstraction channel are competitive over the whole temperature range. For reaction R3, the gamma-hydrogen abstraction channel is preferred at lower temperatures, while the contribution of beta-hydrogen abstraction will become more significant with a temperature increase. The branching ratio to the alpha-hydrogen abstraction channel is found negligible over the whole temperature range.     

Thermochemical properties, DeltafH degrees (298), S degrees (298), and Cp degrees (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.     

Theoretical study and rate constant computation on the reaction HFCO + OH --> CFO + H2O

The potential energy surface, including the geometries and frequencies of the stationary points, of the reaction HFCO + OH is calculated using the MP2 method with 6-31+G(d,p) basis set, which shows that the direct hydrogen abstraction route is the most dominating channel with respect to addition and substitution channels. For the hydrogen abstraction reaction, the single-point energies are refined at the QCISD(T) method with 6-311++G(2df,2pd) basis set. The calculated standard reaction enthalpy and barrier height are -17.1 and 4.9 kcal mol(-1), respectively, at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The reaction rate constants within 250-2500 K are calculated by the improved canonical variational transition state theory (ICVT) with small-curvature tunneling (SCT) correction at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The fitted three-parameter formula is k = 2.875 x 10(-13) (T/1000)1.85 exp(-325.0/T) cm(3) molecule(-1) s(-1). The results indicate that the calculated ICVT/SCT rate constant is in agreement with the experimental data, and the tunneling effect in the lower temperature range plays an important role in computing the reaction rate constants.     

Theoretical study on the hydrogen abstraction reaction of CH3CH2F (HFC-161) with Cl atom

A direct dynamics method is employed to study the hydrogen abstraction reaction of CH₃CH₂F+Cl. Three distinct transition states are located, one for -H abstraction and two for β-H abstraction. The potential energy surface (PES) information is obtained at the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) and G2//MP2/6-311G(d,p) level. Based on the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) results, the rate constants of the three reaction channels are evaluated by using the canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions over the temperature range of 220–2800 K. The calculated results indicate that -H abstraction dominates the total reaction almost over the whole temperature range.     

Kinetics of hydrogen abstraction reactions of butene isomers by OH radical

The rate coefficients of H-abstraction reactions of butene isomers by the OH radical were determined by both canonical variational transition-state theory and transition-state theory, with potential energy surfaces calculated at the CCSD(T)/6-311++G(d,p)//BH&HLYP/6-311G(d,p) level and CCSD(T)/6-311++G(d,p)//BH&HLYP/cc-pVTZ level and quantum mechanical tunneling effect corrected by either the small-curvature tunneling method or the Eckart method. While 1-butene contains allylic, vinylic, and alkyl hydrogens that can be abstracted to form different butene radicals, results reveal that s-allylic H-abstraction channels have low and broad energy barriers, and they are the most dominant channels which can occur via direct and indirect H-abstraction channels. For the indirect H-abstraction s-allylic channel, the reaction can proceed via forming two van der Waals prereactive complexes with energies that are 2.7-2.8 kcal mol(-1) lower than that of the entrance channel at 0 K. Assuming that neither mixing nor crossover occurs between different reaction pathways, the overall rate coefficient was calculated by summing the rate coefficients of the s-allyic, methyl, and vinyl H-abstraction paths and found to agree well with the experimentally measured OH disappearance rate. Furthermore, the rate coefficients of p-allylic H abstraction of cis-2-butene, trans-2-butene, and isobutene by the OH radical were also determined at 300-1500 K, with results analyzed and compared with available experimental data.     

Mechanism and kinetics of the reaction of OH radicals with glyoxal and methylglyoxal: a quantum chemistry + CVT/SCT approach.

A theoretical study of the mechanism and kinetics of the OH hydrogen abstraction from glyoxal and methylglyoxal is presented. Optimum geometries, frequencies, and gradients have been computed at the BHandHLYP/6-311++G(d,p) level of theory for all the stationary points, as well as for 12 additional points along the minimum energy path (MEP). Energies were obtained by single-point calculations at the above geometries using CCSD(T)/ 6-311++G(d,p) to produce the potential energy surface. The rate coefficients were calculated for the temperature range 200-500 K by using canonical variational theory (CVT) with small-curvature tunneling (SCT) corrections. Our analysis suggests a stepwise mechanism, which involves the formation of a reactant complex. The overall agreement between the calculated and experimental kinetic data is very good. This agreement supports the reliability of the Arrhenius parameters of the glyoxal + OH reaction that are proposed in this work for the first time. The Arrhenius expressions that best describe the studied reactions are k1 = (9.63 +/- 0.23) x l0(-13)exp[(517 +/- 7)/T] and k2 = (3.93 +/- 0.11) x 10(-13)exp[(1060 +/- 8)/T]cm3 molecule(-1)s(-1) for glyoxal and methylglyoxal, respectively.     

Theoretical study on the hydrogen abstraction reactions from hydrazine derivatives by H atom

The hydrogen abstraction reactions from hydrazine and its methyl derivatives by the H atom have been investigated theoretically by using CBS-QB3//DSD-BLYP-D3(BJ)/Def2-TZVP quantum chemical calculations and transition state theory calculations coupled with various tunneling correction methods. Both the products and transition state energies of the hydrogen abstraction from the amino group were stabilized by the methyl group substitution. The substitution effect on the ^αN site was two times larger than that on the ^βN site. On the other hand, the substitution effect was negligible on the hydrogen abstraction from the methyl group. The overall rate coefficients of N₂H₄ + H reaction calculated by canonical variational transition state theory with the small-curvature tunneling correction agreed well with previously reported values, but those of CH₃NHNH₂/(CH₃)₂NNH₂ + H were slightly lower than a previous experimental value. The product-specific rate coefficients have been proposed for the kinetics modeling of these fuels’ combustion.     

Computational study of the reaction of CH2(X3B1) with CH3OH

The reaction of triplet methylene with methanol is a key process in alcohol combustion but surprisingly this reaction has never been studied. The reaction mechanism is investigated by using various high-level ab initio methods, including the complete basis set extrapolation (CBS-QB3 and CBS-APNO), the latest Gaussian-n composite method (G4), and the Weizmann-1 method (W1U). A total of five product channels and six transition states are found. The dominant mechanism is direct hydrogen abstraction, and the major product channel is CH(3) + CH(3)O, involving a weak prereactive complex and a 7.4 kcal/mol barrier. The other hydrogen abstraction channel, CH(3) + CH(2)OH, is less important even though it is more exothermic and involves a similar barrier height. The rate coefficients are predicted in the temperature range 200-3000 K. The tunneling effect and the hindered internal rotational freedoms play a key role in the reaction. Moreover, the reaction shows significant kinetic isotope effect.     

Direct Ab Initio Dynamics Study on the Reaction of CH3CHF2 (HFC‐152a) with the Cl Atom

A direct ab initio dynamics method is used to investigate the hydrogen‐abstraction reaction CH₃CHF₂+Cl. One transition state is located for α‐H abstraction, and two are identified for β‐H abstraction. The potential‐energy surface (PES) is obtained at the G3(MP2)//MP2/6‐311G(d, p) level. Furthermore, the rate constants of the three channels are evaluated by using canonical variational transition‐state theory (CVT) with small‐curvature tunneling (SCT) contributions over a wide temperature range of 200–2500 K. The dynamic calculations show that the reaction proceeds mainly by α‐H abstraction over the whole temperature range. The calculated rate constants and branching ratios are both in good agreement with the available experimental values.