Controlling Heterogeneous Catalysis of Water Dissociation Using Cu–Ni Bimetallic Alloy Surfaces: A Quantum Dynamics Study

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.

Mechanistic Insights for Dry Reforming of Methane on Cu/Ni Bimetallic Catalysts: DFT-Assisted Microkinetic Analysis for Coke Resistance

Density functional theory (DFT) calculations have been utilized to evaluate the complete reaction mechanism of methane dry reforming (DRM) over Ni2Cu (111) bimetallic catalyst. The detailed catalytic cycle on Ni2Cu (111) catalyst demonstrated superior coke resistance compared to pure Ni (111) and Ni2Fe (111) reported in the literature. Doping Cu in the Ni–Ni network enhanced the competitive CH oxidation by both atomic O and OH species with the latter having only 0.02 eV higher than the 1.06 eV energy barrier required for CH oxidation by atomic O. Among the C/CH oxidation pathways, C* + O* → CO (g) was the most favorable with an energy barrier of 0.72 eV. This was almost half of the energy barrier required for the rate-limiting step of CH decomposition (1.40 eV) and indicated enhanced coke deposition removal. Finally, we investigated the effect of temperature (800~1000 K) on the carbon deposition and elimination mechanism over Ni2Cu (111) catalyst. Under those realistic DRM conditions, the calculations showed a periodic cycle of simultaneous carbon deposition and elimination resulting in improved catalyst stability.

Polyhedron-Assembled Ternary PtCuCo Nanochains: Integrated Functions Enhance the Electrocatalytic Performance of Methanol Oxidation at Elevated Temperature

Recently, the preparation of a high-performance one-dimensional alloy nanostructure for fuel cells has been given increasing attention due to its smart-structure merits and electronic effect triggered by alloying different kinds of metals at the nanoscale. In this study, unique ternary PtCuCo nanochains assembled with small polyhedra are first achieved and used as high-performance anode electrocatalysts toward methanol oxidation at elevated temperature (60 °C) that is closer to the operating temperature of direct methanol fuel cells than room temperature. The specific activity/mass activity of Pt₄₅Cu₃₅Co₂₀ one-dimensional nanochains can reach up to 18.24 mA cm^-2/4.19 A mg^-1_Pt that is 9.25/10.47 times that of commercial Pt black in sulfuric acid medium. After a 3600 s durability test, the remaining current density of Pt₄₅Cu₃₅Co₂₀ one-dimensional nanochains is 73.3 times that of commercial Pt black. The structure characterizations show that the high density of surface active sites, d-band center of the Pt downshift, moderate strain effect, and synergetic effect are jointly responsible for the enhanced electrocatalytic performance of one-dimensional ternary PtCuCo nanochains.

Coverage dependent structure and energy of water dissociative adsorption on clean and O-pre-covered Ni (100) and Ni(110)

H₂O dissociative adsorption on clean and O pre-covered Ni(100) and Ni(110) surfaces has been computed systematically on the basis of periodic density functional theory and ab initio atomistic thermodynamics.