In-plane and out-of-plane diffraction of H(2) from metal surfaces

Reply to the 'Comment on "Perturbation theory of scattering for grazing-incidence fast-atom diffraction"' by G. A. Bocan, H. Breiss, S. Szilasi, A. Momeni, E. M. S. Casagrande, E. A. Sánchez, M. S. Gravielle and H. Khemliche, Phys. Chem. Chem. Phys., 2023, , DOI: 10.1039/D3CP02486E

In this Reply, we show that criticisms of perturbation theory for grazing-incidence fast-atom diffraction (GIFAD) are ill-founded. We show explicitly that our formulation (W. Allison, S. Miret-Artés and E. Pollak, Phys. Chem. Chem. Phys., 2022, 24, 15851) provides a similar precision in describing the observed phenomena as ab initio potentials. Since that is the main criterion to distinguish between methods, it seems reasonable to conclude that the perturbation approach using a Morse-type potential reproduces the essential aspects of the dynamics correctly. In addition we expand on the historical context and summarize the physical insights provided by our methods.

A setup for grazing incidence fast atom diffraction

We describe a UHV setup for grazing incidence fast atom diffraction (GIFAD) experiments. The overall geometry is simply a source of keV atoms facing an imaging detector. Therefore, it is very similar to the geometry of reflection high energy electron diffraction experiments used to monitor growth at surfaces. Several custom instrumental developments are described making GIFAD operation efficient and straightforward. The difficulties associated with accurately measuring the small scattering angle and the related calibration are carefully analyzed.

Adiabatic models for the quantum dynamics of surface scattering with lattice effects

In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.

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.

Measuring surface phonons using molecular spin-echo

A new method to measure surface phonons with a molecular beam is presented. The method extends the principles of ³He spin-echo spectroscopy, to the more complex case of a molecular beam exchanging energy with the surface. Measurements are presented for inelastic scattering of D₂ from a Cu(111) surface. Similarly to helium spin-echo, experiments can be performed along optimal tilted projections making it possible to resolve energy peaks with a high energy resolution which is not restricted by the spread of energies of the incident beam. Two analysis methods for these molecular spin echo experiments are presented. A classical approach, analogous to that used for helium spin-echo, explains the most dominant excitation peaks measured, whereas a semi-classical approach allows us to identify smaller peaks which are related to the complexity of the multiple spin-rotation states which exist for molecules.

Prominent out-of-plane diffraction in helium scattering from a methyl-terminated Si(111) surface

Due to their electrochemical and oxidative stability, organic-terminated semiconductor surfaces are well suited to applications in, for example, photoelectrodes and electrochemical cells, which explains the lively interest in their detailed characterization. Helium atom scattering (HAS) is a useful tool to carry out such characterization. Here, we have simulated HAS in He/CH3-Si(111) based on density functional theory (DFT) potential energy surfaces (PESs) and multi-configuration time-dependent Hartree (MCTDH) dynamics. Our analysis of HAS shows that most diffraction taking place in this system corresponds to high-order out-of-plane peaks. This is a general trend that does not depend on the specific features of the simulations, such as the inclusion or not of the van der Waals long-range effects. This is the first and only He-surface system for which such huge out-of-plane diffraction has been described. This striking theoretical finding should encourage new experimental developments to confirm this previously unreported effect.

Dynamics of N2 sticking on W(100): the decisive role of van der Waals interactions

The reactive dynamics of N₂ on W(100) has been investigated by means of quasi-classical trajectory calculations using an interpolated six-dimensional potential energy surface (PES) based on density functional theory energies obtained employing the vdW-DF2 functional.

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.

Reactive and nonreactive scattering of N2 from Ru(0001): a six-dimensional adiabatic study

