Atomistic simulation of the self-diffusion in Fe (1 1 1) surface

The formation energies, the intra- and inter-layer self-diffusion activation energies of a single vacancy in Fe (1 1 1) surface have been investigated with the modified analytical embedded-atom method (MAEAM). The results show that the effect of the surface is down to the sixth layer for the formation and intra-layer migration of the vacancy. It is easier for a vacancy to form and to migrate in intra-layer in the first (especially), the second and the third layer. For inter-layer migration, a single vacancy in each of the first six layers is favorable to migrate to the upper layers. On the contrary, a single vacancy in the seventh, the eighth and the ninth layers is favorable to migrate to the lower layers.     

Atomistic simulation of the vacancy in Ni (1 1 0) surface

Both the formation energies and the intra- and inter-layer diffuse activation energies of a vacancy in the first six atomic planes of Ni (1 1 0) surface have been investigated by means of molecular dynamics (MD) in conjunction with the semi-experiential many-body potential of the modified analytical embedded-atom method (MAEAM). The results show that the effect of the surface is only down to the fifth-layer. It is easer for a vacancy in the first or second layer to form and to migrate in intra-layer. For the inter-layer migration, a vacancy in the second or third layer is favorable to migrate to the upper layer, this is not the case for a vacancy in the fourth or fifth layer.     

Vacancy diffusion in Cu∑ = 5[0 0 1] twist grain boundary

Both the formation energies and the diffusive activation energy of a single vacancy migrating intra- and inter-layer in the first four atomic planes near Cu∑ = 5[0 0 1] twist GB have been investigated by means of MD in conjunction with MAEAM. The effects of the GB on the vacancy formation and migration are only to the third layer. The vacancy is favorable to be formed on the un-coincident site in the first, second and third layers near the GB plane and this case is enhanced successively following the third, second and first layers. A single vacancy either on un-coincident site or on coincident site in the forth, third and second layers is favorable to migrate to un-coincident site (its first-nearest-neighbor) in its adjacent layer near the GB. But for the first layer, the favorable migration path of the vacancy on the un-coincident site is between un-coincident sites of the first layer or to its nearest-neighbor of the first layer in the rotating grain, which is not the case for the vacancy on the coincident site ‘1’ that is migrated difficultly. So, there are collective tendency of the vacancy in the GB.     

Structural properties and diffusion processes of the Cu3Au (0 0 1) surface

The surface relaxation and surface energy of both the mixed AuCu and pure Cu terminated Cu₃Au (0 0 1) surfaces are simulated and calculated by using the modified analytical embedded-atom method. We find that the mixed AuCu termination is energetically preferred over the pure Cu termination thereby the mono-vacancy diffusion is also investigated in the topmost few layers of the mixed AuCu terminated Cu₃Au (0 0 1) surface. In the mixed AuCu terminated surface the relaxed Au atoms are raised above Cu atoms for 0.13 Å in the topmost layer. All the surface atoms displace outwards, this effect occurs in the first three layers and changes the first two inter-layer spacing. For mono-vacancy migration in the first layer, the migration energies of Au and Cu mono-vacancy via two-type in-plane displace: the nearest neighbor jump (NNJ) and the second nearest neighbor jump (2NNJ), are calculated and the results show that the NNJ requires a much lower energy than 2NNJ. For the evolution of the energy requirements for successive nearest neighbor jumps (SNNJ) along three different paths: circularity, zigzag and beeline, we find that the circularity path is preferred over the other two paths due to its minimum energy barriers and final energies. In the second layer, the NN jumps in intra- and inter-layer of the Cu mono-vacancy are investigated. The calculated energy barriers and final energies show that the vacancy prefer jump up to a proximate Cu site. This replacement between the Cu vacancy in the second layer and Cu atom in the first layer is remunerative for the Au atoms enrichment in the topmost layer.     

Computer simulation study of self-diffusion in Pd(0 0 1) surface

Both the formation energies and the intra- and inter-layer activation energies of self-diffusion of a single vacancy in the first six planes of Pd(0 0 1) surface have been investigated by means of molecular dynamics (MD) in conjunction with the semi-experiential many-body potential of the modified analytical embedded-atom method (MAEAM). The results show that the effect of the surface on the vacancy is only down to the fifth-layer. It is easer for a single vacancy to form and to migrate in the first layer. Furthermore, the vacancy in the second layer is favorable to migrate to the first layer. This is in agreement with the experimental results that the first layer has the highest concentration of the vacancy.     

Structure and composition of the NiPd(1 1 0) surface

