Determining the Effect of Hot Electron Dissipation on Molecular Scattering Experiments at Metal Surfaces

Vibrational energy transfer in collisions of molecules with metal surfaces

The Born-Oppenheimer approximation (BOA), which serves as the basis for our understanding of chemical bonding, reactivity and dynamics, is routinely violated for vibrationally inelastic scattering of molecules at metal surfaces. The title-field therefore represents a fascinating challenge to our conventional wisdom calling for new concepts that involve explicit electron dynamics occurring in concert with nuclear motion. Here, we review progress made in this field over the last decade, which has witnessed dramatic advances in experimental methods, thereby providing a much more extensive set of diverse observations than has ever before been available. We first review the experimental methods used in this field and then provide a systematic tour of the vast array of observations that are currently available. We show how these observations - taken together and without reference to computational simulations - lead us to a simple and intuitive picture of BOA failure in molecular dynamics at metal surfaces, one where electron transfer between the molecule and the metal plays a preeminent role. We also review recent progress made in the theory of electron transfer mediated BOA failure in molecule-surface interactions, describing the most important methods and their ability to reproduce experimental observation. Finally, we outline future directions for research and important unanswered questions.

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.

Microcanonical treatment of HCl dissociative chemisorption on Au(111): Reactive dampening through inefficient translational energy coupling and an active surface

Microcanonical unimolecular rate theory is applied to Shirhatti and Wodtke's recent supersonic molecular beam experiments examining the activated dissociative chemisorption of HCl on Au(111). A precursor mediated microcanonical trapping (PMMT) model (where the surface vibrates and HCl rotations, vibration, and translation directed along the surface normal are treated as active degrees of freedom) gave dissociative sticking coefficient predictions that are several orders of magnitude higher than experimental values but in good accord with prior quantum and molecular dynamics simulations. Density functional theory (DFT) electronic structure calculations using the Perdew-Burke-Ernzerhof (PBE) functional served to fix the vibrational frequencies of the reactive transition state and the threshold energy for dissociation, E0 = 72.9 kJ/mol. To explore the possibilities of varying threshold energy, coupling to phonons, and dynamics, a three-parameter [E0, s, ɛn] dynamically biased (d-) PMMT model was fit to the experiments. A dynamical bias was introduced using an efficiency, ɛn, of normal translational energy to contribute to the active exchangeable energy capable of promoting reactivity. To achieve the low sticking probabilities observed in experiment, severe normal translational energy dampening (ɛn → 0.26) was imposed, leading to a large vibrational efficacy of ηv = εv/εn = 3.85. The optimal threshold energy for dissociation was E0 = 30.88 kJ/mol, some 40 kJ/mol below the PBE-DFT prediction, and the optimal number of Au surface oscillators was s = 1. The d-PMMT modeling indicates that HCl/Au(111) reactivity can be consistent with electronically adiabatic passage across a relatively low and late transition state that dynamically disfavors normal translational energy.

A Linearized Semiclassical Dynamics Study of the Multiquantum Vibrational Relaxation of NO Scattering from a Au(111) Surface

The vibrational relaxation of NO molecules scattering from a Au(111) surface has served as the focus of efforts to understand nonadiabatic energy transfer at metal-molecule interfaces. Experimental measurements and previous theoretical efforts suggest that multiquantal NO vibrational energy relaxation occurs via electron-hole pair excitations in the metal. Here, using a linearized semiclassical approach, we accurately predict the vibrational relaxation of NO from the νi = 3 state for different incident translational energies. We also accurately capture the central role of transient electron transfer from the metal to the molecule in mediating the vibrational relaxation process but fall short of quantitatively predicting the full extent of multiquantum relaxation for high incident vibrational excitations (νi = 16).

Hydrogen atom scattering at the Al2O3(0001) surface: a combined experimental and theoretical study

Investigating atom-surface interactions is the key to an in-depth understanding of chemical processes at interfaces, which are of central importance in many fields - from heterogeneous catalysis to corrosion. In this work, we present a joint experimental and theoretical effort to gain insights into the atomistic details of hydrogen atom scattering at the α-Al2O3(0001) surface. Surprisingly, this system has been hardly studied to date, although hydrogen atoms as well as α-Al2O3 are omnipresent in catalysis as reactive species and support oxide, respectively. We address this system by performing hydrogen atom beam scattering experiments and molecular dynamics (MD) simulations based on a high-dimensional machine learning potential trained to density functional theory data. Using this combination of methods we are able to probe the properties of the multidimensional potential energy surface governing the scattering process. Specifically, we compare the angular distribution and the kinetic energy loss of the scattered atoms obtained in experiment with a large number of MD trajectories, which, moreover, allow to identify the underlying impact sites at the surface.

