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

Spiers Memorial Lecture: New directions in molecular scattering

The field of molecular scattering is reviewed as it pertains to gas-gas as well as gas-surface chemical reaction dynamics. We emphasize the importance of collaboration of experiment and theory, from which new directions of research are being pursued on increasingly complex problems. We review both experimental and theoretical advances that provide the modern toolbox available to molecular-scattering studies. We distinguish between two classes of work. The first involves simple systems and uses experiment to validate theory so that from the validated theory, one may learn far more than could ever be measured in the laboratory. The second class involves problems of great complexity that would be difficult or impossible to understand without a partnership of experiment and theory. Key topics covered in this review include crossed-beams reactive scattering and scattering at extremely low energies, where quantum effects dominate. They also include scattering from surfaces, reactive scattering and kinetics at surfaces, and scattering work done at liquid surfaces. The review closes with thoughts on future promising directions of research.

First-Principles Nonadiabatic Dynamics of Molecules at Metal Surfaces with Vibrationally Coupled Electron Transfer

Accurate description of nonadiabatic dynamics of molecules at metal surfaces involving electron transfer has been a long-standing challenge for theory. Here, we tackle this problem by first constructing high-dimensional neural network diabatic potentials including state crossings determined by constrained density functional theory, then applying mixed quantum-classical surface hopping simulations to evolve coupled electron-nuclear motion. Our approach accurately describes the nonadiabatic effects in CO scattering from Au(111) without empirical parameters and yields results agreeing well with experiments under various conditions for this benchmark system. We find that both adiabatic and nonadiabatic energy loss channels have important contributions to the vibrational relaxation of highly vibrationally excited CO(v_{i}=17), whereas relaxation of low vibrationally excited states of CO(v_{i}=2) is weak and dominated by nonadiabatic energy loss. The presented approach paves the way for accurate first-principles simulations of electron transfer mediated nonadiabatic dynamics at metal surfaces.

Electronic friction near metal surface: Incorporating nuclear quantum effect with ring polymer molecular dynamics

The molecular dynamics with electronic friction (MDEF) approach can accurately describe nonadiabatic effects at metal surfaces in the weakly nonadiabatic limit. That being said, the MDEF approach treats nuclear motion classically such that the nuclear quantum effects are completely missing in the approach. To address this limitation, we combine Electronic Friction with Ring Polymer Molecular Dynamics (EF-RPMD). In particular, we apply the averaged electronic friction from the metal surface to the centroid mode of the ring polymer. We benchmark our approach against quantum dynamics to show that EF-RPMD can accurately capture zero-point energy as well as transition dynamics. In addition, we show that EF-RPMD can correctly predict the electronic transfer rate near metal surfaces in the tunneling limit as well as the barrier crossing limit. We expect that our approach will be very useful to study nonadiabatic dynamics near metal surfaces when nuclear quantum effects become essential.

Sustained Hydrogen Spillover on Pt/Cu(111) Single-Atom Alloy: Dynamic Insights into Gas-Induced Chemical Processes

Hydrogen spillover, involving the surface migration of dissociated hydrogen atoms from active metal sites to the relatively inert catalyst support, plays a crucial role in hydrogen-involved catalytic processes. However, a comprehensive understanding of how H atoms are driven to spill over from active sites onto the catalyst support is still lacking. Here, we examine the atomic-scale perspective of the H spillover process on a Pt/Cu(111) single atom alloy surface using machine-learning accelerated molecular dynamics calculations based on density functional theory. Our results show that when an impinging H2 dissociates at an active Pt site, the Pt atom undergoes deactivation due to the dissociated hydrogen atoms that attach to it. Interestingly, collisions between H2 and sticking H atoms facilitate H spillover onto the host Cu, leading to the reactivation of the Pt atom and the realization of a continuous H spillover process. This work underscores the importance of the interaction between gas molecules and adsorbates as a driving force in elucidating chemical processes under a gaseous atmosphere, which has so far been underappreciated in thermodynamic studies.

Modeling Equilibration Dynamics of Hyperthermal Products of Surface Reactions Using Scalable Neural Network Potential with First-Principles Accuracy

