Diffusion of oxygen atom in the topmost layer of the Si(1 0 0) surface: Structures and oxidation kinetics

The incorporations and migrations of the atomic oxygen in the topmost layer Si(1 0 0)-p(2 × 2) silicon surface, are investigated theoretically using density functional theory. We show that the diffusion is dependent on the starting and the final surrounding environment and does not simply consist in hops from one silicon-silicon bond to another. The activation energies range from 0.11 eV to 2.59 eV.     

Molecular nano-arches on silicon

The formation of molecular nano-arches on the Si(1 1 1)-7 × 7 surface was modeled using density functional theory (DFT). It has been suggested, based on the calculations, that the arches are formed by molecular dimers of chlorobenzene at near-monolayer coverages. Molecules of the dimer are covalently bound to two silicon adatoms and to each other thereby forming a molecular arch on the surface. The structure of the molecular dimer was calculated at the B3LYP/6-31G(d) level of theory. The dimers were found to be stable at room temperature, and to form a near-monolayer coverage, which has been observed in the experiment [X.H. Chen, Q. Kong, J.C. Polanyi, D. Rogers, S. So, Surf. Sci. 340 (1995) 224; Y. Cao, J.F. Deng, G.Q. Xu, J. Chem. Phys. 112 (2000) 4759].     

Theoretical investigation of the initial steps of the adsorption of N atoms on Si(1 0 0)-2×1

Structural, energetics, and mechanistics aspects of initial steps of the reaction of a N atom with Si(1 0 0)-2×1 modeled by the Si₉H₁₂+N system are reported. Hybrid density functional B3LYP calculations predict a barrierless first step leading to an adsorbate where N is bound to one of the dimer Si. Two possible activated routes for internal rearrangements were found, with that leading to the incorporation of Si below the first layer predicted to be kinetically dominant (98%) under the experimental conditions. This structure and frequency calculations are consistent with the experimental finding of a planar NSi₃ moeity and with the experimental SiN asymmetric stretching frequency of the NSi₃ groups.     

Active oxidation: Silicon etching and oxide decomposition basic mechanisms using density functional theory

The etching of silicon atom from the Si(1 0 0)-p(2 × 2) surface, i.e. the desorption of SiO molecules from this surface, either clean or pre-oxidized, is investigated at the density functional theory level. The reaction paths for desorption are given as a function of the initial oxidation state of the extracted silicon atom. The associated activation energies and the atomic configurations are discussed. Particularly, it is shown that desorption of SiO molecules takes place during conventional thermal oxide growth (∼2 eV activation) via non-oxidized silicon atoms. Further SiO extraction mechanisms of higher silicon oxidation states required higher temperatures. In particular, doubly oxidized silicon atoms (Si²⁺) are able to decompose with an activation of ∼4 eV which corresponds to the actual temperature where decomposition of oxides is observed. This confirms the statement that decomposition of oxide layer nucleates at the interface with silicon where Si²⁺ has been detected thanks to XPS experiment.     

A density functional theory study of formaldehyde adsorption and oxidation on CeO2(1 1 1) surface

The effects of different oxygen species and vacancies on the adsorption and oxidation of formaldehyde over CeO₂(1 1 1) surface were systematically investigated by using density functional theory (DFT) method. On the stoichiometric CeO₂(1 1 1) surface, the C-H bond rupture barriers of chemisorbed formaldehyde are much higher than that of formaldehyde desorption. On the reduced CeO₂(1 1 1) surface, the energy barriers of C-H bond ruptures are less than those on the stoichiometric CeO₂(1 1 1) surface. If the C-H bond rupture occurs, CO and H₂ form quickly with low energy barriers. When O₂ adsorbs on the reduced (1 1 1) surface (O₂/O_v species), the C-H bond rupture barriers of formaldehyde are greatly reduced in comparison with those on the stoichiometric CeO₂(1 1 1) surface. If O₂ adsorbs on oxygen vacancy at sub-layer surface, its oxidative roles on formaldehyde are much similar to that of O₂/O_v species.     

First principles study of Si(3 3 5)-Au surface

The structural and electronic properties of gold decorated Si(335) surface are studied by means of density-functional calculations. The resulting structural model indicates that the Au atoms substitute some of the Si atoms in the middle of the terrace in the surface layer. Calculated electronic band structure near the Fermi energy features two metallic bands, one coming from the step edge Si atoms and the other one having its origin in hybridization between the Au and neighboring Si atoms in the middle of the terrace. The obtained electronic bands remain in good agreement with photoemission data.     

