Structure and magnetism of bcc(1 1 0)-like Fe films grown on Cu(0 0 1) via ion beam triangulation and spin-polarized electron capture

The structure and magnetism of thin epitaxial Fe layers grown on Cu(0 0 1) is investigated by grazing scattering of fast H and He atoms. Information on the atomic structure of the film and substrate surfaces is obtained by making use of ion beam triangulation with protons. The magnetic behavior is studied via the polarization of light emitted after capture of spin-polarized electrons into excited atomic terms during scattering of He atoms. For the formation of bcc(1 1 0)-like Fe films at higher coverages, we detect differences in structural and magnetic properties for room and low temperature growth. We suggest that the crystalline structure depends on the film morphology and that Cu impurities affect the magnetic properties.     

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 study of the CO adsorption on ferromagnetic Co(0 0 0 1) and Co(1 1 1) surfaces

The regular CO overlayers at coverage θ = 1/3 adsorbed on the (0 0 0 1) surface of hcp Co and (1 1 1) surface of fcc Co are studied by first-principles density-functional theory with the exchange-correlation component in the PBE form. Adsorption in atop, bridge, and three-fold hcp or fcc position are considered. The adsorption energies, CO stretching frequencies, geometry, work function, and local magnetic moments are studied, and, when possible, compared with experimental or theoretical data. Particularly, we show that the recently proposed correction to adsorption energy of CO prefers correctly the atop adsorption site, whereas the remaining sites are almost degenerate in energy. The CO molecule lowers magnetization on neighbouring Co atoms, and the effect decreases with the adsorption site coordination. We show, however, that this trend is not the result of the different C-Co separation at different adsorption sites. A very small magnetic moment appears on CO that couples antiferromagnetically to Co. Most results are very similar for the Co(0 0 0 1) and Co(1 1 1) surfaces.     

Study of Si(0 0 1)4 × 2-Ga structure by scanning tunneling microscopy and ab initio calculation

We have studied Si(0 0 1)-Ga surface structures formed at Ga coverages of slightly above 0.50 monolayer (ML) at 250 °C by scanning tunneling microscopy (STM). 4 × 2-, 5 × 2-, and 6 × 2-Ga structures were observed in a local area on the surface. The 4 × 2-Ga structure consists of three protrusions, as observed in filled- and empty-state STM images. The characters of these structures are clearly different from those of other Si(0 0 1)-Ga structures. We also performed an ab initio calculation of the energetics for several possible models for the 4 × 2-Ga structure, and clarified that the three-orthogonal-Ga-dimer model is the most stable. Also, the results of comparing the simulated STM images and observation images at various bias voltages indicate that this structural model is the most favorable.     

A multi-scale Monte Carlo study of oxide structures on the Cu(1 0 0) surface

We present a multi-scale Monte Carlo study of the oxidation of the Cu(1 0 0) surface based on the Bortz-Kalos-Lebowitz model with the equilibrium energetics obtained from ab initio calculations. The radial and island size distribution functions are examined and Cu-O structures are analyzed at different temperatures and coverages. We concentrate on the coverages of 0.3 monolayer O or less, with variable sub-monolayer coverages of Cu. The results show that even though the ab initio calculations yield a higher barrier for O than for Cu adatom diffusion on Cu(1 0 0), the stability of Cu structures causes the O adatoms to be more mobile on the Cu(1 0 0) surface than the Cu adatoms. We are able to reproduce the c(2 × 2)-O domains seen in the experiments. However, we give an alternative explanation based on the repulsive interactions of O that, on one hand, cause the local ordering and, on the other hand, prohibits large well-ordered domains. We also give interpretation on the formation of the R45°-O reconstruction of Cu(1 0 0) above the O coverages of 0.3 monolayer based on the ab initio energetics.     

Coadsorption of methyl radicals and oxygen on Rh(1 1 1)

The chemisorption of CH₃ on Rh(1 1 1) is studied to understand the origin of the weakened symmetric stretch mode. A few different explanations for this weakened mode have been suggested in previous studies. These include C-H bond depletion and donation into C-H anti-bond orbitals either in an upright or tilted geometry. We investigate these possibilities by performing first-principles density functional calculations. Our results show strong adsorption at all high-symmetry sites with methyl in two possible orientations. A thorough analysis of the adsorption geometry shows that C_3v symmetry is preferred over a tilted species, ruling out tilting as a mechanism for C-H mode softening. Evidence of a multi-center bond between methyl and the surface rhodium atoms (similar to the kind shown recently by Michaelides and Hu for methyl on Ni(1 1 1)) is presented, showing that C-H bond depletion is the cause of mode-softening for methyl on Rh(1 1 1). Experimental results have shown that mode-softening diminishes when an electronegative species is coadsorbed, suggesting that donation into C-H anti-bonding orbitals is the mechanism for mode-softening. We therefore examine the coadsorption of oxygen and methyl on Rh(1 1 1). Our results suggest a new model for the effect of O on CH₃. Analysis of charge density differences shows that the dominant initial effects of O coadsorption are the removal of charge from the C-surface bond and the transfer of charge to the C-H bond. Subsequent increase of the H-Rh distance further reduces mode softening.     

