NH3 on Si(0 0 1): Can Gaussian cluster and planewave slab models agree on energetics?

We critically evaluate the use of cluster and periodic slab models in describing the NH₃/Si(0 0 1) molecule-surface reaction system. We show that considerable discrepancies in the relative adsorbate energetics originate in the limitations of the small cluster and slab models commonly used. These limitations in turn are the consequence of the balance that must be struck between the competing demands of cluster/slab size, basis set size and exchange-correlation model. This leads us to consider “cluster compound models” in which the results of several smaller calculations (separately probing the effects of cluster/slab size, basis set size and exchange-correlation model) are combined to estimate the energy of a converged model.     

Study of vacancy on diamond (1 0 0) (2 × 1) surface from first-principles

Structure and energy related properties of neutral and charged vacancies on relaxed diamond (1 0 0) (2 × 1) surface were investigated by means of density functional theory. Calculations indicate that the diffusion of a single vacancy from the top surface layer to the second layer is not energetically favored. Analysis of energies in charged system shows that neutral state is most stable on diamond (1 0 0) (2 × 1) surface. The multiplicity of possible states can exist on diamond (1 0 0) surface in dependence on the surface Fermi level, which supports that surface diffusion of a vacancy is mediated by the change of vacancy charge states. Analysis of density of states shows surface vacancy can be effectively measured by photoelectricity technology.     

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

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

A quantum-mechanical study of the vinyl fluoride adsorbed on the rutile TiO2(1 1 0) surface

The adsorption of vinyl fluoride on the rutile TiO₂(1 1 0) surface has been simulated, on the basis of a recently proposed experimental model, using hybrid-exchange density functional theory. Different surface coverages have been considered and the lateral interaction between adsorbed vinyl fluoride molecules has been quantified through a simple model of nearest and next nearest neighbouring molecules. The vibrational frequencies of the adsorbed molecule have been calculated and are found to be in excellent agreement with those observed providing support for the proposed adsorption model. The effect of the adsorption on the electronic structure of the molecule and the surface have been characterised by computing electrostatic potential maps and the local density of states.     

Optical anisotropy of the In/Si(1 1 1)(4 × 1)/(8 × 2) nanowire array

Ab initio calculations of the reflectance anisotropy of Si(1 1 1)-In surfaces are presented. A very pronounced optical anisotropy around 2 eV is found that is related to In-chain states. The distortion of the indium chains characteristic for the (4 × 1) → (8 × 2) phase transition results in a splitting of the 2 eV peak, as observed experimentally. The splitting occurs irrespective wether the phase transition occurs according to the trimer or hexamer model.     

The dissociative adsorption of borane on the Ge(1 0 0)-2 × 1 surface: A density functional theory study

By means of cluster models coupled with density functional theory, we have studied the hydroboration of the Ge(1 0 0)-2 × 1 surface with BH₃. It was found that the Ge(1 0 0) surface exhibits rather different surface reactivity toward the dissociative adsorption of BH₃ compared to the C(1 0 0) and Si(1 0 0) surfaces. The strong interaction still exists between the as-formed BH₂ and H adspeices although the dissociative adsorption of BH₃ on the Ge(1 0 0) surface occurs readily, which is in distinct contrast to that on the C(1 0 0) and Si(1 0 0) surfaces. This can be understood by the electrophilic nature of the down Ge atom, which makes it unfavourable to form a GeH bond with the dissociating proton-like hydrogen. Alternatively, it can be attributed to the weak proton affinity of the Ge(1 0 0) surface. Nevertheless, the overall dissociative adsorption of BH₃ on group IV semiconductor surfaces is favourable both thermodynamically and kinetically, suggesting the interesting analogy and similar diversity chemistry of solid surface in the same group.     

A density functional theory study of ethylene adsorption on Ni10(1 1 1), Ni13(1 0 0) and Ni10(1 1 0) surface cluster models and Ni13 nanocluster

Ethylene adsorption was studied by use of DFT/B3LYP with basis set 6-31G(d,p) in Gaussian’03 software. It was found that ethylene has adsorbed molecularly on all clusters with π adsorption mode. Relative energy values were calculated to be −50.86 kcal/mol, −20.48 kcal/mol, −32.44 kcal/mol and −39.27 kcal/mol for Ni₁₃ nanocluster, Ni₁₀(1 1 1), Ni₁₃(1 0 0) and Ni₁₀(1 1 0) surface cluster models, respectively. Ethylene adsorption energy is inversely proportional to Ni coordination number when Ni₁₀(1 1 1), Ni₁₃(1 0 0) and Ni₁₀(1 1 0) cluster models and Ni₁₃ nanocluster are compared with each other.     

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

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

