First-principles calculations of hydrogen molecule adsorption on Ti (0 0 0 1)-(2 × 1) surface

Adsorption of H₂ molecule on the Ti (0 0 0 1)-(2 × 1) surface was studied by density functional theory with generalized gradient approximation (GGA). The parallel and vertical absorption cases were investigated in detail by adsorption energy and electronic structure analysis, we obtained three stable configurations of FCC-FCC (the two H atoms adsorption on the two adjacent fcc sites of Ti (0 0 0 1) surface, respectively), HCP-HCP (the two H atoms adsorption on the two adjacent hcp sites of Ti (0 0 0 1) surface, respectively) and FCC-HCP (the one H atom adsorption on the fcc site and the other adsorption on the near hcp site) based on the six different parallel adsorption sites after the H₂ molecule dissociates. However, all the end configurations of four vertical adsorption sites were unstable, H₂ molecule was very easy to desorb from Ti surface. The H-H bond breaking and Ti-H bond forming result from the H₂ molecule dissociation. H-H bond breaking length ranges from 1.9 Å to 2.3 Å for different adsorption configurations due to the strong Ti-H bond forming. The H₂ dissociative approach and the end stable configurations formation in parallel adsorption processes are attributed to the quantum mechanics steering effects.     

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.     

Oxygen vacancy pairs on CeO2(110): A DFT + U study

Oxygen vacancy pairs have been suggested to play a role in the reduction of NO molecules on ceria and for the oxidation processes of reducible rare-earth oxides. The formation energy of the oxygen vacancy pairs and the changes in the structural and electronic properties of the ceria (110) surface with oxygen vacancy pairs are investigated using density-functional theory (DFT + U) methodology within the generalized gradient approximation. It is found that the excess electrons localize on the Ce ions neighbouring the vacancies, and the most stable structure for the oxygen vacancy pairs on the ceria (110) surface is at next-nearest-neighbour site.     

Electronic Structure of Si1-xⅣx/Si Superlattices on Si （001）

We have preformed systematical ab initio studies of the structural and electronic properties of short-period Si1-xⅣx/Si （x = 0.125, 0.25, 0.5, Ⅳ=Ge, Sn） superlattices （SLs） grown along the [001] direction on bulk Si. The present calculations reveal that the Si0.875Ge0.125/Si, Si0.75Ge0.25/Si and Si0.875Sn0.125/Si are the F-point direct bandgap semiconductors. The technological importance lies in the expectation that the direct gap Si1-xⅣx/Si SLs may be used as components in integrated optoelectronic devices, in conjunction with the already well-established and highly advanced silicon technology.     

Adsorption structure of germanium on the Ru(0 0 0 1) surface

Coverage-dependent adsorption energy of the Ge/Ru(0 0 0 1) growth system and the geometrical distortions of the most stable adsorption structure are investigated through first-principles calculations within density functional theory. A local minimum in adsorption energy is found to be at a Ge coverage of 1/7 monolayer with a Ru(0 0 0 1)- symmetry. Based on this stale superstructure, the scanning tunneling microscopy (STM) and scanning tunneling spectroscopy (STS) images are simulated by means of surface local-density of states (LDOS). The results are consistent well with the STM measurements on the phase for Ge overlayer on Ru(0 0 0 1). From this stimulation, the relations between the STM images and the lattice distortion are also clarified.     

Study of arsenic for antimony exchange at the Sb-stabilized GaSb(0 0 1) surface

In this paper we present a first-principle study on the energetics of a single As₂ molecule on GaSb(0 0 1) reconstructed surface. In order to shed light into the mechanisms of anion exchange at the Sb-rich GaSb(0 0 1) surface, we studied firstly As₂ adsorption and then As for Sb exchange. We identify a surface region where both the processes are energetically favored. The results of this twofold analysis can be combined to derive possible reaction paths for the anion exchange process.     

Ruthenium adsorption and diffusion on the GaN(0 0 0 1) surface

