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.     

First principle calculations of benzotriazole adsorption onto clean Cu(1 1 1)

The adsorption of benzotriazole (BTAH or C₆N₃H₅) on a Cu(1 1 1) surface is investigated by using first principle density functional theory calculations (VASP). It is found that BTAH can be physisorbed (<0.1 eV) or weakly chemisorbed (∼0.43 eV) onto Cu(1 1 1), and the chemical bond is formed through nitrogen sp² lone pairs. The weak chemisorption can be stabilized by reaction with neighboring protonphilic radicals, like OH⁻. Furthermore, the geometries and associated energies of intermolecular hydrogen bonds between adsorbates on Cu(1 1 1) are also calculated. A model of the first layer of BTAH/BTA⁻ on Cu(1 1 1) surface is developed based on a hydrogen bond network structure.     

Cu(1 1 1) and Cu(0 0 1) surface electronic states. Comparison between theory and experiment

Recent advances in both the experimental resolution and in the computational capabilities motivate new studies of surface properties. In particular, a detailed comparison between theoretical and experimental data is expected to provide a better insight into surface and image states. In this work we present a joint effort analyzing such features of the Cu(1 1 1) and Cu(0 0 1) surfaces. The experiments are performed by using both linear and non-linear angle-resolved photoemission. From the theoretical point of view, we make use of the Green function embedding technique within density functional theory. We are able to account for the image states by suitably modifying the effective potential in the Kohn-Sham equation and the generalized boundary conditions on the Green function. Comprehensive theoretical and experimental results on the effective mass and the binding energy of the Shockley state and the first image states are reported.     

Quantitative determination of the local structure of thymine on Cu(1 1 0) using scanned-energy mode photoelectron diffraction

The local adsorption structures of the surface species formed by interaction of thymine with a Cu(1 1 0) surface at room temperature, and after heating to ∼530 K, have been investigated. Initial characterisation by soft-X-ray photoelectron spectroscopy and O K-edge near-edge X-ray absorption fine structure (NEXAFS) indicates the effect of sequential dehydrogenation of the NH species and provides information on the molecular orientation. O 1s and N 1s scanned-energy mode photoelectron diffraction shows the species at both temperatures bond to the surface through both carbonyl O atoms and the deprotonated N atom between them, each bonding atom adopting near-atop sites on the outermost Cu surface layer. The associated bondlengths are 1.96 ± 0.03 Å for Cu-N and 1.91 ± 0.03 Å and 2.03 ± 0.03 Å for the two inequivalent Cu--O bonds. The molecular plane lies almost exactly in the close-packed azimuth, but with a tilt relative to the surface normal of approximately 20°. Heating to ∼530 K, or deposition at this temperature, appears to lead to dehydrogenation of the second N atom in the ring, but no significant change in the adsorption geometry.     

Quantum states of hydrogen (H, D, T) atoms on Cu(1 0 0) and (1 1 0) surfaces

We investigate the quantum mechanical behavior of adsorbed hydrogen (H, D, T) on Cu(1 0 0) and (1 1 0) surfaces. We construct potential energy surfaces (PESs) for the motion of the hydrogen H atom on Cu(1 0 0) and (1 1 0) surfaces within the framework of density functional theory. The potential energy takes a minimum value on the hollow site of Cu(1 0 0) and on the short bridge site of Cu(1 1 0). Moreover, we calculate the quantum states of hydrogen atom motion on these calculated PESs. The ground state wave function of the hydrogen atom motion is strongly localized around the hollow site on the Cu(1 0 0) surface. On the other hand, the ground state wave function of the hydrogen atom motion on Cu(1 1 0) is distributed from the short bridge site to two neighboring pseudo-threefold sites. We finally show isotope effects on the quantum states of the motion of hydrogen on both surfaces.     

Adsorption of chlorophenol on the Cu(1 1 1) surface: A first-principles density functional theory study

