Quasi-classical trajectory study of the dynamics of the Cl + CH4→ HCl + CH3 reaction

We present an on-the-fly classical trajectory study of the Cl + CH(4)→ HCl + CH(3) reaction using a specific reaction parameter (SRP) AM1 Hamiltonian that was previously optimized for the Cl + ethane reaction [S. J. Greaves et al., J. Phys Chem A, 2008, 112, 9387]. The SRP-AM1 Hamiltonian is shown to be a good model for the potential energy surface of the title reaction. Calculated differential cross sections, obtained from trajectories propagated with the SRP-AM1 Hamiltonian compare favourably with experimental results for this system. Analysis of the vibrational modes of the methyl radical shows different scattering distributions for ground and vibrationally excited products.     

Dynamics of the O(3P) + CH4 → OH + CH3 reaction is similar to that of a triatomic reaction

The O((3)P) + CH(4) reaction has been investigated using the quasi-classical trajectory (QCT) method and an ab initio pseudotriatomic potential energy surface (PES). This has been mainly motivated by very recent experiments which support the reliability of the triatomic modeling even at high collision energy ( = 64 kcal mol(-1)). The QCT results agree rather well with the experiments (translational and angular distributions of products); i.e., the ab initio pseudotriatomic modeling "captures" the essence of the reaction dynamics, although the PES was not optimized for high E(col). Furthermore, similar experiments on the O((3)P) + CD(4) reaction at moderate E(col) (12.49 kcal mol(-1)) have also been of a large interest here and, under these softer reaction conditions, the QCT method leads to results which are almost in quantitative agreement with experiments. The utility of the ab initio pseudotriatomic modeling has also been recognized for other analogous systems (X + CH(4)) but with very different PESs.     

Ab initio potential energy surface and quantum dynamics for the H + CH4 → H2 + CH3 reaction

A new full-dimensional potential energy surface for the title reaction has been constructed using the modified Shepard interpolation scheme. Energies and derivatives were calculated using the UCCSD(T) method with aug-cc-pVTZ and 6-311++G(3df,2pd) basis sets, respectively. A total number of 30,000 data points were selected from a huge number of molecular configurations sampled by trajectory method. Quantum dynamical calculations showed that the potential energy surface is well converged for the number of data points for collision energy up to 2.5 eV. Total reaction probabilities and integral cross sections were calculated on the present surface, as well as on the ZBB3 and EG-2008 surfaces for the title reaction. Satisfactory agreements were achieved between the present and the ZBB3 potential energy surfaces, indicating we are approaching the final stage to obtain a global potential energy surface of quantitative accuracy for this benchmark polyatomic system. Our calculations also showed that the EG-2008 surface is less accurate than the present and ZBB3 surfaces, particularly in high energy region.     

Bimolecular reaction rates from ring polymer molecular dynamics: application to H + CH4 → H2 + CH3

In a recent paper, we have developed an efficient implementation of the ring polymer molecular dynamics (RPMD) method for calculating bimolecular chemical reaction rates in the gas phase, and illustrated it with applications to some benchmark atom-diatom reactions. In this paper, we show that the same methodology can readily be used to treat more complex polyatomic reactions in their full dimensionality, such as the hydrogen abstraction reaction from methane, H + CH(4) → H(2) + CH(3). The present calculations were carried out using a modified and recalibrated version of the Jordan-Gilbert potential energy surface. The thermal rate coefficients obtained between 200 and 2000 K are presented and compared with previous results for the same potential energy surface. Throughout the temperature range that is available for comparison, the RPMD approximation gives better agreement with accurate quantum mechanical (multiconfigurational time-dependent Hartree) calculations than do either the centroid density version of quantum transition state theory (QTST) or the quantum instanton (QI) model. The RPMD rate coefficients are within a factor of 2 of the exact quantum mechanical rate coefficients at temperatures in the deep tunneling regime. These results indicate that our previous assessment of the accuracy of the RPMD approximation for atom-diatom reactions remains valid for more complex polyatomic reactions. They also suggest that the sensitivity of the QTST and QI rate coefficients to the choice of the transition state dividing surface becomes more of an issue as the dimensionality of the reaction increases.     

Mechanistic and kinetic study of CF3CH═CH2 + OH reaction

