A kinetic Monte Carlo/UBI-QEP model of O2 adsorption on fcc (1 1 1) metal surfaces

The previously developed kinetic Monte Carlo model of molecular oxygen adsorption on fcc (1 0 0) metal surfaces has been extended to fcc (1 1 1) surfaces. The model treats uniformly all elementary steps of the process—O₂ adsorption, dissociation, recombination, desorption, and atomic oxygen hopping—at various coverages and temperatures. The model employs the unity bond index—quadratic exponential potential (UBI-QEP) formalism to calculate coverage-dependent energetics (atomic and molecular binding energies and activation barriers of elementary steps) and a Metropolis-type algorithm including the Arrhenius-type reaction rates to calculate coverage- and temperature-dependent features, particularly the adsorbate distribution over the surface. Optimal values of non-energetic model parameters (the spatial constraint, a travel distance of “hot” atoms, attempt frequencies of elementary steps) have been chosen. Proper modifications of the fcc (1 0 0) model have been made to reflect structural differences in the fcc (1 1 1) surface, in particular the presence of two different hollow sites (fcc and hcp). Detailed simulations were performed for molecular oxygen adsorption on Ni(1 1 1). We found that at very low coverages, only O₂ adsorption and dissociation were effective, while O₂ desorption and O₂ and O diffusion practically did not occur. At a certain O + O₂ coverage, the O₂ dissociation becomes the fastest process with a rate one-two orders of magnitude higher than adsorption. Dissociation continuously slows down due to an increase in the activation energy of dissociation and due to the exhaustion of free sites. The binding energies of both molecular and atomic oxygen decrease with coverage, and this leads to greater mobility of atomic oxygen and more pronounced desorption of molecular oxygen. Saturation is observed when the number of adsorbed molecules becomes approximately equal to the number of desorbed molecules. Simulated coverage dependences of the sticking probability and of the atomic binding energy are in reasonable agreement with experimental data. From comparison with the results of the previous work, it appears that the binding energy profiles for Ni(1 1 1) and Ni(1 0 0) have similar shapes, although at any coverage the absolute values of the oxygen binding energy are higher for the (1 0 0) surface. For metals other than Ni, particularly Pt, the model projections were found to be too parameter-dependent and therefore less certain. In such cases further model developments are needed, and we briefly comment on this situation.     

Room temperature observation of nitric oxide on Ir(1 1 1) by scanning tunneling microscopy

Room temperature (RT) adsorption of nitric oxide (NO) on Ir(1 1 1) was studied by scanning tunneling microscopy (STM). At low exposures, NO molecules can not be imaged by STM, because at RT the diffusion of NO is much faster than the STM scanning speed. At high exposures near the saturation coverage, however, a well-ordered 2 × 2 structure is observed. The coverage of the major 2 × 2 species is 0.25 and they can be assigned to the NO molecules adsorbed on the Ir ontop sites. A small number of less bright spots are assigned to nitrogen atoms produced by dissociation. Their number increases by annealing the NO-saturated surface at 380 K. A small number of another dissociation product, oxygen, are observed as black lines, indicating that the diffusion of oxygen atoms is fast. Scratch-like noise features were also detected by the STM, which suggests that a mobile precursor state exists, which was clearly shown by the effects of electron irradiation from the STM tip. These results are consistent with the previous molecular beam studies. Hopping of the 2 × 2 ordered NO species was frequently observed at the anti-phase domain boundaries and edges of the 2 × 2 islands.     

In situ oxidation study of MgO(1 0 0) supported Pd nanoparticles

We have investigated the oxidation behavior of Pd nanoparticles grown epitaxially on MgO(1 0 0) single crystal substrates. We find that the interaction of oxygen with octahedral Pd nanoparticles at 500 K can be subdivided in three stages: above 10⁻⁶ mbar O₂ pressure, the particles start to flatten; above 10⁻³ mbar, the particles begin to shrink laterally and to be less truncated at the corners. The formation of epitaxial bulk PdO sets in at oxygen pressures above 0.1 mbar, which is accompanied by a continuous shrinkage of the Pd particles. Our results point to a novel nanoparticle oxidation mechanism: the Pd particles act as dissociation centers for O₂ and serve at the same time as source for Pd atoms resulting in epitaxial PdO growth on MgO(1 0 0).     

