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

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.

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.

Stopping molecular rotation using coherent ultra-low-energy magnetic manipulations

Rotational motion lies at the heart of intermolecular, molecule-surface chemistry and cold molecule science, motivating the development of methods to excite and de-excite rotations. Existing schemes involve perturbing the molecules with photons or electrons which supply or remove energy comparable to the rotational level spacing. Here, we study the possibility of de-exciting the molecular rotation of a D₂ molecule, from J = 2 to the non-rotating J = 0 state, without using an energy-matched perturbation. We show that passing the beam through a 1 m long magnetic field, which splits the rotational projection states by only 10⁻¹²eV, can change the probability that a molecule-surface collision will stop a molecule from rotating and lose rotational energy which is 9 orders larger than that of the magnetic manipulation. Calculations confirm that different rotational orientations have different de-excitation probabilities but underestimate rotational flips (∆m_J\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ne$$\end{document}≠0), highlighting the importance of the results as a sensitive benchmark for further developing theoretical models of molecule-surface interactions.

Manipulating the rotational motions of molecules may provide a tool for controlling chemical processes. Here the authors demonstrate that the rotation of a D₂ molecule can be stopped, upon collision with a metal surface, by a magnetic field that affects the rotational levels to a much smaller extent than the energy difference upon de-excitation.

Prominence of Terahertz Acoustic Surface Plasmon Excitation in Gas-Surface Interaction with Metals

The current understanding of the dynamics of gas-surface interactions is that all of the energy lost in the collision is transferred to vibrations of the target. Electronic excitations were shown to play a marginal role except for cases in which the impinging particles have energies of several electronvolts. Here we show that this picture does not hold for metal surfaces supporting acoustic surface plasmons. Such loss, dressed with a vibronic structure, is shown to make up a prominent energy transfer route down to the terahertz region for Ne atoms scattering off Cu(111) and is expected to dominate for most metals. This mechanism determines, e.g., the drag force acting on telecommunication satellites, which are typically gold-plated to reduce overheating by sunshine. The electronic excitations can be unambiguously discerned from the vibrational ones under mild hyperthermal impact conditions.

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.

Chemical reactions on surfaces for applications in catalysis, gas sensing, adsorption-assisted desalination and Li-ion batteries: opportunities and challenges for surface science

The study of chemical processes on solid surfaces is a powerful tool to discover novel physicochemical concepts with direct implications for processes based on chemical reactions at surfaces, largely exploited by industry. Recent upgrades of experimental tools and computational capabilities, as well as the advent of two-dimensional materials, have opened new opportunities and challenges for surface science. In this Perspective, we highlight recent advances in application fields strictly connected to novel concepts emerging in surface science. Specifically, we show for selected case-study examples that surface oxidation can be unexpectedly beneficial for improving the efficiency in electrocatalysis (the hydrogen evolution reaction and oxygen evolution reaction) and photocatalysis, as well as in gas sensing. Moreover, we discuss the adsorption-assisted mechanism in membrane distillation for seawater desalination, as well as the use of surface-science tools in the study of Li-ion batteries. In all these applications, surface-science methodologies (both experimental and theoretical) have unveiled new physicochemical processes, whose efficiency can be further tuned by controlling surface phenomena, thus paving the way for a new era for the investigation of surfaces and interfaces of nanomaterials. In addition, we discuss the role of surface scientists in contemporary condensed matter physics, taking as case-study examples specific controversial debates concerning unexpected phenomena emerging in nanosheets of layered materials, solved by adopting a surface-science approach.

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.

Designing new SRP density functionals including non-local vdW-DF2 correlation for H2 + Cu(111) and their transferability to H2 + Ag(111), Au(111) and Pt(111)

Specific reaction parameter density functionals (SRP-DFs) that can describe dissociative chemisorption molecular beam experiments of hydrogen (H₂) on cold transition metal surfaces with chemical accuracy have so far been shown to be only transferable among different facets of the same metal, but not among different metals. We design new SRP-DFs that include non-local vdW-DF2 correlation for the H₂ + Cu(111) system, and evaluate their transferability to the highly activated H₂ + Ag(111) and H₂ + Au(111) systems and the non-activated H₂ + Pt(111) system. We design our functionals for the H₂ + Cu(111) system since it is the best studied system both theoretically and experimentally. Here we demonstrate that a SRP-DF fitted to reproduce molecular beam sticking experiments for H₂ + Cu(111) with chemical accuracy can also describe such experiments for H₂ + Pt(111) with chemical accuracy, and vice versa. Chemically accurate functionals have been obtained that perform very well with respect to reported van der Waals well geometries, and which improve the description of the metal over current generalized gradient approximation (GGA) based SRP-DFs. From a systematic comparison of our new SRP-DFs that include non-local correlation to previously developed SRP-DFs, for both activated and non-activated systems, we identify non-local correlation as a key ingredient in the construction of transferable SRP-DFs for H₂ interacting with transition metals. Our results are in excellent agreement with experiment when accurately measured observables are available. It is however clear from our analysis that, except for the H₂ + Cu(111) system, there is a need for more, more varied, and more accurately described experiments in order to further improve the design of SRP-DFs. Additionally, we confirm that, when including non-local correlation, the sticking of H₂ on Cu(111) is still well described quasi-classically.

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.

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.