Equilibration dynamics of hot oxygen atoms following the dissociation of O2 on Pd(100) and Pd(111) surfaces are investigated by molecular dynamics simulations based on a scalable neural network potential enabling first-principles description of the interaction of O2 and O interacting with variable Pd supercells. By analyzing hundreds of trajectories with appropriate initial sampling, the measured distance distribution of equilibrated atom pairs on Pd(111) is well reproduced. However, our results on Pd(100) suggest that the ballistic motion of hot atoms predicted previously is a rare event under ideal conditions, while initial molecular orientation and surface thermal fluctuation could significantly affect the overall postdissociation dynamics. On both surfaces, dissociated hyperthermal oxygen atoms primarily locate their nascent positions and experience a similar random walk motion nearby.

Disentangling the role of electrons and phonons in the photoinduced CO desorption and CO oxidation on (O,CO)-Ru(0001)

The role played by electronic and phononic excitations in the femtosecond laser induced desorption and oxidation of CO coadsorbed with O on Ru(0001) is investigated using ab initio molecular dynamics with electronic friction. To this aim, simulations that account for both kind of excitations and that only consider electronic excitations are performed. Results for three different surface coverages are obtained. We unequivocally demonstrate that CO desorption is governed by phononic excitations. In the case of oxidation the low statistics does not allow to give a categorical answer. However, the analysis of the adsorbates kinetic energy gain and displacements strongly suggest that phononic excitations and surface distortion also play an important role in the oxidation process.

Why Ultrafast Photoinduced CO Desorption Dominates over Oxidation on Ru(0001)

CO oxidation on Ru(0001) is a long-standing example of a reaction that, being thermally forbidden in ultrahigh vacuum, can be activated by femtosecond laser pulses. In spite of its relevance, the precise dynamics of the photoinduced oxidation process as well as the reasons behind the dominant role of the competing CO photodesorption remain unclear. Here we use ab initio molecular dynamics with electronic friction that account for the highly excited and nonequilibrated system created by the laser to investigate both reactions. Our simulations successfully reproduce the main experimental findings: the existence of photoinduced oxidation and desorption, the large desorption to oxidation branching ratio, and the changes in the O K-edge X-ray absorption spectra attributed to the initial stage of the oxidation process. Now, we are able to monitor in detail the ultrafast CO desorption and CO oxidation occurring in the highly excited system and to disentangle what causes the unexpected inertness to the otherwise energetically favored oxidation.

Time-dependent density functional theory calculations of electronic friction in non-homogeneous media

The excitation of low-energy electron-hole pairs is one of the most relevant processes in the gas-surface interaction. An efficient tool to account for these excitations in simulations of atomic and molecular dynamics at surfaces is the so-called local density friction approximation (LDFA). The LDFA is based on a strong approximation that simplifies the dynamics of the electronic system: a local friction coefficient is defined using the value of the electronic density for the unperturbed system at each point of the dynamics. In this work, we apply real-time time-dependent density functional theory to the problem of the electronic friction of a negative point charge colliding with spherical jellium metal clusters. Our non-adiabatic, parameter-free results provide a benchmark for the widely used LDFA approximation and allow the discussion of various processes relevant to the electronic response of the system in the presence of the projectile.

H atom scattering from W(110): A benchmark for molecular dynamics with electronic friction

Molecular dynamics with electronic friction (MDEF) at the level of the local density friction approximation (LDFA) has been applied to describe electronically non-adiabatic energy transfer accompanying H atom collisions with many solid metal surfaces. When implemented with full dimensional potential energy and electron density functions, excellent agreement with experiment is found. Here, we compare the performance of a reduced dimensional MDEF approach involving a simplified description of H atom coupling to phonons to that of full dimensional MDEF calculations known to yield accurate results. Both approaches give remarkably similar results for H atom energy loss distributions with a 300 K W(110) surface. At low surface temperature differences are seen; but, quantities like average energy loss are still accurately reproduced. Both models predict similar conditions under which H atoms that have penetrated into the subsurface regions could be observed in scattering experiments.

Molecular dynamics with electronic friction (MDEF) at the level of the local density friction approximation (LDFA) has been applied to describe electronically non-adiabatic energy transfer accompanying H atom collisions with many solid metal surfaces.

Dissipative tunneling rates through the incorporation of first-principles electronic friction in instanton rate theory. I. Theory

