Characterization of diamond (1 0 0) surface with oxygen termination

Ab initio density functional theory (DFT) was employed to study reconstructions of diamond (1 0 0) surfaces in the presence of hydrogen, oxygen and hydroxyl. Clean and (2 × 1):1H surfaces are taken as reference. The properties of oxidization diamond surfaces with several adsorption structures, namely, O-on-top (OT) site, O-bridge (BR) site, hydroxyl (-OH), hydroxyl/hydroxyl, OT/hydroxyl, BR/hydroxyl have been considered. The calculated results indicate that the BR model is much more stable than the OT model, and the most energetically favorable structures of oxygenated surfaces are those with chemisorbed hydroxyl (-OH) group. Furthermore, the stability of the structures is also discussed from the point of HOMO-LUMO gap. Analysis of electronic structures shows that the presence of hydrogen induces surface conductivity whereas oxygen weakens it.     

Study of ammonia molecule adsorbing on diamond (1 0 0) surface

The diamond (1 0 0) surface with amino terminations is investigated based on density function theory within the generalized gradient approximation. Our calculated negative electron affinity of diamond (1 0 0) surface with hydrogen termination provides a necessary condition for initiating radical reaction. The results display that the ammonia molecule can form stable C-N covalent bonds on the diamond surface. In addition, due to the lower adsorption energy of one amino group binding on diamond surface, single amino group (SAG) model is easy to be realized in experiment with the comparison of double amino group (DAG) model. The adsorbed ammonia molecule will induce acceptor-like gap states with little change of the valence and conduction band of diamond in SAG model. The adsorption mechanism in the formation of ammonia monolayer on H-terminated diamond (1 0 0) surface, and two possible adsorption structures (SAG and DAG) were especially studied.     

Olefin bonding mechanism and growth to (1 0 0)(2 × 1):H diamond

The bonding and growth mechanism of photochemically attached olefin molecules to (1 0 0)(2 × 1):H diamond is characterized using atomic force (AFM) and scanning tunneling microscopy (STM) experiments in combination with molecular orbital calculations. To identify growth schemas, diamond surfaces after 10, 40 and 90 min of photo-chemically stimulated growth have been characterized. These data show clearly island formation which is discussed taking into account a growth model from silicon. The island growth shows no directional properties which are attributed to arrangement and geometrical properties of hydrogen terminated carbon bonds at the surface of (1 0 0) oriented (2 × 1) reconstructed diamond.     

Kinetic Monte Carlo simulations of Pd deposition and island growth on MgO(1 0 0)

The deposition and ripening of Pd atoms on the MgO(1 0 0) surface are modeled using kinetic Monte Carlo simulations. The density of Pd islands is obtained by simulating the deposition of 0.1 ML in 3 min. Two sets of kinetic parameters are tested and compared with experiment over a 200-800 K temperature range. One model is based upon parameters obtained by fitting rate equations to experimental data and assuming the Pd monomer is the only diffusing species. The other is based upon transition rates obtained from density functional theory calculations which show that small Pd clusters are also mobile. In both models, oxygen vacancy defects on the MgO surface provide strong traps for Pd monomers and serve as nucleation sites for islands. Kinetic Monte Carlo simulations show that both models reproduce the experimentally observed island density versus temperature, despite large differences in the energetics and different diffusion mechanisms. The low temperature Pd island formation at defects is attributed to fast monomer diffusion to defects in the rate-equation-based model, whereas in the DFT-based model, small clusters form already on terraces and diffuse to defects. In the DFT-based model, the strong dimer and trimer binding energies at charged oxygen vacancy defects prevent island ripening below the experimentally observed onset temperature of 600 K.     

Characteristic energy band values and electron attenuation length of a chemical-vapor-deposition diamond (0 0 1)2 × 1 surface

