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.     

Adsorption site preference of CO2 on the Pt(1 0 0) surface by ab initio calculations

This paper investigates the adsorption sites and electronic structure of the adsorption of CO₂ on the Pt(1 0 0) surface by ab initio periodic density functional theory (DFT) methods. Several parallel and vertical adsorption sites are examined in detail. The results show that the adsorption energy for the atop hollow horizontal (thh) site is 0.34 eV. However, other sites only have small binding energies and these values are very close. For an atop hollow horizontal site, the calculated elecronic interaction is contributed to not only the Pt-O atoms, but also Pt-C atoms. So the results indicate that the thh site is the most favorable and stable.     

Formation of oxygen vacancies on the TiO2(1 1 0) surfaces

Density functional theory was employed to investigate the formation and properties of the oxygen vacancies on the rutile TiO₂(1 1 0) surface. It is found that the formation of the positively charged bridging-oxygen vacancy (BOV⁺, 4.2 eV) is the most favored one, followed by the positively charged in-plane-oxygen vacancy (POV⁺, 4.5 eV). In contrast, the formation of the neutral bridging-oxygen and in-plane-oxygen vacancies (BOV and POV), and their dication oxygen vacancies (BOV²⁺ and POV²⁺) needs much higher energies (7.9 and 8.3 vs. 8.1 and 8.6 eV), respectively.     

First-principles study for the atomic structures and electronic properties of PbTiO3 oxygen-vacancies (0 0 1) surface

The structural and electronic properties of fully-relaxed PbTiO₃ (0 0 1) oxygen-vacancy surface with PbO and TiO₂ terminations are investigated by first-principles calculations. In contrast to the perfect surface, the smaller surface rumples and interlayer distances have been found. The largest relaxation occurs on the second layer atoms not on the surface layer ones, and some in-gap Ti 3d states at about −1.1 eV below the Fermi-level are observed in the TiO₂-terminated surface caused by oxygen-vacancies. For the PbO-terminated surface, some in-gap Ti 3d states and Pb 6p states also move into the bulk midgap region to become partially occupied, and two different chemical states of the Pb 6s states were found. One is attributed to the bulk perovskite Pb atoms and another one is caused by the relaxation of surface Pb atoms. These theoretical results would give a good reference for the future experimental studies.     

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.     

Growth mode and electronic structure of Au nano-clusters on NiO(0 0 1) and TiO2(1 1 0)

The growth mode and electronic structure of Au nano-clusters grown on NiO and TiO₂ were analyzed by reflection high-energy electron diffraction, a field-emission type scanning electron microscope, medium energy ion scattering and photoelectron spectroscopy. Au was deposited on clean NiO(0 0 1)-1 × 1 and TiO₂(1 1 0)-1 × 1 surfaces at room temperature with a Knudsen cell at a rate of 0.25-0.35 ML/min (1 ML = 1.39 × 10¹⁵ atoms/cm²:Au(1 1 1)). Initially two-dimensional (2D) islands with thickness of one Au-atom layer grow epitaxially on NiO(0 0 1) and then neighboring 2D-islands link each other to form three-dimensional (3D)-islands with the c-axis oriented to the [1 1 1] direction. The critical size to form 3D-islands is estimated to be about 5 nm². The shape of the 3D-islands is well approximated by a partial sphere with a diameter d and height h ranging from 2.0 to 11.8 nm and from 0.95 to 4.2 nm, respectively for Au coverage from 0.13 to 4.6 ML. The valence band spectra show that the Au/NiO and Au/TiO₂ surfaces have metallic characters for Au coverage above 0.9 ML. We observed Au 4f spectra and found no binding energy shift for Au/NiO but significant higher binding energy shifts for Au/TiO₂ due to an electron charge transfer from Au to TiO₂. The work function of Au/NiO(0 0 1) gradually increases with increase in Au coverage from 4.4 eV (NiO(0 0 1)) to 5.36 eV (Au(1 1 1)). In contrast, a small Au deposition(0.15 to 1.5 ML) on TiO₂(1 1 0) leads to reduction of the work function, which is correlated with an electron charge transfer from Au to TiO₂ substrate.     

