Comparison of O adsorption on Ni3Al (0 0 1), (0 1 1), and (1 1 1) surfaces through first-principles calculations

First-principles calculations were performed to study the properties of O adsorption on Ni₃Al (0 0 1), (0 1 1), and (1 1 1) surfaces using the Cambridge serial total package (CASTEP) code. Stable adsorption sites are identified. The atomic and electronic structures and adsorption energies are predicted. The adsorption sites for O on the Ni₃Al (0 0 1) surface are at the 2Ni–2Al fourfold hollow site, whereas O prefers to adsorb at the Ni–Al bridge site on (0 1 1) surface and 2Ni–Al threefold hollow site on (1 1 1) surface. It is found that O shows the strongest affinity for Al and the state of O is the most stabilized when O adsorbs on (0 0 1) surface, while the affinity of O for Al on (0 1 1) surface is weaker than (0 0 1) surface, and (1 1 1) surface is the weakest. The stronger O and Al affinity indicates more stable Al₂O₃ when oxidized. The experiment has shown that the oxidation resistance of single crystal superalloy in different orientations improves in the order of (1 1 1), (0 1 1), and (0 0 1) surface, suggesting that the oxidation in different crystallographic orientations may be related to the affinity of O for Al in the surface.     

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.     

Theoretical study of Ti and Fe surface alloys on Al(0 0 1) substrate

Accurate density-functional calculations are performed to investigate the formation of Ti and Fe ultrathin alloys on Al(0 0 1) surface. It is demonstrated that a deposition of Ti monolayer on Al(0 0 1) substrate leads to the formation of Al₃Ti surface alloy with Ti atoms arranged according to the L1₂ stacking, distinct from the D0₂₂ structure characteristic of a bulk Al₃Ti compound. A quest for the reason of this distinct atomic arrangement led us to the study of the surface structure of Al₃Ti(0 0 1) compound. It is concluded that even the Al₃Ti(0 0 1) surface is terminated with three layers assuming a L1₂ stacking and hence this stacking fault can be classified as a surface-induced stacking fault. Several possibilities of Fe atoms distributed in the surface region of Al(0 0 1) have been examined. The most stable configuration is the one with the compact Fe monolayer on Al(0 0 1) and covered by one Al monolayer. Lastly, our calculations show that there is no barrier for the penetration of Fe adatoms below the Al(0 0 1) surface; however, such a barrier is present for a Ti-alloyed Al(0 0 1) surface.     

Simulation of STM images of p(3 × 2)/Na/Ge(0 0 1) surface

Adsorption of Na on the Ge(0 0 1) surface is known to be a cause of surface reconstruction. It is expected to find one Na atom per unit cell of the reconstructed surface, however, the precise atomic configuration of this system is still a matter of controversy. Consequently, the aim of our present theoretical study is to examine the atomic structure of stable p(3 × 2)/Na/Ge(0 0 1) surfaces with and without the possible change of the number of Ge atoms in the surface layer (so-called mass transport). Structural and electronic properties of the considered system are investigated using the local-orbital density functional method. Our considerations are completed by a simulation of STM images of the structures following from molecular dynamics calculations.     

Anharmonic effects on B2-FeAl(1 1 0) surface: A molecular dynamics study

Using molecular dynamics simulations and the analytic embedded-atom method (AEAM), the surface anharmonicity of B2-FeAl(1 1 0) has been studied in the temperature range from 0 K to 1400 K. The temperature dependence of the interlayer spacing, mean square vibrational amplitudes, surface phonon frequencies and line-widths, and layer structure factor have been calculated. The obtained results indicate that the anharmonic effects are small in the temperature range from 0 K to 900 K. The temperature dependences of the interlayer spacing indicates that the rippling effect of the B2-FeAl(1 1 0) surface is exhibited by the contraction of Fe surface atoms and the expansion of Al atoms, which persists at high temperatures. The temperature dependence of the layer structure factors shows that the B2-FeAl(1 1 0) surface does not disorder until the temperature of 1300 K.     

Study of vacancy on diamond (1 0 0) (2 × 1) surface from first-principles

