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.     

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.     

H+ versus D+) transfer from HOD+ to N2: mode- and bond-selective effects

Reactions of HOD(+) with N(2) have been studied for HOD(+) in its ground state and with one quantum of excitation in each of its vibrational modes: (001)--predominately OH stretch, 0.396 eV, (010)--bend, 0.153 eV, and (100)--predominately OD stretch, 0.293 eV. Integral cross sections and product recoil velocities were recorded for collision energies from threshold to 4 eV. The cross sections for both H(+) and D(+) transfer rise slowly from threshold with increasing collision energy; however, all three vibrational modes enhance reaction much more strongly than equivalent amounts of collision energy and the enhancements remain large even at high collision energy, where the vibration contributes less than 10% of the total energy. Excitation of the OH stretch enhances H(+) transfer by a factor of ～5, but the effect on D(+) transfer is only slightly larger than that from an equivalent increase in collision energy, and smaller than the effect from the much lower energy bend excitation. Similarly, OD stretch excitation strongly enhances D(+) transfer, but has essentially no effect beyond that of the additional energy on H(+) transfer. The effects of the two stretch vibrations are consistent with the expectation that stretching the bond that is broken in the reaction puts momentum in the correct coordinate to drive the system into the exit channel. From this perspective it is quite surprising that bend excitation also results in large (factor of 2) enhancements of both H(+) and D(+) transfer channels, such that its effect on the total cross section at collision energies below ～2 eV is comparable to those from the two stretch modes, even though the bend excitation energy is much smaller. For collision energies above ～2 eV, the vibrational effects become approximately proportional to the vibrational energy, though still much larger than the effects of equivalent addition of collision energy. Measurements of the product recoil velocity distributions show that reaction is direct at all collision energies, with roughly half the products in a sharp peak corresponding to stripping dynamics and half with a broad and approximately isotropic recoil velocity distribution. Despite the large effects of vibrational excitation on reactivity, the effects on recoil dynamics are small, indicating that vibrational excitation does not cause qualitative changes in the reaction mechanism or in the distribution of reactive impact parameters.     

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.     