We report first principles calculations to analyze the ruthenium adsorption and diffusion on GaN(0 0 0 1) surface in a 2×2geometry. The calculations were performed using the generalized gradient approximation (GGA) with ultrasoft pseudopotential within the density functional theory (DFT). The surface is modeled using the repeated slabs approach. To study the most favorable ruthenium adsorption model we considered T₁, T₄ and H₃ special sites. We find that the most energetically favorable structure corresponds to the Ru- T₄ model or the ruthenium adatom located at the T₄ site, while the ruthenium adsorption on top of a gallium atom (T₁ position) is totally unfavorable. The ruthenium diffusion on surface shows an energy barrier of 0.612 eV. The resultant reconstruction of the ruthenium adsorption on GaN(0 0 0 1)- 2×2 surface presents a lateral relaxation of some hundredth of Å in the most stable site. The comparison of the density of states and band structure of the GaN(0 0 0 1) surface without ruthenium adatom and with ruthenium adatom is analyzed in detail.     

Study on the adsorption of fluorescein on Ag(1 1 0) substrate

The adsorption of fluorescein on the Ag(1 1 0) surface has been investigated by the first-principles pseudopotential method. Various adsorption geometries have been calculated and the energetically most favorable structure of fluorescein/Ag(1 1 0) was identified. The fluorescein molecule, in most favorable structure, is on hollow site, and the adsorption energy is 2.34 eV. Here the adsorption sites refer to the positions at the first layer of the substrate where the middle carbon atom of the fluorescein molecule is located. The bonding strength of the fluorescein molecule to the Ag substrate is site selective, being determined by electron transfer to the oxygen atoms of the molecule and local electrostatic attraction between the oxygen atoms and the silver atoms.     

Density functional theory study of hydrogenation mechanism in Fe-doped Mg(0 0 0 1) surface

Using density functional theory (DFT) in combination with nudged elastic band (NEB) method, the dissociative chemisorptions and diffusion processes of hydrogen on both pure and Fe-doped Mg(0 0 0 1) surfaces are studied. Firstly, the dissociation pathway of H₂ and the relative barrier were investigated. The calculated dissociation barrier (1.08 eV) of hydrogen molecule on a pure Mg(0 0 0 1) surface is in good agreement with comparable experimental and theoretical studies. For the Fe-doped Mg(0 0 0 1) surface, the activated barrier decreases to 0.101 eV due to the strong interaction between the s orbital of H and the d orbital of Fe. Then, the diffusion processes of atomic hydrogen on pure and Fe-doped Mg(0 0 0 1) are presented. The obtained diffusion barrier to the first subsurface is 0.45 eV and 0.98 eV, respectively. Finally, Chou method was used to investigate the hydrogen sorption kinetic mechanism of pure MgH₂ and Mg mixed with 5 at.% Fe atoms composites. The obtained activation energies are 0.87 ± 0.02 and 0.31 ± 0.01 eV for H₂ dissociation on the pure surface and H atom diffusion in Fe-doped Mg surfaces, respectively. It suggests that the rate-controlling step is dissociation of H₂ on the pure Mg surface while it is diffusion of H atom in the Fe-doped Mg surface. And both of fitting data are matching well with our calculation results.     

Ab initio prediction of half-metallic ferromagnetism in Zn(TM)O2 (TM=Cr,Mn, Fe,Co, Ni) compounds

From the results of first principles tight-binding linear muffin-tin orbital (TB-LMTO) calculations, half-metallic ferromagnetism is proposed in Zn(TM)O₂ with a chalcopyrite structure. The calculated electronic and magnetic property shows that consistent with the integer value for the total magnetic moment, half metallicity is obtained for ZnCrO₂, ZnMnO₂, ZnFeO₂, ZnCoO₂ and ZnNiO₂. A careful analysis of the spin density reveals the ferromagnetic coupling between the p–d states and the cation dangling-bond p states, which is believed to be responsible for the stabilization of the ferromagnetic phase. The calculated heat of formation, bulk modulus and cohesive energy are reported.     

Mn and Fe impurities in MgB2

Based on first principles calculations, we show that Mn impurities are magnetic in MgB₂ due to exchange-splitting of band. Thus, Mn impurities could act as strong magnetic scattering centres leading to pair-breaking effects in MgB₂. In contrast, we find Fe impurities in MgB₂ to be nearly non-magnetic.     

Local environment of nitrogen in GaNyAs1−y epilayers on GaAs (0 0 1) studied using X-ray absorption near edge spectroscopy

