Effect of hydrogenation on B/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 the energetics of substitutional boron on the generic Si(0 0 1)-(1 × 2) surface. For a single B atom substitution corresponding to 0.5 ML coverage, we have considered two different sites: (i) the mixed Si-B dimer structure and (ii) boron substituting for the second-layer Si to form Si-B back-bond structure, which is energetically more favorable than the mixed Si-B dimer by 0.1 eV/dimer. However, when both of these cases are passivated by hydrogen atoms, the situation is reversed and the Si-B back-bond case becomes 0.1 eV/dimer higher in energy than the mixed Si-B dimer case. For the B incorporation corresponding to 1 ML coverage, among the substitutional cases, 100% interdiffusion into the third layer of Si and 50% interdiffusion into the second layer of Si are energetically similar and more favorable than the other cases that are considered. However, when the surface is passivated with hydrogen, the B atoms energetically prefer to stay at the third layer of the Si substrate.     

Methanol adsorption on silicon (0 0 1)

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

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.     

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

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

Density functional study of oxygen adsorption on the Mg3Nd (0 0 1) surface

The adsorption of oxygen atoms on Mg₃Nd (0 0 1) surface was studied based on density function theory (DFT), in which the exchange-correlation potential was chosen as the generalized gradient approximation (GGA) in the Perdew and Wang (PW91). The most preferred adsorption position was at the top-hollow site. Upon the optimization on top-hollow site with different coverage, it was found that the adsorption energy decreased with oxygen coverage. The density of states analysis showed that obvious charge transfer took place between O atom and the nearest Nd atom and chemical bond formed between O atom and the nearest Nd atom after O adsorption. The result of surface energy as a function of chemical potential change of oxygen indicated the clean Mg₃Nd (0 0 1) surface was easy to adsorb oxygen and form 1.00 ML surface.     

The role of carbon impurities on the Si(0 0 1)-c(4 × 4) surface reconstruction: Theoretical calculations

In this work we employ the state-of-the-art pseudopotential method, within a generalized gradient approximation to the density functional theory, combined with a recently developed method for the calculation of HREELS spectra to study a series of different proposed models for carbon incorporation on the silicon (0 0 1) surface. A fully discussion on the geometry, energetics and specially the comparison between experimental and theoretical STM images and electron energy loss spectra indicate that the Si(1 0 0)-c(4 × 4) is probably induced by Si-C surface dimers, in agreement with recent experimental findings.     

A DFT + U description of oxygen vacancies at the TiO2 rutile (1 1 0) surface

Experimental observations indicate that removing bridging oxygen atoms from the TiO₂ rutile (1 1 0) surface produces a localised state approximately 0.7 eV below the conduction band. The corresponding excess electron density is thought to localise on the pair of Ti atoms neighbouring the vacancy; formally giving two Ti³⁺ sites. We consider the electronic structure and geometry of the oxygen deficient TiO₂ rutile (1 1 0) surface using both gradient-corrected density functional theory (GGA DFT) and DFT corrected for on-site Coulomb interactions (GGA + U) to allow a direct comparison of the two methods. We show that GGA fails to predict the experimentally observed electronic structure, in agreement with previous uncorrected DFT calculations on this system. Introducing the +U term encourages localisation of the excess electronic charge, with the qualitative distribution depending on the value of U. For low values of U (?4.0 eV) the charge localises in the sub-surface layers occupied in the GGA solution at arbitrary Ti sites, whereas higher values of U (?4.2 eV) predict strong localisation with the excess electronic charge mainly on the two Ti atoms neighbouring the vacancy. The precise charge distribution for these larger U values is found to differ from that predicted by previous hybrid-DFT calculations.     

A theoretical study of hydrogen adsorption and diffusion on a W(1 1 0) surface