The interaction between a 2-chlorophenol (C₆H₄OHCl) molecule and the Cu(1 1 1) surface has been investigated using density functional theory as an initial step in gaining a better understanding of the catalyzed formation of dioxin compounds on a clean copper surface. The 2-chlorophenol molecule is found to form several weakly bonded, horizontally and vertically oriented configurations. Dissociative modes have also been investigated. For the latter, the formation of phenyl and benzyne fragments is found to be more energetically favourable than the formation of 2-chlorophenoxy radicals.     

Dissociative chemisorption of H2 on Pt(1 1 1): isotope effect and effects of the rotational distribution and energy dispersion

Six-dimensional quantum dynamics calculations on dissociative scattering of H₂ and D₂ from Pt(1 1 1) are performed. The six-dimensional potential energy surface used was generated using density functional theory employing the generalized gradient approximation. The isotope effect, the effect of widening the rotational distribution on the dissociation probability and the effect of the energy dispersion are investigated, as they are possible reasons for a discrepancy between previous theoretical work and molecular beam experiments. It was found that these effects cannot explain the differences between the theoretical and experimental results.     

Optical anisotropy at the W(1 1 0) surface

We report ab initio calculations of the anisotropic dielectric function of tungsten (1 1 0) surface using the linear muffin-tin-orbital method. The calculated anisotropy in the optical spectrum, for polarization of light parallel to the surface, exhibits three dominant broad structures at 3.00, 4.01 and 5.34 eV successively positive, negative and then positive. The first peak is clearly assigned to p → d interband transitions in surface atomic sites whereas the two others have their origin in interband transitions in bulk like atoms. Our results, including the interlayer relaxation effect on the surface optical response, are compared to recent reflectance anisotropy measurements.     

Empty-state scanning tunneling microscopy image of C2H2 on Si(0 0 1)--new evidence for paired end-bridge di-σ configuration

New evidence of the paired end-bridge configuration in the room-temperature adsorption geometries of C₂H₂ molecules on Si(0 0 1) is presented by scanning tunneling microscopy (STM) and ab initio pseudopotential calculations. The distinct four-leaved feature occupying two adjacent Si dimer sites in the experimental empty-state STM images are well reproduced by simulations of the paired end-bridge adsorption configuration. Calculated energetics suggests that the Si(0 0 1) surface is covered by paired end-bridge structures at the saturation coverage of 1 ML, in agreement with the existing experiments.     

Reaction of chlorinated benzenes with Si(1 0 0)2 × 1: a theoretical study

Theoretical (HF + DFT) investigations of the adsorption of chlorobenzene (ClPh), 1,2- and 1,4-dichlorobenzene (1,2-diClPh and 1,4-diClPh) on a silicon (1 0 0) surface are reported for the first time, and are compared with one another and with benzene. Binding energies for various structures with the molecules attached on-top and in-between the surface dimer rows are correlated with the STM experimental data. Novel structures with the molecules linking two dimer rows, stabilised by detachment of Cl (or H)-atoms forming Cl-Si (or H-Si) bonds, are described. For 1,4 and 1,2 binding, these linking structures are predicted to attach the phenyl ring parallel or perpendicular to the Si surface, respectively, while preserving its aromaticity. The potential-energy barriers between several different structures are evaluated, and compared with available experimental evidence. For 1,4-diClPh it is shown that the potential-energy barrier for the second Cl transfer is significantly lower than for the first one in contrast to the gas-phase, and comparable to the barrier for lifting the Bz-ring into a vertical position and forming a singly bonded ‘displaced’ structure. The predicted barrier-heights are consistent with the experimentally observed relative occurrence of the on-top, linking, and displaced structures.     

Local structure determination of a chiral adsorbate: Alanine on Cu(1 1 0)

