Kinetics study of the OH + alkene --> H2O + alkenyl reaction class

In this paper we report on the kinetics of hydrogen abstraction for the OH + alkene reaction class, using the reaction class transition state theory (RC-TST) combined with the linear energy relationship (LER) and the barrier height grouping (BHG) approaches. Parameters for the RC-TST were derived from theoretical calculations using a set of 15 reactions representing the hydrogen abstractions from the terminal and nonterminal carbon sites of the double bond of alkene compounds. Both the RC-TST/LER, where only reaction energy is needed at either density functional theory BH&HLYP or semiempirical AM1 levels, and RC-TST/BHG, where no additional information is required, are found to be promising methods for predicting rate constants for a large number of reactions in this reaction class. Detailed error analyses show that, when compared to explicit theoretical calculations, the averaged systematic errors in the calculated rate constants using both the RC-TST/LER and RC-TST/BHG methods are less than 25% in the temperature range 300-3000 K. The estimated rate constants using these approaches are in good agreement with available data in the literature.     

Kinetics of the C-C bond beta scission reactions in alkyl radical reaction class

Kinetics of the β-scission in alkyl radical reaction class was studied using the reaction class transition state theory (RC-TST) combined with the linear energy relationship (LER) and the barrier height grouping (BHG) approach. All necessary parameters were derived from first-principle density functional calculations for a representative set of 21 reactions. Different error analyses and comparisons with available literature data were made. Direct comparison with available experimental data indicates that the RC-TST/LER, where only reaction energy is needed, can predict rate constants for any reaction in this reaction class with excellent accuracy. Specifically for this reaction class, the RC-TST/LER method has less than 60% systematic errors on average in the predicted rate constants when compared to explicit rate calculations.     

Kinetics of the hydrogen abstraction CHO + Alkane → HCHO + Alkyl reaction class: an application of the reaction class transition state theory

The kinetics of the hydrogen abstraction at alkanes by formyl radicals is investigated using the reaction class transition state theory (RC-TST) approach combined with the linear energy relationship (LER) or the barrier height grouping (BHG). The rate constants of a reaction in this class can be estimated through those of the reference reaction, CHO + C₂H₆, which are obtained from rate constants of the reaction that involves the smallest species, namely CHO + CH₄, using the explicit RC-TST scaling. The thermal rate constants of this smallest reaction are evaluated at the canonical variational transition state theory (CVT) with the corrections from the small-curvature tunneling (SCT) and hindered rotation (HR) treatments. Our analyses indicate that less than 40% systematic errors, on the average, exist in the predicted rate constants using both the LER approach, where only reaction energy is needed, and the BHG approach, where no additional information is needed; while comparing to explicit rate calculations the differences are less than 60%. Contribution to Mark S. Gordon 65th Birthday Festschrift Issue.     

Kinetics of the hydrogen abstraction R?OH + H → R?O• + H2 reaction class

This paper presents an application of the reaction class transition state theory (RC‐TST) to predict thermal rate constants for the hydrogen abstraction R? OH + H → R? O^? + H₂ reaction class, where R is an alkyl group. We have derived all parameters for the RC‐TST method for this reaction class from rate constants of 19 representative reactions, coupling with linear energy relationships (LERs) and the barrier height grouping (BHG) approach. Error analyses indicate that the RC‐TST/LER, where only reaction energy is needed, and RC‐TST/BHG, where no other information is needed, can predict rate constants for any reaction in this reaction class with satisfactory accuracy for combustion modeling. Specifically for this reaction class, the RC‐TST/LER method has less than 25% systematic errors in the predicted rate constants, whereas the RC‐TST/BHG method has less than 35% error when compared to explicit rate calculations. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 42: 414–429, 2010     

Kinetics of 1,5-hydrogen migration in alkyl radical reaction class

Kinetics of the 1,5-intramolecular hydrogen migration in the alkyl radicals reaction class has been studied using the reaction class transition state theory combined with the linear energy relationship (LER) and the barrier height grouping (BHG) approach. The high pressure limits of the rate constants for the reference reaction of 1-pentyl → 1-pentyl, calculated by the Canonical Variational Transition State Theory (CVT) with the Small Curvature Tunneling (SCT), are taken from the literature. Direct comparison with available experimental data indicates that the RC-TST/LER, where only reaction energy is needed, can predict rate constants for any reaction in this reaction class with excellent accuracy. Specifically for this reaction class, the RC-TST/LER method has less than 65% systematic errors in the predicted rate constants when compared to explicit rate calculations.     

