CO Blocking of D2Dissociative Adsorption on Ru(0001)

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 characterization of the surface composition of ruthenium nanoparticles in equilibrium with syngas

A deeper understanding of the relationship between experimental reaction conditions and the surface composition of nanoparticles is crucial in order to elucidate mechanisms involved in nanocatalysis. In the framework of the Fischer-Tropsch synthesis, a resolution of this complex puzzle requires a detailed understanding of the interaction of CO and H with the surface of the catalyst. In this context, the single- and co-adsorption of CO and H to the surface of a 1 nm ruthenium nanoparticle has been investigated with density functional theory. Using several indexes (d-band center, crystal overlap Hamilton population, density of states), a systematic analysis of the bond properties and of the electronic states has also been done, in order to bring an understanding of structure/property relationships at the nanoscale. The H : CO surface composition of this ruthenium nanoparticle exposed to syngas has been evaluated according to a thermodynamic model fed with DFT energies. Such ab initio thermodynamic calculations give access to the optimal H : CO coverage values under a wide range of experimental conditions, through the construction of free energy phase diagrams. Surprisingly, under the Fischer-Tropsch synthesis experimental conditions, and in agreement with new experiments, only CO species are adsorbed at the surface of the nanoparticle. These findings shed new light on the possible reaction pathways underlying the Fischer-Tropsch synthesis, and specifically the initiation of the reaction. It is finally shown that the joint knowledge of the surface composition and energy descriptors can help to identify possible reaction intermediates.

Quantum and classical dynamics of reactive scattering of H2 from metal surfaces

State-of-the-art theoretical models allow nowadays an accurate description of H₂/metal surface systems and phenomena relative to heterogeneous catalysis. Here we review the most relevant ones investigated during the last 10 years.

An ab initio molecular dynamics study of D2 dissociation on CO-precovered Ru(0001)

In dynamics studies of hydrogen dissociation on CO pre-covered Ru(0001) the simulation cell size is important for describing energy exchange.

High coverage adsorption and co-adsorption of CO and H2 on Ru(0001) from DFT and thermodynamics

The adsorption and co-adsorption of CO and H₂ at different coverages on p(4 × 4) Ru(0001) have been computed using periodic density functional theory (GGA-RPBE) and atomistic thermodynamics.

Adsorption features of flavonoids on macroporous adsorption resins functionalized with ionic liquids

A series of macroporous adsorption resins (MARs) with novel structures is synthesized via Friedel-Crafts catalyzed reaction. The adsorption kinetics of the synthetic resins with respect to the purification effect is systematically investigated by means of the response surface methodology (RSM). The kinetic data cannot be fitted to the classical model because it does not take multicompartments and desorption rates into consideration. A new multicompartment louver-tide theory is thus developed considering that adsorption is an indefinite dynamic equilibrium process, which can be divided into innumerable ingredients with different desorption rates. This theory produces much better fits to the experimental data and provides a quantitative explanation with multicompartments and adsorption/desorption rates.

Towards chemically accurate simulation of molecule-surface reactions

This perspective addresses four challenges facing theorists whose aim is to make quantitatively accurate predictions for reactions of molecules on metal surfaces, and suggests ways of meeting these challenges, focusing on dissociative chemisorption reactions of H(2), N(2), and CH(4). Addressing these challenges is ultimately of practical importance to a more accurate description of overall heterogeneously catalysed reactions, which play a role in the production of more than 90% of man-made chemicals. One challenge is to describe the interaction of a molecule with a metal surface with chemical accuracy, i.e., with errors in reaction barrier heights less than 1 kcal mol(-1). In this framework, the potential of a new implementation of specific reaction parameter density functional theory (SRP-DFT) will be discussed, with emphasis on applications to reaction of H(2) with metal surfaces. Two additional challenges are to come up with improved descriptions of the effects of phonons and electron-hole pairs on reaction of molecules like N(2) on metal surfaces. Phonons can be tackled using sudden approximations in quantum dynamics, and through Ab Initio Molecular Dynamics (AIMD) calculations using classical dynamics. To additionally achieve an accurate description of the effect of electron-hole pair excitation on dissociative chemisorption within a classical dynamics framework, it may be possible to combine AIMD with electronic friction. The fourth challenge we will consider is how to achieve an accurate quantum mechanical description of the dissociative chemisorption of a polyatomic molecule, like methane, on a metal surface. A method of potential interest is the Multi-Configuration Time-Dependent Hartree (MCTDH) method.

Evidence of stable high-temperature D(x)-CO intermediates on the Ru(0001) surface

We demonstrate the formation of complexes involving attractive interactions between D and CO on Ru(0001) that are stable at significantly higher temperatures than have previously been reported for such intermediate species on this surface. These complexes are evident by the appearance of new desorption features upon heating of the sample. They decompose in stages as the sample temperature is increased, with the most stable component desorbing at >500 K. The D:CO ratio remaining on the surface during the final stages of desorption tends towards 1:1. The new features are populated during normally incident molecular beam dosing of D(2) on to CO pre-covered Ru(0001) surfaces (180 K) when the CO coverage exceeds 50% of the saturation value. The amount of complex formed decreases somewhat with increasing CO pre-coverage. It is almost absent in the case of dosing on to the fully saturated surface. The results are interpreted in terms of both local and long-range rearrangements of the overlayer that give rise to the observed CO coverage dependence and limit the amount of complex that can be formed.