Atomistic simulation of the vacancy diffusion in (0 0 1) surface of MoTa alloy

The formation and diffusion of a single Mo or Ta vacancy in the (0 0 1) surface of the B₂-type MoTa alloy have been investigated by using modified analytical embedded-atom method (MAEAM). The results show that the effect of the surface on the vacancy is only down to the sixth layer. It is easier for the vacancy to form in the first layer. Comparing the migration energy of the vacancy migrating in the intra-layer, to the upper layer and to the nether layer via 2NN jump, we find that the vacancy in the first or second layer is preferred to migrate in intra-layer, and that in the third or fourth layer is favorable to migrate to the upper layer. Although 1NN jump will result in an anti-site so that a disorder in the order alloy, it may also occur due to the much lower migration energy especially for the vacancy in the second and third layer to migrate to the first and second layer, respectively.     

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.     

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.     

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.     

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.     

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.     

The properties and structures of the mono- and the di- vacancy in Cu crystal

The structures and energies of formation and migration of the mono- and di-vacancy in Cu crystal have been described and calculated with modified analytical embedded atom method (MAEAM). The lattice relaxation is considered with molecular dynamics (MD) method at T=0 K. The results show the FN di-vacancy is the most stable and likely occurs in practice from the energy minimization. Compared with the mono-vacancy, the formation energy of the FN di-vacancy is higher than that of a mono-vacancy, but lower than that of two isolated mono-vacancy. The preferred migration mechanism of the FN di-vacancy is multi-jump of either vacancy (rotating the di-vacancy). The calculated migration energy of the FN di-vacancy is lower than that of a mono-vacancy, so the FN di-vacancy is easier to migrate. All of the calculated results are in good agreement with the experimental values.     

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.     

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.     

Ethylene hydrogenation on Ni(1 0 0) surface

The hydrogenation of ethylene on Ni(1 0 0) surface has been studied by TDS. The decrease in the bonding energy with increasing coverage is revealed for both of adsorbed hydrogen and ethylene by the shift of desorption to lower temperatures. Ethane formation is only observed on the preadsorbed hydrogen coverage exceeding 0.5 monolayer (ML), coupled with the growth of H₂ shoulder peak at lower temperatures. Further increase of H coverage to saturation reduces the bonding energy of subsequently adsorbed ethylene by 15 kJ/mol and decreases the saturation coverage of ethylene to about one-third on the clean surface. This leads to the shift of ethane desorption from 250 to 220 K and an appearance of additional ethane peak at 180 K. The latter ethane formation coincides with the hydrogenation of surface ethyl species derived from ethyl iodide as a precursor. It indicates that the rate of ethyl formation on the surface would be comparable to that of subsequent hydrogen addition to the surface ethyl species in the hydrogenation of ethylene when the preadsorbed hydrogen coverage approaches 1.0 ML.     

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.     

Theoretical study of Ni adsorption on the GaN(0 0 0 1) surface

First-principles pseudo-potential calculations within density-functional theory framework are performed in order to study the structural and electronic properties of nickel adsorption and diffusion on a GaN(0 0 0 1)-2×2 surface. The adsorption energies and potential energy surfaces are investigated for a Ni adatom on the Ga-terminated (0 0 0 1) surface of GaN. This surface is also used to study the effect of the nickel surface coverage. The results show that the most stable positions of a Ni adatom on GaN(0 0 0 1) are at the H₃ sites and T₄ sites, for low and high Ni coverage respectively. In addition, confirming previous experimental results, we have found that the growth of Ni monolayers on the GaN(0 0 0 1) surface is possible.     

Ni adsorbate dynamics on a NiO(0 0 1) surface

We present results concerning the dynamical behavior of a Ni³⁺ adsorbate on a NiO(0 0 1) surface obtained by molecular dynamics simulations. In a first place, we examined at low temperature the position of the Ni³⁺ ion as an adatom on the surface and the corresponding modification of its local environment as reflected on the pair-wise radial distribution function. The calculation of the vibrational properties of the adatom by means of the phonon local density of states (LDOS) shows that there is an anisotropic behavior both in the two principal in-plane directions as well as in the direction normal to the surface in accordance with the structural results. We compare the phonon LDOS of the Ni³⁺ adatom with the corresponding results for the Ni²⁺ adatom and the Ni²⁺ surface cations.Static energetic calculations are indicative that the exchange of the Ni³⁺ ion with a surface Ni²⁺ ion could be favorable. Such a behavior is confirmed by results observed at temperatures higher than 700 K where the Ni³⁺ adsorbate is located on a substitutional position on the surface and not on adatom position. The exchange takes place through simple or double exchange mechanisms. The structural and dynamical behavior of the Ni³⁺ ion at the substitution position was investigated in the temperature range 700-2000 K through the calculation of the pair distribution function, the relaxed interlayer relative position (RIRP), mean-square displacements (MSDs) and phonon LDOS. Results show that in comparison with the Ni²⁺ surface ions the Ni³⁺ ion at substitution position is more tightly bound especially in the direction normal to the surface as is indicated by the local structure and the contraction it presents as well as its phonon LDOS. As temperature increases the binding of the Ni³⁺ ion becomes less important as reflected on the physical properties mentioned above.     

Atomic-level investigation of Al and Ni thin film growth on Ni(1 1 1) surface: Molecular dynamics simulation

Using molecular dynamics (MD) simulation, the structural characteristics of Al and Ni thin film growth on Ni(1 1 1) substrate according to the incident energy of adatoms were investigated. In case of Al on Ni(1 1 1), Al adatoms were grown basically through the layer-by-layer growth mode. On the other hand, Ni thin films on Ni(1 1 1) surface at low incident energy were shown to favor island growth. The steering effect due to atomic attraction, which results in rougher surface, was significantly observed at low incident energy. The growth mode of Ni film was, however, changed to follow layer-by-layer growth mode for the incident energy of 6 eV. The different aspects of surface morphology between Al and Ni deposition on Ni(1 1 1) surface could be successfully explained by the surface diffusion and impact cascade diffusion.     

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

High-resolution angle-resolved photoemission study of Ni(1 1 1) surface state

High-resolution angle-resolved photoemission spectroscopy (ARPES) has been conducted to study the Shockley state (SS) in ferromagnetic Ni(1 1 1) located at the point of the surface Brillouin zone. We have determined the Fermi wave vector and Fermi energy of the state with excitation photon energies of hν = 6.9-27.5 eV. On the basis of ARPES spectral shape analyses, we have found significant electron-electron interaction in the SS.     

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.     

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.     

Mixed layer formation of copper overlayers on Ni(1 1 0) surface

Copper overlayer formation on the Ni(1 1 0) surface was studied by scanning tunneling microscopy (STM) in an ultrahigh vacuum. Atom-resolved STM images showed that initially deposited Cu is replaced with surface Ni atoms forming atom-size depressions on the Ni(1 1 0) terraces and a Ni-rich quasi-one-dimensional island along the direction. Further Cu deposition yields a mosaic structure on the islands, indicating Cu/Ni mixed layer formation. From the quantitative measurement of the Cu/Ni ratio on the substrate and the islands, impinging Cu will be replaced with surface Ni whereas expelled Ni and directly impinging Cu to the island form the mixed island. The number of Cu atoms in the islands, however, more than the directly impinging Cu, indicate significant Cu/Ni replacement at the periphery of the island.