Adsorption behavior and reaction properties of NO and CO on Rh(1 1 1)

Reactions between NO and CO on Rh(1 1 1) surfaces were investigated using infrared reflection absorption spectroscopy, X-ray photoelectron spectroscopy, and temperature-programmed desorption. NO adsorbed on the fcc, atop, and hcp sites in that order, whereas CO adsorbed initially on the atop sites and then on the hollow (fcc + hcp) sites. The results of experiments with NO exposure on CO-preadsorbed Rh(1 1 1) surfaces indicated that the adsorption of NO on the hcp sites was inhibited by preadsorption of CO on the atop sites, and NO adsorption on the atop and fcc sites was inhibited by CO preadsorbed on each type of site, which indicates that NO and CO competitively adsorbed on Rh(1 1 1). From a Rh(1 1 1) surface with coadsorbed NO and CO, N₂ was produced from the dissociation of fcc-NO, and CO₂ was formed by the reaction of adsorbed CO with atomic oxygen from dissociated fcc-NO. The CO₂ production increased remarkably in the presence of hollow-CO. Coverage of fcc-NO and hollow-CO on Rh(1 1 1) depended on the composition ratio of the NO/CO gas mixture, and a gas mixture with NO/CO ? 1/2 was required for the co-existence of fcc-NO and hollow-CO at 273 K.     

Hydrogen adsorption and diffusion on Pd(1 1 1)

The adsorption, diffusion and ordering of hydrogen on Pd(1 1 1) was studied by scanning tunneling microscopy in the temperature range of 37-90 K. At low coverage isolated hydrogen atoms were observed. They formed √3×√3-1H islands as the coverage increased. Above 1/3 monolayer (ML) coverage areas of a new phase with √3×√3-2H structure were formed, with both structures coexisting between 1/3 and 2/3 ML. Finally a 1 × 1 structure was formed after high exposures of hydrogen above 50 K, with a coverage close to 1 ML. Atomically resolved images reveal that H binds to fcc hollow sites.     

Adsorption and dissociation of ammonia on Au(1 1 1) surface: A density functional theory study

Periodic density functional theory (DFT) calculations using plane waves had been performed to systematically investigate the stable adsorption amine and its dehydrogenated reaction on Au(1 1 1) surface. The equilibrium configuration including on top, bridge, and hollow (fcc and hcp) sites had been determined by relaxation of the system. The adsorption both NH₃ on top site and NH₂ on bridge site is favorable on Au(1 1 1) surface, while the adsorption of NH on hollow (fcc) site is preferred. The adsorbates are adsorbed on the gold surface with the interaction between p orbital of adsorbate and the d orbital of gold atoms. The interaction between adsorbate and gold slab is more evident on the first layer than on any others. Furthermore, the dissociation reaction of NH₃ on clean gold surface, as well as on the pre-covered oxygen atom and pre-covered hydroxyl group surface had been investigated. The results show that the dehydrogenated reaction energy barrier on the pre-covered oxygen gold surface is lower. The adsorbed O can promote the dehydrogenation of amine. Additionally, OH as the product of the NH₃ dissociation reaction participates in continuous dehydrogenation reaction, and the reaction energy barrier is the lowest (22.77 kJ/mol). The results indicated that OH_ads play a key role in the dehydrogenated reaction on Au(1 1 1) surface.     

Initial stages of oxidation of Cu(1 1 1)

