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.     

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.     

Diffusion of single adatom Cu on Cu (0 0 1) and (1 1 0) surfaces

With static relaxation, the surface diffusion activation energies of a single Cu adatom migrated by both atomic exchange and hopping mechanisms and the forces acted on the diffusing adatom from other atoms of Cu (0 0 1) or (1 1 0) surface are calculated by using the MAEAM. When adatom migrated on Cu (0 0 1) or (1 1 0) surface, the increment curves of the system energy by hopping mechanism are symmetrical and the saddle points are in the midpoints of the migration path, but the ones by the exchange mechanism are dissymmetrical and the saddle points are always close to the initial hole positions of the adatom and away from the initial equilibrium positions of the exchanged atom. From minimization of both the diffusion activation energy and the force acted on the diffusing adatom from other atoms, we found that, on Cu (0 0 1) surface the favorable diffusion mechanism is hopping mechanism, however, on Cu (1 1 0) surface, hopping via long bridge is easier than the exchange mechanism but the hopping via short bridge is more difficult than the exchange mechanism.     

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.     

Collision-induced diffusion and vacancy migration in alkanethiol monolayers on Au(1 1 1)

Scanning tunneling microscopy is used to characterize the collision-induced migration of molecules within well-ordered octanethiol and nonanethiol self-assembled monolayers. A seeded molecular beam is used to create xenon atoms with a kinetic energy of 1.3 eV, and collisions with these atoms cause measurable changes in alkanethiol monolayer surface structure. Migration rates are calculated and compared for molecules in close-packed domains, at domain-boundary defects, and along the perimeter of vacancy-island defects. The number of nearest-neighbor molecules (within the 5 Å lattice distance) is strongly predictive of molecular stability with respect to xenon bombardment, and the overall dependence of stability on nearest neighbors is well fit by a simple exponential curve. The incident direction of the molecular beam is not observed to influence the direction of molecular motion; however, in some cases, migration correlates strongly to surface lattice directions. Finally, there is no evidence that substrate restructuring or gold-atom diffusion accompanies alkanethiol migration under these non-equilibrium conditions.     

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.     

Kinetic length and step permeability on the Si(1 1 1) (1 × 1) surface

Information on the kinetic regime of step motion and step permeability on the Si(1 1 1) (1 × 1) surface has been obtained from observations of island decay that were made with low energy electron microscopy. Island area during decay exhibits the expected power law dependence on time, with exponent, α, that is a qualitative indicator of the kinetic regime. A new method is presented for determining the kinetic length quantitatively from measurements of the decay exponent in the symmetric island decay geometry on top of a larger concentric circular island. Using this approach, we determine the kinetic length on the Si(1 1 1) (1 × 1) surface at 1163 K to be d ∼ 75a, where a is the lattice constant. It is shown that this result locates step motion firmly in the diffusion limited regime. Mass conservation of decaying island stacks is also observed at this temperature, which indicates that steps are effectively impermeable in the context of diffusion limited step kinetics.     

Site selective atomic chlorine adsorption on the Si(1 1 1)7 × 7 surface

The spontaneous dissociation of trichloroethylene molecules on the Si(1 1 1)7 × 7 surface was investigated using STM. Chlorine atoms were identified by using voltage dependent imaging and by observing voltage dependent tip-induced diffusion. At low coverage, we identify one chlorine that dissociates and binds to an adatom, leaving a nearby chlorovinyl group as the other product bound to the surface. Chlorine atoms show strong site selectivity for corner adatoms and some preference for the faulted half of the unit cell. This result differs significantly from previous studies of chlorine on this surface and a site-selective mobile precursor model is used to explain this discrepancy. The observed site-selectivity is consistent with the high electronegativity value for chlorine.     

First-principles calculations of C diffusion through the surface and subsurface of Ag/Ni(1 0 0) and reconstructed Ag/Ni(1 0 0)

Density functional theory calculations are performed to investigate the C diffusion through the surface and subsurface of Ag/Ni(1 0 0) and reconstructed Ag/Ni(1 0 0). The calculated geometric parameters indicate the center of doped Ag is located above the Ni(1 0 0) surface owing to the size mismatch. The C binding on the alloy surface is substantially weakened, arising from the less attractive interaction between C and Ag atoms, while in the subsurface, the C adsorption is promoted as the Ag coverage is increased. The effect of substitutional Ag on the adsorption property of Ni(1 0 0) is rather short-range, which agrees well with the analysis of the projected density of states. Seven pathways are constructed to explore the C diffusion behavior on the bimetallic surface. Along the most kinetically favorable pathway, a C atom hops between two fourfold hollow sites via an adjacent octahedral site in the subsurface of reconstructed Ag/Ni(1 0 0). The “clock” reconstruction which tends to improve the surface mobility, is more favorable on the alloy surface because the c(2 × 2) symmetry is inherently broken by the Ag impurity. As a consequence, the local lattice strain induced by the C transport is effectively relieved by the Ag-enhanced surface mobility and the C diffusion barrier is lowered from 1.16 to 0.76 eV.     

