Quantum and classical dynamics of reactive scattering of H₂ from metal surfaces

Abstract

We review the state-of-the art in dynamics calculations on the reactive scattering of H2 from metal surfaces, which is an important model system of an elementary reaction that is relevant to heterogeneous catalysis. In many applications, quantum dynamics and classical trajectory calculations are performed within the Born-Oppenheimer static surface model. However, ab initio molecular dynamics (AIMD) is finding increased use in applications aimed at modeling the effect of surface phonons on the dynamics. Molecular dynamics with electronic friction has been used to model the effect of electron-hole pair excitation. Most applications are still based on potential energy surfaces (PESs) or forces computed with density functional theory (DFT), using a density functional within the generalized gradient approximation to the exchange-correlation energy. A new development is the use of a semi-empirical version of DFT (the specific reaction parameter (SRP) approach to DFT). We also discuss the accurate methods that have become available to represent electronic structure data for the molecule-surface interaction in global PESs. It has now become possible to describe highly activated H2 + metal surface reactions with chemical accuracy using the SRP-DFT approach, as has been shown for H2 + Cu(111) and Cu(100). However, chemical accuracy with SRP-DFT has yet to be demonstrated for weakly activated systems like H2 + Ru(0001) and non-activated systems like H2 + Pd(111), for which SRP DFs are not yet available. There is now considerable evidence that electron-hole pair (ehp) excitation does not need to be modeled to achieve the (chemically) accurate calculation of dissociative chemisorption and scattering probabilities. Dynamics calculations show that phonons can be safely neglected in the chemically accurate calculation of sticking probabilities on cold metal surfaces for activated systems, and in the calculation of a number of other observables. However, there is now sufficient evidence to suggest that the decision on whether or not to neglect phonons should be taken with care, with appropriate consideration of the observable to be computed and of the relevant surface temperature. AIMD calculations have provided valuable insights into the mechanisms that are operative in the dissociative adsorption and absorption of hydrogen on/in precovered metal surfaces. Classical and quantum dynamics calculations have shown that the reaction probability of H2 on Pt surfaces consisting of (100) steps and (111) terraces can to a very good approximation be computed as a weighted average of the reactivities on the steps and terraces. Progress obtained with dynamics calculations on the scattering of H2 from alloys and from simple low index metal surfaces is also reported. Insights that may be obtained on the reactivity of a metal surface from the prominent presence of out-of-plane diffraction or, conversely, the complete absence of diffraction, are discussed. A new field has been opened up by experiments on H2 scattering from surfaces at fast graz...

Quantum State–Resolved Studies of Chemisorption Reactions

“…In a recent study using double resonance excitation by RAP, Hundt et al (52) measured different efficacies for the three different symmetry components of the 2ν 3 overtone for CH 4 dissociation on Pt (111). The authors attributed the observed mode specificity to different degrees of C-H stretch vibrational amplitude localization in the dissociating bond (75).…”

“…All these results demonstrate that the intramolecular vibrational energy redistribution (IVR) in the reactant molecule is slow on the several-hundred-femtosecond timescale of the reactive collision (76), in contradiction with the basic assumption of fast and complete IVR made by statistical rate theories (77-80). Killelea et al (81) were the first to demonstrate vibrational bond selectivity in the dissociative chemisorption of methane on Ni (111). They showed that excitation of the C-H stretch fundamental ν 1 in CHD 3 results in a C-H:C-D cleavage branching ratio greater than 30:1.…”

“…The subsurface D atom titration method is limited to studying bond selectivity in CHD 3 on Ni(111). Using RAIRS detection, Chen et al (82) extended the bond selectivity studies to all the partially deuterated isotopologues of methane on Pt (111). The RAIR spectra in Figure 2 show that excitation of a C-H stretch fundamental leads exclusively to C-H bond cleavage for CH 3 D(ν 4 ), CH 2 D 2 (ν 6 ), and CHD 3 (ν 1 ).…”

“…The influence of the surface temperature on the sticking coefficient has been studied by Campbell and coworkers (92,93) for CH 4 (ν 3 ) on Ni (111). Oxygen and hydrogen titrations were used to quantify the adsorbed species.…”

H2/LiF(001) diffractive scattering under fast grazing incidence using a DFT-based potential energy surface

“…Six-dimensional quantum dynamics has already become a standard tool to study diatomic molecule/surface interactions at low incidence energy(<2 eV), from reactive scattering [20][21][22][23][24] to molecular diffraction [24][25][26]. However, its accuracy in reproducing experimental data, disregarding phonons * cristina.diaz@uam.es and electron-hole pair excitations effects, relies on the accuracy of the underlying potential energy surface (PES) on which the dynamics is carried out [23,27].…”