An Ab initio determination of the rate constant for H2 + CN → H + HCN

The reactants, products, and saddle point for the reaction H₂ + CN → H + HCN have been studied by ab initio calculations. The computed structures, frequencies, and energetics are compared directly to available measurements and, indirectly, to experimental rateconstants. The theoretical rate constants used in the comparison are calculated with conventional transition state theory. By reduction of the computed reaction barrier to 4.1 kcal mol,^?1 good agreement with experimental rate constants is obtained over a 3250-K temperature range. This computed rate constant is well represented by the form 4.9 × 10_?18 T^2.45 e^{?1, 126/T} over the temperature range of 250 K–3500 K. Substantial reaction rate curvature is found due to low-frequency bending modes at the saddle point. The results for this reaction are compared to other abstraction reactions involving H atom transfer to identify correlations between reaction exothermicity and both abstraction barriers and reaction rate curvature.     

Direct ab initio dynamics study of the reaction of C2(A3Πu) with CH4

The reaction of C₂(A³Π_u) with CH₄ has been investigated over a wide temperature range 200–3,000 K by direct ab initio dynamics method at the BMC‐CCSD//BB1K/6‐311+G(2d,2p) level of theory. The optimized geometries and frequencies of the stationary points are calculated at the BB1K/6‐311+G(2d,2p) level, and then the energy profiles of the reactions are refined using the BMC‐CCSD method. The activation barrier height for H‐abstraction reaction was calculated to be 4.44 kcal/mol in temperature range (337–605 K), and the electron transfer behavior was also analyzed by quasi‐restricted molecular orbital method in detail. The canonical variational transition‐state theory (CVT) with the small curvature tunneling (SCT) correction method is used to calculate the rate constants over a wide temperature range 200–3,000 K. The theoretical results shows that variational effect is to some extent large in lower temperature range, and small curvature and tunneling effect play important roles to the H‐atom abstraction only at lower temperatures. The CVT/SCT rate constants are in good agreement with the available experimental results. Our theoretical study is expected to provide a direct insight into the reaction mechanism and may be useful for estimating the kinetics of the title reaction over a wide temperature range where no experimental data are available so far. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Kinetic isotope effects for the H2 + C2H ↔ C2H2 + H reaction based on the ab initio calculations and a global potential energy surface

In the present paper, kinetic isotope effects of the title reaction are studied with canonical variational transition state theory on the modified Wang Bowman (MWB) potential energy surface (PES) (Chem Phys Lett 2005, 409, 249) and the ab initio calculations at the quadratic configuration interaction (QCISD (T, full))/aug‐cc‐pVTZ//QCISD (full)/cc‐pVTZ level. The calculated rate constants for the isotopic variants of this title reaction on the MWB PES have good agreement with those of the present ab initio calculations over the temperature range of 20–5000 K for the forward reactions and 800–5000 K for the reverse reactions, respectively. In particular, the forward rate constants for the title reaction and its isotopically substituted reactions have negative temperature dependences at about 40 K. Rate expressions are presented for all the studied reactions. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 42: 289–298, 2010     

Reaction rates for O + HD → OH + D and O + HD → OD + H

We have calculated reaction rates for the reactions O + HD → OH + D and O + DH → OD + H using improved canonical variational transition state theory and least-action ground-state transmission coefficients with an ab initio potential energy surface. The kinetic isotope effects are in good agreement with experiment. The optimized tunneling paths and properties of the variational transition states and the rate enhancement for vibrationally excited reactants are also presented and compared with those for the isotopically unsubstituted reaction O + H₂ → OH + H. The thermal reactions at low and room temperature are predicted to occur by tunneling at extended configurations, i.e., to initiate early on the reaction path and to avoid the saddle point regions. Tunneling also dominates the low and room temperature reactions for excited vibrational states, but in these cases the results are not as sensitive to the nature of the tunneling path. Overbarrier mechanisms dominate for both thermal and excited-vibrational state reactions for T > 600 K. For the excited-state reaction (with initial vibrational quantum number n > 0) a transition state switch occurs for T > 1000 K for the O + HD(n = 1) → OD + H case and for T > 1500 K for the O + DH(n = 1) → OD + H reaction, and this may be a general phenomenon for excited-state reactions at higher temperature. In the present case the switch occurs from an early variational transition state where the vibrationally adiabatic approximation is expected to be valid to a tighter variational transition state where nonadiabatic effects are probably important and should be included.     

