Density functional theory investigation of CN on Cu(1 1 1), Ni(1 1 1) and Ni(1 0 0)

The adsorption of CN on Cu(1 1 1), Ni(1 1 1) and Ni(1 0 0) has been investigated using density functional theory (DFT). While experimental studies of CN on Cu(1 1 1) show the molecular axis to be essentially parallel to the surface, the normally-preferred DFT approach using the generalised gradient approximation (GGA) yields a lowest energy configuration with the C-N axis perpendicular to the surface, although calculations using the local density approximation (LDA) do indicate that the experimental geometry is energetically favoured. The same conclusions are found for CN on Ni(1 1 1); on both surfaces bonding through the N atom is always unfavourable, in contrast to some earlier published results of ab initio calculations for Ni(1 1 1)/CN and Ni(1 0 0)/CN. The different predictions of the GGA and LDA approaches may lie in subtly different relative energies of the CN 5σ and 1π orbitals, a situation somewhat similar to that for CO adsorbed on Pt(1 1 1) which has proved challenging for DFT calculations. On Ni(1 0 0) GGA calculations favour a lying-down species in a hollow site in a geometry rather similar to that found experimentally and in GGA calculations for CN on Ni(1 1 0).     

LCAO study of the Cu(1 1 0)-p(2 × 1)O surface

We have performed semi-empirical LCAO calculations of the electronic structure of the Cu(1 1 0)-p(2 × 1)O surface. This has been done accounting for the Cu-Cu interactions by means of a recently proposed set of parameters, which give very good results for the bulk as well as for the surfaces of lowest Miller indices. Furthermore, the O-O interactions, which have been neglected in the preceding similar studies, have been taken into account. The resulting surface bands are in very good agreement with the overall set of the available experimental data. Several issues concerning the physical properties of this surface are addressed in the present paper: the changes induced on the clean surface bands by the adsorption and the reconstruction; the arrangement of the Cu and O atoms in the added rows; the position of the p_y antibonding band of the oxygen. In particular, we have found that the latter has an energy of −0.2 eV at the point. This result confirms an experimental indication in the same direction previously reported by Courths et al. [R. Courths, S. Hüfner, P. Kemkes, G. Wiesen, Surf. Sci. 376 (1997) 43].     

Density functional study of ethylamine and allylamine on Si(1 0 0)-2 × 1 and Ge(1 0 0)-2 × 1 surfaces

We have investigated the interactions of ethylamine and allylamine with models of the Si(1 0 0)-2 × 1 and Ge(1 0 0)-2 × 1 semiconductor surfaces. Ab initio molecular orbital calculations, along with density functional theory (DFT), are used to examine the interaction of these amines with cluster models of the semiconductor surfaces. The transition states and final adsorption products for adsorption of the molecules are predicted. The DFT calculations show the amines form N-dative bond states with Si(1 0 0)-2 × 1 or Ge(1 0 0)-2 × 1 as the initial adsorption product. The initial dative-bond products can be further activated, resulting in N-H bond cleavage on both surfaces. The overall reaction of a given amine on Si(1 0 0) via N-H dissociation is more exothermic than on the Ge(1 0 0) surface.     

Molecular assembly of tetracene on the (2 × 1)O reconstructed Cu(1 1 0) surface

Using a combination of scanning tunneling microscopy (STM) and density functional theory calculations, we have studied the adsorption of tetracene on the Cu(1 1 0) (2 × 1)O substrate. At monolayer coverage the adsorbed molecules are in the flat-laying geometry with their long axis along the close-packed [0 0 1] direction of the substrate and a long-range ordered structure on the length scale up to 100 nm has been observed. DFT calculation results indicate a stronger interaction between tetracene molecules and Cu(1 1 0) substrate than Cu(1 1 0) (2 × 1)O substrate. The preferential adsorption sites have also been pointed out on both substrates. The observed wavelike structure is explained by the interdigitation of C-H bonds of adjacent molecules.     

Theoretical STM images of alkaline-earth metal adsorbed Si(1 1 1)3 × 2 surfaces

