CO Oxidation Reaction over Oxygen-Rich Ru(0001) Surfaces

Crystal-Phase Engineering in Heterogeneous Catalysis

The performance of a chemical reaction is critically dependent on the electronic and/or geometric structures of a material in heterogeneous catalysis. Over the past century, the Sabatier principle has already provided a conceptual framework for optimal catalyst design by adjusting the electronic structure of the catalytic material via a change in composition. Beyond composition, it is essential to recognize that the geometric atomic structures of a catalyst, encompassing terraces, edges, steps, kinks, and corners, have a substantial impact on the activity and selectivity of a chemical reaction. Crystal-phase engineering has the capacity to bring about substantial alterations in the electronic and geometric configurations of a catalyst, enabling control over coordination numbers, morphological features, and the arrangement of surface atoms. Modulating the crystallographic phase is therefore an important strategy for improving the stability, activity, and selectivity of catalytic materials. Nonetheless, a complete understanding of how the performance depends on the crystal phase of a catalyst remains elusive, primarily due to the absence of a molecular-level view of active sites across various crystal phases. In this review, we primarily focus on assessing the dependence of catalytic performance on crystal phases to elucidate the challenges and complexities inherent in heterogeneous catalysis, ultimately aiming for improved catalyst design.

Theoretical design of a technetium-like alloy and its catalytic properties

Based on the concept of density of states (DOS) engineering, we theoretically designed a pseudo-Tc material (Mo-Ru alloy) and investigated its electronic structure, phase stability and catalytic activity by using density functional theory. Through comparing the DOS shape, peak distribution, and DOS area differences between Tc and the Mo-Ru alloy, we noticed that bcc-Mo8Ru8 and hcp-Mo8Ru8 had the most similar electronic structures to Tc. The excess energies after entropy correction of hcp-Mo8Ru8 and bcc-Mo8Ru8 are stable when the temperature is up to 765 and 745 K, respectively. These results provided the possibility of pseudo-Tc alloy (hcp-Mo8Ru8 and bcc-Mo8Ru8) synthesis. Finally, according to reaction coordinate analysis, the similar catalytic activity between hcp-Mo8Ru8 and Tc have been demonstrated in CO oxidation and N2 dissociation. In N2 dissociation, Tc has a suitable ratio of transition state (TS) barrier to reaction energy which make Tc an efficient catalyst for NH3 synthesis, in addition to our designed pseudo-Tc (hcp-MoRu) because of the similar electronic structures. Our finding provides valuable insight into materials and catalyst design.

Evidence of the O–Pd–O and Pd–O4 structure units as oxide seeds and their origin on Pd(211): revealing the mechanism of surface oxide formation

The formation of surface oxides on metal surfaces is not only important in materials science, but also of significance in heterogeneous catalysis due to the fact that during most oxidation reactions the metal catalysts are inevitably oxidized, which may cause dramatic consequences in the reactions.

Strong Anisotropic Interaction Controls Unusual Sticking and Scattering of CO at Ru(0001)

Complete sticking at low incidence energies and broad angular scattering distributions at higher energies are often observed in molecular beam experiments on gas-surface systems which feature a deep chemisorption well and lack early reaction barriers. Although CO binds strongly on Ru(0001), scattering is characterized by rather narrow angular distributions and sticking is incomplete even at low incidence energies. We perform molecular dynamics simulations, accounting for phononic (and electronic) energy loss channels, on a potential energy surface based on first-principles electronic structure calculations that reproduce the molecular beam experiments. We demonstrate that the mentioned unusual behavior is a consequence of a very strong rotational anisotropy in the molecule-surface interaction potential. Beyond the interpretation of scattering phenomena, we also discuss implications of our results for the recently proposed role of a precursor state for the desorption and scattering of CO from ruthenium.

Ferroelectric triggering of carbon monoxide adsorption on lead zirco-titanate (001) surfaces

Atomically clean lead zirco-titanate PbZr_0.2Ti_0.8O₃ (001) layers exhibit a polarization oriented inwards P⁽⁻⁾, visible by a band bending of all core levels towards lower binding energies, whereas as introduced layers exhibit P⁽⁺⁾ polarization under air or in ultrahigh vacuum. The magnitude of the inwards polarization decreases when the temperature is increased at 700 K. CO adsorption on P⁽⁻⁾ polarized surfaces saturates at about one quarter of a monolayer of carbon, and occurs in both molecular (oxidized) and dissociated (reduced) states of carbon, with a large majority of reduced state. The sticking of CO on the surface in ultrahigh vacuum is found to be directly related to the P⁽⁻⁾ polarization state of the surface. A simple electrostatic mechanism is proposed to explain these dissociation processes and the sticking of carbon on P⁽⁻⁾ polarized areas. Carbon desorbs also when the surface is irradiated with soft X-rays. Carbon desorption when the polarization is lost proceeds most probably in form of CO₂. Upon carbon desorption cycles, the ferroelectric surface is depleted in oxygen and at some point reverses its polarization, owing to electrons provided by oxygen vacancies which are able to screen the depolarization field produced by positive fixed charges at the surface.

