Diffraction of CH4 from a Metal Surface

Characterisation of magnetic atomic and molecular beamlines for the extraction of empirical scattering-matrices

A recently developed magnetic molecular interferometry technique allows the experimental determination of how the amplitudes and phases of the molecular wave-function change during the collision of a gas phase molecule with a surface. This information, quantified by a scattering-matrix, provides a very stringent benchmark for developing accurate theoretical models as they can also be determined from scattering calculations and are particularly sensitive to the underlying interaction potential. However, the value of this comparison is necessarily limited by the accuracy with which an empirical scattering-matrix can be extracted from the experimental data. This paper presents the methods used to analyse the measurements and uses simulations to determine how various uncertainties in modelling the different magnetic elements which make up the beamline of the apparatus affect the accuracy with which the scattering-matrix can be extracted. It is shown that when signals have a noise level which corresponds to on the order of 1% of the oscillation amplitude, the uncertainties in the modelling do not significantly affect the ability to extract the scattering-matrix elements, with the error in the extracted values increasing to a few percent as the noise in the signals is increased to 10% of the oscillation amplitude. This therefore gives an estimate of the accuracy of the parameters that can be obtained from future measurements.

An examination of phonon-inelastic molecule-metal scattering using reduced density matrix and stochastic wave packet methods

We explore the application of reduced density matrix-based approaches to molecules interacting with the lattice vibrations of metals, an interaction responsible for the temperature dependence of many of the fundamental steps of catalysis. We avoid the use of simple models for the bath and instead use density functional theory to compute all molecule-phonon interactions and the properties of the lattice phonons, for methane scattering from Ir(111). We find that while the large metal mass leads to long bath correlation times, these are not significantly longer than the time over which the reduced density matrix changes due to interactions with the bath. We show that the neglect of memory is reasonable and the use of the Redfield equation is justified. We also show how the commonly used rotating wave approximation is far too severe for this scattering problem. A less restrictive approximation that is nearly exact for our system gives an equation of motion in the Lindblad form. As a result, the Monte Carlo wave packet methods can be used to describe gas-phonon scattering, guaranteeing positivity, and with all couplings derived from first-principles.

Temperature dependence in fast-atom diffraction at surfaces

Grazing incidence fast atom diffraction at crystal surfaces (GIFAD or FAD) has demonstrated coherent diffraction both at effective energies close to one eV (λ_⊥ ≈ 14 pm for He) and at elevated surface temperatures offering high topological resolution and real time monitoring of growth processes. This is explained by a favorable Debye-Waller factor specific to the multiple collision regime of grazing incidence. This paper presents the first extensive evaluation of the temperature behavior between 177 and 1017 K on a LiF surface. Similarly to diffraction at thermal energies (TEAS), an exponential attenuation of the elastic intensity is observed but, contrarily to TEAS, the maximum coherence is not directly reduced by the attraction forces that increase the effective impact energy. It is more influenced by the surface stiffness and appears very sensitive to surface defects.

The trapping of methane on Ir(111): A first-principles quantum study

We implement a fully quantum mechanical study of methane trapping on Ir(111), where the phonons, the molecule-surface interaction, and the molecule-phonon coupling are all computed from first-principles. We find that both the surface corrugation and the phonon coupling vary strongly with molecular orientation and that there is a "chemical" aspect to this due to the catalytic nature of the metal. For example, molecules with reactive orientations can approach close to surface sites with low barriers to dissociation. Moreover, lattice motion can lower the barrier to dissociation, leading to unusual behavior for the phonon coupling. We find good agreement with experiment and two recent classical studies if we average our potential energy surface over several orientations of the molecule. We also find reasonable agreement with a recent study of methane diffraction, although we show that diffraction does not play a major role in trapping on the smooth Ir(111) surface and that trapping obeys normal energy scaling, consistent with experiment. We show that the trapping probability can be sensitive to the temperature at both high and low incidence energies. Relaxation and desorption of trapped particles are examined.

Computational approaches to dissociative chemisorption on metals: towards chemical accuracy

