First principle study of SrTiO3 (0 0 1) surface and adsorption of NO on SrTiO3 (0 0 1)

The results of first-principles calculations about the two possible terminations of (0 0 1) surface of SrTiO₃ perovskite and adsorption of NO on SrTiO₃ (0 0 1) surface were presented. Both surface parameters (atomic structures and electronic configurations) and adsorption parameters (bond, energy and charge) of NO on SrTiO₃ (0 0 1) surface, which have never been investigated before as far as we know were investigated using density functional theory calculations with the local-density approximation (DFT-LDA). It was found that the two possible terminations of SrTiO₃ (0 0 1) surface have large surface relaxation, which leads to surface polarization and exhibits different reactivity toward the dissociative adsorption of NO. The electronic states of TiO₂-terminated surface have larger difference than that of bulk, so it is more favorable for adsorption of NO than SrO-terminated surface.     

First-principles study of the (0 0 1) surface of cubic SrZrO3

Density functional calculations have been used to investigate the (0 0 1) surface of cubic SrZrO₃ with both SrO and ZrO₂ termination. Surface structure and electronic structure have been obtained. The SrO surface is found to be similar to its counterpart in SrTiO₃, while there are marked differences between the ZrO₂ and TiO₂ terminations in SrZrO₃ and SrTiO₃, respectively, concerning surface relaxation and rumpling. For the ZrO₂-terminated surface of SrZrO₃, the covalency of the interaction between the outmost Zr and the O beneath is enhanced as a result of their bond contraction. The band gap reduction and the presence of the surface states are also discussed in relation with the behavior of the electrostatic potential.     

DFT study of Pt adsorption on low index SrTiO3 surfaces: SrTiO3(1 0 0), SrTiO3(1 1 1) and SrTiO3(1 1 0)

Density Functional Theory has been used to determine the energetically preferred structures of submonolayer, monolayer, and multilayer Pt films on both ideal terminations of SrTiO₃(1 0 0), SrTiO₃(1 1 1), and SrTiO₃(1 1 0). The strength of the resulting metal/metal oxide interfaces was characterized by the adsorption energy of the film and the film’s work of separation. The two polar surfaces, SrTiO₃(1 1 1) and SrTiO₃(1 1 0), form significantly stronger interfaces than the non-polar SrTiO₃(1 0 0) surface. Approximate criteria were applied to predict the growth mode of Pt on each surface.     

Surface properties of rutile TiO2(1 1 0) from molecular dynamics and lattice dynamics at 300 K: Variable-charge model results

Geometric structure, atomic vibrations and atomic charges and their thermally induced fluctuations have been calculated as a function of depth in, and thickness of, rutile TiO₂(1 1 0) slabs, within the framework of the variable-charge potential of Swamy and Gale [V. Swamy, J.D. Gale, Phys. Rev. B 62 (2000) 5406] at 300 K. Molecular dynamics simulations and lattice dynamics calculations were performed with a 2D periodic slab model for slab thicknesses between 3 and 11 triple layers (approximately 9-35 Å). Odd-even oscillations with respect to the number of slab layers are found for the surface relaxation for very thin slabs, and for the (slowly converging) rumpling in the middle of the slab. The Ti and O atomic charges in the outermost three atomic layers differ from the rest of the slab (they are less ionic); the thermal vibrations do not alter this picture. The atomic mean-square amplitudes are some 50% larger (more for O, less for Ti) at the surface than in the middle of the slab and decay rather slowly to the bulk values. Comparisons with the results of a rigid-ion potential for titania [M. Matsui, M. Akaogi, Mol. Simul. 6 (1991) 238] are presented for non-electronic properties.     

Density functional theory investigation of the structure of SO2 and SO3 on Cu(1 1 1) and Ni(1 1 1)

