Molecular dynamics simulation of low-energy bombardment of Pt (1 1 1) surface by Cu atoms with various incident angles

The low-energy bombardment of Pt (1 1 1) surface by Cu atoms with various incident angles (θ) is studied with MD simulations. In the case of near-normal incidence (θ≤20°), the result of energy deposition is similar to that of θ=0°. In contrast, in the case wherein the incident angles are higher than 60°, the incident atom cannot penetrate through the first layer and is scattered directly on the surface. The low-energy deposition has no obvious effect on the substrate. For 20°≤θ≤60°, the oblique incidence contributes to uniformity of nucleation and layer-by-layer growth of film as well as the layer-by-layer removal of atoms in the surface layers. Based on our MD simulations, the mechanism behind the deposition and thin film formation is related to the horizontal component and the vertical component of the impact momentum.     

低能原子沉积在Pt(111)表面的分子动力学模拟

采用嵌入原子方法的原子间相互作用势,利用分子动力学方法模拟了六种贵金属原子(Ni,Pd,Pt,Cu,Ag,Au)分别在Pt(111)表面低能沉积的动力学过程.结果表明:随着入射能量从0.1eV升高到200eV,基体表面原子是按层迁移的,沉积过程对基体表面的影响和沉积原子在基体表层的作用均存在两个转变能量(ET1≈5eV,ET2≈70eV).当入射能量低于5eV时,基体表面几乎没有吸附原子和空位形成,沉积原子在基体表层几乎没有注入产生;当入射能量在5—70eV范围内时,沉积原子在基体表层有注入产生,其注入深度小于两个原子层,即为亚注入,此时吸附原子主要由基体表层原子形成,基体表面第三层以下没有空位形成;当入射能量高于70eV时,沉积原子的注入深度大于两个原子层,将会导致表面以下第三层形成空位,并且空位产额随入射能量的升高而急剧增加.基于分子动力学模拟的结果,对低能沉积作用下的薄膜生长以及最优沉积参数的选择进行了讨论.     

Structural and electrical analysis of S ion bombarded p-InP(1 0 0)

The chemical state of sulfur and surface structure on low-energy S⁺ ion-treated p-InP(1 0 0) surface have been investigated by high-resolution X-ray photoelectron spectroscopy (XPS) and low-energy electron diffraction (LEED). S⁺ ion energy over the range of 10-100 eV was used to study the effect of ion energy on surface damage and the process of sulfur passivation on p-InP(1 0 0) by S⁺ ion beam bombardment. It was found that sulfur species formed on the S⁺ ion-treated surface. The S⁺ ions with energy above 50 eV were more effective in formation of In-S species, which assisted the InP surface in reconstruction into an ordered (1 × 1) structure upon annealing. After taking into account physical damage due to the process of ion bombardment, we found that 50 eV was the optimal ion energy to form In-S species in the sulfur passivation of p-InP(1 0 0). The subsequent annealing process removed donor states that were introduced during the ion bombardment of p-InP(1 0 0). Results of theoretical simulations by Transport of Ions in Materials (TRIM) are in accordance with those of experiments.     

The energy and structure of (0 0 1) twist grain boundary in noble metals

The relaxed energy and structure of (0 0 1) twist grain boundary (GB) in noble metals Au, Ag and Cu are simulated by the MAEAM. In-boundary translation between two adjacent grains results in a periodic energy variation and the period is a square with the side length L_Σ/Σ. The lowest energy appears when the two grains are translated relatively to either corner or center of the periodic square. The relaxed GB energy increases smoothly for low-angle boundaries and levels off for larger-angle boundaries except a cusp appeared at θ = 36.87° (Σ = 5). After relaxation, the symmetry of the GB structure is not changed but the displacement of the atoms parallel to the GB plane decreases with increasing the distance of the atoms from the GB plane.     

Growth and thermal stability of ultra-thin Ag and Au layers on Mo(1 1 1) surface

The influence of temperature on the growth process of ultra-thin Ag and Au layers on the Mo(1 1 1) surface was investigated. At 300 K growth of the Stranski-Krastanov type was found for Ag; for Au growth of the monolayer plus simultaneous multilayers type was found, where a base layer is one physical layer. The first three geometrical adsorbed layers for Ag are thermally stable. For annealed Au layers triangle features with base side length from 15 to 35 Å were formed for θ < 6 monolayer (ML), and for θ > 6 ML part of the Au formed a flat adlayer with Au atoms grouped in equilateral triangles with side length 7 Å. The presence of Au layers does not cause faceting, layers are not smooth which could be caused by the fact that Au does not wets the substrate. For Ag thick layers reversible wetting/non-wetting transition was observed at 600 K. Ag layers on Mo(1 1 1) surface did not lead to faceting.     

