Hot Adatom Diffusion Following Oxygen Dissociation on Pd(100) and Pd(111): A First-Principles Study of the Equilibration Dynamics of Exothermic Surface Reactions

On the quantum dynamical treatment of surface vibrational modes for reactive scattering of H2 from Cu(111) at 925 K

We construct the effective Hartree potential for H2 on Cu(111) as introduced in our earlier work [Dutta et al., J. Chem. Phys. 154, 104103 (2021), and Dutta et al., J. Chem. Phys. 157, 194112 (2022)] starting from the same gas-metal interaction potential obtained for 0 K. Unlike in that work, we now explicitly account for surface expansion at 925 K and investigate different models to describe the surface vibrational modes: (i) a cluster model yielding harmonic normal modes at 0 K and (ii) slab models resulting in phonons at 0 and 925 K according to the quasi-harmonic approximation-all consistently calculated at the density functional theory level with the same exchange-correlation potential. While performing dynamical calculations for the H2(v = 0, j = 0)-Cu(111) system employing Hartree potential constructed with 925 K phonons and surface temperature, (i) the calculated chemisorption probabilities are the highest compared to the other approaches over the energy domain and (ii) the threshold for the reaction probability is the lowest, in close agreement with the experiment. Although the survival probabilities (v' = 0) depict the expected trend (lower in magnitude), the excitation probabilities (v' = 1) display a higher magnitude since the 925 K phonons and surface temperature are more effective for the excitation process compared to the phonons/normal modes obtained from the other approaches investigated to describe the surface.

Advances in the Dynamics of Adsorbate Diffusion on Metal Surfaces: Focus on Hydrogen and Oxygen

Adsorbates on metal surfaces are typically formed from the dissociative chemisorption of molecules occurring at gas-solid interfaces. These adsorbed species exhibit unique diffusion behaviors on metal surfaces, which are influenced by their translational energy. They play crucial roles in various fields, including heterogeneous catalysis and corrosion. This review examines recent theoretical advancements in understanding the diffusion dynamics of adsorbates on metal surfaces, with a specific emphasis on hydrogen and oxygen atoms. The diffusion processes of adsorbates on metal surfaces involve two energy transfer mechanisms: surface phonons and electron-hole pair excitations. This review also surveys new theoretical methods, including the characterization of the electron-hole pair excitation within electronic friction models, the acceleration of quantum chemistry calculations through machine learning, and the treatment of atomic nuclear motion from both quantum mechanical and classical perspectives. Furthermore, this review offers valuable insights into how energy transfer, nuclear quantum effects, supercell sizes, and the topography of potential energy surfaces impact the diffusion behavior of hydrogen and oxygen species on metal surfaces. Lastly, some preliminary research proposals are presented.

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.

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.

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.