The NiPd(1 1 0) alloy surface was studied using low energy electron diffraction to measure the structure and composition of the first three atomic layers. The surface layer is highly enriched in Pd and has a significantly buckled structure. The second layer is also buckled, with displacements even larger than the surface layer. The second layer also exhibits intralayer segregation (chemical ordering), with alternate close-packed rows of atoms being Ni enriched and Pd enriched. The third layer has a structure and composition close to that of the bulk alloy. These results are compared with results for the other low index faces of NiPd, the extensive literature on NiPt alloy surfaces, and the growing body of theoretical literature for NiPd alloy surfaces.     

Study of vacancy on diamond (1 0 0) (2 × 1) surface from first-principles

Structure and energy related properties of neutral and charged vacancies on relaxed diamond (1 0 0) (2 × 1) surface were investigated by means of density functional theory. Calculations indicate that the diffusion of a single vacancy from the top surface layer to the second layer is not energetically favored. Analysis of energies in charged system shows that neutral state is most stable on diamond (1 0 0) (2 × 1) surface. The multiplicity of possible states can exist on diamond (1 0 0) surface in dependence on the surface Fermi level, which supports that surface diffusion of a vacancy is mediated by the change of vacancy charge states. Analysis of density of states shows surface vacancy can be effectively measured by photoelectricity technology.     

Study on low-energy bombardment of Au (1 1 1) by noble metal atoms with molecular dynamics simulations

The low-energy bombardment of Au (1 1 1) surface by noble metal atoms is studied with molecular dynamics (MD) simulations. With the incident-energy dependence of adatom yields, sputtering yields, and vacancy yields for different projectiles, we find that the implantation of projectiles in shallow layers below surface can be distinguished by subplantation (in the first and second layers) and implantation (deeper than the third layer). The transition from subplantation to implantation occurs at the incident energy of about 45 eV for the low-energy bombardment of noble metal atoms on Au (1 1 1). The incident-energy dependence of defect yields is obviously different for the subplantation and implantation of projectiles. Based on our MD simulations, we discuss the influence of low-energy bombardment on film growth and the guide to the search for optimum deposition parameters.     

Hydrogen-induced metallization of the 3C-SiC(0 0 1)-(3 × 2) surface

A surprising metallization of the SiC(0 0 1)-(3 × 2) surface induced by hydrogen adsorption was discovered in recent experiments. The effect was ascribed to dangling bonds created on the third layer of the surface system by H adsorption and stabilized by steric hindrance. We have investigated the surface metallization by density functional calculations. Our total-energy minimizations show that dangling bonds on the third layer are very unstable. Instead, H adatoms form angular Si-H-Si bonds on the third layer after the asymmetric dimers on the top layer have been saturated by H forming monohydrides. The novel Si-H-Si bonds on the third layer give rise to a metallic surface, indeed. But the mechanism for metallization is very different from the one suggested originally. Likewise, H atoms can also occupy bridge positions in angular Si-H-Si bonds on the second layer and induce metallization, as well. In addition to monohydrides on the top-layer dimers, we have also investigated dihydride surfaces with additional H on the second and/or third layer. The dihydride surface structure with H adsorbed on both the second and third layers is energetically most favorable and is also metallic. In all three cases the new Si-H-Si bonds are the origin of the surface metallization while its nature is somewhat more intricate, as will be discussed.     

An ordered surface ternary alloy of a c(6 × 4) structure formed on Cu(0 0 1) by substitutional coadsorption of Mg and Bi

A c(6 × 4) structure formed on Cu(0 0 1) by the coadsorption of Mg and Bi atoms at room temperature has been determined by a tensor low energy electron diffraction analysis. It is an ordered surface ternary alloy with a thickness of single layer, in which Mg, Bi and Cu atoms are mixed in the top layer. In the primitive unit cell, there are one Mg, four Bi, six Cu atoms and one vacancy in the top layer, and substituted Mg and Bi atoms form MgBi₄ plane clusters being arranged in the c(6 × 4) order. Structural parameters show that Mg-Bi bond distances in the MgBi₄ cluster are 3.01 and 3.07 Å, which are shorter than the summation of metallic radii of Mg and Bi. It is concluded that a direct, attractive interaction between Mg and Bi atoms plays critical role in the formation of the c(6 × 4) structure.     

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.     

Formation and structure of a (√19 × √19)R23.4°-Ge/Pt(1 1 1) surface alloy

