Substantially low desorption barriers in recombinative desorption of deuterium from a Si(1 0 0) surface

We have studied desorption kinetics of deuterium molecules from a Si(1 0 0) surface by means of temperature-programmed desorption (TPD) spectra and isothermal desorptions.Three desorption components, denoted as β_1,A,β_1,B, and C, can be distinguished in semi-logarithmic plots of the TPD spectra.Their peak positions and intensities are strongly affected by the surface preparation methods employed, either with or without annealing to control the initial D coverage .Peak C appears at the leading edge of the TPD peak.It accounts for only about 5% of the TPD peak at and it diminishes rapidly with decreasing , vanishing at .In contrast, together the β_1,A and β_1,B peaks account for the whole TPD peak at any less than 1.0 ML. The maximum of the β_1,A peak is nearly constant at around the maximum temperature of the TPD peak.On the other hand, the β_1,B peak appears on the high-temperature side of the TPD peak and it systematically shifts to higher temperatures with decreasing .These results imply that first- and second-order kinetics are operating for the β_1,A and β_1,B desorptions, respectively.Isothermal desorption experiments confirm the above predicted kinetics for a limited region, namely .From the results for the rate curve analysis, the desorption barriers are evaluated to be 1.6 ± 0.1 eV and 1.8 ± 0.1 eV for the β_1,A and β_1,B desorptions, respectively.These values are substantially lower than the widely accepted value of ∼2.5 eV. To reproduce the measured TPD spectra we take the Arrhenius-type rate equation containing the first- and second-order rate terms for the β_1,A and β_1,B desorptions.The TPD spectra measured for can be reasonably fit with the proposed rate equation when the values given above for E_d,A and E_d,B are used. For , however, the TPD curves are not fit with the same values; rather, the best-fit curves require values for E_d,A and E_d,B larger than those given above. Combining the present kinetics results with those obtained by STM along with the studies, the β_1,A and β_1,B peaks may be attributed to desorption along the 2H path, while peak C may be attributed to desorption along the 4H path. The atomistic desorption mechanism as well as the energy relationship between the desorption barrier and isosteric heat of adsorption are discussed.     

A scanning tunneling microscopy study of PH3 adsorption on Si(1 1 1)-7 × 7 surfaces, P-segregation and thermal desorption

PH₃ adsorption on Si(1 1 1)-7 × 7 was studied after various exposures between 0.3 and 60 L at room temperature by means of scanning tunneling microscopy (STM). PH₃-, PH₂-, H-reacted, and unreacted adatoms can be identified by analyzing empty-state STM images at different sample biases. PH_x-reacted rest-atoms can be observed in empty-state STM images if neighboring adatoms are hydrogen terminated. Most of the PH₃ adsorbs dissociatively on the surface, generating H- and PH₂-adsorbed rest-atom and adatom sites. Dangling-bonds at rest-atom sites are more reactive than adatom sites and the faulted half of the 7 × 7 unit cell is more reactive than the unfaulted half. Center adatoms are overwhelmingly preferred over corner adatoms for PH₂ adsorption. The saturation P coverage is ∼0.18 ML. Annealing of PH₃-reacted 7 × 7 surfaces at 900 K generates disordered, partially P-covered surfaces, but dosing PH₃ at 900 K forms P/Si(1 1 1)- surfaces. Si deposition at 510 K leaves disordered clusters on the surface, which cannot be reordered by annealing up to 800 K. However, annealing above 900 K recreates P/Si(1 1 1)- surfaces. Surface morphologies formed by sequential rapid thermal annealing are also presented.     

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

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

Oxygen atom-induced D2 and D2O desorption on D/Si(1 1 1) surfaces

We studied reaction of oxygen atoms with D-terminated Si(1 1 1) surfaces from a desorption point of view. As the D (1 ML)/Si(1 1 1) surface was exposed to O atoms D₂ and D₂O molecules were found to desorb from the surface. The desorption kinetics of D₂ and D₂O molecules exhibited a feature characterized with a quick rate jump at the very beginning of O exposure, which was followed by a gradual increase with a delayed maximum and then by an exponential decrease. The O-induced D₂ desorption spectra as a function of T_s appeared to be very similar to the H-induced D₂ desorption spectrum from the D/Si(1 1 1) surfaces. Possible mechanisms for the O-induced desorption reactions were discussed.     

Reactions of atomic oxygen with the D-covered Si(1 0 0) surfaces

