SiO adsorption on a p(2 × 2) reconstructed Si(1 0 0) surface

We have investigated the adsorption mechanism of SiO molecule incident on a clean Si(1 0 0) p(2 × 2) reconstructed surface using density functional theory based methods. Stable adsorption geometries of SiO on Si surface, as well as their corresponding activation and adsorption energies are identified. We found that the SiO molecule is adsorbed on the Si(1 0 0) surface with almost no activation energy. An adsorption configuration where the SiO binds on the channel separating the dimer rows, forming a Si-O-Si bridge on the surface, is the energetically most favourable geometry found. A substantial red-shift in the calculated vibrational frequencies of the adsorbed SiO molecule in the bridging configurations is observed. Comparison of adsorption energies shows that SiO adsorption on a Si(1 0 0) surface is energetically less favourable than the comparable O₂ adsorption. However, the role of SiO in the growth of silicon sub-oxides during reactive magnetron plasma deposition is expected to be significant due to the relatively large amount of SiO molecules incident on the deposition surface and its considerable sticking probability. The stable adsorption geometries found here exhibit structural properties similar to the Si/SiO₂ interface and may be used for studying SiO_x growth.     

An ab initio-based approach to adsorption-desorption behavior of Si adatoms on GaAs(1 1 1)A-(2 × 2) surfaces

The adsorption-desorption behavior of Si adatoms on GaAs(1 1 1)A-(2 × 2) surfaces is investigated using our ab initio-based approach, in which adsorption and desorption behavior of Si adatoms is described by comparing the calculated desorption energy obtained by total-energy electronic-structure calculations with the chemical potential estimated by quantum statistical mechanics. We find that the Si adsorption at the Ga-vacancy site on the (2 × 2) surfaces with As adatoms occurs less than 1140-1590 K while the adsorption without As adatom does less than 630-900 K. The change in adsorption temperature of Si adatoms by As adatoms is due to self-surfactant effects of As adatoms: the promotion of the Si adsorption triggered by As adatoms is found to be interpreted in terms of the band-energy stabilization. Furthermore, the stable temperature range for Si adsorbed surfaces with As adatoms agrees with the experimental results. The obtained results provide a firm theoretical framework to clarify n-type doping processes during GaAs epitaxial growth.     

Ga-induced superstructures on the Si(1 1 1) 7 × 7 surface

Monolayer Ga adsorption on Si surfaces has been studied with the aim of forming p-delta doped nanostructures. Ga surface phases on Si can be nitrided by N₂⁺ ion bombardment to form GaN nanostructures with exotic electron confinement properties for novel optoelectronic devices. In this study, we report the adsorption of Ga in the submonolayer regime on 7 × 7 reconstructed Si(1 1 1) surface at room temperature, under controlled ultrahigh vacuum conditions. We use in-situ Auger electron spectroscopy, electron energy loss spectroscopy and low energy electron diffraction to monitor the growth and determine the properties. We observe that Ga grows in the Stranski-Krastanov growth mode, where islands begin to form on two flat monolayers. The variation in the dangling bond density is observed during the interface evolution by monitoring the Si (LVV) line shape. The Ga adsorbed system is subjected to thermal annealing and the residual thermal desorption studied. The difference in the adsorption kinetics and desorption dynamics on the surface morphology is explained in terms of strain relaxation routes and bonding configurations. Due to the presence of an energetic hierarchy of residence sites of adatoms, site we also plot a 2D phase diagram consisting of several surface phases. Our EELS results show that the electronic properties of the surface phases are unique to their respective structural arrangement.     

Coverage effects on phenol adsorption on Si(1 0 0)2 × 1 as: A first principle calculation

