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.

Dissociative adsorption dynamics of nitrogen on a Fe(111) surface

The dissociative adsorption dynamics of N₂ on clean Fe(111) surfaces is theoretically investigated by means of quasi-classical trajectory calculations based on a multidimensional potential energy surface built from density functional theory.

Hydrogen diffusion into the subsurfaces of model metal catalysts from first principles

We present in detail the diffusion pathways of hydrogen on seven different open and closed model catalyst surfaces from first principles calculations.

Adsorption dynamics of molecular nitrogen at an Fe(111) surface

We present an extensive theoretical study of N₂ adsorption mechanisms on an Fe(111) surface.

Methane adsorption and dissociation on iron oxide oxygen carriers: the role of oxygen vacancies

We investigated the interaction between methane and iron oxide oxygen carrier for chemical looping reaction systems and found the oxygen vacancies can facilitate methane conversion.

Molecular dynamics simulation of O2 adsorption on Ag(110) from first principles electronic structure calculations

State of the art simulations show that the physisorption state could be important for O₂/Ag(110) adsorption.

Oxygen vacancy promoted methane partial oxidation over iron oxide oxygen carriers in the chemical looping process

We found that oxygen vacancies can promote CH₄ partial oxidation on iron oxide oxygen carriers during the chemical looping process.

Dissociative dynamics of O2 on Ag(110)

The dissociative dynamics of O₂ on Ag(110) is studied using classical (CT) and quasiclassical (QCT) trajectory calculations performed on an accurate multidimensional DFT potential energy surface.

The dynamics of adsorption and dissociation of N2 in a monolayer of iron on W(110)

Does N₂ adsorption increase on strained Fe monolayers?

Surface strain improves molecular adsorption but hampers dissociation for N2 on the Fe/W(110) surface

We compare the adsorption dynamics of N(2) on the unstrained Fe(110) and on a 10% expanded Fe monolayer grown on W(110) by performing classical molecular dynamics simulations that use potential energy surfaces calculated with density functional theory. Our results allow us to understand why, experimentally, the molecular adsorption of N(2) is observed on the strained layer but not on Fe(110). Surprisingly, we also find that while surface strain favors the molecular adsorption of N(2) it seems, on the contrary, to impede the dissociative adsorption. This result contrasts with previous examples for which strain is found to modify equally the energetics of chemisorption and dissociation.

Role of physisorption states in molecular scattering: a semilocal density-functional theory study on O2/Ag(111)

We simulate the scattering of O2 from Ag(111) with classical dynamics simulations performed on a six-dimensional potential energy surface calculated within semilocal density-functional theory. The enigmatic experimental trends that originally required the conjecture of two types of repulsive walls, arising from a physisorption and chemisorption part of the interaction potential, are fully reproduced. Given the inadequate description of the physisorption properties in semilocal density-functional theory, our work casts severe doubts on the prevalent notion to use molecular scattering data as indirect evidence for the existence of such states.