A time-dependent wave packet quantum scattering study of the reaction HD+ (v = 0 - 3;j0 = 1) + He --> HeH+(HeD+) + D(H)

Time-dependent wave packet quantum scattering (TWQS) calculations are presented for HD(+) (v = 0 - 3;j(0)=1) + He collisions in the center-of-mass collision energy (E(T)) range of 0.0-2.0 eV. The present TWQS approach accounts for Coriolis coupling and uses the ab initio potential energy surface of Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. For a fixed total angular momentum J, the energy dependence of reaction probabilities exhibits quantum resonance structure. The resonances are more pronounced for low J values and for the HeH(+) + D channel than for the HeD(+) + H channel and are particularly prominent near threshold. The quantum effects are no longer discernable in the integral cross sections, which compare closely to quasiclassical trajectory calculations conducted on the same potential energy surface. The integral cross sections also compare well to recent state-selected experimental values over the same reactant and translational energy range. Classical impulsive dynamics and steric arguments can account for the significant isotope effect in favor of the deuteron transfer channel observed for HD(+)(v<3) and low translational energies. At higher reactant energies, angular momentum constraints favor the proton-transfer channel, and isotopic differences in the integral cross sections are no longer significant. The integral cross sections as well as the J dependence of partial cross sections exhibit a significant alignment effect in favor of collisions with the HD(+) rotational angular momentum vector perpendicular to the Jacobi R coordinate. This effect is most pronounced for the proton-transfer channel at low vibrational and translational energies.     

Exact quantum dynamics study of the O(+)+H2(v=0,j=0)-->OH(+)+H ion-molecule reaction and comparison with quasiclassical trajectory calculations