Reactions involving adsorbates on metallic surfaces and impurities in bulk metals are ubiquitous in a wide range of technological applications. The theoretical modeling of such reactions presents a formidable challenge for theory because nuclear quantum effects (NQEs) can play a prominent role and the coupling of the atomic motion with the electrons in the metal gives rise to important non-adiabatic effects (NAEs) that alter atomic dynamics. In this work, we derive a theoretical framework that captures both NQEs and NAEs and, due to its high efficiency, can be applied to first-principles calculations of reaction rates in high-dimensional realistic systems. More specifically, we develop a method that we coin ring polymer instanton with explicit friction (RPI-EF), starting from the ring polymer instanton formalism applied to a system-bath model. We derive general equations that incorporate the spatial and frequency dependence of the friction tensor and then combine this method with the ab initio electronic friction formalism for the calculation of thermal reaction rates. We show that the connection between RPI-EF and the form of the electronic friction tensor presented in this work does not require any further approximations, and it is expected to be valid as long as the approximations of both underlying theories remain valid.

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.

Effective medium theory for bcc metals: electronically non-adiabatic H atom scattering in full dimensions

We report a newly derived Effective Medium Theory (EMT) formalism for bcc metals and apply it for the construction of a full-dimensional PES for H atoms interacting with molybdenum (Mo) and tungsten (W). We construct PESs for the (111) and (110) facets of both metals. The EMT-PESs have the advantage that they automatically provide the background electron density on the fly which makes incorporation of ehp excitation within the framework of electronic friction straightforward. Using molecular dynamics with electronic friction (MDEF) with these new PESs, we simulated 2.76 eV H atoms scattering and adsorption. The large energy losses at a surface temperature of 300 K is very similar those seen for H atom scattering from the late fcc metals and is dominated by ehp excitation. We see significant differences in the scattering from different surface facets of the same metal. For the (110) facet, we see strong evidence of sub-surface scattering, which should be observable in experiment and we predict the best conditions for observing this novel type of scattering process. At low temperatures the MD simulations predict that H atom scattering is surface specific due to the reduced influence of the random force.

We derive a many-body formalism for interaction energies of adsorbates in metals and use it for. electronically non-adiabatic H atom scattering simulations from metal surfaces.

Dynamics of Initial Hydrogen Spillover from a Single Atom Platinum Active Site to the Cu(111) Host Surface: The Impact of Substrate Electron-Hole Pairs

The initial impulsive diffusion of hot hydrogen atoms resulted from the dissociative chemisorption of H₂ at atomically dispersed Pt atoms embedded in Cu(111) is investigated using ab initio molecular dynamics. Upon dissociation, one of the two hydrogen atoms tends to roam away from the dissociation site while the other remains trapped. It is shown that the fraction of diffusion and the average diffusion length increase with the incident energy and H₂ vibrational excitation, due apparently to the increased initial kinetic energy of the hot atoms. Most importantly, the strong interaction with surface electron-hole pairs, modeled using an electronic friction model, is shown to play an important role in rapid energy dissipation and significant retardation of the impulsive diffusion.

Photoinduced Desorption Dynamics of CO from Pd(111): A Neural Network Approach

Modeling the ultrafast photoinduced dynamics and reactivity of adsorbates on metals requires including the effect of the laser-excited electrons and, in many cases, also the effect of the highly excited surface lattice. Although the recent ab initio molecular dynamics with electronic friction and thermostats, (T_e,T_l)-AIMDEF [AlducinM.;Phys. Rev. Lett.2019, 123, 246802]31922860, enables such complex modeling, its computational cost may limit its applicability. Here, we use the new embedded atom neural network (EANN) method [ZhangY.;J. Phys. Chem. Lett.2019, 10, 496231397157] to develop an accurate and extremely complex potential energy surface (PES) that allows us a detailed and reliable description of the photoinduced desorption of CO from the Pd(111) surface with a coverage of 0.75 monolayer. Molecular dynamics simulations performed on this EANN-PES reproduce the (T_e,T_l)-AIMDEF results with a remarkable level of accuracy. This demonstrates the outstanding performance of the obtained EANN-PES that is able to reproduce available density functional theory (DFT) data for an extensive range of surface temperatures (90–1000 K); a large number of degrees of freedom, those corresponding to six CO adsorbates and 24 moving surface atoms; and the varying CO coverage caused by the abundant desorption events.

Multibounce and Subsurface Scattering of H Atoms Colliding with a van der Waals Solid

