SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Review

Self-organised synthesis of Rh nanostructures with tunable chemical reactivity

F Buatier de Mongeot1, A Toma1, A Molle1, S Lizzit2, L Petaccia2 and A Baraldi34*

Author Affiliations

1 Dipartimento di Fisica, Università di Genova and CNISM, Via Dodecaneso 33, 16146, Genova, Italy

2 Sincrotrone Trieste S.C.p.A, S.S. 14 Km 163.5, 34012, Trieste, Italy

3 Physics Department and Center of Excellence for Nanostructured Materials, Trieste University, Via Valerio 2, 34127, Trieste, Italy

4 Laboratorio TASC INFM-CNR, S.S. 14 Km 163.5, 34012, Trieste, Italy

For all author emails, please log on.

Nanoscale Research Letters 2007, 2:251-264  doi:10.1007/s11671-007-9059-3


The electronic version of this article is the complete one and can be found online at:


Received:30 March 2007
Accepted:13 April 2007
Published:22 May 2007

© 2007 to the authors

Abstract

Nonequilibrium periodic nanostructures such as nanoscale ripples, mounds and rhomboidal pyramids formed on Rh(110) are particularly interesting as candidate model systems with enhanced catalytic reactivity, since they are endowed with steep facets running along nonequilibrium low-symmetry directions, exposing a high density of undercoordinated atoms. In this review we report on the formation of these novel nanostructured surfaces, a kinetic process which can be controlled by changing parameters such as temperature, sputtering ion flux and energy. The role of surface morphology with respect to chemical reactivity is investigated by analysing the carbon monoxide dissociation probability on the different nanostructured surfaces.

Keywords:
Nanostructured materials; Nanoscale pattern formation; Rhodium; Surface chemical reactivity; Carbon monoxide

Introduction

The control of the atomic step distribution of clusters and nanostructures is of utmost importance in determining, among others, their magnetic [1], electrical, and catalytic properties. Recent experiments and theoretical models have tried to elucidate the atomistic details underlying the enhanced surface chemical reactivity of these active sites of transition metals (TM). Among these, extensive studies of carbon monoxide chemisorption on TM surfaces have been a valuable resource for the development of surface chemistry. CO, a toxic molecule contained in the automotive exhaust gases, is object of conversion via catalytic oxidation reduction [2]. The carbon monoxide dissociation process is a key step in the syngas reaction which is widely used in the industrial chemistry [3] for methane formation via the reaction or in the Fisher–Tropsch reaction [4] where CO and H2 are transformed in hydrocarbons via the (n > 2) reaction. This, along with the basic interest in understanding the mechanism involved in dissociative adsorption of heteroatomic molecules, has placed CO in the list of the most extensively studied adsorbed molecules [5]. The room temperature interaction of CO with transition metal surfaces can be divided in two main groups. The first includes transition metals from the left side of the periodic table, such as Fe, W and Mo, which adsorb CO dissociatively, while the second is composed by elements from the right side of the periodic table such as Co, Ni, Ru, Rh, Pd, Ir and Pt, which tend to adsorb CO molecularly.

However CO dissociation can also occur on the latter metals, in particular Ru, Rh and Ni, under defined temperature, pressure and surface structural conditions, which allow the molecules to overcome the activation barrier for dissociation. Detailed experimental and theoretical investigations performed in the last 15 years report that the chemical reactivity strongly increases on corrugated surfaces and that CO dissociation is sensitive to the structure of the substrate: steps and kinks drastically modify the reaction paths on solid surfaces and appear to be the most active sites for the C–O bond breaking [6,7].

To this respect a general relation between the chemical reactivity, the d-band center and thus the coordination number of surface atoms, has been established [8-10]: the lower the coordination number of TM surface atoms, the smaller the local bandwidth and the higher the d-band centre position relative to the Fermi level for metals like Rh with a more than half-filled d band. Detailed calculations [11] based on Density Functional Theory (DFT) report reaction barriers for the reaction which strongly decrease when passing from the flat Rh(111) surface (Ea = 1.17 eV) to steps (Ea = 0.30 eV) and kinks (Ea = 0.21 eV). Coordination numbers of these metal atoms range from n = 9, to n = 7 (steps) and n = 6 (kinks). Other DFT calculations by Mavrikakis et al. [12] have shown that the energy of the transition state for CO dissociation on Rh(211) is about 120 kJ/mol lower than on the (111) terrace. In this case the coordination of the atoms at the steps is 7.

Experimental investigations indeed report that CO dissociation is negligible on close packed (111), (110) and (100) Rh surfaces [13] and that it increases on stepped (211) [12], (210) Rh substrates [14]. The dissociation process has been also extensively studied on Rh nanoparticles supported on thin Al2O3 films grown on a NiAl(110) single crystal, as a function of particle size [15,16]. Maximum activity has been measured for particles containing about 1000 atoms, but the nature of the active sites was not explained.

It is a natural consequence of these detailed surface science studies to expect that promotion, enhancement, steering and control of CO dissociation can be reached by simply tuning surface morphology with the purpose of changing the density of reaction centers. Recently, it has been found that it is possible to tune the morphology and step distribution of a crystalline Rh(110) substrate by controlled exposure to a beam of noble gas ions: Xe ion irradiation at few hundreds eV leads to the formation of nonequilibrium periodic nanostructures such as nanoscale ripples, oriented either along [001] or [1-10] directions, mounds, and unexpected rhomboidal pyramids (RP) [17,18]. The latter nanostructures are particularly interesting as candidate model systems for testing catalytic reactivity, since they are endowed with steep facets running along nonequilibrium directions, exposing a high density of undercoordinated atoms.

Morphological characterisation of nanostructures

The irradiation of transition metal surfaces with an energetic noble ion beam yields the self-organization of a great number of nanoscaled patterns. They originate from the surface instability induced by the ion sputtering as well as from the diffusion balance among the removed adatoms.

Here we show how the sputtering parameters, i.e. the substrate temperatureT, the ion fluxf, the impact energy ɛ may influence the structural features of the surface patterns in the case of a Rh(110) surface sputtered with energetic Xe ions at an incidence angle of 15° from the surface normal.

Experimental: Spot Profile Analysis-LEED

Structural characterization of nanoscale surfaces was performed by in situ spot profile analysis low energy electron diffraction (SPA-LEED) which provides information on large surface area by integrating the diffraction signal on the scale of the electron beam size (about 0.1 mm). As a general consideration on the electron diffraction from sample crystal surfaces, constructive (in-phase) and destructive (anti-phase) interference from the exposed terraces are identified respectively by integer and half-integer values of the vertical scattering phase Sz = kzd/2π, kz and d being the vertical momentum transfer and the monoatomic step height [19].