Bi induced step-flow growth in the homoepitaxial growth of Au(1 1 1)

Homoepitaxial growth of Au on Bi-covered Au(1 1 1) was studied at room temperature using reflection high-energy electron diffraction (RHEED) and Auger electron spectroscopy (AES). From observations of RHEED it is found that the Au(1 1 1) (23 × 1) reconstruction structure changes to a (1 × 1) by about 0.16-0.5 ML deposition of Bi and to a (2√3 × 2√3)R30° by about 1.0 ML deposition of Bi, respectively. The surface morphology evolution by Bi deposition leads to a change of Au homoepitaxial growth behavior from layer-by-layer to step flow. This indicates that the surface diffusion distance of Au atoms on the Bi-precovered (1 × 1) and (2√3 × 2√3)R30° surfaces is longer than that on the Au(1 1 1) (23 × 1) clean surfaces. A strong surface segregation of Bi was found at top of surface. It is concluded that Bi atoms acted as an effective surfactant in the Au homoepitaxial growth by promoting Au intralayer mass transport.     

Anharmonic effects on B2-FeAl(1 1 0) surface: A molecular dynamics study

Using molecular dynamics simulations and the analytic embedded-atom method (AEAM), the surface anharmonicity of B2-FeAl(1 1 0) has been studied in the temperature range from 0 K to 1400 K. The temperature dependence of the interlayer spacing, mean square vibrational amplitudes, surface phonon frequencies and line-widths, and layer structure factor have been calculated. The obtained results indicate that the anharmonic effects are small in the temperature range from 0 K to 900 K. The temperature dependences of the interlayer spacing indicates that the rippling effect of the B2-FeAl(1 1 0) surface is exhibited by the contraction of Fe surface atoms and the expansion of Al atoms, which persists at high temperatures. The temperature dependence of the layer structure factors shows that the B2-FeAl(1 1 0) surface does not disorder until the temperature of 1300 K.     

Theoretical and experimental studies of the structure and dynamics of the CaF2(1 1 1) surface

The structure and dynamics of the CaF₂(1 1 1) surface were investigated by means of low-energy electron diffraction (LEED) and molecular dynamics (MD) simulations at 300 K. LEED beam intensities were recorded as a function of electron energy and were analyzed with the tensor LEED approach. Positions as well as mean square amplitudes of the ions in the first layers were fitted to the experimental I(E) curves. According to both LEED and MD, the CaF₂(1 1 1) surface structure is similar to the bulk-terminated structure with only small relaxation of the outermost ions. Moreover, both methods show an enhancement of vibrational amplitudes in the outermost F-Ca-F triple layer.     

Structure and electronic properties of gold adsorbed on Ti(0 0 0 1)

The Au/Ti(0 0 0 1) adsorption system was studied by low energy electron diffraction (LEED) and photoemission spectroscopy with synchrotron radiation after step-wise Au evaporation onto the Ti(0 0 0 1) surface. For adsorption of Au at 300 K, no additional superstructures were observed and the (1 × 1) pattern of the clean surface simply became diffuse. Annealing of gold layers more than 1 ML thick resulted in the formation of an ordered Au-Ti surface alloy. Depending on the temperature and annealing time, three surface reconstructions were observed by LEED: (√3 × √3) R30°, (2 × 2) and a one-dimensional incommensurate (√3 × √3) rectangular pattern. The Au 4f core level and valence band photoemission spectra provided evidence of a strong chemical interaction between gold and titanium. The data indicated formation of an intermetallic interface and associated valence orbital hybridization, together with diffusion of gold into the bulk. Au core-level shifts were found to be dependent on the surface alloy stoichiometry.     

Angular distributions of atoms sputtered from NiAl{1 1 0} and Ni{1 0 0}

Using the laser post-ionization surface analysis technique, I have for the first time studied angular distributions of Ni and Al atoms sputtered from NiAl{1 1 0}. Emission angular distributions from Ni{1 0 0} have also been measured. I have observed preferential emissions of Ni and Al atoms along 〈1 1 1〉 and 〈1 0 0〉crystallographic directions for NiAl{1 1 0} and of Ni atoms along 〈1 1 0〉 and 〈1 0 0〉 directions for Ni{1 0 0}. The observed preferential ejections can be explained in terms of the theory of focusing-collision sequences. Because of the difference in surface binding energy between Al and Ni atoms, preferential ejection angles of Ni atoms are slightly different from those of Al atoms along the 〈1 1 1〉 ejections. For NiAl, the 〈1 1 1〉 preferential ejections were less prominent than the 〈1 0 0〉 preferential ejections and this can be related to the low efficiency of momentum transfer in Ni-Al collision sequences along 〈1 1 1〉 lattice directions. The low efficiency of momentum transfer due to the mass mismatch can also be responsible for the experimental observation that the preferential ejections in the alloy were less prominent than those in the Ni metal.     