Indication of non-thermal contribution to visible femtosecond laser-induced CO oxidation on Ru(0001)

We studied CO oxidation on Ru(0001) induced by 400 nm and 800 nm femtosecond laser pulses where we find a branching ratio between CO oxidation and desorption of 1:9 and 1:31, respectively, showing higher selectivity towards CO oxidation for the shorter wavelength excitation. Activation energies computed with density functional theory show discrepancies with values extracted from the experiments, indicating both a mixture between different adsorbed phases and importance of non-adiabatic effects on the effective barrier for oxidation. We simulated the reactions using kinetic modeling based on the two-temperature model of laser-induced energy transfer in the substrate combined with a friction model for the coupling to adsorbate vibrations. This model gives an overall good agreement with experiment except for the substantial difference in yield ratio between CO oxidation and desorption at 400 nm and 800 nm. However, including also the initial, non-thermal effect of electrons transiently excited into antibonding states of the O-Ru bond yielded good agreement with all experimental results.

Oscillatory Behavior in the CO-oxidation over Bulk Ruthenium Dioxide — the Effect of the CO/O2 Ratio

CO oxidation over polycrystalline ruthenium dioxide was monitored in an in-situ XRD setup. The evolution of the bulk state of the catalyst was followed by in-situ XRD during reaction, while the surface morphology and chemical state before and after reaction were investigated by HRSEM and EDX. The commercial RuO2 powder was calcined prior reaction to ensure the formation of completely oxidized RuO2. This pre-calcined RuO2 is initially inactive in CO oxidation regardless of the CO/O2 feed ratio and requires an induction period, the length of which strongly depends whether the catalyst is diluted with boron nitride or not. After this induction period oscillations in the CO2 yield occur under O2-rich conditions only. These oscillations exhibit two time constants for the diluted catalyst, while the low frequency oscillations were not observed in the case of undiluted RuO2. Furthermore, the state of the catalyst after activation in O2-rich feed conditions differs dramatically from the state after activation in CO-rich feed conditions. Firstly, the catalyst activated in an O2-rich atmosphere remains inactive under CO-rich conditions in contrast to the catalyst activated in CO-rich conditions which is afterwards active under all feed ratios examined. Secondly, the surface morphology of the catalyst is quite different. While the apical surfaces of the RuO2 crystals become roughened upon activation in the CO-rich feed, they become facetted under O2 rich activation conditions. Therefore, we conclude that at least two different active surface states on the bulk RuO2 catalyst exist.

On the CO-Oxidation over Oxygenated Ruthenium

The oxidation of carbon monoxide over polycrystalline ruthenium dioxide (RuO2) powder was studied in a packed-bed reactor and by bulk and surface analytical methods. Activity data were correlated with bulk phases in an in-situ X-ray diffraction (XRD) setup at atmospheric pressure. Ruthenium dioxide was pre-calcined in pure oxygen at 1073 K. At this stage RuO2 is completely inactive in the oxidation of CO. After a long induction period in the feed at 503 K RuO2 becomes active with 100% conversion, while in-situ XRD reveals no changes in the RuO2 diffraction pattern. At this stage selective roughening of apical RuO2 facets was observed by scanning electron microscopy (SEM). Seldom also single lateral facets are roughened. EDX indicated higher oxygen content in the following order: flat lateral facets > rough lateral facets > rough apical facets. Further, experiments in the packed bed reactor indicated oscillations in the CO2 formation rate. At even higher temperatures in reducing feed (533–543 K) the sample reduces to ruthenium metal according to XRD. The reduced particles exhibiting lower ignition temperature are very rough with cracks and deep star-shaped holes. An Arrhenius plot of the CO2 formation rate below the ignition temperature reveals the reduced samples to be significantly more active based on mass unit and shows lower apparent activation energy than the activated oxidized sample. Micro-spot X-ray photoelectron spectroscopy (XPS) and XPS microscopy experiments were carried out on a Ru(0001) single crystal exposed to oxygen at different temperature. Although low energy electron diffraction (LEED) images show a strong 1×1 pattern, the XPS data indicated a wide lateral inhomogeneity with different degree of oxygen dissolved in the subsurface layers. All these and the literature data are discussed in the context of different active states and transport issues, and the metastable nature of a phase mixture under conditions of high catalytic activity.

Reactivity of Oxide Precursor States on Ru(0001)

