Many recent density functionals are numerically ill-behaved.

Most computational studies in chemistry and materials science are based on the use of density functional theory. Although the exact density functional is unknown, several density functional approximations (DFAs) offer a good balance of affordable computational cost and semi-quantitative accuracy for applications. The development of DFAs still continues on many fronts, and several new DFAs aiming for improved accuracy are published every year. However, the numerical behavior of these DFAs is an often-overlooked problem. In this work, we look at all 592 DFAs for three-dimensional systems available in Libxc 5.2.2 and examine the convergence of the density functional total energy based on tabulated atomic Hartree-Fock wave functions. We show that several recent DFAs, including the celebrated SCAN family of functionals, show impractically slow convergence with typically used numerical quadrature schemes, making these functionals unsuitable both for routine applications and high-precision studies, as thousands of radial quadrature points may be required to achieve sub-μE accurate total energies for these functionals, while standard quadrature grids like the SG-3 grid only contain O(100) radial quadrature points. These results are both a warning to users to always check the sufficiency of the quadrature grid when adopting novel functionals, as well as a guideline to the theory community to develop better-behaved density functionals.