Direct dynamics studies on the hydrogen abstraction reactions CF3O+CH4(CD4)→CF3OH(CF3OD)+CH3(CD3)

The dynamics properties of the hydrogen abstraction reaction CF₃O+CH₄→CF₃OH+CH₃ are studied by dual-level direct dynamics method. Optimization calculations are preformed by B3LYP and MP2 with the 6-311G(d,p) basis set, and the single-point calculations are done at the multi-coefficient correction method based on quadratic configuration interaction with single and double excitations (MC-QCISD) method. The rate constants are evaluated by canonical variational transition-state theory with a small-curvature tunneling correction over a wide range of temperature 200–2000 K. The agreement between theoretical and experimental rate constants is good in the measured temperature range. The calculated results show that the variational effect is small and almost neglected over the whole temperature range, whereas, the tunneling correction plays a role in the lower temperature range. The kinetic isotope effect for the reaction is ‘normal’. The value of k_H/k_D is 2.38 at room temperature and it decreases with the temperature increasing.     

Direct dynamics study on hydrogen abstraction reaction of CF3CHOHCF3 with OH radical

Dual-level direct dynamics method is employed to investigate the H-abstraction reaction CF₃CHOHCF₃ with OH radical. Two hydrogen-abstraction reaction channels are possible: one from the methylene (–CH–) position and the other from the hydroxyl (–OH) position. The minimum energy path is calculated at the B3LYP/6-311G(d,p) level, and the energetic information is further refined by a new powerful and inexpensive BMC-CCSD method. To testify the accuracy of the structures and the energies, the recently developed hybrid density functional theory BB1K and higher level MC-QCISD are applied to this system. Hydrogen-bonded complexes are presented at both reactants and products sides of these two channels, which indicating that the reaction may proceed via an indirect mechanism. The rate constants for each reaction channel are evaluated by canonical variational transition state theory (CVT) with a small-curvature tunneling correction (SCT) over a wide range of temperatures from 200 to 2000 K. The calculated CVT/SCT rate constants are in good agreement with the available experimental values in the temperature region 250–430 K. The present results indicate that the two channels are competitive. At lower temperature, the reaction occurs mainly via the hydroxyl-H-abstraction channel, while the methylene-H-abstraction channel is preferred when the temperature is higher than 273 K.     

Kinetic study of the reaction H2O2+H→H2O+OH by ab initio and density functional theory calculations