We investigated the adsorption of a 6-dimers Si(1 0 0)2 × 1 surface as a function of coverage and adsorption type (molecular/dissociative) by first principle calculations. In particular, we performed calculations on models with 2, 3, 4 and 6 phenol molecules, corresponding to coverage Θ = 0.34, 0.5, 0.67 and 1. We found that total adsorption energy, when at least one phenol is in a molecular state is lower than the sum of the corresponding singly adsorbed molecules. The dissociative adsorption of multiple molecules, both in parallel and switched configuration is most favoured for a coverage Θ = 0.34 (2.6 eV per adsorbed molecule). This values decreases to 2.0 eV and remains constant till the coverage 1 is reached.The energy barrier for the molecular-to-dissociated transition of a phenol molecule, in presence of another dissociatively adsorbed molecule is ∼0.008 eV and it is similar to the value in case of single adsorption. Possible hydrogen displacements were also considered.     

Adsorption of Al on the Si(0 0 1) surface

The first-principles calculations have been presented to study the adsorption of aluminum (Al) on the Si(0 0 1)(2×1) surface. We have investigated the optimized geometries and electronic structures of the adatom-substrate system. The adsorption energy of the system has been calculated. The most stable adsorption sites were consequently determined to be HH site and T3+T4. It is shown that the Si-Si dimer is asymmetric on the reconstructed bare surface and become symmetric upon Al adsorption. In addition, the bond length of Si-Si was found to be considerably elongated in the adsorption system. It is found that the work function change obtained in our work is different from other previous results on the adsorption of alkali metals on the Si(0 0 1) surface. In order to investigate the relative stability of phases at different coverages, the surface formation energy of the adsorption system was calculated. To shed light on the nature of the Al-Si bond and the character of silicon surface, the density of states (DOS) and difference charge density of the system were evaluated.     

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.     

First principles calculations of the adsorption and diffusion of Y on the Si(0 0 1)-c(4 × 2) surface

We have studied the adsorption and diffusion of yttrium on the Si(0 0 1)-c(4 × 2) surface in the early stages of growth. Our first principles total energy calculations are based on the density functional theory as implemented in the SIESTA code. The exchange and correlation energies are treated within the generalized gradient approximation according to the Perdew, Burke and Ernzerhof parametrization. Our results demonstrate that the most favorable adsorption site is in the trench between two silicon dimer rows, identified as valley-bridge (V). Our studies show that the diffusion of an Y adatom on Si(0 0 1)-c(4 × 2) surface presents an anisotropic behavior. We found two values for the barriers along the valley (0.54 and 1.07 eV) and one of 1.24 eV in the perpendicular direction, showing that diffusion along the valley is more probable. The analysis of the Mulliken overlap populations shows that the bonding between an Y adatom and the surface is partially covalent. Two Y atoms on the surface do not form dimers instead they are adsorbed as adatoms.     

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.     

First principles calculations of the adsorption of acetylene on the Si(0 0 1) surface at low and full coverage

The adsorption of acetylene on the Si(0 0 1)-c(2 × 4) surface at low and full coverage is studied by first principles density functional calculations using the generalized gradient approximation. For a single acetylene molecule, the most stable configuration corresponds to the di-σ site, on-top of a silicon dimer. This configuration is 0.36 eV more stable than the end-bridge site between two adjacent Si dimers. However, if there are two acetylene molecules, the paired end bridge configuration becomes the most stable. We have also studied the kinetics of the adsorption of a single acetylene molecule. Our calculations show that the reaction is barrier-free for adsorption in the di-σ configuration, while there is an energy barrier of 0.19 eV for adsorption in the end-bridge site. At monolayer coverage, the most stable configuration corresponds to acetylene molecules in the pair-end bridge configuration, in agreement with previous calculations. We have found a noticeable coverage dependence only for the end-bridge site, but not for the di-σ. Our results show that to have an accurate picture of the adsorption of acetylene on the Si(0 0 1) surface, a large unit cell is needed.     

Density functional study of ethylamine and allylamine on Si(1 0 0)-2 × 1 and Ge(1 0 0)-2 × 1 surfaces