Structure and energy related properties of neutral and charged vacancies on relaxed diamond (1 0 0) (2 × 1) surface were investigated by means of density functional theory. Calculations indicate that the diffusion of a single vacancy from the top surface layer to the second layer is not energetically favored. Analysis of energies in charged system shows that neutral state is most stable on diamond (1 0 0) (2 × 1) surface. The multiplicity of possible states can exist on diamond (1 0 0) surface in dependence on the surface Fermi level, which supports that surface diffusion of a vacancy is mediated by the change of vacancy charge states. Analysis of density of states shows surface vacancy can be effectively measured by photoelectricity technology.     

A new structural model for the SiC(0 0 0 1)(3 × 3) surface derived from first principles studies

A new structural model with fluctuant Si-trimers and missing Si-adatom is proposed for Si-terminated 6H-SiC(0 0 0 1)(3 × 3) reconstruction. The atomic and electronic structures of the model are studied using first principles pseudopotential density-functional approach. The calculated surface electronic density of states coincides quantitatively with the experimental results of photoemission and electron energy loss spectroscopy. Based on the calculations, the Patterson map and scanning tunneling microscopic (STM) images simulated for the new model agree more satisfactorily with the experimental X-ray diffraction and STM observations than that for previously proposed models. The calculations of formation energies suggest that the new structure would be formed under the environment of dilute Si vapor around the surface in the preparation process.     

Atomic and electronic structure of Sr/Si(0 0 1)-(2 × 2)

The adsorption of Sr on the Si(0 0 1) surface with the semiantiphase dimer (2 × 2) reconstruction is studied, based upon the ab initio pseudopotential calculations. It is calculated that the semiantiphase dimer (2 × 2) reconstruction (2 dimers per unit cell) is more favorable than the (2 × 1) phase (1 dimer per unit cell) by an energy of about 0.24 eV/dimer. Considering the energetically more stable reconstruction, we have assumed four possible locations for 1/4 monolayer (ML) Sr adsorption on this surface: (i) bridge, (ii) cave, (iii) pedestal, and (iv) valley-bridge. We find that Sr adsorption on the valley-bridge site is energetically more favorable than all other cases studied here. Interestingly, one of the dimers becomes symmetric, but the other one is still asymmetric with the buckling angle reduced from 18° to 14°, when compared with the clean Si(0 0 1)-(2 × 2) surface. The calculated bond length between Sr and Si in the case of valley-bridge adsorption site is 3.05 Å, and in good agreement with other theoretical calculations. We also present and compare the electronic band structures for the clean and covered surfaces as well as the corresponding charge density plots.     

CO methanation on Ni(1 1 1) and modified Ni3Al(1 1 1) surfaces: A first-principle study

The adsorption energies of intermediates in CO methanation on the modified Ni₃Al(1 1 1) surface and the Ni(1 1 1) surface are calculated using density functional theory. A microkinetic analysis based on the calculated adsorption energies is performed to explain the different kinetics of CO methanation catalyzed by Ni₃Al and Ni powders. The electronic structures of different atoms on the modified Ni₃Al surface are also presented, and correlate well with the adsorption energies and geometries.     

Theoretical investigation of the Au/Si(1 1 1)-(5 × 2) surface structure

The atomic structure of the Au/Si(1 1 1)-(5 × 2) surface has been studied by density-functional theory calculations. Two structure models, proposed experimentally by Marks et al. and Hasegawa et al., have been examined on an equal ground. In our total-energy calculations, both models are found to be locally stable and energetically comparable. In our electronic-structure analyses, however, both models fail to reproduce the key features of angle-resolved photoemission spectra and scanning-tunneling-microscopy images, indicating that the considered models need to be modified. Suggestions for the modification are given based on the present calculations.     

Cyclic cluster study of water adsorption structures on the MgO(1 0 0) surface