The characteristic energy band values such as the Fermi-level position with respect to valence band top for a boron-doped p-type hydrogen-terminated chemical-vapor-deposition (CVD) diamond (0 0 1)2 × 1 surface and for a clean CVD diamond (0 0 1)2 × 1 surface have been determined by a new method with an accuracy of ±0.02 eV. The electron attenuation length for the clean diamond (0 0 1)2 × 1 surface for the electron kinetic energy of C 1s X-ray photoemission peak by Mg Kα excitation is experimentally determined to be 2.1-2.2 nm. These values are compared and discussed with the previously reported experimental and simulation values.     

Electron-stimulated desorption from an unexpected source: Internal hot electrons for Br-Si(1 0 0)-(2 × 1)

The desorption of Br adatoms from Br-saturated Si(1 0 0)-(2 × 1) was studied with scanning tunneling microscopy as a function of dopant type, dopant concentration, and temperature for 620-775 K. Analysis yields the activation energies and prefactors for desorption, and the former correspond to the energy separation between the Fermi level and Si-Br antibonding states. Thus, electron capture in long-lived states results in Br expulsion via a Franck-Condon transition. Analysis of the prefactors reveals that optical phonons provide the energy needed for the electronic excitation. These results show that desorption induced by an electronic transition can occur in closed system without external stimulus, and they indicate that thermally-excited charge carriers may play a general role in surface reactions.     

Early stages of vanadium deposition on Si(1 1 1)-7 × 7

We investigate the low-coverage regime of vanadium deposition on the Si(1 1 1)-7 × 7 surface using a combination of scanning tunnelling microscopy (STM) and density-functional theory (DFT) adsorption energy calculations. We theoretically identify the most stable structures in this system: (i) substitutional vanadium atoms at silicon adatom positions; (ii) interstitial vanadium atoms between silicon adatoms and rest atoms; and (iii) interstitial vanadium - silicon adatom vacancy complexes. STM images reveal two simple vanadium-related features near the Si adatom positions: bright spots at both polarities (BB) and dark spots for empty and bright spots for filled states (DB). We relate the BB spots to the interstitial structures and the DB spots to substitutional structures.     

Comparative investigation of the c(2 × 2)-Si/Cu(0 1 1) and (√3 × √3)R30°-Cu2Si/Cu(1 1 1) surface alloys using DFT

The electronic structure of the c(2 × 2)-Si/Cu(0 1 1) surface alloy has been investigated and compared to the structures seen in the three phases of the (√3 × √3)R30°Cu₂Si/Cu(1 1 1) system, using LCAO-DFT. The weighted surface energy increase between the alloyed Cu(0 1 1) and Cu(1 1 1) surfaces is 126.7 meV/Si atom. This increase in energy for the (0 1 1) system when compared to the (1 1 1) system is assigned to the transition from a hexagonal to a rectangular local bonding environment for the Si ion cores, with the hexagonal environment being energetically more favorable. The Si 3s state is shown to interact covalently with the Cu 4s and 4p states whereas the Si 3p state, and to a lesser extent the Si 3d state, forms a mixture of covalent and metallic bonds with the Cu states. The Cu 4s and 4p states are shown to be altered by approximately the same amount by both the removal of Cu ion cores and the inclusion of Si ion cores during the alloying of the Cu(0 1 1) surface. However, the Cu 3d states in the surface and second layers of the alloy are shown to be more significantly altered during the alloying process by the removal of Cu ion cores from the surface layer rather than by the addition of Si ion cores. This is compared to the behavior of the Cu 3d states in the surface and second layers of the each phase of the (√3 × √3)R30°-Cu₂Si/Cu(1 1 1) alloy and consequently the loss of Cu-Cu periodicity during alloying of the Cu(0 1 1) surface is conjectured as the driving force for changes to the Cu 3d states. The accompanying changes to the Cu 4s and 4p states in both the c(2 × 2)-Si/Cu(0 1 1) and (√3 × √3)R30°-Cu₂Si/Cu(1 1 1) alloys are quantified and compared. The study concludes with a brief quantitative study of changes in the bond order of the Cu-Cu bonds during alloying of both Cu(0 1 1) and Cu(1 1 1) surfaces.     