Taking benefit from an instrumental transfer width of about 0.1 μm (i.e. the maximum lateral extension over which the electron wavelength is coherent), the SPA-LEED provides a high resolution magnification on the diffraction spot features. In this way we can investigate the facet distribution of the nano-structured entities and the dominant Fourier modes in the height–height correlation of the surface profile [17]. While a selective faceting can be recognized from the emergence of a couple of Sz dependent Bragg peaks at both sides of the central (0,0) Bragg peak, the presence of a lateral correlation Λ between the surface nanostructures is displayed as an additional, Sz independent, satellite splitting of the near to in-phase spot profile (Sz close to integer values). Figure 1 displays a three dimensional (3D) plot of the (0,0) spot profile for a one dimensional (1D) ripple-like pattern on the parallel momentum transfer plane . The image was recorded in a scattering phase condition such that both couples of satellites can be envisioned (see also the dashed lines in the contour plot of the spot profile reported in the inset). As mentioned above, the inner couple of satellites is due to the average periodicity of the rippled corrugations whereas the outer couple, whose splitting turns out to be linearly dependent on the scattering phase, reflects the diffraction from the slope selected ripple facets. Therefore the out-of-phase spot profiles, i.e. far from integer values of vertical scattering phase Sz, can be effectively regarded as maps of the average facet slopes which have been selected during the surface structuring/re-organization upon ion irradiation.

thumbnailFigure 1. 3D plot of the (0,0) diffraction peak acquired under near-in phase conditions (Sz = 2.1) by means of SPA-LEED for a rippled surface. The multiple satellite splitting indicates the presence of a well defined lateral correlation and facets (see text for details). Inset: 2D contour plot of the spot profile (the dashed lines indicate the position of the satellites). Reproduced from Ref. [18]

The role of the substrate temperature during ion irradiation is reported in Fig. 2. Each spot profile is linked to a peculiar morphology obtained by patterning the initially flat Rh(110) substrate (terrace width above 60 nm) with Xe ion beam with ɛ = 1000 eV at three different temperatures, 450, 500 and 550 K. A ripple corrugation elongated in the direction (Low Temperature Ripple-LTR motif) is inferred from the spot profile in Fig. 2a, measured in a nearly out-of-phase condition (Sz ≈ 2.1), with a lateral periodicity of the ripple nm. The outer satellites are due to the presence of regular step arrays (and facets) along the direction. When moving towards an out-of-phase diffraction point, the outer satellites dominate the diffraction pattern and their splitting increases linearly with Sz giving the facet slope [20]. Increasing the temperature up to 550 K (Fig. 2c) the spot profile apparently presents some qualitative affinities—e.g. twofold symmetry, satellite splitting—with the LTR case in Fig. 2a apart from a rotation of 90° in the reciprocal space which indicates a ripple orientation along the , with facets and dominant steps in the same direction (High Temperature Ripples pattern). The real space morphology corresponding to the HTR pattern is displayed in the surface topography in Fig. 2d taken by Atomic Force Microscopy (AFM). The AFM image elucidates that the High Temperature Ripples (HTR) pattern consists of almost 1D corrugations roughly elongated in the surface direction with a high degree of lateral order (Λ = 60 ± 2 nm) as confirmed by well defined couple of first order correlation peaks resulting in the fast Fourier transform diagram shown in the inset. Therefore the temperature increase from 450 K to 550 K gives rise to a rotation of the ripples by 90°.

thumbnailFigure 2. Out-of-phase Spot Profile (Sz = 2.2) after sputtering at: T = 450 K (a), T = 500 K (b), and T = 550 K (c). The axis scale is expressed in % of Brillouin Zone referred to the short side of the reciprocal surface cell. Panel (b) and (c) are reproduced from Ref. [17]. (d) AFM topography corresponding to an HTR morphology equivalent to panel (c) on a 600 × 600 nm scale. Inset: FFT contour profile of the AFM image. ɛ = 1000 eV

If the intermediate temperature of 500 K is chosen, the out-of-phase diffraction profile in Fig. 2b is modified in a cross-like shape originating from the coexsistence of high symmetry 〈1−10〉 and oriented facets which bound the sides of the rectangular shaped islands constituting the Round Mouth (RM) pattern. The arrangement of the RM pattern actually appears as the interplay between the majority bounding steps in the two extreme cases of the LTR and HTR patterns. The temperature sequence of the surface patterns formed during Rh(110) sputtering is similar to that observed in the case of Ag(110) [21] and Cu(110) [22], apart from a shift to higher temperatures consistent with the larger diffusion barriers for Rh(110). The observation of the RM pattern was also reported in the homoepitaxial growth on Ag(110) [20].

The effect of the impact energy ɛ of the impinging Xe ions on the SPA-LEED pattern is shown in Fig. 3 where the out-of-phase (Sz ≈ 2.2) SPA-LEED plots related to three values of the impact energy, 700, 400 and 220 eV are reported. At ɛ = 700 eV both the two-fold symmetry and the 1D character of the spot profile identifies a LTR motif with = (15.0 ± 0.5) nm and = 10.8° ± 2.5° (Fig. 3a), qualitatively similar to that of Fig. 2a, i.e. with ripple corrugations bound by oriented steps and facets. The LTR pattern undergoes a loss of the lateral correlation between adjacent corrugations and of the structural order in the facet arrangement upon a slight reduction of the impact energy to ɛ = 400 eV as demonstrated by the broadening of the spot profile in Fig. 3b. With ɛ = 220 eV, the (00) diffraction spot splits in fourfold fashion, with satellites along low symmetry directions (see the four-lobe shaped spot profile in Fig. 3c).

thumbnailFigure 3. Out-of-phase (Sz ≈ 2.2) diffraction pattern after irradiation of Rh(110) at different ion energies ɛ. (a) ɛ = 700 eV, (b) ɛ = 400 eV, and (c) ɛ = 220 eV. The other parameters of ion irradiation are T = 450 K, f = 1.5 ML/min, and ion fluence 67.5 ml. The figure is reproduced from Ref. [17]