Smooth and defect-rich Ru(0001) surfaces prepared under ultrahigh-vacuum (UHV) conditions have been loaded with oxygen under high-pressure (p </= 1 bar) and low-temperature (T < 550 K) conditions. On these surfaces the CO oxidation reaction has been investigated by means of thermal desorption spectroscopy (TDS), ultraviolet photoelectron spectroscopy (UPS) and reactive molecular beam scattering (RMBS). Both surfaces are oxide-free and exhibit a high reactivity. The maximum CO/CO(2) conversion probability observed for a defect-rich Ru(0001) surface amounts to 6 x 10(-3) and is comparable to that of a surface covered with rutile RuO(2)(110) domains. RMBS experiments led to the identification of three different reaction channels. The first and second channel is related to CO adsorbing at oxygen-free defect sites and follow the Langmuir-Hinshelwood mechanism. Whereas the first reaction channel is already observed at room temperature, the second is thermally activated, contributing to the CO(2) yield only for reaction temperatures above 400 K. The third channel is due to the recombination of CO molecules with oxygen atoms located in smooth areas of the surface undisturbed by defects. This reaction channel is thermally activated as well.

Reactivity of the 4d Transition Metals toward N Hydrogenation and NH Dissociation: A DFT-Based HSAB Analysis

The density functional theory (DFT) based hard-soft acid-base (HSAB) reactivity indices, including the electrophilicity index, have been successfully applied to many areas of molecular chemistry. In this work we test the applicability of such an approach to fundamental surface chemistry. We have considered, as prototypical surface reactions, both the hydrogenation of atomic nitrogen and the dissociative adsorption of the NH molecular radical. By use of a DFT methodology, the minimum energy reaction pathways, and corresponding reaction barriers, of the above reactions over Zr(001), Nb(110), Mo(110), Tc(001), Ru(001), Rh(111), and Pd(111) have been determined. By consideration of the chemical potential and chemical hardness of the surface metal atoms, and the principle of electronegativity equalization, it is found that the charge transferred to the NH radical during the process of dissociative adsorption correlates very well with that determined by Mulliken population analysis. Furthermore, it is found that the stability of the NH/surface transition state complex relates directly to this charge transfer and that the trend in transition state stability predicted by a HSAB treatment correlates very strongly with that determined by DFT calculations. With regards to N hydrogenation, we find that during the course of the reaction, H loses cohesion to the surface, as it must migrate from a 3-fold hollow site to either a bridge or top site, to react with N. Partial density of states (PDOS) and Mulliken population analysis reveal that this loss of bonding is accompanied by charge transfer from H to the surface metal atoms. Moreover, by simple modeling, we show that the reaction barriers are directly proportional to this mandatory charge transfer. Indeed, it is found that the reaction barriers correlate very well with the electrophilicity index of the metal atoms.

Ultrathin Rh films on Ru(0001): Oxidation in confinement

Ultrathin rhodium films with a thickness ranging from 1 to a few monolayers were deposited on a single-crystal Ru(0001) surface in order to investigate the oxidation behavior of ultrathin epitaxial films on a dissimilar substrate. It is found that rhodium grows on Ru(0001) initially layer by layer, adapting the in-plane lattice parameters of Ru(0001). When exposing Rh films to oxygen environment (approximately 4.8 x 10(6) L O2 exposure) at 660 K, 2-4 ML Rh films form a surface oxide composed of (9 x 9) O-Rh-O trilayers. Quite in contrast, oxidation of the 1 ML RhRu(0001) film leads to a poorly ordered oxide with a rutile structure reminiscent of RuO2(110) on Ru(0001). The oxidized 1 ML RhRu(0001) film contains much more oxygen than the oxidized thicker Rh films. Lower temperatures (535 K) and high doses of oxygen lead to a (1 x 1)-O overlayer on the 1 ML RhRu(0001) surface, whose atomic geometry resembles closely that of the (1 x 1)-O phase on clean Ru(0001).

Ru(0001) Model Catalyst under Oxidizing and Reducing Reaction Conditions: In-Situ High-Pressure Surface X-ray Diffraction Study

With surface X-ray diffraction (SXRD) using a high-pressure reaction chamber we investigated in-situ the oxidation of the Ru(0001) model catalyst under various reaction conditions, starting from a strongly oxidizing environment to reaction conditions typical for CO oxidation. With a mixture of O(2) and CO (stoichiometry, 2:1) the partial pressure of oxygen has to be increased to 20 mbar to form the catalytically active RuO(2)(110) oxide film, while in pure oxygen environment a pressure of 10(-5) mbar is already sufficient to oxidize the Ru(0001) surface. For preparation temperatures in the range of 550-630 K a self-limiting RuO(2)(110) film is produced with a thickness of 1.6 nm. The RuO(2)(110) film grows self-acceleratedly after an induction period. The RuO(2) films on Ru(0001) can readily be reduced by H(2) and CO exposures at 415 K, without an induction period.

