Structure and Stability of a High-Coverage (1 x 1) Oxygen Phase on Ru(0001)

Interfacial catalytic materials; challenge for inorganic synthetic chemistry

Interfacial catalysts are indispensable functional materials in the energy transformation. The traditional empirical search strategies reach their potential. Knowledge-based approaches have not been able to deliver innovative and scalable solutions. Following a short analysis of the origin of these shortcomings a fresh attempt on the material challenge of catalysis is proposed. The approach combines functional understanding of material dynamics derived from operando analysis with digital catalysis science guiding the exploration of non-linear interactions of material genes to catalytic functions. This critically requires the ingenuity of the synthetic inorganic chemist to let us understand the reactivity of well-defined materials under the specific conditions of catalytic operation. It is the understanding of how the kinetics of phase changes brings about and destroys active sites in catalytic materials that forms the basis of realistic material concepts. A rigorous prediction and engineering of these processes may not be possible due to the complexity of options involved.

Exfoliating silica bilayers via intercalation at the silica/transition metal interface

The growth of the silica (SiO₂) bilayer (BL) films on transition metal (TM) surfaces creates a new class of two-dimensional (2D) crystalline, self-contained materials that interact weakly with the TM substrate. The BL-silica/TM heterojunction has shown unique physical and chemical properties that can lead to new chemical reaction mechanisms under the sub-nm confinement and broad potential applications ranging from surface protection, nano transistors, molecular sieves to nuclear waste removal. Novel applications of BL-silica can be further explored as a constituent of van der Waals assembly of 2D materials. Key to these applications is an unmet technical challenge to exfoliate and transfer BL-silica films in a large area from one substrate to another without material damage. In this study, we propose a new exfoliation mechanism based on gas molecule intercalation from density functional theory studies of the BL-silica/TM heterojunction. We found that the intercalation of O atoms and CO molecules at the BL-silica/TM interface weakens the BL-silica-TM hybridization, which results in an exponential decrease of the exfoliation energy against the interface distance as the coverage of interfacial species increases. This new intercalation mechanism opens up the opportunity for non-damaging exfoliation and transfer of large area silica bilayers.

Ultrafast Diffusion at the Onset of Growth: O/Ru(0001)

Nanoscopic clustering in a 2D disordered phase is observed for oxygen on Ru(0001) at low coverages and high temperatures. We study the coexistence of quasistatic clusters (with a characteristic length of ∼9  Å) and highly mobile atomic oxygen which diffuses between the energy-inequivalent, threefold hollow sites of the substrate. We determine a surprisingly low activation energy for diffusion of 385±20  meV. The minimum of the O-O interadsorbate potential appears to be at lower separations than previously reported.

Structural investigation of ternary PdRuM (M = Pt, Rh, or Ir) nanoparticles using first-principles calculations

We perform first-principles calculations and Monte Carlo sampling to investigate the structures of ternary PdRuM (M = Pt, Rh, or Ir) nanoparticles (NPs) with respect to three different spherical shapes. The morphologies include hexagonal close-packed (hcp), truncated-octahedral (fcc), and icosahedral (Ih, fcc) shapes with 57, 55, and 55 atoms, respectively. The calculations show that the atomic position is dominant in determining the stability of the ternary NPs. For bare ternary NPs, Pd and Ru atoms favor a location on the vertex sites and the core, respectively, which can be understood by the surface energy of the corresponding slab models. For single-crystalline NPs, the binary shell could be either a solid solution or a segregation alloy depending on composition and morphology. However, polycrystalline Ih NPs only form segregated binary shells surrounding the Ru core. Such configurations tend to minimize the surface lattice to gain more energy from the d orbital of the transition metals. In addition to the bare NPs, we study the oxidized ternary NPs. The results show that the Ru atoms penetrate outwards from the core to the surface reducing the oxidation formation energy. Furthermore, oxygen adsorption facilitates Pt, Pd, and Pd penetration into the PdRuPt, PdRuRh, and PdRuIr NPs, respectively. Most of the oxide shells are a solid solution, except for the PdRuRh NP with an Ih shape, which is found to be in a segregation shell. The free energy calculation reveals that the pure hcp NPs are thermodynamically unstable under oxygen-rich conditions. This work clearly demonstrates the structural trends of small ternary NPs and their oxidation, unveiling that the structural trends can be understood by the surface formation energy and the interplay between adsorbent and adsorbing oxygen atoms.