A quantitative description of the new structure at the atomistic level is reported in Fig. 4. Figure 4a shows a sequence of diffraction patterns of the (00) spot, acquired at increasing Sz values after patterning the Rh(110) surface with Xe ions with ɛ = 220 eV at T = 450 K. Moving away from the condition of constructive interference (Sz = 2), the (0,0) peak dissolves into four equivalent sub-peaks with a mutual separation linearly proportional to the variation of the scattering phase. Therefore the nature of the sub-peaks can be directly ascribed to the formation of facets with non-conventional bounding step-edges. The orientation of the majority steps has been deduced by plotting the satellite position in the parallel momentum transfer plane (), for different values of (see the diagram in Fig. 4b). The majority steps constituting the facets turn out to be placed along the diagonal axis of the surface unit cell, i.e. the directions of low symmetry which form an angle with respect to the axis (see the inset). The average tilt angle of the facets with respect to the (110) reference plane is derived by linear extrapolation of the satellite splitting in a completely destructive interference condition. The method is shown in Fig. 4c where from the plot of the satellite splitting versus the scattering phase Sz a facet slope α = 11° ± 2° and terrace width Γ ≈ 6.6 Å are derived. As sketched in the hard sphere cartoon model of Fig. 4d, the terrace width is compatible with three lattice constant along the close-packed atomic rows. The four-lobe diffraction profile can thus be related to real-space objects bound by four dominant facets pointing towards low symmetry directions, consistently with the presence of rhomboidal pyramids (RP) on the surface. The lateral periodicity Λ ≈ 14.7 nm, deduced from near-in-phase measurements, definitely supports the identification of the four-lobe spot profiles in Fig. 4a with an ordered pattern of RP mounds in the real space. These observations allow us to unambiguously identify the novel morphological state with the RP state theoretically predicted in Ref. [23]. However, since the latter model is expressed in terms of dimensionless parameters, the experiment should elucidate how the experimental variables determine the transition between the various morphological states.

thumbnailFigure 4. Structural characterization of the RP state: (a) sequence of the (00) spot profile at several vertical scattering phases Sz (the arrows indicate the position of one satellite); (b) satellite position in the parallel momentum transfer plane using Sz as a parameter; (c) plot of the satellite splitting as a function of Sz; (d) simplified cartoon of a RP island with step edges running along well-defined directions. Darker colors correspond to topmost layers. The figure is reproduced from Ref. [17]

The role of the impact energy in the surface nanostructuring and the concomitant formation of the RP pattern have been also generalized to the Cu(110) surface [17,18] provided that the substrate temperature was properly rescaled because of the different diffusion barriers for the two transition metal surfaces. It is just the case of the Cu(110) surface erosion which enables us to discriminate two contrasting trends in the energy dependence of the correlation length Λ and of the facet slope α of the relevant surface patterns within the energy range from 200 eV to 3,000 eV (see the diagram in Fig. 5). For 200 eV eV, Λ follows a steep monotonic decrease with the energy, whereas, above 500 eV, this behaviour is reversed. An opposite trend is observed for the slope of the relevant surface structures: first increasing ( eV) and then decreasing ( eV). The observation of well-distinguished behavior for Λ and α in the considered range of impact energies is strictly related to the structural transition from the RP phase to the LTR phase as the critical energy of 500 eV is crossed. This argument is supported by the observation of different facets which characterize the two structural phases, as reported in the energy dependence of the facet slope in Fig. 5. For clarity, the RP pattern develops facet slopes extended along the real space direction, having majority steps along the , as deduced from the analysis in Fig. 4, and linearly increasing with the energy. Conversely, the LTR morphology consists of 1D corrugations elongated in the direction with roof-top facets through the perpendicular direction with decreasing slope for increasing energy.

thumbnailFigure 5. Energy dependence of the correlation length Λ along the (open circles) and of the facet slope (rhomboidal points refer to RP facets, triangular point to LTR facets) of the surface structures. The vertical line marks the border between the RP and the LTR regime (see text for details). The line between dots is a guide to the eye. The figure is reproduced from Ref. [18]

The emergence of the LTR phase, for ɛ > 500 eV, can be rationalized from an atomistic viewpoint based on Scanning Tunneling Microscopy (STM) investigations of single ion impacts on the Ag(001) and Pt(111) surfaces [24,25]. In those studies vacancy clusters with size of several nanometers are generated by the ion collisions and are coupled to several surrounding clusters which consist of the displaced adatoms. Within the crater width every correlation is reasonably suppressed due to the locally hyperthermal collision transient. The ascending character of the wavelength, observed in our data, can be thus regarded as a consequence of the increasing crater radius with the impact energy in agreement with molecular dynamics simulations [26]. The LTR state arises then from the onset of impact-induced “hot spots” involving a local surface melting in the volume around the collision point [25,27]. This assumption is corroborated by the slope relaxation intervening at higher impact energy (see Fig. 5), i.e. when the thermal spike affects wider areas of the surface [28]. In addition, the loss of correlation for higher energy can be reasonably associated either to the stronger excitation transferred to the surface atoms from the impinging ions or to an increase of the lateral extension of the impact crater.

Let’s now consider the decrease of the correlation length Λ of the RP state, when ɛ ranges from 200 eV to 500 eV. Again, from the atomistic approach, this behavior can be understood in terms of the actual damage produced by the ion impact. Contrary to the “thermal spike” picture applied to the LTR case, decreasing the impact energy may have a critical role in the concentration of mobile defects at the topmost surface layer. This argument is supported by the observation of a monotonic increase of the adatom yield with the energy on the Pt(111) surface in the energy range 40–10,000 eV, as follows from STM analysis of single ion impacts [25] and from molecular dynamics simulations [26]. The role of the impact energy is effectively analogous to that of the deposition flux in homoepitaxial growth: increasing the energy yields a higher concentration of mobile defects (mainly adatoms), which rearrange in stable nuclei after the ion impact. The higher density of stable nuclei on the surface and the lower correlation length between them, can explain the behavior of Λ in the RP regime (Fig. 5).

These arguments—especially the tendency of the impact energy to act as the deposition flux—find a phenomenological confirmation in the comparative study shown in Fig. 6, where the surface patterns induced on Cu(110) surface by Xe ion beam are investigated as a function of the substrate temperature T, the ion flux f and impact energy ɛ, as for Rh(110). In Fig. 6a–c the Xe sputtering temperature dependence of the surface morphology is explored at fixed energy ɛ = 400 eV: the panels show a set of out-of-phase spot profiles for the three temperatures of 200, 230 and 290 K after exposing the Cu(110) surface to a Xe ion beam. Sputtering at T = 200 K (Fig. 6a) generates a LTR pattern. At the intermediate temperature T = 230 K (Fig. 6b) the spot profile acquires the fourfold symmetry, indicating that the RP state is formed. Finally, increasing the temperature to 290 K (Fig. 6c), the transition towards the HTR configuration is observed. This sequence confirms that also for Cu(110) the reduction of the energy below 400 eV is mandatory to access the RP pattern. On the other hand, the strong dependence of the pattern on the substrate temperature suggests that the evolution of the RP state is mainly dictated by the kinetics of thermally activated diffusion processes involving the adatoms displaced after the ion collision, rather than by the details of the hyperthermal ion-surface collision transient.

