New analytical potential energy surface for the F(2P)+CH4 hydrogen abstraction reaction: kinetics and dynamics

A new potential energy surface for the gas-phase F(2P)+CH4 reaction and its deuterated analogues is reported, and its kinetics and dynamics are studied exhaustively. This semiempirical surface is completely symmetric with respect to the permutation of the four methane hydrogen atoms, and it is calibrated to reproduce the topology of the reaction and the experimental thermal rate constants. For the kinetics, the thermal rate constants were calculated using variational transition-state theory with semiclassical transmission coefficients over a wide temperature range, 180-500 K. The theoretical results reproduce the experimental variation with temperature. The influence of the tunneling factor is negligible, due to the flattening of the surface in the entrance valley, and we found a direct dependence on temperature, and therefore positive and small activation energies, in agreement with experiment. Two sets of kinetic isotope effects were calculated, and they show good agreement with the sparse experimental data. The coupling between the reaction coordinate and the vibrational modes shows qualitatively that the FH stretching and the CH3 umbrella bending modes in the products appear vibrationally excited. The dynamics study was performed using quasi-classical trajectory calculations, including corrections to avoid zero-point energy leakage along the trajectories. First, we found that the FH(nu',j') rovibrational distributions agree with experiment. Second, the excitation function presents an oscillatory pattern, reminiscent of a reactive resonance. Third, the state specific scattering distributions present reasonable agreement with experiment, and as the FH(nu') vibrational state increases the scattering angle becomes more forward. These kinetics and dynamics results seem to indicate that a single, adiabatic potential energy surface is adequate to describe this reaction, and the reasonable agreement with experiment (always qualitative and sometimes quantitative) lends confidence to the new surface.     

Potential energy surface, kinetics, and dynamics study of the Cl+CH4-->HCl+CH3 reaction

A modified and recalibrated potential energy surface for the gas-phase Cl+CH4-->HCl+CH3 reaction is reported and tested. It is completely symmetric with respect to the permutation of the four methane hydrogen atoms and is calibrated with respect to updated experimental and theoretical stationary point properties and experimental forward thermal rate constants. From the kinetics point of view, the forward and reverse thermal rate constants and the activation energies were calculated using the variational transition-state theory with semiclassical transmission coefficients over a wide temperature range of 150-2500 K. The theoretical results reproduce the available experimental data, with a small curvature of the Arrhenius plot which indicates the role of tunneling in this hydrogen abstraction reaction. A dynamics study was also performed on this PES using quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories. First, we found a noticeable internal energy in the coproduct methyl radical, both in the ground-state [CH4 (v=0)] and vibrationally excited [CH4 (v=1)] reactions. This CH3 internal energy was directly precluded in some experiments or oversimplified in previous theoretical studies using pseudotriatomic models. Second, our QCT calculations give HCl rotational distributions slightly hotter than those in experiment, but correctly describing the experimental trend of decreasing the HCl product rotation excitation in going from HCl (v'=0) to HCl (v'=1) for the CH4 (v=1) reaction. Third, the state specific scattering distributions present a reasonable agreement with experiment, although they tend to make the reaction more forward and backward scattered than found experimentally probably because of the hotter rotational distribution and the deficiencies of the QCT methods.     

Potential energy surface for the CCl4 + H --> CCl3 + ClH reaction: kinetics and dynamics study

An analytical potential energy surface for the gas-phase CCl4 + H --> CCl3 + ClH reaction was constructed with suitable functional forms to represent vibrational modes. This surface is completely symmetric with respect to the permutation of the four chlorine atoms and is calibrated with respect to experimental thermal rate constants available over the temperature range 297-904 K. On this surface, the thermal rate constants were calculated using variational transition-state theory with semiclassical transmission coefficients over a wider temperature range 300-2500 K, therefore obtaining kinetics information at higher temperatures than are experimentally available. This surface was also used to analyze dynamical features, such as tunneling and reaction-path curvature. In the first case, the influence of the tunneling factor is very small since a heavy chlorine atom has to pass through the barrier. In the second, it was found that vibrational excitation of the Cl-H stretching mode can be expected in the exit channel.     

Quasiclassical trajectory study of the F + CH4 reaction dynamics on a dual-level interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the F + CH4 reactive system has been constructed using the interpolation method of Collins and co-workers. The ab initio calculations have been performed using second-order M?ller-Plesset (MP2) perturbation theory to build the initial PES. Scaling all correlation (SAC) methodology has been employed to improve the ab initio calculations and to construct a dual-level PES. Using this PES, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations and internal energy distributions has been carried out for the F + CH4 and F + CD4 reactions and the theoretical results have been compared with the available experimental data.     

Direct dynamic study on the hydrogen abstraction reaction CH3CN + OH-->CH2CN + H2O

The reaction of acetonitrile with hydroxyl has been studied using the direct ab initio dynamics methods. The geometries, vibrational frequencies of the stationary points, as well as the minimum energy paths were computed at the BHandHLYP and MP2 levels of theory with the 6-311G(d, p) basis set. The energies were further refined at the PMP4/6-311+G(2df, 2pd) and QCISD(T)/6-311+G(2df, 2pd) levels of theory based on the structures optimized at BHandHLYP/6-311G(d, p) and MP2/6-311G(d, p) levels of theory. The Polyrate 8.2 program was employed to predict the thermal rate constants using the canonical variational transition state theory incorporating a small-curvature tunneling correction. The computed rate constants are in good agreement with the available experimental data.     

Nonadiabatic reactant-product decoupling calculation for the F(2P(1/2)) + H2 reaction

In this paper we present a theoretical study using time-dependent nonadiabatic reactant-product decoupling method for the state-to-state reactive scattering calculation of F((2)P(1/2))+H(2) (nu=j=0) reaction on the Alexander-Stark-Werner potential energy surface. In this nonadiabatic state-to-state calculation, the full wave function is partitioned into reactant component and a sum of all product components. The reactant and product components of the wave function are solved independently. For the excited state reaction, the state-to-state reaction probabilities for J=0.5 are calculated. Comparing the state-to-state reaction probabilities, it is found that the vibrational population of the HF product is dominated by vibrational levels nu=2 and 3. The rotation specific reaction probabilities of HF product in j=1 and 2 are larger than those in other rotational levels. As the rotation quantum number j increases, the positions of the peak in the rotational reaction probability of HF product in nu=3 shift to higher collision energy.     

An ab initio spin-orbit-corrected potential energy surface and dynamics for the F + CH4 and F + CHD3 reactions

We report an analytical ab initio three degrees of freedom (3D) spin-orbit-correction surface for the entrance channel of the F + methane reaction obtained by fitting the differences between the spin-orbit (SO) and non-relativistic electronic ground state energies computed at the MRCI+Q/aug-cc-pVTZ level of theory. The 3D model surface is given in terms of the distance, R(C-F), and relative orientation, Euler angles ? and θ, of the reactants treating CH(4) as a rigid rotor. The full-dimensional (12D) "hybrid" SO-corrected potential energy surface (PES) is obtained from the 3D SO-correction surface and a 12D non-SO PES. The SO interaction has a significant effect in the entrance-channel van der Waals region, whereas the effect on the energy at the early saddle point is only ～5% of that at the reactant asymptote; thus, the SO correction increases the barrier height by ～122 cm(-1). The 12D quasiclassical trajectory calculations for the F + CH(4) and F + CHD(3) reactions show that the SO effects decrease the cross sections by a factor of 2-4 at low collision energies and the effects are less significant as the collision energy increases. The inclusion of the SO correction in the PES does not change the product state distributions.     

Direct dynamics studies on hydrogen abstraction reactions of CH3CHFCH3 and CH3CH2CH2F with OH radicals

The hydrogen abstraction reactions of CH3CHFCH3 and CH3CH2CH2F with the OH radicals have been studied theoretically by a dual-level direct dynamics method. The geometries and frequencies of all the stationary points are optimized by means of the DFT calculation. There are complexes at the reactant side or exit route, indicating these reactions may proceed via indirect mechanisms. To improve the reaction enthalpy and potential barrier of each reaction channel, the single point energy calculation is performed by the MC-QCISD/3 method. The rate constants are evaluated by canonical variational transition state theory (CVT) with the small-curvature tunneling correction method (SCT) over a wide temperature range 200-2000 K. The canculated CVT/SCT rate constants are consistent with available experimental data. The results show that both the variation effect and the SCT contribution play an important role in the calculation of the rate constants. For reactions CH3CHFCH3 and CH3CH2CH2F with OH radicals, the channels of H-abstraction from -CHF- and -CH2- groups are the major reaction channels, respectively, at lower temperature. Furthermore, to further reveal the thermodynamics properties, the enthalpies of formation of reactants CH3CHFCH3, CH3CH2CH2F, and the product radicals CH3CFCH3, CH3CHFCH2, CH3CH2CHF, CH3CHCH2F, and CH2CH2CH2F are studied using isodesmic reactions.     

Direct dynamics study on the hydrogen abstraction reaction CH2O + HO2 --> CHO + H2O2

We present a direct ab initio dynamics study on the hydrogen abstraction reaction CH2O + HO2 --> CHO + H2O2, which is predicted to have four possible reaction channels caused by different attacking orientations of HO2 radical to CH2O. The structures and frequencies at the stationary points and the points along the minimum energy paths (MEPs) of the four reaction channels are calculated at the B3LYP/cc-pVTZ level of theory. Energetic information of stationary points and the points along the MEPs is further refined by means of some single-point multilevel energy calculations (HL). The rate constants of these channels are calculated using the improved canonical variational transition-state theory with the small-curvature tunneling correction (ICVT/SCT) method. The calculated results show that, in the whole temperature range, the more favorable reaction channels are Channels 1 and 3. The total ICVT/SCT rate constants of the four channels at the HL//B3LYP/cc-pVTZ level of theory are in good agreement with the available experiment data over the measured temperature ranges, and the corresponding three-parameter expression is k(ICVT/SCT) = 3.13 x 10(-20) T(2.70) exp(-11.52/RT) cm3 mole(-1) s(-1) in the temperature range of 250-3000 K. Additionally, the flexibility of the dihedral angle of H2O2 is also discussed to explain the different experimental values.     

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     

Accurate potential energy surface and quantum reaction rate calculations for the H+CH4-->H2+CH3 reaction

Calculations for the cumulative reaction probability N(E) (for J=0) and the thermal rate constant k(T) of the H+CH(4)-->H(2)+CH(3) reaction are presented. Accurate electronic structure calculations and a converged Shepard-interpolation approach are used to construct a potential energy surface which is specifically designed to allow the precise calculation of k(T) and N(E). Accurate quantum dynamics calculations employing flux correlation functions and multiconfigurational time-dependent Hartree wave packet propagation compute N(E) and k(T) based on this potential energy surface. The present work describes in detail the various convergence test performed to investigate the accuracy of the calculations at each step. These tests demonstrate the predictive power of the present calculations. In addition, approximate approaches for reaction rate calculations are discussed. A quite accurate approximation can be obtained from a potential energy surface which includes only interpolation points on the minimum energy path.     

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

The dual-level direct dynamics method has been employed to investigate the H-abstraction reaction of CF(3)CF(2)CH(2)OH with OH radical, which is predicted to have two classes of possible reaction channels caused by different positions of hydrogen atom attack. The minimum-energy path is calculated at the B3LYP/6-311G(d,p) level, and the energetic information is further refined by the MC-QCISD method. To compare the structures, the other method MPW1K/6-311G(d,p) is also applied to this system. Hydrogen-bonded complexes are presented in the reactant and product sides of the three channels, indicating that each reaction may proceed via an indirect mechanism. The rate constants for each reaction channel are evaluated by canonical variational transition-state theory (CVT) with the small-curvature tunneling correction (SCT) over a wide range of temperatures from 200 to 2000 K. The calculated CVT/SCT rate constants are found to be in good agreement with the available experimental values. The result shows that the variational effect is small, and in the lower-temperature range, the SCT effect is important for each reaction. It is shown that hydrogen abstracted from the -CH(2)- position is the major channel, while H-abstraction from the -OH position may be neglected with the temperature increasing.     

Direct dynamics study on the hydrogen abstraction reactions N2H4 + R-->N2H3 + RH (R=NH2, CH3)

We present a direct ab initio dynamics study on the hydrogen abstraction reactions N(2)H(4)+R-->N(2)H(3)+RH (R=NH(2),CH(3)), which are predicted to have six possible reaction channels for NH(2) abstraction and four for CH(3) abstraction caused by the different N(2)H(4) isomers and various attacking orientations of foreign radical to N(2)H(4). The structures and frequencies at the stationary points and the points along the minimum energy paths (MEPs) of all reaction channels are obtained at the UMP2(full)6-31+G(d,p) level of theory. Energetic information of stationary points and the points along the MEPs is further refined by means of MC-QCISD method. The rate constants of these channels are calculated using the improved canonical variational transition-state theory with the small-curvature tunneling correction (ICVT/SCT) method. The calculated results show that the favorable reaction channels are channels (n1) and (n4) as well as (c1) and (c3) (refer to Scheme 1) in the whole temperature range. The total ICVT/SCT rate constants of all channels for the two reactions at the MC-QCISDUMP2(full)6-31+G(d,p) level are both in good agreement with the available experimental data, and corresponding three-parameter expressions of k(ICVTSCT) in 220-3000 K are fitted as 6.46 x 10(-15)(T298)(3.60) exp(-386T) cm(3) mol(-1) s(-1) for NH(2) abstraction and 1.04 x 10(-14)(T298)(4.00) exp(-2037T) cm(3) mol(-1) s(-1) for CH(3) abstraction. Additionally, the long range interaction between the H atom of X-H bond in foreign radicals and the lone pair on the nonreactive N atom of the transition states is further discussed to explain the various transition-state numbers of the two similar hydrogen abstraction reactions.     

Seven-dimensional quantum dynamics study of the O(3P)+CH4 reaction

The initial state selected time-dependent wave packet calculations have been carried out to study the title reaction with seven degrees of freedom included by restricting the nonreacting CH(3) group under C(3V) symmetry and the CH bond length in the group. Total reaction probabilities as well as integral cross sections were calculated for the ground and four vibrationally excited reagent states. Our calculation shows that the reactivity is very small for the reaction for collision energy up to 1.0 eV for all the initial states. Initial vibration excitation of CH(4), in particular, the CH stretch excitation, enhances the reactivity, but only part of the excitation energy deposited can be used to reduce the reaction threshold. The rate constant for the ground initial state agrees rather well with that from a recent quasiclassical trajectory study and is larger than that from the semirigid vibrating rotor target calculations, in particular, in the low temperature region. On the other hand, the thermal rate constant calculated from the integral cross sections for these five vibrational states is about a factor of 20 smaller than that from the multiconfiguration time-dependent Hartree study.     

Imaging photon-initiated reactions: a study of the Cl(2P(3/2))+CH4-->HCl+CH3 reaction

The hydrogen or deuterium atom abstraction reactions between Cl((2)P(3/2)) and methane, or its deuterated analogues CD(4) and CH(2)D(2), have been studied at mean collision energies around 0.34 eV. The experiments were performed in a coexpansion of molecular chlorine and methane in helium, with the atomic Cl reactants generated by polarized laser photodissociation of Cl(2) at 308 nm. The Cl-atom reactants and the methyl radical products were detected using (2+1) resonantly enhanced multiphoton ionization, coupled with velocity-map ion imaging. Analysis of the ion images reveals that in single-beam experiments of this type, careful consideration must be given to the spread of reagent velocities and collision energies. Using the reactions of Cl with CH(4), CD(4), and CH(2)D(2), as examples, it is shown that the data can be fitted well if the reagent motion is correctly described, and the angular scattering distributions can be obtained with confidence. New evidence is also provided that the CD(3) radicals from the Cl+CD(4) reaction possess significant rotational alignment under the conditions of the present study. The results are compared with previous experimental and theoretical works, where these are available.     

Nonadiabatic dynamics in the CH3+HCl-->CH4+Cl(2PJ) reaction

Nonadiabatic dynamics in the title reaction have been investigated by 2+1 REMPI detection of the Cl(2P(3/2)) and Cl*(2P(1/2)) products. Reaction was initiated by photodissociation of CH(3)I at 266 nm within a single expansion of a dilute mixture of CH(3)I and HCl in argon, giving a mean collision energy of 7800 cm(-1) in the center-of-mass frame. Significant production of Cl* was observed, with careful checks made to ensure that no additional photochemical or inelastic scattering sources of Cl* perturbed the measurements. The fraction of the total yield of Cl(2P(J)) atoms formed in the J=1/2 level at this collision energy was 0.150+/-0.024, and must arise from nonadiabatic dynamics because the ground potential energy surface correlates to CH(4)+Cl(2P(3/2)) products.     

Ab initio potential energy surface of CH and reaction dynamics of H + CH+

This work presents a new ground state potential energy surface (PES) for CH. The potential is tested using quasi classical trajectory (QCT) and quantum reactive scattering methods for the H + CH(+) reaction. Cross sections and rate coefficients for all reaction channels up to 300 K are calculated. The abstraction rate coefficients follow the expected slightly decreasing behaviour above 90 K, but have a positive gradient with lower temperatures. The inelastic collision and exchange reaction rate constants are increasing monotonically with temperature. The rate coefficients of the exchange reaction differ significantly between QCT and quantum reactive scattering, due to intrinsic shortcomings of the QCT final state distributions.     

Theoretical investigation of the hydrogen abstraction reaction of the OH radical with CH2FCH2F (HFC-152): a dual-level direct dynamics study

The hydrogen abstraction reaction of the OH radical with CH(2)FCH(2)F (HFC-152) is studied theoretically over the 150-3000 K temperature range. In this study, the two most recently developed hybrid density functional theories, namely, BB1K and MPWB1K, are applied, and their efficiency in reaction dynamics calculation is discussed. The BB1K/6-31+G(d,p) method gives the best result for the potential energy surface (PES) calculations, including barrier heights, reaction path information (the first and second derivatives of PES), geometry of transition state structures, and even weak hydrogen bond orientations. The rate constants were obtained by the dual-level direct dynamics with the interpolated single-point energy method (VTST-ISPE) using the BB1K/MG3S//BB1K/6-31+G(d,p) quantum model. The canonical variational transition state theory (CVT) with the small-curvature tunneling correction methods are used to calculate the rate constants in comparison to the experimental data. The total rate constant and its temperature dependency in the form of a fitted three-parameter Arrhenius expression is k(T) = 5.4 x 10(-13)(T/298)3.13 exp{-322/T} cm3 molecule(-1) s(-1). A significant variational effect, which is not common generally for hydrogen-transfer reactions, is reported and analyzed.     

H原子与CH3NH2抽氢反应的动力学计算

用QC ISD(T)/6-311 G(3DF,3PD)/MP2/6-311G(D,P)方法研究了H原子与CH3NH2的抽氢反应过程。该反应包含两个反应通道:H分别从CH3基团(R1)和NH2(R2)基团上抽氢。R1势垒比R2势垒低3.42kJ/mol,表明R1是主反应通道。在从头算的基础上,用变分过渡态理论(CVT)加小曲率隧道效应(SCT)研究了各反应温度范围为200~4000K内的速率常数,所得结果与实验值符合的很好。动力计算表明,在所研究的温度范围内,变分效应对速率常数的计算影响不大,而在低温范围内,隧道效应起了很重要的作用。