Multiconfiguration time-dependent Hartree method applied to molecular dissociation on surfaces: H2 + Pt(111).

Four-dimensional quantum dynamics calculations are performed on the dissociative chemisorption of H(2) on Pt(111) using the multiconfiguration time-dependent Hartree method. The aim of this work is to study the performance of the multiconfiguration time-dependent Hartree method for a gas-surface reaction by comparison with the standard time-dependent wave-packet propagation method. The initial-state resolved dissociation probability of H(2) is calculated within two four-dimensional models. The first four-dimensional model treats explicitly the rotational motion of the molecule and the H(2) dissociation is studied above two different sites of the Pt(111). For this model, only a potential-energy surface of general form was available. This potential was refitted to a sum of product form to allow efficient calculations with the multiconfiguration time-dependent Hartree method. The second model focuses on the description of the center-of-mass motion parallel to the surface, the rotational motion of the molecule being frozen. These four-dimensional quantum dynamics calculations yield important insights which can help with performing full six-dimensional calculations on H(2) dissociating on Pt(111). The multiconfiguration time-dependent Hartree method is shown to be particularly efficient for computing initial-state selective dissociation probabilities for the system studied, with a good accuracy and a reduced amount of memory and computational time when compared to the standard time-dependent wave-packet method.

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.