Thermal and state-selected rate constant calculations for O(3p) + H2 → OH + H and isotopic analogs

We present a new parametrization (based on ab initio calculations) of the bending potentials for the two lowest potential energy surfaces of the reaction O(³P) + H₂, and we use it for rate constant calculations by variational transition-state theory with multidimensional semiclassical tunneling corrections. We present results for the temperature range 250–2400 K for both the rate constants and the intermolecular kinetic isotope effects for the reactions of O(³P) with D₂ and HD. In general, the calculated rate constants for the thermal reactions are in excellent agreement with available experiments. We also calculate the enhancement effect for exciting H₂ to the first excited vibrational state. The calculations also provide information on which aspects of the potential energy surfaces are important for determining the predicted rate constants.     

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.     

Quantifying the non-RRKM effect in the H + O2 ⇄ OH + O reaction

In the present investigation the non-RRKM behavior in the title reaction is quantified in two different ways: (1) Quasiclassical trajectory calculations of the thermal rate coefficient are compared with results from a microcanonical variational transition-state theory/RRKM model. Results on both the Varandas DMBE IV and Melius-Blint potentials indicate that the non-RRKM behavior acts to reduce the thermal rate coefficient by about a factor of two, independent of temperature from 250 K to 5500 K. The QCT thermal rate coefficients on the two potentials are in remarkably good agreement with experiment and with each other over the entire temperature range. (2) The non-RRKM behavior as a classical phenomenon is demonstrated and quantified on both potentials by a direct test of the fundamental assumption. Complex-forming classical trajectories, started as either O + OH or H + O₂, are shown preferentially to return to the region of configuration space from which they were started. This test is discussed in detail in the text. The transition of the non-RRKM behavior from classical to quantum mechanics is also discussed. © 1997 John Wiley & Sons, Inc. Int J Chem Kinet 29: 275–287, 1997.     

Direct dynamics study of the hydrogen abstraction reaction of CF3CH2Cl + Cl → CF3CHCl + HCl

The hydrogen abstraction reaction of Cl atoms with CF₃CH₂Cl (HCFC‐133a) is investigated by using density function theory and ab initio approach, and the rate constants are calculated by using the dual‐level direct dynamics method. Optimized geometries and frequencies of reactants, transition state, and products are computed at the B3LYP/6‐311+G(2d,2p) level. To refine the energetic information along the minimum energy path, single‐point energy calculations are carried out at the G3(MP2) level of theory. The interpolated single‐point energy method is employed to correct the energy profiles for the title reaction. The rate constants are evaluated by using the canonical variational transition state theory with a small‐curvature tunneling correction over a wide range of temperature, 200–2000 K. The variational effect for the reaction is moderate at low temperatures and very small at high temperatures. However, the tunneling correction has an important contribution in the lower temperature range. The agreement between calculated rate constants and available experimental values is good at lower temperatures but diverges significantly at higher temperatures. © 2012 Wiley Periodicals, Inc. Int J Chem Kinet 44: 661–667, 2012     

Ab initio study on the kinetics of hydrogen abstraction for the H+alkene-->H2+alkenyl reaction class

Kinetics of the hydrogen abstraction reaction class of the H+alkene has been studied using the reaction class transition state theory (RC-TST) combined with the linear energy relationship (LER) and the barrier height grouping (BHG) approach. The rate constants for the reference reaction, H+C2H4, were obtained by the canonical variational transition state theory (CVT) with the small curvature tunneling (SCT) correction in the temperature range of 300-3000 K. Combined with these data, both the RC-TST/LER, where only reaction energy is needed, and RC-TST/BHG, where no other information is needed, are found to be promising methods for predicting rate constants for a large number of reactions in this reaction class. Our analysis indicates that less than 50% systematic errors on the average exist in the predicted rate constants using the RC-TST/LER or RC-TST/BHG method while in comparison to explicit rate calculations the differences are less than 100% or a factor of 2 on the average.     