We have performed total-energy calculations to study theoretical scanning tunneling microscopy (STM) images of the Si(1 1 1)3 × 2 surfaces induced by the adsorption of alkaline-earth metals (AEMs). Previously, in a series of works on Ba/Si(1 1 1) system, we have found that the observed Si(1 1 1)3 × 1-Ba LEED phase indeed has a 3 × 2 periodicity with a Ba coverage of 1/6 ML and the HCC substrate structure. Based on results of the Ba case, we proposed that the HCC structure is also adopted for other AEM atoms, which was confirmed by our recent work. In this paper, we mainly report the STM simulations for different AEM systems to compare with existing experimental data. We discuss the difference in the detailed STM images for different AEM adsorbates. Especially, the difference in filled-state images between Mg and other AEM atoms is attributed to the strong Mg-Si interaction.     

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.     

Structure of the SiC (0 0 0 1) 3 × 3 reconstruction studied by surface X-ray diffraction

The surface structure of the 3 × 3 reconstruction of the 4H-SiC (0 0 0 1) surface was investigated with surface X-ray diffraction (SXRD).Of the studied models, the twist model proposed by Starke et al. [U. Starke, J. Schardt, J. Bernhardt, M. Franke, K. Reuter, H. Wedler, K. Heinz, J. Furthmuller, P. Kackell, F. Bechstedt, Phys. Rev. Lett. 80 (1998) 758] gave the best fit to the experimental data. The structural parameters were determined accurately.     

Reconstruction and de-reconstruction of the Ir(1 0 0) surface and ultrathin Fe/Ir(1 0 0) films

The structure, energetics and magnetic properties of the quasihexagonal reconstruction of the Ir(1 0 0) surface and nanostructures formed by Fe atoms on this surface have been investigated using first-principles density functional theory with generalized gradient corrections. We find the reconstructed (1 × 5) surface to be 0.10 eV/(1 × 1) area lower in energy than the unreconstructed surface and we demonstrate that first-principles calculations can achieve quantitative agreement with experiment even for such long-period and deep-going reconstructions. For Fe coverage of 0.4 monolayers (ML) we have studied the stripe-like structure with biatomic Fe rows placed in the troughs of the (1 × 5)-reconstructed surface. Results of nonmagnetic calculations agree well with the structure inferred from STM data. Higher Fe coverages lead to a de-reconstruction of the Ir substrate. At 0.8 ML coverage a surface compound with composition Fe₄Ir is formed, which shows an appreciable buckling. In this case, a ferromagnetic calculation leads to good agreement with the low-temperature LEED data. We predict that the (1 × 5) periodicity of the mixed interface layer will persist also in thicker films with a pure Fe surface. Films with 1-4 ML Fe are predicted to be tetragonally distorted and ferromagnetic, with an axial ratio corresponding well to an elastic distortion of the Fe lattice.     

Me-PTCDI thin film deposition on InSb(1 1 1)A

The evolution of N,N′-dimethylperylene-3,4,9,10-dicarboxyimide (Me-PTCDI) thin films formed by vapour deposition on InSb(1 1 1)A substrates has been studied by X-ray photoelectron spectroscopy (XPS), high-resolution electron energy loss spectroscopy (HREELS) and low energy electron diffraction (LEED). XPS studies of the Me-PTCDI covered surface indicate that no significant interaction occurs at sub-monolayer coverage when compared to multilayer Me-PTCDI films. HREELS studies suggest only a weak interaction as evidenced by very small changes in the frequencies of several molecular vibrational modes. LEED patterns show the Me-PTCDI overlayer adopts a structure commensurate with the underlying InSb(1 1 1)A substrate surface and that can be rationalised by van der Waals intermolecular energy calculations for the Me-PTCDI unit cell. The results are consistent with a weak interaction at the Me-PTCDI/InSb interface, the formation of the commensurate structure being sufficient to overcome the small energetic penalty associated with deviation from the calculated intermolecular interaction energy minimum.     

On-surface and sub-surface oxygen on ideal and reconstructed Cu(1 0 0)

In order to understand the first steps of the Cu(1 0 0) oxidation we performed first principles calculations for on-surface and sub-surface oxygen on this surface. According to our calculations, the adsorption energies for all on-surface site oxygen atoms increase, whereas the energies of the sub-surface atoms decrease with the increasing oxygen coverage. At coverage 1 ML and higher on the reconstructed surface, structures including both on- and sub-surface atoms are energetically more favourable than structures consisting only of on-surface adsorbates. On the ideal (1 0 0) surface this change can be perceived at coverage 0.75 ML.     

Vibrational entropy-driven dealloying of Mo(1 0 0) and W(1 0 0) surface alloys

