Reproducibility of density functional approximations: how new functionals should be reported

Abstract

Density functional theory is the workhorse of chemistry and materials science, and novel density functional approximations (DFAs) are published every year. To become available in program packages, the novel DFAs need to be (re)implemented. However, according to our experience as developers of Libxc [Lehtola et al, SoftwareX 7, 1 (2018)], a constant problem in this task is verification, due to the lack of reliable reference data. As we discuss in this work, this lack has lead to several non-equivalent implementations of functionals such as BP86, PW91, PBE, and B3LYP across various program packages, yielding different total energies. Through careful verification, we have also found many issues with incorrect functional forms in recent DFAs. The goal of this work is to ensure the reproducibility of DFAs: DFAs must be verifiable in order to prevent reappearances of the abovementioned errors and incompatibilities. A common framework for verification and testing is therefore needed. We suggest several ways in which reference energies can be produced with free and open source software, either with non-self-consistent calculations with tabulated atomic densities or via self-consistent calculations with various program packages. The employed numerical parameters -- especially, the quadrature grid -- need to be converged to guarantee a $\lesssim 0.1\mu E_h$ precision in the total energy, which is nowadays routinely achievable in fully numerical calculations. Moreover, as such sub-$\mu E_h$ level agreement can only be achieved when fully equivalent implementations of the DFA are used, also the source code of the reference implementation should be made available in any publication describing a new DFA.

arXiv comment16 pages, 1 figure

References in corpus (27)

• Climbing the Density Functional Ladder: Non-Empirical Meta-Generalized Gradient Approximation Designed for Molecules and Solids
• Quantum ESPRESSO toward the exascale
• Development of Novel Density Functionals for Thermochemical Kinetics
• Recent developments in the PySCF program package
• Long-Range Corrected Hybrid Density Functionals with Improved Dispersion Corrections
• Regularized SCAN functional
• The Role of the Basis Set: Assessing Density Functional Theory
• On the integration accuracy in molecular density functional theory calculations using Gaussian basis sets
• The Elephant in the Room of Density Functional Theory Calculations
• An overview of self-consistent field calculations within finite basis sets
• A review on non-relativistic fully numerical electronic structure calculations on atoms and diatomic molecules
• A self-interaction-free local hybrid functional: Accurate binding energies vis-à-vis accurate ionization potentials from Kohn-Sham eigenvalues
• Many recent density functionals are numerically ill-behaved
• Fully numerical calculations on atoms with fractional occupations. Range-separated exchange functionals
• Evolving symbolic density functionals
• Fully numerical Hartree-Fock and density functional calculations. I. Atoms
• Meta-GGA density functional calculations on atoms with spherically symmetric densities in the finite element formalism
• Communication: Curing basis set overcompleteness with pivoted Cholesky decompositions
• Improving the performance of Tao-Mo non-empirical density functional with broader applicability in quantum chemistry and material sciences
• Atomic electronic structure calculations with Hermite interpolating polynomials
• Fully numerical Hartree-Fock and density functional calculations. II. Diatomic molecules
• Polarized Gaussian basis sets from one-electron ions
• Generalized gradient approximation correlation energy functionals based on the uniform electron gas with gap model
• The MRChem multiresolution analysis code for molecular electronic structure calculations: performance and scaling properties
• Beyond the local approximation to exchange and correlation: the role of the Laplacian of the density in the energy density of Si
• Improving the applicability of the Pauli kinetic energy density based semilocal functional for solids
• Meta-local density functionals: a new rung on Jacob's ladder

Cited by in corpus (7)

• JARVIS-Leaderboard: A Large Scale Benchmark of Materials Design Methods
• A call to arms: making the case for more reusable libraries
• Predicting The One-Particle Density Matrix With Machine Learning
• Density functional benchmark for quadruple hydrogen bonds
• Ensemble Generalization of the Perdew-Zunger Self-Interaction Correction: a Way Out of Multiple Minima and Symmetry Breaking
• Review of the finite difference Hartree-Fock method for atoms and diatomic molecules, and its implementation in the x2dhf program
• libwignernj: a reusable C/C++/Fortran/Python library for exact Wigner symbols and related coefficients