Electronic properties of Au nano-particles supported on stoichiometric and reduced TiO2(1 1 0) substrates

Growth modes and electronic properties were analyzed for Au nano-particles grown on stoichiometric and reduced TiO₂(1 1 0) substrates by medium energy ion scattering (MEIS) and photoelectron spectroscopy(PES) using synchrotron-radiation-light. Initially, two-dimensional islands (2D) with a height of one and two atomic layers grow and higher coverage increases the islands height to form three-dimensional (3D) islands for the stoichiometric TiO₂(1 1 0) substrate. In contrast, 3D islands start to grow from initial stage with a small Au coverage (?0.1 ML, 1 ML = 1.39 × 10¹⁵  atoms/cm²: Au(1 1 1)) probably due to O-vacancies acting as a nucleation site. Above 0.7 ML, all the islands become 3D ones taking a shape of a partial sphere and the Au clusters change to metal for both substrates. We observed the Au 4f and Ti 3p core level shifts together with the valence band spectra. The Ti 3p peak for the O-deficient surface shifts to higher binding energy by 0.25 ± 0.05 eV compared to that for the stoichiometric surface, indicating downward band bending by an electron charge transfer from an O-vacancy induced surface state band to n-type TiO₂ substrate. Higher binding energy shifts of Au 4f peaks observed for both substrates reveal an electron charge transfer from Au to TiO₂ substrates. The work functions of Au nano-particles supported on the stoichiometric and reduced TiO₂ substrates were also determined as a function of Au coverage and explained clearly by the above surface and interface dipoles.     

The small unit cell reconstructions of SrTiO3(1 1 1)

We analyze the basic structural units of simple reconstructions of the (1 1 1) surface of SrTiO₃ using density functional calculations. The prime focus is to answer three questions: what is the most appropriate functional to use; how accurate are the energies; what are the dominant low-energy structures and where do they lie on the surface phase diagram. Using test calculations of representative small molecules we compare conventional PBE-GGA with higher-order methods such as the TPSS meta-GGA and on-site hybrid methods PBE0 and TPSSh, the later being the most accurate. There are large effects due to reduction of the metal d oxygen sp hybridization when using the hybrid methods which are equivalent to a dynamical GGA + U, which leads to rather substantial improvements in the atomization energies of simple calibration molecules, even though the d-electron density for titanium compounds is rather small. By comparing the errors of the different methods we are able to generate an estimate of the theoretical error, which is about 0.25 eV per 1 × 1 unit cell, with changes of 0.5-1.0 eV per 1 × 1 cell with the more accurate method relative to conventional GGA. An analysis of the plausible structures reveals a new low-energy TiO₂-rich configuration with octahedral co-ordination. This structure can act as a template for layers of either TiO or Ti₂O₃, consistent with experimental results. The results also suggest that both the fracture surface and the stoichiometric SrTiO₃(1 1 1) surface should spontaneously disproportionate into SrO and TiO₂ rich domains.     

TiO2-rich reconstructions of SrTiO3(0 0 1): a theoretical study of structural patterns

We have recently reported structure solutions for the (2 × 1) and c(4 × 2) reconstructions of SrTiO₃(0 0 1) based on high-resolution electron microscopy, direct methods analysis of diffraction data and density functional theory. Both reconstructions were found to be TiO₂-rich and feature a single overlayer of TiO₂ stoichiometry on top of a bulk-like TiO₂ layer. Qualitatively, the two reconstruction geometries differ in the cation sub-lattice of the overlayer only, where Ti atoms occupy half of the fivefold cation sites. In the present work we use density functional theory to generate a number of variations of this structural motif in search of patterns of stability. We find a reliable predictor for the reconstruction energy in the ability of oxygen atoms to relax vertically out of the overlayer plane to minimize non-bonded oxygen-oxygen repulsions. Out-of-plane relaxation of oxygen atoms in turn is modulated by the number and relative position of coordinating Ti atoms, which yields simple empirical rules as to how cations are distributed in low energy reconstructions.     

Defects at the TiO2(1 0 0) surface probed by resonant photoelectron diffraction

