Effect of hydrogenation on P/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 energetics of substitutional phosphorus (P) on the generic Si(0 0 1)-(1 × 2) surface. For the 0.5 ML coverage of P, we have considered three different substitutional sites: (i) the mixed Si-P dimer structure (i.e., the P-nondiffused case), (ii) P-interdiffused to the second layer Si (i.e., intermixed P-Si bond structure) and (iii) P-interdiffused to the third layer Si. We have found that the mixed Si-P dimer structure is 0.79 eV/dimer energetically more favorable than the P-interdiffused case. However, for the hydrogenation of above cases, we have found that the situation is reversed and the interdiffused case is 0.3 eV/dimer energetically more favorable than the P-nondiffused case. Reductions in the number of P-Si is identified as a contributing factor which determines energetically the stable structures during P on Si(0 0 1).     

Mixed PbSi dimer chains on Si(1 0 0): a first-principles study

It has been a common belief that the one-dimensional structures observed by STM at low coverage of Pb on Si(1 0 0) are buckled Pb-Pb dimer chains. However, using first-principles density functional calculations, we found that it is energetically favorable for Pb adatoms to intermix with Si atoms to form mixed dimer chains on Si(1 0 0), instead of Pb-Pb dimer chains as assumed in previous studies. Up to a Pb coverage of 0.125 ML, mixed PbSi dimer chain is 0.19 eV per Pb atom lower in energy than Pb dimer chain.     

The interaction of carbon with Si(1 1 1):As and Si(1 1 1):H surfaces, a theoretical study

Density functional theory calculations have been performed to determine the adsorption site of carbon at the Si(1 1 1):As and Si(1 1 1):H surfaces at different coverages. The As- and H-passivated surfaces were simulated by replacing the topmost Si layer by As or by saturating the Si dangling bonds with hydrogen atoms, respectively. Different high symmetry sites were considered. Carbon was placed successively in the fourfold (T₄) or threefold coordinated (H₃), the ontop (T₁) sites or substituted for a Si atom in the S₅ position located underneath the Si adatom in the T₄ site. We found that the preferred carbon adsorption site depends on the coverage of the passivated surfaces. At low coverages i.e. at 1/16 ML and 1/3 ML, it prefers a distorted T₄ position whereas at 1 ML, it occupies an H₃ site. This contrasts with the clean surface where the most energetically favored site is the S₅ at all coverages. Carbon adsorption induces a significant change in the structural geometry of the surface atoms, leading to a charge re-arrangement in the surface layers.     

Atomic and electronic structure of Sr/Si(0 0 1)-(2 × 2)

The adsorption of Sr on the Si(0 0 1) surface with the semiantiphase dimer (2 × 2) reconstruction is studied, based upon the ab initio pseudopotential calculations. It is calculated that the semiantiphase dimer (2 × 2) reconstruction (2 dimers per unit cell) is more favorable than the (2 × 1) phase (1 dimer per unit cell) by an energy of about 0.24 eV/dimer. Considering the energetically more stable reconstruction, we have assumed four possible locations for 1/4 monolayer (ML) Sr adsorption on this surface: (i) bridge, (ii) cave, (iii) pedestal, and (iv) valley-bridge. We find that Sr adsorption on the valley-bridge site is energetically more favorable than all other cases studied here. Interestingly, one of the dimers becomes symmetric, but the other one is still asymmetric with the buckling angle reduced from 18° to 14°, when compared with the clean Si(0 0 1)-(2 × 2) surface. The calculated bond length between Sr and Si in the case of valley-bridge adsorption site is 3.05 Å, and in good agreement with other theoretical calculations. We also present and compare the electronic band structures for the clean and covered surfaces as well as the corresponding charge density plots.     

Effect of hydrogenation on the electronic structure of the P/Si(0 0 1)-(1 × 2) surface

Ab initio calculations, based on pseudopotentials and density functional theory (DFT), have been performed to investigate the effect of hydrogenation on the electronic properties of P/Si(0 0 1)-(1 × 2) surface. In parallel with this, the electronic band structure of the hydrogenated and non-hydrogenated P/Si(0 0 1)-(1 × 2) surface have been calculated for half- and full-monolayer P. For the mixed Si-P dimer structure, we have identified two occupied and one unoccupied surface state, which correspond to 0.5 ML coverage of P. When this surface is terminated with H, we see that two occupied states completely disappeared and that one unoccupied state is shifted towards the conduction band. A similar calculations for the 1 ML coverage of P have been also carried out. It is seen that the unoccupied state C₁ appeared in the P/Si(0 0 1)-(1 × 2) surface is passivated when this surface is terminated with the H atoms. To explain the nature of the surface states, we have also plotted the total and partial charge densities at the point of the Surface Brillouin Zone (SBZ).     