N 1s and O 1s scanned-energy mode photoelectron diffraction (PhD) has been used to investigate the local structure of a single enantiomer of deprotonated alanine, alaninate, NH₂CH₃CHCOO-, on Cu(1 1 0) in the (3 × 2) phase. The local site is found to be similar to that of glycinate on Cu(1 1 0), with the N atoms in near-atop sites and the O atoms sites consistent with bonding to single surface Cu atoms but substantially off-atop. Unlike the Cu(1 1 0)(3 × 2)pg-glycinate phase, however, in which the two molecular species per unit mesh are mirror images of one another in identical local sites, the intrinsic chirality of l-alaninate means that the two molecules per unit mesh of the (3 × 2) surface phase occupy slightly different local sites. However, an excellent fit to the PhD data can be achieved by a minor modification of the structure found in DFT calculations by R.B. Rankin and D.S. Sholl [Surf. Sci. 574 (2005) L1] in which the heights of the N and O atoms above the surface are reduced by approximately 0.1 Å. The resulting average N-Cu and O-Cu values are 2.02 and 1.98 Å, respectively, with an estimated precision of ±0.03 Å. These bondlengths are shorter than those obtained from DFT by 0.08 and 0.10 Å, respectively.     

LCAO study of the Cu(1 1 0)-p(2 × 1)O surface

We have performed semi-empirical LCAO calculations of the electronic structure of the Cu(1 1 0)-p(2 × 1)O surface. This has been done accounting for the Cu-Cu interactions by means of a recently proposed set of parameters, which give very good results for the bulk as well as for the surfaces of lowest Miller indices. Furthermore, the O-O interactions, which have been neglected in the preceding similar studies, have been taken into account. The resulting surface bands are in very good agreement with the overall set of the available experimental data. Several issues concerning the physical properties of this surface are addressed in the present paper: the changes induced on the clean surface bands by the adsorption and the reconstruction; the arrangement of the Cu and O atoms in the added rows; the position of the p_y antibonding band of the oxygen. In particular, we have found that the latter has an energy of −0.2 eV at the point. This result confirms an experimental indication in the same direction previously reported by Courths et al. [R. Courths, S. Hüfner, P. Kemkes, G. Wiesen, Surf. Sci. 376 (1997) 43].     

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.     

Dependence of surface properties on adsorbate-substrate distance: Work function changes and binding energy shifts for I/Pt(1 1 1)

This paper presents an analysis of the character of the bond of I adsorbed at on-top and 3-fold sites of Pt(1 1 1). At both sites, the bonding is dominated by an ionic interaction supplemented with some covalent character due to donation from the adsorbed I anion to the Pt surface. The way in which the I-Pt interaction affects observed properties has been established. In particular, the origins of the anomalous work function changes induced by the adsorption of I and the shifts of I core level binding energies are explained. It is shown that the magnitudes of the changes in these properties can be directly correlated with the distance of the I from the Pt surface. Thus, these shifts can be interpreted to indicate adsorbate height. The fact that the negatively charged I adsorbate leads to a work function decrease, rather than the increase expected due to the charge of the adsorbate, may appear to be an anomaly. However, it is shown that this decrease arises from electronic reorganizations that cancel the dipole due to the charge of the adsorbate. Furthermore, the electronic terms that contribute to a lowering of the work function are larger as the adsorbate moves closer to the surface.     

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.     

Cu adsorption on pyrite (1 0 0): Ab initio and spectroscopic studies

Surface optimised S 2p photoelectron spectra show that both surface S²⁻ monomers and (S-S)²⁻ dimers are present at pyrite (1 0 0) fracture surfaces. In order to determine which sulfur species are involved in Cu adsorption, fresh pyrite surfaces were exposed to Cu²⁺ in solution. The S 2p spectra suggest that both types of S surface species are involved in the mechanism of Cu adsorption (activation). Ab initio density functional theory was used to model Cu adsorbed onto pyrite (1 0 0) to support the interpretation of the spectroscopy. Mulliken population analysis confirms the charge distribution suggested by the core line shifts as observed in the photoelectron spectra. The ab initio calculations were consistent with a two-coordinate bond between Cu(I), a surface S monomer and a surface S 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.     

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.