The formation and stability of Cu, Ag and Au-induced c(2 × 2) alloys at the Mo(1 0 0) and W(1 0 0) surfaces have been investigated with low-energy electron microscopy and diffraction. The ordered alloys transform to disordered overlayer structures at elevated temperature. Comparison of the transformation temperatures with energetics obtained from first principles calculations reveals the vibrational entropic contribution to the system free energy that defines alloy thermal stability. Effective Debye temperatures for metal adatoms are determined that exhibit the expected mass and bond strength dependence.     

Surface and subsurface oxide formation on Ni(1 0 0) and Ni(1 1 1)

The oxidation of Ni(1 0 0) and Ni(1 1 1) at elevated temperatures and large oxygen exposures, typical of the methods used in the preparation of NiO(1 0 0) films for surface studies, has been investigated by medium energy ion scattering (MEIS) using 100 keV H⁺ incident ions. Oxide film growth proceeds significantly faster on Ni(1 1 1) than on Ni(1 0 0), but on both surfaces oxide penetration occurs to depths significantly greater than 100 Å with total exposures of 1200 and 6000 L respectively. The metal/oxide interface is extremely rough, with metallic Ni extending to the surface, even for much thicker oxide films on Ni(1 1 1). On Ni(1 1 1), NiO growth occurs with the (1 0 0) face parallel to the Ni(1 1 1) surface and the close-packed 〈1 1 0〉 directions parallel. On Ni(1 0 0) the MEIS blocking curves cannot be reconciled with a single orientation of NiO(1 0 0) (with the 〈1 1 0〉 directions parallel) on the surface, but is consistent with the substantial orientational disorder (including tilt) previously identified by spot-profile analysis LEED.     

Structural and statistical analysis of Yb/Si(1 1 1) and Eu/Si(1 1 1) reconstructions

We have performed the structural and statistical analysis of Yb/Si(1 1 1) and Eu/Si(1 1 1) surfaces in the submonolayer regime utilizing low-energy electron diffraction and scanning tunneling microscopy (STM). The almost identical series of one-dimensional chain structures (e.g., 3 × 2/3 × 1, 5 × 1, 7 × 1, 9 × 1, and 2 × 1 phases) are found in order of increasing metal coverage for both adsorbed systems, however, only the Eu/Si system reveals the ‘√3’-like reconstruction before the 2 × 1 endpoint phase. The atomic models of chain structures are proposed and discussed. In particular, our results suggest the odd-order n×1 (n=5,7,9,…) intermediate reconstructions to incorporate the Seiwatz chains and honeycomb chains with the proportion of m:1, where . The statistical analysis of STM images is carried out to examine the correlation of atomic rows on Eu/Si and Yb/Si surfaces. It is found that Eu stabilizes more ordered row configuration compared to Yb, which can be explained in terms of indirect electronic interaction of atomic chains or/and different magnetic properties of adsorbed species.     

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.     

Effect of hydrogenation on B/Si(0 0 1)-(1 × 2)

Ab initio calculations, based on pseudopotentials and density functional theory, have been performed to investigate the effect of hydrogenation on the atomic geometries and the energetics of substitutional boron on the generic Si(0 0 1)-(1 × 2) surface. For a single B atom substitution corresponding to 0.5 ML coverage, we have considered two different sites: (i) the mixed Si-B dimer structure and (ii) boron substituting for the second-layer Si to form Si-B back-bond structure, which is energetically more favorable than the mixed Si-B dimer by 0.1 eV/dimer. However, when both of these cases are passivated by hydrogen atoms, the situation is reversed and the Si-B back-bond case becomes 0.1 eV/dimer higher in energy than the mixed Si-B dimer case. For the B incorporation corresponding to 1 ML coverage, among the substitutional cases, 100% interdiffusion into the third layer of Si and 50% interdiffusion into the second layer of Si are energetically similar and more favorable than the other cases that are considered. However, when the surface is passivated with hydrogen, the B atoms energetically prefer to stay at the third layer of the Si substrate.     

NO dissociation pathways on Rh(1 0 0), (1 1 0), and (1 1 1) surfaces: A comparative density functional theory study

