Wave packet dynamics of H2(v1=8-14)+H2(v2=0-2): the role of the potential energy surface on different reactive and dissociative processes

A time-dependent wave packet method has been used to study different competing products of H(2)+H(2) collisions: four center reaction, collision induced dissociation, reactive dissociation, and three-body complex formation. A three-degree-of-freedom reduced dimensionality model has been used for five different geometries of the colliding complex (parallel H, crossed X, collinear L, and two T-shaped geometries T(I) and T(II)), with reactants in selected vibrational states with one diatom vibrationally "hot" and the other one vibrationally "cold." Product probabilities have been calculated using two potential energy surfaces [J. Chem. Phys. 101, 4004 (1994); J. Chem. Phys. 116, 666 (2002)] in order to compare their performance in the dynamics. The regions of the potential energy surfaces responsible of the threshold behavior of the probabilities have been identified. Overall, we have found that the most recent potential energy surface is less anisotropic, provides a smaller propensity for insertion-type processes, and gives lower energy thresholds.     

Probing Feshbach resonances in F+H2(j=1)-->HF+H: dynamical effect of single quantum H2-rotation

Full quantum state resolved scattering of the F atom reaction with H(2)(j=0) and H(2)(j=1) was investigated at the collision energies of 0.19 and 0.56 kcalmol. Dramatic difference between the dynamics for the F+H(2)(j=0,1) reactions at both collision energies have been observed. Forward scattering HF(v(')=2) products have been observed unambiguously for the F+H(2)(j=1) reaction at low collision energies, which was attributed to the Feshbach resonances. This study provides a unique case of reaction resonances involving a rotationally excited reagent.     

Communication: Experimental and theoretical investigations of the effects of the reactant bending excitations in the F + CHD3 reaction

The effects of the reactant bending excitations in the F+CHD(3) reaction are investigated by crossed molecular beam experiments and quasiclassical trajectory (QCT) calculations using a high-quality ab initio potential energy surface. The collision energy (E(c)) dependence of the cross sections of the F+CHD(3)(v(b)=0,1) reactions for the correlated product pairs HF(v('))+CD(3)(v(2)=0,1) and DF(v('))+CHD(2)(v(4)=0,1) is obtained. Both experiment and theory show that the bending excitation activates the reaction at low E(c) and begins to inactivate at higher E(c). The experimental F+CHD(3)(v(b)=1) excitation functions display surprising peak features, especially for the HF(v(')=3)+CD(3)(v(2)=0,1) channels, indicating reactive resonances (quantum effects), which cannot be captured by quasiclassical calculations. The reactant state-specific QCT calculations predict that the v(5)(e) bending mode excitation is the most efficient to drive the reaction and the v(6)(e) and v(5)(e) modes enhance the DF and HF channels, respectively.     

Differential cross section for the H+D(2)-->HD(v(')=1,j(')=2,6,10)+D reaction as a function of collision energy

We have measured differential cross sections (DCSs) for the HD (v(')=1,j(')=2,6,10) products of the H+D(2) exchange reaction at five different collision energies in the range 1.48< or =E(coll)< or =1.94 eV. The contribution from the less energetic H atoms formed upon spin-orbit excitation of Br in the photolysis of the HBr precursor is taken into account for two collision energies, E(coll)=1.84 and 1.94 eV, allowing us to disentangle the two different channels. The measured DCSs agree well with new time-dependent quantum-mechanical calculations. As the product rotational excitation increases, the DCSs shift from backward to sideward scattering, as expected. We also find that the shapes of the DCSs show only a small overall dependence on the collision energy, with a notable exception occurring for HD (v(')=1,j(')=2), which appears bimodal at high collision energies. We suggest that this feature results from both direct recoil and indirect scattering from the conical intersection.     

Nonadiabatic time-dependent wave packet study of the D+ + H2 reaction system

A theoretical investigation on the nonadiabatic processes of the D(+) + H(2) reaction system has been carried out by means of exact three-dimensional nonadiabatic time-dependent wave packet calculations with an extended split operator scheme (XSOS). The diabatic potential energy surface newly constructed by Kamisaka et al. (J. Chem. Phys. 2002, 116, 654) was employed in the calculations. This study provided quantum cross sections for three competing channels of the reactive charge transfer, the nonreactive charge transfer, and the reactive noncharge transfer, which contrasted markedly to many previous quantum theoretical reports on the (DH(2))(+) system restricted to the total angular momentum J = 0. These quantum theoretical cross sections derived from the ground rovibrational state of H(2) show wiggling structures and an increasing trend for both the reactive charge transfer and the nonreactive charge transfer but a decreasing trend for the reactive noncharge transfer throughout the investigated collision energy range 1.7-2.5 eV. The results also show that the channel of the reactive noncharge transfer with the largest cross section is the dominant one. A further investigation of the v-dependent behavior of the probabilities for the three channels revealed an interesting dominant trend for the reactive charge transfer and the nonreactive charge transfer at vibrational excitation v = 4 of H(2). In addition, the comparison between the centrifugal sudden (CS) and exact calculations showed the importance of the Coriolis coupling for the reactive system. The computed quantum cross sections are also compared with the experimental measurement results.     

Dynamical regimes on the Cl + H2 collisions: inelastic rainbow scattering

While Cl + H(2) reactive collisions have been a subject of numerous experimental and theoretical studies, inelastic collisions leading to rotational energy transfer and/or vibrational excitation have been largely ignored. In this work, extensive quantum mechanical calculations covering the 0.5-1.5 eV total energy range and various initial rovibrational states have been carried out and used to perform a joint study of inelastic and reactive Cl + H(2) collisions. Quasiclassical trajectories calculations complement the quantum mechanical results. The analysis of the inelastic transition probabilities has revealed the existence of two distinct dynamical regimes that correlate with low and high impact parameters, b, and are neatly separated by glory scattering. It has been found that while high-b collisions are mainly responsible for |Δj| = 2 transitions which dominate the inelastic scattering, they are very inefficient in promoting higher |Δj| transitions. The effectiveness of this type of collision also drops with rotational excitation of H(2). In contrast, reactive scattering, that competes with |Δj|?> 2 inelastic transitions, is exclusively caused by low-b collisions, and it is greatly favored when the reactants get rotationally excited. Previous studies focusing on the reactivity of the Cl + H(2) system established that the van der Waals well located in the entrance channel play a key role in determining the mechanism of the collisions. Our results prove this to be also a case for inelastic processes, where the origin of the double dynamical regime can be traced back to the influence exerted by this well that shapes the topology of the entrance channel of the Cl-H(2) system.     

Exact quantum scattering study of the H + HS reaction on a new ab initio potential energy surface H2S (3A")

We present an exact quantum dynamical study and quasi-classical trajectory (QCT) calculations for the exchange and abstraction processes for the H + HS reaction. These calculations were based on a newly constructed high-quality potential energy surface for the lowest triplet state of H(2)S ((3)A"). The ab initio single-point energies were computed using complete active space self-consistent field and multi-reference configuration interaction method with a basis set of aug-cc-pV5Z. The time-dependent wave packet (TDWP) method was used to calculate the total reaction probabilities and integral cross sections over the collision energy (E(col)) range of 0.0-2.0 eV for the reactant HS initially at the ground state and the first vibrationally excited state. It was found that the initial vibrational excitation of HS enhances both abstraction and exchange processes. In addition, a good agreement is found between QCT and TDWP reaction probabilities at the total momentum J = 0 as a function of collision energy for the H + HS (v = 0, j = 0) reaction.     

Bimolecular reactions of vibrationally excited molecules. Roaming atom mechanism at low kinetic energies

Quasiclassical trajectory calculations have been performed for the H + H'X(v) → X + HH' abstraction and H + H'X(v) → XH + H' (X = Cl, F) exchange reactions of the vibrationally excited diatomic reactant at a wide collision energy range extending to ultracold temperatures. Vibrational excitation of the reactant increases the abstraction cross sections significantly. If the vibrational excitation is larger than the height of the potential barrier for reaction, the reactive cross sections diverge at very low collision energies, similarly to capture reactions. The divergence is quenched by rotational excitation but returns if the reactant rotates fast. The thermal rate coefficients for vibrationally excited reactants are very large, approach or exceed the gas kinetic limit because of the capture-type divergence at low collision energies. The Arrhenius activation energies assume small negative values at and below room temperature, if the vibrational quantum number is larger than 1 for HCl and larger than 3 for HF. The exchange reaction also exhibits capture-type divergence, but the rate coefficients are larger. Comparisons are presented between classical and quantum mechanical results at low collision energies. At low collision energies the importance of the exchange reaction is enhanced by a roaming atom mechanism, namely, collisions leading to H atom exchange but bypassing the exchange barrier. Such collisions probably have a large role under ultracold conditions. The calculations indicate that for roaming to occur, long-range attractive interaction and small relative kinetic energy in the chemical reaction at the first encounter are necessary, which ensures that the partners can not leave the attractive well. Large orbital angular momentum of the primary products (equivalent to large rotational excitation in a unimolecular reaction) is favorable for roaming.     

Study of the H+O2 reaction by means of quantum mechanical and statistical approaches: the dynamics on two different potential energy surfaces

The possible existence of a complex-forming pathway for the H+O(2) reaction has been investigated by means of both quantum mechanical and statistical techniques. Reaction probabilities, integral cross sections, and differential cross sections have been obtained with a statistical quantum method and the mean potential phase space theory. The statistical predictions are compared to exact results calculated by means of time dependent wave packet methods and a previously reported time independent exact quantum mechanical approach using the double many-body expansion (DMBE IV) potential energy surface (PES) [Pastrana et al., J. Phys. Chem. 94, 8073 (1990)] and the recently developed surface (denoted XXZLG) by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. The statistical approaches are found to reproduce only some of the exact total reaction probabilities for low total angular momenta obtained with the DMBE IV PES and some of the cross sections calculated at energy values close to the reaction threshold for the XXZLG surface. Serious discrepancies with the exact integral cross sections at higher energy put into question the possible statistical nature of the title reaction. However, at a collision energy of 1.6 eV, statistical rotationally resolved cross sections managed to reproduce the experimental cross sections for the H+O(2)(v=0,j=1)-->OH(v(')=1,j('))+O process reasonably well.     

A quantum wave packet dynamics study of the N(2D) + H2 reaction

We report a dynamics study of the reaction N((2)D) + H(2) (v=0, j=0-5) --> NH + H using the time-dependent quantum wave packet method and a recently reported single-sheeted double many-body expansion potential energy surface for NH(2)(1(2)A' ') which has been modeled from accurate ab initio multireference configuration-interaction calculations. The calculated probabilities for (v=0, j=0-5) are shown to display resonance structures, a feature also visible to some extent in the calculated total cross sections for (v=0, j=0). A comparison between the calculated centrifugal-sudden and coupled-channel reaction probabilities validate the former approximation for the title system. Rate constants calculated using a uniform J-shifting scheme and averaged over a Boltzmann distribution of rotational states are shown to be in good agreement with the available experimental values. Comparisons with other theoretical results are also made.     

Dynamics of (H(-), H(2)) collisions: a time-dependent quantum mechanical investigation on a new ab initio potential energy surface

A global analytical potential energy surface for the ground state of H(3)(-) has been constructed by fitting an analytic function to the ab initio potential energy values computed using coupled cluster singles and doubles with perturbative triples [CCSD(T)] method and Dunning's augmented correlation consistent polarized valence triple zeta basis set. Using this potential energy surface, time-dependent quantum mechanical wave packet calculations were carried out to calculate the reaction probabilities (P(R)) for the exchange reaction H(-)+H(2)(v, j)-->H(2)+H(-), for different initial vibrational (v) and rotational (j) states of H(2), for total angular momentum equal to zero. With increase in v, the number of oscillations in the P(R)(E) plot increases and the oscillations become more pronounced. While P(R) increases with increase in rotational excitation from j=0 to 1, it decreases with further increase in j to 2 over a wide range of energies. In addition, rotational excitation quenches the oscillations in P(R)(E) plots.     

Vibrational deexcitation and rotational excitation of H(2) and D(2) scattered from Cu(111): Adiabatic versus non-adiabatic dynamics

We have studied survival and rotational excitation probabilities of H(2)(v(i) = 1, J(i) = 1) and D(2)(v(i) = 1, J(i) = 2) upon scattering from Cu(111) using six-dimensional (6D) adiabatic (quantum and quasi-classical) and non-adiabatic (quasi-classical) dynamics. Non-adiabatic dynamics, based on a friction model, has been used to analyze the role of electron-hole pair excitations. Comparison between adiabatic and non-adiabatic calculations reveals a smaller influence of non-adiabatic effects on the energy dependence of the vibrational deexcitation mechanism than previously suggested by low-dimensional dynamics calculations. Specifically, we show that 6D adiabatic dynamics can account for the increase of vibrational deexcitation as a function of the incidence energy, as well as for the isotope effect observed experimentally in the energy dependence for H(2)(D(2))/Cu(100). Furthermore, a detailed analysis, based on classical trajectories, reveals that in trajectories leading to vibrational deexcitation, the minimum classical turning point is close to the top site, reflecting the multidimensionally of this mechanism. On this site, the reaction path curvature favors vibrational inelastic scattering. Finally, we show that the probability for a molecule to get close to the top site is higher for H(2) than for D(2), which explains the isotope effect found experimentally.     

Effect of surface modes on the six-dimensional molecule-surface scattering dynamics of H2-Cu(100) and D2-Cu(111) systems

We include the phonon modes originating from the three layers of Cu(100)/Cu(111) surface atoms on the dynamics of molecular [H(2)(v,j)/D(2)(v,j)] degrees of freedom (DOFs) through a mean field approach, where the surface temperature is incorporated into the effective Hamiltonian (potential) either by considering Boltzmann probability (BP) or by including the Bose-Einstein probability (BEP) factor for the initial state distribution of the surface modes. The formulation of effective potential has been carried out by invoking the expression of transition probabilities for phonon modes known from the "stochastic" treatment of linearly forced harmonic oscillator (LFHO). We perform four-dimensional (4D?2D) as well as six-dimensional (6D) quantum dynamics on a parametrically time and temperature-dependent effective Hamiltonian to calculate elastic/inelastic scattering cross-section of the scattered molecule for the H(2)(v,j)-Cu(100) system, and dissociative chemisorption-physisorption for both H(2)(v,j)-Cu(100) and D(2)(v,j)-Cu(111) systems. Calculated sticking probabilities by either 4D?2D or 6D quantum dynamics on an effective potential constructed by using BP factor for the initial state distribution of the phonon modes could not show any surface temperature dependence. In the BEP case, (a) both 4D?2D and 6D quantum dynamics demonstrate that the phonon modes of the Cu(100) surface affect the state-to-state transition probabilities of the scattered H(2) molecule substantially, and (b) the sticking probabilities due to the collision of H(2) on Cu(100) and D(2) on Cu(111) surfaces show noticeable and substantial change, respectively, as function of surface temperature only when the quantum dynamics of all six molecular DOFs are treated in a fully correlated manner (6D).     

Dissociative chemisorption of H2 on the Cu(110) surface: a quantum and quasiclassical dynamical study

Six-dimensional quantum dynamical and quasiclassical trajectory (QCT) calculations are reported for the reaction and vibrationally inelastic scattering of (v = 0,1,j = 0) H(2) scattering from Cu(110), and for the reaction and rovibrationally elastic and inelastic scattering of (v = 1,j = 1) H(2) scattering from Cu(110). The dynamics results were obtained using a potential energy surface obtained with density functional theory using the PW91 functional. The reaction probabilities computed with quantum dynamics for (v = 0,1,j = 0) were in excellent agreement with the QCT results obtained earlier for these states, thereby validating the QCT approach to sticking of hydrogen on Cu(110). The vibrational de-excitation probability P(v=1,j = 0 --> v = 0) computed with the QCT method is in remarkably good agreement with the quantum dynamical results for normal incidence energies E(n) between 0.2 and 0.6 eV. The QCT result for the vibrational excitation probability P(v = 0,j = 0 --> v = 1) is likewise accurate for E(n) between 0.8 and 1 eV, but the QCT method overestimates vibrational excitation for lower E(n). The QCT method gives probabilities for rovibrationally (in)elastic scattering, P(v = 1,j = 1 --> v('),j(')), which are in remarkably good agreement with quantum dynamical results. The rotationally averaged, initial vibrational state-selective reaction probability obtained with QCT agrees well with the initial vibrational state-selective reaction probability extracted from molecular beam experiments for v = 1, for the range of collision energies for which the v=1 contribution to the measured total sticking probability dominates. The quantum dynamical probabilities for rovibrationally elastic scattering of (v = 1,j = 1) H(2) from Cu(110) are in good agreement with experiment for E(n) between 0.08 and 0.25 eV.     

A full dimensional time-dependent wave packet study for the H4 four-center, collision induced dissociation, and single exchange reactions: reaction probabilities for J=0

A time-dependent initial state selected wave packet method has been developed to study the H2(v(1)=10-11,j1=0)+H2'(v2=0,j2=0)-->HH'+HH' four-center (4C) reaction, and two other competing reactions: the H2+H2'-->H+H+H2' collision induced dissociation (CID) and the H2+H2'-->H+HH'+H' single exchange (SE) reaction, in full six dimensions. Initial state-specific total reaction probabilities for these three competing reactions are presented for total angular momentum J=0 and the effects of reagent vibration on reactions are examined. It is found that (a) the CID process is the dominant process over the whole energy range considered in this study, but the 4C and SE processes also have non-negligible probabilities; (b) the SE process has a lower threshold energy than the 4C process, but the SE probability increases slower than the 4C probability as collision energy increases; (c) the vibrational excitation of H2(v1) is much more efficient than translational motion for promoting these processes, in particular to the CID process.     

Constraints at the transition state of the D + H2 reaction: quantum bottlenecks vs. stereodynamics

This article presents a quasiclassical trajectory method for the calculation of cumulative reaction probabilities by sampling of the helicity quantum number of the reagents (k). The method is applied to the D + H(2) reaction at various total angular momentum (J) values, and the helicity-resolved quasiclassical cumulative reaction probabilities are compared to their quantum mechanical counterparts. The agreement between the two sets of results is fairly good. In particular, k-dependent, J-independent reaction thresholds found with quantum methods are reproduced by the quasiclassical calculations. The shift of these thresholds with increasing k, which has been previously attributed to the quantum bottleneck states taking part in the reaction, is revisited and discussed also in terms of the reaction stereodynamics.     

Quantum wave-packet calculation of reaction probabilities, cross sections, and rate constants for Li + H2+ reaction

The Li + H2+(upsilon,j) --> LiH(upsilon',j') + H+ reactive scattering has been studied by using quantum real wave-packet method. The state-to-state and state-to-all reaction probabilities for the entitled collision have been calculated. The probabilities show a smooth variation for all initial rotational quantum states. The J-shifting approximation has been employed to estimate the integral cross sections and thermal rate constants have been calculated.     

The state-to-state-to-state model for direct chemical reactions: application to D+H2-->HD+H

A simple theoretical model is developed to predict the state-to-state dynamics of direct chemical reactions. Motivated by traditional ideas from transition state theory, expressions are derived for the reactive S matrix that may be computed using the local transition state dynamics. The key approximation involves the use of quantum bottleneck states to represent the near separable dynamics taking place near the transition state. Explicit expressions for the S matrix are obtained using a Franck-Condon treatment for the inelastic coupling between internal states of the collision complex. It is demonstrated that the energetic thresholds for various initial reagent states of the D+H(2) reaction can be understood in terms of our theory. Specifically, the helicity of the reagent states are found to correlate directly to the symmetry of the quantum bottleneck states, which thus possess very different thresholds. Furthermore, the rotational product state distributions for D+H(2) are found to be associated with interfering pathways through the quantum bottleneck states.