Canonical variational transition-state theory study of the CF3CH2CH3 + OH reaction

Variational transition-state theory rate constants with multidimensional tunneling contributions using the small curvature method have been calculated for the CF3CH2CH3 (HFC-263fb) + OH reaction over a temperature range from 200 to 373 K. The mPW1B95-41.0 hybrid functional, parametrized by Albu and Swaminathan to generate theoretical rate constants nearly identical to the experimental values for the CH3F + OH reaction, has been used in conjunction with the 6-31+G** basis set to explore the potential energy surface of the title reaction. The good agreement found between theoretical predictions and the experimental data available suggests that the present approach is an excellent option to obtain high-quality results at low computational cost for direct dynamics studies of hydrogen abstraction reactions from complex hydrofluorocarbons. The reliability of the structure activity relationship used to estimate rate constant values for OH reactions with hydrofluorocarbons is also discussed in detail.     

Theoretical study on the rate constants for the C2H5 + HBr --> C2H6 + Br reaction

The reaction C(2)H(5) + HBr --> C(2)H(6) + Br has been theoretically studied over the temperature range from 200 to 1400 K. The electronic structure information is calculated at the BHLYP/6-311+G(d,p) and QCISD/6-31+G(d) levels. With the aid of intrinsic reaction coordinate theory, the minimum energy paths (MEPs) are obtained at the both levels, and the energies along the MEP are further refined by performing the single-point calculations at the PMP4(SDTQ)/6-311+G(3df,2p)//BHLYP and QCISD(T)/6-311++G(2df,2pd)//QCISD levels. The calculated ICVT/SCT rate constants are in good agreement with available experimental values, and the calculate results further indicate that the C(2)H(5) + HBr reaction has negative temperature dependence at T < 850 K, but clearly shows positive temperature dependence at T > 850 K. The current work predicts that the kinetic isotope effect for the title reaction is inverse in the temperature range from 200 to 482 K, i.e., k(HBr)/k(DBr) < 1.     

Comprehensive theoretical studies on the CF3H dissociation mechanism and the reactions of CF3H with OH and H free radicals

The potential energy surfaces for the CF3H unimolecular dissociation reaction and reactions of CF3H with free radical OH and H were investigated at the B3LYP6-311++G(**) and QCISD(T)6-311++G(**) levels and by the G3B3 theory. All the possible stationary and first-order saddle points along the reaction paths were verified by the vibrational analysis. The calculations account for all the product channels. The reaction enthalpies obtained at the G3B3 level are in good agreement with the available experiments. Canonical transition-state theory with Wigner tunneling correction was used to predict the rate constants for the temperature range of 298-2500 K without any artificial adjustment, and tshe computed rate constants for elementary channels can be accurately fitted with three-parameter Arrhenius expressions. The theoretical rate constants of the CF3H+H reaction agree with the available experimental data very well. The theoretical and experimental rate constants for the CF3H+OH reaction are in reasonable agreement. The H abstraction of CF3H by OH is found to be the main reaction channel for the CF3H fire extinguishing reactions while the CF3H unimolecular dissociation reaction plays a negligible role.     

Isotope branching and tunneling in O(3P)+HD-->OH+D; OD+H reactions

The O((3)P)+HD and O((3)P)+D(2) reactions are studied using quantum scattering calculations and chemically accurate potential energy surfaces developed for the O((3)P)+H(2) system by Rogers et al. [J. Phys. Chem. A 104, 2308 (2000)]. Cross sections and rate coefficients for OH and OD products are calculated using accurate quantum methods as well as the J-shifting approximation. The J-shifting approach is found to work remarkably well for both O+HD and O+D(2) collisions. The reactions are dominated by tunneling at low temperatures and for the O+HD reaction the hydrogen atom transfer leading to the OH product dominates at low temperatures. Our result for the OH/OD branching ratio is in close agreement with previous calculations over a wide range of temperatures. The computed OH/OD branching ratios are also in close agreement with experimental results of Robie et al. [Chem. Phys. Lett. 134, 579 (1987)] at temperatures above 400 K but the theoretical results do not reproduce the rapid rise in the experimental values of the branching ratio for temperatures lower than 350 K. We believe that new measurements could resolve the long-standing discrepancy between experiment and theory for this benchmark reaction.     