The chemisorption and dissociation pathways of NO on the Rh(1 0 0), (1 1 0), and (1 1 1) surfaces are studied by the plane-wave density functional theory (DFT) with CASTEP program. In addition, the electronic and geometrical effects that affect the NO dissociation reactions have been investigated in detail. The calculation results are presented as following: The effective activation energies of the best NO dissociation pathways on the Rh(1 0 0), the Rh(1 1 0), and the Rh(1 1 1) are 0.63, 0.66 and 1.77 eV, respectively. The activity of the Rh planes for NO dissociation is in the order of Rh(1 0 0) ≈ Rh(1 1 0) > Rh(1 1 1). The low dissociation barrier for Rh(1 0 0) and Rh(1 1 0) is associated with the existence of a lying-down NO structure which acts as a precursor for dissociation. By Mulliken population analysis and structure analysis, both electronic and geometrical effects are found to affect the NO dissociation reactions, but the geometrical effect exceed the electronic. The energy decomposition scheme has been used to provide further insight into the NO dissociation reactions. Based on the calculations, the interaction energy between N and O in the transition state on the Rh(1 1 1) is found much larger than that on the Rh(1 0 0) and the Rh(1 1 0). The major differences of should originate from the variation of the bonding competition effect.     

Vacancy charging on Si(1 1 1)-(7 × 7) investigated by density functional theory

The structure and energetics of charged vacancies on Si(1 1 1)-(7 × 7) are investigated using density functional theory calculations supplemented by estimates of ionization entropy. The calculations predict multiple possible charge states for the unfaulted edge vacancy in the adatom layer, although the −2 state is most stable on real Si(1 1 1) surfaces for which the Fermi level lies near the middle of the band gap.     

Bistability in a H-terminated Si(1 0 0)2 × 1 surface obtained by ab initio transport calculations

We have analyzed electron tunneling due to the electric field from a hydrogen-terminated Si(1 0 0)2 × 1 ultrathin film on a metal substrate by density functional transport calculations. We have obtained a hysteresis loop in the tunneling current, which comes from the existence of two electronic structures. Furthermore, we have clarified that, as a condition of bistable electron transport, a double-barrier potential structure is not necessarily required for zero field, because it can be induced by the electric field.     

Vertical bonding distances of PTCDA on Au(1 1 1) and Ag(1 1 1): Relation to the bonding type

The vertical bonding distance of 3,4,9,10-perylene-tetracarboxylic-dianhydride (PTCDA) above the Au(1 1 1) surface has been measured by the normal incidence X-ray standing wave (NIXSW) technique. The carbon skeleton of PTCDA has a vertical distance of D = (3.27 ± 0.02) Å to the Au(1 1 1) substrate. This distance corresponds very nearly to the sum of the van der Waals radii of carbon and gold, suggesting the adsorption to be a physisorptive one. In contrast, the PTCDA/Ag(1 1 1) interface which according to spectroscopic data follows the standard model of chemisorption very closely, shows a considerably smaller bonding distance of D = (2.86 ± 0.01) Å [A. Hauschild, K. Karki, B.C.C. Cowie, M. Rohlfing, F.S. Tautz, M. Sokolowski, Phys. Rev. Lett. 94 (2005) 036106, comment: Rurali et al., Phys. Lett. 95 (2005) 209205, reply: Phys. Rev. Lett. 95 (2005) 209206]. The different vertical adsorption heights of PTCDA on gold and silver are discussed in relation to the different bonding mechanisms on both noble metal surfaces.     

Electron and lattice structure of ultra thin Ag films on Si(1 1 1) and Si(0 0 1)

We studied the low temperature (T ? 130 K) growth of Ag on Si(0 0 1) and Si(1 1 1) flat surfaces prepared by Si homo epitaxy with the aim to achieve thin metallic films. The band structure and morphology of the Ag overlayers have been investigated by means of XPS, UPS, LEED, STM and STS. Surprisingly a (√3 × √3)R30° LEED structure for Ag films has been observed after deposition of 2-6 ML Ag onto a Si(1 1 1)(√3 × √3)R30°Ag surface at low temperatures. XPS investigations showed that these films are solid, and UPS measurements indicate that they are metallic. However, after closer STM studies we found that these films consists of sharp Ag islands and (√3 × √3)R30°Ag flat terraces in between. On Si(0 0 1) the low-temperature deposition yields an epitaxial growth of Ag on clean Si(0 0 1)-2 × 1 with a twinned Ag(1 1 1) structure at coverage’s as low as 10 ML. Furthermore the conductivity of few monolayer Ag films on Si(1 0 0) surfaces has been studied as a function of temperature (40-300 K).