A classical dynamics method for H2 diffraction from metal surfaces

Grazing incidence fast atom and molecule diffraction: theoretical challenges

This perspective article reviews the state-of-the-art of grazing incidence fast atom and molecule diffraction (GIFAD and GIFMD) simulations and addresses the main challenges that theorists, aiming to provide useful inputs in this topic, are facing. We first discuss briefly the methods used to build accurate potential energy surfaces describing the interaction between the projectile and the surface. Subsequently, we focus on the dynamics simulation methods for GIFAD, a phenomenon that has received a lot of experimental attention since 2007, when the first measurements were published. Following this experimental effort, theorists have developed and adapted a bunch of methods able to simulate, analyze and extract information from the experimental outputs. We review these methods, from the very simple ones based on classical dynamics to the full quantum ones, paying special attention to more versatile semiclassical approaches, which include quantum ingredients in the dynamics at a computational cost only slightly higher than that required in classical dynamics. Within the semiclassical framework it is possible, for example, to include in the dynamics the surface phonons and the projectile coherence, two factors that may have a relevant influence on the experimental measurements, at a reasonable computational cost. Finally, we address GIFMD, a phenomenon that has received much less attention and for which there is still a lot of room for research. We review the few examples of GIFMD available in the literature, and we discuss new phenomena associated with the molecular internal degrees of freedom, which may have some impact in other closely related fields, such as molecular reactivity on metal surfaces. Finally, we point out opened questions, raised from the comparisons between theoretical and experimental results, which claim for further experimental efforts.

When classical trajectories get to quantum accuracy: II. The scattering of rotationally excited H2 on Pd(111)

The classical trajectory method in a quantum spirit assigns statistical weights to classical paths on the basis of two semiclassical corrections: Gaussian binning and the adiabaticity correction.

When Classical Trajectories Get to Quantum Accuracy: The Scattering of H2 on Pd(111)

When elementary reactive processes occur at such low energies that only a few states of reactants and/or products are available, quantum effects strongly manifest and the standard description of the dynamics within the classical framework fails. We show here, for H₂ scattering on Pd(111), that by pseudoquantizing in the spirit of Bohr the relevant final actions of the system, along with adequately treating the diffraction-mediated trapping of the incoming wave, classical simulations achieve an unprecedented agreement with state-of-the-art quantum dynamics calculations.

Experimental and theoretical study of rotationally inelastic diffraction of H2(D2) from methyl-terminated Si(111)

Fundamental details concerning the interaction between H2 and CH3-Si(111) have been elucidated by the combination of diffractive scattering experiments and electronic structure and scattering calculations. Rotationally inelastic diffraction (RID) of H2 and D2 from this model hydrocarbon-decorated semiconductor interface has been confirmed for the first time via both time-of-flight and diffraction measurements, with modest j = 0 → 2 RID intensities for H2 compared to the strong RID features observed for D2 over a large range of kinematic scattering conditions along two high-symmetry azimuthal directions. The Debye-Waller model was applied to the thermal attenuation of diffraction peaks, allowing for precise determination of the RID probabilities by accounting for incoherent motion of the CH3-Si(111) surface atoms. The probabilities of rotationally inelastic diffraction of H2 and D2 have been quantitatively evaluated as a function of beam energy and scattering angle, and have been compared with complementary electronic structure and scattering calculations to provide insight into the interaction potential between H2 (D2) and hence the surface charge density distribution. Specifically, a six-dimensional potential energy surface (PES), describing the electronic structure of the H2(D2)/CH3-Si(111) system, has been computed based on interpolation of density functional theory energies. Quantum and classical dynamics simulations have allowed for an assessment of the accuracy of the PES, and subsequently for identification of the features of the PES that serve as classical turning points. A close scrutiny of the PES reveals the highly anisotropic character of the interaction potential at these turning points. This combination of experiment and theory provides new and important details about the interaction of H2 with a hybrid organic-semiconductor interface, which can be used to further investigate energy flow in technologically relevant systems.

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.

Dynamics simulation of N(2) scattering onto W(100,110) surfaces: A stringent test for the recently developed flexible periodic London-Eyring-Polanyi-Sato potential energy surface

An efficient method to construct the six dimensional global potential energy surface (PES) for two atoms interacting with a periodic rigid surface, the flexible periodic London-Eyring-Polanyi-Sato model, has been proposed recently. The main advantages of this model, compared to state-of-the-art interpolated ab initio PESs developed in the past, reside in its global nature along with the small number of electronic structure calculations required for its construction. In this work, we investigate to which extent this global representation is able to reproduce the fine details of the scattering dynamics of N(2) onto W(100,110) surfaces reported in previous dynamics simulations based on locally interpolated PESs. The N(2)/W(100) and N(2)/W(110) systems are chosen as benchmarks as they exhibit very unusual and distinct dissociative adsorption dynamics although chemically similar. The reaction pathways as well as the role of dynamic trapping are scrutinized. Besides, elastic/inelastic scattering dynamics including internal state and angular distributions of reflected molecules are also investigated. The results are shown to be in fair agreement with previous theoretical predictions.

Molecular effects in H2 scattering from metal surfaces at grazing incidence

Collisions of fast atoms with surfaces at grazing incidence have been recently proposed as a promising new tool to determine surface parameters with unprecedented accuracy. Here we show, by means of classical dynamics calculations performed with first-principles six-dimensional potential energy surfaces for H2/NiAl(110) and H2/Pd(111) that, under grazing incidence conditions, fast light molecular projectiles are also useful to determine sticking probabilities at thermal energies, from the threshold up to the saturation limit. Thus they are the ideal complement to traditional experiments at thermal energies to determine sticking curves up to the saturation limit.

Experimental evidence of dynamic trapping in the scattering of H2 from Pd(110)

We have performed H2(D2) diffraction experiments on a Pd(110) surface using two different high-sensitivity set-ups. We have found that, although the total reflectivity of Pd(110) is comparable to that observed in other reactive systems, the corresponding H2(D2) diffraction patterns are quite different: no diffraction peak, including the specular one, is observed on Pd(110). This unexpected result is the consequence of dynamic trapping. Such interpretation is supported by classical dynamics calculations based on accurate ab initio potential energy surfaces.

Adsorption and scattering of H2 and D2 by NiAl(110)

We present quasiclassical dynamics calculations of H2 and D2 scattering by the NiAl(110) surface using a recently proposed six-dimensional potential-energy surface (PES) obtained from density-functional theory calculations. The results for dissociative adsorption confirm several experimental predictions using (rotationally hot) D2 beams, namely, the existence of a dissociation barrier, the small isotopic effect, the importance of vibrational enhancement, and the existence of normal energy scaling. The latter conclusion shows that normal energy scaling is not necessarily associated with weak corrugated surfaces. The results for rotationally elastic and inelastic diffractions are also in reasonable agreement with experiment, but they show that many more diffractive transitions are responsible for the observed structures than previously assumed. This points to the validity of the PES recently proposed [P. Riviere, H. F. Busnengo, and F. Martin, J. Chem. Phys. 121, 751 (2004)] to describe dissociative adsorption as well as rotationally elastic and inelastic diffractions in the H2NiAl(110) system.