Density functional theory (DFT) slab calculations, mainly using the generalised gradient approximation, have been used to investigate the minimum energy structures of molecular SO₂ and SO₃ on Cu(1 1 1) and Ni(1 1 1) surfaces. On Ni(1 1 1) the optimal local adsorption structures are in close agreement with experimental results for both molecular species obtained using the X-ray standing wavefield technique, although for adsorbed SO₂ the energetic difference between two alternative lateral positions of the lying-down molecule on the surface is marginally significant. On Cu(1 1 1) the results for adsorbed SO₂, in particular, were sensitive to the DFT functional used in the calculations, but in all cases failed to reproduce the experimentally-established preference for adsorption with the molecular plane perpendicular to the surface. This result is discussed in the context of previously published DFT results for these species adsorbed on Cu(1 0 0). The optimal geometry found for SO₃ on Cu(1 1 1) is similar to that on Ni(1 1 1), providing agreement with experiment regarding the molecular orientation but not the adsorption site.     

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.     

Structure and electronic properties calculation of ultrathin α-Al2O3 films on (0 0 0 1) α-Cr2O3 templates

Ab initio total energy Hartree-Fock calculations of ultrathin films of α-Al₂O₃ on (0 0 0 1) α-Cr₂O₃ templates are presented. The surface relaxation, the in-plane reconstruction and the surface and strain energies of the slabs are studied as a function of alumina film thickness. The surface Al layer is found to relax inwards considerably, with the magnitude of the inwards relaxation depending on the thickness of the ultrathin alumina film in a non-linear manner. The calculations also reveal that ultrathin films of alumina lower the surface energy of (0 0 0 1) α-chromia substrates. This indicates that the (0 0 0 1) α-chromia surface provides favourable conditions for the templated growth of α-alumina. However, increasing the alumina film thickness is found to give rise to a significant increase in strain energy. Finally, the electronic properties at the surface of the (0 0 0 1) α-Al₂O₃/α-Cr₂O₃ slabs are investigated. Here it is found that the alumina coating gives rise to an increase in the covalency of the bonds at the surface of the slabs. In contrast, the influence of an alumina layer on the electrostatic potential at the surface of the chromia slab is relatively minor, which should also be beneficial for the templated growth of α-alumina on (0 0 0 1) α-chromia substrates.     

Structural relaxation and Jahn-Teller distortion of LaMnO3 (0 0 1) surface

We studied in detail the structural relaxation and Jahn-Teller distortion in LaMnO₃ (0 0 1) surface of the orthorhombic phase by means of classical atomistic simulation. It is found that MnO₂-terminated surface is more energetically favorable than LaO-terminated surface by 0.34 eV. The standard deviation of Mn-O bond lengths of MnO₆ octahedra and Jahn-Teller distortion oscillate in LaMnO₃ (0 0 1) surface. Our simulated atomic displacements in the surface are compared with some ab initio studies.     

Molecular dynamics simulation of the (0 0 0 1) α-Al2O3 and α-Cr2O3 surfaces

A simple, rigid pair-potential model is applied to investigate the dynamics of the (0 0 0 1) α-Al₂O₃ and α-Cr₂O₃ surfaces using the molecular dynamics technique. The simulations employ a two-stage equilibration process: in the first stage the simulation-cell size is determined via the constant-stress ensemble, and in the second stage the equilibration of the size-corrected simulation cell is continued in the canonical ensemble. The thermal expansion coefficients of bulk alumina and chromia are evaluated as a function of temperature. Furthermore, the surface relaxation and mean-square displacement of the atoms versus depth into the slab are calculated, and their behaviour in the surface region analysed in detail. The calculations show that even moderate temperatures (∼400 °C) give rise to displacements of the atoms at the surface which are similar to the lattice mismatch between α-alumina and chromia. This will help in the initial nucleation stage during thin film growth, and thus facilitate the deposition of α-Al₂O₃ on (0 0 0 1) α-Cr₂O₃ templates.     

Ab initio simulation of the BaZrO3 (0 0 1) surface structure

Electronic and atomic structures of different terminations of the (0 0 1) non-polar orientation of BaZrO₃ surfaces have been studied using first-principles calculations. We found that surface energies at both possible surface terminations, BaO and ZrO₂, were very close. The (0 0 1)-BaO and (0 0 1)-ZrO₂ terminated surfaces have bandgap values smaller than that of a bulk BaZrO₃ crystal. In addition, the relative surface stability has been analyzed as a function of chemical environment.     

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.     

Effects of N adsorption on the structural and electronic properties of SrTiO3(0 0 1) surface