Germanium growth on Br-terminated Si(1 0 0)

The consequences of Ge deposition on Br-terminated Si(1 0 0) were studied with scanning tunneling microscopy at ambient temperature after annealing at 650 K. One monolayer of Br was sufficient to prevent the formation of Ge huts beyond the critical thickness of 3 ML. This is possible because Br acts as a surfactant whose presence lowered the diffusivity of Ge adatoms. Hindered mobility was manifest at low coverage through the formation of short Ge chains. Further deposition resulted in the extension and connection of the Ge chains and gave rise to the buildup of incomplete layers. The deposition of 7 ML of Ge resulted in a rough surface characterized by irregularly shaped clusters. A short 800 K anneal desorbed the Br and allowed Ge atoms to reorganize into the more energetically favorable “hut” structures produced by conventional Ge overlayer growth on Si(1 0 0).     

Progressive changes in atomic structure and gap states on Si(0 0 1) by Bi adsorption

From ab initio studies employing the pseudopotential method and the density functional scheme, we report on progressive changes in geometry, electronic states, and atomic orbitals on Si(0 0 1) by adsorption of different amounts of Bi coverage. For the 1/4 ML coverage, uncovered Si dimers retain the characteristic asymmetric (tilted) geometry of the clean Si(0 0 1) surface and the Si dimers underneath the Bi dimer have become symmetric (untilted) and elongated. For this geometry, occupied as well as unoccupied surface states are found to lie in the silicon band gap, both sets originating mainly from the uncovered and tilted silicon dimers. For the 1/2 ML coverage, there are still both occupied and unoccupied surface states in the band gap. The highest occupied state originates from an elaborate mixture of the p_z orbital at the Si and Bi dimer atoms, and the lowest unoccupied state has a ppσ* antibonding character derived from the Bi dimer atoms. For 1 ML coverage, there are no surface states in the fundamental bulk band gap. The highest occupied and the lowest unoccupied states, lying close to band edges, show a linear combination of the p_z orbitals and ppσ* antibonding orbital characters, respectively, derived from the Bi dimer atoms.     

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.     

A theoretical study of Fe adsorption along Bi-nanolines on the H/Si(0 0 1) surface

We have investigated the energetic stability and equilibrium geometry of the adsorption of transition metal Fe atoms near the self-organized Bi lines on hydrogen passivated Si(0 0 1) surface. Our total energy results show that there is an attractive interaction between Fe adatoms along the Bi-nanolines. For the energetically most stable configuration, the Fe adatoms are seven-fold coordinated, occupying the subsurface interstitial sites aside the Bi-nanolines. With increased coverage, Fe atoms are predicted to form two parallel lines, symmetrically on both sides of the Bi line. Within our local spin-density functional calculations, we find that for the most stable geometries the Fe adatoms exhibit an antiferromagnetic coupling.     

Diffusion of Pb adatom and ad-dimer on Si(1 0 0) from ab initio studies

The diffusion pathways of Pb adatoms and ad-dimers on Si(1 0 0) are investigated by first-principles calculations. Pb adatoms are found to diffuse on top of the Si(1 0 0) dimer row with an energy barrier of 0.31 eV. However, Pb dimers are energetically more stable. Pb dimers on top of the dimer row have a high energy barrier (0.95 eV) to rotate from the lowest energy configuration to the orientation parallel to the underlying Si(1 0 0) dimer row. Once the ad-dimer is oriented parallel to Si(1 0 0) dimer row, they can diffuse along the dimer row with an energy barrier of only 0.32 eV.     

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.     

First-principles study of indium on silicon (1 0 0)—the structure, defects and interdiffusion

The atomic structures of indium (In) on silicon (Si) (1 0 0)-(2 × 1) surface are investigated by the local density approximation using first-principles pseudopotentials. Total energy optimizations show that the energetically favored structure is the parallel ad-dimer model. The adsorption energy of In on ideal Si(1 0 0)-(1 × 1) surface is significantly higher than that on reconstructed Si(1 0 0)-(2 × 1) surface, suggesting that In adsorption does not break the Si-Si dimer bond of the substrate. When Si surface contains single dimer vacancy defects, In chain will be interrupted, leading to disconnected In nanowires. Displacive adsorption of In on Si(1 0 0) is also considered, and the calculation suggests that interdiffusion of In into Si substrate will not be favorable under equilibrium conditions.     

Spectroscopic interpretation for Na-induced phase transitions on Na/Si(1 1 1)3 × 1