thumbnailFigure 6. Out-of-phase diffraction maps (Sz = 1.88) after ion irradiation of Cu(110). (a), (b), (c) Temperature dependence after dosing at ɛ = 400 eV, f = 1.5 ML/min for a total fluence of 21 ML. (d), (e), (f) Flux dependence after dosing at ɛ = 400 eV and T = 230 K for a total fluence of 21 ML. (g), (h), (i) Energy dependence for irradiation at the flux f = 1 ML/min and at T = 230 K for a total fluence of 45 ML. The figure is reproduced from Ref. [17]

Since the ion fluxfdefines the relaxation time between two subsequent collision events, a morphological variation of the surface structure similar to that observed in the T dependence of Fig. 6a–c, is expected even for different ion flux f. This is shown in Fig. 6d–f: here the out-of-phase diffraction maps of the (00) spot are reported for three different fluxes at ɛ = 400 eV and T = 230 K. Sputtering at a relatively low flux (f = 0.1 ML/min—Fig. 6f) results in a faint HTR pattern; in the intermediate range (f = 0.3–1 ML/min—Fig. 6e) the RP state emerges, while further increasing the flux (f = 3 ML/min—Fig. 6d) the transition to a well-resolved LTR state occurs. From this it can be concluded that, the decrease of the ion flux on the surface morphology corresponds to the increase of the substrate temperature, further confirming the kinetic and diffusive mechanisms underlying the pattern formation.

According to the atomistic approach proposed for the discussion of Fig. 5, a further question which has to be addressed is the role of the ion impact energy in the formation of the surface pattern. In Fig. 6g–i we show the dependence of the surface morphology of the Cu(110) surface on ion energy, for fixed f = 1 ML/min and T = 230 K. As already observed in Fig. 3 and 5, the LTR state transforms into the RP pattern when decreasing the energy from ɛ = 600 eV (Fig. 6g) to ɛ = 400 eV (Fig. 6h), whereas a further decrease of the impact energy down to 200 eV allows to revert the surface morphology into an HTR pattern (Fig. 6i). The sequence of Fig. 6g–i is consistent with the sequence of panels Fig. 6d–f, suggesting that an increase of the ion flux f is equivalent to an increase of the ion energy ɛ. Such behavior can be rationalised if we recall that both STM experiments [25] and molecular dynamics simulations of a single ion impact [26] show a monotonic increase of the adatom production yields when the impact energy is increased in the range 0.1–10 keV. Therefore an increase of the average production rate of adatoms can be achieved either through an increase of ɛ (which affects the number of adatoms produced per impact event) or through an increase of f (which modifies the rate of single ion impact events). This picture is fairly compatible with the atomistic discussion on Fig. 5 according to which the role of the energy can be regarded as that of the deposition flux, i.e. as source of mobile defects.

Furthermore, from the data in Fig. 5 and 6g–i, we can also conclude that the impact energy affects the formation process of the RP state not through a selective anisotropic etching of the surface, but by controlling the total number of adatoms emitted per collision event which constitute the mobile species that enter the destabilizing massive transport at the base of the self-diffusion in fcc(110) terminated metal surfaces.

Theory: the continuum model

The diffusion of the mobile species can be treated in a unified scheme both under erosion as well as growth conditions by considering a non-equilibrium, tilt-dependent flux of defects Jup as responsible for the surface instability which gives rise to the pattern formation. In a continuum approach, the evolution of the surface profile h(x,y,t) obeys the conservation law described in terms of the total adatom current density J = Jsd + Jup(m), Jsd being the curvature dependent surface diffusion term (vanishing in flat regions such as facets) and Jup(m) the destabilising contribution which depends on the local slope vector and biases diffusion uphill, towards ascending step edges [29]; η(x,y,t) accounts for the randomness of the adatom (ion) arrivals.

Metastable structures are reached whenever steps rearrange by selecting local surface slopes m* which minimize the Jup current, i.e. for stable zero solutions of Jup[29,30]. Following Ref. [23,31], the arrangement of each surface pattern is dictated by the vectorial character of the m* solutions. In more detail, non-equivalent pairs of doublet-like solutions (m* = (±m1 , 0) or (0 , ±m2)) identify two kinds of rippled structures, i.e. the LTR and the HTR pattern, depending on whether the faceting (and the periodicity) is extended along the or the direction, respectively [32]. The coexistence of both pairs of doublet-like solutions (m* = (±m1, 0) and (0, ±m2)) expresses the occurrence of rectangular faceting which identifies the RM pattern. Quartet-like solutions of the form m*=(±m1, ±m2) reflect the presence of 2D pyramid structures having a rhomboidal shaped contour lines formed by step-edges (the RP state). Consecutive transitions among these surface configurations can be achieved upon variation of the intrinsic details of the Jup current as results in the morphological phase diagram of Fig. 7a obtained from a continuum approach to the surface evolution under erosion [23,31]. It is also interesting to stress how the same theoretical study allows to single out the formation of surface patterns with structural arrangements similar to those observed in our experiments. For instance, surface models elaborated in Fig. 7b discriminate a universal rippled state from peculiar RP pattern and from a RM pattern. The strong similarity between the nanopatterning of the Rh(110) and the Cu(110) surfaces suggests that the predictions of the continuum approach and the consecutive morphological transformations can be extended to a broader class of unreconstructed fcc(110) substrates. Finally, the observation of the scaling law of the lateral correlation and of the interface width (rms roughness) of the rhomboidal pyramids and of the LTR ripples as a function of the irradiation dose (Fig. 8a and 8b) provides a good agreement with the predictions of the continuum models of Ref. [23,31] further reinforcing the assignments made. Additionally the data shown in Fig. 8 demonstrate that it is possible to tune the separation of the nanopyramids and of the LTR ripples as well as their roughness within significant range, which provides a valuable tool for tuning their morphological properties in view of applications.

thumbnailFigure 7. (a) Kinetic phase diagram depicting various interfacial states in terms of their basic properties (as obtained from simulations of ref. [31]): surface contour plots, magnitudes of interfacial height Fourier transforms (FT), corresponding to near in-phase diffraction patterns, and slope distributions (SD) in the slope space, corresponding to out-of-phase diffraction patterns. R1 and R2 are the two rippled states, RhP is the rhomboidal pyramid state, R1rec and R2rec are the two rectangular rippled states, and R1buc(hatched domain) is the buckled rippled . See Ref. [31] for details. (b) Three-dimensional (3D) views of various interfacial states from the simulations of Ref. [31]: (a) rippled state with the period λ, (b) rhomboidal pyramid state with the periods λ1 and λ2, and (c) rectangular rippled state, with the period λ, which motif is a rooflike pyramid (hut), with a long rooftop edge of the length ξ. The figures are reproduced from Ref. [31]

thumbnailFigure 8. Evolution of the correlation length through the direction (a) and of the estimated surface roughness West (b) for the RP and LTR surface patterns at different fluences. In panel (a) the scaling behavior of the LTR wavelength is also reported. The figure is reproduced from Ref. [18]

