Rate coefficients of the CF3CHFCF3 + H → CF3CFCF3 + H2 reaction at different temperatures calculated by transition state theory with ab initio and DFT reaction paths

The minimum energy path (MEP) of the reaction, CF₃CHFCF₃ + H → transition state (TS) → CF₃CFCF₃ + H₂, has been computed at different ab initio levels and with density functional theory (DFT) using different functionals. The computed B3LYP/6‐31++G**, BH&HLYP/cc‐pVDZ, BMK/6‐31++G**, M05/6‐31+G**, M05‐2X/6‐31+G**, UMP2/6‐31++G**, PUMP2/6‐31++G**//UMP2/6‐31++G**, RCCSD(T)/aug‐cc‐pVDZ//UMP2/6‐31++G**, RCCSD(T)/aug‐cc‐pVTZ(spd,sp)//UMP2//6‐31++G**, RCCSD(T)/CBS//M05/6‐31+G**, and RCCSD(T)/CBS//UMP2/6‐31++G** MEPs, and associated gradients and Hessians, were used in reaction rate coefficient calculations based on the transition state theory (TST). Reaction rate coefficients were computed between 300 and 1500 K at various levels of TST, which include conventional TST, canonical variational TST (CVT) and improved CVT (ICVT), and with different tunneling corrections, namely, Wigner, zero‐curvature, and small‐curvature (SCT). The computed rate coefficients obtained at different ab initio, DFT and TST levels are compared with experimental values available in the 1000–1200 K temperature range. Based on the rate coefficients computed at the ICVT/SCT level, the highest TST level used in this study, the BH&HLYP functional performs best among all the functionals used, while the RCCSD(T)/CBS//MP2/6‐31++G** level is the best among all the ab initio levels used. Comparing computed reaction rate coefficients obtained at different levels of theory shows that, the computed barrier height has the strongest effect on the computed reaction rate coefficients as expected. Variational effects on the computed rate coefficients are found to be negligibly small. Although tunneling effects are relatively small at high temperatures (～1500 K), SCT corrections are significant at low temperatures (～300 K), and both barrier heights and the magnitudes of the imaginary frequencies affect SCT corrections. © 2012 Wiley Periodicals, Inc.     

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.     

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     

Rate constants from the reaction path Hamiltonian. II. Nonseparable semiclassical transition state theory

For proton transfer reactions, the tunneling contributions to the rates are often much larger than thermally activated rates at temperatures of interest. A number of separable tunneling corrections have been proposed that capture the dependence of tunneling rates on barrier height and imaginary frequency size. However, the effects of reaction pathway curvature and barrier anharmonicity are more difficult to quantify. The nonseparable semiclassical transition state theory (TST) of Hernandez and Miller [Chem. Phys. Lett. 214, 129 (1993)] accounts for curvature and barrier anharmonicity, but it requires prohibitively expensive cubic and quartic derivatives of the potential energy surface at the transition state. This paper shows how the reaction path Hamiltonian can be used to approximate the cubic and quartic derivatives used in nonseparable semiclassical transition state theory. This enables tunneling corrections that include curvature and barrier anharmonicity effects with just three frequency calculations as required by a conventional harmonic transition state theory calculation. The tunneling correction developed here is nonseparable, but can be expressed as a thermal average to enable efficient Monte Carlo calculations. For the proton exchange reaction NH2 + CH4 <==> NH3 + CH3, the nonseparable rates are very accurate at temperatures from 300 K up to about 1000 K where the TST rate itself begins to diverge from the experimental results.     

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 Calculation and Kinetic Study for the Abstraction Reaction of H with Sihcl3

ＴｒｉｃｈｌｏｒｏｓｉｌａｎｅｉｓａｎｉｍｐｏｒｔａｎｔｍａｔｅｒｉａｌｉｎｐｌａｓｍａＣｈｅｍｉｃａｌＶａｐｏｒＤｅｐｏｓｉｔｉｏｎ (ＣＶＤ)ａｎｄｉｎｓｅｍｉｃｏｎｄｕｃｔｏｒｄｅｖｉｃｅｐｒｏｃｅｓｓ .1 4 Ｔｈｅｒｅａｃｔｉｏｎｏｆｔｒｉｃｈｌｏｒｏｓｉｌａｎｅｗｉｔｈａｔｏｍｉｃｈｙｄｒｏｇｅｎ ,ｔｈｅｓｉｍｐｌｅｓｔｆｒｅｅ ｒａｄｉｃａｌｓｐｅｃｉｅｓ,ｈａｓｄｒａｗｎｃｏｎｓｉｄｅｒａｂｌｅａｔｔｅｎｔｉｏｎ :ｋｉｎｅｔｉｃｐａｒａｍｅｔｅｒｓｆｏｒＨ ａｔｏｍｒｅａｃｔｉｏｎａｒｅｄ…     

