First-principle Wannier functions and effective lattice fermion models for narrow-band compounds

We propose a systematic procedure for constructing effective lattice fermion models for narrow-band compounds on the basis of first-principles electronic structure calculations. The method is illustrated for the series of transition-metal (TM) oxides: SrVO$_3$, YTiO$_3$, V$_2$O$_3$, and Y$_2$Mo$_2$O$_7$. It consists of three parts, starting from LDA. (i) construction of the kinetic energy Hamiltonian using downfolding method. (ii) solution of an inverse problem and construction of the Wannier functions (WFs) for the given kinetic energy Hamiltonian. (iii) calculation of screened Coulomb interactions in the basis of \textit{auxiliary} WFs, for which the kinetic-energy term is set to be zero. The last step is necessary in order to avoid the double counting of the kinetic-energy term, which is included explicitly into the model. The screened Coulomb interactions are calculated in a hybrid scheme. First, we evaluate the screening caused by the change of occupation numbers and the relaxation of the LMTO basis functions, using the conventional constraint-LDA approach, where all matrix elements of hybridization involving the TM $d$ orbitals are set to be zero. Then, we switch on the hybridization and evaluate the screening associated with the change of this hybridization in RPA. The second channel of screening is very important, and results in a relatively small value of the effective Coulomb interaction for isolated $t_{2g}$ bands. We discuss details of this screening and consider its band-filling dependence, frequency dependence, influence of the lattice distortion, proximity of other bands, and the dimensionality of the model Hamiltonian.Comment: 35 pages, 25 figure

Towards a first-principles theory of surface thermodynamics and kinetics

Understanding of the complex behavior of particles at surfaces requires detailed knowledge of both macroscopic and microscopic processes that take place; also certain processes depend critically on temperature and gas pressure. To link these processes we combine state-of-the-art microscopic, and macroscopic phenomenological, theories. We apply our theory to the O/Ru(0001) system and calculate thermal desorption spectra, heat of adsorption, and the surface phase diagram. The agreement with experiment provides validity for our approach which thus identifies the way for a predictive simulation of surface thermodynamics and kinetics.Comment: 4 pages including 3 figures. Related publications can be found at http://www.fhi-berlin.mpg.de/th/paper.htm

Orbital densities functional

Local density approximation (LDA) to the density functional theory (DFT) has continuous derivative of total energy as a number of electrons function and continuous exchange-correlation potential, while in exact DFT both should be discontinuous as number of electrons goes through an integer value. We propose orbital densities functional (ODF) (with orbitals defined as Wannier functions) that by construction obeys this discontinuity condition. By its variation one-electron equations are obtained with potential in the form of projection operator. The operator increases a separation between occupied and empty bands thus curing LDA deficiency of energy gap value systematic underestimation. Orbital densities functional minimization gives ground state orbital and total electron densities. The ODF expression for the energy of orbital densities fluctuations around the ground state values defines ODF fluctuation Hamiltonian that allows to treat correlation effects. Dynamical mean-field theory (DMFT) was used to solve this Hamiltonian with quantum Monte Carlo (QMC) method for effective impurity problem. We have applied ODF method to the problem of metal-insulator transition in lanthanum trihydride LaH_{3-x}. In LDA calculations ground state of this material is metallic for all values of hydrogen nonstoichiometry x while experimentally the system is insulating for x < 0.3. ODF method gave paramagnetic insulator solution for LaH_3 and LaH_{2.75} but metallic state for LaH_{2.5}.Comment: 35 pages, 5 figure

Lattice Distortion and Magnetism of 3d-t2gt_{2g} Perovskite Oxides

Several puzzling aspects of interplay of the experimental lattice distortion and the the magnetic properties of four narrow $t_{2g}$-band perovskite oxides (YTiO$_3$, LaTiO$_3$, YVO$_3$, and LaVO$_3$) are clarified using results of first-principles electronic structure calculations. First, we derive parameters of the effective Hubbard-type Hamiltonian for the isolated $t_{2g}$ bands using newly developed downfolding method for the kinetic-energy part and a hybrid approach, based on the combination of the random-phase approximation and the constraint local-density approximation, for the screened Coulomb interaction part. Then, we solve the obtained Hamiltonian using a number of techniques, including the mean-field Hartree-Fock (HF) approximation, the second-order perturbation theory for the correlation energy, and a variational superexchange theory. Even though the crystal-field splitting is not particularly large to quench the orbital degrees of freedom, the crystal distortion imposes a severe constraint on the form of the possible orbital states, which favor the formation of the experimentally observed magnetic structures in YTiO$_3$, YVO_, and LaVO$_3$ even at the HF level. Beyond the HF approximation, the correlations effects systematically improve the agreement with the experimental data. Using the same type of approximations we could not reproduce the correct magnetic ground state of LaTiO$_3$. However, we expect that the situation may change by systematically improving the level of approximations for dealing with the correlation effects.Comment: 30 pages, 17 figures, 8 tables, high-quality figures are available via e-mai