We have used density functional theory to investigate hydrogen adsorption and diffusion on a W(1 1 0) surface. Hydrogen adsorption structures were examined from low coverage to one monolayer, and a threefold hollow site was found to be the most stable site at all coverages. In contrast to previous assertions, the work function decrease is not due to electron transfer from the hydrogen atoms to the W surface, but due to electron depletion at the vacuum region above the hydrogen atoms. Hydrogen atoms can diffuse via short-bridge sites and long-bridge sites at a coverage of θ = 1.0. Although the calculated activation energy for hydrogen diffusion via a short-bridge site is as small as 0.05 eV, field emission microscope experiments have shown that the activation energy for hydrogen diffusion is about 0.20 eV, which agrees fairly well with our calculated value of the activation energy via a long-bridge site. This discrepancy can be related to hydrogen delocalization on the W(1 1 0) surface, which has been suggested by electron energy loss spectroscopy experiments.     

Structural stability and electronic structure of iron adsorption on the GaN(0 0 0 1) surface

In this work, we have investigated by means of first-principles spin-polarized calculations, the electronic and magnetic properties of iron (Fe) adsorption and diffusion on the GaN(0 0 0 1) surface using density functional theory (DFT) within a plane-wave pseudopotential scheme. In the surface adsorption study, results show that the most stable positions of a Fe adatom on GaN(0 0 0 1) surface are the H3 sites and T4 sites, for low and high Fe coverage respectively. We found that the Fe-H3 2 × 2 surface reconstruction exhibits a half-metallic behavior with a spin band gap and stable ferromagnetism ordering, which is a desirable property for high-efficiency magnetoelectronic devices. In addition, confirming previous experimental results, we found that the iron monolayers present a ferromagnetic order and a large thermal stability. This is interesting from a theoretical point of view and for its technological applications.     

Atomic and electronic structures of thin NaCl films grown on a Ge(0 0 1) surface

Density functional theory (DFT) with LDA and GGA have been employed to study the interface and thin film properties of NaCl on a Ge(0 0 1) surface. The atomic and electronic structures of thin NaCl films from one to ten monolayers were analyzed. The layer adsorption energies show that a quasi-crystalline (0 0 1) fcc NaCl film is built up via a layer-by-layer growth mode with NaCl thickness above 2 ML. Simulated STM images show a well-resolved (1 × 1) NaCl atomic structure for sample bias voltage V_s < −2.5 V and the bright protrusions should be assigned to the Cl⁻ ions of the NaCl film. The Ge substrate dimer is reserved and buckled like a clean Ge(0 0 1)-p(2 × 2) surface as the result of weak interface interaction between the dangling bonds coming from valence π states of the Ge substrate and the 3p states of the interfacial Cl⁻ ion. These results are consistent with the experiments of STM, LEED and EELS.     

Density functional study of hypophosphite adsorption on Ni (1 1 1) and Cu (1 1 1) surfaces

Surface structures and electronic properties of hypophosphite, H₂PO₂⁻, molecularly adsorbed on Ni(1 1 1) and Cu(1 1 1) surfaces are investigated in this work by density functional theory at B3LYP/6-31++g(d, p) level. We employ a four-metal-atom cluster as the simplified model for the surface and have fully optimized the geometry and orientation of H₂PO₂⁻ on the metal cluster. Six stable orientations have been discovered on both Ni (1 1 1) and Cu (1 1 1) surfaces. The most stable orientation of H₂PO₂⁻ was found to have its two oxygen atoms interact the surface with two PO bonds pointing downward. Results of the Mulliken population analysis showed that the back donation from 3d orbitals of the transition metal substrate to the unfilled 3d orbital of the phosphorus atom in H₂PO₂⁻ and 4s orbital's acceptance of electron donation from one lone pair of the oxygen atom in H₂PO₂⁻ play very important roles in the H₂PO₂⁻ adsorption on the transition metals. The averaged electron configuration of Ni in Ni₄ cluster is 4s^0.634p^0.023d^9.35 and that of Cu in Cu₄ cluster is 4s^1.004p^0.033d^9.97. Because of this subtle difference of electron configuration, the adsorption energy is larger on the Ni surface than on the Cu surface. The amount of charge transfers due to above two donations is larger from H₂PO₂⁻ to the Ni surface than to the Cu surface, leading to a more positively charged P atom in Ni_nH₂PO₂⁻ than in Cu_nH₂PO₂⁻. These results indicate that the phosphorus atom in Ni_nH₂PO₂⁻ complex is easier to be attacked by a nucleophile such as OH⁻ and subsequent oxidation of H₂PO₂⁻ can take place more favorably on Ni substrate than on Cu substrate.     