Several surface analysis techniques were combined to study the initial stages of oxidation of Cu(1 1 1) surfaces exposed to O₂ at low pressure (<5 × 10⁻⁶ mbar) and room temperature. Scanning tunneling microscopy (STM) results show that the reactivity is governed by the restructuring of the Cu(1 1 1) surface. On the terraces, oxygen dissociative adsorption leads to the formation of isolated O adatoms and clusters weakly bound to the surface. The O adatoms are located in the fcc threefold hollow sites of the unrestructured terraces. Friedel oscillations with an amplitude lower than 5 pm have been measured around the adatoms. At step edges, surface restructuring is initiated and leads to the nucleation and growth of a two-dimensional disordered layer of oxide precursor. The electronic structure of this oxide layer is characterised by a band gap measured by scanning tunneling spectroscopy to be ∼1.5 eV wide. The growth of the oxide islands progresses by consumption of the upper metal terraces to form triangular indents. The extraction of the Cu atoms at this interface generates a preferential orientation of the interface along the close-packed directions of the metal. A second growth front corresponds to the step edges of the oxide islands and progresses above the lower metal terraces. This is where the excess Cu atoms extracted at the first growth front are incorporated. STM shows that the growing disordered oxide layer consists of units of hexagonal structure with a first nearest neighbour distance characteristic of a relaxed Cu-Cu distance (∼0.3 nm), consistent with local Cu₂O(1 1 1)-like elements. Exposure at 300 °C is necessary to form an ordered two-dimensional layer of oxide precursor. It forms the so-called “29” superstructure assigned to a periodic distorted Cu₂O(1 1 1)-like structure.     

First-principles calculations of hydrogen molecule adsorption on Ti (0 0 0 1)-(2 × 1) surface

Adsorption of H₂ molecule on the Ti (0 0 0 1)-(2 × 1) surface was studied by density functional theory with generalized gradient approximation (GGA). The parallel and vertical absorption cases were investigated in detail by adsorption energy and electronic structure analysis, we obtained three stable configurations of FCC-FCC (the two H atoms adsorption on the two adjacent fcc sites of Ti (0 0 0 1) surface, respectively), HCP-HCP (the two H atoms adsorption on the two adjacent hcp sites of Ti (0 0 0 1) surface, respectively) and FCC-HCP (the one H atom adsorption on the fcc site and the other adsorption on the near hcp site) based on the six different parallel adsorption sites after the H₂ molecule dissociates. However, all the end configurations of four vertical adsorption sites were unstable, H₂ molecule was very easy to desorb from Ti surface. The H-H bond breaking and Ti-H bond forming result from the H₂ molecule dissociation. H-H bond breaking length ranges from 1.9 Å to 2.3 Å for different adsorption configurations due to the strong Ti-H bond forming. The H₂ dissociative approach and the end stable configurations formation in parallel adsorption processes are attributed to the quantum mechanics steering effects.     

Initial stages of the autocatalytic oxidation of the InAs(0 0 1)-(4 × 2)/c(8 × 2) surface by molecular oxygen

The initial stages of oxidation of the In-rich InAs(0 0 1)-(4 × 2)/c(8 × 2) surface by molecular oxygen (O₂) were studied using scanning tunneling microscopy (STM) and density functional theory (DFT). It was shown that the O₂ dissociatively chemisorbs along the rows in the [1 1 0] direction on the InAs surface either by displacing the row-edge As atoms or by inserting between In atoms on the rows. The dissociative chemisorption is consistent with being autocatalytic: there is a high tendency to form oxygen chemisorption sites which grow in length along the rows in the [1 1 0] direction at preexisting oxygen chemisorption sites. The most common site size is about 21-24 Å in length at ∼25% ML coverage, representing 2-3 unit cell lengths in the [1 1 0] direction (the length of ∼5-6 In atoms on the row). The autocatalysis was confirmed by modeling the site distribution as non-Poisson. The autocatalysis and the low sticking probability (∼10⁻⁴) of O₂ on the InAs(0 0 1)-(4 × 2)/c(8 × 2) are consistent with activated dissociative chemisorption. The results show that is it critical to protect the InAs surface from oxygen during subsequent atomic layer deposition (ALD) or molecular beam epitaxy (MBE) oxide growth since oxygen will displace As atoms.     

Formation, stability and CO adsorption properties of PdAg/Pd(1 1 1) surface alloys