X-ray absorption near-edge spectroscopy (XANES) is used to study the N environment in bulk GaN and in GaN_yAs_1−y epilayers on GaAs (0 0 1), for y∼5%. Density-functional optimized structures were used to predict XANES via multiple-scattering theory. We obtain striking agreement for pure GaN. An alloy model with nitrogen pairs on Ga accurately predicts the threshold energy, the width of the XANES ‘white line’, and features above threshold, for the given X-ray polarization. The presence of large quantitities of N-pairs may point to a role for molecular N₂ in epitaxial growth kinetics.     

Theoretical study of CuxAg1−xI alloys

We have performed first-principles calculations using full potential linearized augmented plane wave (FP-LAPW) method within density functional theory (DFT) to investigate the fundamental properties of Cu_xAg_1−xI alloys. We used both GGA96 [J.P. Perdew, K. Burke, M. Ernzerhof, Phys. Rev. Lett. 77 (1996) 3865.] and EVGGA [E. Engel, S.H. Vosko, Phys. Rev. B. 47 (1993) 13164.] generalized gradient approximations of the exchange-correlation energy that are based on the optimization of total energy and corresponding potential. Quantities such as lattice constants, bulk modulus, band gap, density of occupied states and effective mass were calculated as a function of copper molar fraction x. These parameters were found to depend non-linearly on alloy composition x, except the lattice parameter, which follows Vegard's law. The microscopic origins of the gap bowing were explained using the approach of Zunger and co-workers; we have concluded that the band-gap energy bowing was mainly caused by the chemical charge-transfer effect and the volume deformation , while the structural relaxation contribute to the gap bowing parameter at smaller magnitude. The calculated phase diagram shows a broad miscibility gap for this alloy with a high critical temperature.     

First principles study of structural phase transition of YSb and ScSb compounds

High pressure induced phase transition of YSb and ScSb compounds have been studied using Density Functional Theory method within Generalized Gradient Approximation. It was found that the phase transition from the NaCl-type (B1) to a CsCl-type structure (B2) began to occur at around 29 GPa for YSb compound, agreeing well with available experiments and theoretical calculations. For ScSb compound it was suggested that structural phase transition from B1 to B2 will occur at about 40 GPa, differing greatly with experimental and theoretical results. The finding that the transition pressures increase with decreasing lattice constant in the NaCl-type structure for YSb and ScSb compounds was found to be similar to the phenomena observed for LnSb (Ln: lanthanide) compounds. Mulliken charge and overlap population analysis revealed that YSb and ScSb compounds in B1 structure show similar interaction between anion and cation, while in B2 structure a higher degree of covalency was found for ScSb than that in YSb. Also, DOS and band structure of these two compounds in B1 and B2 structures were presented and analyzed.     

First-principles prediction of half-metallic ferromagnetism in Cd(TM)O2 (TM=Cr, Mn, Fe, Co, Ni) compounds

We report the electronic structure of Cd(TM)O₂ (TM=Cr, Mn, Fe, Co, Ni) in the chalcopyrite structures. From this study we find that Cd(TM)O₂ is a half-metallic ferromagnetic compound. From the energy consideration we find that Cd(TM)O₂ is more stable in chalcopyrite structure rather than in rock salt structure. A careful analysis of the spin density reveals the ferromagnetic coupling between the p-d states and the cation dangling-bond p states, which is believed to be responsible for the stabilization of the ferromagnetic phase. The calculated heat of formation, bulk modulus and cohesive energy are reported.     

Structural, elastic and electronic properties and formation energies for hexagonal (W0.5Al0.5)C in comparison with binary carbides WC and Al4C3 from first-principles calculations