We have studied D abstraction by O on the D/Si(1 0 0) surfaces using a continuous as well as pulsed O-beams. Both D₂ and D₂O molecules are detected during O-exposure. The D₂ desorption is found to take place more efficiently on the monodeuteride/dideuteride surface than on the monodeuteride surface. The pulsed beam experiments exhibit occurrence of both a slow and a fast D₂ desorption. The D₂ desorption is found to obey the second-order rate law in on the monodeuteride surfaces and 3.5th-order rate law on the monodeuteride/dideuteride surfaces. The D₂O desorption is found to be governed also by the second-order rate law, however regardless of D coverage even on the monodeuteride/dideuteride surfaces. Possible mechanisms for the O-induced desorption from the D/Si(1 0 0) surfaces are discussed.     

D abstraction by H at a D-saturated Ru(0 0 1) surface

D abstraction (ABS) by H at Ru(0 0 1) surfaces initially saturated with D adatoms has been investigated using in situ mass spectrometry. HD and D₂ desorption rates are measured at various surface temperatures T as a function of H exposure time. Yield of D₂ desorption increases with T, while that of HD is little affected. Analyzing the measured rate curves, HD and D₂ desorption orders are evaluated to be 1.7 ± 0.1 and 2.5 ± 0.1, respectively, with respect to D coverage θ_D. To pursue the origin of the derived non-integral reaction orders the rate curves are further analyzed with the rate equations constructed to involve several ABS channels. Consequently, we find that the HD desorption is mainly governed by a second-order rate law in θ_D rather than the conventional hot atom-mediated ABS reaction even when it is corrected to include an isotope effect on ABS. We argue that such second-order ABS kinetics becomes important when the H atoms in excited state of chemisorption have energetically relaxed to some extent, and thereby tend to reside at, e.g. hexagonal closed packed hollow sites, interacting with nearby adatoms. On the other hand, the D₂ rate curves can be fit with third-order kinetics, consistent with the Langmuir-Hinshelwood mechanism in a super-saturation state. The isotope effect plays an essential role in the ABS reaction of D abstraction by H which competes with H abstraction by H as D adatoms are replaced by H atoms.     

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.     

O2 dissociation on Pd(2 1 1) and Cu(2 1 1) surfaces

We have performed ab initio Density Functional Theory (DFT) based calculations to observe the reactivity of the Pd(2 1 1) and Cu(2 1 1) surfaces towards O₂. In order to properly address the adsorption dynamics, the static potential energy surface calculations have been complemented with first principles molecular dynamics calculations, which reveal interesting steering effects that complicate the dissociation dynamics. We have found that on both surfaces the step microfacets are very reactive and the dissociation of the O₂ molecule at room temperature occurs mostly on those sites.     

Initial stages of H2O adsorption and hydroxylation of Fe-terminated α-Fe2O3(0 0 0 1) surface

The adsorption modes of H₂O on a Fe-terminated hematite(0 0 0 1) surface have been investigated by first principles Density Functional theory within a periodic slab model and the generalized gradient approximation. Molecular adsorption and dissociative adsorption in monolayer coverage, one H₂O per surface Fe, were both considered. Five plausible orientations were studied to determine the most favorable adsorption position. Molecular adsorption is shown to have a small effect on the underlying surface structure, while hydroxylation has a strong effect on the surface geometry. Electronic densities of state calculations reveal details of these different interactions. The heterolytic dissociation, which produces two types of surface hydroxyls, is the preferable adsorption mode, being slightly favored energetically over the molecular adsorption. Homolytic dissociative adsorption, forming a single hydroxyl on surface Fe, is energetically unfavored, even though strong binding interaction (∼3 eV) is found between the OH radical and surface. Dissociative adsorption on an oxidized ferryl site was also studied to investigate suggested local reactivity enhancement.     

D2 layers on MgO(0 0 1): Simulation study

In response to recent helium atom scattering (HAS) and neutron scattering results, Monte Carlo simulations and perturbation theory calculations have been performed for D₂ on MgO(0 0 1). Monte Carlo simulations predict that D₂ molecules form a series of interesting structures, p(2 × 2) → p(4 × 2) → p(6 × 2), with coverages Θ = 0.5, 0.75, 0.83 respectively, and followed by a formation of a top layer of p(6 × 2) unit cell symmetry. The three types of mono-layers are stable up to 13 K, whereas the top layer still exists up to 10 K. This is in partial agreement with the neutron scattering and HAS results that report c(2 × 2), c(4 × 2) and c(6 × 2); they agree in terms of coverage and stability, but disagree in terms of symmetry. A quantum mechanical examination of the D₂ molecules’ rotational motion shows the molecular axes are azimuthally delocalized and hence the simulated structures are c-type rather than p-type. These calculations also indicate that ortho-D₂ and helicoptering para-D₂ prefer cationic sites, while cartwheeling para-D₂ prefers anionic sites.     