Investigations of meta-stable and post-spinel silicon nitrides

Since the discovery of post-spinel Si₃N₄, its fundamental physical properties are highly required. In this paper, theoretical calculations are performed to investigate the structural and elastic properties of the β-, γ-, wII- and post-spinel Si₃N₄ polymorphs. The calculated ground-state properties compare well with available experiments. The phase transformations of the β-, γ-, wII- and post-spinel phases are investigated by the famous plane-wave pseudo-potential density functional theory. From the elastic constants obtained, we find that β-, γ- and wII-Si₃N₄ are stable at 0 GPa and the post-spinel phase is unstable/stable at 0 GPa/160 GPa. When the high-temperature β→γ transformation is bypassed due to kinetic reasons, β-Si₃N₄ is predicted to undergo a first-order phase transition to a new phase (wII-Si₃N₄). It is found that the transition pressures of β→wII and γ→post-spinel transitions are 20.8 GPa and 152.5 GPa, respectively. The phase boundary of the γ→post-spinel transition can be described as P=152.3631−6.39×10⁻³T+2.01062×10⁻⁵T²−1.93962×10⁻⁹T³. Through the quasi-harmonic approximation, the dependences of heat capacity, entropy, thermal expansion coefficient and the Debye temperature on temperature, are also successfully predicted.     

Density functional study of initial HfCl4 adsorption and decomposition reactions on silicon surfaces with SiON interfacial layer

The initial adsorption and decomposition of HfCl₄ on silicon surfaces with different types of SiON interfacial layers are investigated using density functional theory. We find that the reactions of HfCl₄ on both the hydroxylated and nitrided silicon surfaces proceed through similar reaction pathways. By comparison of the reaction energies of HfCl₄ with the hydroxyl and amino surface sites, we find that it is both kinetically and thermodynamically favorable for the reactions of HfCl₄ on hydroxyl site of silicon substrates. Comparing with the adjacent bridging oxygen, we also find that the neighboring hydroxyl can facilitate the adsorption of HfCl₄ on the amido surface site. Also, it is more kinetically and thermodynamically favorable for the reaction of HfCl₄ with bridging NH site than that with NH₂ site.     

A review on techniques to fabricate silicon oxide arrays for biomolecules patterning

The development of methods for pattern proteins and other functional molecules on the surfaces with nanoscale accuracy is indispensable to take advantage of their properties in ultrasensitive and/or high-density devices. Several methods are used to fabricate organized micro/nanohierarchical structures on a surface and give the ability of molecules to self-assemble by using the mutual recognition properties. However, the supramolecular organization is difficult to extend from nano- to mesoscopic length scales or does not allow accurate placement of the desired structures on a specific region of an inhomogeneous surface. This paper reviews the different techniques used to fabricate nano/millimeter range patterns on SiO₂/Si substrates and local chemical grafting to perform successful attachment of biomolecules on predetermined areas.     

Vibrational study of silicon oxidation: H2O on Si(100)