Structural analysis of the c(4 × 2) reconstruction in Si(0 0 1) and Ge(0 0 1) surfaces by low-energy electron diffraction

The c(4 × 2) structures in (0 0 1) surfaces of Si and Ge have been studied by low-energy electron diffraction (LEED). Using a proper cleaning method for the Si surface, we were able to observe clear c(4 × 2) LEED patterns up to incident energy of ∼400 eV as well as the Ge surface. Extensive experimental intensity-voltage curves allowed us to optimize the asymmetric dimer model up to the eighth layer (including the dimer layer) in depth in the dynamical LEED calculation. Optimized structural parameters are almost the same for the Si and Ge except for the height of the buckled-up atom of the asymmetric dimer. For the Ge surface, the structural parameters are in excellent agreement with those obtained by a previous theoretical calculation. The tilt angle and bond length of the dimer are 18 ± 1 (19 ± 1)° and 2.4 ± 0.1 (2.5 ± 0.1) Å for the Si(0 0 1) (Ge(0 0 1)), respectively.     

Vacancy charging on Si(1 1 1)-(7 × 7) investigated by density functional theory

The structure and energetics of charged vacancies on Si(1 1 1)-(7 × 7) are investigated using density functional theory calculations supplemented by estimates of ionization entropy. The calculations predict multiple possible charge states for the unfaulted edge vacancy in the adatom layer, although the −2 state is most stable on real Si(1 1 1) surfaces for which the Fermi level lies near the middle of the band gap.     

Density functional theory study of selenium adsorption on Fe(1 1 0)

The adsorption of atomic Se on a Fe(1 1 0) surface is examined using the density functional theory (DFT). Selenium is adsorbed in high-symmetry adsorption sites: the -short and long-bridge, and atop sites at 1/2, 1/4, and 1 monolayer (ML) coverages. The long bridge (LB) site is found to be the most stable, followed by the short bridge (SB) and top sites (T). The following overlayer structures were examined, p(2 × 2), c(2 × 2), and p(1 × 1), which correspond to 1/4 ML, 1/2 ML, and 1 ML respectively. Adsorption energy is −5.23 eV at 1/4 ML. Se adsorption results in surface reconstruction, being more extensive for adsorption in the long bridge site at 1/2 ML, with vertical displacements between +8.63 and −6.69% -with regard to the original Fe position-, affecting the 1st and 2nd neighbours. The largest displacement in x or y-directions was determined to be 0.011, 0.030, and 0.021 Å for atop and bridge sites. Comparisons between Se-adsorbed and pure Fe surfaces revealed reductions in the magnetic moments of surface-layer Fe atoms in the vicinity of the Se. At the long bridge site, the presence of Se causes a decrease in the surface Fe d-orbital density of states between 4 and 5 eV below Fermi level. The density of states present a contribution of Se states at −3.1 eV and −12.9 eV. stabilized after adsorption. The Fe-Fe overlap population decrease and a Fe-Se bond are formed at the expense of the metallic bond.     

Growth of (√3 × √3)-Ag and (1 1 1) oriented Ag islands on Ge/Si(1 1 1) surfaces

We have studied the growth of Ag on Ge/Si(1 1 1) substrates. The Ge/Si(1 1 1) substrates were prepared by depositing one monolayer (ML) of Ge on Si(1 1 1)-(7 × 7) surfaces. Following Ge deposition the reflection high energy electron diffraction (RHEED) pattern changed to a (1 × 1) pattern. Ge as well as Ag deposition was carried out at 550 °C. Ag deposition on Ge/Si(1 1 1) substrates up to 10 ML has shown a prominent (√3 × √3)-R30° RHEED pattern along with a streak structure from Ag(1 1 1) surface. Scanning electron microscopy (SEM) shows the formation of Ag islands along with a large fraction of open area, which presumably has the Ag-induced (√3 × √3)-R30° structure on the Ge/Si(1 1 1) surface. X-ray diffraction (XRD) experiments show the presence of only (1 1 1) peak of Ag indicating epitaxial growth of Ag on Ge/Si(1 1 1) surfaces. The possibility of growing a strain-tuned (tensile to compressive) Ag(1 1 1) layer on Ge/Si(1 1 1) substrates is discussed.     