First-principle study of NO adsorption on the LaFeO3 (0 1 0) surface

The adsorption of NO molecule on the LaFeO₃ (0 1 0) surface was studied using first-principle calculations based on density functional theory. The calculated results indicate that the Fe-top site is the most favorable for NO adsorption. The N-O bond length, Mulliken charge, and the N-O vibration frequency of the NO molecule are discussed after adsorption. The analysis results of the density of the states show that when NO is adsorbed with the Fe-NO configuration, the bonding mechanism is mainly from the interaction between the NO and the Fe d orbit.     

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

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

Density functional theory calculations of surface properties and H2 adsorption on the Cu2O (1 1 1) surface

First-principles calculation on the basis of the density functional theory (DFT) and generalized gradient approximation have been applied to study the adsorption of H₂ on the stoichiometric O-terminated Cu₂O (1 1 1), Cu₂O (1 1 1)-Cu_CUS and Cu-terminated Cu₂O (1 1 1) surfaces. The optimal adsorption position and orientation of H₂ on the stoichiometric O-terminated Cu₂O (1 1 1) surface and Cu-terminated Cu₂O (1 1 1) surface were determined and electronic structural changes upon adsorption were investigated by calculating the Local Density of States (LDOS) of the Cu_CUS 3d and Cu_CUS 4s of stoichiometric O-terminated Cu₂O (1 1 1) surface. These results showed that H₂ molecule adsorption on Cu_CUS site parallel to stoichiometric O-terminated Cu₂O (1 1 1) surface and H₂ molecule adsorption on Cu₂ site parallel to Cu-terminated Cu₂O (1 1 1) surface were the most favored, respectively. The presence of surface copper vacancy has a little influence on the structures when H₂ molecule adsorbs on Cu_CSA, O_CUS and O_CSA atoms and the H₂ molecule is only very weakly bound to the Cu₂O (1 1 1)-Cu_CUS surface. From the analysis of stoichiometric O-terminated Cu₂O (1 1 1) Local Density of States, it is observed that Cu_CUS 3d orbital has moved to a lower energy and the sharp band of Cu_CUS 4s is delocalized when compared to that before H₂ molecule adsorption, and overlapped substantially with bands due to adsorbed H₂ molecule. The Mulliken charges of H₂ adsorption on Cu_CUS site showed that H₂ molecule obtained electron from Cu_CUS which was consistent with the calculated electronic structural changes upon H₂ adsorption.     

CO2 adsorption on Cr(1 1 0) and Cr2O3(0 0 0 1)/Cr(1 1 0)

We attempt to correlate qualitatively the surface structure with the chemical activity for a metal surface, Cr(1 1 0), and one of its surface oxides, Cr₂O₃(0 0 0 1)/Cr(1 1 0). The kinetics and dynamics of CO₂ adsorption have been studied by low energy electron diffraction (LEED), Aug er electron spectroscopy (AES), and thermal desorption spectroscopy (TDS), as well as adsorption probability measurements conducted for impact energies of E_i = 0.1-1.1 eV and adsorption temperatures of T_s = 92-135 K. The Cr(1 1 0) surface is characterized by a square shaped LEED pattern, contamination free Cr AES, and a single dominant TDS peak (binding energy E_d = 33.3 kJ/mol, first order pre-exponential 1 × 10¹³ s⁻¹). The oxide exhibits a hexagonal shaped LEED pattern, Cr AES with an additional O-line, and two TDS peaks (E_d = 39.5 and 30.5 kJ/mol). The initial adsorption probability, S₀, is independent of T_s for both systems and decreases exponentially from 0.69 to 0.22 for Cr(1 1 0) with increasing E_i, with S₀ smaller by ∼0.15 for the surface oxide. The coverage dependence of the adsorption probability, S(Θ), at low E_i is approx. independent of coverage (Kisliuk-shape) and increases initially at large E_i with coverage (adsorbate-assisted adsorption). CO₂ physisorbs on both systems and the adsorption is non-activated and precursor mediated. Monte Carlo simulations (MCS) have been used to parameterize the beam scattering data. The coverage dependence of E_d has been obtained by means of a Redhead analysis of the TDS curves.     