CH(^4∑^-)+H2O→CH2(^3B1)+OH的直接动态学研究

用从头算方法计算了反应CH(^4∑^-)+H2O→CH2(^3B1)+OH的反应途径。在此基础上, 计算沿反应途径的动态学性质和正则变分过渡态理论的速率常数, 并进行隧道效应校正。结果表明, 电子相关能对反应活化位垒影响较大; 反应存在返回效应和隧道效应, 用正则变分过渡态方法和小曲率近似的隧道校正是有效的。     

Theoretical studies of the reactions of Cl atoms with CF3CH2OCH n F(3−n) (n = 1, 2, 3)

This paper presents the theoretical studies of the reactions of Cl atoms with CF₃CH₂OCH₃, CF₃CH₂OCH₂F and CF₃CH₂OCHF₂ using an ab initio direct dynamics theory. The geometries and vibrational frequencies of the reactants, complexes, transition states and products are calculated at the MP2/6-31+(d,p) level. The minimum energy path is also calculated at same level. The MC-QCISD method is carried out for further refining the energetic information. The rate constants are evaluated with the canonical variational transition state theory (CVT) and CVT with small curvature tunneling contributions in the temperature range 200–1,500 K. The results are in good agreement with experimental values.     

Theoretical study of reactions between AlH(1Σ) and HF molecule

Reaction mechanisms between AlH (¹Σ) and HF molecule are theoretically investigated. Ab initio calculations demonstrate that there are two parallel reaction channels: one is an addition reaction to give H₂AlF via the three‐membered ring transition state (TS) and the other is a dehydrogenation reaction to give AlF+H₂ via the four‐membered ring TS. The addition reaction is thermodynamically favorable and the dehydrogenation reaction is kinetically favorable. Thermodynamics and Eyring transition state theory (TST) with the Wigner correction are also used to compute the thermodynamic functions, the equilibrium constants, A factors, and the rate constants of these reaction channels at 200–1000 K. From the thermodynamics and TST calculations, it is valuable to point out that consideration on the entropy and thermal enthalpy is quite important in the study of chemical reactions on the basis of ab initio method. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 417–424, 1999     

TheRate: Program for ab initio direct dynamics calculations of thermal and vibrational-state-selected rate constants

We introduce TheRate (THEoretical RATEs), a complete application program with a graphical user interface (GUI) for calculating rate constants from first principles. It is based on canonical variational transition-state theory (CVT) augmented by multidimensional semiclassical zero and small curvature tunneling approximations. Conventional transition-state theory (TST) with one-dimensional Wigner or Eckart tunneling corrections is also available. Potential energy information needed for the rate calculations are obtained from ab initio molecular orbital and/or density functional electronic structure theory. Vibrational-state-selected rate constants may be calculated using a diabetic model. TheRate also introduces several technical advancements, namely the focusing technique and energy interpolation procedure. The focusing technique minimizes the number of Hessian calculations required by distributing more Hessian grid points in regions that are critical to the CVT and tunneling calculations and fewer Hessian grid points elsewhere. The energy interpolation procedure allows the use of a computationally less demanding electronic structure theory such as DFT to calculate the Hessians and geometries, while the energetics can be improved by performing a small number of single-point energy calculations along the MEP at a more accurate level of theory. The CH₄+H↔CH₃+H₂ reaction is used as a model to demonstrate usage of the program, and the convergence of the rate constants with respect to the number of electronic structure calculations. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1039–1052, 1998     

Hydrogen abstraction reactions of OH radicals with CH₃CH₂CH₂Cl and CH₃CHClCH₃: a mechanistic and kinetic study

The hydrogen abstraction reactions of OH radicals with CH₃CH₂CH₂Cl (R1) and CH₃CHClCH₃ (R2) have been investigated theoretically by a dual‐level direct dynamics method. The optimized geometries and frequencies of the stationary points are calculated at the B3LYP/6‐311G(d,p) level. To improve the reaction enthalpy and potential barrier of each reaction channel, the single point energy calculation is performed by the BMC‐CCSD method. Using canonical variational transition‐state theory (CVT) with the small‐curvature tunneling correction, the rate constants are evaluated over a wide temperature range of 200–2000 K at the BMC‐CCSD//B3LYP/6‐311G(d,p) level. For the reaction channels with the negative barrier heights, the rate constants are calculated by using the CVT. The calculated total rate constants are consistent with available experimental data. The results show that at lower temperatures, the tunneling correction has an important contribution in the calculation of rate constants for all the reaction channels with the positive barrier heights, while the variational effect is found negligible for some reaction channels. For reactions OH radicals with CH₃CH₂CH₂Cl (R1) and CH₃CHClCH₃ (R2), the channels of H‐abstraction from –CH₂– and –CHCl groups are the major reaction channels, respectively, at lower temperatures. With temperature increasing, contributions from other channels should be taken into account. Finally, the total rate constants are fitted by two models, i.e., three‐parameter and four‐parameter expressions. The enthalpies of formation of the species CH₃CHClCH₂, CH₃CHCH₂Cl, and CH₂CH₂CH₂Cl are evaluated by isodesmic reactions. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Direct dynamics study on the reaction of acetaldehyde with ozone