The potential energy surfaces of the CF(3)CH═CH(2) + OH reaction have been investigated at the BMC-CCSD level based on the geometric parameters optimized at the MP2/6-311++G(d,p) level. Various possible H (or F)-abstraction and addition/elimination pathways are considered. Temperature- and pressure-dependent rate constants have been determined using Rice-Ramsperger-Kassel-Marcus theory with tunneling correction. It is shown that IM1 (CF(3)CHCH(2)OH) and IM2 (CF(3)CHOHCH(2)) formed by collisional stabilization are major products at 100 Torr pressure of Ar and in the temperature range of T < 700 K (at P = 700 Torr with N(2) as bath gas, T ≤ 900 K), whereas CH(2)═CHOH and CF(3) produced by the addition/elimination pathway are the dominant end products at 700-2000 K. The production of CF(3)CHCH and CF(3)CCH(2) produced by hydrogen abstractions become important at T ≥ 2000 K. The calculated results are in good agreement with available experimental data. The present theoretical study is helpful for the understanding the characteristics of the reaction of CF(3)CH═CH(2) + OH.     

A transition state view on reactive scattering: initial state-selected reaction probabilities for the H + CH4 → H2 + CH3 reaction studied in full dimensionality

Initial state-selected reaction probabilities for the H+CH(4)→H(2)+CH(3) reaction are computed for vanishing total angular momentum by full-dimensional calculations employing the multiconfigurational time-dependent Hartree approach. An ensemble of wave packets completely describing reactivity for total energies up to 0.58 eV is constructed in the transition state region by diagonalization of the thermal flux operator. These wave packets are then propagated into the reactant asymptotic region to obtain the initial state-selected reaction probabilities. Reaction probabilities for reactants in all rotational states of the vibrational 1A(1), 1F(2), and 1E levels of methane are presented. Vibrational excitation is found to decrease reactivity when reaction probabilities at equivalent total energies are compared but to increase reaction probabilities when the comparison is done at the basis of equivalent collision energies. Only a fraction of the initial vibrational energy can be utilized to promote the reaction. The effect of rotational excitation on the reactivity differs depending on the initial vibrational state of methane. For the 1A(1) and 1F(2) vibrational states of methane, rotational excitation decreases the reaction probability even when comparing reaction probabilities at equivalent collision energies. In contrast, rotational energy is even more efficient than translational energy in increasing the reaction probability when the reaction starts from the 1E vibrational state of methane. All findings can be explained employing a transition state based interpretation of the reaction process.     

Direct ab initio dynamics study of radical C4H (X̃2Σ+) + CH4 reaction

The methane (CH(4)) hydrogen abstraction reaction by linear butadiynyl radical C(4)H (CCCCH) has been investigated by direct ab initio dynamics over a wide temperature range of 100-3000 K, theoretically. The potential energy surfaces (PESs) have been constructed at the CCSD(T)/aug-cc-pVTZ//BB1K/6-311G(d,p) levels of theory. Two different hydrogen abstraction channels by C(1) and C(4) of C(4)H (C(1)C(2)C(3)C(4)H) have been considered. The results indicate that the C(1) position of C(4)H is a more reactive site. The electron transfer behaviors of two possible channels are also analyzed by quasi-restricted orbital (QRO) in detail. The rate constants calculated by canonical variational transition-state theory (CVT) with the small-curvature tunneling correction (SCT) are in excellent agreement with available experimental values. The normal and three-parameter expressions of Arrhenius rate constants are also provided within 100-3000 K. It is expected to be helpful for further studies on the reaction dynamics behaviors over a wide temperature range where no experimental data is available so far.     

Semiempirical calculation of the CH4 + CH3 · → CH3 · + CH4 radical-reaction potential surfaces in a basis of frontier orbitals