Adsorption of CO, CO2 and NO molecules on a BaTiO3 (0 0 1) surface

Since the development of Scanning Tunnelling Microscopy (STM) technique, considerable attention has been devoted to various molecules adsorbed on various surfaces. Also, a new concept emerged with molecules on surfaces considered as nano machines by themselves. In this context, a thorough knowledge of surfaces and adsorbed molecules at an atomic scale are thus particularly invaluable. The present work describes the first Density Functional Theory (DFT) study of adsorption of CO, CO₂ and NO molecules on a BaTiO₃ surface following a first preliminary calculation of O and O₂ adsorption on the same surface. In the previously considered work, we found that a (0 0 1) surface with BaO termination is more stable than the one with TiO₂-termination. Consequently, we extended our study to CO, CO₂ and NO molecules adsorbed on a (0 0 1) surface with BaO termination. The present calculation was performed on a (1 × 1) cell with one monolayer of adsorbed molecules. Especially, a series of cases implying CO molecules adsorbed in various geometrical configurations has been examined. The corresponding adsorption energy varies in the range of −0.17 to −0.10 eV. The adsorption energy of a CO₂ molecule directly located above an O surface atom (called O_s) is of the order of −0.18 eV. The O-C distance length is then 1.24 Å and the O-C-O and O-C-O_s angles are 134.0° and 113.0°, respectively. For NO adsorption, the most important induced structural changes are the followings: (i) the N-O bond is broken when a NO molecule is absorbed on a Ba-O_s bridge site. In that case, N and O atoms are located above an O and a Ba surface atom, respectively, whereas the O-Ba-O_s and N-O_s-Ba angles are 106.5° and 63.0°, respectively. The N-O distance is as large as 2.58 Å and the adsorption energy is as much as −2.28 eV. (ii) In the second stable position, the NO molecule has its N atom adsorbed above an O_s atom, the N-O axis being tilted toward the Ba atom. The N-O_s-Ba angle is then 41.1° while the adsorption energy is only −0.10 eV. At last, the local densities of states around C, O as well as N atoms of the considered adsorbed molecules have also been discussed.     

Intramolecular structures of C60 and C84 molecules on Si(1 1 1)-7 × 7 surfaces by scanning tunneling microscopy

The intramolecular features of carbon 60 and carbon 84 molecules on Si(1 1 1)-7 × 7 surfaces were studied under a UHV-scanning tunneling microscope. Carbon molecules preferentially appear in faulted halves, rather than in unfaulted halves and corner holes; they are embedded in silicon substrates. The orientation and details of the structure of carbon molecules are determined by applying various sample biases to the silicon substrate. As compared with other fullerenes, a bright pentagonal ring with nebulous clusters which represents the cage structure is clearly observed on top of carbon 60 molecules. The bright stripes associated with partitioned curves which depict eight features of asymmetrical C₈₄ molecules are also investigated on Si(1 1 1)-7 × 7 surfaces. The orientations and possible configurations of C₆₀ and C₈₄ are considered in this work. The energy differences for various features of C₆₀ and C₈₄ molecules are estimated and discussed. The corresponding models with respect to each intramolecular feature are proposed and compared with recent theoretical calculation.     

CO2 dissociation on Ni(2 1 1)

The direct and H-mediated dissociation of CO₂ on Ni(2 1 1) were investigated at the level of density functional theory. Although formate (HCOO) formation via CO₂ hydrogenation was widely reported for CO₂ adsorption on metal surfaces, it is found that on Ni(2 1 1) HCOO dissociation into CHO and O is much difficult, while direct dissociation of adsorbed CO₂ into CO and O is more favorable. It is also found that the degree of electron transfer from surface to adsorbed CO₂ correlates with the elongation of C-O bond lengths and the reduction of the CO₂ dissociation barrier.     

Interaction between chemisorbed N2 and Li promoter atoms: A comparison between the stepped Ru(1 0 9) and the atomically smooth Ru(0 0 1) surfaces