We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (E_b, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate E_b can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol^-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for E_b for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH₄. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of E_b for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate E_b for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O₂, H₂O, NH₃, CO₂, and CH₃OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.

Time-of-flight measurements of the low-energy scattering of CH4 from Ir(111)

We have measured high-resolution time-of-flight (TOF) spectra of methane scattered from an Ir(111) surface at an incident energy of 81 meV. The angular distributions of scattered CH₄ reveal the presence of a sharp and intense specular peak in addition to sharp features corresponding to rotationally inelastic diffraction (RID) peaks along the two main symmetry directions of Ir(111). TOF spectra have been recorded at several RID positions for the two high-symmetry directions. The data show that the scattering dynamics of CH₄ is more complex than the one reported for H₂/D₂, where energy losses in TOF correspond to the expected excitation/deexcitation RID energy transitions. For CH₄, this is the case only for RID peaks showing up far from the specular peak, whereas those appearing close to the specular position present different behaviors, depending on the incident direction. The results are compared with Ne scattering TOF data, which allows to assess the relevance of multiphonon scattering in the energy-exchange process. Finally, we report experimental evidence of selective adsorption resonances detected with CH₄ beams. This will allow characterizing the CH₄-metal surface physisorption well by measuring angular distributions with CH₄ beams.

Material properties particularly suited to be measured with helium scattering: selected examples from 2D materials, van der Waals heterostructures, glassy materials, catalytic substrates, topological insulators and superconducting radio frequency materials

Helium Atom Scattering (HAS) and Helium Spin-Echo scattering (HeSE), together helium scattering, are well established, but non-commercial surface science techniques. They are characterised by the beam inertness and very low beam energy (<0.1 eV) which allows essentially all materials and adsorbates, including fragile and/or insulating materials and light adsorbates such as hydrogen to be investigated on the atomic scale. At present there only exist an estimated less than 15 helium and helium spin-echo scattering instruments in total, spread across the world. This means that up till now the techniques have not been readily available for a broad scientific community. Efforts are ongoing to change this by establishing a central helium scattering facility, possibly in connection with a neutron or synchrotron facility. In this context it is important to clarify what information can be obtained from helium scattering that cannot be obtained with other surface science techniques. Here we present a non-exclusive overview of a range of material properties particularly suited to be measured with helium scattering: (i) high precision, direct measurements of bending rigidity and substrate coupling strength of a range of 2D materials and van der Waals heterostructures as a function of temperature, (ii) direct measurements of the electron-phonon coupling constant λ exclusively in the low energy range (<0.1 eV, tuneable) for 2D materials and van der Waals heterostructures (iii) direct measurements of the surface boson peak in glassy materials, (iv) aspects of polymer chain surface dynamics under nano-confinement (v) certain aspects of nanoscale surface topography, (vi) central properties of surface dynamics and surface diffusion of adsorbates (HeSE) and (vii) two specific science case examples - topological insulators and superconducting radio frequency materials, illustrating how combined HAS and HeSE are necessary to understand the properties of quantum materials. The paper finishes with (viii) examples of molecular surface scattering experiments and other atom surface scattering experiments which can be performed using HAS and HeSE instruments.

Towards bridging the structure gap in heterogeneous catalysis: the impact of defects in dissociative chemisorption of methane on Ir surfaces

The negatively activated region in CH₄ dissociation is attributed to a precursor-mediated mechanism involving surface defects.

Setting benchmarks for modelling gas-surface interactions using coherent control of rotational orientation states

The coherent evolution of a molecular quantum state during a molecule-surface collision is a detailed descriptor of the interaction potential which was so far inaccessible to measurements. Here we use a magnetically controlled molecular beam technique to study the collision of rotationally oriented ground state hydrogen molecules with a lithium fluoride surface. The coherent control nature of the technique allows us to measure the changes in the complex amplitudes of the rotational projection quantum states, and express them using a scattering matrix formalism. The quantum state-to-state transition probabilities we extract reveal a strong dependency of the molecule-surface interaction on the rotational orientation of the molecules, and a remarkably high probability of the collision flipping the rotational orientation. The scattering matrix we obtain from the experimental data delivers an ultra-sensitive benchmark for theory to reproduce, guiding the development of accurate theoretical models for the interaction of H₂ with a solid surface.

A fundamental and predictive understanding of molecule-surface interactions is challenging to obtain. Here the authors report an experimental technique allowing direct measurement of the scattering matrix, which reports on the coherent evolution of quantum states of a molecule scattering from a surface.