Quasi-classical trajectory-gaussian binning study of the OH + D2 → HOD(v1',v2',v3') + D angle-velocity and vibrational distributions at a collision energy of 0.28 eV

The angle-velocity and product vibrational state distributions for the OH + D(2) reaction at a collision energy of 0.28 eV have been calculated using the quasi-classical trajectory-gaussian binning (QCT-GB) method and the Wu-Schatz-Lendvay-Fang-Harding (WSLFH) analytical potential energy surface. Comparison with high resolution molecular beam experiments shows that, differing from what happens when using the standard QCT method (i.e., histogram binning), very good results are obtained for both distributions. Hence, the strong differences previously observed between QCT and experimental results mainly come from an inadequate pseudoquantization of HOD rather than from other quantum effects. This is probably the first time that such a high level of agreement between theory and high resolution experimental data has been found in polyatomic reaction dynamics.     

H+ versus D+ transfer from HOD+ to CO2: bond-selective chemistry and the anomalous effect of bending excitation

Reactions of HOD(+) with CO(2) have been studied for HOD(+) in its ground state, and with one quantum of excitation in each of its vibrational modes: (001)--predominantly OH stretch, 0.396 eV; (010)--bend, 0.153 eV; and (100)--predominantly OD stretch, 0.293 eV. Integral cross sections and product recoil velocities were recorded for collision energies from threshold to 3 eV. The cross sections for both H(+) and D(+) transfer rise with increasing collision energy from threshold to ～1 eV, then become weakly dependent of the collision energy. All three vibrational modes enhance the total reactivity, but quite mode specifically. The H(+) transfer reaction is enhanced by OH stretch excitation, whereas OD stretch excitation has little effect. Conversely, the D(+) transfer reaction is enhanced by OD stretch excitation, while the OH stretch has little effect. Excitation of the bend strongly enhances both channels. The effects of the stretch excitations are consistent with previous studies of neutral HOD mode-selective chemistry, and can be at least qualitatively understood in terms of a late barrier to product formation. The fact that bend excitation produces the largest overall enhancement is surprising, because this is the lowest energy excitation, and is not obviously connected with the reaction coordinates for either H(+) or D(+) transfer. A rationalization in terms of the effects of water distortion on the potential surface is proposed.     

Dynamics of the D+ + H2 and H+ + D2 reactions: a detailed comparison between theory and experiment

An extensive set of experimental measurements on the dynamics of the H(+) + D(2) and D(+) + H(2) ion-molecule reactions is compared with the results of quantum mechanical (QM), quasiclassical trajectory (QCT), and statistical quasiclassical trajectory (SQCT) calculations. The dynamical observables considered include specific rate coefficients as a function of the translational energy, E(T), thermal rate coefficients in the 100-500 K temperature range. In addition, kinetic energy spectra (KES) of the D(+) ions reactively scattered in H(+) + D(2) collisions are also presented for translational energies between 0.4 eV and 2.0 eV. For the two reactions, the best global agreement between experiment and theory over the whole energy range corresponds to the QCT calculations using a gaussian binning (GB) procedure, which gives more weight to trajectories whose product vibrational action is closer to the actual integer QM values. The QM calculations also perform well, although somewhat worse over the more limited range of translational energies where they are available (E(T) < 0.6 eV and E(T) < 0.2 eV for the H(+) + D(2) and D(+) + H(2) reactions, respectively). The worst agreement is obtained with the SQCT method, which is only adequate for low translational energies. The comparison between theory and experiment also suggests that the most reliable rate coefficient measurements are those obtained with the merged beams technique. It is worth noting that none of the theoretical approaches can account satisfactorily for the experimental specific rate coefficients of H(+) + D(2) for E(T)≤ 0.2 eV although there is a considerable scatter in the existing measurements. On the whole, the best agreement with the experimental laboratory KES is obtained with the simulations carried out using the state resolved differential cross sections (DCSs) calculated with the QCT-GB method, which seems to account for most of the observed features. In contrast, the simulations with the SQCT data predict kinetic energy spectra (KES) considerably cooler than those experimentally determined.     

Vibrational inelastic and charge transfer processes in H(+)+H2 system: an ab initio study

State-resolved differential cross sections, total and integral cross sections, average vibrational energy transfer, and the relative probabilities are computed for the H(+)+H2 system using the newly obtained ab initio potential energy surfaces at the full CI/cc-pVQZ level of accuracy which allow for both the direct vibrational inelastic and the charge transfer processes. The quantum dynamics is treated within the vibrational close-coupling infinite-order-sudden approximation approach using the two ab initio quasidiabatic potential energy surfaces. The computed collision attributes for both the processes are compared with the available state-to-state scattering experiments at E(c.m.)=20 eV. The results are in overall good agreement with most of the observed scattering features such as rainbow positions, integral cross sections, and relative vibrational energy transfers. A comparison with the earlier theoretical study carried out on the semiempirical surfaces (diatomics in molecules) is also made to illustrate the reliability of the potential energy surfaces used in the present work.     

New results for the OH (nu = 0,j = 0) + CO (nu = 0,j = 0) --> H + CO2 reaction: Five- and full-dimensional quantum dynamical study on several potential energy surfaces

Full- [six-dimensional (6-D)] and reduced-dimensional [five-dimensional (5-D)] quantum wave packet calculations have been performed for the title reaction to obtain reaction probabilities deriving from the ground rovibrational states of OH and CO with total angular momentum J = 0. Three potential energy surfaces (PES) are studied, namely, those of Bradley and Schatz (BS), Yu, Muckerman, and Sears (YMS), and Lakin, Troya, Schatz, and Harding (LTSH). 6-D calculations are performed only for the BS PES, while 5-D results are reported for all three PES'. The 6-D results obtained in the present work improve on those previously reported, since a larger vibrational basis and a better representation of the OH and CO bonds has been introduced. In particular, we now employ a generalized Lanczos-Morse discrete variable representation for both the OH and CO vibrations. In a further improvement, the generalized discrete variable representation of the CO vibration is based on different CO intramolecular potentials for the asymptotic and product grids employed in our projection formalism. This new treatment of the vibrational bases allows for a large reduction in computation time with respect to our previous implementation of the wave packet method, for a given level of accuracy. As a result, we have been able to extend the range of collision energies for which we can obtain converged 6-D results to a higher energy (0.8 eV) than was possible before (0.5 eV). The comparison of the new 6-D and previous 5-D results for the BS PES shows good agreement of the general trend in the reaction probabilities over all collision energies considered (0.1-0.8 eV), while our previous 6-D calculation showed reaction probabilities that differed from the 5-D results by up to 10% between 0.5 and 0.8 eV. The 5-D reaction probabilities reveal interesting trends for the different PES'. In particular, at low energies (< 0.2 eV) the LTSH PES gives rise to much larger reactivity than the other PES', while at high energies (> 0.3 eV) its reaction probability decreases with respect to the BS and YMS PES', being more than a factor of 2 smaller at 0.8 eV. A 5-D calculation on a modified version of the LTSH surface shows that the van der Waals interaction in the entrance channel, which is not correctly described in the other PES' is largely responsible for its larger reactivity at low energies. The large difference between the 5-D reaction probabilities for the YMS and LTSH PES' serves to emphasize the importance of the van der Waals interaction for the reactivity at low energies, because most of the stationary point energies on the YMS and LTSH PES are rather similar, being in line with high-level ab initio information.     

Influence of nuclear exchange on nonadiabatic electron processes in H(+)+H2 collisions

H(+)+H(2) collisions are studied by means of a semiclassical approach that explicitly accounts for nuclear rearrangement channels in nonadiabatic electron processes. A set of classical trajectories is used to describe the nuclear motion, while the electronic degrees of freedom are treated quantum mechanically in terms of a three-state expansion of the collision wavefunction. We describe electron capture and vibrational excitation, which can also involve nuclear exchange and dissociation, in the E = 2-1000 eV impact energy range. We compare dynamical results obtained with two parametrizations of the potential energy surface of H(3)(+) ground electronic state. Total cross sections for E > 10 eV agree with previous results using a vibronic close-coupling expansion, and with experimental data for E < 10 eV. Additionally, some prototypical features of both nuclear and electron dynamics at low E are discussed.     

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.     

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.     

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.     

Molecular beam scattering studies of the reaction D + H2 (v = 0) and D + H2 (v = 1) → HD + H

The reaction D + H₂ → HD + H has been investigated in two molecular beam scattering experiments. Angular and time-of-flight distributions have been measured for the initial vibrational ground state (v = 0) at a most probable collision energy of E_cm = 1.5 eV and for the first vibrational excited state (v = 1) at E_cm = 0.28 eV with the same apparatus. Results for the ground-state experiment are compared with quasiclassical trajectory calculations(QCT) on the LSTH-hypersurface transformed into the laboratory system and averaged over the apparatus distributions. The agreement isquite satisfactory. At this high collision energy the HD products are no longer scattered in a backward direction but in a wide angular region concentrated about θ = 90° in the center-of-mass system. The absolute reactive cross section has been determined and the agreement with the theoretical value from QCT calculations is within the experimental error. The high sensitivity of the experiment to different properties of the doubly differential cross section has also been demonstrated. A preliminary evaluation of the experiment with initial vibrational excitation (v = 1) shows that the HD-product molecules are preferably backward scattered and the change of internal energy is small supporting the concept of a reaction which is adiabatic with respect to the internal degrees of freedom.     

A pulsed-field ionization photoelectron secondary ion coincidence study of the H2+ (X,upsilon+=0-15,N+=1)+He proton transfer reaction

The endothermic proton transfer reaction, H2+(upsilon+)+He-->HeH+ + H(DeltaE=0.806 eV), is investigated over a broad range of reactant vibrational levels using high-resolution vacuum ultraviolet to prepare reactant ions either through excitation of autoionization resonances, or using the pulsed-field ionization-photoelectron-secondary ion coincidence (PFI-PESICO) approach. In the former case, the translational energy dependence of the integral reaction cross sections are measured for upsilon+=0-3 with high signal-to-noise using the guided-ion beam technique. PFI-PESICO cross sections are reported for upsilon+=1-15 and upsilon+=0-12 at center-of-mass collision energies of 0.6 and 3.1 eV, respectively. All ion reactant states selected by the PFI-PESICO scheme are in the N+=1 rotational level. The experimental cross sections are complemented with quasiclassical trajectory (QCT) calculations performed on the ab initio potential energy surface provided by Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. The QCT cross sections are significantly lower than the experimental results near threshold, consistent with important contributions due to resonances observed in quantum scattering studies. At total energies above 2 eV, the QCT calculations are in excellent agreement with the present results. PFI-PESICO time-of-flight (TOF) measurements are also reported for upsilon+=3 and 4 at a collision energy of 0.6 eV. The velocity inverted TOF spectra are consistent with the prevalence of a spectator-stripping mechanism.     

Quasi-classical trajectory calculations analyzing the reactivity and dynamics of asymmetric stretch mode excitations of methane in the H + CH4 reaction

An exhaustive dynamics study was performed at two collision energies, 1.52 and 2.20 eV, analyzing the effects of the asymmetric (nu3) stretch mode excitation in the reactivity and dynamics of the gas-phase H + CH4 reaction. Quasi-classical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were performed on an analytical potential energy surface previously developed by our group. First, strong coupling between different vibrational modes in the entry channel was observed, indicating that energy can flow between these modes, and therefore that they do not preserve their adiabatic character along the reaction path; i.e., the reaction is nonadiabatic. Second, we found that the reactant vibrational excitation has a significant influence on the vibrational and rotational product distributions. With respect to the vibrational distribution, our results confirm the purely qualitative experimental evidence, although the theoretical results presented here are also quantitative. The rotational distributions are predictive, because no experimental data have been reported. Third, with respect to the reactivity, we found that the nu3 mode excitation by one quantum is more reactive than the ground state by a factor of about 2, independently of the collision energy, and in agreement with the experimental measurement of 3.0 +/- 1.5. Fourth, the state-to-state angular distributions of the products reproduce the experimental behavior at 1.52 eV, where the CH3 products scatter sideways and backward. At 2.20 eV this experimental information is not available, and therefore the results reported here are again predictive. The satisfactory reproduction of a great variety of experimental data by the present QCT study lends confidence to the potential energy surface constructed by our group and to those results whose accuracy cannot be checked by comparison with experiment.     

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.     

Ab initio potential energy surfaces and nonadiabatic collision dynamics in H(+)+O(2) system

The adiabatic potential energy surfaces for the lowest five electronic states of (3)A" symmetry for the H(+)+O(2) collision system have been obtained at the multireference configuration interaction level of accuracy using Dunning's correlation consistent polarized valence triple zeta basis set. The radial nonadiabatic coupling terms and the mixing angle between the lowest two electronic states (1 (3)A" and 2 (3)A"), which adiabatically correlate in the asymptotic limit to H((2)S)+O(2) (+)(X (2)Pi(g)) and H(+)+O(2)(X (3)Sigma(g)(-)), respectively, have been computed using ab initio procedures at the same level of accuracy to yield the corresponding quasidiabatic potential energy matrix. The computed strengths of the vibrational coupling matrix elements reflect the trend observed for inelastic vibrational excitations of O(2) in the experiments at collision energy of 9.5 eV. The quantum dynamics has been preformed on the newly obtained coupled quasidiabatic potential energy surfaces under the vibrational close-coupling rotational infinite-order sudden framework at the experimental collision energy of 9.5 eV. The present theoretical results for vibrational elastic/inelastic excitations of O(2) are in overall good agreement with the available experimental data obtained from the proton energy-loss spectra in molecular beam experiments [F. A. Gianturco et al., J. Phys. B 14, 667 (1981)]. The results for the complementary charge transfer processes are also presented at this collision energy.     

Dynamically biased statistical model for the ortho/para conversion in the H(2)+H(3) (+) → H(3) (+)+ H(2) reaction

In this work we present a dynamically biased statistical model to describe the evolution of the title reaction from statistical to a more direct mechanism, using quasi-classical trajectories (QCT). The method is based on the one previously proposed by Park and Light [J. Chem. Phys. 126, 044305 (2007)]. A recent global potential energy surface is used here to calculate the capture probabilities, instead of the long-range ion-induced dipole interactions. The dynamical constraints are introduced by considering a scrambling matrix which depends on energy and determine the probability of the identity/hop/exchange mechanisms. These probabilities are calculated using QCT. It is found that the high zero-point energy of the fragments is transferred to the rest of the degrees of freedom, what shortens the lifetime of H(5) (+) complexes and, as a consequence, the exchange mechanism is produced with lower proportion. The zero-point energy (ZPE) is not properly described in quasi-classical trajectory calculations and an approximation is done in which the initial ZPE of the reactants is reduced in QCT calculations to obtain a new ZPE-biased scrambling matrix. This reduction of the ZPE is explained by the need of correcting the pure classical level number of the H(5) (+) complex, as done in classical simulations of unimolecular processes and to get equivalent quantum and classical rate constants using Rice-Ramsperger-Kassel-Marcus theory. This matrix allows to obtain a ratio of hop/exchange mechanisms, α(T), in rather good agreement with recent experimental results by Crabtree et al. [J. Chem. Phys. 134, 194311 (2011)] at room temperature. At lower temperatures, however, the present simulations predict too high ratios because the biased scrambling matrix is not statistical enough. This demonstrates the importance of applying quantum methods to simulate this reaction at the low temperatures of astrophysical interest.