We have studied the dissociative chemisorption and scattering of N(2) on and from Ru(0001), using a six-dimensional quasiclassical trajectory method. The potential energy surface, which depends on all the molecular degrees of freedom, has been built applying a modified Shepard interpolation method to a data set of results from density functional theory, employing the RPBE generalized gradient approximation. The frozen surface and Born-Oppenheimer [Ann. Phys. (Leipzig) 84, 457 (1927)] approximations were used, neglecting phonons and electron-hole pair excitations. Dissociative chemisorption probabilities are found to be very small even for translational energies much higher than the minimum reaction barrier, in good agreement with experiment. A comparison to previous low dimensional calculations shows the importance of taking into account the multidimensional effects of N(2) rotation and translation parallel to the surface. The new calculations strongly suggest a much smaller role of nonadiabatic effects than previously assumed on the basis of a comparison between low dimensional results and experiments [J. Chem. Phys. 115, 9028 (2001)]. Also in agreement with experiment, our theoretical results show a strong dependence of reaction on the initial vibrational state. Computed angular scattering distributions and parallel translation energy distributions are in good agreement with experiments on scattering, but the theory overestimates vibrational and rotational excitations in scattering.

Quantum scattering of fast atoms and molecules on surfaces

We present evidence for the diffraction of light keV atoms and molecules grazingly scattered on LiF(001) and NaCl(001) surfaces. At such energies, the de Broglie wavelength is 2 orders of magnitude smaller that the mean thermal atomic displacement in the crystal. Thus, no coherent scattering was expected and interaction of keV atoms with surfaces is routinely treated with classical mechanics. We show here that well-defined diffraction patterns can be observed indicating that, for grazing scattering, the pertinent wavelength is that associated with the slow motion perpendicular to the surface. The experimental data are well reproduced by an ab initio calculation.

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.

Reactive and Nonreactive Scattering of H2 from a Metal Surface Is Electronically Adiabatic

The Born-Oppenheimer approximation of uncoupled electronic and nuclear motion is a standard tool of the computational chemist. However, its validity for molecule-metal surface reactions, which are important to heterogeneous catalysis, has been questioned because of the possibility of electron-hole pair excitations. We have performed experiments and calculations on the scattering of molecular hydrogen from a catalytically relevant metal surface, obtaining absolute probabilities for changes in the molecule's velocity parallel to the representative Pt(111) surface. The comparison for in-plane and out-of-plane scattering and results for dissociative chemisorption in the same system show that for hydrogen-metal systems, reaction and diffractive scattering can be accurately described using the Born-Oppenheimer approximation.

The role of molecular rotation in activated dissociative adsorption on metal surfaces

The role of molecular rotation in dissociative adsorption of H2 on the activated NiAl(110) metal surface is systematically investigated by means of classical dynamics calculations performed on ab initio six-dimensional potential energy surfaces. The calculations show that molecules rotate abruptly when they are close to the surface and that this rotation allows the molecules to adopt the orientation that is more convenient for dissociation (i.e., nearly parallel to the surface). Also, in reactive sectors of the NiAl(110) unit cell, there is an "angular threshold" below which molecules cannot dissociate. This angular threshold goes down as the incidence energy increases, which explains the rise of the dissociation probability and the fact that it reaches a value close to 1 at incidence energies of the order of 2 eV. The fact that switching on molecular rotation favors dissociation establishes a competition between dissociation and rotational excitation of reflected molecules above the dissociation threshold. Measurements on rotational excitation might thus bring indirect evidence on the dissociation dynamics. Sample calculations for nonactivated Pd(111) and activated Cu(110) metal surfaces suggest that some of these conclusions may be of general validity.

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.

A classical dynamics method for H2 diffraction from metal surfaces

We present a discretization method that allows one to interpret measurements on diffraction of diatomic molecules from solid surfaces using six-dimensional (6D) classical trajectory calculations. It has been applied to the D2NiAl(110) and H2Pd(111) systems (which are models for activated and nonactivated dissociative chemisorption, respectively) using realistic potential energy surfaces obtained from first principles. Comparisons with experimental results and 6D quantum dynamical calculations show that, in general, the method is able to predict the relative intensity of the most important diffraction peaks. We therefore conclude that classical mechanics can be an efficient guide for experimentalists in the search for the most significant diffraction channels.