The formation, stability and CO adsorption properties of PdAg/Pd(1 1 1) surface alloys were investigated by X-ray photoelectron spectroscopy (XPS) and by adsorption of CO probe molecules, which was characterized by temperature-programmed desorption (TPD) and high resolution electron energy loss spectroscopy (HREELS). The PdAg/Pd(1 1 1) surface alloys were prepared by annealing (partly) Ag film covered Pd(1 1 1) surfaces, where the Ag films were deposited at room temperature. Surface alloy formation leads to a modification of the electronic properties, evidenced by core-level shifts (CLSs) of both the Pd(3d) and Ag(3d) signal, with the extent of the CLSs depending on both initial Ag coverage and annealing temperature. The role of Ag pre-coverage and annealing temperature on surface alloy formation is elucidated. For a monolayer Ag covered Pd(1 1 1) surface, surface alloy formation starts at ∼450 K, and the resulting surface alloy is stable upon annealing at temperatures between 600 and 800 K. CO TPD and HREELS measurements demonstrate that at 120 K CO is exclusively adsorbed on Pd surface atoms/Pd sites of the bimetallic surfaces, and that the CO adsorption behavior is dominated by geometric ensemble effects, with adsorption on threefold hollow Pd₃ sites being more stable than on Pd₂ bridge sites and finally Pd₁ a-top sites.     

Growth and thermal evolution of submonolayer Pt films on Ru(0 0 0 1) studied by STM

The growth of submonolayer Pt on Ru(0 0 0 1) has been studied with scanning tunneling microscopy. We focus on the island evolution depending on Pt coverage θ_Pt, growth temperature T_G and post-growth annealing temperature T_A. Dendritic trigonal Pt islands with atomically rough borders are observed at room temperature and moderate deposition rates of about 5 × 10⁻⁴ ML/s. Two types of orientation, rotated by 180° and strongly influenced by minute amounts of oxygen are observed which is ascribed to nucleation starting at either hcp or fcc hollow sites. The preference for fcc sites changes to hcp in the presence of about one percent of oxygen. At lower growth temperatures Pt islands show a more fractal shape. Generally, atomically rough island borders smooth down at elevated growth temperatures higher than 300 K, or equivalent annealing temperatures. Dendritic Pt islands, for example, transform into compact, almost hexagonal islands, indicating similar step energies of A- and B-type of steps. Depending on the Pt coverage the thermal evolution differs somewhat: While regular islands on Ru(0 0 0 1) are formed at low coverages, vacancy islands are observed close to completion of the Pt layer.     

Interaction of O2 with Pd single crystals in the range 1-150 Torr: Surface morphology transformations

The interaction of O₂ with Pd single crystals including Pd(1 1 1), Pd(1 1 0) and Pd(1 0 0) in the pressure range 1-150 Torr was studied using scanning tunneling microscopy (STM). The Pd single crystal surface morphologies were determined by the oxidation conditions: O₂ pressure, exposure time and treatment temperature. Oxygen dissolution into Pd metal followed by the formation of bulk oxide was observed. The dissolution of oxygen resulted in the increase of the inter-planar spacing between the first two layers, 9-14% increase after an exposure of Pd(1 1 1) to 10-25 Torr O₂ at 600 K for 10 min, and 10-20% increase after exposing Pd(1 1 0) and Pd(1 0 0) to 1 Torr O₂ at 600 K for 10 min. Elongated or semi-spherical oxide agglomerates along the steps nucleated and grew on both Pd(1 1 0) and Pd(1 0 0) surfaces after oxidation in 5-25 Torr O₂ at 600 K. When bulk PdO was formed, the single crystal surface was covered with semi-spherical agglomerates 2-4 nm in size, which tended to aggregate to form a “cauliflower-like” structure. The single crystal surface area also increased during oxidation.     

Ethylene dissociation on flat and stepped Ni(1 1 1): A combined STM and DFT study

