Quantum dynamics of the dissociation of H2 on Cu(100): Dependence of the site-reactivity on initial rovibrational state

We perform six-dimensional (6D) quantum wavepacket calculations for H2 dissociatively adsorbing on Cu(100) from a variety of rovibrational initial states. The calculations are performed on a new potential energy surface (PES), the construction of which is also detailed. Reaction probabilities are in good agreement with experimental findings. Using a new flux analysis method, we calculate the reaction probability density as a function of surface site and collision energy, for a variety of initial states. This approach is used to study the effects of rotation and vibration on reaction at specific surface sites. The results are explained in terms of characteristics of the PES and intrinsically dynamic effects. An important observation is that, even at low collision energies, reaction does not necessarily proceed predominantly in the region of the minimum potential barrier, but can occur almost exclusively at a site with a higher barrier. This suggests that experimental control of initial conditions could be used to selectively induce reaction at particular surface sites. Our predictions for site-reactivity could be tested using contemporary experimental methods: The calculations predict that, for reacting molecules, there will be a dependence of the quadrupole alignment of j on the incident vibrational state, v. This is a direct result of PES topography in the vicinity of the preferred reaction sites of v = 0 and v = 1 molecules. Invoking detailed balance, evidence for this difference in preferred reaction site of v = 0 and 1 molecules could be obtained through associative desorption experiments.