Identification of Subsurface Oxygen Species Created during Oxidation of Ru(0001)

The oxidation states formed during low-temperature oxidation (T < 500 K) of a Ru(0001) surface are identified with photoelectron spectromicroscopy and thermal desorption (TD) spectroscopy. Adsorption and consecutive incorporation of oxygen are studied following the distinct chemical shifts of the Ru 3d(5/2) core levels of the two topmost Ru layers. The evolution of the Ru 3d(5/2) spectra with oxygen exposure at 475 K and the corresponding O2 desorption spectra reveal that about 2 ML of oxygen incorporate into the subsurface region, residing between the first and second Ru layer. Our results suggest that the subsurface oxygen binds to the first and second layer Ru atoms, yielding a metastable surface "oxide", which represents the oxidation state of an atomically well ordered Ru(0001) surface under low-temperature oxidation conditions. Accumulation of more than 3 ML of oxygen is possible via defect-promoted penetration below the second layer when the initial Ru(0001) surface is disordered. Despite its higher capacity for oxygen accumulation, also the disordered Ru surface does not show features characteristic for the crystalline RuO2 islands. Development of lateral heterogeneity in the oxygen concentration is evidenced by the Ru 3d(5/2) images and microspot spectra after the onset of oxygen incorporation, which becomes very pronounced when the oxidation is carried out at T > 550 K. This is attributed to facilitated O incorporation and oxide nucleation in microregions with a high density of defects.

Oxide-free oxygen incorporation into Ru(0001)

A smooth Ru(0001) surface prepared under ultra-high vacuum conditions has been loaded with oxygen under high-pressure (p approximately 1 bar) and low-temperature (T < 600 K) conditions. Oxygen phases created in this way have been investigated by means of thermal desorption spectroscopy, low-energy electron diffraction, and ultraviolet photoelectron spectroscopy. The exposure procedures applied lead to oxygen incorporation into the subsurface region without creation of RuO2 domains. For oxygen exposures ranging from 10(11) to 10(14) L oxygen contents up to about 4 monolayer equivalent could be achieved. The oxygen incorporation is thermally activated. The CO oxidation reaction conducted at mild temperatures (T < 500 K) at a sample loaded with subsurface oxygen reaches CO --> CO2 conversion probabilities of 10(-3).

A Systematic Study of CO Oxidation on Metals and Metal Oxides: Density Functional Theory Calculations

CO oxidation on Ru(0001), Rh(111), Pd(111), Os(0001), Ir(111), Pt(111), and their corresponding metal oxides is studied using density functional theory. It is found that (i) the reactivity of metal oxide is generally higher than that of the corresponding metal, and (ii) on both metals and metal oxides, the higher the chemisorption energy is in the initial state, the larger the reaction barrier. The barriers are further analyzed by decomposing them into electronic and geometric effects, and the higher reactivity of metal oxides is attributed mainly to the surface geometric effect. Moreover, the electronic effect on both metals and metal oxides follows the same pattern: the shorter the OC-O bond distance in the TS, the higher the barrier.

First-principles atomistic thermodynamics for oxidation catalysis: surface phase diagrams and catalytically interesting regions

Present knowledge of the function of materials is largely based on studies (experimental and theoretical) that are performed at low temperatures and ultralow pressures. However, the majority of everyday applications, like, e.g., catalysis, operate at atmospheric pressures and temperatures at or higher than 300 K. Here we employ ab initio, atomistic thermodynamics to construct a phase diagram of surface structures in the (T,p) space from ultrahigh vacuum to technically relevant pressures and temperatures. We emphasize the value of such phase diagrams as well as the importance of the reaction kinetics that may be crucial, e.g., close to phase boundaries.

Heterogeneous catalysis on the atomic scale

Information about the elementary processes underlying heterogeneous catalysis may be obtained by investigating well-defined single crystal surfaces. The success of this "surface science" approach for "'real" catalysis can be demonstrated, for example, with ammonia synthesis. The progress of catalytic reactions can be observed on an atomic scale by applying scanning tunneling microscopy and other surface physical techniques, as is shown with different examples in this paper: CO oxidation on a Pt(111) surface proceeds preferentially along the boundaries between adsorbed O and CO patches. Ru is practically inactive for the same reaction under lower pressure conditions but is transformed into RuO2 under atmospheric pressure conditions, where part of the surface Ru atoms function as coordinatively unsaturated sites (cus). In contrast, in the hydrogen oxidation reaction on Pt(111), an autocatalytic reaction step comes into prominence, and is responsible for the formation of propagating concentration patterns on the surface as a characteristic of nonlinear dynamics. Additionally, the limits of the concept of thermal equilibrium in surface rate processes are explored by applying ultrafast (femtosecond) laser techniques.