Interaction of small vacancy clusters with (1 1 4) twin-boundary in gold

The molecular dynamics simulation technique with many-body and semi-empirical potentials is used to calculate the (1 1 4) twin-boundary in gold at different temperatures. Relaxations are found on both sides of the interface with the same magnitude and the phenomenon of coalescence is observed near the interface. The interactions of single-, di- and tri-vacancies with twin-interface at 300 K on mirror and off-mirror sites are calculated. Off-mirror arrangements are favorable for all vacancy clusters, except for the single-vacancy cluster, which is less repulsive on the mirror site. Vacancy clusters energetically prefer to lie at planes closest to the (1 1 4) interface rather than away from it. The effect of temperature on interaction behavior is also calculated.     

Formate-induced destabilization of the Cu(1 1 0) surface

Scanning tunneling microscopy (STM) images show that adsorbed formate has a profound affect on the step edges of Cu(1 1 0) surfaces at room temperature. For low exposures, the presence of formate enhances step fluctuations as confirmed by a correlation function analysis. For formate coverages approaching 0.5 monolayers, drastic restructuring of step edges is observed. Quantum chemical calculations help to explain this behavior.     

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.     

Analysis of energy and force of Pt adatom on Pt (0 0 1) surface by MAEAM

The energy and perpendicular force of a Pt adatom on Pt (0 0 1) surface have been calculated by MAEAM. With increasing the distance of the adatom from the surface, the energy and force maps can be classified into four regions: repulsive region, transformed region, strongly attractive region and weakly attractive region. In repulsive region, the maximum (minimum) values of the energy and repulsive force appear on the top (hole) of the first layer atoms of Pt (0 0 1) surface due to stronger pair-potential interaction. In other regions, the energy and force maps are more complicated than those in repulsive region due to the effects of the many body interactions and nonspherical distribution of the electrons of the atoms in crystal. The most stable position is 0.1664 nm above the hole of the first layer atoms for a Pt adatom on Pt (0 0 1) surface.     

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].     

Surface mass diffusion and step stiffness on an anisotropic surface; Mo(0 1 1)

Results of step fluctuation experiments for Mo(0 1 1), using low-energy electron microscopy, are re-examined using recently developed procedures that offer accurate coefficients of surface mass diffusion. By these means, surface diffusion D_s is documented at T/T_m ∼ 0.5, while the crossover to relaxation driven by bulk vacancy diffusion is inferred for T/T_m ∼ 0.6. Here, T_m is the melting temperature T_m = 2896 K. We obtain D_s = 4 × 10⁻⁴ exp(−1.13 eV/k_BT) cm²/s for the temperature interval 1080-1680 K. Possible indications of diffusion along step edges appear for T/T_m ∼ 0.4. The same measurements of step fluctuation amplitudes determine also the step stiffness, which by symmetry is anisotropic on Mo(0 1 1). It is shown that three independent procedures yield mutually consistent step stiffness anisotropies. These are (1) step fluctuation amplitudes; (2) step relaxation rate anisotropies; and (3) the observed anisotropies of islands in equilibrium on the Mo(0 1 1) surface. The magnitude of the step stiffness obtained from step edge relaxation is consistent with earlier measurements that determine diffusion from grain boundary grooving.     

The surface stress of the (1 1 0) and (1 0 0) surfaces of rutile and the effect of water adsorbents

The surface stress on clean TiO₂ (1 1 0) and (1 0 0) surfaces, and those with four types of adsorbent - (i) molecularly adsorbed water, (ii) dissociatively adsorbed water, (iii) dissociatively adsorbed water at an oxygen vacancy, and (iv) adsorbed hydrogen - was investigated in the framework of density functional theory using a slab model. The calculations were intended to rationalize the effect of the artificially introduced stress that occurs in experimentally photoinduced hydrophilicity. Tensile stress was observed for a clean (1 1 0) surface, and a mixture of tensile and compressive stress for a clean (1 0 0) surface. The adsorbate-induced surface stresses were analyzed in terms of the sixfold coordinated character of the surface titanium atoms, hydrogen bonds between the adsorbents and the bridging oxygen atoms, and the change in electron density in the vicinity of the surface.     

Faceting of nanoscale fingers on the (1 1 1) surface of gold

Gold fingers, one atomic layer (0.25 nm) high, 4-5 nm wide, and several hundred nm long are formed on the (1 1 1) surface of gold at room temperature by a combination of atomic manipulation and surface self-organisation. Each finger has two parallel edges (type A and type B, respectively) running along its length. The type A step is found to have higher step energy and become nanofaceted when disturbed by either thermal energy or the electric field under the STM tip, leading to the transformation of fingers to “nano-knives”. Our findings reveal the important role of step energy in the process of nanostructure fabrication on surfaces. The gold fingers also provide an ideal system for the investigation of meta-stable nanostructures.     

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.     

Hydrogen adsorption and diffusion on Pd(1 1 1)

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