Chemical reactivity

We have studied CO interaction with the Rh nanostructured surfaces described in Sect. 2.1 by using High-Energy Resolution Core Level Spectroscopy with synchrotron radiation, to probe the coverage evolution and the molecular dissociation process. The photoemission studies were performed at the SuperESCA beamline [33,34] of the Elettra third generation synchrotron radiation source in Trieste, Italy. The experimental chamber is equipped with a double pass hemispherical electron energy analyser with 96 channels detector [35]. During the measurements the background pressure in the main chamber was always better than 2 × 10−10 mbar. The Rh(110) single crystal was cleaned by Xe ion sputtering at room temperature (E = 1 keV), flash to 1300 K, oxygen cycles in order to remove residual carbon (in the range 570–1070 K at PO2 = 5 × 10−8 mbar) and finally, hydrogen reduction to remove residual oxygen traces (PH2 = 1 × 10−7 mbar, T = 470–770 K). Surface cleanliness prior to nanostructures preparation was checked by measuring C1s, S2p and O1s signals. C1s and O1s spectra were recorded always at a sample temperature of 200 K in order to reduce temperature broadening of the peaks and in normal emission conditions. Photon energies of 400 and 650 eV were used for C1s and O1s spectra, with an overall energy resolution (X-ray monochromator and electron energy analyser) of 150 and 300 meV, respectively. In these conditions typical data acquisition time was 5 min/spectrum. Core level spectra binding energies have always been calibrated with respect to the Fermi level.

The XPS analysis was done by fitting the core level spectra with a Doniach-Šunjić (DS) function [36], characterized by two parameters: the singularity index α (describing the asymmetry of the core level spectra due to electron–hole pairs excitations) and the Lorentzian width Γ (because of the natural core-hole lifetime), convoluted with a Gaussian, which takes into account the broadening due to unresolved vibrations, many-body effects and the instrumental resolution. A linear background was also included in the fit.

Chemical reactivity: experimental results on RP

It is well established that both oxygen and carbon 1s core-level signals are strongly sensitive to the local molecular and atomic adsorption sites, and can be used to determine the CO adsorption geometry. In particular, for a large number of carbon monoxide adsorption systems it was found that the binding energy (BE) decreases with increasing CO coordination to the substrate atoms, i.e. in the order BE(on-top) > BE(bridge) > BE(hollow), with a shift which is about twice as large for O1s than for C1s [37]. The O1s BE on different TM single-crystal surfaces varies in the range 531.6–532.6 eV for on-top bonded CO and between 530.5 and 531.6 eV for bridge-bonded CO [38]. The reason of this trend can be understood from total energy considerations, the major contribution to the shift originating from the changes in the energy of the core ionized final state. Indeed the difference of the CO adsorption energies between different adsorption sites for the neutral initial state is very small (∼100 meV) [37]. Carbon and oxygen atomic species when chemisorbed on transition metal surfaces, usually produce core level components at lower BE. In particular carbon species are found at about 284 eV, while chemisorbed oxygen at about 530 eV.

The CO adsorption and the temperature evolution of the chemisorbed layer have been measured for HTR, LTR and RP nanostructures, produced using the procedures reported above.

The growth of the RP nanostructures was characterized in situ by low energy electron diffraction (LEED). The appearance of a fourfold splitting of the (00) diffraction peak along diagonal directions demonstrates the formation of the RP facets. The diffraction pattern is in agreement with the SPA-LEED results reported in Fig. 4. Carbon monoxide was firstly dosed on the Rh(110) nanostructured surface at T = 200 K, i.e. well before the CO desorption onset on the clean (1 × 1) Rh(110) surface [39-41], at different initial coverage, ranging from ∼0.03 ML for the RP to saturation.

As for the (111), (100) and (110) Rh surfaces, CO adsorption on RPs at saturation (Fig. 9) leads to the occupation of two different adsorption sites. By analogy with the (110) flat surface [41], the higher BE C1s (Fig. 9a) and O1s (Fig. 9b) components, at 286 and 531.7 eV respectively, are assigned to molecules sitting in on-top sites, while the lower BE peaks, to CO placed in bridge sites (BE of 285.55 and 530.5 eV). The spectra reported in Fig. 9 show that, for both signals, the Gaussian width of the on-top component is always larger than that of the bridge-bonded CO. Besides the presence of unresolved vibrational excitations [40,42], this behavior can be explained in terms of occupation of local inequivalent on-top configurations. Indeed the small terraces present on the RPs can be occupied by CO sitting just below or above the steps (TA and TC sites in Fig. 9c) or in the terrace (TB site). At the contrary the Gaussian width of C1s and O1s signals relative to bridge-bonded CO is much lower and this can be tentatively interpreted as due to single site occupation.

thumbnailFigure 9. (a) C1s and (b) O1s core-level spectra showing the components corresponding to molecular CO adsorbed at saturation on the Rh nano-pyramids. (c) Possibile inequivalent on-top (TA,TB,TC) and bridge (BA, BB) adsorption sites are indicated

Figures 10 and 11 show the evolution of the C1s and O1s spectra after annealing of the CO saturated layer, which evidence the different behavior of the two CO bonding configurations, as well as the onset of CO dissociation. The two CO related components behave similarly in the C1s and O1s regions. The C1s spectra drastically change with increasing the annealing temperature: they lose intensity due to CO desorption while the ratio on-top to bridge population clearly increases with increasing the temperature. This behavior parallels that observed on the flat (110) surface using Temperature Programmed XPS where a conversion from bridge to on-top sites takes place during heating [41].

thumbnailFigure 10. Evolution of the C1s core level spectra after annealing of the CO saturated layer prepared on the Rhomboidal Pyramids at different temperatures and quenching at 250 K. The peak at 286 and 285.55 eV correspond to CO in on-top (orange) and bridge (green) sites, respectively. The component at 283.55 eV (blue) is due to atomic carbon. hν = 400 eV

thumbnailFigure 11. Evolution of the O1s core level spectra after annealing of the CO saturated layer prepared on the Rhomboidal Pyramids at different temperatures and quenching at 250 K. The peak at 531.7 and 530.5 eV correspond to CO in on-top (orange) and bridge (green) sites, respectively. The component at ∼530 eV (blue) is due to atomic oxygen. hν = 650 eV

The evolution of the intensity of the C1s components is reported in Fig. 12. The relative intensities have been rescaled because of photoelectron diffraction effects which are known to be relevant at this photoelectron kinetic energies. In order to do this, high energy spectra (Ekin > 250 eV) have been acquired. This is the reason also of the different on-top to bridge bonded CO population observed in the O1s experiments. Up to about 450 K the decrease of the total CO coverage is mainly dominated by the CO bridge-bonded depopulation. The changes in the relative population of the on-top and bridge sites is in agreement with the predominant occupation of on-top sites at CO coverage > 0.3 ML. The shift to lower binding energy of the on-top component with increasing the annealing temperature (ΔE = −80 meV) can be explained in different ways: (i) preferential occupation of one of the three available on-top adsorption sites (see Fig. 9c), (ii) changes of interatomic molecular interactions between CO molecules at different coverage or (iii) partial occupation of three-fold sites in the (111) facets of the nano-pyramids.

