Localization and delocalization errors in density functional theory and implications for band-gap prediction

The band-gap problem and other systematic failures of approximate functionals are explained from an analysis of total energy for fractional charges. The deviation from the correct intrinsic linear behavior in finite systems leads to delocalization and localization errors in large or bulk systems. Functionals whose energy is convex for fractional charges such as LDA display an incorrect apparent linearity in the bulk limit, due to the delocalization error. Concave functionals also have an incorrect apparent linearity in the bulk calculation, due to the localization error and imposed symmetry. This resolves an important paradox and opens the possibility to obtain accurate band-gaps from DFT.Comment: 4 pages 4 figure

Nonempirical Density Functionals Investigated for Jellium: Spin-Polarized Surfaces, Spherical Clusters, and Bulk Linear Response

Earlier tests show that the Tao-Perdew-Staroverov-Scuseria (TPSS) nonempirical meta-generalized gradient approximation (meta-GGA) for the exchange-correlation energy yields more accurate surface energies than the local spin density (LSD) approximation for spin-unpolarized jellium. In this study, work functions and surface energies of a jellium metal in the presence of ``internal'' and external magnetic fields are calculated with LSD, Perdew-Burke-Ernzerhof (PBE) GGA, and TPSS meta-GGA and its predecessor, the nearly nonempirical Perdew-Kurth-Zupan-Blaha (PKZB) meta-GGA, using self-consistent LSD orbitals and densities. The results show that: (i) For normal bulk densities, the surface correlation energy is the same in TPSS as in PBE, as it should be since TPSS strives to represent a self-correlation correction to PBE; (ii) Normal surface density profiles can be scaled uniformly to the low-density or strong-interaction limit, and TPSS provides an estimate for that limit that is consistent with (but probably more accurate than) other estimates; (iii) For both normal and low densities, TPSS provides the same description of surface magnetism as PBE, suggesting that these approximations may be generally equivalent for magnetism. The energies of jellium spheres with up to 106 electrons are calculated using density functionals and compared to those obtained with Diffusion Quantum Monte Carlo data, including our estimate for the fixed-node correction. Finally we calculate the linear response of bulk jellium using these density functionals, and find that not only LSD but also PBE GGA and TPSS meta-GGA yield a linear-response in good agreement with that of the Quantum Monte Carlo method, for wavevectors of the perturbing external potential up to twice the Fermi wavevector.Comment: 14 pages, 9 figure

Climbing the Density Functional Ladder: Non-Empirical Meta-Generalized Gradient Approximation Designed for Molecules and Solids

The electron density, its gradient, and the Kohn-Sham orbital kinetic energy density are the local ingredients of a meta-generalized gradient approximation (meta-GGA). We construct a meta-GGA density functional for the exchange-correlation energy that satisfies exact constraints without empirical parameters. The exchange and correlation terms respect {\it two} paradigms: one- or two-electron densities and slowly-varying densities, and so describe both molecules and solids with high accuracy, as shown by extensive numerical tests. This functional completes the third rung of ``Jacob's ladder'' of approximations, above the local spin density and GGA rungs.Comment: 4 pages, 1 figure, 1 table. updated with minor and yet necessary corrections. New references are adde

A natural orbital functional for the many-electron problem

The exchange-correlation energy in Kohn-Sham density functional theory is expressed as a functional of the electronic density and the Kohn-Sham orbitals. An alternative to Kohn-Sham theory is to express the energy as a functional of the reduced first-order density matrix or equivalently the natural orbitals. In the former approach the unknown part of the functional contains both a kinetic and a potential contribution whereas in the latter approach it contains only a potential energy and consequently has simpler scaling properties. We present an approximate, simple and parameter-free functional of the natural orbitals, based solely on scaling arguments and the near satisfaction of a sum rule. Our tests on atoms show that it yields on average more accurate energies and charge densities than the Hartree Fock method, the local density approximation and the generalized gradient approximations

Two Avenues to Self-Interaction Correction within Kohn-Sham Theory: Unitary Invariance is the Shortcut

The most widely-used density functionals for the exchange-correlation energy are inexact for one-electron systems. Their self-interaction errors can be severe in some applications. The problem is not only to correct the self-interaction error, but to do so in a way that will not violate size-consistency and will not go outside the standard Kohn-Sham density functional theory. The solution via the optimized effective potential (OEP) method will be discussed, first for the Perdew-Zunger self-interaction correction (whose performance for molecules is briefly summarized) and then for the more modern self-interaction corrections based upon unitarily-invariant indicators of iso-orbital regions. For the latter approaches, the OEP construction is greatly simplified. The kinetic-energy-based iso-orbital indicator \tau^W_\sigma(\re)/\tau_\sigma(\re) will be discussed and plotted, along with an alternative exchange-based indicator

Ab initio pseudopotential study of Fe, Co, and Ni employing the spin-polarized LAPW approach

The ground-state properties of Fe, Co, and Ni are studied with the linear-augmented-plane-wave (LAPW) method and norm-conserving pseudopotentials. The calculated lattice constant, bulk modulus, and magnetic moment with both the local-spin-density approximation (LSDA) and the generalized gradient approximation (GGA) are in good agreement with those of all-electron calculations, respectively. The GGA results show a substantial improvement over the LSDA results, i.e., better agreement with experiment. The accurate treatment of the nonlinear core-valence exchange and correlation interaction is found to be essential for the determination of the magnetic properties of 3d transition metals. The present study demonstrates the successful application of the LAPW pseudopotential approach to the calculation of ground-state properties of magnetic 3d transition metals.Comment: RevTeX, 14 pages, 2 figures in uufiles for

Ab initio study of shock compressed oxygen

Quantum molecular dynamic simulations are introduced to study the shock compressed oxygen. The principal Hugoniot points derived from the equation of state agree well with the available experimental data. With the increase of pressure, molecular dissociation is observed. Electron spin polarization determines the electronic structure of the system under low pressure, while it is suppressed around 30 $\sim$ 50 GPa. Particularly, nonmetal-metal transition is taken into account, which also occurs at about 30 $\sim$ 50 GPa. In addition, the optical properties of shock compressed oxygen are also discussed.Comment: 5 pages, 5 figure

Intrinsic hole localization mechanism in magnetic semiconductors

The interplay between clustering and exchange coupling in magnetic semiconductors for the prototype (Ga_{1-x},Mn_x)As with manganese concentrations x of 1/16 and 1/32 in the interesting experimental range is investigated. For x ~ 6 %, when all possible arrangements of two atoms within a large supercell are considered, the clustering of Mn atoms at nearest-neighbour Ga sites is energetically preferred. As shown by spin density analysis, this minimum energy configuration localizes further one hole and reduces the effective charge carrier concentration. Also the exchange coupling constant increases to a value corresponding to lower Mn concentrations with decreasing inter Mn distance.Comment: Accepted for publication in Journal of Physics: Condensed Matte

A local density functional for the short-range part of the electron-electron interaction

Motivated by recent suggestions --to split the electron-electron interaction into a short-range part, to be treated within the density functional theory, and a long-range part, to be handled by other techniques-- we compute, with a diffusion Monte Carlo method, the ground-state energy of a uniform electron gas with a modified, short-range-only electron-electron interaction \erfc(\mu r)/r, for different values of the cutoff parameter $\mu$ and of the electron density. After deriving some exact limits, we propose an analytic representation of the correlation energy which accurately fits our Monte Carlo data and also includes, by construction, these exact limits, thus providing a reliable ``short-range local-density functional''.Comment: 7 pages, 3 figure