A semiempirical approach is suggested to describe potential-energy surfaces (PESs) of some radical reactions in the zero-differential overlap (ZDO) approximation. An incomplete basis set is used including only frontier (single-filled) radical molecular orbitals (MOs) depending on geometrical parameters. All possible configurations are taken into account. The parameter selection techniques are analyzed. The approach is applied to PES calculation of the CH₄ + CH₃ · CH₃ · + CH₄ reaction. Some points of the PES are verified by a nonempirical method using the perturbation theory and taking into account the correlation energy. The relaxation energies are calculated. The one-center parameters are determined nonempirically from the CH₃ + CH₃ ·, and CH₃ ^– energetics. The two-center parameters are found by modeling the CH₄ CH₃ · + H and C₂H₆ 2CH₃ reactions in the same single-orbital approximation. The energy parameters of the reactions considered are overestimated by 10%, whereas the geometrical parameters are under-estimated by 15%. Further, a comparative analysis of the Hartree-Fock solutions and those including correlation interactions (CIs) is given. The variations in the spin and charge densities on the reaction centers are considered.Institute of Chemical Kinetics and Combustion, Academy of Sciences of the USSR, Siberian Branch, Novosibirsk. Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 27, No. 4, pp. 499–506, July–August, 1991. Original article submitted February 10, 1989.     

CF3CH3 --> HF + CF2CH2: a non-RRKM reaction?

The experimental shock tube data recently reported by Kiefer et al. [J. Phys. Chem. A 2004, 108, 2443-2450] for the title reaction at temperatures between 1600 and 2400 K have been compared to master equation simulations using three models: (a) standard RRKM theory, (b) RRKM theory modified by local random matrix theory, which introduces dynamical corrections arising from slow intramolecular vibrational energy randomization, and (c) an ad hoc empirical non-RRKM model. Only the third model provides a good fit of the Kiefer et al. unimolecular reaction rate data. In separate simulations, all three models accurately reproduce the experimental 300 K chemical activation data of Marcoux and Setser [J. Phys. Chem. 1978, 82, 97-108] when the energy transfer parameters are freely varied to fit the data. When experimental energy transfer parameters for a geometrical isomer (1,1,2-trifluoroethane) are used, the standard RRKM model fits the chemical activation data better than the other models, but if energy transfer in the 1,1,1-trifluoroethane is significantly reduced in comparison to the 1,1,2 isomer, then the empirical ad hoc non-RRKM model also gives a good fit. While the ad hoc empirical non-RRKM model can be made to fit the data, it is not based on theory, and we argue that it is physically unrealistic. We also show that the master equation simulations can mimic the Kiefer et al. vibrational relaxation data, which was the first shock tube observation of double-exponential relaxation. We conclude that, until more data on the trifluoroethanes become available, the current evidence is insufficient to decide with confidence whether non-RRKM effects are important in this reaction, or whether the Kiefer et al. data can be explained in some other way.     

Quantum dynamics of heavy light heavy reactions: Application to (F + CH4 → FCH3 + H) reaction

We apply the time-dependent theory to the collinear exchange reaction F + CH₄ → FCH₃ + H. We have performed detailed calculations on two-dimensional potential surfaces representing the ground electronic potentials of the collinear F + CH₄ → FCH₃ + H reaction, at incident energies. Transmission coefficients range from zero to unity, depending upon the incident energy. Normal modes of vibrations are displayed along the reaction path.     

MCSCF and VTST studies of reaction CH + H_2→CH_3

The reaction path of CH+H2→CH3 is traced with Fukui's IRC theory and MCSCF/6-31G method. On this basis. the dynamical properties along the reaction path and variational transition state theory, (VTST. including CVT, ICVT,μVT and US methods) rate constants are investigated. The results show that the recrossing effect is small, but the curvature effect on rate constant is considerable and it must be taken into account in theoretical study.     

Crossed-beams studies of the dynamics of the H-atom abstraction reaction, O(3P) + CH4 → OH + CH3, at hyperthermal collision energies

The H-atom abstraction reaction, O((3)P) + CH(4) → OH + CH(3), has been studied at a hyperthermal collision energy of 64 kcal mol(-1) by two crossed-molecular-beams techniques. The OH products were detected with a rotatable mass spectrometer employing electron-impact ionization, and the CH(3) products were detected with the combination of resonance-enhanced multiphoton ionization (REMPI) and time-sliced ion velocity-map imaging. The OH products are mainly formed through a stripping mechanism, in which the reagent O atom approaches the CH(4) molecule at large impact parameters and the OH product is scattered in the forward direction: roughly the same direction as the reagent O atoms. Most of the available energy is partitioned into product translation. The dominance of the stripping mechanism is a unique feature of such H-atom abstraction reactions at hyperthermal collision energies. In the hyperthermal reaction of O((3)P) with CH(4), the H-atom abstraction reaction pathway accounts for 70% of the reactive collisions, while the H-atom elimination pathway to produce OCH(3) + H accounts for the other 30%.     