thumbnailFigure 12. Plot of the normalized C1s intensity versus annealing temperature for a CO saturated layer prepared on the Rhomboidal Pyramids at T = 200 K relative to the overall CO coverage (black curve) and to CO species in on-top (orange) and bridge (green) sites. Atomic carbon is shown in blue

However, the most interesting result of the heating process is that not all the CO desorbs but a minor fraction converts into atomic species already at ∼450 K, as evidenced by the increase of lower binding energy components in both, C1s (283.55 eV) and O1s spectra (∼530 eV). Indeed, after the removal of bridge-bonded CO, the on-top sites are gradually depopulated and the surface remains completely free of CO for T > 525 K. After heating to 563 K, 9.4 ± 0.5% of the initial CO has converted into atomic carbon. O1s spectra show the presence of a residual amount of atomic oxygen. Atomic oxygen species are expected to desorb as molecular oxygen at temperatures higher than 750 K. The lower amount of atomic oxygen is therefore interpreted as due to the CO + O→CO2 reaction followed by CO2 desorption [43].

The heating experiment has been repeated with a lower CO coverage (0.21 ML) adsorbed on the RP at 200 K (Fig. 13). In this case the initial layer is formed mainly by CO molecules sitting in on-top sites with a residual population of bridge-bonded molecules (see Fig. 13a top spectrum). We observe a similar temperature behavior of the surface population, with the bridge-bonded CO molecules disappearing before the on-top moieties and the atomic carbon species increasing already at temperatures above 370 K. We interpret the increase of carbon as due to CO dissociation in on-top sites. In this case the fraction of dissociated CO after annealing at 563 K represent 22 ± 3% of the total CO initial intensity, i.e. much larger than the value obtained for the CO saturated layer.

thumbnailFigure 13. (a) Evolution of the C1s core level spectra after annealing at different temperatures of the 0.21 ML CO layer prepared on the Rhomboidal Pyramids and quenching at 250 K. (b) Same as in (a) but with a very small initial CO coverage (Θ = 0.03 ML) hν = 400 eV

The increase of CO dissociation probability on RP with decreasing the initial CO coverage is unambiguously confirmed by the temperature behavior of a very low CO coverage layer (0.03 ML) adsorbed at 200 K (Fig. 13b). In this case 80 ± 14% of the CO molecules adsorbed in on-top sites undergo dissociation.

Chemical reactivity: experimental results on HT and LT Ripples

The large set of inequivalent Rh nanostructured surfaces prepared under different experimental conditions, allowed us to investigate the CO dissociation process also on HTR structures.

Figure 14 shows the temperature evolution of the C1s spectra after CO adsorption at 200 K on the HTR structures and Fig. 15 shows the behavior of the intensity of the different C1s components. As for the RP case, the low temperature spectra can be analyzed by using two fitting components, corresponding to CO in on-top and bridge sites. The saturation coverage is very similar to the RP case, being only 8% larger. For the HTR the normalized top/bridge intensity ratio is slightly smaller, indicating a higher tendency to form CO bridge bonds with the Rh substrate. Up to 350 K the overall CO coverage is almost constant, as is the ratio of on-top to bridge CO populations. CO desorption starts at T > 350 K, while atomic carbon formation takes place only at 475 K. The onset of CO dissociation on the HTR is therefore shifted by about 85 K to higher temperature, compared to CO on the RP nanostructures (T = 375 K). Our results show that for the HTR only 4.2% of the initial CO population converts into atomic carbon, while the rest of CO molecules desorb. The same experiment performed on the LTR nanostructures shows a similar behavior, with the onset of CO dissociation shifted at even higher temperatures with only 3.8% of the CO saturated layer undergoing dissociation.

thumbnailFigure 14. Evolution of the C1s core level spectra after annealing at different temperatures of the CO saturated layer prepared on the High Temperature Ripples and quenching at 250 K. The peak at 286 and 285.55 eV correspond to CO in on-top (orange) and bridge (green) sites, respectively. The component at 283.55 eV (blue) is due to atomic carbon. hν = 400 eV

thumbnailFigure 15. Plot of the normalized C1s intensity versus annealing temperature for a CO saturated layer prepared on the High Temperature Ripples at T = 200 K relative to the overall CO coverage (black curve) and to CO species in on-top (orange) and bridge (green) sites. Atomic carbon is shown in blue

The same annealing experiments have been performed after CO saturation of the (1 × 1) Rh(110) in order to make a detailed comparison with the nanostructured surfaces. We find evidence that in this condition only 2.4% of CO undergoes dissociation. This result can be explained as due to a residual density of surface defects on the (1 × 1) substrate. Indeed, even on a well prepared surface with a typical miscut of 0.2°, about 1% of the surface atoms reside at the steps.

Chemical reactivity: discussion

Among the factors that govern the Rh surface reactivity, the density of under-coordinated atoms involved in the reaction plays a very important role. In order to compare the CO dissociation behavior on Rh nanostructured surfaces with that on other Rh crystallographic planes, we have to consider the different atomic coordination on Rh nanostructures. At variance with the (111), (100) and (110) flat Rh surfaces, which present first layer atoms with coordination 9, 8 and 7 respectively, the nanopyramids are composed by a large number of under-coordinated atoms, as shown in Fig. 16a.

thumbnailFigure 16. Model for Rhomboidal Pyramids (a) determined from the SPA-LEED data and for the High Temperature Ripples (b) and Low Temperature Ripples (c). Atomic coordination nl are also reported, with n = number of nearest neighbors andl = number of next-nearest neighbors