The hydrogen abstraction reaction of ozone with acetaldehyde has been studied theoretically over the temperature range 250-2500 K. Two different reactive sites of acetaldehyde molecule, CH(3) and CHO groups have been investigated, and results confirm that the CHO group is a highly reactive site. In this study, the geometries and harmonic vibrational frequencies of all stationary points are calculated at the MPW1K, BHandHLYP, and MPWB1K levels of theory. The minimum energy paths (MEPs) were obtained at the MPW1K/6-31+G(d,p) level of theory. To refine the energies along the MEPs of each channel, single-point energy calculations were performed by a higher-level energy calculation method (denoted as HL). The rate constants were evaluated based on the MEPs from the HL method in the temperature range 250-2500 K by using the conventional transition state theory (TST), the canonical variational transition state theory (CVT), the microcanonical variational transition state theory (muVT), the CVT coupled with small-curvature tunneling (SCT) correction (CVT/SCT), and the muVT coupled with Eckart tunneling correction (muVT/Eckart). The fitted three-parameter Arrhenius expressions of the calculated CVT/SCT and muVT/Eckart rate constants of the H abstraction from CHO group are k CVT/SCT(T) = 4.92 x 10(-27).T 3.77.e(-7867.0/T) and k muVT/Eckart(T) = 2.10 x 10(-27).T(3.90).e(-7706.2/T), respectively. The fitted three-parameter Arrhenius expressions of the calculated CVT/SCT and muVT/Eckart rate constants of the H abstraction from CH3 group are k(CVT/SCT)(T) = 1.27 x 10(-27).T(3.94).e(-14554.1/T) and k muVT/Eckart(T) = 1.62 x 10(-26).T(3.66).e(-15459.8/T), respectively.     

A dual‐level direct dynamics study on the hydrogen abstraction reaction of oxygen atom with methylhydrazine

The mechanism of the multichannel reaction CH₃NHNH₂ (SC₁ and SC₂) + O → products is investigated theoretically using ab initio and density functional theory, and dynamics properties are explored by a dual‐level direct dynamics method. The calculation of the potential energy surface is carried out at the BMC‐CCSD//MPW1K/6‐311G(d,p) level. Using canonical variational transition state theory with a small‐curvature tunneling correction, the rate constants of each channel are evaluated over a wide temperature range of 200–2000 K on the basis of obtained electronic structures and energy information. The total rate constants are calculated from the sum of the individual rate constants taking into account the Boltzmann distribution of two conformers. The reactivity of the H atom located in different groups is compared. © 2013 Wiley Periodicals, Inc.     

Direct Ab Initio Dynamics Study for the Hydrogen Abstraction Reaction: CH _{2} ({^{3}}B_{1}) + H2CO \rightarrow CH3 + CHO

We present a direct ab initio dynamics study of thermal rate constants of the hydrogen abstraction reaction of CH $_{2}(^{3}B_{1})$ + H₂CO $\rightarrow $ CH₃ + CHO. The MP2/cc-pVDZ method is employed to optimize the geometries of stationary points as well as the points on the minimum energy path. The energies of all the points were further refined at the CCSD(T)/cc-pVTZ level of theory based on the Moller– Plesset perturbation theory (MP2) optimized geometries. The rate constants were evaluated using the conventional transition state theory, the canonical variational TST, and the improved canonical variational TST, also both including small-curvature tunneling correction in the temperature range of 300–2,500 K. The calculated results show that the rate constants have positive temperature dependence in the calculated temperature range. The calculated results show that the tunneling effect is important at low temperature region.     

Criegee自由基CH2O2与H2O反应机理及动力学的理论研究

在6-311＋G(2d,2p)水平下, 采用密度泛函理论(DFT)的B3LYP方法, 研究了Criegee 自由基CH₂O₂与H₂O的反应. 结果表明反应存在三个通道: CH₂O₂＋H₂O®HOCH₂OOH (R1); CH₂O₂＋H₂O®HCO＋OH＋H₂O (R2); CH₂O₂＋H₂O®HCHO＋H₂O₂ (R3), 各通道的势垒高度分别为43.35, 85.30和125.85 kJ/mol. 298 K下主反应通道(R1)的经典过渡态理论(TST)与变分过渡态理论(CVT)的速率常数k^TST与k^CVT均为2.47×10^－17 cm³•molecule^－1•s^－1, 而经小曲率隧道效应模型(SCT)校正后的速率常数k^CVT/SCT为 5.22×10^－17 cm³•molecule^－1•s^－1. 另外, 还给出了200～2000 K 温度范围内拟合得到的速率常数随温度变化的三参数Arrhenius方程.     