Assessment of Two Problems of Specific Reaction Parameter Density Functional Theory: Sticking and Diffraction of H2 on Pt(111)

It is important that theory is able to accurately describe dissociative chemisorption reactions on metal surfaces, as such reactions are often rate-controlling in heterogeneously catalyzed processes. Chemically accurate theoretical descriptions have recently been obtained on the basis of the specific reaction parameter (SRP) approach to density functional (DF) theory (DFT), allowing reaction barriers to be obtained with chemical accuracy. However, being semiempirical, this approach suffers from two basic problems. The first is that sticking probabilities (to which SRP density functionals (DFs) are usually fitted) might show differences across experiments, of which the origins are not always clear. The second is that it has proven hard to use experiments on diffractive scattering of H₂ from metals for validation purposes, as dynamics calculations using a SRP-DF may yield a rather poor description of the measured data, especially if the potential used contains a van der Waals well. We address the first problem by performing dynamics calculations on three sets of molecular beam experiments on D₂ + Pt(111), using four sets of molecular beam parameters to obtain sticking probabilities, and the SRP-DF recently fitted to one set of experiments on D₂ + Pt(111). It is possible to reproduce all three sets of experiments with chemical accuracy with the aid of two sets of molecular beam parameters. The theoretical simulations with the four different sets of beam parameters allow one to determine for which range of incidence conditions the experiments should agree well and for which conditions they should show specific differences. This allows one to arrive at conclusions about the quality of the experiments and about problems that might affect the experiments. Our calculations on diffraction of H₂ scattering from Pt(111) show both quantitative and qualitative differences with previously measured diffraction probabilities, which were Debye-Waller (DW)-extrapolated to 0 K. We suggest that DW extrapolation, which is appropriate for direct scattering, might fail if the scattering is affected by the presence of a van der Waals well and that theory should attempt to model surface atom motion for reproducing diffraction experiments performed for surface temperatures of 500 K and higher.

Diffraction of CH4 from a Metal Surface

Diffraction with matter waves has been reported since the beginning of quantum mechanics. In free space, diffraction effects have been observed even with objects as large as C₆₀ molecules. However, in scattering from a solid surface, pure elastic diffraction features have never been observed with molecules larger than D₂. Here we report the observation of pure molecular diffraction for CH₄ scattered off of an Ir(111) surface. These results prove that quantum coherence is preserved, despite the small separation between rotational levels and the interaction with surface phonons. Density functional theory calculations of the potential energy surface provide some clues to understand the larger corrugation sampled by CH₄ molecules in comparison to Ne atoms. Our results show that isotope separation of polyatomic molecules may be possible using gas-surface diffraction.

Performance of van der Waals DFT approaches for helium diffraction on metal surfaces

The ability of the different approaches proposed to date to include the effects of van der Waals (vdW) dispersion forces in density functional theory (DFT) is currently under debate. Here, we used the diffraction of He on a Ru(0 0 0 1) surface as a challenging benchmark system to analyze the suitability of several representative approaches, from the ones correcting the exchange-correlation generalized gradient approximation (GGA) functional, to the ones correcting the DFT energies through pairwise-based methods. To perform our analysis, we have built seven continuous potential energy surfaces (PESs) and carried out quantum dynamics simulations using a multi-configuration time-dependent Hartree method. Our analysis reveals that standard DFT within the PBE-GGA framework, although it overestimates diffraction probabilities, yields the best results in comparison with available experimental measurements. On the other hand, although several of the existing vdW DFT approaches yield physisorption wells in very good agreement with experiment, they all seem to overestimate the long-distance corrugation of the PES, the region probed by He scattering, resulting in a large overestimation of diffraction probabilities.

Test of the Transferability of the Specific Reaction Parameter Functional for H2 + Cu(111) to D2 + Ag(111)