Aiming at a comparison in surface reactivity between different nanostructured surfaces we firstly evaluate the density of under-coordinated surface atoms on RP. Figure 16a shows the local model structure of the top of a RP: it is evident that, together with a small amount of Rh atoms with coordination n = 5 (n = number of nearest neighbors) (yellow atoms in Fig. 16), there is a large number of atoms with coordination 6 (orange). These are the atoms either forming the steps oriented along the [1–12] direction or sitting at the top of the nano-pyramid. Also the LTR and the HTR nanostructures, reported in Fig. 16b and c are formed by a considerable amount of 6-fold coordinated atoms. However the main difference among these structures resides in the next-nearest (l) neighbors. In Table 1 we have estimated the density of local atomic configurations for the different morphologies that we have produced. The typical average dimension of the RP clusters are 15 nm along the [1-10] direction, corresponding to an average number of 12 exposed layers. The data have been reported according to the number nl of nearest (n) and next-nearest (l) neighbors. What appears peculiar to the RP, with respect to the LTR and HTR morphology, is the very high density of kinked microfacets, produced by the steps running along [1-12], which provide a high density of low-coordinated atoms labeled (6) in Fig. 16. The local atomic configuration of the RP microfacets represents an open morphology [the nearest neighbor number n is 6 with only 1 next n.n. (l = 1)], with an increase of ≈73% of the distance between Rh atoms along the steps compared to the nearest neighbor distance of bulk Rh. Indeed, with decreasing the initial CO coverage the availability of 61 Rh atoms, with CO sitting in on-top sites (estimated surface densities between 3% and 34%, respectively), increases, thus resulting in a higher fraction of dissociated CO. We expect that the activation energy for CO dissociation on these sites can be even lower than on the Rh(211) steps (n = 7) or kinks (n = 6) for which calculations [11] have found activation energies of 0.30 and 0.21 eV, respectively.

Table 1. Occupacy (in %) of the different sites labelednl, according to the number of nearestnand next-nearest neighborslfor the different nanostructures (RP, HTR and LTR) and for the (110) flat terraces. The sites can be identified in the following way: 51 = upper step corner atom, 61 = upper step (100) facet, 71 = upper step (111) facet, 72 = (110) terrace atoms, 81 = lower step RP facets, 92 = lower step (111) facets, 92 = lower step (100) facets

Another interesting feature of our experimental results is the increased CO dissociation rate with smaller initial CO coverage on the RP. This can be correlated with the increased availability of undercoordinated atoms, with CO occupying on-top sites at the steps, i.e. the most probable adsorption site. The on-top CO adsorption geometries on the small (110) microterraces and steps can be very similar, but we expect that, as for the (211) stepped surface [12], CO is more stable at the steps. We associate this with the fact that the Rh 4d electrons are higher in energy and thus interact more strongly with the CO valence states as for CO on Rh stepped surfaces. With increasing the CO coverage only a residual amount of CO molecules can stay on the atoms at the steps. Moreover, the fact that also the adsorbed C and O strongly bind to the steps indicates that during the annealing experiment the products can poison the active sites and we expect that the reaction should stop when all step sites are occupied.

The increased chemical reactivity on nanoparticle catalysts supported on surface oxide is also connected to the presence of special defect configurations. This has been observed in core-level photoemission experiments on alumina-supported Rh particles as a function of particle size [15,16]. The maximum dissociation rate of 50% was found for particles containing about 1000 Rh atoms, with the reactivity decreasing for larger clusters. The RP clusters produced in our experiments, formed by ∼3000 atoms, result in a dissociation fraction ranging from 10% up to 80%, which nicely compares with the 30% of Rh nanoclusters of similar dimension. These findings suggest that the role of the support in CO dissociation is marginal and that chemical reactivity is mainly governed by the density of undercoordinated surface atoms.

Conclusions

The growth of out-of-equilibrium morphological structures formed during Xenon ion irradiation on Rh(110) substrate is reported. The formation of rhomboydal pyramids and the transition to high temperature and low temperture ripple phases result from the interplay between different kinetic processes depending on temperature, flux and impact energy. Our results illustrate the close correlation between dissociation probability and surface morphology of nanostructured surfaces, a fundamental issue in surface chemistry and nano-catalysis. The RP strongly increase the CO dissociation probability because of the high density of under-coordinated Rh atoms. Detailed comparison of dissociation probability of different nanostructured surfaces shows that CO dissociation takes place in on-top configuration at the kinked step edges. This study shows that it is possibile to change in a controlled way the catalytic properties by changing the density of low coordinated atoms on nanostructured surfaces.

Acknowledgements

We thank U. Valbusa for useful discussion and suggestions. Financial support from the MIUR under the Programs PRIN 2003, FIRB 2001, from Fondazione CARIGE and from MAE under bilateral program Italia-Slovenia is acknowledged.

