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.     

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.     

Quasi-classical trajectory study of the Ne + H2(+) → NeH(+) + H reaction based on global potential energy surface

Using the multireference configuration interaction method with a Davidson correction and a large orbital basis set (aug-cc-pVQZ), we obtain an energy grid that includes 32 038 points for the construction of a new analytical potential energy surface (APES) for the Ne + H(2)(+) → NeH(+) + H reaction. The APES is represented as a many-body expansion containing 142 parameters, which are fitted from 31?000 ab initio energies using an adaptive nonlinear least-squares algorithm. The geometric characteristics of the reported APES and the one presented here are also compared. On the basis of the APES we obtained, reaction cross sections are computed by means of quasi-classical trajectory (QCT) calculations and compared with the experimental and theoretical data in the literature.     

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.     

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.     

Theoretical studies on the reaction path and dynamics of the reaction CH_3+H_2O→CH_4+OH

The reaction path, the dynamical properties along the reaction path and CVT rate constants are computed by the ab initio MO method, the reaction path Hamiltonian theory and the variational transition state theory. The results show that the effect of the electron correlation energy on activation barrier is large, the recrossing and tunneling effects exist in the reaction.     

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.     

Theoretical Studies on CH3SiH3+H→CH3SiH2+H2 Reaction with the Variational Transitional State Theory

The abstract reaction of CH3SiH3 with H has been studied by using the“direct dynamics”method of variational transition-state theory,which is based on the information on geometries,frequencies and energies calculated by ab initio along the minimum energy path.The rate constants and transmission coefficients were calculated for the temperature range 298～1700K.The result indicates that the variational effect on this reaction is great and the tunneling effect is very obvious at room temperature.The rate constants calculated match well with the experimental value.     

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.     

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...     

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.     

Quantum dynamics calculation of reaction probability for H + Cl2 → HCl + Cl

We present in this paper a time-dependent quantum wave packet calculation of the initial state selected reaction probability for H + Cl₂ based on the GHNS potential energy surface with total angular momentumJ = 0. The effects of the translational, vibrational and rotational excitation of Cl₂ on the reaction probability have been investigated. In a broad region of the translational energy, the rotational excitation enhances the reaction probability while the vibrational excitation depresses the reaction probability. The theoretical results agree well with the fact that it is an early down-hill reaction.     

Lack of salt effect on the kinetics of the reaction SO2−4 + H3O+ → HSO−4 + H2O

It is found that the broadening of the 1100-cm⁻¹ line of SO⁻²₄, caused by increasing [H₃O⁺], is unaffected by addition of 4 M LiCl, NaBr, KCl and NH₄Cl. This finding is in line with the lack of influence of NaCl reported earlier. The significance of these findings, in terms of the reaction mechanism, is discussed.     

Comment on the possible role of the reaction H+H2O→H2+OH in the radiolysis of water at high temperatures

In a recent paper (Radiation Physics and Chemistry, 2005, vol. 74, pp. 210) it was suggested that the anomalous increase of molecular hydrogen radiolysis yields observed in high-temperature water is explained by a high activation energy for the reaction H+H₂O→H₂+OH. In this comment we present thermodynamic arguments to demonstrate that this reaction cannot be as fast as suggested. A best estimate for the rate constant is 2.2×10³ M⁻¹ s⁻¹ at 300 °C. Central to this argument is an estimate of the OH radical hydration free energy vs. temperature, ΔG_hyd(OH)=0.0278t−18.4 kJ/mole (t in °C, equidensity standard states), which is based on analogy with the hydration free energy of water and of hydrogen peroxide.     

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.     

Quantum dynamics calculation of reaction probability for H+Cl_2→HC1+Cl

We present in this paper a time-dependent quantum wave packet calculation of the initial state selected reaction probability for H + CI2 based on the GHNS potential energy surface with total angular momentum J= 0. The effects of the translational, vibrational and rotational excitation of CI2 on the reaction probability have been investigated. In a broad region of the translational energy, the rotational excitation enhances the reaction probability while the vibrational excitation depresses the reaction probability. The theoretical results agree well with the fact that it is an early down-hill reaction.     

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%.     

Quantum instanton calculation of rate constants for the C2H6 + H → C2H5 + H2 reaction: anharmonicity and kinetic isotope effects

Thermal rate constants and kinetic isotope effects for the title reaction are calculated by using the quantum instanton approximation within the full dimensional Cartesian coordinates. The obtained results are in good agreement with experimental measurements at high temperatures. The detailed investigation reveals that the anharmonicity of the hindered internal rotation motion does not influence the rate too much compared to its harmonic oscillator approximation. However, the motion of the nonreactive methyl group in C(2)H(6) significantly enhances the rates compared to its rigid case, which makes conventional reduced-dimensionality calculations a challenge. In addition, the temperature dependence of kinetic isotope effects is also revealed.     

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.