invDFT: A CPU-GPU massively parallel tool to find exact exchange-correlation potentials from groundstate densities

Density functional theory (DFT) remains the most widely used electronic structure method. Although exact in principle, in practice, it relies on approximations to the exchange-correlation (XC) functional, which is known to be a unique functional of the electron density. Despite 50 years of active research, existing XC approximations remain far from general purpose chemical accuracy of various thermochemical and materials properties. In that light, the inverse DFT problem, of finding the exact XC potential corresponding to an accurate groundstate density, offers an insightful tool to understand the nature of the XC functional as well as aid in the development of more accurate functionals. However, solving the inverse DFT problem is fraught with several numerical challenges, such as non-uniqueness or spurious oscillations in the solution and non-convergence. We present invDFT as an open-source framework to address the outstanding challenges in inverse DFT and computed XC potentials solely from a target density. We do so by use of a systematically convergent finite-element basis and asymptotic corrections to the target density. We also employ several numerical and high-performance computing (HPC) advances that affords both efficiency and parallel scalability, on CPU-GPU hybrid architectures. We demonstrate the accuracy and scalability of invDFT using accurate full-configuration interaction (FCI) densities as well as model densities, ranging up to 100 electrons and spanning both weakly and strongly correlated molecules.