Best-of-Both-Worlds Predictive Approach to Dissociative Chemisorption on Metals

Predictive capability, accuracy, and affordability are essential features of a theory that is capable of describing dissociative chemisorption on a metal surface. This type of reaction is important for heterogeneous catalysis. Here we present an approach in which we use diffusion Monte Carlo (DMC) to pin the minimum barrier height and construct a density functional that reproduces this value. This predictive approach allows the construction of a potential energy surface at the cost of density functional theory while retaining near DMC accuracy. Scrutinizing effects of energy dissipation and quantum tunneling, dynamics calculations suggest the approach to be of near chemical accuracy, reproducing molecular beam sticking experiments for the showcase H2 + Al(110) system to ∼1.4 kcal/mol.

Nonadiabatic dynamics near metal surfaces under Floquet engineering: Floquet electronic friction vs Floquet surface hopping

In the previous study Wang and Dou [J. Chem. Phys. 158, 224109 (2023)], we have derived a Floquet classical master equation (FCME) to treat nonadiabatic dynamics near metal surfaces under Floquet engineering. We have also proposed a trajectory surface hopping algorithm to solve the FCME. In this study, we map the FCME into a Floquet Fokker-Planck equation in the limit of fast Floquet driving and fast electron motion as compared to nuclear motion. The Fokker-Planck equation is then being solved using Langevin dynamics with explicit friction and random force from the nonadiabatic effects of hybridized electrons and Floquet states. We benchmark the Floquet electronic friction dynamics against Floquet quantum master equation and Floquet surface hopping. We find that Floquet driving results in a violation of the second fluctuation-dissipation theorem, which further gives rise to heating effects.

Assessing Mixed Quantum-Classical Molecular Dynamics Methods for Nonadiabatic Dynamics of Molecules on Metal Surfaces

Mixed quantum-classical (MQC) methods for simulating the dynamics of molecules at metal surfaces have the potential to accurately and efficiently provide mechanistic insight into reactive processes. Here, we introduce simple two-dimensional models for the scattering of diatomic molecules at metal surfaces based on recently published electronic structure data. We apply several MQC methods to investigate their ability to capture how nonadiabatic effects influence molecule-metal energy transfer during the scattering process. Specifically, we compare molecular dynamics with electronic friction, Ehrenfest dynamics, independent electron surface hopping, and the broadened classical master equation approach. In the case of independent electron surface hopping, we implement a simple decoherence correction approach and assess its impact on vibrationally inelastic scattering. Our results show that simple, low-dimensional models can be used to qualitatively capture experimentally observed vibrational energy transfer and provide insight into the relative performance of different MQC schemes. We observe that all approaches predict similar kinetic energy dependence but return different vibrational energy distributions. Finally, by varying the molecule-metal coupling, we can assess the coupling regime in which some MQC methods become unsuitable.

First-principles surface reaction rates by ring polymer molecular dynamics and neural network potential: role of anharmonicity and lattice motion

Elementary gas-surface processes are essential steps in heterogeneous catalysis. A predictive understanding of catalytic mechanisms remains challenging due largely to difficulties in accurately characterizing the kinetics of such steps. Experimentally, thermal rates for elementary surface reactions can now be measured using a novel velocity imaging technique, providing a stringent testing ground for ab initio rate theories. Here, we propose to combine ring polymer molecular dynamics (RPMD) rate theory with state-of-the-art first-principles-determined neural network potential to calculate surface reaction rates. Taking NO desorption from Pd(111) as an example, we show that the harmonic approximation and the neglect of lattice motion in the commonly-used transition state theory overestimates and underestimates the entropy change during the desorption process, respectively, leading to opposite errors in rate coefficient predictions and artificial error cancellations. Including anharmonicity and lattice motion, our results reveal a generally neglected surface entropy change due to significant local structural change during desorption and obtain the right answer for the right reasons. Although quantum effects are found to be less important in this system, the proposed approach establishes a more reliable theoretical benchmark for accurately predicting the kinetics of elementary gas-surface processes.