First-principles studies of the oxygen adsorption on unreconstructed and reconstructed Ni(1 1 0) surfaces

First-principles calculations have been performed to investigate the adsorption of oxygen on unreconstructed and reconstructed Ni(1 1 0) surfaces. The energetics, structural, electronic and magnetic properties are given in detail. For oxygen adsorption on unreconstructed surface, (n×1)(n=2,3) substrate with oxygen atom on short-bridge site is found to be the most stable adsorption configuration. Whereas energetically most favorable adsorption phase of reconstructed surface is p(n×1) substrate with oxygen atom located at long-bridge site. Our calculations suggest that the surface reconstruction is induced by the oxygen adsorption. We also find there are redistributions of electronic structure and electron transfer from the substrate to adsorbate. Our calculations also indicate surface magnetic moment is enhanced on clean surfaces and oxygen atoms are magnetized weakly after oxygen adsorption. Interestingly, adsorption on unreconstructed surface does not change surface magnetic moment. However, adsorbate leads to reduction of surface magnetic moment in reconstructed system remarkably.     

Function of subsurface boron on Si(0 0 1)-2 × 1: water adsorption

Using density functional theory calculations we investigate the function of subsurface boron in determining surface properties of Si(0 0 1). To demonstrate its effect on surface reactivity we compare the behaviors of water adsorption on the clean and B-modified surfaces. We find that subsurface boron brings about a significant change in surface chemical properties by altering charge polarization of Si(0 0 1) locally. As a consequence, water adsorption on the B-modified surface shows a distinctively different feature from that on the clean surface.     

Nitrogen adsorption and thin surface nitrides on Cu(1 1 1) from first-principles

Experimental studies of nitrogen adsorbed on a Cu(1 1 1) surface show that the surface layer undergoes a reconstruction to form a pseudo-(1 0 0) structure. We use ab initio techniques to demonstrate the theoretical stability of this reconstructed surface phase over a range of conditions. We systematically investigate the chemisorption of N on the Cu(1 1 1) surface, from 0.06 to 1 ML coverage. A peculiar atomic relaxation of N atoms for 0.75 ML is identified, which results in the formation of a (metastable) “N-trimer cluster” on the surface. We have also investigated surface nitride formation, as suggested from experiments. A surface nitride-like structure similar to the reported pseudo-(1 0 0) reconstruction is found to be highly energetically favored. Using concepts from “ab initio atomistic thermodynamics”, we predict that this surface nitride exists for a narrow range of nitrogen chemical potential before the formation of bulk Cu₃N.     

Surface core level shift observed on NiAl(1 1 0)

Using high resolution core level spectroscopy, a surface core level shift towards lower binding energy of −0.13 eV is determined for the 2p level of the outwardly relaxed Al surface atoms on NiAl(1 1 0). Density functional theory based calculations with inclusion of final state effects yield a value of −0.14 eV for this shift in excellent agreement with experiment. We show that the initial state approximation yields a value of +0.09 eV, i.e. the inclusion of final state relaxation effects is vital not only to obtain the correct value but even the correct sign for this shift.     

DFT study on dissociative adsorption of SiH4 and GeH4 on SiGe(1 0 0)-2 × 1 surface

SiH₄ and GeH₄ dissociative adsorptions on a buckled SiGe(1 0 0)-2 × 1 surface have been analyzed using density functional theory (DFT) at the B3LYP level. The Ge alloying in the Si(1 0 0)-2 × 1 surface affects the dimer buckling and its surface reactivity. Systematic Ge influences on the reaction energetics are found in SiH₄ and GeH₄ reactions with four dimers of Si^∗-Si, Ge^∗-Si, Ge^∗-Ge, and Si^∗-Ge (∗ denotes the protruded atom). On a half H-covered surface, the energy barriers for silane and germane adsorption are higher than those on the pristine surface. The energy barrier for silane adsorption is higher than the corresponding barrier for germane adsorption. Rate constants are also calculated using the transition-state theory. We conclude that the SiGe surface reactivity in adsorption reaction depends on the Ge presence in dimer form. If the surface Ge is present in form of Ge^∗-Ge, the surface reactivity decreases as the Ge^∗-Ge content increases. If the surface Ge prefers to be in form of Ge^∗-Si at low Ge contents, the surface reactivity increases first, then decreases at high Ge surface contents when Ge^∗-Ge prevails. The calculated rate constant ratio of GeH₄ adsorption on Si^∗-Si over Ge^∗-Ge at 650 °C is 2.1, which agrees with the experimental ratio of GeH₄ adsorption probability on Si(1 0 0) over Si(1 0 0) covered by one monolayer Ge. The experimental ratio is 1.7 measured through supersonic molecular beam techniques. This consistency between calculation and experimental results supports that one monolayer of Ge on Si(1 0 0) exists in form of Ge^∗-Ge dimer.     

