NO adsorption and diffusion on unreconstructed Pt{1 0 0} surface. A density functional theory investigation

Ab initio density functional theory was used to investigate the adsorption and diffusion of a single NO molecule on the unreconstructed Pt{1 0 0}-(1 × 1) surface. To our knowledge this is the first theoretical study of the NO diffusion activation energy on the Pt{1 0 0} surface. The most stable adsorption position for NO corresponds to the bridge site with the axis of the molecule perpendicular to the surface. The bond of the NO molecule to the surface is through the N-atom. We found that there is a low adsorption energy when the NO molecule is bonded through the O-atom and the axis is perpendicular to the surface, for the three high symmetry sites investigated. NO diffusion between bridge-hollow sites, bridge-atop sites, and hollow-atop sites was also investigated. The barrier for NO diffusion is 0.41 eV, which corresponds to the energy difference between the bridge and hollow sites. This value is around 15% of the highest adsorption energy found on this surface. NO stretch frequencies are also calculated for the three high symmetry sites investigated.     

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.     

Atomic and molecular adsorption on Pt(1 1 1)

The adsorption of several atomic (H, O, N, S, and C) and molecular (N₂, HCN, CO, NO, and NH₃) species and molecular fragments (CN, CNH₂, NH_2, NH, CH₃, CH₂, CH, HNO, NOH, and OH) on the (1 1 1) facet of platinum, an important industrial and fuel cell catalyst, was studied using self-consistent periodic density functional theory (DFT-GGA) calculations at a coverage of 1/4 ML. The best binding site, energy, and position, as well as an estimated diffusion barrier, of each species were determined. The binding strength for all the species can be ordered as follows: N₂ < NH₃ < HCN < NO < CO < CH₃ < OH < NH₂ < H < CN < NH < O < HNO < CH₂ < NOH < CNH₂ < N < S < CH < C. Although the atomic species generally preferred fcc sites, there was no clear trend in site preference by the molecular species or molecular fragments. The vibrational frequencies of all the stable adsorbates in their best and second best adsorption sites were calculated and found to be in good agreement with experimental values reported in the literature. Finally, the decomposition thermochemistry of NOH, HNO, NO, NH₃, N₂, CO, and CH₃ was analyzed.     

Density functional theory investigation of CN on Cu(1 1 1), Ni(1 1 1) and Ni(1 0 0)

The adsorption of CN on Cu(1 1 1), Ni(1 1 1) and Ni(1 0 0) has been investigated using density functional theory (DFT). While experimental studies of CN on Cu(1 1 1) show the molecular axis to be essentially parallel to the surface, the normally-preferred DFT approach using the generalised gradient approximation (GGA) yields a lowest energy configuration with the C-N axis perpendicular to the surface, although calculations using the local density approximation (LDA) do indicate that the experimental geometry is energetically favoured. The same conclusions are found for CN on Ni(1 1 1); on both surfaces bonding through the N atom is always unfavourable, in contrast to some earlier published results of ab initio calculations for Ni(1 1 1)/CN and Ni(1 0 0)/CN. The different predictions of the GGA and LDA approaches may lie in subtly different relative energies of the CN 5σ and 1π orbitals, a situation somewhat similar to that for CO adsorbed on Pt(1 1 1) which has proved challenging for DFT calculations. On Ni(1 0 0) GGA calculations favour a lying-down species in a hollow site in a geometry rather similar to that found experimentally and in GGA calculations for CN on Ni(1 1 0).     

Adsorbate influence on the coercive field of ultrathin Co/Cu{1 1 0}

We report on experimental coercive field measurements for ultrathin films of Co on Cu{1 1 0}, revealing a complex dependence upon the surface coverage of oxygen and hydrogen adatoms. With reference to first-principles calculations, we rationalise this behaviour in terms of initial reaction with surface contaminants and/or adsorption at defects, followed by (in the case of oxygen) subsequent surface ordering and eventual formation of surface oxide.     

Reconstruction and de-reconstruction of the Ir(1 0 0) surface and ultrathin Fe/Ir(1 0 0) films

The structure, energetics and magnetic properties of the quasihexagonal reconstruction of the Ir(1 0 0) surface and nanostructures formed by Fe atoms on this surface have been investigated using first-principles density functional theory with generalized gradient corrections. We find the reconstructed (1 × 5) surface to be 0.10 eV/(1 × 1) area lower in energy than the unreconstructed surface and we demonstrate that first-principles calculations can achieve quantitative agreement with experiment even for such long-period and deep-going reconstructions. For Fe coverage of 0.4 monolayers (ML) we have studied the stripe-like structure with biatomic Fe rows placed in the troughs of the (1 × 5)-reconstructed surface. Results of nonmagnetic calculations agree well with the structure inferred from STM data. Higher Fe coverages lead to a de-reconstruction of the Ir substrate. At 0.8 ML coverage a surface compound with composition Fe₄Ir is formed, which shows an appreciable buckling. In this case, a ferromagnetic calculation leads to good agreement with the low-temperature LEED data. We predict that the (1 × 5) periodicity of the mixed interface layer will persist also in thicker films with a pure Fe surface. Films with 1-4 ML Fe are predicted to be tetragonally distorted and ferromagnetic, with an axial ratio corresponding well to an elastic distortion of the Fe lattice.     