First principle calculations have been performed with the purpose to understand the peculiarities of the structural, elastic parameters and electronic properties and interatomic bonding for novel hexagonal carbide (W_0.5Al_0.5)C in comparison with binary phases WC and Al₄C₃. The geometries of all phases were optimized and their structural, elastic parameters and theoretical density were established. Besides, we have evaluated the formation energies (E_form) of W_0.5Al_0.5C for different possible preparation routes (namely for the reactions with the participation of simple substances (metallic W, Al and graphite, binary W or Al carbides and metallic Al and W, or binary W and Al carbides). The results show that the synthesis of the ternary carbide from simple substances is more favorable in comparison with the reactions with participation of W and Al carbides. Moreover, band structures, total and partial densities of states were obtained and analyzed systematically for (W_0.5Al_0.5)C, WC and Al₄C₃ phases in comparison with available theoretical and experimental data. The bonding picture in W_0.5Al_0.5C was described as a mixture of metallic, ionic and covalent contributions with the high anisotropy for the covalent W-C and Al-C bonds, where p-p like Al-C bonds become weaker than p-d like W-C bonds.     

Electronic structure and magnetism of MnFeP1−xSix alloys from first-principles calculations

First-principles calculations of electronic structure and magnetic properties based on density-functional theory were performed for MnFeP_1−xSi_x (0.44?x?0.60) alloys which are considered as promising magnetocaloric refrigerants. We used the full-potential APW+lo method and treated the random order of P(Si) atoms in the ZrNiAl-type structure in a virtual-crystal approximation. A non-monotonic behavior of the alloy magnetization as a function of x was obtained, in qualitative agreement with experiment, and explained in terms of the spin-polarized densities of states.     

On the propensity of magnetism in 3d transition-metal-MgCNi3 alloys

The effects of disorder and incipient magnetism in MgC(_Ni1-xTx)3$\mathrm{MgC}({\mathrm{Ni}}_{1-x}{\mathrm{T}}_x{)}_3$ (T≡Fe$\mathrm{T}\equiv \mathrm{Fe}$, Co or Cu) alloys are studied using coherent-potential approximation and Ginzburg–Landau coefficients. The first-principles, local-density-functional-based calculations for substitutionally disordered Fe and Co impurities in the Ni sub-lattice of MgCNi₃, in low concentrations, show that incipient magnetism resides in these materials. The overestimation of the calculated magnetic properties points to the limitations of the local-density approximation. However, using a phenomenological approach based on Ginzburg–Landau coefficients and the fixed-spin moment method, we show that MgC(_Ni1-xTx)3$\mathrm{MgC}({\mathrm{Ni}}_{1-x}{\mathrm{T}}_x{)}_3$ alloys remain paramagnetic. At expanded volumes, we also find the possibility of a ferromagnetic state for MgC(_{Ni0.95Fe0.05)3}$\mathrm{MgC}({\mathrm{Ni}}_{0.95}{\mathrm{Fe}}_{0.05}{)}_3$ and MgC(_{Ni0.90Co0.10)3}$\mathrm{MgC}({\mathrm{Ni}}_{0.90}{\mathrm{Co}}_{0.10}{)}_3$ alloys.     

Theoretical investigations on KClxBr1−x, KClxI1−x and KBrxI1−x: A first-principles study

Using first-principles total energy calculations within the full-potential linearized augmented plane wave (FP-LAPW) method, we have investigated the structural, electronic and thermodynamic properties of potassium halides (KCl_xBr_1−x, KCl_xI_1−x and KBr_xI_1−x), with x concentrations varying from 0% up to 100%. The effect of composition on lattice constants, bulk modulus, band gap and dielectric function was investigated. Deviations of the lattice constants from Vegard's law and the bulk modulus from linear concentration dependence (LCD) were observed for the three alloys. The microscopic origins of the gap bowing were explained by using the approach of Zunger and coworkers. On the other hand, the thermodynamic stability of these alloys was investigated by calculating the excess enthalpy of mixing ΔH_m as well as the phase diagram.     

Ab initio study of pressure-induced magnetic transition in manganese pnictides

We report a density functional calculation on the NiAs-type Mn-based pnictides. The total energy as a function of volume is obtained by means of self-consistent tight-binding linear muffin–tin orbital method by performing spin and non-spin polarized calculation. From the present study, we predict a magnetic-phase transition from ferromagnetic (FM) to non-magnetic (NM) around 49 and 35.7 GPa for MnAs and MnSb, respectively. The pressure-induced transition is found to be a second-order transition. The band structure and density of states (DOS) are plotted for FM and NM states. Apart from this the ground-state properties like magnetic moment, lattice parameter and bulk modulus are calculated and are compared with the available results. Under large volume expansion these compounds exist in zinc-blende (ZB) structure, which shows half metallicity. The magnetic moment and equilibrium lattice constants for ZB structure are obtained as well as band structure and DOS are presented.