Water adsorption and dissociation on BeO(0 0 1) and (1 0 0) surfaces

Plateaus in water adsorption isotherms on hydroxylated BeO surfaces suggest significant differences between the hydroxylated (1 0 0) and (0 0 1) surface structures and reactivities. Density functional theory structures and energies clarify these differences. Using relaxed surface energies, a Wulff construction yields a prism crystal shape exposing long (1 0 0) sides and much smaller (0 0 1) faces. This is consistent with the BeO prisms observed when beryllium metal is oxidized. A water oxygen atom binds to a single surface beryllium ion in the preferred adsorption geometry on either surface. The water oxygen/beryllium bonding is stronger on the surface with greater beryllium atom exposure, namely the less-stable (0 0 1) surface. Water/beryllium coordination facilitates water dissociation. On the (0 0 1) surface, the dissociation products are a hydroxide bridging two beryllium ions and a metal-coordinated hydride with some surface charge depletion. On the (1 0 0) surface, water dissociates into a hydroxide ligating a Be atom and a proton coordinated to a surface oxygen but the lowest energy water state on the (1 0 0) surface is the undissociated metal-coordinated water. The (1 0 0) fully hydroxylated surface structure has a hydrogen bonding network which facilitates rapid proton shuffling within the network. The corresponding (0 0 1) hydroxylated surface is fairly open and lacks internal hydrogen bonding. This supports previous experimental interpretations of the step in water adsorption isotherms. Further, when the (1 0 0) surface is heated to 1000 K, hydroxides and protons associate and water desorbs. The more open (0 0 1) hydroxylated surface is stable at 1000 K. This is consistent with the experimental disappearance of the isotherm step when heating to 973 K.     

CO methanation on Ni(1 1 1) and modified Ni3Al(1 1 1) surfaces: A first-principle study

The adsorption energies of intermediates in CO methanation on the modified Ni₃Al(1 1 1) surface and the Ni(1 1 1) surface are calculated using density functional theory. A microkinetic analysis based on the calculated adsorption energies is performed to explain the different kinetics of CO methanation catalyzed by Ni₃Al and Ni powders. The electronic structures of different atoms on the modified Ni₃Al surface are also presented, and correlate well with the adsorption energies and geometries.     

A DFT periodic study of the vanadyl pyrophosphate (1 0 0) surface

The unique ability of the vanadyl pyrophosphate (1 0 0) surface to activate n-butane and then selectively oxidise the hydrocarbon to maleic anhydride was investigated using modern quantum chemical methods. Bulk (VO)₂P₂O₇, together with stoichiometric and phosphorus-enriched (1 0 0) surfaces, were analysed using periodic density functional theory calculations. Also simulated was surface ionic relaxation from bulk geometry, and surface hydration. Density of states (DOS) plots show that, whether stoichiometric or phosphorus-enriched, bulk terminated or relaxed, bare or hydrated, local covalent reactivity at the (1 0 0) surface is controlled by vanadium species. Terminal P-O oxygen species are the most nucleophilic surface oxygens, as indicated by their predominance of sub-vanadium high-lying valence band levels. A periodic treatment of (VO)₂P₂O₇(1 0 0) hence gives results qualitatively identical to those obtained from earlier cluster calculations. Simulation of surface ionic relaxation shows that in-plane P-O-V oxygens may also be involved in rupture of substrate C-H bonds for mild oxidation, while surface hydration calculations indicate that dissociative chemisorption of water may play a key role in perpetuation of the selective oxidation cycle.     

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.