The dissociative adsorption of ethylene (C₂H₄) on Ni(1 1 1) was studied by scanning tunneling microscopy (STM) and density functional theory (DFT) calculations. The STM studies reveal that ethylene decomposes exclusively at the step edges at room temperature. However, the step edge sites are poisoned by the reaction products and thus only a small brim of decomposed ethylene is formed. At 500 K decomposition on the (1 1 1) facets leads to a continuous growth of carbidic islands, which nucleate along the step edges.DFT calculations were performed for several intermediate steps in the decomposition of ethylene on both Ni(1 1 1) and the stepped Ni(2 1 1) surface. In general the Ni(2 1 1) surface is found to have a higher reactivity than the Ni(1 1 1) surface. Furthermore, the calculations show that the influence of step edge atoms is very different for the different reaction pathways. In particular the barrier for dissociation is lowered significantly more than the barrier for dehydrogenation, and this is of great importance for the bond-breaking selectivity of Ni surfaces.The influence of step edges was also probed by evaporating Ag onto the Ni(1 1 1) surface. STM shows that the room temperature evaporation leads to a step flow growth of Ag islands, and a subsequent annealing at 800 K causes the Ag atoms to completely wet the step edges of Ni(1 1 1). The blocking of the step edges is shown to prevent all decomposition of ethylene at room temperature, whereas the terrace site decomposition at 500 K is confirmed to be unaffected by the Ag atoms.Finally a high surface area NiAg alloy catalyst supported on MgAl₂O₄ was synthesized and tested in flow reactor measurements. The NiAg catalyst has a much lower activity for ethane hydrogenolysis than a similar Ni catalyst, which can be rationalized by the STM and DFT results.     

Sum frequency generation and density functional studies of CO-H interaction and hydrogen bulk dissolution on Pd(1 1 1)

CO-H interaction and H bulk dissolution on Pd(1 1 1) were studied by sum frequency generation (SFG) vibrational spectroscopy and density functional theory (DFT). The theoretical findings are particularly important to rationalize the experimentally observed mutual site blocking of CO and H and the effect of H dissolution on coadsorbate structures. Dissociative hydrogen adsorption on CO-precovered Pd(1 1 1) is impeded due to an activation barrier of ∼2.5 eV for a CO coverage of 0.75 ML, an effect which is maintained down to 0.33 ML CO. Preadsorbed hydrogen prevented CO adsorption at 100 K, while hydrogen was replaced from the surface by CO above 125 K. The temperature-dependent site blocking of hydrogen originates from the onset of hydrogen diffusion into the Pd bulk around 125 K, as shown by SFG and theoretical calculations using various approaches. When Pd(1 1 1) was exposed to 1:1 CO/H₂ mixtures at 100 K, on-top CO was absent in the SFG spectra although hydrogen occupies only threefold hollow sites on Pd(1 1 1). DFT attributes the absence of on-top CO to H atoms diffusing between hollow sites via bridge sites, thereby destabilizing neighboring on-top CO molecules. According to the calculations, the stretching frequency of bridge-bonded CO with a neighboring bridge-bonded hydrogen atom is redshifted by 16 cm⁻¹ when compared to bridging CO on the clean surface. Implications of the observed effects on hydrogenation reactions are discussed and compared to the C₂H₄-H coadsorption system.     

O2 dissociation on Pd(2 1 1) and Cu(2 1 1) surfaces

We have performed ab initio Density Functional Theory (DFT) based calculations to observe the reactivity of the Pd(2 1 1) and Cu(2 1 1) surfaces towards O₂. In order to properly address the adsorption dynamics, the static potential energy surface calculations have been complemented with first principles molecular dynamics calculations, which reveal interesting steering effects that complicate the dissociation dynamics. We have found that on both surfaces the step microfacets are very reactive and the dissociation of the O₂ molecule at room temperature occurs mostly on those sites.     

Structure analysis of W(0 0 1)2 × 1-O surface at room and liquid nitrogen temperatures

Surface structure of O-adsorbed W(0 0 1) surface after annealing to 1200 K has been analyzed by low energy electron diffraction at 77 K as well as at room temperature. The optimum structure has tungsten missing rows and oxygen double rows. Furthermore, the R-factor is minimized at the structure that O atoms are adsorbed on one of the two different threefold hollow sites of the (1 1 0) facet appearing on the W(0 0 1)2 × 1 with missing row. However, the results suggest that two domains of O atoms adsorbed on both the two different threefold hollow sites coexist. Then, I-V curves have been analyzed as a function of the mixing ratio of the two domains having different O adsorption sites at room and low temperatures. The energy difference between these two sites has been estimated to be 6.5 meV from the temperature dependence of the mixing ratio.     