Ab initio study of atomic geometry and electronic states of GaSb(0 0 1) reconstructions

We have employed the pseudopotential method and the density functional scheme to study the atomic geometry and electronic states of the GaSb(0 0 1) surface such as (1 × 3), c(2 × 6) and (4 × 3) reconstructions. It is found that both of (1 × 3) and c(2 × 6) reconstructions are characterised by metallic band structures, and thus violate the so-called electron counting rule, one of the main building principles of the stability of compound semiconductor surfaces. We establish that the stability of these reconstructions results from significant elastic deformation in the top atomic layers of the surface, a process which overcomes the penality incurred by the violation of the electron counting rule. The atomic geometry and electronic states for the two reconstructions are compared and contrasted with each other. The α and β phases of the (4 × 3) reconstruction also show large atomic relaxations but are semiconducting and obey the electron counting rule.     

Halogen chemisorption, the pairwise diffusion of I, and trapping by defects on Si(1 0 0)

Halogen molecules dissociatively chemisorb on Si(1 0 0)-(2 × 1), and the bonding structures that they adopt can be elucidated with scanning tunneling microscopy. Of the Cl, Br, and I group, Cl has the highest single atom diffusion barrier, and both single and paired adatoms are observed at 295 K. The barrier is smaller for Br, and the adatoms can interrogate the surface until they form pairs, which are then immobile, or are trapped at C-type defects. The barrier is smallest for I, allowing the formation of pairs and trapped states, but the pairs are mobile at ambient temperature. Their motion is thermally activated, the events are random, and the diffusivities along and across the dimer row are ∼0.42 and ∼0.17 Å²/s at 295 K. The respective energy barriers for pairwise diffusion are ∼0.76 and ∼0.82 eV, assuming an attempt frequency of 10¹² s⁻¹. Studies over long times reveal that pairwise diffusion at low coverage is ultimately quenched by the increasing density of C-type defects, i.e. the increasing amounts of dissociated H₂O.     

Effects of electronic confinement and substrate on the low-temperature growth of Pb islands on Si(1 0 0)-2 × 1 surfaces

The growth of Pb films on the Si(1 0 0)-2 × 1 surface has been investigated at low temperature using scanning tunneling microscopy. Although the orientation of the substrate is (1 0 0), flat-top Pb islands with (1 1 1) surface can be observed. The island thickness is confined within four to nine atomic layers at low coverage. Among these islands, those with a thickness of six layers are most abundant. Quantum-well states in Pb(1 1 1) islands of different thickness are acquired by scanning tunneling spectroscopy. They are found to be identical to those taken on the Pb(1 1 1) islands grown on the Si(1 1 1)7 × 7 surface. Besides Pb(1 1 1) islands, two additional types of Pb islands are formed: rectangular flat-top Pb(1 0 0) islands and rectangular three-dimensional (3D) Pb islands, and both their orientations rotate by 90° from a terrace to the adjacent one. This phenomenon implies that the structures of Pb(1 0 0) and 3D islands are influenced by the Si(1 0 0)-2 × 1 substrate.     

Layer-resolved photoelectron diffraction from Si(0 0 1) and GaAs(0 0 1)

Photoelectron diffraction in the layer-resolved mode brings more detailed information about local atomic arrangement than is obtained in the standard mode. This is demonstrated in crystals with diamond and zinc-blende structures, both for unpolarized photon excitation as well as for circularly polarized excitation. The full angular distributions of photoemission intensities are evaluated for large atomic clusters representing ideally truncated surfaces of Si(0 0 1) and GaAs(0 0 1). Highly structured layer-resolved patterns enable a more detailed understanding of the standard mode outcomes. Photoelectron intensities from atomic layers placed at different depths under the crystal surface provide direct evidence about electron attenuation and its anisotropy in crystals.     