First-principles calculations and Monte Carlo sampling was used to investigate the structures of PdRuM (M = Pt, Rh, or Ir) nanoparticles with three different spherical shapes. Oxygen adsorption facilitates Ru penetration from core to surface sites.

Quantifying Confidence in DFT-Predicted Surface Pourbaix Diagrams of Transition-Metal Electrode-Electrolyte Interfaces

Density functional theory (DFT) calculations have been widely used to predict the activity of catalysts based on the free energies of reaction intermediates. The incorporation of the state of the catalyst surface under the electrochemical operating conditions while constructing the free-energy diagram is crucial, without which even trends in activity predictions could be imprecisely captured. Surface Pourbaix diagrams indicate the surface state as a function of the pH and the potential. In this work, we utilize error-estimation capabilities within the Bayesian ensemble error functional with van der Waals correlations exchange correlation functional as an ensemble approach to propagate the uncertainty associated with the adsorption energetics in the construction of Pourbaix diagrams. Within this approach, surface-transition phase boundaries are no longer sharp and are therefore associated with a finite width. We determine the surface phase diagram for several transition metals under reaction conditions and electrode potentials relevant for the oxygen reduction reaction. We observe that our surface phase predictions for most predominant species are in good agreement with cyclic voltammetry experiments and prior DFT studies. We use the OH* intermediate for comparing adsorption characteristics on Pt(111), Pt(100), Pd(111), Ir(111), Rh(111), and Ru(0001) since it has been shown to have a higher prediction efficiency relative to O*, and find the trend Ru > Rh > Ir > Pt > Pd for (111) metal facets, where Ru binds OH* the strongest. We robustly predict the likely surface phase as a function of reaction conditions by associating confidence values for quantifying the confidence in predictions within the Pourbaix diagram. We define a confidence quantifying metric, using which certain experimentally observed surface phases and peak assignments can be better rationalized. The probabilistic approach enables a more accurate determination of the surface structure and can readily be incorporated in computational studies for better understanding the catalyst surface under operating conditions.

The Nanostructuring of Atomically Flat Ru(0001) upon Oxidation and Reduction

The O/Ru(0001) system is widely studied due to its rich phase variety of various stoichiometry and atomic arrangements, including the formation of a RuO₂/Ru(0001) oxide layer. Apart from homogeneous ruthenium surfaces in certain oxidation states, also strongly heterogeneous surfaces can exist due to oxidation state's variation at the nanoscale. We report on a scanning tunneling microscopy (STM) study of the nanostructuring of the oxidized Ru(0001) surface as a result of its interaction with molecular oxygen at elevated temperatures and subsequent reduction of a resulting RuO₂ film by CO or HCl molecules from the gas phase in high-vacuum environment.

Factors controlling the CO intercalation of h-BN overlayers on Ru(0001)

Critical factors influencing the CO intercalation of h-BN were investigated including CO partial pressure, h-BN coverage, and oxygen pre-adsorption on Ru.

Role of RuO2(100) in surface oxidation and CO oxidation catalysis on Ru(0001)

Oxidation of Ru(0001) induces the simultaneous formation of RuO₂(100) and RuO₂(110) and a structure-sensitive oxygen spillover during CO oxidation.

Kinetic analysis of interaction between N atoms and O-covered Ru(0001)

Eley-Rideal (ER) reactions involving neutral atoms heavier than hydrogen reacting with adsorbed atoms of similar mass were first observed in recent molecular beam experiments by Zaharia et al. [Phys. Rev. Lett. 113, 053201 (2014)]. Through analysis of two types of measurements, they obtained different estimations for the N-O ER reaction cross section, one of which is unexpectedly high. This was qualitatively accounted for by invoking a secondary effect whereby the presence of N adatoms on the surface acted to "shield" O adatoms from prompt recombinative desorption. We apply a rate equation model that includes two ER processes involving different adsorbed species (N-Oad and N-Nad) and an N-adsorption process to the full-beam exposure subset of the experimental data in order to study the reaction kinetics. Values for the individual reaction cross sections are derived. The measured N2 response can be well described by the model, but it is insufficient to completely describe the NO response. Modeling of different exposures is used to evaluate the qualitative picture presented by Zaharia et al.

