Energetics and dynamics of CH4 and H2O dissociation on metal surfaces

Lowering the C–H bond activation barrier of methane by means of SAC@Cu(111): periodic DFT investigations

Methane has long been in the world's spotlight as the simplest yet one of the most notorious hydrocarbons; here, we study the efficiency of single-atom catalysts (SACs) for methane activation using density functional theory (DFT).

Theoretical study on water behavior on the copper surfaces

We calculated the adsorption of H, O, OH, and H₂O and the dissociation of H₂O molecule on the Cu(111), Cu(100), and Cu(110) surfaces using density functional theory. H, O, and OH tend to adsorb stably at the highly coordinated dh and h sites on the Cu(111) and Cu(100) surfaces. OH and H tend to adsorb on sb site on the Cu(110) surface. The more charge transfer of the adsorbed substance, the more stable the adsorption. The dissociation product is O+H on the Cu(111) surface, while the dissociation product is OH+H on the Cu(100) and Cu(110) surfaces. Due to the different geometric structures of initial state (IS), transition state (TS), and final state (FS) in the dissociation reaction, the dissociation of water on the copper surface does not establish a linear Brønsted-Evans-Polanyi (BEP) relationship. These results provide theoretical support for the understanding of the interaction between water and metals as well as the behavior of water molecules.

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.

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.