We report the results of inelastic differential scattering experiments and full-dimensional molecular dynamics trajectory simulations for 2.76 eV H atoms colliding at a surface of solid xenon. The interaction potential is based on an effective medium theory (EMT) fit to density functional theory (DFT) energies. The translational energy-loss distributions derived from experiment and theory are in excellent agreement. By analyzing trajectories, we find that only a minority of the scattering results from simple single-bounce dynamics. The majority comes from multibounce collisions including subsurface scattering where the H atoms penetrate below the first layer of Xe atoms and subsequently re-emerge to the gas phase. This behavior leads to observable energy-losses as large as 0.5 eV, much larger than a prediction of the binary collision model (0.082 eV), which is often used to estimate the highest possible energy-loss in direct inelastic surface scattering. The sticking probability computed with the EMT-PES (0.15) is dramatically reduced (5 × 10-6) if we employ a full-dimensional potential energy surface (PES) based on Lennard-Jones (LJ) pairwise interactions. Although the LJ-PES accurately describes the interactions near the H-Xe and Xe-Xe energy minima, it drastically overestimates the effective size of the Xe atom seen by the colliding H atom at incidence energies above about 0.1 eV.

Random Force in Molecular Dynamics with Electronic Friction

Originally conceived to describe thermal diffusion, the Langevin equation includes both a frictional drag and a random force, the latter representing thermal fluctuations first seen as Brownian motion. The random force is crucial for the diffusion problem as it explains why friction does not simply bring the system to a standstill. When using the Langevin equation to describe ballistic motion, the importance of the random force is less obvious and it is often omitted, for example, in theoretical treatments of hot ions and atoms interacting with metals. Here, friction results from electronic nonadiabaticity (electronic friction), and the random force arises from thermal electron-hole pairs. We show the consequences of omitting the random force in the dynamics of H-atom scattering from metals. We compare molecular dynamics simulations based on the Langevin equation to experimentally derived energy loss distributions. Despite the fact that the incidence energy is much larger than the thermal energy and the scattering time is only about 25 fs, the energy loss distribution fails to reproduce the experiment if the random force is neglected. Neglecting the random force is an even more severe approximation than freezing the positions of the metal atoms or modelling the lattice vibrations as a generalized Langevin oscillator. This behavior can be understood by considering analytic solutions to the Ornstein-Uhlenbeck process, where a ballistic particle experiencing friction decelerates under the influence of thermal fluctuations.

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.

Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications

This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.

Electrons and Phonons Cooperate in the Laser-Induced Desorption of CO from Pd(111)

Femtosecond laser induced desorption of CO from a CO-covered Pd(111) surface is investigated with ab initio molecular dynamics with electronic friction that incorporates effects due to the excited electronic and phononic systems, as well as out-of-phase coadsorbate interactions. Our simulations show evidence of an important electron-phonon synergy in promoting CO desorption that has largely been neglected in other similar systems. At the saturated coverage of 0.75 ML, effects due to CO-CO interadsorbate energy exchange are also important. Our dynamics simulations, in concert with site-specific desorption energy calculations, allow us to understand the large coverage dependence of the desorption yields observed in experiments.

Quantum Dynamics of Dissociative Chemisorption of H2 on the Stepped Cu(211) Surface

Reactions on stepped surfaces are relevant to heterogeneous catalysis, in which a reaction often takes place at the edges of nanoparticles where the edges resemble steps on single-crystal stepped surfaces. Previous results on H₂ + Cu(211) showed that, in this system, steps do not enhance the reactivity and raised the question of whether this effect could be, in any way, related to the neglect of quantum dynamical effects in the theory. To investigate this, we present full quantum dynamical molecular beam simulations of sticking of H₂ on Cu(211), in which all important rovibrational states populated in a molecular beam experiment are taken into account. We find that the reaction of H₂ with Cu(211) is very well described with quasi-classical dynamics when simulating molecular beam sticking experiments, in which averaging takes place over a large number of rovibrational states and over translational energy distributions. Our results show that the stepped Cu(211) surface is distinct from its component Cu(111) terraces and Cu(100) steps and cannot be described as a combination of its component parts with respect to the reaction dynamics when considering the orientational dependence. Specifically, we present evidence that, at translational energies close to the reaction threshold, vibrationally excited molecules show a negative rotational quadrupole alignment parameter on Cu(211), which is not found on Cu(111) and Cu(100). The effect arises because these molecules react with a site-specific reaction mechanism at the step, that is, inelastic rotational enhancement, which is only effective for molecules with a small absolute value of the magnetic rotation quantum number. From a comparison to recent associative desorption experiments as well as Born-Oppenheimer molecular dynamics calculations, it follows that the effects of surface atom motion and electron-hole pair excitation on the reactivity fall within chemical accuracy, that is, modeling these effect shifts extracted reaction probability curves by less than 1 kcal/mol translational energy. We found no evidence in our fully state-resolved calculations for the "slow" reaction channel that was recently reported for associative desorption of H₂ from Cu(111) and Cu(211), but our results for the fast channel are in good agreement with the experiments on H₂ + Cu(211).

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.