Predicting Promoter-Induced Bond Activation on Solid Catalysts Using Elementary Bond Orders

In this Letter, we examine bond activation induced by nonmetal surface promoters in the context of dehydrogenation reactions. We use C-H bond activation in methane dehydrogenation on transition metals as an example to understand the origin of the promoting or poisoning effect of nonmetals. The electronic structure of the surface and the bond order of the promoter are found to establish all trends in bond activation. On the basis of these results, we develop a predictive model that successfully describes the energetics of C-H, O-H, and N-H bond activation across a range of reactions. For a given reaction step, a single data point determines whether a nonmetal will promote bond activation or poison the surface and by how much. We show how our model leads to general insights that can be directly used to predict bond activation energetics on transition metal sulfides and oxides, which can be perceived as promoted surfaces. These results can then be directly used in studies on full catalytic pathways.

Surface resonances in electron reflection from overlayers

Electron scattering by oxygen monolayers on the Ru(0 0 0 1) surface is studied both experimentally and theoretically. Sharp transmission resonances at low energies are revealed and established to originate from critical points of a special kind in the complex band structure of the substrate. Electron reflection from the clean and oxidized Ru(0 0 0 1) is measured for kinetic energies up to 40 eV at normal incidence for oxygen coverages of 1/4, 1/2, 3/4, and one monolayer. The reflection spectra R(E) are analyzed using a Bloch-waves based ab initio scattering theory. In addition to the substrate-induced resonances the reconstructed (2 × 1) and (2 × 2) surfaces show surface resonances due to pre-emergent secondary diffraction beams. The R(E) spectra are shown to give unambiguous evidence of the hcp stacking of the oxygen layer.

Insights into the gas phase oxidation of Ru(0001) on the mesoscopic scale using molecular oxygen

We present an extensive mesoscale study of the initial gas phase oxidation of Ru(0001), employing in situ low-energy electron microscopy (LEEM), micro low-energy electron diffraction (μ-LEED) and scanning tunneling microscopy (STM).

On the role of the surface oxygen species during A–H (A = C, N, O) bond activation: a density functional theory study

During A–H (A = C, N, O) bond cleavage on O* or OH* pre-covered (111) surfaces, the oxygen species play the role of modifying the reaction energy by changing the species involved in the initial and final states of the reaction.

Initial stages of FeO growth on Ru(0001)

We study how FeO wüstite films on Ru(0001) grow by oxygen-assisted molecular beam epitaxy at elevated temperatures (800–900 K). The nucleation and growth of FeO islands are observed in real time by low-energy electron microscopy (LEEM). When the growth is performed in an oxygen pressure of 10(−6) Torr, the islands are of bilayer thickness (Fe–O–Fe–O). In contrast, under a pressure of 10(−8) Torr, the islands are a single FeO layer thick. We propose that the film thickness is controlled by the concentration of oxygen adsorbed on the Ru. More specifically, when monolayer growth increases the adsorbed oxygen concentration above a limiting value, its growth is suppressed. Increasing the temperature at a fixed oxygen pressure decreases the density of FeO islands. However, the nucleation density is not a monotonic function of oxygen pressure.

Real-time microscopy of graphene growth on epitaxial metal films: role of template thickness and strain

Epitaxial transition metal films have recently been introduced as substrates for the scalable synthesis of transferable graphene. Here, real-time microscopy is used to study graphene growth on epitaxial Ru films on sapphire. At high temperatures, high-quality graphene grows in macroscopic (>100 μm) domains to full surface coverage. Graphene nucleation and growth characteristics on thin (100 nm) Ru films are consistent with a pure surface chemical vapor deposition process, without detectable contributions from C segregation. Experiments on thicker (1 μm) films show a systematic suppression of the C uptake into the metal to levels substantially below those expected from bulk C solubility data, consistent with a strain-induced reduction of the C solubility due to gas bubbles acting as stressors in the epitaxial Ru films. The results identify two powerful approaches--i) limiting the template thickness and ii) tuning the interstitial C solubility via strain--for controlling graphene growth on metals with high C solubility, such as Ru, Pt, Rh, Co, and Ni.