H3PO→H2POH异构化反应的直接动力学研究

在QCISD(T)/6-311C++G(2df,2pd)//QCISD/6-311C++G(d,p)+ZPE水平上,对H3PO的异构化反应H3PO→(1)H2POH(trans)→(2)H2POH(cis)进行了计算研究.结果表明,H原子由P原子向O原子迁移反应(1)的能垒为250.0kJ/mol,是反应速率控制步骤,而O_H键绕P_O键旋转的构型转化反应(2)的能垒只为12.3kJ/mol.利用经典过渡态理论(TST)与变分过渡态理论(CVT)分别计算了反应(1)在200～2000K温度区间内的速率常数kTST和kCVT,获得了经小曲率隧道效应(SCT)及Eckart模型校正后的速率常数kTST/Eckart和kCVT/SCT.对只涉及H原子迁移的反应(1),量子力学隧道效应的影响在低温段非常明显,而变分效应对反应速率常数的影响很小.     

Ab initio investigation on the reaction path and rate for the gas-phase reaction of HO + H2O <--> H2O + OH

This article describes an ab initio investigation on the potential surfaces for one of the simplest hydrogen atom abstraction reactions, that is, HO + H2O <--> H2O + OH. In accord with the findings in the previously reported theoretical investigations, two types of the hydrogen-bonding complexes [HOH--OH] and [H2O--HO] were located on the potential energy surface. The water molecule acts as a hydrogen donor in the [HOH--OH] complex, while the OH radical acts as a hydrogen donor in the [H2O--HO] complex. The energy evaluations at the MP2(FC) basis set limit, as well as those through the CBS-APNO procedure, have provided estimates for enthalpies of association for these complexes at 298 K as -2.1 approximately -2.3 and -4.1 approximately -4.3 kcal/mol, respectively. The IRC calculations have suggested that the [H2O--HO] complex should be located along the reaction coordinate for the hydrogen abstraction. Our best estimate for the classical barrier height for the hydrogen abstraction is 7.8 kcal/mol, which was obtained from the CBS-APNO energy evaluations. After fitting the CBS-APNO potential energy curve to a symmetrical Eckart function, the rate constants were calculated by using the transition state theory including the tunneling correction. Our estimates for the Arrhenius parameters in the temperature region from 300 to 420 K show quite reasonable agreement with the experimentally derived values.     

H原子和(CH3)2SiH2抽提反应的理论研究

采用UMP2/6-31G(d)理论水平优化了H原子和(CH3)2SiH2抽提反应势能面上的所有驻点,并在此水平基础上进行了内禀反应坐标(IRC)的计算,得到该反应的反应途径(MEP)。应用变分过渡态理论及最小能量途径半经典绝热基态隧道效应校正(MEPSAG)、小曲率半经典绝热基态隧道效应校正(SCSAG)等方法对上述反应进行了动力学研究,期望从理论上提供一套温度范围较宽、精度较高的动力学数据,为阐明反应机理和解释实验结果提供理论依据。     

Direct dynamic study on the hydrogen abstraction reaction C2(Πu)+H2 → C2H+H

The ab initio direct dynamics method at the G2//UQCISD/6-311 + G(d,p) level is employed to study the hydrogen abstraction reaction C₂(³Π_u)+H₂ → C₂H+H over a wide temperature range 100–4650 K. The barrier heights obtained for the forward and reverse reactions are 7.78 and 17.53 kcal/mol, respectively. Comparing with one recent experiment, the calculated forward rate constants over the temperature range 2580–4650 K are about 4.4–13.5 times greater and show a steeper temperature-dependent effect. This indicates that further experimental investigation on this simple radical reaction may still be desired. Finally, G2//UQCISD/6-311 + G(2df,2p) calculations are performed to test the reliability of the G2//UQCISD/6-311 + G(d,p) results.     

