Potential energy surface and reaction mechanism for the ion-molecule reaction: CH4 + CH4+ → CH3 + CH5+

The potential energy surface for the CH₄+CH₄⁺ reaction system has been calculated with the ab initio method. A stable complex, responsible for the complex mechanism, has been found but is hard to reach. Each of the two direct mechanisms, hydrogen transfer and proton transfer, has been shown to consist of a combination of electron transfer and hydrogen atom transfer processes.     

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.     

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.     

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.     

Ten-Dimensional Quantum Dynamics Study of H+CH3D→H2+CH2D Reaction

The mode selectivity of the H+CH3D→H2+CH2D reaction was studied using a recently developed ten-dimensional time-dependent wave packet method.The reac-tion dynam...     

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.     

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.     

Quantum dynamics of the S+OH→SO+H reaction

First accurate quantum mechanical scattering calculations have been carried out for the S((3)P)+OH(X?(2)Π)→SO(X?(3)Σ(-))+H((2)S) reaction using a recent ab initio potential energy surface for the ground electronic state, X?(2)A("), of HSO. Total and state-to-state reaction probabilities for a total angular momentum J=0 have been determined for collision energies up to 0.5 eV. A rate constant has been calculated by means of the J-shifting approach in the 10-400 K temperature range. Vibrational and rotational product distributions show no specific behavior and are consistent with a mixture of direct and indirect reaction mechanisms.     

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.     

H+CH3NO2H2+CH2NO2反应途径和变分速率常数计算研究

采用MP2(FULL)/6-311G**从头算方法, 优化了H+CH3NO2H2+ CH2NO2反应的过渡态结构, 得出该反应的正逆反应的活化位垒分别是82.73和57.14 kJ*mol-1. 沿IRC分析指出该反应是一个H-H键生成和C-H键断裂的协同反应, 而且在反应途径上存在一个引导反应进行的振动模式, 这一反应模式引导反应进行的区间在-0.7～0.2( amu)1/2*a0之间; 在1 000～1 400 K温度范围内, 运用变分过渡态理论(CVT), 计算了该反应的速率常数, 计算结果与实验相一致.     

用变分过渡态理论研究CH3SiH3+H→CH3SiH2+H2反应动力学

用变分过渡态理论对CH3SiH3与H的抽提反应进行了理论研究;利用从头算计算了反应体系的构型、振动频率和能量等信息;计算了温度在298 ～1700K内反应的速率常数和穿透系数。结果表明,在室温下,变分对于此反应影响较大,隧道效应特别明显,计算得到的速率常数和实验值符合得很好。     

The TST-CEQ calculations on the light-heavy-light H+ClH reaction

The calculated TST-CEQ state-selected reaction cross sections and rate constants forthe exchange reaction of H+ClH are reported.It is found that,although the collinear proba-bilities exhibit the oscillating behaviour,the TST-CEQ cross sections do not oscillate.Compa-rison of the different state-selected cross sections at a constant total energy indicates thattranslational energy is more effective than vibrational energy in promoting the reaction H+ClH.We also find that,the larger the vibrational quantum number of the reactant,the larger thereactive rate constant k~(TST-CEQ)(T,v)is at a given temperature.The state-selected rate constantsare then averaged to yield the thermal rate constants which are shown to be in good agreementwith our VTST calculations and those by Schwenke et al..     

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.     

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.     

Dynamics and mode-selected reaction of NCO+H_2→HNCO+H

The reaction path of the reaction NCO+H2→HNCO + H has been traced by Fukui's theory and the ab initio method. On this basis, the dynamical properties along the reaction path, canonical variational theory (CVT) rate constants and vibrational-mode-selected rate constants have been computed. The results show that the effect of the electron correlation energy on the activation barrier is large, and tiros the correction by MP4 method is effective; the results also show that the recrossing and tunneling effects exist, and thus the corrections by the variational transition state theory (VTST) and the small curvature (SC) approximation method are also effective. In the reaction, the coupling and energy transfer between mode 8(7) and reaction path are strong, so the rate is effectively enhanced while these modes, especially H2 stretching, are vibrationally excited.     

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.     

Mean potential statistical theory of the H+ + D2 → HD + D+ reaction

The reactive collision process H(+) + D(2)(ν = 0, j = 0) → HD + D(+) is theoretically analyzed for collision energies ranging from threshold up to 1.3 eV. It is assumed that the reaction takes place via formation of a collision complex. In calculations, a statistical theory is used, based on a mean isotropic potential deduced from a full potential energy surface. Calculated integral cross sections, opacity functions, and rotational distributions of the HD products are compared with recent statistical and quantum mechanical calculations performed using a full potential energy surface. Satisfactory agreement between the results obtained using the two statistical methods is found, both of which however overestimate the existing quantum mechanical predictions. The effects due to the presence of identical particles are also discussed.