The proton‐coupled proton transfer mechanism,H2O catalysis,and hydrogen tunneling effects in the reaction of HNCH2 with HCOOH in the interstellar medium

The reaction HNCH₂ + HCOOH → H₂NCH₂COOH is supposed to be an important reaction related to the possible origin of amino acids on the early Earth. We find that it has an energy barrier of 87.37 kcal mol⁻¹ obtained with MP2/6‐311+G** in the gas phase, but it is likely enhanced to occur in the interstellar medium (ISM) through a proton‐coupled proton transfer reaction, initiated by HNCH₂ coupled with H₂⁺, H₃⁺, or H₃⁺O. H₂⁺, H₃⁺, and H₃⁺O serve as a donor of energy in the coupled reactions. H⁺, which is a key species to the coupled reactions, further, plays a catalytic role in reducing a barrier up to 14.14 kcal mol⁻¹. In the coupled reaction with H₃⁺O, H₂O, which can seize, transport, and deliver a proton from HCOOH to H₂NCH₂⁺, reduces a barrier up to 14.96 kcal mol⁻¹. A significant hydrogen‐tunneling pathway is predicted by the temperature dependences of k_H^CVT/SCT, calculated using the small curvature tunneling (SCT) approximation and canonical variational transition state theory (CVT). Hydrogen tunneling is another important mechanism to make the reaction happen in the ISM. The achieved results can be applied to discuss the origin of amino acids from the materials of the Earth itself. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Reaction path dynamics and theoretical rate constants for the SiH3Cl+H→SiH2Cl+H2 reaction by ab initio direct dynamics method

Ab initio direct dynamics method has been used to study the title reaction. Electronic structure information including geometries, gradients and force constants (Hessians) are calculated at the UQCISD/6-311+G^** level. Energies along the minimum energy path are improved by a series of single-point G2//QCISD calculations. The changes of the geometries, vibratioanal frequencies, potential energies and total curvature along the reaction path are discussed. The rate constants in the temperature range 200–3000 K are calculated by canonical variational transition state theory with small-curvature tunneling correction (CVT/SCT) method. The results show that the variational effect is small and in the lower temperature range, the small curvature tunneling effect is important for the reaction.     

On the kinetic mechanism of the hydrogen abstraction reactions of the hydroxyl radical with CH3CF2Cl and CH3CFCl2: A dual level direct dynamics study

By means of the dual‐level direct dynamics method, the mechanisms of the reactions, CH₃CF₂Cl + OH → products (R1) and CH₃CFCl₂ + OH → products (R2), are studied over a wide temperature range 200–2000 K. The optimized geometries and frequencies of the stationary points are calculated at the MP2/6‐311G(d,p) level, and then the energy profiles of the reactions are refined with the interpolated single‐point energy method at the G3(MP2) level. The canonical variational transition‐state theory with the small‐curvature tunneling (SCT) correction method is used to calculate the rate constants. For the title reactions, three reaction channels are identified and the H‐abstraction channel is the major pathway. The results indicate that F substitution has a significant (reductive) effect on hydrochlorofluorocarbon reactivity. Also, for all H‐abstraction reaction channels the variational effect is small and the SCT effect is only important in the lower temperature range on the rate constants calculation. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

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),量子力学隧道效应的影响在低温段非常明显,而变分效应对反应速率常数的影响很小.     

Direct dynamics studies for the reactions of CF3CHFCF3 and CF3CF2CHF2 with H atoms

The rate constants of the hydrogen abstraction reactions of CF₃CHFCF₃ + H (R1) and CF₃CF₂CHF₂ + H (R2) have been calculated by means of the dual-level direct dynamics method. Optimized geometries and frequencies of stationary points and extra points along the minimum-energy path (MEP) are obtained at the MPW1K/6-311+G(d,p) level, and the classical energetic information is further corrected with the interpolated single-point energy (ISPE) approach by the G3(MP2) level of theory. Using the canonical variational transition state theory (CVT) with small-curvature tunneling corrections (SCT), the rate constants are evaluated over a wide temperature range of 200-2000 K. The calculated CVT/SCT rate constants are in good agreement with available experimental values. It is found that the variational effect is very small and almost negligible over the whole temperature region. However, the small-curvature tunneling correction plays an important role in the lower temperature range. Furthermore, the heats of formation of species CF₃CF₂CHF₂ (SC₁ or SC₂) and CF₃CF₂CF₂ are studied using isodesmic reactions to further elucidate the thermodynamic properties.