We present a direct side-by-side comparison of the interaction of Li atoms and N₂ molecules on the atomically stepped Ru(1 0 9) single crystal surface and on the atomically smooth Ru(0 0 1) single crystal surface using infrared reflection absorption spectroscopy (IRAS) and temperature programmed desorption (TPD). At low adsorbate coverages there is spectroscopic evidence for the formation of a Li_x(N₂)_y complex on the Ru(1 0 9) surface, whereas no such complex is observed on the Ru(0 0 1) surface. This complex is due to local interactions between an adsorbed Li atom and N₂ adsorbed on the atomic steps of Ru(1 0 9). The short range interaction near the atomic steps is characterized by the development of several highly red-shifted ν(N₂) modes in the region of ∼2130 cm⁻¹ in the IR spectra. Adsorbed N₂ molecules on both Ru(1 0 9) and Ru(0 0 1) also are influenced by the long range electrostatic field produced by Li adsorbate atoms, causing a red shift in the uncomplexed N₂ species, which monotonically increases as the Li coverage in increased. On the Ru(0 0 1) surface, small coverages of N₂ influenced by the long range effect of Li are initially chemisorbed parallel to the surface resulting in the absence of infrared activity. In addition we have also found that Li does not cause N-N bond scission on Ru(0 0 1) below 250 K.     

Si 2p and O 1s photoemission from oxidized Si(0 0 1) surfaces depending on translational kinetic energy of incident O2 molecules

The influence of translational kinetic energy of incident O₂ molecules for the passive oxidation of the partially oxidized Si(0 0 1) surface has been studied by photoemission spectroscopy. The incident energy of O₂ molecules was controlled up to 3 eV by a supersonic molecular beam technique. Two incident energy thresholds (1.0 and 2.6 eV) were found out in accordance with the first-principle calculations. Si 2p and O 1s photoemission spectra measured at representative incident energies showed the incident energy induced oxidation at the backbonds of the dimer and the second layer (subsurface) Si atoms. Moreover, the difference of oxygen chemical bonds was found out to be as the low and the high binding energy components in the O 1s photoemission spectra. They were assigned to bridge sites oxygen and dangling bond sites oxygen, respectively.     

Theoretical study of oxygen adsorption at the Fe(1 1 0) and (1 0 0) surfaces

Electronic, magnetic and structural properties of atomic oxygen adsorbed in on-surface and subsurface sites at the two most densely packed iron surfaces are investigated using density functional theory combined with a thermodynamics formalism. Oxygen coverages varying from a quarter to two monolayers (MLs) are considered. At a 1/4 ML coverage, the most stable on-surface adsorption sites are the twofold long bridge sites on the (1 1 0), and the fourfold-hollow sites on the (1 0 0) surface. The presence of on-surface oxygen atoms enhances the magnetic moments of the atoms of the two topmost Fe layers. Detailed results on the surface magnetic properties, due to O incorporation, are presented as well. Subsurface adsorption is found unfavored. The most stable subsurface O, in tetrahedral positions at the (1 0 0) and octahedral ones at the (1 1 0) surface, are characterized by substantially lower binding than that in the on-surface sites. Subsurface oxygen increases the interplanar distance between the uppermost Fe layers. The preadsorbed oxygen overlayer enhances binding of subsurface O atoms, particularly for tetrahedral sites beneath the (1 1 0) surface.     

Hopping motion of chlorine atoms on Si(1 0 0)-(2 × 1) surfaces induced by carrier injection from scanning tunneling microscope tips

Injection of tunneling electrons and holes from the probe tips of a scanning tunneling microscope was found to enhance the hopping motion of Cl atoms between neighboring dangling-bond sites of Si dimers on Si(1 0 0)-(2 × 1) surfaces, featured by the rate of hopping linearly dependent on the injection current. The hopping rate formed peaks at sample biases of V_S∼+1.25 and −0.85 V, which agree with the peaks in the local density of states spectrum measured by scanning tunneling spectroscopy. The Cl hopping was enhanced at Cl-adsorbed sites even remote from the injection point. The Cl hopping by hole injection was more efficiently enhanced by sweeping the tip along the Si dimer row than by tip-sweeping along the perpendicular direction. Such anisotropy, on the other hand, was insignificant in the electron injection case. All of these findings can be interpreted by the model that the holes injected primarily into a surface band originated from the dangling bonds of Si dimers propagate quite anisotropically along the surface, and become localized at Cl sites somehow to destabilize the Si-Cl bonds causing hopping of the Cl atoms. The electrons injected into a bulk band propagate in an isotropic manner and then get resonantly trapped at Si-Cl antibonding orbitals, resulting in bond destabilization and hopping of the Cl atoms.