Kinetics of the hydrogen abstraction OH + alkane --> H2O + alkyl reaction class: an application of the reaction class transition state theory

This paper presents an application of the reaction class transition state theory (RC-TST) to predict thermal rate constants for hydrogen abstraction reactions of the type OH + alkane --> HOH + alkyl. We have derived all parameters for the RC-TST method for this reaction class from rate constants of 19 representative reactions, coupling with linear energy relationships (LERs), so that rate constants for any reaction in this class can be predicted from its reaction energy calculated at either the AM1 semiempirical or BH&HLYP/cc-pVDZ level of theory. The RC-TST/LER thermal rate constants for selected reactions are in good agreement with those available in the literature. Detailed analyses of the results show that the RC-TST/LER method is an efficient method for accurately estimating rate constants for a large number of reactions in this class. Analysis of the LERs leads to the discovery of the beta-carbon radical stabilization effect that stabilizes the transition state of any reaction in this class that yields products having one or more beta-carbons, and thus leads to the lower barrier for such a reaction.     

Kinetics of the hydrogen abstraction R−OH + H → R•−OH + H2 reaction class: An application of the reaction class transition state theory

This paper presents an application of the reaction class transition state theory (RC‐TST) to predict thermal rate constants for hydrogen abstraction reactions of the type R‐OH + H → R^?‐OH + H₂. We have derived all parameters for the RC‐TST method with linear energy relationships (LERs) and the barrier height grouping (BHG) approach for this reaction class from rate constants of 37 representative reactions divided in two types of hydrogen abstraction, namely from α carbon sites and non‐α carbon sites two training sets. Error analyses indicate that the RC‐TST/LER, where only reaction energy is needed, and RC‐TST/ BHG, where no other information is needed, can predict rate constants for any reaction in this reaction class with satisfactory accuracy for combustion modeling. Specifically for this reaction class, the RC‐TST/LER and RC‐TST/BHG methods have, respectively, less than 40% and 90% systematic errors in the predicted rate constants, when compared to the explicit full TST/Eckart method. The branching ratio analysis shows that in the low‐temperature regime α abstractions are dominant, whereas, for T > 1500 K, abstractions at other sites become more important. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 43: 78–98, 2011     

Kinetics of the hydrogen abstraction *CH3 + alkane --> CH4 + alkyl reaction class: an application of the reaction class transition state theory

Kinetics of the hydrogen abstraction reaction (*)CH(3) + CH(4) --> CH(4) + (*)CH(3) is studied by a direct dynamics method. Thermal rate constants in the temperature range of 300-2500 K are evaluated by the canonical variational transition state theory (CVT) incorporating corrections from tunneling using the multidimensional semiclassical small-curvature tunneling (SCT) method and from the hindered rotations. These results are used in conjunction with the Reaction Class Transition State Theory/Linear Energy Relationship (RC-TST/LER) to predict thermal rate constants of any reaction in the hydrogen abstraction class of (*)CH(3) + alkanes. Our analyses indicate that less than 40% systematic errors on the average exist in the predicted rate constants using the RC-TST/LER method while comparing to explicit rate calculations the differences are less than 100% or a factor of 2 on the average.     

烷基环己烷夺氢反应类的动力学研究

应用反应类过渡态理论(RC-TST)研究了氢自由基抽取烷基环己烷氢的反应类,给出了该类反应的动力学参数.通过研究该反应类中的14个代表反应,分析得到了RC-TST方法所需要的各种参数,建立了该类反应的线性能量关系式(LER).应用RC-TST/LER方法计算得到速率常数与TST/Eckart方法计算结果吻合很好,分析表明RC-TST/LER方法对于准确预测夺氢反应类的速率常数十分有效,可以节约大量计算成本.     

