Analysis of over-magnetization of elemental transition metal solids from the SCAN density functional

• Published in: Physical review. B 2019-07, Vol.100 (4), p.041113(R), Article 041113
• Main Authors: Mejía-Rodríguez, Daniel, Trickey, S. B.
• Format: Article
• Language: English
• Subjects: Kinetic energy; Magnetization; Nickel; Switching; Transition metals
• ISSN: 2469-9950; 2469-9969
• DOI: 10.1103/PhysRevB.100.041113

Recent investigations have found that the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation exchange-correlation functional significantly over-magnetizes elemental Fe, Co, and Ni solids. For the paradigmatic case, bcc Fe, the error relative to experiment is ≳20%. A comparative analysis of magnetization results from SCAN and its deorbitalized counterpart, SCAN-L, leads to identification of the source of the discrepancy. It is not from the difference between Kohn-Sham (SCAN-L) and generalized Kohn-Sham (SCAN) procedures. The key is the iso-orbital indicator α (the ratio of the local Pauli and Thomas-Fermi kinetic energy densities). Its deorbitalized counterpart αL has more dispersion in both spin channels with respect to magnetization in an approximate region between 0.6 and 1.2 bohrs around an Fe nucleus. The overall effect is that the SCAN switching function evaluated with αL reduces the energetic disadvantage of the down channel with respect to up compared to the original α, which in turn reduces the magnetization. This identifies the cause of the SCAN magnetization error as insensitivity of the SCAN switching function to α values in the approximate range 0.5≲α≲0.8 and oversensitivity for α≳0.8.