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

Ring Polymer Molecular Dynamics Calculations of Thermal Rate Constants for the O((3)P) + CH4 → OH + CH3 Reaction: Contributions of Quantum Effects

The thermal rate constant of the O((3)P) + CH4 → OH + CH3 reaction is investigated with ring polymer molecular dynamics on a full-dimensional potential energy surface. Good agreement with experimental and full-dimensional quantum multiconfiguration time-dependent Hartree results between 300 and 1500 K was obtained. It is shown that quantum effects, for example, tunneling and zero-point energy, can be effectively and efficiently included in this path-integral based approach implemented with classical trajectories. Convergence with respect to the number of beads is rapid, suggesting wide applicability for other reactions involving polyatomic molecules.

Mode Selectivity for a "Central" Barrier Reaction: Eight-Dimensional Quantum Studies of the O((3)P) + CH4 → OH + CH3 Reaction on an Ab Initio Potential Energy Surface

The dynamics of a combustion reaction, namely, O((3)P) + CH4 → OH + CH3, is investigated with an eight-dimensional quantum model that includes representatives of all vibrational modes of CH4 and with a full-dimensional quasi-classical trajectory (QCT) method. The calculated excitation functions for the ground vibrational state CH4 agree well with experiment. Both quantum and QCT results suggest that excitation of the stretching modes of CH4 enhances the reaction, while the bending and umbrella modes have a smaller impact on reactivity, again consistent with experimental findings. However, none of the vibrational excitations has comparable efficiency in promoting the reaction as translational energy.

Theoretical Study of the Validity of the Polanyi Rules for the Late-Barrier Cl + CHD3 Reaction

The Polanyi rules, which state that vibrational energy is more efficient in promoting a late-barrier reaction than translational energy, were questioned recently by an experimental unexpected finding that the CH stretch excitation is no more effective in promoting the late-barrier Cl + CHD3 reaction than the translational energy. However, the present quantum dynamics study on the best-available potential energy surface for the title reaction reveals that the CH stretch excitation does promote the reaction significantly, except at low collision energies. Further studies should be carried out to solve the disagreements between theory and experiment on the reaction.

Dynamics of the O(3P) + CHD3(vCH = 0,1) reactions on an accurate ab initio potential energy surface

Recent experimental and theoretical studies on the dynamics of the reactions of methane with F and Cl atoms have modified our understanding of mode-selective chemical reactivity. The O + methane reaction is also an important candidate to extend our knowledge on the rules of reactivity. Here, we report a unique full-dimensional ab initio potential energy surface for the O((3)P) + methane reaction, which opens the door for accurate dynamics calculations using this surface. Quasiclassical trajectory calculations of the angular and vibrational distributions for the ground state and CH stretching excited O + CHD(3)(v(1) = 0,1) → OH + CD(3) reactions are in excellent agreement with the experiment. Our theory confirms what was proposed experimentally: The mechanistic origin of the vibrational enhancement is that the CH-stretching excitation enlarges the reactive cone of acceptance.

A quantum reaction dynamics study of the translational, vibrational, and rotational motion effects on the HD + H3+ reaction

Time-dependent, quantum reaction dynamics wavepacket approach is employed to investigate the impacts of the translational, vibrational, and rotational motion on the HD+H(3)(+) → H(2)D(+) + H(2) reaction using the Xie-Braams-Bowman potential energy surface [Z. Xie, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 122, 224307 (2005)]. We treat this five atom reaction with a seven-degree-of-freedom model by fixing one Jacobi and one torsion angle related to H(3) (+) at the lowest saddle point geometry of the potential energy surface. The initial state selected reaction probabilities show that the rotational excitations of H(+)-H(2) greatly enhance the reactivity with the reaction probabilities increased double at high rotational states compared to the ground state. However, the vibrational excitations of H(3) (+) hinder the reactivity. The ground state reaction probability shows no reaction threshold for this exoergic reaction, and as the translational energy increases, the reaction probability decreases. Furthermore, reactive resonances and zero point energy play very important roles on the reaction dynamics. The obtained integral cross section has the character of an exoergic reaction without a threshold: it decreases with the translational energy increasing. The calculated thermal rate constants using this seven-degree-of-freedom model are in agreement with a later experiment measurement.

The method of Gaussian weighted trajectories. V. On the 1GB procedure for polyatomic processes

In recent years, many chemical reactions have been studied by means of the quasiclassical trajectory (QCT) method within the Gaussian binning (GB) procedure. The latter consists of "quantizing" the final vibrational actions in Bohr spirit by putting strong emphasis on the trajectories reaching the products with vibrational actions close to integer values. A major drawback of this procedure is that if N is the number of product vibrational modes, the amount of trajectories necessary to converge the calculations is ∼10(N)×larger than with the standard QCT method. Applying it to polyatomic processes is thus problematic. In a recent paper, however, Czakó and Bowman propose to quantize the total vibrational energy instead of the vibrational actions [G. Czakó and J. M. Bowman, J. Chem. Phys. 131, 244302 (2009)], a procedure called 1GB here. The calculations are then only ∼10 times more time consuming than with the standard QCT method, allowing thereby for considerable numerical saving. In this paper, we propose some theoretical arguments supporting the 1GB procedure and check its validity on model test cases as well as the prototype four-atom reaction OH+D(2)→HOD+D.

Seven-dimensional quantum dynamics study of the H+NH3-->H2+NH2 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H+NH3-->H2+NH2 reaction using a seven-dimensional model and an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values. The total reaction probabilities are calculated for the initial ground and seven excited states of NH3 with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants are calculated for the temperature range 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the total reaction probabilities are overall very small, (b) the symmetric and asymmetric NH stretch excitations enhance the reaction significantly and almost all of the excited energy deposited was used to reduce the reaction threshold, (c) the excitation of the umbrella and bending motion have a smaller contribution to the enhancement of reactivity, (d) the main contribution to the thermal rate constants is thought to come from the ground state at low temperatures and from the stretch excited states at high temperatures, and (e) the calculated thermal rate constants are three to ten times smaller than the experimental data and transition state theory results.