Cyclic cluster calculations were performed with the quantum chemical method MSINDO to elucidate the relative stabilities of c(4 × 2), p(3 × 2) and (1 × 1) overlayer structures of water molecules on the MgO(1 0 0) surface. For the c(4 × 2) and p(3 × 2) structures both molecular adsorption and partially dissociated adsorption were considered. In agreement with earlier theoretical studies partial dissociation was found to be more stable than molecular adsorption. For the c(4 × 2) structure both monolayer and double layer coverage were studied. Adsorption was found to be more stabilized with increasing degree of dissociation until 50% of the water molecules were dissociated. In the case of 50% dissociated water molecules we found that one quarter of the Mg atoms were pulled out of the MgO surface when surface relaxation was taken into account. A new structure for the fully dissociated (1 × 1) water monolayer was found which is considerably more stable than previously studied arrangements. In all cases surface relaxation was found to be important. The most stable structures of c(4 × 2), p(3 × 2) and (1 × 1) symmetry have adsorption energies which differ by no more than 13 kJ/mol. This offers an explanation of phase transitions of overlayer structures found by experiments between 180 and 300 K.     

Surface energy of M2AC(0 0 0 1) determined by density functional theory (M = Ti, V, Cr; A = Al, Ga, Ge)

We have studied the correlation between the valence electron configuration and the electronic structure of M₂AC(0 0 0 1) surfaces (M = Ti, V, Cr; A = Al, Ga, Ge) by density functional theory. The A surface termination is the most stable configuration for all systems studied according to our surface energy data. As the M valence electron population is increased, the surface energy increases by 22% and 12% for A = Al and Ga, respectively, while it decreases by 29% for A = Ge. This can be understood by evaluating the valence electron concentration induced changes in the surface density of states. Antibonding surface Md-Ap states are present as Ti is substituted by Cr in M₂AC(0 0 0 1) for A = Al and Ga, while antibonding surface Md-Ap states are not present as Ti is substituted by Cr in M₂GeC(0 0 0 1).     

Effect of hydrogenation on P/Si(0 0 1)-(1 × 2)

Ab initio calculations, based on pseudopotentials and density functional theory, have been performed to investigate the effect of hydrogenation on the atomic geometries and energetics of substitutional phosphorus (P) on the generic Si(0 0 1)-(1 × 2) surface. For the 0.5 ML coverage of P, we have considered three different substitutional sites: (i) the mixed Si-P dimer structure (i.e., the P-nondiffused case), (ii) P-interdiffused to the second layer Si (i.e., intermixed P-Si bond structure) and (iii) P-interdiffused to the third layer Si. We have found that the mixed Si-P dimer structure is 0.79 eV/dimer energetically more favorable than the P-interdiffused case. However, for the hydrogenation of above cases, we have found that the situation is reversed and the interdiffused case is 0.3 eV/dimer energetically more favorable than the P-nondiffused case. Reductions in the number of P-Si is identified as a contributing factor which determines energetically the stable structures during P on Si(0 0 1).     

Surface segregation of aluminum atoms on Cu-9 at.% Al(1 1 1) studied by Auger electron spectroscopy and low energy electron diffraction

Auger electron spectroscopy (AES) and low energy electron diffraction (LEED) were applied to investigate the segregation of aluminum atoms on a Cu-9 at.% Al(1 1 1) surface. We observed that the Al concentration in the top layer ranged between about 9 and 36 at.% after the sample we used was annealed at different temperatures. The phenomenon of Al atoms segregating on the surfaces was explained well by considering the diffusion length of Al atoms in bulk Cu. LEED measurements showed that R30° structures grew as the concentration of Al atoms increased. The segregation phenomena on surfaces resulted in a stable two-dimensional Cu₆₇Al₃₃ alloy phase in the top layer.     

Simulations of iodine adsorbed Ge(0 0 1) surface and its STM images

A layer of iodine at Ge(0 0 1) surface develops an ordered structure of iodine atoms bound to Ge dimers. Here are discussed atomic structures of Ge(0 0 1) surface covered by 0.25 monolayer of iodine. The p(2×4), p(2×2), c(2×4) and p(1×4) surface structures are found in calculations. The structure with two iodine atoms of the dissociated I₂ molecule adsorbed at both ends of the same germanium dimer is found to be energetically favourable over iodine adsorption at neighbouring dimers. Simulated STM images of the obtained surface structures are presented and compared with experimental data.     