The accurate description of the dissociative chemisorption of a molecule on a metal surface requires a chemically accurate description of the molecule-surface interaction. Previously, it was shown that the specific reaction parameter approach to density functional theory (SRP-DFT) enables accurate descriptions of the reaction of dihydrogen with metal surfaces in, for instance, H₂ + Pt(111), H₂ + Cu(111), and H₂ + Cu(100). SRP-DFT likewise allowed a chemically accurate description of dissociation of methane on Ni(111) and Pt(111), and the SRP functional for CH₄ + Ni(111) was transferable to CH₄ + Pt(111), where Ni and Pt belong to the same group. Here, we investigate whether the SRP density functional derived for H₂ + Cu(111) also gives chemically accurate results for H₂ + Ag(111), where Ag belongs to the same group as Cu. To do this, we have performed quasi-classical trajectory calculations using the six-dimensional potential energy surface of H₂ + Ag(111) within the Born-Oppenheimer static surface approximation. The computed reaction probabilities are compared with both state-resolved associative desorption and molecular beam sticking experiments. Our results do not yet show transferability, as the computed sticking probabilities and initial-state selected reaction probabilities are shifted relative to experiment to higher energies by about 2-3 kcal/mol. The lack of transferability may be due to the different character of the SRP functionals for H₂ + Cu and CH₄ + group 10 metals, the latter containing a van der Waals correlation functional and the former not.

Editorial for PCCP themed issue "Developments in Density Functional Theory"

This issue provides an overview of the state-of-the-art of DFT, ranging from mathematical and software developments, via topics in chemical bonding theory, to all kinds of molecular and material properties. Through this issue, we also celebrate the enormous contributions that Evert Jan Baerends has made to this field.

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.

Energy dissipation to tungsten surfaces upon hot-atom and Eley–Rideal recombination of H2

Adiabatic and nonadiabatic quasi-classical molecular dynamics simulations are performed to investigate the role of electron–hole pair excitations in hot-atom and Eley–Rideal H₂ recombination mechanisms on H-covered W(100). The influence of the surface structure is analyzed by comparing with previous results for W(110).

Misconceptions in the Exploding Flask Demonstration Resolved through Students' Critical Thinking

As it connects to a large set of important fundamental ideas in chemistry and analytical techniques discussed in high school chemistry curricula, we review the exploding flask demonstration. In this demonstration, methanol vapor is catalytically oxidized by a Pt wire catalyst in an open container. The exothermicity of reactions occurring at the catalytic surface heats the metal to the extent that it glows. When restricting reactant and product gas flow, conditions may favor repetitive occurrence of a small explosion. We show how mass spectrometry and infrared spectroscopy allow for unravelling the chemical background of this demonstration and discuss various ideas on how to use it in a classroom setting to engage students' critical thinking about chemical research. Along the way, we show that two commonly published ideas about the chemical background of this demonstration are incorrect, and we suggest simple tests that may be performed in a high school setting either as an addition to the demonstration or as a student research project.

Reactivity of O2 on Pd/Ru(0001) and PdRu/Ru(0001) surface alloys

The reactivity of a Pd monolayer epitaxially grown on Ru(0001) toward O2 has been investigated by molecular beam techniques. O2 initial sticking coefficients were determined using the King and Wells method in the incident energy range of 40-450 meV and for sample temperatures of 100 K and 300 K, and compared to the corresponding values measured on the clean Ru(0001) and Pd(111) surfaces. In contrast to the high reactivity shown by Ru(0001) at 100 K, the Pd/Ru(0001) system exhibits a monotonic decrease in the sticking probability of O2 as a function of normal incident energy. At room temperature, the system was found to be inert. Thermal desorption measurements show that O2 is adsorbed molecularly at 100 K. A completely different behaviour has been measured for the Pd0.95Ru0.05/Ru(0001) surface alloy. On this surface, the O2 sticking probability increases with incident energy and resembles the one observed on the clean Ru(0001) surface, even at 300 K. Thermal desorption measurements point to dissociative adsorption of O2 in this system. Both the charge transfer from the Pd to the Ru substrate and the compressive strain on the Pd monolayer contribute to decrease in the reactivity of the Pd/Ru(0001) system well below those of both Ru(0001) and Pd(111).

Dissipation of the excess energy of the adsorbate-thermalization via electron transfer

A new scenario for the thermalization process of adsorbates at solid surfaces is proposed. The scenario is based on the existence of an electric dipole layer in which the electron wavefunctions extend over the positive ions, creating a strong local electric field which drags the electrons into the solid interior and repels the positive ions. During adsorption the electrons tunnel into the solid interior, conveying the excess energy. The positive ions are retarded by the field, losing the excess kinetic energy, and are located smoothly into the adsorption sites. In such a scheme, the excess energy is not dissipated locally, avoiding melting or the creation of defects which is in accordance with experiments. The scenario is supported by ab initio calculation results, including density function theory of the slabs representing the AlN surface and the Schrodinger equation for the time evolution of hydrogen-like atoms at the solid surface.