Efficient implementation and performance analysis of the independent electron surface hopping method for dynamics at metal surfaces

Independent electron surface hopping (IESH) is a computational algorithm for simulating the mixed quantum-classical molecular dynamics of adsorbate atoms and molecules interacting with metal surfaces. It is capable of modeling the nonadiabatic effects of electron-hole pair excitations on molecular dynamics. Here, we present a transparent, reliable, and efficient implementation of IESH, demonstrating its ability to predict scattering and desorption probabilities across a variety of systems, ranging from model Hamiltonians to full dimensional atomistic systems. We further show how the algorithm can be modified to account for the application of an external bias potential, comparing its accuracy to results obtained using the hierarchical quantum master equation. Our results show that IESH is a practical method for modeling coupled electron-nuclear dynamics at metal surfaces, especially for highly energetic scattering events.

Stereodynamics of adiabatic and non-adiabatic energy transfer in a molecule surface encounter

Molecular energy transfer and reactions at solid surfaces depend on the molecular orientation relative to the surface. While such steric effects have been largely understood in electronically adiabatic processes, the orientation-dependent energy transfer in NO scattering from Au(111) was complicated by electron-mediated nonadiabatic effects, thus lacking a clear interpretation and posing a great challenge for theories. Herein, we investigate the stereodynamics of adiabatic and nonadiabatic energy transfer via molecular dynamics simulations of NO(v = 3) scattering from Au(111) using realistic initial orientation distributions based on accurate neural network fitted adiabatic potential energy surface and electronic friction tensor. Our results reproduce the observed stronger vibrational relaxation for N-first orientation and enhanced rotational rainbow for O-first orientation, and demonstrate how adiabatic anisotropic interactions steer molecules into the more attractive N-first orientation to experience more significant energy transfer. Remaining disagreements with experiment suggest the direction for further developments of nonadiabatic theories for gas-surface scattering.

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.

Quantum Dynamics with Electronic Friction

A theory of electronic friction is developed using the exact factorization of the electronic-nuclear wave function. No assumption is made regarding the electronic bath, which can be made of independent or interacting electrons, and the nuclei are treated quantally. The ensuing equation of motion for the nuclear wave function is a nonlinear Schrödinger equation including a friction term. The resulting friction kernel agrees with a previously derived mixed quantum-classical result by Dou et al., [Phys. Rev. Lett. 119, 046001 (2017)]PRLTAO0031-900710.1103/PhysRevLett.119.046001, except for a pseudomagnetic contribution in the latter that is here removed. More specifically, it is shown that the electron dynamics generally washes out the gauge fields appearing in the adiabatic dynamics. However, these are fully re-established in the typical situation where the electrons respond rapidly on the slow time scale of the nuclear dynamics (Markov limit). Hence, we predict Berry's phase effects to be observable also in the presence of electronic friction. Application to a model vibrational relaxation problem proves that the proposed approach represents a viable way to account for electronic friction in a fully quantum setting for the nuclear dynamics.

NQCDynamics.jl: A Julia package for nonadiabatic quantum classical molecular dynamics in the condensed phase

Accurate and efficient methods to simulate nonadiabatic and quantum nuclear effects in high-dimensional and dissipative systems are crucial for the prediction of chemical dynamics in the condensed phase. To facilitate effective development, code sharing, and uptake of newly developed dynamics methods, it is important that software implementations can be easily accessed and built upon. Using the Julia programming language, we have developed the NQCDynamics.jl package, which provides a framework for established and emerging methods for performing semiclassical and mixed quantum-classical dynamics in the condensed phase. The code provides several interfaces to existing atomistic simulation frameworks, electronic structure codes, and machine learning representations. In addition to the existing methods, the package provides infrastructure for developing and deploying new dynamics methods, which we hope will benefit reproducibility and code sharing in the field of condensed phase quantum dynamics. Herein, we present our code design choices and the specific Julia programming features from which they benefit. We further demonstrate the capabilities of the package on two examples of chemical dynamics in the condensed phase: the population dynamics of the spin-boson model as described by a wide variety of semiclassical and mixed quantum-classical nonadiabatic methods and the reactive scattering of H₂ on Ag(111) using the molecular dynamics with electronic friction method. Together, they exemplify the broad scope of the package to study effective model Hamiltonians and realistic atomistic systems.

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.

REANN: A PyTorch-based end-to-end multi-functional deep neural network package for molecular, reactive, and periodic systems