NO dissociation pathways on Rh(1 0 0), (1 1 0), and (1 1 1) surfaces: A comparative density functional theory study

The chemisorption and dissociation pathways of NO on the Rh(1 0 0), (1 1 0), and (1 1 1) surfaces are studied by the plane-wave density functional theory (DFT) with CASTEP program. In addition, the electronic and geometrical effects that affect the NO dissociation reactions have been investigated in detail. The calculation results are presented as following: The effective activation energies of the best NO dissociation pathways on the Rh(1 0 0), the Rh(1 1 0), and the Rh(1 1 1) are 0.63, 0.66 and 1.77 eV, respectively. The activity of the Rh planes for NO dissociation is in the order of Rh(1 0 0) ≈ Rh(1 1 0) > Rh(1 1 1). The low dissociation barrier for Rh(1 0 0) and Rh(1 1 0) is associated with the existence of a lying-down NO structure which acts as a precursor for dissociation. By Mulliken population analysis and structure analysis, both electronic and geometrical effects are found to affect the NO dissociation reactions, but the geometrical effect exceed the electronic. The energy decomposition scheme has been used to provide further insight into the NO dissociation reactions. Based on the calculations, the interaction energy between N and O in the transition state on the Rh(1 1 1) is found much larger than that on the Rh(1 0 0) and the Rh(1 1 0). The major differences of should originate from the variation of the bonding competition effect.     

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.     

Diffusion of Pb adatom and ad-dimer on Si(1 0 0) from ab initio studies

The diffusion pathways of Pb adatoms and ad-dimers on Si(1 0 0) are investigated by first-principles calculations. Pb adatoms are found to diffuse on top of the Si(1 0 0) dimer row with an energy barrier of 0.31 eV. However, Pb dimers are energetically more stable. Pb dimers on top of the dimer row have a high energy barrier (0.95 eV) to rotate from the lowest energy configuration to the orientation parallel to the underlying Si(1 0 0) dimer row. Once the ad-dimer is oriented parallel to Si(1 0 0) dimer row, they can diffuse along the dimer row with an energy barrier of only 0.32 eV.     

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.     

A surface X-ray diffraction study of TiO2(1 1 0)(3 × 1)-S

Surface X-ray diffraction has been used to investigate the structure of TiO₂(1 1 0)(3 × 1)-S. In concert with existing STM and photoemission data it is shown that on formation of a (3 × 1)-S overlayer, sulphur adsorbs in a position bridging 6-fold titanium atoms, and all bridging oxygens are lost. Sulphur adsorption gives rise to significant restructuring of the substrate, detected as deep as the fourth layer of the selvedge. The replacement of a bridging oxygen atom with sulphur gives rise to a significant motion of 6-fold co-ordinated titanium atoms away from the adsorbate, along with a concomitant rumpling of the second substrate layer.     

Surface reconstruction of clean bcc-Fe{1 1 0}: A quasi-hexagonal top-layer with periodic height modulation

Hexagonal-pillar shaped pure Fe single crystal whiskers with six {1 1 0} side planes were obtained by means of chemical vapor deposition. Atomically resolved scanning tunneling microscopy images obtained on the {1 1 0} surface showed a quasi-hexagonal atomic array with mesoscopic-range periodic height modulation of about 1/3 of an atomic step. This height modulation was found to be a result of an interference between the quasi-hexagonal top-layer and the sub-surface bcc-Fe{1 1 0} layer. Unit vectors of the mesoscopic-range modulation turned out to be expressed as , where and are the primitive vectors of the two-dimensional atomic array in the top-layer and those in the sub-surface layer, respectively. The two-dimensional density of atoms in the top-layer is slightly higher by 0.46% than that in the sub-surface layer.     

Driving force for the WO3(0 0 1) surface relaxation

The optimized structure of the WO₃(0 0 1) surface with various types of termination ((1 × 1)O, (1 × 1)WO₂, and c(2 × 2)O) has been simulated using density functional theory with the Perdew-Wang 91 gradient corrected exchange-correlation functional. While the energy of bulk WO₃ depends weakly on the distortions and tilting of the WO₆ octahedra, relaxation of the (0 0 1) surface results in a significant decrease of surface energy (from 10.2 × 10⁻² eV/Ų for the cubic ReO₃-like, c(2 × 2)O-terminated surface to 2.2 × 10⁻² eV/Ų for the relaxed surface). This feature illustrates a potential role of surface relaxation in formation of crystalline nano-size clusters of WO₃. The surface relaxation is accompanied by a dramatic redistribution of the density of states near the Fermi level, in particular a transformation of surface electronic states. This redistribution is responsible for the decrease of electronic energy and therefore is suggested to be the driving force for surface relaxation of the WO₃(0 0 1) surface and, presumably, similar surfaces of other transition metal oxides.     

