Competition between electron and phonon excitations in the scattering of nitrogen atoms and molecules off tungsten and silver metal surfaces

First-Principles Insights into Adiabatic and Nonadiabatic Vibrational Energy-Transfer Dynamics during Molecular Scattering from Metal Surfaces: The Importance of Surface Reactivity

Energy transfer is ubiquitous during molecular collisions and reactions at gas-surface interfaces. Of particular importance is vibrational energy transfer because of its relevance to bond forming and breaking. In this Perspective, we review recent first-principles studies on vibrational energy-transfer dynamics during molecular scattering from metal surfaces at the state-to-state level. Taking several representative systems as examples, we highlight the intrinsic correlation between vibrational energy transfer in nonreactive scattering and surface reactivity and how it operates in both electronically adiabatic and nonadiabatic pathways. Adiabatically, the presence of a dissociation barrier softens a bond in the impinging molecule and increases its couplings with other molecular modes and surface phonons. In the meantime, the stronger interaction between the molecule and the surface also changes the electronic structure at the barrier, resulting in an increase of nonadiabatic effects. We further discuss future prospects toward a more quantitative understanding of this important surface dynamical process.

Mode Selective Chemistry for the Dissociation of Methane on Efficient Ni/Pt-Bimetallic Alloy Catalysts

The mode selectivity of methane dissociation is studied on three different Ni/Pt-bimetallic alloy surfaces using a fully quantum approach based on reaction path Hamiltonian. Dissociative sticking probability depends on the...

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.

Mechanical Vibrational Relaxation of NO Scattering from Metal and Insulator Surfaces: When and Why They Are Different

NO scattering from metallic and insulating surfaces represents contrasting benchmark systems for understanding energy transfer at gas-surface interface. Strikingly different behaviors of highly vibrationally excited NO scattered from Au(111) and LiF(001) were observed and attributed to disparate electronic structures between metals and insulators. Here, we reveal an alternative mechanical origin of this discrepancy by comparative molecular dynamics simulations with globally accurate adiabatic neural network potentials of both systems. We find that highly vibrating NO can reach for the high-dissociation barrier on Au(111), by which vibrational energy can largely transfer to translation or rotation and further dissipate into substrate phonons. This mechanical energy transfer channel is forbidden in the purely repulsive NO/LiF(001) system or for low-vibrating NO on Au(111), where molecular vibration is barely coupled to other degrees of freedom. Our results emphasize that the initial state and potential energy landscape concurrently influence the mechanical energy transfer dynamics of gas-surface scattering.

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.

Dynamics in reactions on metal surfaces: A theoretical perspective

Recent advances in theoretical characterization of reaction dynamics on metal surfaces are reviewed. It is shown that the widely available density functional theory of metals and their interactions with molecules have enabled first principles theoretical models for treating surface reaction dynamics. The new theoretical tools include methods to construct high-dimensional adiabatic potential energy surfaces, to characterize nonadiabatic processes within the electronic friction models, and to describe dynamics both quantum mechanically and classically. Three prototypical surface reactions, namely, dissociative chemisorption, Eley-Rideal reactions, and recombinative desorption, are surveyed with a focus on some representative examples. While principles governing gas phase reaction dynamics may still be applicable, the presence of the surface introduces a higher level of complexity due to strong interaction between the molecular species and metal substrate. Furthermore, most of these reactive processes are impacted by energy exchange with surface phonons and/or electron-hole pair excitations. These theoretical studies help to interpret and rationalize experimental observations and, in some cases, guide experimental explorations. Knowledge acquired in these fundamental studies is expected to impact many practical problems in a wide range of interfacial processes.

Stochastic wave packet approach to nonadiabatic scattering of diatomic molecules from metals

In this contribution, we present a quantum dynamical approach to study inelastic scattering of diatomic molecules from metal surfaces at normal incidence. The dissipative dynamics obeys a stochastic Schrödinger equation describing the time-evolution of the system as a piecewise deterministic process. Energy exchange between the molecular vibrational degrees of freedom and the metal electrons is represented using operators in tensor product form, which are coupled via anharmonic transition rates calculated from first-order perturbation theory. Full dimensional observables are obtained by averaging over simulations in 4D-including the internal stretch, the distance to the surface, and the orientation angles-at different surface sites. The method is applied to the state-resolved scattering of vibrationally excited NO from Au(111), revealing important channels for quantized energy relaxation.

Hot-electron effects during reactive scattering of H2 from Ag(111): the interplay between mode-specific electronic friction and the potential energy landscape