Seven dimensional quantum dynamics study of the H2+NH2-->H+NH3 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H2+NH2-->H+NH3 reaction using a seven dimensional model on an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values and nonreactive NH2 group keeps C2v symmetry and the rotation-vibration coupling in NH2 is neglected. The total reaction probabilities are calculated when the two reactants are initially at their ground states, when the NH2 bending mode is excited, and when H2 is on its first vibrational excited state, with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants and equilibrium constants are calculated for the temperature range of 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the reaction is dominated by ground-state reactivity and the main contribution to the thermal rate constants is thought to come from this state, (b) the excitation energy of H2 was used to enhance reactivity while the excitation of the NH2 bending mode hampers the reaction, (c) the calculated thermal rate constants are very close to the experimental data and transition state theory results at high and middle temperature, while they are ten times higher than that of transition state theory at low temperature (T=200 K), and (d) the equilibrium constants results indicate that the approximations applied may have different roles in the forward and reverse reactions.     

Prediction of product branching ratios in the C(3P)+C2H2-->l-C3H+H/c-C3H+H/C3+H2 reaction using ab initio coupled clusters calculations extrapolated to the complete basis set combined with Rice-Ramsperger-Kassel-Marcus and radiationless transition theories

Ab initio CCSD(T) calculations of intermediates and transition states on the singlet and triplet C3H2 potential energy surfaces extrapolated to the complete basis set limit are combined with statistical computations of energy-dependent rate constants of the C(3P)+C2H2 reaction under crossed molecular beam conditions. Rice-Ramsperger-Kassel-Marcus theory is applied for isomerization and dissociation steps within the same multiplicity and radiationless transition and nonadiabatic transition state theories are used for singlet-triplet intersystem crossing rates. The calculated rate constants are utilized to predict product branching ratios. The results demonstrate that, in qualitative agreement with available experimental data, c-C3H+H and C3+H2 are the most probable products at low collision energies, whereas l-C3H+H becomes dominant at higher Ec above approximately 25 kJ/mol.     

Seven-dimensional quantum dynamics study of the H+NH3-->H2+NH2 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H+NH3-->H2+NH2 reaction using a seven-dimensional model and an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values. The total reaction probabilities are calculated for the initial ground and seven excited states of NH3 with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants are calculated for the temperature range 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the total reaction probabilities are overall very small, (b) the symmetric and asymmetric NH stretch excitations enhance the reaction significantly and almost all of the excited energy deposited was used to reduce the reaction threshold, (c) the excitation of the umbrella and bending motion have a smaller contribution to the enhancement of reactivity, (d) the main contribution to the thermal rate constants is thought to come from the ground state at low temperatures and from the stretch excited states at high temperatures, and (e) the calculated thermal rate constants are three to ten times smaller than the experimental data and transition state theory results.     

Ab initio rate constants from hyperspherical quantum scattering: application to H+C2H6 and H+CH3OH

The dynamics and kinetics of the abstraction reactions of H atoms with ethane and methanol have been studied using a quantum mechanical procedure. Bonds being broken and formed are treated with explicit hyperspherical quantum dynamics. The ab initio potential energy surfaces for these reactions have been developed from a minimal number of grid points (average of 48 points) and are given by analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimized using the second order perturbation theory method with a correlation consistent polarized valence triple zeta basis set. Single point energies are calculated on the optimized geometries with the coupled cluster theory and the same basis set. The reaction of H with C2H6 is endothermic by 1.5 kcal/mol and has a vibrationally adiabatic barrier of 12 kcal/mol. The reaction of H with CH3OH presents two reactive channels: the methoxy and the hydroxymethyl channels. The former is endothermic by 0.24 kcal/mol and has a vibrationally adiabatic barrier of 13.29 kcal/mol, the latter reaction is exothermic by 7.87 kcal/mol and has a vibrationally adiabatic barrier of 8.56 kcal/mol. We report state-to-state and state-selected cross sections together with state-to-state rate constants for the title reactions. Thermal rate constants for these reactions exhibit large quantum tunneling effects when compared to conventional transition state theory results. For H+CH3OH, it is found that the CH2OH product is the dominant channel, and that the CH3O channel contributes just 2% at 500 K. For both reactions, rate constants are in good agreement with some measurements.     

变分过渡态理论对H＋H2及其同位素选态反应的研究