Structure analysis of stoichiometric LiNbO3(0 0 0 1) surfaces using low-energy neutral scattering spectroscopy

The atomic structure of LiNbO₃(0 0 0 1) surface was investigated by low-energy neutral scattering spectroscopy (LENS). Poled stoichiometric LiNbO₃ (SLN) samples were prepared for the measurements. The LENS was developed for surface structure and composition analysis particularly of highly insulating materials and was successfully applied to the structure analysis of the SLN(0 0 0 1) surface. The polar angle dependences of intensity of scattered He⁰ from the poled SLN surfaces indicate obvious differences between the negatively and the positively charged surfaces. It is suggested that O atoms cover the surfaces, and the first metal layers underneath the O layer consist of Li and Nb for negatively and positively charged surfaces, respectively, parallel to the applied electric field.     

Formation of linear structures of Sr and Gd films on the Mo(1 1 2) and W(1 1 2) surfaces

At low coverages, alkaline earth and rare earth layers adsorbed on furrowed transition metal surfaces--such as W(1 1 2) and Mo(1 1 2)--tend to form commensurate linear structures built of monoatomic chains directed across the furrows. Presented Monte Carlo simulations provide insight into parameters of the long-range indirect interaction that, in competition with a dipole-dipole interaction, leads to formation of the linear structures. The most impressive result of the simulations is the revealed repulsive effective interaction between adjacent atoms in linear Sr chains on W(1 1 2) surface. In this case, formation of the linear chains is accomplished due to pronounced minima in the potential of indirect interaction along the surface furrows. This result predicts that a single Sr chain cannot exist on the W(1 1 2) surface.     

Study of buried Si(1 1 1)-5 × 2-Au by surface X-ray diffraction

The structure of buried Si(1 1 1)-5 × 2-Au capped with amorphous Si was investigated using surface X-ray diffraction. It was found that the 5 × 2 structural periodicity is kept under the amorphous Si from the in-plane measurement. Furthermore, the intensity variation along the fractional-order rod indicates that Au atoms are located almost on the same plane.     

The stability of fcc (1 1 1) metal surfaces

The results of a theoretical study on the stability of fcc (1 1 1) metal surfaces to certain commensurate-incommensurate reconstructive phase transformations is presented. Specifically, we have performed computer simulation studies of the 22×√3 surface reconstruction of Au(1 1 1). This reconstruction involves a uniaxial contraction of the top monolayer corresponding to a surface strain of about 4.3% and has been observed to be the stable structure for the clean surface at low temperatures. The driving force for the reconstruction has been associated with the quantity (f−γ), where f is the surface stress and γ is the surface free energy, while the opposing force is due to the disregistry with the underlying lattice. A continuum model yields a stability criterion that depends on the knowledge of a small number of physical quantities: f, γ, the equilibrium nearest-neighbor spacing r₁ and the shear modulus G. We have performed molecular dynamics simulations as a general stability analysis of these types of reconstructions. The results are in excellent agreement with the continuum model. The simulations using embedded-atom method potentials also accurately reproduce many observed features of the reconstruction on Au(1 1 1).     

Surface diffusion of carbon atom and carbon dimer on Si(0 0 1) surface

Carbon (C) atom and carbon dimer (C2) are known to be the main projectiles in the deposition of diamond-like carbon (DLC) films. The adsorption and diffusion of the C adatom and addimer (C2) on the fully relaxed Si(0 0 1)-(2 × 1) surface was studied by a combination of the molecular dynamics (MD) and Monte Carlo (MC) simulation. The adsorption sites of the C and C2 on the surface and the potential barriers between these sites were first determined using the semi-empirical many-body Brenner and Tersoff potential. We then estimated their hopping rates and traced their pathways. It is found that the diffusion of both C and C2 is strongly anisotropic in nature. In addition, the C adatom can diffuse a long distance on the surface while the adsorbed C2 is more likely to be confined in a local region. Thus we can expect that smoother films will be formed on the Si(0 0 1) surface with single C atoms as projectile at moderate temperature, while with C2 the films will grow in two-dimensional islands. In addition, relatively higher kinetic energy of the projectile, say, a few tens of eV, is needed to grow DLC films of higher quality. This is consistent with experimental findings.     