The atomic configurations, bonding characteristics, and electronic structures of the N-adsorbed (directly and substitutionally) SrTiO₃(0 0 1) surface are studied by using first-principles method based on the density functional theory. From the analysis of the energetics and density of states, it is found that the stability of the directly adsorbed N depends on the relative position of N atom to the surface. To better understand the effects of the substitutionally adsorbed N on the surface, as an example, the behavior of Au atoms adsorbed on the N-substituted surface is discussed in detail. There is clearly a synergy effect between the substitution of N for O_s(I) and the adsorption of Au atoms on the SrTiO₃(0 0 1) surface.     

Adsorption of water on thin V2O3(0 0 0 1) films

V₂O₃(0 0 0 1) films have been grown epitaxially on Au(1 1 1) and W(1 1 0). Under typical UHV conditions these films are terminated by a layer of vanadyl groups as has been shown previously [A.-C. Dupuis, M. Abu Haija, B. Richter, H. Kuhlenbeck, H.-J. Freund, V₂O₃(0 0 0 1) on Au(1 1 1) and W(1 1 0): growth, termination and electronic structure, Surf. Sci. 539 (2003) 99]. Electron irradiation may remove the oxygen atoms of this layer. H₂O adsorption on the vanadyl terminated surface and on the reduced surface has been studied with thermal desorption spectroscopy (TDS), vibrational spectroscopy (IRAS) and electron spectroscopy (XPS) using light from the BESSY II electron storage ring in Berlin. It is shown that water molecules interact only weakly with the vanadyl terminated surface: water is adsorbed molecularly and desorbs below room temperature. On the reduced surface water partially dissociates and forms a layer of hydroxyl groups which may be detected on the surface up to T ∼ 600 K. Below ∼330 K also co-adsorbed molecular water is detected. The water dissociation products desorb as molecular water which means that they recombine before desorption. No sign of surface re-oxidation could be detected after desorption, indicating that the dissociation products desorb completely.     

Structural and electronic properties of cubic SrHfO3 surface: First-principles calculations

Structural, electronic and chemical bonding properties of the (0 0 1) surface of cubic SrHfO₃ have been investigated with both SrO and HfO₂ termination using the plane-wave ultrasoft pseudopotential technique based on the first-principles density-functional theory. The relaxed structures of two slabs have been analyzed, which shows the interplanar distance of two slabs has the same changed trend. The electronic band structures and density of states of two slabs have been discussed, showing the reduced band gaps by comparison with those of bulk system. The chemical bonding between Sr and O between the surface layer and subsurface layer as well as Hf and O has been increased. The surface energy, work function and stability have been calculated, which indicates SrO-terminated slab is more stable.     

Hydroxylated α-Al2O3 (0 0 0 1) surfaces and metal/α-Al2O3 (0 0 0 1) interfaces

X-ray photoelectron spectroscopy was applied to study the hydroxylation of α-Al₂O₃ (0 0 0 1) surfaces and the stability of surface OH groups. The evolution of interfacial chemistry of the α-Al₂O₃ (0 0 0 1) surfaces and metal/α-Al₂O₃ (0 0 0 1) interfaces are well illustrated via modifications of the surface O1s spectra. Clean hydroxylated surfaces are obtained through water- and oxygen plasma treatment at room temperature. The surface OH groups of the hydroxylated surface are very sensitive to electron beam illumination, Ar⁺ sputtering, UHV heating, and adsorption of reactive metals. The transformation of a hydroxylated surface to an Al-terminated surface occurs by high temperature annealing or Al deposition.     

Theoretical study of the electronic structure of the Si3N4(0 0 0 1) surface

Density functional theory has been applied to a study of the electronic structure of the ideally-terminated, relaxed and H-saturated (0 0 0 1) surfaces of β-Si₃N₄ and to that of the bulk material. For the bulk, the lattice constants and atom positions and the valence band density of states are all in good agreement with experimental results. A band gap of 6.7 eV is found which is in fair accord with the experimental value of 5.1-5.3 eV for H-free Si₃N₄. Using a two-dimensionally-periodic slab model, a π-bonding interaction is found between threefold-coordinated Si and twofold-coordinated N atoms in the surface plane leading to π and π* surface-state bands in the gap. A surface-state band derived from s-orbitals is also found in the gap between the upper and lower parts of the valence band. Relaxation results in displacements of surface and first-underlayer atoms and to a stronger π-bonding interaction which increases the π-π* gap. The relaxed surface shows no occupied surface states above the valence band maximum, in agreement with recent photoemission data for a thin Si₃N₄ film. The π* band, however, remains well below the conduction band minimum (but well above the Fermi level). Adsorbing H at all dangling-bond sites on the ideally-terminated surface and then relaxing the surface and first underlayer leads to smaller, but still finite, displacements in comparison to the clean relaxed surface. This surface is more stable, by about 3.67 eV per H, than the clean relaxed surface.     