将选态速度常数的计算推广到任意指定反应物、过渡态的振动激发态.用此法计算了H+H_2(v)及其同位素经不同振动激发过渡态时的速度常数,发现弯曲振动模激发所得结果与实验值更符合,并且在给定能量下,过渡态的弯曲振动模激发比其对称伸缩模激发更有利于反应进行.     

Variational analysis of the phenyl + O2 and phenoxy + O reactions

Variational transition state analysis was performed on the barrierless phenyl + O2 and phenoxy + O association reactions. In addition, we also calculated rate constants for the related vinyl radical (C2H3) + O2 and vinoxy radical (C2H3O) + O reactions and provided rate constant estimates for analogous reactions in substituted aromatic systems. Potential energy scans along the dissociating C-OO and CO-O bonds (with consideration of C-OO internal rotation) were obtained at the O3LYP/6-31G(d) density functional theory level. The CO-O and C-OO bond scission reactions were observed to be barrierless, in both phenyl and vinyl systems. Potential energy wells were scaled by G3B3 reaction enthalpies to obtain accurate activation enthalpies. Frequency calculations were performed for all reactants and products and at points along the potential energy surfaces, allowing us to evaluate thermochemical properties as a function of temperature according to the principles of statistical mechanics and the rigid rotor harmonic oscillator (RRHO) approximation. The low-frequency vibrational modes corresponding to R-OO internal rotation were omitted from the RRHO analysis and replaced with a hindered internal rotor analysis using O3LYP/6-31G(d) rotor potentials. Rate constants were calculated as a function of temperature (300-2000 K) and position from activation entropies and enthalpies, according to canonical transition state theory; these rate constants were minimized with respect to position to obtain variational rate constants as a function of temperature. For the phenyl + O2 reaction, we identified the transition state to be located at a C-OO bond length of between 2.56 and 2.16 A (300-2000 K), while for the phenoxy + O reaction, the transition state was located at a CO-O bond length of 2.00-1.90 A. Variational rate constants were fit to a three-parameter form of the Arrhenius equation, and for the phenyl + O2 association reaction, we found k(T) = 1.860 x 1013T-0.217 exp(0.358/T) (with k in cm3 mol-1 s-1 and T in K); this rate equation provides good agreement with low-temperature experimental measurements of the phenyl + O2 rate constant. Preliminary results were presented for a correlation between activation energy (or reaction enthalpy) and pre-exponential factor for heterolytic O-O bond scission reactions.     

Quantum dynamics study of H+NH3-->H2+NH2 reaction

We report in this paper a quantum dynamics study for the reaction H+NH3-->NH2+H2 on the potential energy surface of Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The quantum dynamics calculation employs the semirigid vibrating rotor target model [J. Z. H. Zhang, J. Chem. Phys. 111, 3929 (1999)] and time-dependent wave packet method to propagate the wave function. Initial state-specific reaction probabilities are obtained, and an energy correction scheme is employed to account for zero point energy changes for the neglected degrees of freedom in the dynamics treatment. Tunneling effect is observed in the energy dependency of reaction probability, similar to those found in H+CH4 reaction. The influence of rovibrational excitation on reaction probability and stereodynamical effect are investigated. Reaction rate constants from the initial ground state are calculated and are compared to those from the transition state theory and experimental measurement.     

Association rate constants for reactions between resonance-stabilized radicals: C3H3 + C3H3, C3H3 + C3H5, and C3H5 + C3H5

