Rotational effects in the dissociative adsorption of H2 on the Pt211 stepped surface

Quantum Dynamics of Dissociative Chemisorption of H2 on the Stepped Cu(211) Surface

Reactions on stepped surfaces are relevant to heterogeneous catalysis, in which a reaction often takes place at the edges of nanoparticles where the edges resemble steps on single-crystal stepped surfaces. Previous results on H₂ + Cu(211) showed that, in this system, steps do not enhance the reactivity and raised the question of whether this effect could be, in any way, related to the neglect of quantum dynamical effects in the theory. To investigate this, we present full quantum dynamical molecular beam simulations of sticking of H₂ on Cu(211), in which all important rovibrational states populated in a molecular beam experiment are taken into account. We find that the reaction of H₂ with Cu(211) is very well described with quasi-classical dynamics when simulating molecular beam sticking experiments, in which averaging takes place over a large number of rovibrational states and over translational energy distributions. Our results show that the stepped Cu(211) surface is distinct from its component Cu(111) terraces and Cu(100) steps and cannot be described as a combination of its component parts with respect to the reaction dynamics when considering the orientational dependence. Specifically, we present evidence that, at translational energies close to the reaction threshold, vibrationally excited molecules show a negative rotational quadrupole alignment parameter on Cu(211), which is not found on Cu(111) and Cu(100). The effect arises because these molecules react with a site-specific reaction mechanism at the step, that is, inelastic rotational enhancement, which is only effective for molecules with a small absolute value of the magnetic rotation quantum number. From a comparison to recent associative desorption experiments as well as Born-Oppenheimer molecular dynamics calculations, it follows that the effects of surface atom motion and electron-hole pair excitation on the reactivity fall within chemical accuracy, that is, modeling these effect shifts extracted reaction probability curves by less than 1 kcal/mol translational energy. We found no evidence in our fully state-resolved calculations for the "slow" reaction channel that was recently reported for associative desorption of H₂ from Cu(111) and Cu(211), but our results for the fast channel are in good agreement with the experiments on H₂ + Cu(211).

Transferability of the Specific Reaction Parameter Density Functional for H2 + Pt(111) to H2 + Pt(211)

The accurate description of heterogeneously catalyzed reactions may require chemically accurate evaluation of barriers for reactions of molecules at the edges of metal nanoparticles. It was recently shown that a semiempirical density functional describing the interaction of a molecule dissociating on a flat metal surface (CHD₃ + Pt(111)) is transferable to the same molecule reacting on a stepped surface of the same metal (Pt(211)). However, validation of the method for additional systems is desirable. To address the question whether the specific reaction parameter (SRP) functional that describes H₂ + Pt(111) with chemical accuracy is also capable of accurately describing H₂ + Pt(211), we have performed molecular beam simulations with the quasi-classical trajectory (QCT) method, using the SRP functional developed for H₂ + Pt(111). Our calculations used the Born-Oppenheimer static surface model. The accuracy of the QCT method was assessed by comparison with quantum dynamics results for reaction of the ro-vibrational ground state of H₂. The theoretical results for sticking of H₂ and D₂ on Pt(211) are in quite good agreement with the experiment, but uncertainties remain because of a lack of accuracy of the QCT simulations at low incidence energies and possible inaccuracies in the reported experimental incidence energies at high energies. We also investigated the nonadiabatic effect of electron-hole pair excitation on the reactivity using the molecular dynamics with the electron friction (MDEF) method, employing the local density friction approximation (LDFA). Only small effects of electron-hole pair excitation on sticking are found.

Site-specific reactivity of molecules with surface defects—the case of H2 dissociation on Pt

The classic system that describes weakly activated dissociation in heterogeneous catalysis has been explained by two dynamical models that are fundamentally at odds. Whereas one model for hydrogen dissociation on platinum(111) invokes a preequilibrium and diffusion toward defects, the other is based on direct and local reaction. We resolve this dispute by quantifying site-specific reactivity using a curved platinum single-crystal surface. Reactivity is step-type dependent and varies linearly with step density. Only the model that relies on localized dissociation is consistent with our results. Our approach provides absolute, site-specific reaction cross sections.

Anomalous Dependence of the Reactivity on the Presence of Steps: Dissociation of D2 on Cu(211)

Stepped metal surfaces are usually assumed to exhibit an increased catalytic activity for bond cleavage of small molecules over their flat single-crystal counterparts. We present experimental and theoretical data on the dissociative adsorption of molecular hydrogen on copper that contradicts this notion. We observe hydrogen molecules to be more reactive on the flat Cu(111) than on the stepped Cu(211) surface. We suggest that this exceptional behavior is due to a geometric effect, that is, that bond cleavage on the flat surface does not occur preferentially over a top site.

Quantum and classical dynamics of reactive scattering of H2 from metal surfaces

State-of-the-art theoretical models allow nowadays an accurate description of H₂/metal surface systems and phenomena relative to heterogeneous catalysis. Here we review the most relevant ones investigated during the last 10 years.

Dynamical pruning of static localized basis sets in time-dependent quantum dynamics

We investigate the viability of dynamical pruning of localized basis sets in time-dependent quantum wave packet methods. Basis functions that have a very small population at any given time are removed from the active set. The basis functions themselves are time independent, but the set of active functions changes in time. Two different types of localized basis functions are tested: discrete variable representation (DVR) functions, which are localized in position space, and phase-space localized (PSL) functions, which are localized in both position and momentum. The number of functions active at each point in time can be as much as an order of magnitude less for dynamical pruning than for static pruning, in reactive scattering calculations of H2 on the Pt(211) stepped surface. Scaling of the dynamically pruned PSL (DP-PSL) bases with dimension is considerably more favorable than for either the primitive (direct product) or DVR bases, and the DP-PSL basis set is predicted to be three orders of magnitude smaller than the primitive basis set in the current state-of-the-art six-dimensional reactive scattering calculations.