The vibrational spectrum of water dissociatively adsorbed on Si(100) surfaces is obtained with surface infrared absorption spectroscopy. Low frequency spectra (< 1450 cm⁻¹ are acquired using a buried CoSi₂ layer as an internal mirror to perform external reflection spectroscopy. On clean Si(100), water dissociates into H and OH surface species as evidenced by EELS results [1] in the literature which show a Si---H stretching vibration (2082 cm⁻¹), and SiO---H vibrations (O---H stretch at 3660 cm⁻¹ and the Si---O---H bend and Si---O stretch of the hydroxyl group centered around 820 cm⁻¹). In this paper, infrared (IR) measurements are presented which confirm and resolve the issue of a puzzling isotopic shift for the Si---O mode of the surface hydroxyl group, namely, that the Si---O stretch of the O---H surface species formed upon H₂O exposure occurs at 825 cm⁻¹, while the Si---O stretch of the ---OD surface species formed upon D₂O exposure shifts to 840 cm⁻¹, contrary to what is expected for simple reduced mass arguments. The higher resolution of IR measurements versus typical EELS measurements makes it possible to identify a new mode at 898 cm⁻¹, which is an important piece of evidence in understanding the anomalous frequency shift. By comparing the results of measurements for adsorption of H¹⁶₂O, H¹⁸₂O and D₂O with the results from recently performed first-principles calculations, it can be shown that a strong vibrational interaction between the Si---O stretching and Si---O---H bending functional group vibrations of the hydroxyl group accounts for the observed isotopic shifts.     

A theoretical study of structural and electronic properties of a missing dimer defect on Si- and C-terminated SiC(0 0 1)

We present a theoretical study of the geometrical and electronic properties of the C- and Si-terminated -SiC(0 0 1) surfaces in the vicinity of the missing dimer defect. The experimental results suggest that the atomic structures of these two surfaces may be considerably modified by external stress. In our present study we have considered the possible influence of this factor on the surface geometry of both systems. We have shown that the structural differences between the C- and Si-terminated surfaces lead to their different behaviour in the presence of a missing dimer and applied stress. In the case of the C-terminated c(2×2) surface, the missing dimer defect causes the buckling of the adjacent carbon dimers lying in the line of the defect (dimer atoms adjacent to the defect have vertical positions lower by 0.18 Å). This effect becomes more pronounced in the presence of compressive stress — the stress of 8% leads to the buckling of these two dimers of around 0.5 Å. The vertical positions of silicon atoms located directly below the defect were increased by 0.2 Å. We have also found that the missing dimer influences the structure of the carbon dimers on the neighbouring lines of dimers. Contrary to the C-terminated surface, the missing dimer defect on the Si-terminated SiC(0 0 1)-p(2×1) surface remains neutral for silicon dimers located in the line of defect, i.e. the dimers do not change their geometrical properties in unstrained structure nor in the presence of a tensile stress. On the other hand, this defect modifies considerably the geometry of the dimers from the two neighbouring lines of dimers by reducing their bond lengths and vertical positions. Changes in the geometrical properties of the second neighbour dimers (with respect to the defect) in these two lines are also noticeable. Moreover, we have found that the presence of a missing dimer modifies significantly the positions of the adjacent subsurface carbon atoms.     

Stress development and impurity segregation during oxidation of the Si(1 0 0) surface

We have studied the segregation of P and B impurities during oxidation of the Si(1 0 0) surface by means of combined static and dynamical first-principles simulations based on density functional theory. In the bare surface, dopants segregate to chemically stable surface sites or to locally compressed subsurface sites. Surface oxidation is accompanied by development of tensile surface stress up to 2.9 Nm⁻¹ at a coverage of 1.5 monolayers of oxygen and by formation of oxidised Si species with charges increasing approximately linearly with the number of neighbouring oxygen atoms. Substitutional P and B defects are energetically unstable within the native oxide layer, and are preferentially located at or beneath the Si/SiO_x interface. Consistently, first-principles molecular dynamics simulations of native oxide formation on doped surfaces reveal that dopants avoid the formation of P-O and B-O bonds, suggesting a surface oxidation mechanism whereby impurities remain trapped at the Si/SiO_x interface. This seems to preclude a direct influence of impurities on the surface electrostatics and, hence, on the interactions with an external environment.     

Deep reduction behavior of iron oxide and its effect on direct CO oxidation

Reduction of metal oxide oxygen carrier has been attractive for direct CO oxidation and CO₂ separation. To investigate the reduction behaviors of iron oxide prepared by supporting Fe₂O₃ on γ-Al₂O₃ and its effect on CO oxidation, fluidized-bed combustion experiments, thermogravimetric analyzer (TGA) experiments, and density functional theory (DFT) calculations were carried out. Gas yield (γCO₂) increases significantly with the increase of temperature from 693 K to 1203 K, while carbon deposition decreases with the increase of temperature from 743 K to 1203 K, where temperature is a very important factor for CO oxidation by iron oxide. Further, it were quantitatively detected that the interaction between CO and Fe₂O₃, breakage of O-Fe bonds and formation of new C-O bonds, and effect of reduction degree were quantitatively detected. Based on adsorptions under different temperatures and reducing processes from Fe³⁺ into Fe²⁺, Fe⁺ and then into Fe, it was found that Fe²⁺ → Fe⁺ was the reaction-controlling step and the high oxidation state of iron is active for CO oxidation, where efficient partial reduction of Fe₂O₃ into FeO rather than complete reduction into iron may be more energy-saving for CO oxidation.     

Initial incorporation of sulfur into the Pd(1 1 1) surface: A theoretical study

Density functional theory is used to investigate the initial inclusion of sulfur into the subsurface interstitial sites of Pd(1 1 1) surface. Pure subsurface adsorption is found to be less energetically favorable than on-surface adsorption. The incorporation of sulfur into the metal becomes more favorable than continuous adsorption on the surface after a critical on-surface sulfur coverage. We find subsurface sulfur occupation to be energetically favorable after adsorption of more than half a monolayer on the surface. Occupation of subsurface sites induces a pronounced structural distortion of the Pd(1 1 1) surface. We find significant expansion of interplanar spacing between the uppermost surface metal layers and rearrangement of the S overlayer. The interplay between the energy cost due to structural distortion of Pd(1 1 1) and the energy gain due to bond formation for different structures is discussed.     

Density functional study of the reaction of O2 with a single site on the zigzag edge of graphene

The reaction between molecular oxygen and an isolated zigzag graphene edge has been studied using density functional theory at the B3LYP/6-31G(d) level of theory. The initial reaction forms a peroxide, ΔH = −135 kJ mol⁻¹. If the graphene edge is pre-oxidised, the dangling peroxy atom can (E_a = 91 kJ mol⁻¹) migrate across contiguous ketone groups until finding another vacant site and stabilizing as a ketone. However, if no further vacant sites are available, the peroxy oxygen has a number of other possibilities open to it, including desorption of an oxygen atom (E_a = 140 kJ mol⁻¹), migration via the basal plane to form a lactone (E_a = 147 kJ mol⁻¹), and direct interaction with an adjacent oxide to form the lactone or a carbonate (E_a = 146 kJ mol⁻¹). The combination of thermal energy and the heat released in the initial formation of the peroxy adduct is likely to be sufficient to overcome these secondary barriers at modest temperatures.Transfer of the dangling peroxy O to the basal plane produces an epoxide that is mobile on the basal surface (E_a = 40–80 kJ mol⁻¹) but that is transferred back to the edge upon coming into proximity of either a vacant edge site or ketone. The instability of the edge epoxide structure implies that it cannot play a significant role in carbon gasification through promoting the reactivity of ketones, contrary to earlier suggestions.The desorption of an oxygen atom creates a very active species capable of reacting with basal or edge sites as well as with oxygen complexes. The reaction of ketone + O has been reported to yield a five-membered ring + CO₂, leading to an overall stoichiometry which is consistent with the observed oxyreactivity of carbon surface oxides identified in isotopic labelling studies in which one O atom is gasified and the other forms a new surface oxide.     

First-principles study on the magnetic properties of transition-metal atoms doped (ZnS)12 cluster

A first-principles density functional investigation has been performed to evaluate the structural, electronic, and magnetic properties of (ZnS)₁₂ doped with one or two transition-metal (TM) atoms (Fe, Co, and Ni). Substitutional- and interstitial-doping are considered. The substitutional isomers are found to be most favorable for Fe-doped clusters, while the interstitial isomers are found to be most favorable for Co- and Ni-doped clusters. Magnetic coupling between the TM atoms at the nearest neighbor position is mainly governed by the competition between direct ferromagnetic and antiferromagnetic interactions between two TM atoms via the S atom due to strong p-d hybridization. The coupling is short-ranged. Most importantly, we demonstrate that the Fe and Ni endohedral bi-doped (ZnS)₁₂ clusters favor the ferromagnetic state, which has potential applications in nanoscale quantum devices.     

Ab initio density functional study of O on the Ag(0 0 1) surface

The adsorption of oxygen on the Ag(1 0 0) is investigated by means of density functional techniques. Starting from a characterization of the clean silver surfaces oxygen adsorption in several modifications (molecularly, on-surface, sub-surface, Ag₂O) for varying coverage was studied. Besides structural parameters and adsorption energies also work-function changes, vibrational frequencies and core level energies were calculated for a better characterization of the adsorption structures and an easier comparison to the rich experimental data.     

First principles investigation on the stabilization mechanisms of the polar copper terminated Cu2O(1 1 1) surface

The stabilization of the unstable, polar copper terminated Cu₂O(1 1 1) surface by reconstruction and hydroxylation was studied theoretically with static and molecular dynamics calculations at ab initio density functional theory (DFT) level. Surface reconstruction was investigated using extensive finite temperature molecular dynamics (MD) combined with a simulated annealing technique. Both the global minimum energy structure obtained during annealing the system at higher temperature (300 K) and the final ‘quenched’ structure which was obtained after cooling the system to 0 K show the expected reconstruction of the adsorbate-free surface. The copper atoms in the first layer and oxygen atoms in the second and third layers are markedly displaced, and the atomic planes merge together to form a uniform mixed layer, thereby minimizing the polarity of the surface. Surface hydroxylation by adsorption of OH⁻ or dissociated water was investigated using static optimization at 0 K. The results show that adsorption is exothermic and that the reconstruction characterizing the annealed OH-free surface does not occur in the presence of adsorbed OH. A surface coverage of 50% results in the surface structure that is the closest to the unrelaxed bulk terminated surface.     

Reactivities of ultrathin alumina films exposed to intermediate pressures of H2O: Substrate-mediated mechanism for growth and loss of surface order

The response of ordered ultrathin Al₂O₃ films on NiAl(1 1 0) and Ni₃Al(1 1 0) substrates to sequential exposures at varying pressures of H₂O between 10⁻⁷ Torr and 10⁻³ Torr, ambient temperature, was characterized by LEED, AES and density functional theory (DFT) calculations. In all cases, an increase in average oxide thickness, as determined by AES, was observed, consistent with a field-induced oxide growth mechanism. Ordered oxide films of initial average thicknesses of 7 Å and 12 Å grown on NiAl(1 1 0) achieved a limiting thickness of 17(1) Å, while films of initial thickness of 7 Å and 11 Å grown on Ni₃Al(1 1 0) achieved a limiting thickness of 12(1) Å. The LEED patterns for the thinner (7 Å) films were not observed after exposure to 10⁻⁵ Torr (NiAl(1 1 0)), or 10⁻⁴ Torr (Ni₃Al(1 1 0)). In contrast, LEED patterns for the films of greater initial thickness persisted after exposures to 10⁻³ Torr UHV. DFT calculations indicate an Al vacancy formation energy that is significantly greater (by ∼0.5 eV) on the surface that has the thicker oxide film, directly opposite to what may be naively expected. A simple coordination argument supports these numerical results. Therefore, the greater limiting oxide thickness observed on NiAl(1 1 0) demonstrates that the rate determining step in the oxide growth process is not Al removal from the metal substrate and transport across the oxide/metal interface. Instead, the results indicate that the determining factor in the oxide growth mechanism is the kinetic barrier to Al diffusion from the substrate bulk to the oxide/metal interface. The persistence of the LEED patterns observed for the films of greater initial oxide thickness indicates that the surface disorder generally observed for alumina films grown on aluminide substrates and exposed to intermediate pressures of H₂O is due to the growth of a disordered alumina layer over an ordered substrate, rather than to direct H₂O interaction with terrace sites.     

Pd diffusion on MgO(1 0 0): The role of defects and small cluster mobility

Density functional theory is used to explore the energy landscape of Pd atoms adsorbed on the terrace of MgO(1 0 0) and at oxygen vacancy sites. Saddle point finding methods reveal that small Pd clusters diffuse on the terrace in interesting ways. The monomer and dimer diffuse via single atom hops between oxygen sites with barriers of 0.34 eV and 0.43 eV respectively. The trimer and tetramer, however, form 3D clusters by overcoming a 2D-3D transition barrier of less than 60 meV. The trimer diffuses along the surface either by a walking or flipping motion, with comparable barriers of ca. 0.5 eV. The tetramer rolls along the terrace with a lower barrier of 0.42 eV. Soft rotational modes at the saddle point lead to an anomalously high prefactor of 1.3 × 10¹⁴ s⁻¹ for tetramer diffusion. This prefactor is two order of magnitude higher than for monomer diffusion, making the tetramer the fastest diffusing species on the terrace at all temperatures for which diffusion is active (above 200 K). Neutral oxygen vacancy sites are found to bind Pd monomers with a 2.63 eV stronger binding energy than the terrace. A second Pd atom, however, binds to this trapped monomer with a smaller energy of 0.56 eV, so that dimers at defects dissociate on a time scale of milliseconds at room temperature. Larger clusters bind more strongly at defects. Trimers and tetramers dissociate from monomer-bound-defects at elevated temperatures of ca. 600 K. These species are also mobile on the terrace, suggesting they are important for the ripening observed at ?600 K during Pd vapor deposition on MgO(1 0 0) by Haas et al. [G. Haas, A. Menck, H. Brune, J.V. Barth, J.A. Venables, K. Kern, Phys. Rev. B 61 (2000) 11105].