The interaction of carbon with Si(1 1 1):As and Si(1 1 1):H surfaces, a theoretical study

Density functional theory calculations have been performed to determine the adsorption site of carbon at the Si(1 1 1):As and Si(1 1 1):H surfaces at different coverages. The As- and H-passivated surfaces were simulated by replacing the topmost Si layer by As or by saturating the Si dangling bonds with hydrogen atoms, respectively. Different high symmetry sites were considered. Carbon was placed successively in the fourfold (T₄) or threefold coordinated (H₃), the ontop (T₁) sites or substituted for a Si atom in the S₅ position located underneath the Si adatom in the T₄ site. We found that the preferred carbon adsorption site depends on the coverage of the passivated surfaces. At low coverages i.e. at 1/16 ML and 1/3 ML, it prefers a distorted T₄ position whereas at 1 ML, it occupies an H₃ site. This contrasts with the clean surface where the most energetically favored site is the S₅ at all coverages. Carbon adsorption induces a significant change in the structural geometry of the surface atoms, leading to a charge re-arrangement in the surface layers.     

Adsorption of Si on Ag(0 0 1) from ab initio study

We have systematically investigated the adsorption of Si on the Ag(0 0 1) surface employing density-functional theory. Various adsorption geometries have been considered for Si coverages up to 2.0 monolayers. Our results show that the behaviors of Si at the early stages of growth on the Ag(0 0 1) surface are governed by a competition between the Si-Si and Si-Ag interactions. From the calculated results, we presented alternative models for the observed 3 × 3 and 4 × 7 structures. Our results provide a reasonable explanation for the experimental findings in a previous work.     

Dissociative adsorption and adsorbate-induced reconstruction on Fe{2 1 1}

The products of CO, NO, O₂ and N₂ dissociation on Fe{2 1 1} have been studied by means of first-principles density functional theory. Preferred adsorption sites for adatoms C, N and O are identified, and trends in charge transfer and surface magnetism described. An experimentally observed (2 × 1) reconstruction induced by O is confirmed to be energetically stable, and a similar reconstruction induced by N is tentatively predicted. It is argued that these reconstructions may be important not only in the context of the catalytic reactivity of the Fe{2 1 1} surface, but also for the initial stages of surface nitridation and oxidation.     

Stabilization of the ferromagnetic configuration through alloying for monolayers of V, Cr, Mn on Pt(0 0 1) substrate

The magnetism is investigated for pseudomorphic monolayers of (V, Cr and Mn) and their substitutional c(2 × 2) alloyed structures adsorbed on Pt(0 0 1) semi-infinite surfaces using the semiempirical tight-binding method. Vanadium pseudomorphic monolayer exhibits a c(2 × 2) magnetic ordering, whereas Cr and Mn overlayers exhibit both c(2 ×  2) and p(1 × 1) magnetic structures. The substitutional alloyed overlayers exhibit ferrimagnetic structures that leads to appreciable net surface magnetization.     

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.     

A theoretical study of hydrogen adsorption and diffusion on a W(1 1 0) surface

We have used density functional theory to investigate hydrogen adsorption and diffusion on a W(1 1 0) surface. Hydrogen adsorption structures were examined from low coverage to one monolayer, and a threefold hollow site was found to be the most stable site at all coverages. In contrast to previous assertions, the work function decrease is not due to electron transfer from the hydrogen atoms to the W surface, but due to electron depletion at the vacuum region above the hydrogen atoms. Hydrogen atoms can diffuse via short-bridge sites and long-bridge sites at a coverage of θ = 1.0. Although the calculated activation energy for hydrogen diffusion via a short-bridge site is as small as 0.05 eV, field emission microscope experiments have shown that the activation energy for hydrogen diffusion is about 0.20 eV, which agrees fairly well with our calculated value of the activation energy via a long-bridge site. This discrepancy can be related to hydrogen delocalization on the W(1 1 0) surface, which has been suggested by electron energy loss spectroscopy experiments.     

Perpendicular magnetic anisotropy of ultrathin FeCo alloy films on Pd(0 0 1) surface: First principles study

Using the full potential linearized augmented plane wave (FLAPW) method, thickness dependent magnetic anisotropy of ultrathin FeCo alloy films in the range of 1 monolayer (ML) to 5 ML coverage on Pd(0 0 1) surface has been explored. We have found that the FeCo alloy films have close to half metallic state and well-known surface enhancement in thin film magnetism is observed in Fe atom, whereas the Co has rather stable magnetic moment. However, the largest magnetic moment in Fe and Co is found at 1 ML thickness. Interestingly, it has been observed that the interface magnetic moments of Fe and Co are almost the same as those of surface elements. The similar trend exists in orbital magnetic moment. This indicates that the strong hybridization between interface FeCo alloy and Pd gives rise to the large magnetic moment. Theoretically calculated magnetic anisotropy shows that the 1 ML FeCo alloy has in-plane magnetization, but the spin reorientation transition (SRT) from in-plane to perpendicular magnetization is observed above 2 ML thickness with huge magnetic anisotropy energy. The maximum magnetic anisotropy energy for perpendicular magnetization is as large as 0.3 meV/atom at 3 ML film thickness with saturation magnetization of . Besides, the calculated X-ray magnetic circular dichroism (XMCD) has been presented.     

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.     