Vibrationally-assisted dissociative adsorption of oxygen on Ru(0 0 0 1)-p(2 × 1)-O

Supersonic molecular beam technique combined with high resolution X-ray photoelectron spectroscopy using synchrotron radiation was applied to the study of the dynamics of dissociative adsorption of oxygen on Ru(0 0 0 1) surface in high coverage region. The Ru(0 0 0 1) surface pre-covered with oxygen atoms of 0.5 monolayer, which corresponds to the p(2 × 1)-O structure, was dosed to oxygen molecules with translational energy of 0.5 eV. Oxygen uptake was compared between the cases with and without the beam source heated in order to verify the effects of internal energy of oxygen. We found drastic enhancement in initial sticking probability of oxygen when the beam source was heated to 1400 K. We concluded that the enhancement of sticking probability is mainly caused by molecular vibrational excitation, indicating that dissociation barrier is located in the exit channel on potential energy surface.     

Photoelectron spectroscopy study of oxygen adsorption on Mo2C(0 0 0 1)

Oxygen adsorption on the α-Mo₂C(0 0 0 1) surface has been investigated with X-ray photoelectron spectroscopy and valence photoelectron spectroscopy utilizing synchrotron radiation. It is found that oxygen adsorbs dissociatively at room temperature, and the adsorbed oxygen atoms interact with both Mo and C atoms to form an oxycarbide layer. As the O-adsorbed surface is heated at ≧800 K, the C-O bonds are broken and the adsorbed oxygen atoms are bound only to Mo atoms. Valence PES study shows that the oxygen adsorption induces a peculiar state around the Fermi level, which enhances the emission intensity at the Fermi edge in PES spectra.     

Electron-induced H atom desorption patterns created with a scanning tunneling microscope: Implications for controlled atomic-scale patterning on H-Si(1 0 0)

We have used scanning tunneling microscopy (STM) to explore the details of single and multiple H atom desorption from the H-Si(1 0 0)-2 × 1 surface induced by the inelastic scattering of electrons from an STM tip. The desorption of pairs of H atoms from individual Si dimers is rarely observed. Two-H atom desorption most often involves pairs of dimers, in the same or adjacent rows. This suggests that recombinative H₂ desorption via an interdimer reaction pathway, like that observed recently under nanosecond laser heating, may also be operative for electron-induced excitation using STM. Repeatable fabrication of desired size-selected dangling bond (DB) clusters is also achieved. The single atomic precision of the fabrication is a result of the intrinsically unfavorable paired H atom desorption from a single dimer, but does not result from the spatial localization of excitation energy of the Si-H bond under the STM tip as suggested in previous studies.     

A density functional study of atomic oxygen and water molecule adsorption on Ni(1 1 1) and chromium-substituted Ni(1 1 1) surfaces

Density functional theory (DFT) for generalized gradient approximation calculations has been used to study the adsorption of atomic oxygen and water molecules on Ni(1 1 1) and different kind of Ni-Cr(1 1 1) surfaces. The fcc hollow site is energetically the most favorable for atomic oxygen adsorption and on top site is favorable for water adsorption. The Ni-Cr surface has the highest absorption energy for oxygen at 6.86 eV, followed by the hcp site, whereas the absorption energy is 5.56 eV for the Ni surface. The Ni-O bond distance is 1.85 Å for the Ni surface. On the other hand, the result concerning the Ni-Cr surface implies that the bond distances are 1.93-1.95 Å and 1.75 Å for Ni-O and Cr-O, respectively. The surface adsorption energy for water on top site for two Cr atom substituted Ni-Cr surface is 0.85 eV. Oxygen atoms prefer to bond with Cr rather than Ni atoms. Atomic charge analysis demonstrates that charge transfer increases due to the addition of Cr. Moreover, a local density of states (LDOS) study examines the hybridization occurring between the metal d orbital and the oxygen p orbital; the bonding is mainly ionic, and water bonds weakly in both cases.