We report photoelectron diffraction (PED) experiments of weakly sub-stoichiometric TiO₂(1 0 0) rutile surfaces. Apart from standard core-level PED from the Ti-2p3/2 line, we have studied valence band PED from the defect induced Ti-3d states in the insulating band gap. For maximum yield, the latter were resonantly excited at the Ti-2p absorption edge. The PED patterns have been analyzed within the forward scattering approximation as well as by comparison with simulated PED patterns obtained in multiple scattering calculations. The analysis shows that the defect induced Ti-3d charge is mainly located on the second layer Ti atoms.     

First-principles study on the half-metallicity of zinc-blende CrP(1 1 0) surface

We have investigated the half-metallicity and magnetism at the (1 1 0) surface of CrP by using the all-electron full-potential linearized augmented plane wave (FLAPW) method within the generalized gradient approximation (GGA). From the calculated local density of states (LDOS), we found that the (1 1 0) surface of CrP preserves the half-metallicity, but the band gaps (∼1.1 eV) of the minority states for the surface Cr and P atoms are much reduced from the bulk value (∼1.9 eV). The magnetic moment of the P is coupled antiferromagnetically to that of the Cr. The magnetic moment of surface Cr atom is calculated to be 3.31μ_B which is increased by 10% compared to the bulk value, 3.00μ_B.     

Electronic structures and optical properties of rutile TiO2 with different point defects from DFT + U calculations

The electronic states and formation energies of four types of lattice point defects in rutile TiO₂ are studied using the first-principles calculations. The existence of oxygen vacancy leads to a deep donor defect level in the forbidden band, while the Ti interstitial forms two local states. It is predicted that oxygen vacancy prefers to combine with Ti-interstitial to form V_O–Ti_i dimer by a partial 3d electron transfer from the Ti_i to its neighboring V_O. The charge distribution between a Ti interstitial and its neighboring Ti ions partially shields the Coulomb interactions. Lastly, optical properties of these defective lattices are discussed.     

The surface stress of the (1 1 0) and (1 0 0) surfaces of rutile and the effect of water adsorbents

The surface stress on clean TiO₂ (1 1 0) and (1 0 0) surfaces, and those with four types of adsorbent - (i) molecularly adsorbed water, (ii) dissociatively adsorbed water, (iii) dissociatively adsorbed water at an oxygen vacancy, and (iv) adsorbed hydrogen - was investigated in the framework of density functional theory using a slab model. The calculations were intended to rationalize the effect of the artificially introduced stress that occurs in experimentally photoinduced hydrophilicity. Tensile stress was observed for a clean (1 1 0) surface, and a mixture of tensile and compressive stress for a clean (1 0 0) surface. The adsorbate-induced surface stresses were analyzed in terms of the sixfold coordinated character of the surface titanium atoms, hydrogen bonds between the adsorbents and the bridging oxygen atoms, and the change in electron density in the vicinity of the surface.     

N incorporation and electronic structure in N-doped TiO2(1 1 0) rutile

We describe the growth and properties of well-defined epitaxial TiO_2−xN_x rutile for the first time. A mixed beam of atomic N and O radicals was prepared in an electron cyclotron resonance plasma source and Ti was supplied from a high-temperature effusion cell or an electron beam evaporator, depending on the required flux. A very high degree of structural quality is generally observed for films grown under optimized anion-rich conditions. N substitutes for O in the lattice, but only at the ∼1 at.% level, and is present as N³⁻. Epitaxial growth of TiO₂ and TiO_2−xN_x rutile prepared under anion-rich conditions is accompanied by Ti indiffusion, leading to interstitial Ti (Ti_i), which is a shallow donor in rutile. Our data strongly suggest that Ti_i donor electrons compensate holes associated with substitutional N²⁻ (i.e., Ti(III) + N²⁻ → Ti(IV) + N³⁻), leading to highly resistive or weakly n-type, but not p-type material. Ti 2p core-level line shape analysis reveals hybridization of N and Ti, as expected for substitutional N. Ti-N hybridized states fall in the gap just above the VBM, and extend the optical absorption well into the visible.     

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.     