Variational transition‐state theory study of the atmospheric reaction OH + O3 → HO2 + O2

We report variational transition‐state theory calculations for the OH + O₃→ HO₂ + O₂ reaction based on the recently reported double many‐body expansion potential energy surface for ground‐state HO₄ [Chem Phys Lett 2000, 331, 474]. The barrier height of 1.884 kcal mol^?1 is comparable to the value of 1.77–2.0 kcal mol^?1 suggested by experimental measurements, both much smaller than the value of 2.16–5.11 kcal mol^?1 predicted by previous ab initio calculations. The calculated rate constant shows good agreement with available experimental results and a previous theoretical dynamics prediction, thus implying that the previous ab initio calculations will significantly underestimate the rate constant. Variational and tunneling effects are found to be negligible over the temperature range 100–2000 K. The O₁? O₂ bond is shown to be spectator like during the reactive process, which confirms a previous theoretical dynamics prediction. © 2007 Wiley Periodicals, Inc. 39: 148–153, 2007     

Reflected shock tube studies of high-temperature rate constants for OH + CH4 --> CH3 + H2O and CH3 + NO2 --> CH3O + NO

The reflected shock tube technique with multipass absorption spectrometric detection of OH radicals at 308 nm has been used to study the reactions OH + CH(4) --> CH(3) + H(2)O and CH(3) + NO(2) --> CH(3)O + NO. Over the temperature range 840-2025 K, the rate constants for the first reaction can be represented by the Arrhenius expression k = (9.52 +/- 1.62) x 10(-11) exp[(-4134 +/- 222 K)/T] cm(3) molecule(-1) s(-1). Since this reaction is important in both combustion and atmospheric chemistry, there have been many prior investigations with a variety of techniques. The present results extend the temperature range by 500 K and have been combined with the most accurate earlier studies to derive an evaluation over the extended temperature range 195-2025 K. A three-parameter expression describes the rate behavior over this temperature range, k = (1.66 x 10(-18))T(2.182) exp[(-1231 K)/T] cm(3) molecule(-1) s(-1). Previous theoretical studies are discussed, and the present evaluation is compared to earlier theoretical estimates. Since CH(3) radicals are a product of the reaction and could cause secondary perturbations in rate constant determinations, the second reaction was studied by OH radical production from the fast reactions CH(3)O --> CH(2)O + H and H + NO(2) --> OH + NO. The measured rate constant is 2.26 x 10(-11) cm(3) molecule(-1) s(-1) and is not dependent on temperature from 233 to 1700 K within experimental error.     

Full dimensional time-dependent quantum dynamics study of the H + NH3 --> H2 + NH2 reaction

A rigorous full dimensional time-dependent wave packet method has been developed for the reactive scattering between an atom and a tetra-atomic molecule. The method has been applied to the hydrogen abstraction reaction H+NH(3)-->H(2)+NH(2). Initial state-selected total reaction probabilities are investigated for the reactions from the ground vibrational state and from four excited vibrational states of ammonia. The total reaction probabilities from two lowest "tunneling doublets" due to the inversion barrier for the umbrella bending motion of NH(3) and from two pairs of doubly degenerate vibrational states of NH(3) are also inspected. Integral cross sections and rate constants are calculated for the reaction from the ground state with the centrifugal-sudden approximation. The calculated results are compared with those from the previous seven dimensional calculations [M. Yang and J. C. Corchado, J. Chem. Phys. 126, 214312 (2007)]. This work shows that the full dimensional rate constants are a factor of 3 larger than the corresponding seven dimensional calculated values at T=200 K and are overall smaller than those obtained from the variational transition state theory in the whole temperature region. The work also reveals that nonreactive NH bonds of NH(3) cannot be treated as spectators due to the fact that three NH bonds are coupled with each other during the reaction process.