Anomalous hybridization in the In-rich InAs(0 0 1) reconstruction

The surface bonding arrangement in nearly all the confirmed reconstructions of InAs(0 0 1) and GaAs(0 0 1) have only two types of hybridization present. Either the bonds are similar to those in the bulk and the surface atoms are sp³ hybridized or the surface atoms are in a tricoordinated bonding arrangement and are sp² hybridized. However, dicoordinated In atoms with sp hybridization are observed on the InAs(0 0 1), In-rich, room temperature and low temperature surfaces. Scanning tunneling microscopy (STM) images of the room temperature (300 K) InAs(0 0 1) surface reveal that the In-rich surface reconstruction consists of single-atom rows with areas of high electron density that are separated by ∼4.3 Å. The separation in electron density is consistent with rows of undimerized, sp hybridized, In atoms, denoted as the β3′(4 × 2) reconstruction. As the sample is cooled to 77 K, the reconstruction spontaneously changes. STM images of the low temperature surface reveal that the areas of high electron density are no longer separated by ∼4.3 Å but instead by ∼17 Å. In addition, the LEED pattern changes from a (4 × 2) pattern to a (4 × 4) pattern at 77 K. The 77 K reconstruction is consistent with two (4 × 2) subunit cells; one that contains In dimers on the row and another subunit cell that contains undimerized, sp hybridized, In atoms on the row. This combination of dimerized and undimerized subunit cells results in a new unit cell with (4 × 4) periodicity, denoted as the β3(4 × 4) reconstruction. Density functional theory (DFT) and STM simulations were used to confirm the experimental findings.     

First-principles study of oxidized Nb(1 0 0) surface structures

The atomic positions of the oxygen-induced c(2 × 2)-O, (3 × 1)-O and (4 × 1)-O surface structures on Nb(1 0 0) are determined by first-principles electronic structure calculations within the density functional theory comparing experimentally observed scanning tunneling microscopy (STM) images. STM images of these surfaces are calculated on the basis of the theory of Tersoff and Hamann. The theoretical and experimental STM images of the oxygen-chemisorbed c(2 × 2)-O structural model agree well. However, only the oxide-covered (3 × 1)-O and (4 × 1)-O structural models with two layers of NbO and contraction of the unit length along longitudinal 〈1 0 0〉 direction by 10% result in the theoretical STM images that agree with the experimental ones.     

Electronic structures of the suboxide films formed on TiC(1 0 0) and ZrC(1 0 0) surfaces: Density functional theory studies

Density functional theory (DFT) calculations have been performed to elucidate the electronic structures of the TiO-like film on TiC(1 0 0) and the ZrO-like film on ZrC(1 0 0), which are assumed to be monolayers of suboxide films with (1 × 1) periodicity with respect to the substrate (1 0 0) surfaces. It was revealed that the electronic structures of both films were characterized by the existence of a band around 6 eV and a band around the Fermi level. The former and latter bands were mostly composed of O 2p and metal d orbitals, respectively, indicating the substantial ionic nature of the film. The calculated DOS well reproduced the previously obtained photoelectron spectra. From the inspection of the optimized structures, it was found that the both suboxide films have rippled structures; the metal and oxygen atoms are displaced vertically downward and upward, respectively, maintaining the (1 × 1) structures.     

Photoelectron spectroscopy study of oxygen adsorption on Mo2C(0 0 0 1)

Oxygen adsorption on the α-Mo₂C(0 0 0 1) surface has been investigated with X-ray photoelectron spectroscopy and valence photoelectron spectroscopy utilizing synchrotron radiation. It is found that oxygen adsorbs dissociatively at room temperature, and the adsorbed oxygen atoms interact with both Mo and C atoms to form an oxycarbide layer. As the O-adsorbed surface is heated at ≧800 K, the C-O bonds are broken and the adsorbed oxygen atoms are bound only to Mo atoms. Valence PES study shows that the oxygen adsorption induces a peculiar state around the Fermi level, which enhances the emission intensity at the Fermi edge in PES spectra.     

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.     

Determining the GaSb/GaAs-(2 × 8) reconstruction

Highly strained thin layers of GaSb/GaAs possess a (2 × 4) reconstruction at low Sb overpressures, and a (2 × 8) reconstruction at high Sb overpressures. While the atomic details of the Sb/GaAs-(2 × 4) are well known, the details of the (2 × 8) are not understood. In this paper, we use density functional theory to analyze possible (2 × 8) structures. Comparing scanning tunneling microscope images from both simulation and experiment and examining the relative energies of possible (2 × 8) structures, we show the α(2 × 8) and β(2 × 8) are the thermodynamically stable surface reconstructions for high Sb content films strained to the GaAs lattice parameter. The α and β(2 × 8) reconstructions are related to the GaAs-α2(2 × 4) and GaAs-β2(2 × 4) through the addition of 2 cations and 8 anions into the trench between adjacent (2 × 4) unit cells.