CH3 CH2 +N(4S)反应机理的理论研究

The reaction for CH3CH2+N(4S) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single point calculations for all the stationary points were carried out at the QCISD(T)/ 6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2CH2+3NH and H2CN＋CH3, and the minor products are the CH3CHN+H in the reaction. The majority of the products CH2CH2+3NH are formed via a direct hydrogen abstraction channel. The products H2CN＋CH3 are produced via an addition/dissociation channel. The products CH3CHN+H are produced via an addition/dissociation channel.     

Reduced dimensionality spin-orbit dynamics of CH3 + HCl ⇌ CH4 + Cl on ab initio surfaces

A reduced dimensionality quantum scattering method is extended to the study of spin-orbit nonadiabatic transitions in the CH(3) + HCl ? CH(4) + Cl((2)P(J)) reaction. Three two-dimensional potential energy surfaces are developed by fitting a 29 parameter double-Morse function to CCSD(T)/IB//MP2/cc-pV(T+d)Z-dk ab initio data; interaction between surfaces is described by geometry-dependent spin-orbit coupling functions fit to MCSCF/cc-pV(T+d)Z-dk ab initio data. Spectator modes are treated adiabatically via inclusion of curvilinear projected frequencies. The total scattering wave function is expanded in a vibronic basis set and close-coupled equations are solved via R-matrix propagation. Ground state thermal rate constants for forward and reverse reactions agree well with experiment. Multi-surface reaction probabilities, integral cross sections, and initial-state selected branching ratios all highlight the importance of vibrational energy in mediating nonadiabatic transition. Electronically excited state dynamics are seen to play a small but significant role as consistent with experimental conclusions.     

The dissociation/recombination reaction CH4 (+M) ⇔ CH3 + H (+M): a case study for unimolecular rate theory

The dissociation/recombination reaction CH(4) (+M) ? CH(3) + H (+M) is modeled by statistical unimolecular rate theory completely based on dynamical information using ab initio potentials. The results are compared with experimental data. Minor discrepancies are removed by fine-tuning theoretical energy transfer data. The treatment accounts for transitional mode dynamics, adequate centrifugal barriers, anharmonicity of vibrational densities of states, weak collision and other effects, thus being "complete" from a theoretical point of view. Equilibrium constants between 300 and 5000 K are expressed as K(c) = k(rec)/k(dis) = exp(52,044 K/T) [10(-24.65) (T/300 K)(-1.76) + 10(-26.38) (T/300 K)(0.67)] cm(3) molecule(-1), high pressure recombination rate constants between 130 and 3000 K as k(rec,∞) = 3.34 × 10(-10) (T/300 K)(0.186) exp(-T/25,200 K) cm(3) molecule(-1) s(-1). Low pressure recombination rate constants for M = Ar are represented by k(rec,0) = [Ar] 10(-26.19) exp[-(T/21.22 K)(0.5)] cm(6) molecule(-2) s(-1), for M = N(2) by k(rec,0) = [N(2)] 10(-26.04) exp[-(T/21.91 K)(0.5)] cm(6) molecule(-2) s(-1) between 100 and 5000 K. Weak collision falloff curves are approximated by asymmetric broadening factors [J. Troe and V. G. Ushakov, J. Chem. Phys. 135, 054304 (2011)] with center broadening factors of F(c) ≈ 0.262 + [(T - 2950 K)/6100 K](2) for M = Ar. Expressions for other bath gases can also be obtained.     

Kinetic isotope effects for Cl + CH4 ⇌ HCl + CH3 calculated using ab initio semiclassical transition state theory