Oxygen impurities in NiAl: Relaxation effects

We have used a full-potential linear muffin-tin orbital method to calculate the effects of oxygen impurities on the electronic structure of NiAl. Using the supercell method with a 16-atom supercell we have investigated the cases where an oxygen atom is substitutionally placed at either a nickel or an aluminum site. Full relaxation of the atoms within the supercell was allowed. We found that oxygen prefers to occupy a nickel site over an aluminum site with a site selection energy of 138 mRy (21,370 K). An oxygen atom placed at an aluminum site is found to cause a substantial relaxation of its nickel neighbors away from it. In contrast, this steric repulsion is hardly present when the oxygen atom occupies the nickel site and is surrounded by aluminum neighbors. We comment on the possible relation of this effect to the pesting degradation phenomenon (essentially spontaneous disintegration in air) in nickel aluminides.Comment: To appear in Phys. Rev. B (Aug. 15, 2001

Adlayer core-level shifts of random metal overlayers on transition-metal substrates

We calculate the difference of the ionization energies of a core-electron of a surface alloy, i.e., a B-atom in a A_(1-x) B_x overlayer on a fcc-B(001)-substrate, and a core-electron of the clean fcc-B(001) surface using density-functional-theory. We analyze the initial-state contributions and the screening effects induced by the core hole, and study the influence of the alloy composition for a number of noble metal-transition metal systems. Data are presented for Cu_(1-x)Pd_x/Pd(001), Ag_(1-x) Pd_x/Pd(001), Pd_(1-x) Cu_x/Cu(001), and Pd_(1-x) Ag_x/Ag(001), changing x from 0 to 100 %. Our analysis clearly indicates the importance of final-state screening effects for the interpretation of measured core-level shifts. Calculated deviations from the initial-state trends are explained in terms of the change of inter- and intra-atomic screening upon alloying. A possible role of alloying on the chemical reactivity of metal surfaces is discussed.Comment: 4 pages, 2 figures, Phys. Rev. Letters, to appear in Feb. 199

Screened Coulomb interactions in metallic alloys: II Screening beyond the single-site and atomic sphere approximations

A quantitative description of the configurational part of the total energy of metallic alloys with substantial atomic size difference cannot be achieved in the atomic sphere approximation: It needs to be corrected at least for the multipole moment interactions in the Madelung part of the one-electron potential and energy. In the case of a random alloy such interactions can be accounted for only by lifting the atomic sphere and single-site approximations, in order to include the polarization due to local environment effects. Nevertheless a simple parameterization of the screened Coulomb interactions for the ordinary single-site methods, including the generalized perturbation method, is still possible. We obtained such a parameterization for bulk and surface NiPt alloys, which allows one to obtain quantitatively accurate effective interactions in this system.Comment: 24 pages, 2 figure

Onset of magnetism in B2 transition metals aluminides

Ab initio calculation results for the electronic structure of disordered bcc Fe(x)Al(1-x) (0.4<x<0.75), Co(x)Al(1-x) and Ni(x)Al(1-x) (x=0.4; 0.5; 0.6) alloys near the 1:1 stoichiometry, as well as of the ordered B2 (FeAl, CoAl, NiAl) phases with point defects are presented. The calculations were performed using the coherent potential approximation within the Korringa-Kohn-Rostoker method (KKR-CPA) for the disordered case and the tight-binding linear muffin-tin orbital (TB-LMTO) method for the intermetallic compounds. We studied in particular the onset of magnetism in Fe-Al and Co-Al systems as a function of the defect structure. We found the appearance of large local magnetic moments associated with the transition metal (TM) antisite defect in FeAl and CoAl compounds, in agreement with the experimental findings. Moreover, we found that any vacancies on both sublattices enhance the magnetic moments via reducing the charge transfer to a TM atom. Disordered Fe-Al alloys are ferromagnetically ordered for the whole range of composition studied, whereas Co-Al becomes magnetic only for Co concentration >0.5.Comment: 11 pages with 9 embedded postscript figures, to be published in Phys.Rev.