Bridging the Gap between Direct Dynamics and Globally Accurate Reactive Potential Energy Surfaces Using Neural Networks

Direct dynamics simulations become increasingly popular in studying reaction dynamics for complex systems where analytical potential energy surfaces (PESs) are unavailable. Yet, the number and/or the propagation time of trajectories are often limited by high computational costs, and numerous energies and forces generated on-the-fly become wasted after simulations. We demonstrate here an example of reusing only a very small portion of existing direct dynamics data to reconstruct a 90-dimensional globally accurate reactive PES describing the interaction of CO₂ with a movable Ni(100) surface based on a machine learning approach. In addition to reproducing previous results with much better statistics, we predict scattering probabilities of CO₂ at the state-to-state level, which is extremely demanding for direct dynamics. We propose this unified way to investigate gaseous and gas-surface reactions of medium size, initiating with hundreds of preliminary direct dynamics trajectories, followed by low-cost and high-quality simulations on full-dimensional analytical PESs.

Reactive and Nonreactive Scattering of HCl from Au(111): An Ab Initio Molecular Dynamics Study

The HCl + Au(111) system has recently become a benchmark for highly activated dissociative chemisorption, which presumably is strongly affected by electron-hole pair excitation. Previous dynamics calculations, which were based on density functional theory at the generalized gradient approximation level (GGA-DFT) for the molecule-surface interaction, have all overestimated measured reaction probabilities by at least an order of magnitude. Here, we perform ab initio molecular dynamics (AIMD) and AIMD with electronic friction (AIMDEF) calculations employing a density functional that includes the attractive van der Waals interaction. Our calculations model the simultaneous and possibly synergistic effects of surface temperature, surface atom motion, electron-hole pair excitation, the molecular beam conditions of the experiments, and the van der Waals interaction on the reactivity. We find that reaction probabilities computed with AIMDEF and the SRP32-vdW functional still overestimate the measured reaction probabilities, by a factor 18 for the highest incidence energy at which measurements were performed (≈2.5 eV). Even granting that the experiment could have underestimated the sticking probability by about a factor three, this still translates into a considerable overestimation of the reactivity by the current theory. Likewise, scaled transition probabilities for vibrational excitation from ν = 1, j = 1 to ν = 2 are overestimated by the AIMDEF theory, by factors 3-8 depending on the initial conditions modeled. Energy losses to the surface and translational energy losses are, however, in good agreement with experimental values.

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.

Hot electron effects during reactive scattering of H2 from Ag(111): assessing the sensitivity to initial conditions, coupling magnitude, and electronic temperature

We use an analytical representation of electronic friction for H₂ on Ag(111) to assess the validity and robustness of the MDEF method based on TDPT.

Electron-Mediated Phonon-Phonon Coupling Drives the Vibrational Relaxation of CO on Cu(100)

We bring forth a consistent theory for the electron-mediated vibrational intermode coupling that clarifies the microscopic mechanism behind the vibrational relaxation of adsorbates on metal surfaces. Our analysis points out the inability of state-of-the-art nonadiabatic theories to quantitatively reproduce the experimental linewidth of the CO internal stretch mode on Cu(100) and it emphasizes the crucial role of the electron-mediated phonon-phonon coupling in this regard. The results demonstrate a strong electron-mediated coupling between the internal stretch and low-energy CO modes, but also a significant role of surface motion. Our nonadiabatic theory is also able to explain the temperature dependence of the internal stretch phonon linewidth, thus far considered a sign of the direct anharmonic coupling.

Characterization of the Platinum-Hydrogen Bond by Surface-Sensitive Time-Resolved Infrared Spectroscopy

The vibrational dynamics of Pt-H on a nanostructured platinum surface has been examined by ultrafast infrared spectroscopy. Three bands are observed at 1800, 2000, and 2090 cm^-1, which are assigned to Pt-CO in a bridged and linear configuration and Pt-H, respectively. Lifetime analysis revealed a time constant of (0.8 ± 0.1) ps for the Pt-H mode, considerably shorter than that of Pt-CO because of its stronger coupling to the metal substrate. Two-dimensional attenuated total reflection infrared spectroscopy provided additional evidence for the assignment based on the anharmonic shift, which is large in the case of Pt-H (90 cm^-1), in agreement with the density functional theory calculations. The absorption cross section of Pt-H is smaller than that of the very strong Pt-CO vibration by only a modest factor of ∼1.5-3. Because Pt-H is transiently involved in catalytic water splitting on Pt, the present spectroscopic characterization paves the way for in-operando kinetic studies of such 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.

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