The breakdown of the Born-Oppenheimer approximation gives rise to nonadiabatic effects in gas-surface reactions at metal surfaces. However, for a given reaction, it remains unclear which factors quantitatively determine whether these effects measurably contribute to surface reactivity in catalysis and photo/electrochemistry. Here, we systematically investigate hot electron effects during H2 scattering from Ag(111) using electronic friction theory. We combine first-principles calculations of tensorial friction by time-dependent perturbation theory based on density functional theory and an analytical neural network representation, to overcome the limitations of existing approximations and explicitly simulate mode-specific nonadiabatic energy loss during molecular dynamics. Despite sizable hot-electron-induced energy loss, no measurable nonadiabatic effects can be found for H2 scattering on Ag(111). This is in stark contrast to previous reports for vibrationally excited H2 scattering on Cu(111). By detailed analysis of the two systems, we attribute this discrepancy to a subtle interplay between the magnitude of electronic friction along intramolecular vibration and the shape of the potential energy landscape that controls the molecular velocity at impact. On the basis of this characterization, we offer guidance for the search of highly nonadiabatic surface reactions.

Hydrogen abstraction from metal surfaces: when electron-hole pair excitations strongly affect hot-atom recombination

Using molecular dynamics simulations, we predict that the inclusion of nonadiabatic electronic excitations influences the dynamics of preadsorbed hydrogen abstraction from the W(110) surface by hydrogen scattering. The hot-atom recombination, which involves hyperthermal diffusion of the impinging atom on the surface, is significantly affected by the dissipation of energy mediated by electron-hole pair excitations at low coverage and low incidence energy. This issue is of importance as this abstraction mechanism is thought to largely contribute to molecular hydrogen formation from metal surfaces.

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).

Strong Anisotropic Interaction Controls Unusual Sticking and Scattering of CO at Ru(0001)

Complete sticking at low incidence energies and broad angular scattering distributions at higher energies are often observed in molecular beam experiments on gas-surface systems which feature a deep chemisorption well and lack early reaction barriers. Although CO binds strongly on Ru(0001), scattering is characterized by rather narrow angular distributions and sticking is incomplete even at low incidence energies. We perform molecular dynamics simulations, accounting for phononic (and electronic) energy loss channels, on a potential energy surface based on first-principles electronic structure calculations that reproduce the molecular beam experiments. We demonstrate that the mentioned unusual behavior is a consequence of a very strong rotational anisotropy in the molecule-surface interaction potential. Beyond the interpretation of scattering phenomena, we also discuss implications of our results for the recently proposed role of a precursor state for the desorption and scattering of CO from ruthenium.

Vibrational Excitation of H2 Scattering from Cu(111): Effects of Surface Temperature and of Allowing Energy Exchange with the Surface

