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

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

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

A More Accurate Generalized Gradient Approximation for Solids

We present a new nonempirical density functional generalized gradient approximation (GGA) that gives significant improvements for lattice constants, crystal structures, and metal surface energies over the most popular Perdew-Burke-Ernzerhof (PBE) GGA. The new functional is based on a diffuse radial cutoff for the exchange-hole in real space, and the analytic gradient expansion of the exchange energy for small gradients. There are no adjustable parameters, the constraining conditions of PBE are maintained, and the functional is easily implemented in existing codes.Comment: 5 pages, corrected the errors of the sublimation energy of Ih ic

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