Energy Dissipation during Diffusion at Metal Surfaces: Disentangling the Role of Phonons versus Electron-Hole Pairs

Helium spin echo experiments combined with ab initio based Langevin molecular dynamics simulations are used to quantify the adsorbate-substrate coupling during the thermal diffusion of Na atoms on Cu(111). An analysis of trajectories within the local density friction approximation allows the contribution from electron-hole pair excitations to be separated from the total energy dissipation. Despite the minimal electronic friction coefficient of Na and the relatively small mass mismatch to Cu promoting efficient phononic dissipation, about (20±5)% of the total energy loss is attributable to electronic friction. The results suggest a significant role of electronic nonadiabaticity in the rapid thermalization generally relied upon in adiabatic diffusion theories.

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.

Two distinctive energy migration pathways of monolayer molecules on metal nanoparticle surfaces

Energy migrations at metal nanomaterial surfaces are fundamentally important to heterogeneous reactions. Here we report two distinctive energy migration pathways of monolayer adsorbate molecules on differently sized metal nanoparticle surfaces investigated with ultrafast vibrational spectroscopy. On a 5 nm platinum particle, within a few picoseconds the vibrational energy of a carbon monoxide adsorbate rapidly dissipates into the particle through electron/hole pair excitations, generating heat that quickly migrates on surface. In contrast, the lack of vibration-electron coupling on approximately 1 nm particles results in vibrational energy migration among adsorbates that occurs on a twenty times slower timescale. Further investigations reveal that the rapid carbon monoxide energy relaxation is also affected by the adsorption sites and the nature of the metal but to a lesser extent. These findings reflect the dependence of electron/vibration coupling on the metallic nature, size and surface site of nanoparticles and its significance in mediating energy relaxations and migrations on nanoparticle surfaces.

Energy migrations at metal nanomaterial surfaces are fundamentally important to heterogeneous reactions. Here, the authors employ ultrafast vibrational spectroscopy to show two distinctive energy migration pathways of monolayer adsorbate molecules on differently sized metal nanoparticle surfaces.

Electron-Hole Pair Effects in Polyatomic Dissociative Chemisorption: Water on Ni(111)

The influence of electron-hole pairs in dissociative chemisorption of a polyatomic molecule (water) on metal surfaces is assessed for the first time using a friction approach. The atomic local density dependent friction coefficients computed based on a free electron gas embedding model are employed in classical molecular dynamics simulations of the water dissociation dynamics on rigid Ni(111) using a recently developed nine dimensional interaction potential energy surface for the system. The results indicate that nonadiabatic effects are relatively small and they do not qualitatively alter the mode specificity in the dissociation.

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.

Chemical energy dissipation at surfaces under UHV and high pressure conditions studied using metal–insulator–metal and similar devices

Thin film metal heterostructures have allowed new light to be shed on the dissipation of chemical energy into electric excitations on metal surfaces.

Quantum dynamics of polyatomic dissociative chemisorption on transition metal surfaces: mode specificity and bond selectivity

Recent advances in quantum dynamical characterization of polyatomic dissociative chemisorption on accurate global potential energy surfaces are critically reviewed.

Toward a Database of Chemically Accurate Barrier Heights for Reactions of Molecules with Metal Surfaces

Being able to calculate reaction barrier heights to within chemical accuracy (errors < 1 kcal/mol) is crucial to the accurate modeling of chemical reactions. Although accurate databases exist that can help theorists with benchmarking new electronic structure theories on gas-phase chemical reactions, no such databases exist for reactions of molecules with metal surfaces. Nonetheless, most chemicals are made in heterogeneously catalyzed processes, of which many take place over metal particles. Presently, barrier heights for molecule-metal surface reactions have been determined with chemical accuracy for only two systems, that is, H2 + Cu(111) and H2 + Cu(100). This has been done with semiempirically determined density functionals, which were fitted through comparisons of dynamics results with molecular beam sticking probabilities. The prospects of extending the database with chemically accurate data for other molecule-metal reactions, either with the use of semiempirical density functional theory or with first-principles theory, are discussed.

The effect of phonon modes and electron–hole pair couplings on molecule–surface scattering processes