We have investigated the interactions of ethylamine and allylamine with models of the Si(1 0 0)-2 × 1 and Ge(1 0 0)-2 × 1 semiconductor surfaces. Ab initio molecular orbital calculations, along with density functional theory (DFT), are used to examine the interaction of these amines with cluster models of the semiconductor surfaces. The transition states and final adsorption products for adsorption of the molecules are predicted. The DFT calculations show the amines form N-dative bond states with Si(1 0 0)-2 × 1 or Ge(1 0 0)-2 × 1 as the initial adsorption product. The initial dative-bond products can be further activated, resulting in N-H bond cleavage on both surfaces. The overall reaction of a given amine on Si(1 0 0) via N-H dissociation is more exothermic than on the Ge(1 0 0) surface.     

Interaction of copper with sulfur on the sulfur-terminated Si(1 1 1)-(7 × 7) surface

The adsorption of S₂ on the Si(1 1 1)-(7 × 7) surface and the interaction of copper and sulfur on this sulfur-terminated Si(1 1 1) surface have been studied using synchrotron irradiation photoemission spectroscopy and scanning tunneling microscopy. The adsorption of S₂ at room temperature results in the passivation of silicon dangling bonds of Si(1 1 1)-(7 × 7) surface. Excessive sulfur forms S_n species on the surface. Copper atoms deposited at room temperature directly interact with S-adatoms through the formations of Cu-S bonds. Upon annealing the sample at 300 °C, CuS_x nanocrystals were produced on the sulfur-terminated Si(1 1 1) surface.     

A theoretical study of Ge adsorption on Si(0 0 1) covered with Bi nanolines

The energetic stability and the equilibrium geometry of Ge adsorption on the Si(0 0 1) surface covered with Bi nanolines were examined by ab initio total energy calculations. We find that there is an energetic preference of Ge atoms lying on the Si(0 0 1) terraces, forming Si^down-Ge^up mixed dimers. Further investigations reveal a repulsive interaction between the mixed dimers and the Bi nanolines, suggesting that the formation of Si^down-Ge^up dimers can be tailored by the presence of the Bi nanolines.     

Formation of vinylidene and vinylene by selective reactivity of Si(1 0 0)2 × 1 towards iso, cis and trans isomers of dichloroethylene

The room temperature (RT) chemisorption of three (iso, cis and trans) isomers of dichloroethylene (DCE) on Si(1 0 0)2 × 1 have been investigated by X-ray photoelectron spectroscopy (XPS) and temperature programmed desorption (TPD). Unlike ethylene, the lack of molecular desorption features in the TPD data effectively rules out the cycloaddition adsorption mechanism for all three isomers. XPS spectra show that cis- and trans-DCE adsorb dissociatively on the 2 × 1 surface in equal proportion as mono-σ bonded 2-chlorovinyl and di-σ bonded vinylene adspecies, which could be produced by dechlorination mechanisms involving the proposed tri-atom π-complex and diradical intermediates, respectively. Acetylene (m/z 26) evolution from 2-chlorovinyl adspecies at 590 K and vinylene at 750 K are also observed for both cis- and trans-DCE, further confirming the common adsorption mechanisms for these geometrical isomers and the relative stabilities of the adspecies. In contrast, only vinylidene adspecies is found for iso-DCE, which indicates that the high ionicity of the CCl₂ group favours the diradical dechlorination mechanism. The single m/z 26 desorption peak for iso-DCE adspecies observed at a higher temperature (780 K) than cis and trans isomers is consistent with the higher adsorption energy of vinylidene than vinylene on Si(1 0 0) obtained in our ab initio calculations. The different relative locations of the Cl atoms in these isomers therefore play a crucial role in controlling the adsorption and thermal evolution on Si(1 0 0)2 × 1. The selective reactivity of the 2 × 1 surface towards these isomers can be used to generate vinylene or vinylidene templates from their corresponding adspecies.     

First-principles study of adsorption and diffusion on Ge/Si(0 0 1)-(2 × 8) and Ge/Si(1 0 5)-(1 × 2) surfaces

