Stereodynamics and rovibrational effect for H+CH4(v,j,K,n)→H2+CH3 reaction

Full-dimensional quantum mechanical calculations of the reaction probability of the H + CH4 reaction based on a mixed Jacobi and Radau description

A full-dimensional time-dependent wave packet study using mixed polyspherical Jacobi and Radau coordinates for the title reaction has been reported. The non-reactive moiety CH₃ has been described using three Radau vectors, whereas two Jacobi vectors have been used for the bond breaking/formation process. A potential-optimized discrete variable representation basis has been employed to describe the vibrational coordinates of the reagent CH₄. About one hundred billion basis functions have been necessary to achieve converged results. The reaction probabilities for some initial vibrational states are given. A comparison between the present approach and other methods, including reduced and full-dimensional ones, is also presented.

Testing the Palma-Clary Reduced Dimensionality Model Using Classical Mechanics on the CH4 + H → CH3 + H2 Reaction

Application of exact quantum scattering methods in theoretical reaction dynamics of bimolecular reactions is limited by the complexity of the equations of nuclear motion to be solved. Simplification is often achieved by reducing the number of degrees of freedom to be explicitly handled by freezing the less important spectator modes. The reaction cross sections obtained in reduced-dimensionality (RD) quantum scattering methods can be used in the calculation of rate coefficients, but their physical meaning is limited. The accurate test of the performance of a reduced-dimensionality method would be a comparison of the RD cross sections with those obtained in accurate full-dimensional (FD) calculations, which is not feasible because of the lack of complete full-dimensional results. However, classical mechanics allows one to perform reaction dynamics calculations using both the RD and the FD model. In this paper, an RD versus FD comparison is made for the 8-dimensional Palma-Clary model on the example of four isotopologs of the CH4 + H → CH3 + H2 reaction, which has 12 internal dimensions. In the Palma-Clary model, the only restriction is that the methyl group is confined to maintain C3v symmetry. Both RD and FD opacity and excitation functions as well as differential cross sections were calculated using the quasiclassical trajectory method. The initial reactant separation has been handled according to our one-period averaging method [ Nagy et al. J. Chem. Phys. 2016, 144, 014104 ]. The RD and FD excitation functions were found to be close to each other for some isotopologs, but in general, the RD reactivity parameters are lower than the FD reactivity parameters beyond statistical error, and for one of the isotopologs, the deviation is significant. This indicates that the goodness of RD cross sections cannot be taken for granted.

Control of chemical reactivity by transition-state and beyond

It has been long established that the transition state for an activated reaction controls the overall reactivity, serving as the bottleneck for reaction flux. However, the role of the transition state in regulating quantum state resolved reactivity has only been addressed more recently, thanks to advances in both experimental and theoretical techniques. In this perspective, we discuss some recent advances in understanding mode-specific reaction dynamics in bimolecular reactions, mainly focusing on the X + H2O/CH4 (X = H, F, Cl, and O(3P)) systems, extensively studied in our groups. These advances shed valuable light on the importance of the transition state in mode-specific and steric dynamics of these prototypical reactions. It is shown that many mode-specific phenomena can be understood in terms of a transition-state based model, which assumes in the sudden limit that the ability of a reactant mode for promoting the reaction stems from its coupling with the reaction coordinate at the transition state. Yet, in some cases the long-range anisotropic interactions in the entrance (or exit) valley, which govern how the trajectories reach (or leave) the transition state, also come into play, thus modifying the reactive outcomes.

QUANTITATIVE QUANTUM DYNAMICS CALCULATION OF H2 + CH3 → H + CH4 REACTION

We report in this paper quantum dynamics calculation of state-selected reaction probabilities for a benchmark chemical reaction H 2 + CH 3 → H + CH 4 on an ab initio potential energy surface. The quantum dynamics calculation is based on the recently developed semirigid vibrating rotor target (SVRT) model and involves six degrees of freedom. The present result is the first such high-level quantum dynamics calculation of microscopic reaction probability for a chemical reaction between two molecules with at least one of the reagents being larger than a diatomic molecule.

Multidimensional molecular steric opacity function for XeCl*(B, C) formation in the oriented Xe* (³P₂, MJ = 2) + oriented CCl₃F reaction