Active oxidation of Cu3Au(1 1 0) using hyperthermal O2 molecular beam

Oxidation of Cu₃Au(1 1 0) using a hyperthermal O₂ molecular beam (HOMB) was investigated by X-ray photoemission spectroscopy in conjunction with a synchrotron light source. From the incident energy dependence of the O-uptake curve, the precursor-mediated dissociative adsorption occurs, where the trapped O₂ molecule can migrate and dissociate at the lower activation-barrier sites, dominantly at thermal O₂ exposures. Dissociative adsorption of O₂ on Cu₃Au(1 1 0) is as effective at the thermal O₂ exposure as on Cu(1 1 0). On the other hand, at the incident energies of HOMB where the direct dissociative adsorption is dominant, it was determined that the dissociative adsorption of O₂ implies a higher activation barrier and therefore less reactivity due to the Au alloying in comparison with the HOMB oxidation of Cu(1 1 0). The dissociative adsorption progresses with the Cu segregation on Cu₃Au(1 1 0) similarly as on Cu₃Au(1 0 0). The growth of Cu₂O for 2 eV HOMB suggests that the diffusion of Cu atoms also contribute to the oxidation process through the open face, which makes the difference from Cu₃Au(1 0 0).     

Adsorption of gold on hydrogen terminated Si(0 0 1): Formation of chain structure

Possible formation of stable Au atomic wire on the hydrogen terminated Si(0 0 1): 3×1 surface is investigated under the density functional formalism. The hydrogen terminated Si(0 0 1): 3×1 surface is patterned in two different ways by removing selective hydrogen atoms from the surface. The adsorption of Au on such surfaces is studied at different sub-monolayer coverages. At 4/9 monolayer (ML) coverage, zigzag continuous Au chains are found to be stable on the patterned hydrogen terminated Si(0 0 1): 3×1 surface. The reason for the stability of the wire structures at 4/9 ML coverage is explained. It is to be noted that beyond 4/9 ML coverage, the additional Au atoms may introduce clusters on the surface. The continuous atomic gold chains on the substrate may be useful for the fabrication of atomic scale devices.     

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.     

Study of S ion-assisted sulfurization of n-GaAs (1 0 0) surface

The chemical structure and site location of sulfur atoms on n-GaAs (1 0 0) surface treated by bombardment of S⁺ ions over their energy range from 10 to 100 eV have been studied by X-ray photoelectron spectroscopy and low energy electron diffraction. The formation of Ga-S and As-S species on the S⁺ ion bombarded n-GaAs surface is observed. An apparent donor doping effect is observed for the n-GaAs by the 100 eV S⁺ ion bombardment. It is found that the S⁺ ions with higher energy are more effective in the formation of Ga-S species, which assists the n-GaAs (1 0 0) surface in reconstruction into an ordered (1 × 1) structure upon subsequent annealing. The treatment is further extended to repair Ar⁺ ion damaged n-GaAs (1 0 0) surface. It is found that after a n-GaAs (1 0 0) sample is damaged by 150 eV Ar⁺ ion bombardment, and followed by 50 eV S⁺ ion treatment and subsequent annealing process, finally an (1 × 1) ordering GaAs (1 0 0) surface with low surface states is obtained.     

Molecular dynamics investigation on the atomic-scale friction behaviors between copper(0 0 1) and diamond(1 1 1) surfaces

Classical molecular dynamics (MD) simulations are conducted to examine the atomic-scale friction behavior of an infinite flat-flat contact between copper(0 0 1) and diamond(1 1 1) surfaces. Two types of diamond surface, namely H-free and hydrogenated, are constructed and on each of them the copper counterface is brought to slide along the [1 1 −2] and [1 −1 0] crystallographic directions with a variety of loads. The simulation results demonstrate that the hydrogen atoms chemisorbed to the diamond surface can to large extent eliminate the directional dependency of its friction behavior with copper. Under pressures less than 30 GPa, the sliding between copper and hydrogenated is wearless. In this period, the shear stress of them just slightly increases to 0.6 GPa. Between 30 GPa and 32 GPa, copper atoms near the interface begin to be worn and incorporate into the diamond substrate and this causes a sharp shift from 0.6 GPa to 2.7 GPa in their shear stress. In contrast, the sliding process between copper and H-free diamond is always wearless even under pressure beyond 40 GPa. The H-free [1 −1 0] model exhibits much higher shear stress than H-free [1 1 −2] under pressures less than 35 GPa. Beyond 35 GPa, they present nearly consistent shear stress evolution. Moreover, the simulations for hydrogenated diamond models suggest that their friction behavior is independent on sliding velocity only under wearless sliding regime.