Growth and morphology of the epitaxial Fe(1 1 0)/MgO(1 1 1)/Fe(1 1 0) Trilayers

Growth and surface morphology of epitaxial Fe(1 1 0)/MgO(1 1 1)/Fe(1 1 0) trilayers constituting a magnetic tunnel junction were investigated by low-energy electron diffraction (LEED) and scanning tunneling microscopy (STM). STM reveals a grain-like growth mode of MgO on Fe(1 1 0) resulting in dense MgO(1 1 1) films at room temperature as well as at 250 °C. As observed by STM, initial deposition of MgO leads to a partial oxidation of the Fe(1 1 0) surface which is confirmed by Auger electron spectroscopy. The top Fe layer deposited on MgO(1 1 1) at room temperature is relatively rough consisting of clusters which can be transformed by annealing to an atomically flat epitaxial Fe(1 1 0) film.     

A density-functional study on the atomic geometry and adsorption of the Cu(1 0 0) c(2 × 2)/N Surface

In this work we have performed total-energy calculations on the geometric structure and adsorption properties of Cu(1 0 0) c(2 × 2)/N surface by using the density-functional theory and the projector-augmented wave method. It is concluded that nitrogen atom was adsorbed on a FFH site with a vertical distance of 0.2 Å towards from surface Cu layer. The bond length of the shortest Cu-N bonding is calculated to be 1.83 Å. Geometry optimization calculations exclude out the possibilities of adsorbate induced reconstruction mode suggested by Driver and Woodruff and the atop structural model. The calculated workfunction for this absorbate-adsorbent system is 4.63 eV which is quite close to that of a clean Cu(1 0 0) surface. The total-energy calculations showed that the average adsorption energy per nitrogen in the case of Cu(1 0 0) c(2 × 2)-N is about 4.88 eV with respect to an isolated N atom. The absorption of nitrogen on Cu(1 0 0) surface yields the hybridization between surface Cu atoms and N, and generates the localized surface states at −1.0 eV relative to Fermi energy E_F. The stretch mode of the adsorbed nitrogen at FFH site is about 30.8 meV. The present study provides a strong criterion to account for the local surface geometry in Cu(1 0 0) c(2 × 2)/N surface.     

First-principles calculations on electronic structures of Fe-vacancy-codoped TiO2 anatase (1 0 1) surface

The electronic structures of Fe-doped TiO₂ anatase (1 0 1) surfaces have been investigated by all spin-polarized density functional theory (DFT) plane-wave pseudopotential method. The general gradient approximation (GGA)+U (Hubbard coefficient) method has been adopted to describe the exchange-correlation effects. Through the density functional calculations for the formation energies of various configurations, the complex of a substitutional Fe plus an O vacancy was found to form easily in the most range of O chemical potential. The calculated density of the states of the system of Fe-doped surface with a surface oxygen vacancy shows a band gap narrowing from 2.8 to 1.9 eV comparing with the pure surface due to the synergistic effects of surface Fe impurities with O vacancies. The system processes high visible light sensitivity and photocatalytic ability by decreasing extrinsic absorption energy. By comparing the partial DOS of some O and Ti atoms lying in the outermost and bottom layers of Fe-doped surfaces, it was found that the influence of Fe impurities on the electronic structure of the system is localized.     

Growth kinetics of metastable (3 3 1) nanofacet on Au and Pt(1 1 0) surfaces

A theoretical epitaxial growth model with realistic barriers for surface diffusion is investigated by means of kinetic Monte Carlo simulations to study the growth modes of metastable (3 3 1) nanofacets on Au and Pt(1 1 0) surfaces. The results show that under experimental atomic fluxes, the (3 3 1) nanofacets grow by 2D nucleation at low temperature in the submonolayer regime. A metastable growth phase diagram that can be useful to experimentalists is presented and looks similar to the one found for the stationary growth of the bcc(0 0 1) surface in the kinetic 6-vertex model.