TY - JOUR

T1 - First-principles calculation of the magnetocrystalline anisotropy energy of iron, cobalt, and nickel

AU - Daalderop, G. H.O.

AU - Kelly, P. J.

AU - Schuurmans, M. F.H.

PY - 1990/1/1

Y1 - 1990/1/1

N2 - The magnetocrystalline anisotropy energies of the elements iron, cobalt, and nickel have been calculated by means of the linear muffin-tin orbital (LMTO) method in the atomic-sphere approximation (ASA) within the framework of the local-spin-density approximation (LSDA). The so-called force theorem is used to express the total-energy difference, when spin-orbit coupling is included, as a difference in sums of Kohn-Sham single-particle eigenvalues. The results depend strongly on the location and dispersion of degenerate energy bands near the Fermi surface, and particular attention must be paid to the convergence of the Brillouin-zone integral of the single-particle eigenvalues. The calculated values of the anisotropy energy are too small by comparison with experiment, and we do not predict the correct easy axis for cobalt and nickel. We find that the variation of the anisotropy energy with changes in strain, in the magnitude of the spin-orbit coupling, for different choices of the exchange-correlation potential and for varying numbers of valence electrons are not capable of explaining these incorrect results. By comparing our calculated energy bands with those obtained by a full-potential linear augmented plane wave (FLAPW) method we conclude that the discrepancy is not attributable to terms in the potential that are neglected in the ASA.

AB - The magnetocrystalline anisotropy energies of the elements iron, cobalt, and nickel have been calculated by means of the linear muffin-tin orbital (LMTO) method in the atomic-sphere approximation (ASA) within the framework of the local-spin-density approximation (LSDA). The so-called force theorem is used to express the total-energy difference, when spin-orbit coupling is included, as a difference in sums of Kohn-Sham single-particle eigenvalues. The results depend strongly on the location and dispersion of degenerate energy bands near the Fermi surface, and particular attention must be paid to the convergence of the Brillouin-zone integral of the single-particle eigenvalues. The calculated values of the anisotropy energy are too small by comparison with experiment, and we do not predict the correct easy axis for cobalt and nickel. We find that the variation of the anisotropy energy with changes in strain, in the magnitude of the spin-orbit coupling, for different choices of the exchange-correlation potential and for varying numbers of valence electrons are not capable of explaining these incorrect results. By comparing our calculated energy bands with those obtained by a full-potential linear augmented plane wave (FLAPW) method we conclude that the discrepancy is not attributable to terms in the potential that are neglected in the ASA.

UR - http://www.scopus.com/inward/record.url?scp=0001718417&partnerID=8YFLogxK

U2 - 10.1103/PhysRevB.41.11919

DO - 10.1103/PhysRevB.41.11919

M3 - Article

AN - SCOPUS:0001718417

VL - 41

SP - 11919

EP - 11937

JO - Physical review B: Covering condensed matter and materials physics

JF - Physical review B: Covering condensed matter and materials physics

SN - 2469-9950

IS - 17

ER -