Na adsorption at room temperature causes the Na/Si(1 1 1)3 × 1 surface with Na coverage of 1/3 monolayer (ML) to transit into the Na/Si(1 1 1)6 × 1 surface at 1/2 ML and sequentially into the Na/Si(1 1 1)3 × 1 surface at 2/3 ML. The phase transition was studied by Si 2p core-level photoemission spectroscopy. The detailed line shape analysis of the Si 2p core-level spectrum of the Na/Si(1 1 1)3 × 1 surface (2/3 ML) is presented and compared to the Na/Si(1 1 1)3 × 1 surface (1/3 ML) which is composed of Si honeycomb chain-channel structures. This suggests that as additional Na atoms form atomic chains resulting in the Na/Si(1 1 1)3 × 1 surface (2/3 ML), the inner atoms of the Si honeycomb chain-channel structure is buckled due to the additional Na atoms.     

Selforganized nanopatterning of Si(0 0 1)

The (2 × n) superstructure of Si(0 0 1) consists of elongated (2 × 1) reconstructed stripes separated by a dimer-vacancy line every few nanometers, thus offering a means to obtain a nanopattered Si(0 0 1) surface. Scanning tunneling microscopy (STM) investigations of Si(0 0 1) substrates that were deoxidized at 880-920 °C reveal that the formation of the (2 × n) depends strongly on the Si coverage of the topmost surface layer. It forms only in a narrow coverage window ranging from 0.6 to 0.8 ML. Systematic Monte Carlo simulations by an algorithm that combines the diffusion of monomers and dimers with the simultaneous deposition of Si onto the Si(0 0 1) surface, corroborate the STM results and suggest Si deposition as a viable alternative for introducing dimer vacancies in a well-defined manner.     

Chemisorption of 4-(4-amino-phenylazo) benzoic acid molecule on the Si(0 0 1)-(4 × 2) surface

Structural and electronic properties of self-assembled monolayer with 4-(4-amino-phenylazo) benzoic acid (APABA) on the Si(0 0 1)-(4 × 2) surface are investigated by ab initio calculation based on density functional theory. For the APABA chemisorption on the silicon surface, we have assumed two different binding sites: (i) amino group of molecule and (ii) carboxyl group of molecule. Considering amino-site, we have assumed two possible models for the chemisorption of molecules on the Si(0 0 1)-(4 × 2) surface: (i) an intrarow position between two neighboring Si dimers in the same dimer row (Model I), (ii) on-dimer position (Model II). We have found that Model II is 1.10 eV energetically more favorable than Model I. The Si-N bond length was calculated as 1.85 Å which is in excellent agreement with the sum of the corresponding covalent radii of 1.87 Å. Considering carboxyl-site, we have assumed exactly the same model as mentioned above. Again we have found that Model II is energetically favorable than Model I. The calculated bond lengths for Si-O and O-C are 1.76 and 1.35 Å, respectively.     

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.     

Study on the adsorption of fluorescein on Ag(1 1 0) substrate

The adsorption of fluorescein on the Ag(1 1 0) surface has been investigated by the first-principles pseudopotential method. Various adsorption geometries have been calculated and the energetically most favorable structure of fluorescein/Ag(1 1 0) was identified. The fluorescein molecule, in most favorable structure, is on hollow site, and the adsorption energy is 2.34 eV. Here the adsorption sites refer to the positions at the first layer of the substrate where the middle carbon atom of the fluorescein molecule is located. The bonding strength of the fluorescein molecule to the Ag substrate is site selective, being determined by electron transfer to the oxygen atoms of the molecule and local electrostatic attraction between the oxygen atoms and the silver atoms.     

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.     

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.     

Methanol adsorption on silicon (0 0 1)

Using a first-principles pseudopotential technique, we have investigated the adsorption of CH₃OH on the Si(0 0 1) surface. We have found that, in agreement with the overall experimental picture, the most probable chemisorption path for methanol adsorption on silicon (0 0 1) is as follows: the gas phase CH₃OH adsorbs molecularly to the electrophilic surface Si atom via the oxygen atom and then dissociates into Si-OCH₃ and H, bonded to the electrophilic and nucleophilic surface silicon dimer atoms, respectively. Other possible adsorption models and dissociation paths are also discussed. Our calculations also suggest that the most probable methanol coverage is 0.5 ML, i.e., one molecule per Si-Si dimer, in agreement with experimental evidences. The surface atomic and electronic structures are discussed and compared to available theoretical and experimental data. In addition, we propose that a comparison of our theoretical STM images and calculated vibrational modes for the adsorbed systems with detailed experimental investigations could possibly confirm the presented adsorption picture.