The effect of phonon modes and electron–hole pair (elhp) couplings at different surface temperature on D 2(v = 0, 1; j = 0)– Cu (111) collision has been explored by assuming weakly correlated interactions between molecular Degrees of Freedoms (DOFs) with surface modes and elhp excitations through a Hartree product type wavefunction, where the initial state distributions for the phonon modes and the elhp couplings are incorporated by using Bose–Einstein and Fermi–Dirac probability factors, respectively. We carry out four (4D⊗2D)- and six (6D)- dimensional quantum dynamics on such an effective Hamiltonian, and depict the calculated sticking/transition probabilities and energy transfer from molecule to the surface. The phonon modes slightly affect the sticking probability by broadening the profile, but the transition probability are substantially changed with respect to the rigid surface. On the contrary, the inclusion of elhp coupling along with phonon modes does not change the results much compared to the only phonon case.

Semiclassical multi-phonon theory for atom-surface scattering: Application to the Cu(111) system

The semiclassical perturbation theory of Hubbard and Miller [J. Chem. Phys. 80, 5827 (1984)] is further developed to include the full multi-phonon transitions in atom-surface scattering. A practically applicable expression is developed for the angular scattering distribution by utilising a discretized bath of oscillators, instead of the continuum limit. At sufficiently low surface temperature good agreement is found between the present multi-phonon theory and the previous one-, and two-phonon theory derived in the continuum limit in our previous study [Daon, Pollak, and Miret-Artés, J. Chem. Phys. 137, 201103 (2012)]. The theory is applied to the measured angular distributions of Ne, Ar, and Kr scattered from a Cu(111) surface. We find that the present multi-phonon theory substantially improves the agreement between experiment and theory, especially at the higher surface temperatures. This provides evidence for the importance of multi-phonon transitions in determining the angular distribution as the surface temperature is increased.

Unveiling mode-selected electron-phonon interactions in metal films by helium atom scattering

The quasi two-dimensional electron gas on a metal film can transmit to the surface even minute mechanical disturbances occurring in the depth, thus allowing the gentlest of all surface probes, helium atoms, to perceive the vibrations of the deepest atoms via the induced surface-charge density oscillations. A density functional perturbation theory (DFPT) and a helium atom scattering study of the phonon dispersion curves in lead films of up to 7 mono-layers on a copper substrate show that: (a) the electron-phonon interaction is responsible for the coupling of He atoms to in-depth phonon modes; and (b) the inelastic HAS intensity from a given phonon mode is proportional to its electron-phonon coupling. The direct determination of mode-selected electron-phonon coupling strengths has great relevance for understanding superconductivity in thin films and two-dimensional systems.

Role of physisorption states in molecular scattering: a semilocal density-functional theory study on O2/Ag(111)

We simulate the scattering of O2 from Ag(111) with classical dynamics simulations performed on a six-dimensional potential energy surface calculated within semilocal density-functional theory. The enigmatic experimental trends that originally required the conjecture of two types of repulsive walls, arising from a physisorption and chemisorption part of the interaction potential, are fully reproduced. Given the inadequate description of the physisorption properties in semilocal density-functional theory, our work casts severe doubts on the prevalent notion to use molecular scattering data as indirect evidence for the existence of such states.

Electronic friction dominates hydrogen hot-atom relaxation on Pd(100)

We study the dynamics of transient hot H atoms on Pd(100) that originated from dissociative adsorption of H2. The methodology developed here, denoted AIMDEF, consists of ab initio molecular dynamics simulations that include a friction force to account for the energy transfer to the electronic system. We find that the excitation of electron-hole pairs is the main channel for energy dissipation, which happens at a rate that is five times faster than energy transfer into Pd lattice motion. Our results show that electronic excitations may constitute the dominant dissipation channel in the relaxation of hot atoms on surfaces.

Transient quantum trapping of fast atoms at surfaces

We report on the experimental observation and theoretical study of the bound state resonances in fast atom diffraction at surfaces. In our studies, the 4He atom beam has been scattered from a high-quality LiF(001) surface at very small grazing incidence angles. In this regime, the reciprocal lattice vector exchange with the surface allows transient trapping of the 0.3-0.5 keV projectiles into the quasistationary states bound by the attractive atom-surface potential well which is only 10 meV deep. Analysis of the linewidths of the calculated and measured resonances reveals that prior to their release, the trapped projectiles preserve their coherence over travel distances along the surface as large as 0.2  μm, while being in average only at some angstroms in front of the last atomic plane.

Six-dimensional quantum dynamics for dissociative chemisorption of H2 and D2 on Ag(111) on a permutation invariant potential energy surface

Quantum dynamics on a permutation invariant potential energy surface for H₂ dissociation on Ag(111) yield satisfactory agreement with experiment.