The steric effect for the XeCl*(B, C) formations in the oriented Xe* (³P₂, MJ = 2) + oriented CCl₃F reaction has been observed as a function of the mutual configuration between the molecular orientation and the atomic orbital alignment in the collision frame. Molecular steric opacity functions have been determined as a function of the atomic orbital alignment (M(L)') in the collision frame. The XeCl*(B, C) channels show similar molecular steric opacity functions at M(L)' = 0 but not at |M(L)'| = 1. The large molecular alignment dependence (i.e., the reactivity of the Cl₃ end and the F end is comparable, but a very poor reactivity for the sideway) is recognized for the XeCl*(B, C) channels except for the XeCl*(C) channel at |M(L)'| = 1, which shows an almost isotropic molecular orientation dependence. The M(L)' selectivity is different between the XeCl*(B, C) channels. At the molecular axis direction, the XeCl*(B) channel has little M(L)' selectivity whereas the XeCl*(C) channel is significantly favorable at M(L)' = 0. On the other hand, |M(L)'| = 1 is favorable at the sideway for the XeCl*(B, C) channels.

Quantum dynamics study of H+NH3→H2+NH2 reaction

We report in this paper a quantum dynamics study for the reaction H+NH3-->NH2+H2 on the potential energy surface of Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The quantum dynamics calculation employs the semirigid vibrating rotor target model [J. Z. H. Zhang, J. Chem. Phys. 111, 3929 (1999)] and time-dependent wave packet method to propagate the wave function. Initial state-specific reaction probabilities are obtained, and an energy correction scheme is employed to account for zero point energy changes for the neglected degrees of freedom in the dynamics treatment. Tunneling effect is observed in the energy dependency of reaction probability, similar to those found in H+CH4 reaction. The influence of rovibrational excitation on reaction probability and stereodynamical effect are investigated. Reaction rate constants from the initial ground state are calculated and are compared to those from the transition state theory and experimental measurement.

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.

Accurate quantum calculations of the reaction rates for H∕D+CH4

In previous work [T. Wu, H. J. Werner, and U. Manthe, Science 306, 2227 (2004)], accurate quantum reaction rate calculations of the rate constant for the H+CH4-->CH3+H2 reaction have been presented. Both the electronic structure calculations and the nuclear dynamics calculations are converged with respect to the basis sets employed. In this paper, the authors apply the same methodology to an isotopic variant of this reaction: D+CH4-->CH3+HD. Accurate rate constants are presented for temperatures between 250 and 400 K. For temperatures between 400 and 800 K, they use a harmonic extrapolation to obtain approximate rate constants for H/D+CH4. The calculations suggest that the experimentally reported rate constants for D+CH4 are about a factor of 10-20 too high. For H+CH4, more accurate experiments are available and agreement is much better: the difference is less than a factor of 2.6. The kinetic isotope effect for the H/D+CH4 reactions is studied and compared with experiment and transition state theory (TST) calculations. Harmonic TST was found to provide a good description of the kinetic isotope effect.

A full dimensional, nine-degree-of-freedom, time-dependent quantum dynamics study for the H2+C2H reaction

A full dimensional, nine-degree-of-freedom (9DOF), time-dependent quantum dynamics wave packet approach is presented for the study of the H2+C2H-->H+C2H2 reaction system. This is the first full dimensional quantum dynamics study for a diatom-triatom reaction system. The effects of the initial vibrational and rotational excitations of the reactants on the reactivity of this reaction are investigated. This study shows that vibrational excitations of H2 enhance the reactivity; whereas, the vibrational excitations of C2H only have a small effect on the reaction probability. In addition, the bending excitations of C2H, compared to the ground state reaction probability, hinder the reactivity. Comparison of the ground state reaction probabilities of the 9DOF and 8DOF shows the reaction probability from the full dimensional calculation is larger, with more prominent resonance features.

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.

An eight-degree-of-freedom quantum dynamics study for the H2+C2H system

An eight-degree-of-freedom (8DOF) time-dependent wave-packet approach has been developed to study the H(2)+C(2)H-->H+C(2)H(2) reaction system. The 8DOF model is obtained by fixing one of the Jacobi torsion angle in the nine-degree-of-freedom AB+CDE reaction system. This study is an extension of the previous seven-degree-of-freedom (7DOF) computation [J. Chem. Phys. 119, 12057 (2003)] of this reaction system. This study shows that vibrational excitations of H(2) enhance the reaction probability, whereas the stretching vibrational excitations of C(2)H have only a small effect on the reactivity. Furthermore, the bending excitation of C(2)H, compared to the ground-state reaction probability, hinders the reactivity. A comparison of the rate constant between the 7DOF calculation and the present 8DOF results has been made. The theoretical and experimental results agree with each other very well when the present 8DOF results are adjusted to account for the lower transition state barrier heights found in recent ab initio calculations.

Effects of C–H stretch excitation on the H+CH4 reaction

We have investigated the effects of C-H stretching excitation on the H+CH4-->CH3+H2 reaction dynamics using the photo-LOC technique. The CH3 product vibrational state and angular distribution are measured for the reaction of fast H atoms with methane excited in either the antisymmetric stretching fundamental (nu3=1) or first overtone (nu3=2) with a center-of-mass collision energy of Ecoll ranging from 1.52 to 2.20 eV. We find that vibrational excitation of the nu3=1 mode enhances the overall reaction cross section by a factor of 3.0+/-1.5 for Ecoll=1.52 eV, and this enhancement factor is approximately constant over the 1.52-2.20-eV collision energy range. A local-mode description of the CH4 stretching vibration, in which the C-H oscillators are uncoupled, is used to describe the observed state distributions. In this model, the interaction of the incident H atom with either a stretched or an unstretched C-H oscillator determines the vibrational state of the CH3 product. We also compare these results to the similar quantities obtained previously for the Cl+CH4-->CH3+HCl reaction at Ecoll=0.16 eV [Z. H. Kim, H. A. Bechtel, and R. N. Zare, J. Chem. Phys. 117, 3232 (2002); H. A. Bechtel, J. P. Camden, D. J. A. Brown, and R. N. Zare, ibid. 120, 5096 (2004)] in an attempt to elucidate the differences in reactivity for the same initially prepared vibration.

Calculation of converged rovibrational energies and partition function for methane using vibrational–rotational configuration interaction

The rovibration partition function of CH4 was calculated in the temperature range of 100-1000 K using well-converged energy levels that were calculated by vibrational-rotational configuration interaction using the Watson Hamiltonian for total angular momenta J = 0-50 and the MULTIMODE computer program. The configuration state functions are products of ground-state occupied and virtual modals obtained using the vibrational self-consistent field method. The Gilbert and Jordan potential energy surface was used for the calculations. The resulting partition function was used to test the harmonic oscillator approximation and the separable-rotation approximation. The harmonic oscillator, rigid-rotator approximation is in error by a factor of 2.3 at 300 K, but we also propose a separable-rotation approximation that is accurate within 2% from 100 to 1000 K.