The 2 × 1 reconstruction of the rutile TiO2(0 1 1) surface: A combined density functional theory, X-ray diffraction, and scanning tunneling microscopy study

An extensive search for possible structural models of the (2 × 1)-reconstructed rutile TiO₂(0 1 1) surface was carried out by means of density functional theory (DFT) calculations. A number of models were identified that have much lower surface energies than the previously-proposed ‘titanyl’ and ‘microfaceting’ models. These new structures were tested with surface X-ray diffraction (SXRD) and voltage-dependent STM measurements. The model that is (by far) energetically most stable shows also the best agreement with SXRD data. Calculated STM images agree with the experimental ones for appropriate tunneling conditions. In contrast to previously-proposed models, this structure is not of missing-row type; because of its similarity to the fully optimized brookite TiO₂(0 0 1) surface, we call it the ‘brookite (0 0 1)-like’ model. The new surface structure exhibits two different types of undercoordinated oxygen and titanium atoms, and is, in its stoichiometric form, predicted to be rather inert towards the adsorption of probe molecules.     

Island coalescence and diffusion along kinked steps on Cu(0 0 1): Evidence for a large kink Ehrlich-Schwoebel barrier

Using temperature-variable scanning tunneling microscopy, we studied the coalescence of vacancy islands on Cu(0 0 1) in ultra-high vacuum. From the temperature dependence of the relaxation of merged vacancy islands to the equilibrium shape we obtain an activation energy of the island coalescence process of 0.76 eV. From that value we deduce an activation energy for the atomic hopping coefficient of E_Γh=0.89 eV. Comparing our result with previous STM data on step fluctuations with dominant diffusion along straight step segments (E_Γh=0.68 eV; [M. Giesen, S. Dieluweit, J. Mol. Catal. A: Chem. 216 (2004) 263]) and step fluctuations with kink crossing (E_Γh=0.9 eV; [M. Giesen-Seibert, F. Schmitz, R. Jentjens, H. Ibach, Surf. Sci. 329 (1995) 47]), we conclude that there is a large extra barrier for diffusion of atoms across kinks on Cu(0 0 1) of the order of 0.23 eV. This is the first direct experimental evidence for the existence of a large kink Ehrlich-Schwoebel barrier on Cu(0 0 1).     

Molecular dynamics study on surface structure and surface energy of rutile TiO2 (1 1 0)

The formula for surface energy was modified in accordance with the slab model of molecular dynamics (MDs) simulations, and MD simulations were performed to investigate the relaxed structure and surface energy of perfect and pit rutile TiO₂(1 1 0). Simulation results indicate that the slab with a surface more than four layers away from the fixed layer expresses well the surface characteristics of rutile TiO₂ (1 1 0) surface; and the surface energy of perfect rutile TiO₂ (1 1 0) surface converges to 1.801±0.001 J m⁻². The study on perfect and pit slab models proves the effectiveness of the modified formula for surface energy. Moreover, the surface energy of pit surface is higher than that of perfect surface and exhibits an upper-concave parabolic increase and a step-like increase with increasing the number of units deleted along [0 0 1] and [1 1 0], respectively. Therefore, in order to obtain a higher surface energy, the direction along which atoms are cut out should be chosen in accordance with the pit sizes: [] direction for a small pit size and [0 0 1] direction for a big pit size; or alternatively the odd units of atoms along [1 1 0] direction are removed.     

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.     

Surface electronic structure of α-Mo2C(0 0 0 1)

The electronic structure of α-Mo₂C(0 0 0 1) has been investigated by angle-resolved photoemission spectroscopy utilizing synchrotron radiation. A sharp peak is observed at 3.3 eV in normal-emission spectra. Since the peak shows no dispersion as a function of photon energy and is sensitively attenuated by oxygen adsorption, the initial state of the peak is attributed to a surface state. Resonant photoemission study shows that the state includes substantial contribution of 4d orbitals of the Mo atoms in the second layer. The emissions with constant kinetic energies of 22 and 31 eV above the Fermi level (E_F) are found in normal-emission spectra, and these emissions are interpreted as originating from the Mo N₁N₂₃V and N₂₃VV Auger transitions, respectively.