Flat plane based double-counting free and parameter free many-body DFT+U

Aug 15, 2024
22 pages
Published in:
  • Phys.Rev.B 110 (2024) 20, 205150
  • Published: Nov 15, 2024
e-Print:
DOI:

Citations per year

0 Citations
Abstract: (APS)
The recently introduced BLOR corrective exchange-correlation functional [A. C. Burgess , Phys. Rev. B 107, L121115 (2023)] is, by construction, the unique simplified rotationally invariant DFT+U functional that enforces the flat plane condition separately on each effective orbital of a localized subspace. Detached from the Hubbard model, functionals of this type are both double-counting correction free and, when optimized in situ using appropriate error quantifiers, effectively parameter free. They are as computationally undemanding as conventional DFT+U functionals. In this work, the extension of the BLOR functional to address many-body errors (mBLOR) is derived. The mBLOR functional is built to enforce the flat plane condition on the entire subspace, rather than on each orbital individually. It depends solely on the total subspace occupancy and spin magnetization, bringing consistency with how Hubbard U and Hund J values are typically calculated, and very low complexity. In this way, interorbital errors are corrected on the same footing as the single-particle ones. Focusing on exact test cases with strong interorbital interactions, the BLOR and mBLOR functionals were benchmarked against contemporary DFT+U functionals using the total-energy extensivity condition on stretched homonuclear p-block dimers that represent various self-interaction and static-correlation error regimes, namely, singlet N2 and F2, non-spin polarized O2, and doublet Ne2+. The BLOR functional outperformed all other DFT+U functionals tested, which often act to increase total-energy errors, yet it still yielded large errors in some systems. mBLOR instead yielded low-energy errors across all four strongly correlated dimers, while being constructed using only semilocal approximation ingredients. As mBLOR would not otherwise introduce a band-gap correction in the manner that is a desirable feature of DFT+U, we developed a cost-free technique to reintroduce it automatically by moving the functional's unusual explicit derivative discontinuity into the potential. With this in place, mBLOR is the only known DFT+U functional that opens the band gap of stretched neutral homonuclear dimers without the aid of unphysical spin-symmetry breaking.
Note:
  • As accepted for publication in Phys. Rev. B. 22 pages, 9 figures, and 4 tables