The close-coupling hyperspherical (CCH) exact quantum method was used to study the title barrierless reaction up to a collision energy (E(T)) of 0.75 eV, and the results compared with quasiclassical trajectory (QCT) calculations to determine the importance of quantum effects. The CCH integral cross section decreased with E(T) and, although the QCT results were in general quite similar to the CCH ones, they presented a significant deviation from the CCH data within the 0.2-0.6 eV collision energy range, where the QCT method did not correctly describe the reaction probability. A very good accord between both methods was obtained for the OH(+) vibrational distribution, where no inversion of population was found. For the OH(+) rotational distributions, the agreement between the CCH and QCT results was not as good as in the vibrational case, but it was satisfactory in many conditions. The kk(') angular distribution showed a preferential forward character, and the CCH method produced higher forward peaks than the QCT one. All the results were interpreted considering the potential energy surface and plots of a representative sampling of reactive trajectories.     

A time-dependent wave-packet quantum scattering study of the reaction H2+(v = 0-2,4,6;j = 1) + He--> HeH+ + H

The quantum scattering dynamics calculation was carried out for the titled reaction in the collision energy range of 0.0-2.4 eV with reactant H(2) (+) in the rotational state j = 1 and vibrational states v = 0-2, 4, and 6. The present time-dependent wave-packet calculation takes into account the Coriolis coupling (CC) and uses the accurate ab initio potential-energy surface of Palmieri et al. [Mol. Phys. 98, 1835 (2000)]. The importance of including the CC quantum scattering calculation has been revealed by the comparison between the CC calculation and the previous coupled state (CS) calculation. The CC total cross sections for the v = 2, 4, and 6 states show collision energy-dependent behaviors different from those based on the CS calculation. Furthermore, the collision energy dependence of the total cross sections obtained in the present CC calculation only exhibits minor oscillations, indicating that the chance is slim for reactive resonances in total cross sections to survive through the partial-wave averaging. The magnitude and profile of the CC total cross sections for v = 0-2 in the collision energy range of 0.0-2.5 eV are found to be consistent with experimental cross sections obtained recently by Tang et al. [J. Chem. Phys. 122, 164301 (2005)] after taking into account the experimental uncertainties.     

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.     

Collision-induced dissociation in (He, H2(+)(v = 0-2; j = 0-3)) system: a time-dependent quantum mechanical investigation

The collision-induced process He + H(2)(+)(v = 0-2; j = 0-3) → He + H + H(+) has been investigated using a time-dependent quantum mechanical wave packet approach, within the centrifugal sudden approximation. The exchange reaction He + H(2)(+) → HeH(+) + H, which has a lower threshold, dominates over the dissociation process over the entire energy range considered in this study. The reaction cross section for both the exchange and dissociation channels and the branching ratio between the two channels have been computed on the McLaughlin-Thompson-Joseph-Sathyamurthy potential-energy surface and compared with the available experimental and quasiclassical trajectory results.     

Quantum scattering LCAC-SW theory studies on reaction probabilities of three-dimensional H + H2 (v, j) → H2 (v′, j′) + H reaction

Upon the Liu, Siegbahn, Truhlar, Horowitz (LSTH) potential energy surface, the reaction probabilities of the three-dimensional (3-D) state-to-state H + H₂ (v, j) →H ₂(v′, j′) + H reaction are calculated with the linear combination of arrangement channels-scattering wavefunction (LCAC-SW) method. In the calculation, the vibration function of H₂ and the radial propagating wave functions are expanded by the real Gauss functions. The calculated threshold energy and the resonating structure are consistent with the results of the accurate quantum scattering calculations, which shows the accuration, simplicity and practicability of the LCAC-SW method. Project supported by the National Natural Science Fondation of China and the Doctoral Foundation of the State Education Commission of China.     

State-to-state reaction probabilities for the H+O2(v,j)-->O+OH(v',j') reaction on three potential energy surfaces

We report state-to-state and total reaction probabilities for J=0 and total reaction probabilities for J=2 and 4 for the title reaction, both for ground-state and initially rovibrationally excited reactants. The results for three different potential energy surfaces are compared and contrasted. The potential energy surfaces employed are the DMBE IV surface by Pastrana et al. [J. Phys. Chem. 94, 8073 (1990)], the surface by Troe and Ushakov (TU) [J. Chem. Phys. 115, 3621 (2001)], and the new XXZLG ab initio surface by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. Our results show that the total reaction probabilities from both the TU and XXZLG surfaces are much smaller in magnitude for collision energies above 1.2 eV compared to the DMBE IV surface. The three surfaces also show different behavior with regards to the effect of initial state excitation. The reactivity is increased on the XXZLG and the TU surfaces and decreased on the DMBE IV surface. Vibrational and rotational product state distributions for the XXZLG and the DMBE IV surface show different behaviors for both types of distributions. Our results show that for energies above 1.25 eV the dynamics on the DMBE IV surface are not statistical. However, there is also evidence that the dynamics on the XXZLG surface are not purely statistical for energies above the onset of the first excited product vibrational state v'=1. The magnitude of the total reaction probability is decreased for J>0 for the DMBE IV and the XXZLG surfaces for ground-state reactants. However, for initially rovibrationally excited reactants, the total reaction probability does not decrease as expected for both surfaces. As a result the total cross section averaged over all Boltzmann accessible rotational states may well be larger than the cross section reported in the literature for j=1.     

Renner-Teller quantum dynamics of the N(2D) + H2-->NH + H reaction

We present the Born-Oppenheimer (BO) and Renner-Teller (RT) quantum dynamics of the reaction (14)N((2)D)+(1)H(2)(X (1)Sigma(g) (+))-->NH(X (3)Sigma(-))+H((2)S), considering the NH(2) electronic states X (2)B(1) and A (2)A(1). These states correlate to the same (2)Pi(u) linear species, are coupled by RT nonadiabatic effects, and give NH(X (3)Sigma(-))+H and NH(a (1)Delta)+H, respectively. We develop the Hamiltonian matrix elements in the R embedding of the Jacobi coordinates and in the adiabatic electronic representation, using the permutation-inversion symmetry, and taking into account the nuclear-spin statistics. Collision observables are calculated via the real wave-packet (WP) and flux methods, using the potential-energy surfaces of Santoro et al. [J. Phys. Chem. A 106, 8276 (2002)]. WP snapshots show that the reaction proceeds via an insertion mechanism, and that the RT-WP avoids the A (2)A(1) potential barrier, jumping from the excited to the ground surface and giving mainly the NH(X (3)Sigma(-)) products. X (2)B(1) BO probabilities and cross sections show large tunnel effects and are approximately four to ten times larger than the A (2)A(1) ones. This implies a BO rate-constant ratio k(X (2)B(1))k(A (2)A(1)) approximately 10(5) at 300 K, i.e., a negligible BO formation of NH(a (1)Delta). When H(2) is rotationally excited, RT couplings reduce slightly the X (2)B(1) reaction observables, but enhance strongly the A (2)A(1) reactivity. These couplings are important at all collision energies, reduce the collision threshold, and increase remarkably reaction probabilities and cross sections. The RT k(A (2)A(1)) is thus approximately 3.3 order of magnitude larger than the BO value, and degeneracy-averaged, initial-state-resolved rate constants increase by approximately 13% and by approximately 47% at 300 and 500 K, respectively. Owing to an overestimation of the X (2)B(1) potential barrier, the calculated thermal rate is too low with respect to that observed, but we obtain a good agreement by shifting down the calculated cross section.     

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.     

Time dependent quantum dynamics study of the Ne + H2+(v=0-4)-->NeH+ + H proton transfer reaction

The Ne + H2+-->NeH+ + H proton transfer reaction was studied using the time dependent real wave packet quantum dynamics method at the helicity decoupling level, considering the H2+ molecular ion in the (v=0-4, j=0) vibrorotational states and a wide collision energy interval. The calculated reaction probabilities and reaction cross sections were in a rather good agreement with reanalyzed previous exact quantum dynamics results, where a much smaller collision energy interval was considered. Also, a quite good agreement with experimental data was found. These results suggested the adequacy of the approach used here to describe this and related systems.     

Time dependent quantum dynamics study of the O++H2(v=0,j=0)-->OH++H ion-molecule reaction and isotopic variants (D2,HD)

The time dependent real wave packet method using the helicity decoupling approximation was used to calculate the cross section evolution with collision energy (excitation function) of the O++H2(v=0,j=0)-->OH++H reaction and its isotopic variants with D2 and HD, using the best available ab initio analytical potential energy surface. The comparison of the calculated excitation functions with exact quantum results and experimental data showed that the present quantum dynamics approach is a very useful tool for the study of the selected and related systems, in a quite wide collision energy interval (approximately 0.0-1.1 eV), involving a much lower computational cost than the quantum exact methods and without a significant loss of accuracy in the cross sections.     

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

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H2+NH2-->H+NH3 reaction using a seven dimensional model on 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 and nonreactive NH2 group keeps C2v symmetry and the rotation-vibration coupling in NH2 is neglected. The total reaction probabilities are calculated when the two reactants are initially at their ground states, when the NH2 bending mode is excited, and when H2 is on its first vibrational excited state, with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants and equilibrium constants are calculated for the temperature range of 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the reaction is dominated by ground-state reactivity and the main contribution to the thermal rate constants is thought to come from this state, (b) the excitation energy of H2 was used to enhance reactivity while the excitation of the NH2 bending mode hampers the reaction, (c) the calculated thermal rate constants are very close to the experimental data and transition state theory results at high and middle temperature, while they are ten times higher than that of transition state theory at low temperature (T=200 K), and (d) the equilibrium constants results indicate that the approximations applied may have different roles in the forward and reverse reactions.     

Quantum dynamics study of the K+HF(v=0-2,j=0)-->KF+H reaction and comparison with quasiclassical trajectory results

Extensive quantum real wave packet calculations within the helicity decoupling approximation are used to analyze the influence of the HF vibrational excitation on the K+HF(v=0-2,j=0)-->KF+H reaction. Quantum reaction probabilities P and reaction cross sections sigma are compared with corresponding quasiclassical trajectory (QCT) results. Disregarding threshold regions for v=0 and 1 (v=2 has no threshold), both approaches lead to remarkably similar results, particularly for sigma, validating the use of the QCT method for this system. When moving from v=0 to v=1 there is a large increase in P and sigma, as expected for a late barrier system. For v=2 the reaction becomes exoergic and P approximately 0.95 (with the exception of large total angular momenta where centrifugal barriers play a role). While substantial vibrational enhancement of the reactivity is thus seen, it is still quite less than that inferred from experimental data in the intermediate and high collision energy ranges. The origin of this discrepancy is unclear.     

Full dimensional time-dependent quantum dynamics study of the H + NH3 --> H2 + NH2 reaction

A rigorous full dimensional time-dependent wave packet method has been developed for the reactive scattering between an atom and a tetra-atomic molecule. The method has been applied to the hydrogen abstraction reaction H+NH(3)-->H(2)+NH(2). Initial state-selected total reaction probabilities are investigated for the reactions from the ground vibrational state and from four excited vibrational states of ammonia. The total reaction probabilities from two lowest "tunneling doublets" due to the inversion barrier for the umbrella bending motion of NH(3) and from two pairs of doubly degenerate vibrational states of NH(3) are also inspected. Integral cross sections and rate constants are calculated for the reaction from the ground state with the centrifugal-sudden approximation. The calculated results are compared with those from the previous seven dimensional calculations [M. Yang and J. C. Corchado, J. Chem. Phys. 126, 214312 (2007)]. This work shows that the full dimensional rate constants are a factor of 3 larger than the corresponding seven dimensional calculated values at T=200 K and are overall smaller than those obtained from the variational transition state theory in the whole temperature region. The work also reveals that nonreactive NH bonds of NH(3) cannot be treated as spectators due to the fact that three NH bonds are coupled with each other during the reaction process.     

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.     

Quantum and classical studies of the O(3P) + H2(v = 0-3,j = 0) --> OH + H reaction using benchmark potential surfaces

We present results of time-dependent quantum mechanics (TDQM) and quasiclassical trajectory (QCT) studies of the excitation function for O(3P) + H2(v = 0-3,j = 0) --> OH + H from threshold to 30 kcal/mol collision energy using benchmark potential energy surfaces [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)]. For H2(v = 0) there is excellent agreement between quantum and classical results. The TDQM results show that the reactive threshold drops from 10 kcal/mol for v = 0 to 6 for v = 1, 5 for v = 2 and 4 for v = 3, suggesting a much slower increase in rate constant with vibrational excitation above v = 1 than below. For H2(v > 0), the classical results are larger than the quantum results by a factor approximately 2 near threshold, but the agreement monotonically improves until they are within approximately 10% near 30 kcal/mol collision energy. We believe these differences arise from stronger vibrational adiabaticity in the quantum dynamics, an effect examined before for this system at lower energies. We have also computed QCT OH(v',j') state-resolved cross sections and angular distributions. The QCT state-resolved OH(v') cross sections peak at the same vibrational quantum number as the H2 reagent. The OH rotational distributions are also quite hot and tend to cluster around high rotational quantum numbers. However, the dynamics seem to dictate a cutoff in the energy going into OH rotation indicating an angular momentum constraint. The state-resolved OH distributions were fit to probability functions based on conventional information theory extended to include an energy gap law for product vibrations.     

On the dynamics of the H+ +D2(v=0,j=0)-->HD+D + reaction: a comparison between theory and experiment

The H+ +D2(v=0,j=0)-->HD+D + reaction has been theoretically investigated by means of a time independent exact quantum mechanical approach, a quantum wave packet calculation within an adiabatic centrifugal sudden approximation, a statistical quantum model, and a quasiclassical trajectory calculation. Besides reaction probabilities as a function of collision energy at different values of the total angular momentum, J, special emphasis has been made at two specific collision energies, 0.1 and 0.524 eV. The occurrence of distinctive dynamical behavior at these two energies is analyzed in some detail. An extensive comparison with previous experimental measurements on the Rydberg H atom with D2 molecules has been carried out at the higher collision energy. In particular, the present theoretical results have been employed to perform simulations of the experimental kinetic energy spectra.     

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.     

Time-dependent quantum dynamics of reactive scattering and the calculation of product quantum state distributions — A study of the collinear F+H2(v=0) → HF(v′)+H reaction

Summary A new method for the calculation of partial cross sections in the time-dependent quantum theory of molecular reactive scattering processes is discussed. Preliminary calculations are presented which clearly illustrate the power of the method. They show how all the partial cross sections associated with a single initial quantum state may be computed over a very wide energy range from a single propagation of a prepared wavepacket. The resonance behaviour is obtained without difficulty and the energies of the reactive scattering resonances are exactly reproduced.