An ordered (√19 × √19)R23.4°-Ge/Pt(1 1 1) surface alloy can be formed by vapor depositing one-monolayer Ge on a Pt(1 1 1) substrate at room temperature and subsequently annealing at 900-1200 K. The long-range order of this structure was observed by low energy electron diffraction (LEED) and confirmed by scanning tunneling microscopy (STM). The local structure and alloying of vapor-deposited Ge on Pt(1 1 1) at 300 K was investigated by using X-ray Photoelectron Diffraction (XPD) and low energy alkali ion scattering spectroscopy (ALISS). XPS indicates that Ge adatoms are incorporated to form an alloy surface layer at ∼900 K. Results from XPD and ALISS establish that Ge atoms are substitutionally incorporated into the Pt surface layer and reside exclusively in the topmost layer, with excess Ge diffusing deep into the bulk of the crystal. The incorporated Ge atoms at the surface are located very close to substitutional Pt atomic positions, without any corrugation or “buckling”. Temperature Programmed Desorption (TPD) shows that both CO and NO adsorb more weakly on the Ge/Pt(1 1 1) surface alloy compared to that on the clean Pt(1 1 1) surface.     

Growth of Ce on Rh(1 1 1)

We have studied the growth of cerium films on Rh(1 1 1) using STM (scanning tunneling microscopy), LEED (low energy electron diffraction), XPS (X-ray photoelectron spectroscopy) and AES (Auger electron spectroscopy). Measurements of the Ce films after room temperature deposition showed that Ce is initially forming nanoclusters in the low coverage regime. These clusters consist of 12 Ce atoms and have the shape of pinwheels. At a coverage of 0.25 ML (monolayer, ML) an adatom layer with a (2 × 2) superstructure is observed. Above 0.4 ML, Rh is diffusing through pinholes into the film, forming an unstructured mixed layer. Annealing at 250 °C leads to the formation of ordered Ce-Rh compounds based on the bulk compound CeRh₃. At a coverage of 0.1 ML, small ordered (2 × 2) surface alloy domains are observed. The exchanged Rh atoms form additional alloy islands situated on the pure Rh(1 1 1) surface, showing the same (2 × 2) superstructure as the surface alloy. At a coverage of 0.25 ML, the surface is completely covered by the surface alloy and alloy islands. The (2 × 2) structure is equivalent to a (1 1 1)-plane of CeRh₃, contracted by 6%. Annealing a 1 ML thick Ce layer leads to a flat surface consisting of different rotational domains of CeRh₃(1 0 0). The Rh needed for alloy formation comes from 50 Å deep pits in the substrate. Finally we show that LEIS (low energy ion scattering) is not suitable for the characterization of Ce and CeRh films due to strong effects of neutralization.     

Structure and stability of Sb/Au(1 1 0)-c(2 × 2) surface phase

Adsorption of 0.5 monolayers (ML) of Sb on the Au(1 1 0) surface resulted in the formation of a c(2 × 2) surface reconstruction. Analysis of surface X-ray diffraction data by a direct method revealed the existence of an ordered substitutional surface alloy, with every other hollow site occupied by Au and Sb atoms. Quantitative conventional χ² refinement showed a contraction of 0.12 ± 0.03 Å in the spacing of the first Au layer to the second, an expansion of 0.13 ± 0.03 Å in the second-to-third layer distance, and an inward Sb displacement (rumpling) of 0.21 ± 0.04 Å. This surface phase proved to be extremely robust, with the long-range order of this arrangement remaining up to substrate temperatures of 900 K.     

Adsorbate-substrate bonding and the growth of naphthalene thin films on Ag(1 1 1)

The adsorption kinetics, energetics and growth of naphthalene thin films, from submonolayer to about 10 layers, on a Ag(1 1 1) surface at low temperature in a ultrahigh vacuum chamber are examined by using temperature programmed desorption spectroscopy. The first layer adsorption occurs with a desorption energy of 85 ± 5 kJ/mole and results in an interface dipole of 5 ± 1 D, from charge transfer of approximately 0.2 e from naphthalene to Ag. The surface dipole induced inter-adsorbate repulsion causes the lowering of the adsorption energy within the first layer near the saturation coverage so that the second layer deposition begins before the completion of the first layer. The second layer is a metastable phase with desorption energy, 74 ± 3 kJ/mole, smaller than the multilayer desorption energy of 79 ± 5 kJ/mole. Fractional order desorption kinetics were found for both the metastable and the multilayer phases, suggesting desorption from 2-D islanding and 3-D islanding, respectively.     

Collision-induced migration of CO on Pt(9 9 7) surface