Thin silica films on Ru(0001): monolayer, bilayer and three-dimensional networks of [SiO4] tetrahedra

The atomic structure of thin silica films grown over a Ru(0001) substrate was studied by X-ray photoelectron spectroscopy, infrared reflection absorption spectroscopy, low energy electron diffraction, helium ion scattering spectroscopy, CO temperature programmed desorption, and scanning tunneling microscopy in combination with density functional theory calculations. The films were prepared by Si vapor deposition and subsequent oxidation at high temperatures. The silica film first grows as a monolayer of corner-sharing [SiO(4)] tetrahedra strongly bonded to the Ru(0001) surface through the Si-O-Ru linkages. At increasing amounts of Si, the film forms a bilayer of corner-sharing [SiO(4)] tetrahedra which is weakly bonded to Ru(0001). The bilayer film can be grown in either the crystalline or vitreous state, or both coexisting. Further increasing the film thickness leads to the formation of vitreous silica exhibiting a three-dimensional network of [SiO(4)]. The principal structure of the films can be monitored by infrared spectroscopy, as each structure shows a characteristic vibrational band, i.e., ∼1135 cm(-1) for a monolayer film, ∼1300 cm(-1) for the bilayer structures, and ∼1250 cm(-1) for the bulk-like vitreous silica.

Chemical vapor deposition and etching of high-quality monolayer hexagonal boron nitride films

The growth of large-area hexagonal boron nitride (h-BN) monolayers on catalytic metal substrates is a topic of scientific and technological interest. We have used real-time microscopy during the growth process to study h-BN chemical vapor deposition (CVD) from borazine on Ru(0001) single crystals and thin films. At low borazine pressures, individual h-BN domains nucleate sparsely, grow to macroscopic dimensions, and coalescence to form a closed monolayer film. A quantitative analysis shows borazine adsorption and dissociation predominantly on Ru, with the h-BN covered areas being at least 100 times less reactive. We establish strong effects of hydrogen added to the CVD precursor gas in controlling the in-plane expansion and morphology of the growing h-BN domains. High-temperature exposure of h-BN/Ru to pure hydrogen causes the controlled edge detachment of B and N and can be used as a clean etching process for h-BN on metals.

Reaction kinetics and polarization modulation infrared reflection absorption spectroscopy investigations of CO oxidation over planar Pt-group model catalysts

Microscopic and spectroscopic techniques at near-atmospheric pressures have been used in recent years to investigate CO oxidation over Pt-group metals. New insights have been obtained that allow broadening of the understanding of this reaction beyond the ultrahigh vacuum regime where it is well-understood. However, new issues also have arisen that need clarification. In this paper, we review recent work in our laboratory, using polarization modulation infrared reflection absorption spectroscopy (PM-IRAS) and reaction kinetics measurements from ultrahigh vacuum to near-atmospheric pressures. These studies reveal a continuity of this reaction with respect to pressure over Pt, Pd, and Rh; that is, Langmuir-Hinshelwood kinetics is exhibited over a wide pressure range with no apparent "pressure gap". The difference between Ru(0001) and other noble metals is well-understood with respect to higher oxygen binding energies and reduced CO inhibition. It is concluded that for all Pt-group metals the most active phase is one saturated with chemisorbed oxygen and with low CO coverages. The significance of oxide phases under most industrially relevant catalytic conditions suggested recently in the literature is not consistent with the experimental data.

STM, SECPM, AFM and Electrochemistry on Single Crystalline Surfaces