Using first-principles total-energy calculations, we have investigated the adsorption and diffusion of Si and Ge adatoms on Ge/Si(0 0 1)-(2 × 8) and Ge/Si(1 0 5)-(1 × 2) surfaces. The dimer vacancy lines on Ge/Si(0 0 1)-(2 × 8) and the alternate S_A and rebonded S_B steps on Ge/Si(1 0 5)-(1 × 2) are found to strongly influence the adatom kinetics. On Ge/Si(0 0 1)-(2 × 8) surface, the fast diffusion path is found to be along the dimer vacancy line (DVL), reversing the diffusion anisotropy on Si(0 0 1). Also, there exists a repulsion between the adatom and the DVL, which is expected to increase the adatom density and hence island nucleation rate in between the DVLs. On Ge/Si(1 0 5)-(1 × 2) surface, the overall diffusion barrier of Si(Ge) along direction is relative fast with a barrier of ∼0.83(0.61) eV, despite of the large surface undulation. This indicates that the adatoms can rapidly diffuse up and down the (1 0 5)-faceted Ge hut island. The diffusion is also almost isotropic along [0 1 0] and directions.     

Ge(5 × 5) as a wetting layer on Si(1 1 1)

Ab initio total energy methods are used to investigate the effects on a Ge(1 1 1)-5 × 5 surface of the lateral compressive stress that would be due to a Si substrate, and the effects of intermixing at the interface with the substrate. The effects of stress due to the lattice mismatch between Si and Ge are studied on a Ge slab by changing the lattice constant in the surface plane from that of experimental bulk diamond Ge to that of Si. When this is done the height difference of the Ge adatoms in the faulted half-cell from those in the unfaulted half is accentuated. Effects on the Ge surface due to the presence of the Si-Ge interface were studied using a thin Ge layer on a Si substrate. The presence of the substrate leads to corrugations with significant height differences appearing among the faulted adatoms. The energetics of intermixing were investigated for Si-Ge single atom interchanges. Additional corrugations resulted from the shortened bondlengths due to the Si impurity in the wetting layer.     

DFT study of NH3 dissociation on Si(1 1 1)-7 × 7. The role of intermolecular interactions

The adsorption of NH₃ molecule on the Si(1 1 1)-7 × 7 surface modelled with a cluster has been studied using density functional theory (DFT). The results indicate the existence of a precursor state for the non-dissociative chemisorption. The active site for the molecular chemisorption is the adatom; while the NH₃ molecule adsorbs on the Si restatom via this preadsorbed state, the adsorption on the Si adatom is produced practically without an energy barrier. The ammonia adsorption on the adatom induces an electron transfer from the dangling bond of this atom to the dangling bond of the adjacent Si restatom, hindering this site for the adsorption of a second NH₃ incoming molecule. However, this second molecule links strongly by means of two H-bonds. The dissociative chemisorption process was studied considering one and two ammonia molecules. For the dissociation of a lonely NH₃ molecule an energy barrier of ∼0.3 eV was calculated, yielding NH₂ on the adatom and H on the restatom. When two molecules are adsorbed, the NH₃-NH₃ interaction yields the weakening of a N-H bond of the ammonia molecule adsorbed closer the Si surface. As a consequence, the dissociation barrier practically disappears. Thus, the presence of a second NH₃ molecule at the adatom-restatom pair of the Si(1 1 1)-7 × 7 surface makes the dissociative reaction self-assisted, the total adsorption process elapsing with a negligible activation barrier (less than 0.01 eV).     

Adsorption of iso-butane on ZnO(0 0 0 1)-Zn

Presented are thermal desorption spectroscopy (TDS) and adsorption probability measurements of iso-butane on the Zn-terminated surface of ZnO. The initial adsorption probability, S₀, decreases linearly from 0.57 to 0.22 (±0.02) with impact energy, E_i = 0.74-1.92 eV, and is independent of adsorption temperature, T_s = 91-114 K (±5 K), indicating non-activated molecular adsorption. The coverage, Θ, dependent adsorption probabilities, S(Θ), show a cross-over from adsorbate-assisted adsorption (S increases with Θ) to Kisliuk-like dynamics at about the desorption temperature of iso-butane bi-layers (∼110 K). Thus, the adsorption dynamics are precursor-mediated. The enhanced (gas-surface) mass-match, caused by forming a second layer of the alkane, leads to adsorbate-assisted adsorption. A direct fitting procedure of the TDS data yields a pre-exponential factor of 2.5 × 10¹³/s and a coverage dependent heat of adsorption of E_d(Θ) = 39 − 6 ∗ Θ + 2.5 ∗ exp(−Θ/0.07) kJ/mol.     