References

  1. Moroni R, Sekiba D, Buatier de Mongeot F, Gonella G, Boragno C, Mattera L, Valbusa U:

    Phys. Rev. Lett.. 2003, 91:167207.

    COI number [1:STN:280:DC%2BD3srjsVeksw%3D%3D]

    Publisher Full Text OpenURL

  2. Crucq A. (Ed): Catalysis and Automotive Pollution Control II, Studies in Surface Science and Catalysis. Elsevier, Amsterdam; 1991. OpenURL

  3. Newsome DS:

    Cat. Rev.-Sci. Eng. 1980, 21:275.

    COI number [1:CAS:528:DyaL3cXitVShuro%3D]

    Publisher Full Text OpenURL

  4. Zhang Y, Jacobs G, Sparks DE, Dry ME, Davis BH:

    Catal. Today. 2002, 71:411.

    COI number [1:CAS:528:DC%2BD38XhtFGqsLw%3D]

    Publisher Full Text OpenURL

  5. Somorjai GA: Introduction to Surface Chemistry and Catalysis. Wiley, New York; 1994. OpenURL

  6. Zubkov T, Morgan GA, Yates JT, Kuhlert O, Lisowski M, Schillinger R, Fick D, Jansch HJ:

    Surf. Sci.. 2003, 526:57.

    COI number [1:CAS:528:DC%2BD3sXhtVagsb0%3D]

    Publisher Full Text OpenURL

  7. Zubkov T, Morgan GA, Yates JT:

    Chem. Phys. Lett.. 2002, 362:181.

    COI number [1:CAS:528:DC%2BD38XmtlWgsr0%3D]

    Publisher Full Text OpenURL

  8. Hammer B, Nørskov JK:

    Surf. Sci.. 1995, 343:211.

    COI number [1:CAS:528:DyaK28XnsFGq]

    Publisher Full Text OpenURL

  9. Hammer B, Nielsen OH, Nørskov JK:

    Catal. Lett.. 1997, 46:31.

    COI number [1:CAS:528:DyaK2sXksFyrsb0%3D]

    Publisher Full Text OpenURL

  10. Dahl S, Logadottier A, Egeberg RC, Larsen JH, Chorkendorff I, Tornqvist E, Nørskov JK:

    Phys. Rev. Lett.. 1999, 83:1814. Publisher Full Text OpenURL

  11. Liu Z-P, Hu P:

    J. Am. Chem. Soc.. 2003, 125:1958.

    COI number [1:CAS:528:DC%2BD3sXkvFKmtQ%3D%3D]

    Publisher Full Text OpenURL

  12. Mavrikakis M, Bäumer M, Freund H-J, Nørskov JK:

    Catal. Lett.. 2002, 81:153.

    COI number [1:CAS:528:DC%2BD38XlsVerur0%3D]

    Publisher Full Text OpenURL

  13. Campuzano JC: The Chemical Physics of Solid Surfaces and Heterogeneous Catalysis.. Edited by King D.A., Woodruff D.P.. Elsevier, New York; 1990:389. OpenURL

  14. Rebholz M, Prins R, Kruse N:

    Surf. Sci.. 1991, 259:L797.

    COI number [1:CAS:528:DyaK38Xks12jtQ%3D%3D]

    Publisher Full Text OpenURL

  15. Andersson S, Frank M, Sandell A, Giertz A, Brena B, Bruhwiler PA, Martensson N, Libuda J, Baumer M, Freund H-J:

    J. Chem. Phys.. 1998, 108:2967.

    COI number [1:CAS:528:DyaK1cXps1KnsQ%3D%3D]

    Publisher Full Text OpenURL

  16. Ertl G, Freund H-J:

    Phys. Today. 1999, 52:32.

    COI number [1:CAS:528:DyaK1MXos1aktw%3D%3D]

    OpenURL

  17. Molle A, Buatierde Mongeot F, Molinari A, Fuerkaiti F, Boragno C, Valbusa U:

    Phys. Rev. Lett.. 2004, 93:256103.

    COI number [1:STN:280:DC%2BD2M%2FnsFaqsQ%3D%3D]

    Publisher Full Text OpenURL

  18. Molle A, Buatier de Mongeot F, Molinari A, Boragno C, Valbusa U:

    Phys. Rev. B. 2006, 73:155418. Publisher Full Text OpenURL

  19. Rusponi S, Costantini G, Buatierde Mongeot F, Boragno C, Valbusa U:

    Appl. Phys. Lett.. 1999, 75:3318.

    COI number [1:CAS:528:DyaK1MXnsVakt74%3D]

    Publisher Full Text OpenURL

  20. Buatierde Mongeot F, Costantini G, Boragno C, Valbusa U:

    Phys. Rev. Lett.. 2000, 84:2445.

    COI number [1:CAS:528:DC%2BD3cXhs1Cms74%3D]

    Publisher Full Text OpenURL

  21. Rusponi S, Boragno C, Valbusa U:

    Phys. Rev. Lett.. 1997, 78:2795.

    COI number [1:CAS:528:DyaK2sXisVWjtL0%3D]

    Publisher Full Text OpenURL

  22. Golubovic’ L, Levandovsky A, Moldovan D:

    Phys. Rev. Lett.. 2002, 89:266104. Publisher Full Text OpenURL

  23. Costantini G, Buatierde Mongeot F, Boragno C, Valbusa U:

    Phys. Rev. Lett.. 2000, 86:838. Publisher Full Text OpenURL

  24. Michely Th, Teichert C:

    Phys. Rev. B. 1994, 50:11156.

    COI number [1:CAS:528:DyaK2MXhslSgur4%3D]

    Publisher Full Text OpenURL

  25. Bringa EM, Nordlund K, Keinonen J:

    Phys. Rev. B. 2001, 64:235426. Publisher Full Text OpenURL

  26. van Dijken S, de Bruin D, Poelsema B:

    Phys. Rev. Lett.. 2000, 86:4608. Publisher Full Text OpenURL

  27. Ghaly M, Averback RS:

    Phys. Rev. Lett.. 1994, 72:364.

    COI number [1:CAS:528:DyaK2cXhslCitLg%3D]

    Publisher Full Text OpenURL

  28. Villain J:

    J. Phys. (France). 1991, 1:19. Publisher Full Text OpenURL

  29. Levandovsky A, Golubović L, Moldovan D:

    Phys. Rev. E. 2007, 74:061601. Publisher Full Text OpenURL

  30. Boragno C, Buatier de Mongeot F, Costantini G, Molle A, De Sanctis D, Valbusa U, Felici R, Ferrer S:

    Phys. Rev. B. 2003, 68:094102. Publisher Full Text OpenURL

  31. Baraldi A, Barnaba M, Brena B, Cocco D, Comelli G, Lizzit S, Paolucci G, Rosei R:

    Electron Spectrosc. Relat. Phenom.. 1995, 76:145.

    COI number [1:CAS:528:DyaK28XisVSmsrc%3D]

    Publisher Full Text OpenURL

  32. Baraldi A, Comelli G, Lizzit S, Kiskinova M, Paolucci G:

    Surf. Sci. Rep.. 2003, 49:169.

    COI number [1:CAS:528:DC%2BD3sXivVKgurk%3D]

    Publisher Full Text OpenURL

  33. Baraldi A, Dhanak VR:

    J. Electron Spectrosc. Relat. Phenom.. 1994, 67:211.

    COI number [1:CAS:528:DyaK2cXjtVOrsbY%3D]

    Publisher Full Text OpenURL

  34. Doniach S, Šunjić M:

    J. Phys. C. 1970, 3:185. Publisher Full Text OpenURL

  35. Mårtensson N, Nilsson A: Springer Series in Surface Science. Edited by Eberhardt W.. Springer-Verlag, Berlin; 1994. OpenURL

  36. Dhanak VR, Baraldi A, Comelli G, Paolucci G, Kiskinova M, Rosei R:

    Surf. Sci.. 1993, 295:287.

    COI number [1:CAS:528:DyaK2cXhsFygtg%3D%3D]

    Publisher Full Text OpenURL

  37. Baraldi A, Dhanak VR, Comelli G, Prince KC, Rosei R:

    Surf. Sci.. 1993, 293:246.

    COI number [1:CAS:528:DyaK3sXls1yhurs%3D]

    Publisher Full Text OpenURL

  38. Dhanak VR, Baraldi A, Comelli G, Paolucci G, Kiskinova M, Rosei R:

    Surf. Sci.. 1993, 295:287.

    COI number [1:CAS:528:DyaK2cXhsFygtg%3D%3D]

    Publisher Full Text OpenURL

  39. Baraldi A, Comelli G, Lizzit S, Cocco D, Paolucci G, Rosei R:

    Surf. Sci.. 1996, 367:L67.

    COI number [1:CAS:528:DyaK28XnsVyrsrY%3D]

    Publisher Full Text OpenURL

  40. Smedh M, Beutler A, Ramsvik T, Nyholm R, Borg M, Andersen JN, Duschek R, Sock M, Netzer FP, Ramsey MG:

    Surf. Sci.. 2001, 491:99.

    COI number [1:CAS:528:DC%2BD3MXntVyktb8%3D]

    Publisher Full Text OpenURL

  41. Baraldi A, Lizzit S, Cocco D, Comelli G, Paolucci G, Rosei R, Kiskinova M:

    Surf. Sci.. 1997, 385:376.

    COI number [1:CAS:528:DyaK2sXmvVGru78%3D]

    Publisher Full Text OpenURL