Adsorption and dissociation of O2 on CuCl(1 1 1) surface: A density functional theory study

The adsorption and dissociation of O₂ on CuCl(1 1 1) surface have been systematically studied by the density functional theory (DFT) slab calculations. Different kinds of possible modes of atomic O and molecular O₂ adsorbed on CuCl(1 1 1) surface and possible dissociation pathways are identified, and the optimized geometry, adsorption energy, vibrational frequency and Mulliken charge are obtained. The calculated results show that the favorable adsorption occurs at hollow site for O atom, and molecular O₂ lying flatly on the surface with one O atom binding with top Cu atom is the most stable adsorption configuration. The O-O stretching vibrational frequencies are significantly red-shifted, and the charges transferred from CuCl to oxygen. Upon O₂ adsorption, the oxygen species adsorbed on CuCl(1 1 1) surface mainly shows the characteristic of the superoxo (O₂⁻), which primarily contributes to improving the catalytic activity of CuCl, meanwhile, a small quantity of O₂ dissociation into atomic O also occur, which need to overcome very large activation barrier. Our results can provide some microscopic information for the catalytic mechanism of DMC synthesis over CuCl catalyst from oxidative carbonylation of methanol.     

Comparison of S, Pt, and Hf adsorption on NiAl(1 1 0)

First-principles periodic slab density-functional theory (DFT) calculations with a plane-wave basis are used to predict the properties of S, Pt, and Hf adsorption on NiAl(1 1 0). Stable adsorption sites are identified, and adsorbate binding energies and structures are predicted. We find that while S adsorbs in a threefold site, the metals prefer to adsorb in the Ni-Ni twofold bridge site. The latter finding is consistent with scanning tunneling microscopy experiments for adsorption of various transition metals on NiAl(1 1 0) by Ho and coworkers. S is predicted to easily diffuse between threefold sites. We find that Pt and Hf both induce significant changes in the local surface structure, changing twofold bridge sites into fourfold coordination sites by drawing next-nearest-neighbor atoms nearly equidistant with the nearest-neighbor atoms. We find Pt favors interaction with Al slightly more than Ni, while Hf shows a particularly strong affinity for Ni compared to Al. We also predict that Hf may diffuse one-dimensionally along Ni rows with a barrier of ≈0.6 eV.     

Adsorption and reactivity of SO2 on Ir(1 1 1) and Rh(1 1 1)

The adsorption and reactivity of SO₂ on the Ir(1 1 1) and Rh(1 1 1) surfaces were studied by surface science techniques. X-ray photoelectron spectroscopy measurements showed that SO₂ was molecularly adsorbed on both the Ir(1 1 1) surface and the Rh(1 1 1) surface at 200 K. Adsorbed SO₂ on the Ir(1 1 1) surface disproportionated to atomic sulfur and SO₃ at 300 K, whereas adsorbed SO₂ on the Rh(1 1 1) surface dissociated to atomic sulfur and oxygen above 250 K. Only atomic sulfur was present on both surfaces above 500 K, but the formation process and structure of the adsorbed atomic sulfur on Ir(1 1 1) were different from those on Rh(1 1 1). On Ir(1 1 1), atomic sulfur reacted with surface oxygen and was completely removed from the surface, whereas on Rh(1 1 1), sulfur did not react with oxygen.     

A study of the interaction between perylene and the TiO2(1 1 0)-(1 × 1) surface-based on XPS, UPS and NEXAFS measurements

The interaction between a semi-large aromatic hydrocarbon compound (perylene) and the TiO₂(1 1 0)-(1 × 1) surface under ultra high vacuum conditions has been probed by X-ray photoemission spectroscopy (XPS), ultraviolet photoemission spectroscopy (UPS) and near-edge X-ray absorption fine structure (NEXAFS) methods. UPS measurements of the adsorbate system have been compared with an experimental UPS spectrum of perylene in the gas phase and a calculated spectrum obtained by means of density functional theory (DFT) methods. NEXAFS results of perylene molecules adsorbed on TiO₂(1 1 0)-(1 × 1) were compared with data from an α-phase perylene single crystal. A novel analysis of the valence data has been employed to show that no strong chemical interaction takes place between perylene and the TiO₂(1 1 0)-(1 × 1) surface. Furthermore, angle-dependent NEXAFS measurements and the growth curve results suggest that the perylene molecules are oriented flat down onto the TiO₂ substrate due to weak van der Waals interactions.     

