Local-spin-density calculations for iron: Effect of spin interpolation on ground-state properties

Abstract
Scalar-relativistic self-consistent linear muffin-tin orbital (LMTO) calculations for bcc and fcc Fe have been performed with several different local approximations to the exchange and correlation energy density and potential. Overall, in contrast to the conclusions of previous studies, we find that the local-spin-density approximation to exchange and correlation can provide an adequate description of bulk Fe provided that a proper parametrization of the correlation energy density and potential of the homogeneous electron gas over both spin and density is used. Lattice constants, found from the position of the minimum of the total energy as a function of Wigner-Seitz radius, agree to within 1% (for s,p,d LMTO’s only) and within 1–2% (for s,p,d,f LMTO’s) of the experimental lattice constants for all forms used for the local correlation. The best agreement, however, was obtained using a local correlation potential derived from the Vosko-Wilk-Nusair form for the spin dependence of the correlation energy density. The calculation performed with this correlation potential was also the only calculation to correctly predict a bcc ferromagnetic ground state.