Water adsorption on hydrogenated Si(1 1 1) surfaces

The adsorption of water on the hydrogen terminated Si(1 1 1) surface is studied by means of first-principles calculations as well as contact angle measurements. Possible initial adsorption configurations for single water molecules and the potential energy surface are calculated. Only small adsorption energies of the order of meV are predicted. Calculations for higher coverage show that the water-water interactions are stronger than the water-surface bonding. The contact angle formed between a water droplet on the surface approximated from the total-energy calculations amounts to 88°, while our measured value is 91°.     

Density-functional study of the CO adsorption on ferromagnetic Co(0 0 0 1) and Co(1 1 1) surfaces

The regular CO overlayers at coverage θ = 1/3 adsorbed on the (0 0 0 1) surface of hcp Co and (1 1 1) surface of fcc Co are studied by first-principles density-functional theory with the exchange-correlation component in the PBE form. Adsorption in atop, bridge, and three-fold hcp or fcc position are considered. The adsorption energies, CO stretching frequencies, geometry, work function, and local magnetic moments are studied, and, when possible, compared with experimental or theoretical data. Particularly, we show that the recently proposed correction to adsorption energy of CO prefers correctly the atop adsorption site, whereas the remaining sites are almost degenerate in energy. The CO molecule lowers magnetization on neighbouring Co atoms, and the effect decreases with the adsorption site coordination. We show, however, that this trend is not the result of the different C-Co separation at different adsorption sites. A very small magnetic moment appears on CO that couples antiferromagnetically to Co. Most results are very similar for the Co(0 0 0 1) and Co(1 1 1) surfaces.     

A structural parallel between Ge- and Si-induced 4 × 4 and 3 × 3 reconstructions on SiC(0 0 0 1) drawn from comparative RHEED oscillations

The initial Ge growth stages on a (√3 × √3)R30°-reconstructed SiC(0 0 0 1) surface (√3) have been studied using a complete set of surface techniques such as reflection high energy electron diffraction (RHEED), low energy electron diffraction (LEED), atomic force microscopy (AFM) and photoemission and compared with similar Si surface enrichments in place of Ge. The investigations essentially focus on the wetting growth-regimes that are favoured by the use of the √3 surface as a starting substrate, this surface being the closest to a smooth and ideally truncated Si-terminated face of hexagonal SiC(0 0 0 1). Depending on temperature and Ge or Si coverages, varying surface organizations are obtained. They range from unorganized layer by layer growths to relaxed Ge(1 1 1) or Si(1 1 1) island growths, through intermediate attempts of coherent and strained Ge or Si surface layers, characterized by 4 × 4 and 3 × 3 surface reconstructions, respectively. RHEED intensity oscillation recordings, as a function of Ge or Si deposited amounts, have been particularly helpful to pinpoint the limited (by the high lattice mismatch) existence domains of these interesting coherent phases, either in terms of formation temperature or surface coverages. Prominently comparable data for these two Ge- and Si-related reconstructions allow us to propose an atomic model for the still unexplained Ge-4 × 4 one. It is based on a same local organization in trimer and ad-atom units for the Ge excess as admitted for the Si-excess of the 3 × 3 surface, the higher strain nevertheless favouring arrangements, for the Ge-units, in 4 × 4 arrays instead of 3 × 3 Si ones. Admitting such models, 1.25 and 1.44 monolayers of Ge and Si, should, respectively, be able to lie coherently on SiC, with respective lattice mismatches near 30% and 25%. The experimental RHEED-oscillations values are compatible with such theoretical ones. Moreover, these RHEED coverage determinations (for layer completion, for instance) inform us in turn about the initial Si richness of the starting √3 reconstruction and help us to discriminate between earlier contradictory atomic models proposed in the literature.