Water dissociating on rigid Ni(100): A quantum dynamics study on a full-dimensional potential energy surface

Six-dimensional quantum dynamics study for the dissociative chemisorption of H2 on pure and alloyed AgAu surfaces

The 6D time-dependent wave packet calculations were performed to explore H2 dissociation on Ag, Au, and two AgAu alloy surfaces, using four newly fitted potential energy surfaces based on the neural network fitting to density functional theory energy points. The ligand effect resulting from the Ag-Au interaction causes a reduction in the barrier height for H2+Ag/Au(111) compared to H2+Ag(111). However, the scenario is reversed for H2+Au/Ag(111) and H2+Au(111). The 6D dissociation probabilities of H2 on Ag/Au(111) surfaces are significantly higher than those on the pure Ag(111) surface, but the corresponding results for H2 on Au/Ag(111) surfaces are substantially lower than those on the pure Au(111) surface. The reactivity of H2 on Au(111) is larger than that on Ag(111), despite Ag(111) having a slightly lower static barrier height. This can be attributed to the exceptionally small dissociation probabilities at the hcp and fcc regions, which are at least 100 times smaller compared to those at the bridge or top site for H2+Ag(111). Due to the late barrier being more pronounced, the vibrational excitation of H2 on Ag(111) is more effective in promoting the reaction than on Au(111). Moreover, a high degree of alignment dependence is detected for the four reactions, where the H2 dissociation has the highest probability at the helicopter alignment, as opposed to the cartwheel alignment.

Charge-Transfer-Controlled Quantum Dynamics of HCl Dissociation on the Ag/Au(111) Bimetallic Alloy Surface

Understanding polar molecule dynamics on bimetallic surfaces, especially electropositivity and electronegativity, remains a challenge. Here, we report the reactivity of HCl on a strained Ag monolayer on Au(111) using six-dimensional quantum dynamics with a new machine-learning-based potential energy surface. Surprisingly, HCl reactivity is significantly suppressed by the Ag-Au interaction despite a lower HCl+Ag/Au(111) barrier than pure Ag(111). This arises from charge transfer between Ag and Au, where electronegative Au makes the top Ag layer on Ag/Au(111) electropositive, unlike that on pure Ag(111). Electropositive Ag in HCl+Ag/Au(111) attracts Cl, yielding an unfavorable H-Cl configuration and reduced reactivity. These findings deepen our understanding of polar molecule interactions on bimetallic surfaces, highlighting the role of charge transfer in dissociative chemisorption and the implications for catalyst design in heterogeneous catalysis.

Quantum dynamics reveal different ligand effects by vibrational excitation in the dissociative chemisorption of HCl on the Au/Ag(111) surface

The reactivity and selectivity of bimetallic surfaces are of fundamental importance in industrial applications. Here, we report the first six-dimensional (6D) quantum dynamics study for the role of surface strain and ligand effects on the reactivity of HCl on a strained pseudomorphic monolayer of Au deposited onto a Ag(111) substrate, with the aid of accurate machine learning-based potential energy surfaces. The substitute of Au into Ag changes the location of the transition state; however, the static barrier height remains roughly the same as pure Au(111). The 6D quantum dynamics calculations reveal that the surface strain due to lattice expansion slightly enhances the reactivity. The ligand effect due to electronic structure interactions between Au and Ag substantially suppresses the reactivity of HCl in the ground vibrational state but promotes the reactivity via vibrational excitation at high kinetic energies. This finding can be attributed to more close interaction with Ag atoms at the transition state close to the fcc site, as well as the tight transition-state region, making the vibrational excitation highly efficient in enhancing the reactivity. Our study quantitatively unravels the dynamical origin of reactivity control by two metals, which will ultimately provide valuable insight into the selectivity of the catalyst.

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.

Theoretical study on water adsorption and dissociation on the nickel surfaces

Using density functional theory methods, H₂O dissociation was investigated on the Ni(111), Ni(100), and Ni(110) surfaces. H and O atom as well as OH species adsorb stably at the high coordination sites. While on the Ni(110) surface, the OH species prefers at the twofold short bridge site because the adsorption on the fourfold hollow site is less feasible due to the increased distances between the nickel atoms. The amount of charge transfer is related to the adsorption stability. The more charge transfer, the more stable the adsorption. The charge transfer decreases in the order of O > OH > H. H₂O molecule adsorbs at the top site in a configuration parallel to the surface. The final products are different for H₂O dissociation due to the different mechanisms. On the Ni(111) surface, the final product is the O atom. On the Ni(100) and Ni(110) surfaces, the most abundant species are OH and H, but the reaction mechanisms were different. It is not necessary to linear BEP relationship for a given reaction on different surfaces. These results could provide fundamental insights into water behaviors and a favorable theoretical basis for further understanding and research on the interaction between water and metal surfaces.

Physisorbed State Regulates the Dissociation Mechanism of H2O on Ni(100)

Water dissociation is a key step in many industrial catalytic processes. The dissociation of H₂O on a rigid Ni(100) surface was investigated by the quantum instanton method with a full-dimensional potential energy surface. The calculated free-energy barrier maps showed that the free-energy barrier varied dramatically with the surface site. The free-energy well map demonstrated that the physisorption well of H₂O was existent at all of the surface sites, and H₂O could be dissociated by both the direct and steady-state processes. The calculated direct dissociation rate constants at different surface sites decreased rapidly in the order transition state (TS) > bridge > top > hollow. The steady-state dissociation rate constants had the same trend as that of the direct process but the steady-state dissociation rate constant at the top site became the largest at high temperatures. The direct dissociation rate constants were always larger than those of the steady-state process at a given temperature. The calculated kinetic isotope effects for the direct and steady-state processes were extremely large at low temperatures, which was caused by the zero-point energy correction and remarkable quantum tunneling. From low temperature to high temperature, H₂O would undergo stable molecular adsorption at the top site, steady-state dissociation at the TS site, direct rupture at the TS site, and direct decomposition at the impact site.

High-Fidelity Potential Energy Surfaces for Gas-Phase and Gas-Surface Scattering Processes from Machine Learning

In this Perspective, we review recent advances in constructing high-fidelity potential energy surfaces (PESs) from discrete ab initio points, using machine learning tools. Such PESs, albeit with substantial initial investments, provide significantly higher efficiency than direct dynamics methods and/or high accuracy at a level that is not affordable by on-the-fly approaches. These PESs not only are a necessity for quantum dynamical studies because of delocalization of wave packets but also enable the study of low-probability and long-time events in (quasi-)classical treatments. Our focus here is on inelastic and reactive scattering processes, which are more challenging than bound systems because of the involvement of continua. Relevant applications and developments for dynamical processes in both the gas phase and at gas-surface interfaces are discussed.

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.