Adsorption of CO, NH3 and O2 on Fe site of La0.875Sr0.125FeO3 (0 1 0) surface by density function theory calculation

The adsorption of CO, NH₃ and O₂ gas molecules on Fe site of La_0.875Sr_0.125FeO₃ (0 1 0) surface has been investigated based on the density functional theory (DFT) with the spin polarized generalized gradient approximation (GGA). The optimal adsorption orientations as well as the adsorption energies for these molecules adsorption on Fe site of La_0.875Sr_0.125FeO₃ (0 1 0) surface are derived. Results show that adsorption configurations with C-down for CO, N-down for NH₃ are stable. For the O₂ molecule adsorption, the mode with an angle about 120° between the oxygen molecule and La_0.875Sr_0.125FeO₃ (0 1 0) surface is stable.     

Adsorption and reactions of HN3 on Si(1 0 0)-2 × 1: A computational study

We have studied the adsorption and decomposition of HN₃ on Si(1 0 0)-2 × 1 surface using the hybrid density functional B3LYP method and effective core potential basis, LanL2DZ, with Si₁₅H₁₆ as a double dimer surface model for cluster calculations. The result shows that the barriers for the dissociative adsorption of HN₃ forming HN(a) + N₂(g) are quite low by stepwise dissociation processes occurring either on a dimer or across the dimers. The low activation energies are consistent with previous experimental observations that the molecularly adsorbed HN₃(a) can undergo decomposition producing HN(a) at low surface temperatures. On the other hand, the predicted activation energies for the N₃(a) + H(a) formation processes are all relatively higher. These results also explain the absence of the N₃(a) species in HREELS measurements following each annealing experiment. Several selected reaction paths were also confirmed with slab model calculations using an optimization approach coupling the energy and gradient calculations by the slab model with the geometrical optimization using Berny algorithm.In addition, the adsorbate effect was examined for the end-on and side-on molecular configurations. For the side-on adsorption configuration, all possible combinations with 1-4 adsorbates can exist on the four surface Si sites of the double dimers, with adsorption energies lying closely to the multiples of that of a single side-on adsorbate (LM2); i.e., adsorption energies are nearly additive. Interestingly, for the end-on adsorption, only 1 and 2 HN₃ molecules can adsorb on a dimer due to the presence of the negative charges on the remaining Si sites in the neighboring dimer. For the two end-on adsorbates on the same dimer, the total adsorption energy is close to two times that of HN₃(a) or LM1. For the mixed end-on/side-on configurations, only one of each type can co-exist on a single dimer pair (Si₁-Si₂ or Si₃-Si₄) sites with adsorption energy close to the sum of those of one end-on and one side-on adsorbate. Finally, the predicted reaction routes and vibrational frequencies showed good agreement with previous experimental results. The stabilities of many ad-species involved in these reactions with end-on and/or side-on configurations have been predicted together with the transition states connecting those species.     

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.     

First-principles study of CaTiO3 oxygen-vacancies (0 0 1) surface

The structural and electronic properties of fully relaxed CaTiO₃ oxygen-vacancies (0 0 1) surface with CaO and TiO₂ terminations are investigated by first-principles plane waves ultrasoft pseudopotential method based on local density approximation. The present results show that the direction and the magnitudes of the atomic relaxations for oxygen-vacancies surface are different from that of the perfect surface. Compared with the TiO₂-terminated oxygen-vacancy surface, the CaO-terminated oxygen-vacancy surface are likely to be observed in oxygen environment conditions. Much different from the perfect surface, the oxygen-vacancy surface becomes metallic caused by some states in the conduction band are lowered and pulled down in the band-gap region.