Nonadiabatic Vibrational Damping of Molecular Adsorbates: Insights into Electronic Friction and the Role of Electronic Coherence

We present a perturbation approach rooted in time-dependent density-functional theory to calculate electron-hole (e-h) pair excitation spectra during the nonadiabatic vibrational damping of adsorbates on metal surfaces. Our analysis for the benchmark systems CO on Cu(100) and Pt(111) elucidates the surprisingly strong influence of rather short electronic coherence times. We demonstrate how in the limit of short electronic coherence times, as implicitly assumed in prevalent quantum nuclear theories for the vibrational lifetimes as well as electronic friction, band structure effects are washed out. Our results suggest that more accurate lifetime or chemicurrentlike experimental measurements could characterize the electronic coherence.

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.

Mode Specific Electronic Friction in Dissociative Chemisorption on Metal Surfaces: H_{2} on Ag(111)

Electronic friction and the ensuing nonadiabatic energy loss play an important role in chemical reaction dynamics at metal surfaces. Using molecular dynamics with electronic friction evaluated on the fly from density functional theory, we find strong mode dependence and a dominance of nonadiabatic energy loss along the bond stretch coordinate for scattering and dissociative chemisorption of H_{2} on the Ag(111) surface. Exemplary trajectories with varying initial conditions indicate that this mode specificity translates into modulated energy loss during a dissociative chemisorption event. Despite minor nonadiabatic energy loss of about 5%, the directionality of friction forces induces dynamical steering that affects individual reaction outcomes, specifically for low-incidence energies and vibrationally excited molecules. Mode-specific friction induces enhanced loss of rovibrational rather than translational energy and will be most visible in its effect on final energy distributions in molecular scattering experiments.

Mixed Quantum-Classical Dynamics Using Collective Electronic Variables: A Better Alternative to Electronic Friction Theories

An accurate description of nonadiabatic dynamics of molecular species on metallic surfaces poses a serious computational challenge associated with a multitude of closely spaced electronic states. We propose a mixed quantum-classical scheme that addresses this challenge by introducing collective electronic variables. These variables are defined through analytic block-diagonalization applied to the time-dependent Hamiltonian matrix governing the electronic dynamics. We compare our scheme with a simplified Ehrenfest approach and with a full-memory electronic friction model on a 1D "adatom + atomic chain" model. Our simulations demonstrate that collective-mode dynamics with only a few (two to three) electronic variables is robust and can describe a variety of situations: from a chemisorbed atom on an insulator to an atom on a metallic surface. Our molecular model also reveals that the friction approach is prone to unpredictable and catastrophic failures.

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.

Liquid-phase catalytic reactor combined with measurement of hot electron flux and chemiluminescence

Understanding the role of electronically nonadiabatic interactions during chemical reactions on metal surfaces in liquid media is of great importance for a variety of applications including catalysis, electrochemistry, and environmental science. Here, we report the design of an experimental apparatus for detection of the highly excited (hot) electrons created as a result of nonadiabatic energy transfer during the catalytic decomposition of hydrogen peroxide on thin-film metal-semiconductor nanodiodes. The apparatus enables the measurement of hot electron flows and related phenomena (e.g., surface chemiluminescence) as well as the corresponding reaction rates at different temperatures. The products of the chemical reaction can be characterized in the gaseous phase by means of gas chromatography. The combined measurement of hot electron flux, catalytic activity, and light emission can lead to a fundamental understanding of the elementary processes occurring during the heterogeneous catalytic reaction.

Chemically Accurate Simulation of a Polyatomic Molecule-Metal Surface Reaction

Although important to heterogeneous catalysis, the ability to accurately model reactions of polyatomic molecules with metal surfaces has not kept pace with developments in gas phase dynamics. Partnering the specific reaction parameter (SRP) approach to density functional theory with ab initio molecular dynamics (AIMD) extends our ability to model reactions with metals with quantitative accuracy from only the lightest reactant, H2, to essentially all molecules. This is demonstrated with AIMD calculations on CHD3 + Ni(111) in which the SRP functional is fitted to supersonic beam experiments, and validated by showing that AIMD with the resulting functional reproduces initial-state selected sticking measurements with chemical accuracy (4.2 kJ/mol ≈ 1 kcal/mol). The need for only semilocal exchange makes our scheme computationally tractable for dissociation on transition metals.