Growth kinetics of metastable (3 3 1) nanofacet on Au and Pt(1 1 0) surfaces

A theoretical epitaxial growth model with realistic barriers for surface diffusion is investigated by means of kinetic Monte Carlo simulations to study the growth modes of metastable (3 3 1) nanofacets on Au and Pt(1 1 0) surfaces. The results show that under experimental atomic fluxes, the (3 3 1) nanofacets grow by 2D nucleation at low temperature in the submonolayer regime. A metastable growth phase diagram that can be useful to experimentalists is presented and looks similar to the one found for the stationary growth of the bcc(0 0 1) surface in the kinetic 6-vertex model.     

Structural, electronic properties and stability of the (1 × 1) PbTiO3 (1 1 1) polar surfaces by first-principles calculations

Under GGA, the structural, electronic properties and stabilities of four different (1 × 1) terminations of cubic PbTiO₃ (1 1 1) surface, the directly cleaved (1 1 1)-Ti and (1 1 1)-PbO₃ terminations and the constructed (1 1 1)-TiO and (1 1 1)-PbO₂ ones, have been systematically studied by using projector-augmented wave method implemented in VASP. For (1 1 1)-Ti and (1 1 1)-PbO₃ terminations, Ti-O bonds between the outermost two layers are enhanced after relaxation, while those between the second and the third layers are weakened. In addition, a contraction of O-O distance in surface PbO₃ layer is also found for (1 1 1)-PbO₃ termination. Moreover, electronic structures of both (1 1 1)-Ti and (1 1 1)-PbO₂ terminations are significantly influenced by structure relaxations, and the effects of the surface on the DOS are dominantly on the Ti layers, especially the CB. For a constructed (1 1 1)-TiO termination, the relaxation results show both Ti-O bonds between the outermost two layers and those between the second PbO₃ layer and the third Ti layers are enhanced. For a constructed (1 1 1)-PbO₂ termination, Ti-O bonds between the outermost two layers are also enhanced as in the (1 1 1)-TiO termination, however, inequivalent Ti-O bonds between the second layer Ti atom and the third layer O atoms are found, with two bonds expanding and the other one contracting. Results of electronic structure calculations show these two constructed terminations are all insulating and changes of DOS originate dominantly from modifications of surface compositions. Furthermore, it is found that for all four different (1 1 1) terminations, the movements of the cation and/or anion on the outermost layer along the surface normal direction after relaxation all result in a reduction of the space electric field. In O and Pb external environments, it is predicted that (1 1 1)-PbO₂ termination is the most stable one in O- and Pb-rich environments, however, the (1 1 1)-Ti termination is stable one in O- and Pb-poor conditions. The (1 1 1)-TiO termination also shows a stability domain in moderate O and Pb environments.     

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.     

Structure and reactivity of the hydrated hematite (0 0 0 1) surface

The structure of the hydroxylated hematite (0 0 0 1) surface was investigated using crystal truncation rod diffraction and density functional theory. The combined experimental and theoretical results suggest that the surface is dominated by two hydroxyl moieties—hydroxyls that are singly coordinated and doubly coordinated with Fe. The results are consistent with the formation of distinct domains of these surface species; one corresponding to the hydroxylation of the surface Fe-cation predicted to be most stable under UHV conditions, and the second a complete removal of this surface Fe species leaving the hydroxylated oxygen layer. Furthermore, our results indicate that the hydroxylated hematite surface structures are significantly more stable than their dehydroxylated counterparts at high water partial pressures, and this transition in stability occurs at water pressures orders of magnitude below the same transition for α-alumina. These results explain the observed differences in reactivity of hematite and alumina (0 0 0 1) surfaces with respect to water and binding of aqueous metal cations.     

Vacancy charging on Si(1 1 1)-(7 × 7) investigated by density functional theory

The structure and energetics of charged vacancies on Si(1 1 1)-(7 × 7) are investigated using density functional theory calculations supplemented by estimates of ionization entropy. The calculations predict multiple possible charge states for the unfaulted edge vacancy in the adatom layer, although the −2 state is most stable on real Si(1 1 1) surfaces for which the Fermi level lies near the middle of the band gap.