Low temperature adsorption of oxygen on reduced V2O3(0 0 0 1) surfaces

Well ordered V₂O₃(0 0 0 1) films were prepared on Au(1 1 1) and W(1 1 0) substrates. These films are terminated by a layer of vanadyl groups under typical UHV conditions. Reduction by electron bombardment may remove the oxygen atoms of the vanadyl layer, leading to a surface terminated by vanadium atoms. The interaction of oxygen with the reduced V₂O₃(0 0 0 1) surface has been studied in the temperature range from 80 to 610 K. Thermal desorption spectroscopy (TDS), infrared reflection absorption spectroscopy (IRAS), high resolution electron energy loss spectroscopy (HREELS), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) were used to study the adsorbed oxygen species. Low temperature adsorption of oxygen on reduced V₂O₃(0 0 0 1) occurs both dissociatively and molecularly. At 90 K a negatively charged molecular oxygen species is observed. Upon annealing the adsorbed oxygen species dissociates, re-oxidizing the reduced surface by the formation of vanadyl species. Density functional theory was employed to calculate the structure and the vibrational frequencies of the O₂ species on the surface. Using both cluster and periodic models, the surface species could be identified as η²-peroxo () lying flat on surface, bonded to the surface vanadium atoms. Although the O-O vibrational normal mode involves motions almost parallel to the surface, it can be detected by infrared spectroscopy because it is connected with a change of the dipole moment perpendicular to the surface.     

Nanostructure formation aided by self-organised Bi nanolines on Si(0 0 1)

We provide a mini review of recent theoretical investigations of nanostructure formation aided by self-organised Bi nanolines on the Si(0 0 1) surface. It is suggested that hydrogen-passivated single-domain Si(0 0 1) produced by the formation of defect-free, hundreds of nm long, and 1.2-1.5 nm wide Bi nanolines provides an appealing template with preferential sites for adsorption of other elements. Based on ab initio pseudopotential calculations it is suggested that using the Bi nanoline template it should be possible to grow the following structures: mixed Ge-Si dimer structures on the Si(0 0 1) terrace between two neighbouring Bi nanolines; small In nanoclusters along the nanoline; and line and cluster structures of Fe atoms with novel electronic and magnetic properties.     

CO2 adsorption on the bimetallic Zn-on-Cu(1 1 0) system

Adsorption probability measurements (molecular beam scattering) have been conducted to examine the adsorption dynamics (i.e. the gas-surface energy transfer processes) of CO₂ adsorption on the Zn-on-Cu(1 1 0) bimetallic system. The results indicate surface alloy formation, which is in agreement with prior studies. Depositing Zn at 300 K on Cu(1 1 0), above the condensation temperature of CO₂, leads to a “blocking” of CO₂ adsorption sites by Zn which is incorporated in the Cu(1 1 0) surface. This apparent site blocking effect indicates a lowering of the CO₂ binding energy on the alloyed surface as compared with the clean Cu(1 1 0) support. The Zn coverage has been calibrated by Auger electron spectroscopy and thermal desorption spectroscopy.     

The surface structure of BaO on Pt(1 1 1): (2 × 2)-reconstructed BaO(1 1 1)

The surface structure of BaO(1 1 1) has been determined using STM and computer modelling. The BaO(1 1 1) surface was prepared in thin film form on Pt(1 1 1) and presents a surface with twice the lattice parameter expected for that of the bulk termination, i.e. a (2 × 2) reconstruction. Computer modelling indicates that the bulk termination is unstable, but that the (2 × 2) reconstructed BaO(1 1 1) surface has a low surface energy and is hence a stable surface reconstruction. The (2 × 2) reconstruction consists of small, three-sided pyramids with (1 0 0) oriented sides and either oxygen or barium ions at the apices. Less regular surface reconstructions containing the same pyramids are almost equally stable, indicating that we may also expect less regular regions to appear with a fairly random distribution of these surface species. The simulations further suggest that a regular (4 × 4) reconstruction built up of bigger pyramids is even more energetically favourable, and some evidence is found for such a structure in the STM.