Theoretical investigations on the kinetics of the elementary reaction H₂O₂+H→H₂O+OH were performed using the transition state theory (TST). Ab initio (MP2//CASSCF) and density functional theory (B3LYP) methods were used with large basis set to predict the kinetic parameters; the classical barrier height and the pre-exponential factor. The ZPE and BSSE corrected value of the classical barrier height was predicted to be 4.1 kcal mol⁻¹ for MP2//CASSCF and 4.3 kcal mol⁻¹ for B3LYP calculations. The experimental value fitted from Arrhenius expressions ranges from 3.6 to 3.9 kcal mol⁻¹. Thermal rate constants of the title reaction, based on the ab initio and DFT calculations, was evaluated for temperature ranging from 200 to 2500 K assuming a direct reaction mechanism. The modeled ab initio-TST and DFT–TST rate constants calculated without tunneling were found to be in reasonable agreement with the observed ones indicating that the contribution of the tunneling effect to the reaction was predicted to be unimportant at ambient temperature.     

Direct ab initio dynamics calculations on the rate constants for the hydrogen-abstraction reaction of C2H5F with O (3P)

The hydrogen-abstraction reaction C₂H₅F+O → C₂H₄F+OH has been studied by a dual-level direct dynamics method. For the reaction, three reaction channels, one for α-abstraction and two for β-abstraction, have been identified. The potential-energy surface information is obtained at the MP2(full)/6-311G(d,p) and PMP2(full)/6-311G(3df,3pd) (single-point) levels. By canonical variational transition-state theory, rate constants for each reaction channel are calculated with a small-curvature tunneling correction. The total rate constant is calculated from the sum of the individual rate constants and the temperature dependence of the branching ratios is obtained over a wide range of temperatures from 300 to 5,000 K. The agreement of the rate constants with experiment is good in the experimental temperature range from 1,000 to 1,250 K. The calculated results indicate that at low temperatures α-abstraction is most likely to be the major reaction channel, while β-abstraction channels will significantly contribute to the whole reaction rate as the temperature increases. Received: 23 January 2002 / Accepted: 23 June 2002 / Published online: 20 September 2002     

Ab initio molecular orbital study on the gas phase SN2 reaction F + CH3Cl → CH3F + Cl

Ab-initio molecular orbital (MO) and direct ab initio dynamics calculations have been applied to the gas phase S_N2 reaction F⁻ + CH₃Cl → CH₃F + Cl⁻. Several basis sets were examined in order to select the most convenient and best fitted basis set to that of high-quality calculations. The Hartree–Fock (HF) 3−21+G(d) calculation reasonably represents a potential energy surface calculated at the MP2/6−311++G(2df,2pd) level. A direct ab initio dynamics calculation at the HF/3−21+G(d) level was carried out for the S_N2 reaction. A full dimensional ab initio potential energy surface including all degrees of freedom was used in the dynamics calculation. Total energies and gradients were calculated at each time step. Two initial configurations at time zero were examined in the direct dynamics calculations: one is a near collinear collision, and the other is a side-attack collision. It was found that in the near collinear collision almost all total available energy is partitioned into two modes: the relative translational mode between the products (40%) and the C − F stretching mode (60%). The other internal modes of CH₃F were still in the ground state. The lifetimes of the early- and late-complexes F⁻ … CH₃Cl and FCH₃ … Cl⁻ are significantly short enough to dissociate directly to the products. On the other hand, in the side-attack collision, the relative translation energy was about 20% of total available energy.     

An ab initio study of the potential surface for the reaction H2COH → HCO + H2

The potential surface for the reaction H₂CO⁺H → HCO⁺ + H₂ has been studied by ab initio SCF calculations, using gaussian-type basis functions. A saddle point on the surface has been found, and a reaction path is proposed to explain the observed release of kinetic energy. The energy of activation and ΔE for the reaction have been estimated.     

Ab initio and density function theory computational studies of the CH4 + H → CH3 + H2 reaction

The activation barrier for the CH₄ + H → CH₃ + H₂ reaction was evaluated with traditional ab initio and Density Functional Theory (DFT) methods. None of the applied ab initio and DFT methods was able to reproduce the experimental activation barrier of 11.0-12.0 kcal/mol. All ab initio methods (HF, MP2, MP3, MP4, QCISD, QCISD(T), G1, G2, and G2MP2) overestimated the activation energy. The best results were obtained with the G2 and G2MP2 ab initio computational approaches. The zero-point corrected energy was 14.4 kcal mol⁻¹. Some of the exchange DFT methods (HFB) computed energies which were similar to the highly accurate ab initio methods, while the B3LYP hybrid DFT methods underestimated the activation barrier by 3 kcal mol⁻¹. Gradient-corrected DFT methods underestimated the barrier even more. The gradient-corrected DFT method that incorporated the PW91 correlational functional even generated a negative reaction barrier. The suitability of some computational methods for accurately predicting the potential energy surface for this hydrogen radical abstraction reaction was discussed.     

HCN + OH→CN + H2O反应理论研究

The reaction path of the reaction HCN + OH→ CN + H2O was traced with Fukui's theory of intrinsic reaction coordinate by using ab initio MO method (at UMP4/6-31G** level) with gradient technique. On this basis, the dynamics properties along the reaction path was investigated by reaction path Hamiltonian theory. The rate constants of this reaction at different temperatures were calculated by conventional and variational transition state theory with tunneling correction. The theoretically calculated rate constants are in good agreement with experimental results, this shows that the title reaction is an one step, direct reaction.     

CCl2自由基与H2O分子反应动力学研究

用213 nm激光光解CCl4产生CCl2自由基，用LP LIF技术测定了室温下基态CCl2自由基与H2O分子的反应速率常数为(5.45±0.95)×10－14 cm3•molecule－1•s－1.在G2MP2理论水平上计算了CCl2＋H2O反应的最低单重态势能面，揭示了插入与加成 消除两种反应机理，得到了三个可能的产物通道:HCl＋HClCO、HCl＋trans ClCOH以及HCl＋cis ClCOH.并用RRKM TST和传统过渡态理论计算了这三个通道的分支比及其温度效应.结果说明在低温下（273 K），插入机理的产物通道的分支比远大于加成 消除机理的产物通道， HCl＋HClCO是主要产物，分支比为77.4％，其次是HCl＋cis ClCOH，分支比为22.6％.而在高温下（3000 K），加成 消除机理的反应通道大于插入机理， HCl＋trans ClCOH分支比为82.3％.     

Theoretical studies and rate constant calculations of the reactions C2F5CHO with OH radicals and Cl atoms

The hydrogen abstraction reactions of C₂F₅CHO with OH radicals and Cl atoms have been investigated theoretically by a dual-level direct dynamics method. In this study, the optimized geometries and frequencies of the stationary points are calculated at the MP2/cc-pVDZ level of theory. The energies of the stationery points and the selected points along the minimum energy paths are further refined at the MC-QCISD level using the MP2 geometries. Complexes with energies less than those of the reactants or products are located at the entrance or the exit channels of the two reactions. This result indicates that both of reactions proceed via indirect reaction mechanisms. The enthalpies of formation for the reactant C₂F₅CHO and the product radical C₂F₅CO are estimated by isodesmic reactions at the MC-QCISD//MP2/cc-pVDZ level. At the same level, the rate constants are calculated by canonical variational transition state theory (CVT) incorporating with the small-curvature tunneling correction (SCT) in the temperature range 200–1000 K. Good agreement between the calculated and experimental rate constants is obtained at the room temperature. Due to the lack of the kinetic data of these reactions, the fitted three-parameter expressions based on the CVT/SCT rate constants within 200–1000 K are k₁ = 1.64 × 10⁻²⁴ T^4.33 exp (−566.1/T) and k₂ = 6.33 × 10⁻¹⁵ T^1.35 exp (550.3/T) cm³ molecule⁻¹ s⁻¹, respectively.     

Reaction‐path dynamics and theoretical rate constants for the reaction CH4 + O3 → HOOO + CH3

We present a direct ab initio dynamics study of thermal rate constants of the hydrogen abstraction reaction of CH₄ + O₃ → HOOO +CH₃. The geometries of all the stationary points are optimized at MPW1K/6‐31+G(d,p), MPWB1K/6‐31+G(d,p), and BHandHLYP/6‐31+G(d,p) levels of theory. The energies are refined at a multi‐high‐level method. The extended Arrhenius expression fitted from the CVT/SCT and μVT/Eckart rate constants of ozonolysis of methane in the temperature range 200–2500 K are k^CVT/SCT(T) = 5.96 × 10^?29T^4.49e^(?17321.3/T) and k^μVT/Eckart(T) = 7.92 × 10^?29T^4.46e^(?17301.7/T), respectively. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

CF_2ClC(O)OCH_2CH_3+OH的反应机理及动力学性质的理论研究

利用双水平直接动力学方法,在MCG3-MPWB//M06-2X/aug-cc-pVDZ水平上研究了CF_2ClC(0)OCH_2CH_3+OH的微观反应机理.得到了反应物CF_2ClC(O)OCH_2CH_3的5种稳定构象(RCl~RC5),并对每一构象考察了发生在-CH_3-和-CH_2-基团上的所有可能氢提取反应通道.利用改进的变分过渡态理论(ICVT)结合小曲率隧道效应校正(SCT)计算了各反应通道的速率常数,分析了各构象反应位点选择性.结果表明,对于构象RCl和RC2,低温时氢提取反应主要发生在-CH_2-基团上;而对于构象RC3RC4和RC5,发生在-CH_3基团上的氢提取反应通道在整个温度区间内占绝对优势.根据Boltzmann配分函数计算总包反应速率常数,在298 K温度下计算的体系总包反应速率常数与实验值相符,进而给出200~1000 K温度范围内拟合了速率常数的三参数Arrhenius表达式:k_(overall)=5.45×10~(25)T~(4.54)exp(-685/T).     

An ab initio study of the Cl(P)+C2H6→C2H5+HCl abstraction reaction

Electronic energies, geometries, and harmonic vibration frequencies for the reactants, products, and transition state for the Cl(³P)+C₂H₆→C₂H₅+HCl abstraction reaction were evaluated at the HF and MP2 levels using several correlation consistent polarized-valence basis sets. Single-point calculations at PMP2, MP4, QCISD(T), and CCSD(T) levels were also carried out. The values of the forward activation energies obtained at the MP4/cc-pVTZ, QCISD(T)/cc-pVTZ, and CCSD(T)/cc-pVTZ levels using the MP2/cc-pVTZ structures are equal to −0.1, −0.4, and −0.3 kcal/mol, respectively. The experimental value is equal to 0.3±0.2 kcal/mol. We found that the MP2/aug-cc-pVTZ adiabatic vibration energy for the reaction (−2.4 kcal/mol) agrees well with the experimental value −(2.2–2.6) kcal/mol. Rate constants calculated with the zeroth-order interpolated variational transition state (IVTST-0) method are in good agreement with experiment. In general, the theoretical rate constants differ from experiment by, at most, a factor of 2.6.     

光引发BrC2F4Br+C2F4调聚反应的光强影响

实验发现，在光引发BrC2F4Br＋C2F4调聚反应中，光强（或光功率密度）能影响产品分布。提出了反应机理：此反应由加成反应与复合反应组成，而链转移反应可忽略．由此进行了动力学计算，为与实验结果吻合，拟合得的加成反应BrC2F4＋C2F4的速率常数为（2±1）×107cm3•mol-1•s-1，Br（C2F4）n≧2＋C2F4的速率常数为（1．2±0．4）×107Cm3•mol-1•s-1     

Ab initio direct dynamics studies on the hydrogen abstraction reactions of CF2H2 and CF3H with F atom

The dynamic properties of the hydrogen abstraction reactions of CF₂H₂ and CF₃H with F atom are investigated in the temperature range of 182–2000 K. The minimum-energy path (MEP) is optimized at MP2/6-311 G(d, p) level, then the energy profiles are refined at the CCSD(T)/6-311++G(3df, 2pd) level (single-point). The theoretical rate constants, which are calculated by the variational transition state theory (VTST) including the small curvature tunneling (SCT) correction, are in good agreement with the experimental ones. It is found that the rate constant of the CF₂H₂ + F reaction are larger than that of the CF₃H + F reaction and the activation energies exhibit in the just opposite order. This phenomenon can be rationalized by the hardness η of the halomethane molecules. The comparison of the two reactions with the CFH₃ + F reaction is made. It is found that the rate constants decrease in the order of CFH₃ + F > CF₂H₂ + F > CF₃H + F. The effect of fluorine substitution leads to a dramatic increase in the activation energy and a decrease in the preexponential factor. We hope that present theoretical studies for these compounds can give further information concerning how fluorine substitution affects the rate constants of hydrogen abstraction reactions.     

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

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

CH与H2分子反应动力学及选态反应的理论研究

次甲基作为化学反应源曾引起广泛的兴趣.Schaefer 及其合作者于1977年对反应CH(~4Σ~-)+H_2→CH_2(~3B_1)+H 进行过量子化学研究,但是计算中限制了一些自由度.近年来,由于能量梯度方法的发展,反应途径哈密顿理论和变分过渡态理论的提出,有可能进一步对该反应进行分子反应动力学性质的研究.本文用从头算UHF/6-31G 方法和能量梯度方法首先优化出上述反应(原子编号为CH_a+H_bH_c→H_bCH_a+H_c)的过渡态;再用     

用直接动力学方法确定CF3CHFCF3与Cl原子反应的速率常数

采用直接动力学方法研究了CF3CHFCF3与氯原子反应的动力学特性,利用包含小曲率隧道效应的正则变分过渡态理论确定了200～2 000 K温度区间内反应的速率常数.     

An ab initio potential energy surface for the reaction N+ + H2 → NH+ + H

Preliminary results of ab initio unrestricted Hartree-Fock calculations for the potential energy surface for the reaction N⁺ + H₂ → NH⁺ + H are reported. For the collinear approach of N⁺ to H₂, the ³Σ^? surface has no activation barrier and has a shallow well (ca. 1 eV). For perpendicular approach (C_2v symmetry) the ³B₂ state is of high energy, the ³A₂ state has a shallow well but as the bond angle increases the ³B₁ state decreases in energy to become the state of lowest energy. Neither the collinear nor the perpendicular approaches give adiabatic pathways to the deep potential well of ³B₁ (HNH)⁺.     

Theoretical investigation on the reaction of N2O and CO catalyzed by Fe(C6H6)

The reaction of N₂O with CO, catalyzed by Fe⁺(C₆H₆) and producing N₂ and CO₂, has been investigated at the UB3LYP/6-311+G(d) level. The computation results revealed that the reaction of Fe⁺(C₆H₆), N₂O and CO, is an O-atom abstraction mechanism. For the reaction channels, the geometries and the vibrational frequencies of all species have been calculated and the frequency modes analysis also have been given to elucidate the reaction mechanism. On the basis for geometry optimizations, the thermodynamic data of these reactions channels have been calculated using the statistical theory at 295.15 K and pressure of 0.35 Torr. Using Eyring transition state theory with Wigner correction, the activation thermodynamic data, rate constant and frequency factors for the these reaction channels also have been given. The results showed that CO and N₂O do not react without catalyst and Fe⁺(C₆H₆) can excellently mediate the reaction of N₂O and CO.