The dehydrogenation of CH4 on Rh(1 1 1), Rh(1 1 0) and Rh(1 0 0) surfaces: A density functional theory study

CH₄ dehydrogenation on Rh(1 1 1), Rh(1 1 0) and Rh(1 0 0) surfaces has been investigated by using density functional theory (DFT) slab calculations. On the basis of energy analysis, the preferred adsorption sites of CH_x (x = 0-4) and H species on Rh(1 1 1), Rh(1 1 0) and Rh(1 0 0) surfaces are located, respectively. Then, the stable co-adsorption configurations of CH_x (x = 0-3) and H are obtained. Further, the kinetic results of CH₄ dehydrogenation show that on Rh(1 1 1) and Rh(1 0 0) surfaces, CH is the most abundant species for CH₄ dissociation; on Rh(1 1 0) surface, CH₂ is the most abundant species, our results suggest that Rh catalyst can resist the carbon deposition in the CH₄ dehydrogenation. Finally, results of thermodynamic and kinetic show that CH₄ dehydrogenation on Rh(1 0 0) surface is the most preferable reaction pathway in comparison with that on Rh(1 1 1) and Rh(1 1 0) surfaces.     

Adsorption of Al on the Si(0 0 1) surface

The first-principles calculations have been presented to study the adsorption of aluminum (Al) on the Si(0 0 1)(2×1) surface. We have investigated the optimized geometries and electronic structures of the adatom-substrate system. The adsorption energy of the system has been calculated. The most stable adsorption sites were consequently determined to be HH site and T3+T4. It is shown that the Si-Si dimer is asymmetric on the reconstructed bare surface and become symmetric upon Al adsorption. In addition, the bond length of Si-Si was found to be considerably elongated in the adsorption system. It is found that the work function change obtained in our work is different from other previous results on the adsorption of alkali metals on the Si(0 0 1) surface. In order to investigate the relative stability of phases at different coverages, the surface formation energy of the adsorption system was calculated. To shed light on the nature of the Al-Si bond and the character of silicon surface, the density of states (DOS) and difference charge density of the system were evaluated.     

Scaling behavior of island density in submonolayer growth of CaF2 on vicinal Si(1 1 1)

We have studied the scaling behavior of two-dimensional island density during submonolayer growth of CaF₂ on vicinal Si(1 1 1) surfaces using scanning tunneling microscopy. We have analyzed the morphology of the Si(1 1 1) surfaces where CaF₂ partial monolayers with coverages of about 0.1 monolayer are deposited at ∼600 °C. The number density of terrace nucleated islands increases with substrate terrace width l as ∼l⁴ in a low island density regime. This scaling behavior is consistent with predictions for the case of the irreversible growth of islands.     

Ground-state properties of arsenic deposited on the Ge(0 0 1)(2 × 1) surface

We present results of ab initio calculations of structural, electronic and vibrational properties of the Ge(0 0 1) surface covered with a monolayer of arsenic. The fully occupied π_u bonding and π_g antibonding electronic states due to the As-As dimer formation are quite close in energy and their ordering is same as that found on the Si(0 0 1) surface. Using our calculated atomic and electronic structures, surface lattice dynamics was studied by employing a linear response approach based on density functional perturbation theory. A comparison of the phonon spectrum of the Ge(0 0 1)/As(2 × 1) surface with that of the clean Ge(0 0 1)(2 × 1) surface indicates the presence of several new characteristic phonon modes due to adsorption of As atoms.     

Robust half-metallicity at the zincblende CrTe(0 0 1) surfaces and its interface with ZnTe(0 0 1)

All electron full potential calculations based on spin density functional theory are performed to study cubic zincblende (ZB) and hexagonal NiAs structures of bulk CrTe, free (0 0 1) surfaces of ZB CrTe, and interface of ZB CrTe with ZnTe(0 0 1). The ferromagnetic NiAs structure is reported to be about 0.26 eV more stable than the ferromagnetic ZB phase while ZB CrTe is found to be a half-metallic ferromagnet with a half-metallic gap of about 2.90 eV. Thermodynamic stability of CrTe(0 0 1) surfaces are studied in the framework of ab-initio thermodynamic. The obtained phase diagram evidences more stability of the Te terminated surface compared with the Cr termination. We discuss that both Te and Cr ideal terminations of CrTe(0 0 1) retain bulk-like half-metallic property but with a reduced half-metallic gap compared with bulk value. The structural, electronic, magnetic, and band alignment properties of the ZB CrTe/ZnTe(0 0 1) interface are computed and a rather large minority valence band offset of about 1.09 eV is observed in this heterojunction.     

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.