Scanning probe microscopy (SPM) techniques have had a great impact on research fields of surface science and nanotechnology during the last decades. They are used to investigate surfaces with scanning ranges between several 100 μm down to atomic resolution. Depending on experimental conditions, and the interaction forces between probe and sample, different SPM techniques allow mapping of different surface properties. In this work, scanning tunneling microscopy (STM) in air and under electrochemical conditions (EC-STM), atomic force microscopy (AFM) in air and scanning electrochemical potential microscopy (SECPM) under electrochemical conditions, were used to study different single crystalline surfaces in electrochemistry. Especially SECPM offers potentially new insights into the solid-liquid interface by providing the possibility to image the potential distribution of the surface, with a resolution that is comparable to STM. In electrocatalysis, nanostructured catalysts supported on different electrode materials often show behavior different from their bulk electrodes. This was experimentally and theoretically shown for several combinations and recently on Pt on Au(111) towards fuel cell relevant reactions. For these investigations single crystals often provide accurate and well defined reference and support systems. We will show heteroepitaxially grown Ru, Ir and Rh single crystalline surface films and bulk Au single crystals with different orientations under electrochemical conditions. Image studies from all three different SPM methods will be presented and compared to electrochemical data obtained by cyclic voltammetry in acidic media. The quality of the single crystalline supports will be verified by the SPM images and the cyclic voltammograms. Furthermore, an outlook will be presented on how such supports can be used in electrocatalytic studies.

Oxygen-coverage effects on molecular dissociations at a Pt metal surface

The effects of adsorbate coverage on catalytic surface reactions are not well understood. Here, we contrast the rates of O2 and NO2 dissociations, two competing reactions in NO oxidation catalysis, versus oxygen coverage at a Pt(111) surface. In situ x-ray photoelectron spectroscopy experiments show that the NO2 dissociation rate is less sensitive to O coverage than is O2. Density-functional theory simulations reveal an NO2 reaction pathway that is more adaptable to an increasingly crowded surface than is O2 dissociation. While the rates are comparable at low coverage, NO2 dissociation is many orders of magnitude faster at O coverages typical of NO oxidation catalysis.

Monomer adsorption on equilateral triangular lattices with repulsive first-neighbor interactions

A model of monomer adsorption on infinitely long, finite-width M equilateral triangular lattices with nonperiodic boundaries is presented. The study includes adsorbate-adsorbate first- and second-neighbor interactions with results obtained for repulsive first neighbors. The matrix method and numerical algorithms presented here allow determination of the occupational characteristics of the adsorption crystallization phases, which fit exact analytic expressions in the width M of the lattice. The limit as M approaches infinity provides the complete energy phase diagram for the infinite two-dimensional surface and recovers the results obtained by different methods that were often applied only in restricted energy regions of the phase diagram. The ordered phases are (2x1), (2x2), (3x1), (square root of 3xsquare root of 3) R30 degrees, and the complementary phases of (2x2) and (square root of 3xsquare root of 3) R30 degrees. Comparison is made with other theoretical studies and with experimental observations on adsorption systems consistent with the limitations of the model. In some cases, comparison with experimental data yields bounds on the interaction energies between adsorbates. On the basis of the model, suggestions are made on the manner in which to conduct relatively low temperature experiments to allow determination of most, if not all, of the interaction energies from the knowledge of the sequences of phases and the conditions prevailing at the transitions between phases.

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.

Initial Growth of the Water Layer on (1 × 1)-Oxygen-Covered Ru(0001) in Comparison with that on Bare Ru(0001)

The initial growth of a water (D2O) layer on (1 x 1)-oxygen-covered Ru(0001) has been studied in comparison with that on bare Ru(0001) by means of temperature-programmed desorption (TPD) and infrared reflection absorption spectroscopy (IRAS). Although water molecules adsorbed on both bare and (1 x 1)-oxygen-covered Ru(0001) commonly tend to form hydrogen bonds with each other when mobility occurs upon heating, the TPD and IRAS measurements for the two surfaces exhibit distinct differences. On (1 x 1)-oxygen-covered Ru(0001), most of the D2O molecules were desorbed with a peak at 160 K, even at submonolayer coverage, as condensed water desorption. The vibration spectra of adsorbed D2O also showed broad peaks such as a condensed water phase, from the beginning of low coverage. For submonolayer coverage, in addition, we found a characteristic O-D stretching mode at around 2650 cm(-1), which is never clearly observed for D2O on bare Ru(0001). Thus, we propose a distinctive water adsorption structure on (1 x 1)-oxygen-covered Ru(0001) and discuss its influence on water layer growth in comparison with the case of D2O on bare Ru(0001).