Reactions between resonance-stabilized radicals play an important role in combustion chemistry. The theoretical prediction of rate coefficients and product distributions for such reactions is complicated by the fact that the initial complex-formation steps and some dissociation steps are barrierless. In this paper direct variable reaction coordinate transition state theory (VRC-TST) is used to predict accurately the association rate constants for the self and cross reactions of propargyl and allyl radicals. For each reaction, a set of multifaceted dividing surfaces is used to account for the multiple possible addition channels. Because of their resonant nature the geometric relaxation of the radicals is important. Here, the effect of this relaxation is explicitly calculated with the UB3LYP/cc-pvdz method for each mutual orientation encountered in the configurational integrals over the transition state dividing surfaces. The final energies are obtained from CASPT2/cc-pvdz calculations with all pi-orbitals in the active space. Evaluations along the minimum energy path suggest that basis set corrections are negligible. The VRC-TST approach was also used to calculate the association rate constant and the corresponding number of states for the C(6)H(5) + H --> C(6)H(6) exit channel of the C(3)H(3) + C(3)H(3) reaction, which is also barrierless. For this reaction, the interaction energies were evaluated with the CASPT2(2e,2o)/cc-pvdz method and a 1-D correction is included on the basis of CAS+1+2+QC/aug-cc-pvtz calculations for the CH(3) + H reference system. For the C(3)H(3) + C(3)H(3) reaction, the VRC-TST results for the energy and angular momentum resolved numbers of states in the entrance channels and in the C(6)H(5) + H exit channel are incorporated in a master equation simulation to determine the temperature and pressure dependence of the phenomenological rate coefficients. The rate constants for the C(3)H(3) + C(3)H(3) and C(3)H(5) + C(3)H(5) self-reactions compare favorably with the available experimental data. To our knowledge there are no experimental rate data for the C(3)H(3) + C(3)H(5) reaction.     

HNCO+OH->H2O+NCO的反应机理

采用从头算分子轨道法 (UHF/6 31G 水平 ,并用MP4加以相关能校正 )研究了HNCO OHH2 O NCO反应机理 .同时用Morokuma数值法获得了反应途径即内禀反应坐标 (IRC) .沿着IRC ,运用反应途径哈密顿理论 ,获得反应途径动态学信息 .在此基础上 ,根据过渡态理论和相应隧道效应校正 ,计算了在不同温度下的反应速率常数 ,得到了和实验相一致的结果 .计算结果表明 ,此反应是一步直接型的抽提H反应 .     

A coupled states calculation of accurate quantum rate constants for H + H2

In this article we present the results of converged quantum reactivescattering calculations of thermal rate constants for H + H₂ using the Liu-Siegbahn-Truhlar-Horowitz (LSTH) potential energy surface. These calculations are based on the coupled states (CS) approximation wherein rotational states having different body fixed angular momentum projection quantum numbers are decoupled. By comparision with accurate coupled channel results on the Porter-Karplus No. 2 (PK2) potential surface, we estimate that the maximumerror in thermal rate constants arising from both this approximation and from other numerical approximations in the calculation is less than 25%. We also show that the sum over projection quantum numbers Ω associated with the CS calculation may be approximated quite accurately in terms of the Ω = 0 rate constants by assuming that the |Ω| > 0 rate constants differ from Ω = 0 by a shift in activation energy, which reflects the vibrationally adiabatic bending energy associated with each Ω. Comparison of the LSTH rate constants with experiment indicates average errors of 16% and 24% relative to the two modern measurementsof the rate constants for H + H₂. These errors are reduced to 18% and 9% if the CS rate constants are multiplied by exp(0.0065 eV/kT). The expected error (based on recent quantum structure calculations) associated with the 0.425 eV barrier of theLSTH potential surface is 0.0065 eV. Overall, the agreement of either the LSTH or modified LSTH rate constants with experiment is within the 32% maximum disagreement between the two experimental measurements at all butthe lowest temperature that has been studied. Comparison of our CS rate constants with the results from simpler theories is considered using both the LSTH and PK2 potential surfaces. The best overall agreement is with transition state theories that use accurate dynamical methods to calculate tunnelling factors. These include reduced dimensionality quantum dynamics methods and variational transition state theory using either the Marcus-Coltrin or least action ground-state tunnelling paths. Comparison with the results of quasiclassical trajectory calculations indicates substantial errors at low temperatures.     

H2 + CN(n=0,1)→H + HCN振动选态反应

用从头算方法获得了Ｈ２＋ＣＮ反应的内禀反应坐标（ＩＲＣ），沿着ＩＲＣ，计算了各垂直于ＩＲＣ的简正模所对应的频率（Ｗ）以及沿ＩＲＣ运动与垂直ＩＲＣ运动的简正模之间的耦合常数（ＢＫＦ），根据传统过渡态，变分过渡态理论和选态公式，计算了ｎＣＮ＝０及ｎＣＮ＝１时反应的速率常数，并得到了实验相一致的结果，还计算了ｎＣＨ＝１及ｎＣＮ＝１的Ｈ＋ＨＣＮ→Ｈ２＋ＣＮ反应速率常数，可供实验工作者参考。