Lateral displacement of adsorbates induced by collisions with energy-controlled rare gas atoms was examined in an ultra high vacuum chamber using Fourier-transform infrared (FTIR) spectroscopy and a supersonic molecular beam apparatus. A stepped Pt(9 9 7) surface was exposed to CO molecules and subsequently to energy-controlled Ne or Ar atoms. There was no change in the CO stretching mode region of the FTIR spectrum of the Pt(9 9 7) surface after Ne atoms having an average translational energy of 0.23 eV were collided with it. However, when Ne atoms having an average translational energy of 0.56 eV were collided with the surface, the intensity of the peak assigned to the CO stretching mode at terrace sites decreased, while that at step sites increased with increasing the exposure to the Ne atoms. This is the demonstration of collision-induced migration, showing that CO molecules adsorbed at the terrace sites migrate laterally to the step sites upon collision with high-energy Ne atoms. In addition, the experimental results demonstrate the existence of an additional energy barrier for jumps across the steps. This investigation demonstrates an advantage of using a molecular beam for studying adsorbate migration.     

Enhanced reactivity and selectivity in oxidation of Cu(1 0 0) and α-Cu-Al(5 at.%)(1 0 0) surfaces studied by electron and ion spectroscopies

In this study we investigate the influence of alloying on the reactivity and bonding of oxygen on α-Cu-Al(5 at.%)(1 0 0) oriented single crystal surfaces by X-ray photoelectron spectroscopy (XPS), ultra-violet spectroscopy (UPS) and low energy ion scattering (LEIS) spectroscopy, at room temperature. It was found that alloying results in an enhanced reactivity of both Cu and Al sites in comparison with the pure metals. According to adsorption curves calculated from XPS, saturation of the alloy surface occurs for exposures of ∼15 L. At saturation the total amount of adsorbed oxygen is similar for the alloy and pure copper surfaces. It was determined that first mostly Al sites are oxidized, followed by simultaneous oxidation of Cu and Al sites. At saturation the amount of oxygen bonded to Cu sites is ∼1.7 larger then that bonded to Al sites. From a comparison of the XPS and LEIS data analysis as a function of oxygen exposure it was found that oxidation of α-Cu-Al(5 at.%)(1 0 0) alloy is a multi-stage process with fast and slow stages. These stages involve an interplay of chemisorption, sub-surface diffusion of oxygen and Al segregation. UPS measurements show an increase in the work function of the alloy surface with oxygen adsorption. This is a contrast to pure Cu surfaces where the work function decreases at the initial stages of oxidation followed by an increase with oxygen exposure. Annealing to 400 °C drives the oxidized alloy surface into its thermodynamic state resulting in the formation of an aluminum oxide layer. Possible mechanisms to explain the enhanced reactivity of the alloy surface compared to that of pure copper are suggested and discussed.     

Structural properties of oxygen on InN(0 0 0 1) surface

First-principles calculations are performed to study the various structures of oxygen (O) adsorbed on InN(0 0 0 1) surfaces. It is found that the formation energy of O on InN(0 0 0 1) decreases with decreasing oxygen coverage. Of all the adsorbate induced surface structures examined, the structure of InN(0 0 0 1)-(2 × 2) as caused by O adsorption at the H₃ sites with 0.25 monolayers coverage is most energetically favorable. Meanwhile, nitrogen (N) vacancy can form spontaneously. Oxygen atoms may also substitute N atoms, or accumulate at the voids inside InN film or simply stay on the surface during growth. The oxygen impurity then acts as a potential source for the n-type conductivity of InN as well as the large energy band gap measured.     

The structure of Sb(1 1 1) determined by photoelectron diffraction

The Sb(1 1 1) 4d_5/2 core level is found to contain three components. Using photoelectron diffraction, these are assigned to photoemission from the first layer, the second layer and the bulk, respectively. The binding energy for the first and second layer atoms is found to be 120 meV lower and 330 meV higher than for the bulk atoms, respectively. As a by-product of this assignment, the geometric structure of the surface is determined. No substantial relaxations are found.     

Atomic geometry determination of FeO(0 0 1) grown on Ag(0 0 1) by low energy electron diffraction

A low energy electron diffraction (LEED) investigation of the structure of the surface of an FeO(0 0 1) thin film grown on Ag(0 0 1) is presented. The results show that this surface has an almost bulk termination structure with a very small rumple on the first layer, which agrees with the structure found in other studies carried out on the (0 0 1) surface of oxides that have rock-salt structure. Evidences that may support a linear behaviour of the topmost layer rumple with the oxide lattice constant are also discussed.