Calculations were carried out for 25 isotopologues of the title reaction for various combinations of (35)Cl, (37)Cl, (12)C, (13)C, (14)C, H, and D. The computed rate constants are based on harmonic vibrational frequencies calculated at the CCSD(T)/aug-cc-pVTZ level of theory and X(ij) vibrational anharmonicity coefficients calculated at the CCSD(T) /aug-cc-pVDZ level of theory. For some reactions, anharmonicity coefficients were also computed at the CCSD(T)/aug-cc-pVTZ level of theory. The classical reaction barrier was taken from Eskola et al. [J. Phys. Chem. A 2008, 112, 7391-7401], who extrapolated CCSD(T) calculations to the complete basis set limit. Rate constants were calculated for temperatures from ～100 to ～2000 K. The computed ab initio rate constant for the normal isotopologue is in good agreement with experiments over the entire temperature range (～10% lower than the recommended experimental value at 298 K). The ab initio H/D kinetic isotope effects (KIEs) for CH(3)D, CH(2)D(2), CHD(3), and CD(4) are in very good agreement with literature experimental data. The ab initio (12)C/(13)C KIE is in error by ～2% at 298 K for calculations using X(ij) coefficients computed with the aug-cc-pVDZ basis set, but the error is reduced to ～1% when X(ij) coefficients computed with the larger aug-cc-pVTZ basis set are used. Systematic improvements appear to be possible. The present SCTST results are found to be more accurate than those from other theoretical calculations. Overall, this is a very promising method for computing ab initio kinetic isotope effects.     

Capture and dissociation in the complex-forming CH + H2 → CH2 + H, CH + H2 reactions

The rate coefficients for the capture process CH + H(2)→ CH(3) and the reactions CH + H(2)→ CH(2) + H (abstraction), CH + H(2) (exchange) have been calculated in the 200-800 K temperature range, using the quasiclassical trajectory (QCT) method and the most recent global potential energy surface. The reactions, which are of interest in combustion and in astrochemistry, proceed via the formation of long-lived CH(3) collision complexes, and the three H atoms become equivalent. QCT rate coefficients for capture are in quite good agreement with experiments. However, an important zero point energy (ZPE) leakage problem occurs in the QCT calculations for the abstraction, exchange and inelastic exit channels. To account for this issue, a pragmatic but accurate approach has been applied, leading to a good agreement with experimental abstraction rate coefficients. Exchange rate coefficients have also been calculated using this approach. Finally, calculations employing QCT capture/phase space theory (PST) models have been carried out, leading to similar values for the abstraction rate coefficients as the QCT and previous quantum mechanical capture/PST methods. This suggests that QCT capture/PST models are a good alternative to the QCT method for this and similar systems.     

CH4+O(3P)→CH3+OH反应的准经典轨线研究

用准经典轨线方法研究了O（3P）与CH4的反应，计算结果表明，CH4（υ＝0，j＝0）与O（3P）的反应在低及高的碰撞参数下都是直接反应，无短寿命的碰撞复合物生成，产物OH以向后散射为主，基本上处于振转基态.CH4（υ＝1，j＝1）与O（3P）的反应在低及高的碰撞参数下反应机理不一样。在低碰撞参数下是直接反应，无短寿命的碰撞复合物生成，产物OH以向后散射为主，主要处于振动基态，转动基本上是冷的，但比高碰撞参数下的热.在高的碰撞参数下则生成短寿命的碰撞复合物，产物OH以向前散射为主，表现出明显的周边动力学反应的特征，主要处于振动激发态（υ＝1），但转动仍然是较冷的。     

Electronic state dependence of the ion-molecule reaction CH(3)CN(+) + CH(3)CN → CH(4)CN(+) + CH(2)CN: threshold electron-secondary ion coincidence (TESICO) and direct ab initio molecular dynamics study

The ion-molecule reaction, CH(3)CN(+) + CH(3)CN → CH(3)CNH(+) + CH(2)CN, has been investigated using the threshold electron-secondary ion coincidence (TESICO) technique. Relative reaction cross sections for two microscopic reaction mechanisms, i.e., proton transfer (PT) from the acetonitrile ion CH(3)CN(+) to neutral acetonitrile CH(3)CN and hydrogen atom abstraction (HA) by CH(3)CN(+) from CH(3)CN, have been determined for two low-lying electronic states, (2)E and (2)A(1) of the CH(3)CN(+) primary ion. The cross section for PT of the (2)A(1) state was smaller than that of the (2)E state, whereas that of HA are almost the same in the two states. Ab initio calculations showed that the dissociation of the C-H(+) bond of CH(3)CN(+) is easier in the (2)E state than that in the (2)A(1) state. The direct ab initio molecular dynamics (MD) calculations showed that two mechanisms, direct proton transfer and complex formation, contribute the reaction dynamics.