In scattering of H₂ from Cu(111), vibrational excitation has so far defied an accurate theoretical description. To expose the causes of the large discrepancies with experiment, we investigate how the feature due to vibrational excitation (the "gain peak") in the simulated time-of-flight spectrum of (v = 1, j = 3) H₂ scattering from Cu(111) depends on the surface temperature (T_s) and the possibility of energy exchange with surface phonons and electron-hole pairs (ehp's). Quasi-classical dynamics calculations are performed on the basis of accurate semiempirical density functionals for the interaction with H₂ + Cu(111). The methods used include the quasi-classical trajectory method within the Born-Oppenheimer static surface model, the generalized Langevin oscillator (GLO) method incorporating energy transfer to surface phonons, the GLO + friction (GLO+F) method also incorporating energy exchange with ehp's, and ab initio molecular dynamics with electronic friction (AIMDEF). Of the quasi-classical methods tested, comparison with AIMDEF suggests that the GLO+F method is accurate enough to describe vibrational excitation as measured in the experiments. The GLO+F calculations also suggest that the promoting effect of raising T_s on the measured vibrational excitation is due to an electronically nonadiabatic mechanism. However, by itself, enabling energy exchange with the surface by modeling surface phonons and ehp's leads to reduced vibrational excitation, further decreasing the agreement with experiment. The simulated gain peak is quite sensitive to energy shifts in calculated vibrational excitation probabilities and to shifts in a specific experimental parameter (the chopper opening time). While the GLO+F calculations allow important qualitative conclusions, comparison to quantum dynamics results suggests that, with the quasi-classical way of describing nuclear motion and the present box quantization method for assigning the final vibrational state, the gain peak is not yet described with quantitative accuracy. Ways in which this problem might be resolved in the future are discussed.

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.

Dissociative and non-dissociative adsorption of O2 on Cu(111) and CuML/Ru(0001) surfaces: adiabaticity takes over

Adiabatic molecular spin-quenching during the approach of O₂ to Cu(111) and Cu_ML/Ru(0001) surfaces.

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.

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.

First principles calculations on nitrogen reactivity on tungsten surfaces

Using the spin-polarized gradient-corrected density functional theory we investigate the adsorption and dissociation of nitrogen molecule and hydrogen nitride radicals on two model surfaces of tungsten, W{1 0 0} and W{1 1 0}. The goal of the investigations is to predict the number and nature of the nitrogenized moieties that could be found on a tungsten surface after nitrogen-seeded hydrogen (deuterium) plasma discharge. The surfaces are considered clean or saturated in hydrogen or nitrogen atoms as they are expected to be after nitrogen seeded hydrogen plasma irradiation. We find that the radicals NH, NH2 and NH3 are dissociated by the catalytic action of the metallic surfaces if they are animated with an initial kinetic energy of less than 2.5 eV. The products of the reaction are a nitrogen atom incorporated into the metal surface and an adsorbed proton. The general conclusion is that nitrogen reduces hydrogen retention on the surface and more generally that the nitrided surfaces should be less reactive than the clean tungsten even if the two surfaces we consider behave quite differently.

Molecular dynamics simulation of O2 adsorption on Ag(110) from first principles electronic structure calculations

State of the art simulations show that the physisorption state could be important for O₂/Ag(110) adsorption.

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.

Electronic Friction-Based Vibrational Lifetimes of Molecular Adsorbates: Beyond the Independent-Atom Approximation

We assess the accuracy of vibrational damping rates of diatomic adsorbates on metal surfaces as calculated within the local-density friction approximation (LDFA). An atoms-in-molecules (AIM) type charge partitioning scheme accounts for intramolecular contributions and overcomes the systematic underestimation of the nonadiabatic losses obtained within the prevalent independent-atom approximation. The quantitative agreement obtained with theoretical and experimental benchmark data suggests the LDFA-AIM scheme as an efficient and reliable approach to account for electronic dissipation in ab initio molecular dynamics simulations of surface chemical reactions.

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.

Electronic energy dissipation during scattering of vibrationally excited molecules at metal surfaces: ab initio simulations for HCl/Al(111)

In this Letter, ab initio molecular dynamics simulations based on time-dependent density functional theory for the electrons and Ehrenfest dynamics for the nuclei are reported that detail the interaction of a vibrating HCl molecule with an Al(111) substrate. The mechanism responsible for the strong electron-hole-pair (EHP)-vibrational coupling in case of highly vibrationally excited molecules is traced back to a large eigenenergy shift of the spz*-like antibonding HCl lowest unoccupied molecular orbital with the bond length. As a consequence of this mechanism, the electronic excitation spectra turn out to be highly asymmetric. The simulations suggest an explanation of how to reconcile a strong EHP-vibrational coupling in case of highly vibrationally excited molecules with the small, but clearly evident, electronic contribution to the v=0 → v=1 vibrational excitation observed experimentally during the scattering of HCl molecules at a hot Au surface by Ran et al. [Phys. Rev. Lett. 98 237601 (2007)].

Surface temperature effects on the dynamics of N2 Eley-Rideal recombination on W(100)

Quasiclassical trajectories simulations are performed to study the influence of surface temperature on the dynamics of a N atom colliding a N-preadsorbed W(100) surface under normal incidence. A generalized Langevin surface oscillator scheme is used to allow energy transfer between the nitrogen atoms and the surface. The influence of the surface temperature on the N(2) formed molecules via Eley-Rideal recombination is analyzed at T = 300, 800, and 1500 K. Ro-vibrational distributions of the N(2) molecules are only slightly affected by the presence of the thermal bath whereas kinetic energy is rather strongly decreased when going from a static surface model to a moving surface one. In terms of reactivity, the moving surface model leads to an increase of atomic trapping cross section yielding to an increase of the so-called hot atoms population and a decrease of the direct Eley-Rideal cross section. The energy exchange between the surface and the nitrogen atoms is semi-quantitatively interpreted by a simple binary collision model.

Towards chemically accurate simulation of molecule-surface reactions

This perspective addresses four challenges facing theorists whose aim is to make quantitatively accurate predictions for reactions of molecules on metal surfaces, and suggests ways of meeting these challenges, focusing on dissociative chemisorption reactions of H(2), N(2), and CH(4). Addressing these challenges is ultimately of practical importance to a more accurate description of overall heterogeneously catalysed reactions, which play a role in the production of more than 90% of man-made chemicals. One challenge is to describe the interaction of a molecule with a metal surface with chemical accuracy, i.e., with errors in reaction barrier heights less than 1 kcal mol(-1). In this framework, the potential of a new implementation of specific reaction parameter density functional theory (SRP-DFT) will be discussed, with emphasis on applications to reaction of H(2) with metal surfaces. Two additional challenges are to come up with improved descriptions of the effects of phonons and electron-hole pairs on reaction of molecules like N(2) on metal surfaces. Phonons can be tackled using sudden approximations in quantum dynamics, and through Ab Initio Molecular Dynamics (AIMD) calculations using classical dynamics. To additionally achieve an accurate description of the effect of electron-hole pair excitation on dissociative chemisorption within a classical dynamics framework, it may be possible to combine AIMD with electronic friction. The fourth challenge we will consider is how to achieve an accurate quantum mechanical description of the dissociative chemisorption of a polyatomic molecule, like methane, on a metal surface. A method of potential interest is the Multi-Configuration Time-Dependent Hartree (MCTDH) method.