In this work, we present a general purpose deep neural network package for representing energies, forces, dipole moments, and polarizabilities of atomistic systems. This so-called recursively embedded atom neural network model takes advantages of both the physically inspired atomic descriptor based neural networks and the message-passing based neural networks. Implemented in the PyTorch framework, the training process is parallelized on both the central processing unit and the graphics processing unit with high efficiency and low memory in which all hyperparameters can be optimized automatically. We demonstrate the state-of-the-art accuracy, high efficiency, scalability, and universality of this package by learning not only energies (with or without forces) but also dipole moment vectors and polarizability tensors in various molecular, reactive, and periodic systems. An interface between a trained model and LAMMPs is provided for large scale molecular dynamics simulations. We hope that this open-source toolbox will allow for future method development and applications of machine learned potential energy surfaces and quantum-chemical properties of molecules, reactions, and materials.

Plasmonic enhancement of molecular hydrogen dissociation on metallic magnesium nanoclusters

Light-driven plasmonic enhancement of chemical reactions on metal catalysts is a promising strategy to achieve highly selective and efficient chemical transformations. The study of plasmonic catalyst materials has traditionally focused on late transition metals such as Au, Ag, and Cu. In recent years, there has been increasing interest in the plasmonic properties of a set of earth-abundant elements such as Mg, which exhibit interesting hydrogenation chemistry with potential applications in hydrogen storage. This work explores the optical, electronic, and catalytic properties of a set of metallic Mg nanoclusters with up to 2057 atoms using time-dependent density functional tight-binding and density functional theory calculations. Our results show that Mg nanoclusters are able to produce highly energetic hot electrons with energies of up to 4 eV. By electronic structure analysis, we find that these hot electrons energetically align with electronic states of physisorbed molecular hydrogen, occupation of which by hot electrons can promote the hydrogen dissociation reaction. We also find that the reverse reaction, hydrogen evolution on metallic Mg, can potentially be promoted by hot electrons, but following a different mechanism. Thus, from a theoretical perspective, Mg nanoclusters display very promising behaviour for their use in light promoted storage and release of hydrogen.

Perspective on integrating machine learning into computational chemistry and materials science

Machine learning (ML) methods are being used in almost every conceivable area of electronic structure theory and molecular simulation. In particular, ML has become firmly established in the construction of high-dimensional interatomic potentials. Not a day goes by without another proof of principle being published on how ML methods can represent and predict quantum mechanical properties-be they observable, such as molecular polarizabilities, or not, such as atomic charges. As ML is becoming pervasive in electronic structure theory and molecular simulation, we provide an overview of how atomistic computational modeling is being transformed by the incorporation of ML approaches. From the perspective of the practitioner in the field, we assess how common workflows to predict structure, dynamics, and spectroscopy are affected by ML. Finally, we discuss how a tighter and lasting integration of ML methods with computational chemistry and materials science can be achieved and what it will mean for research practice, software development, and postgraduate training.

Searching Configurations in Uncertainty Space: Active Learning of High-Dimensional Neural Network Reactive Potentials

Neural network (NN) potential energy surfaces (PESs) have been widely used in atomistic simulations with ab initio accuracy. While constructing NN PESs, their training data points are often sampled by molecular dynamics trajectories. This strategy can be however inefficient for reactive systems involving rare events. Here, we develop an uncertainty-driven active learning strategy to automatically and efficiently generate high-dimensional NN-based reactive potentials, taking a gas-surface reaction as an example. The difference between two independent NN models is used as a simple and differentiable uncertainty metric, allowing us to quickly search in the uncertainty space and place new samples at which the PES is less reliable. By interfacing this algorithm with the first-principles simulation package, we demonstrate that a globally accurate NN potential of the H₂ + Ag(111) system can be constructed with merely ∼150 data points. This PES can be further refined to describe H₂ dissociation on Ag(100) by adding ∼130 more configurations on this facet. The entire process is completely automatic and self-terminated once the relative error criterion is fulfilled. Impressively, data points sampled by this uncertainty-driven strategy are substantially fewer than by the traditional trajectory-based sampling. The final NN PES not only converges well the quantum dissociation probability of the molecule but also well-reproduces the phonon properties of the substrate and is capable of describing surface temperature effects. These results show the potential of this active learning approach in developing high-dimensional NN reactive potentials in gas and condensed phases.

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.