Abstract
We create accurate exchange-correlation (XC) functionals for density functional theory using neural network models and grid-based density descriptors. We investigate the impact of imposing exact constraints based on established physical laws to improve results, using a physics-informed approach. Additionally, we explore how neural networks can be trained to emulate existing XC functionals, such as PBE, by utilizing both the functionals and their derivatives. To facilitate this work, we have developed a new machine learning framework called xcquinox. This framework employs neural network architectures implemented in JAX, a library for automatically differentiable mathematical operations, allowing us to leverage PySCF-AD, an extension of the PySCF package that incorporates automatic differentiation capabilities. We emphasize the importance of incorporating either exact potentials or exact densities in the training process, as training an energy functional without information about functional derivatives—essential for determining the XC potential—can lead to models that predict correct energies but yield incorrect densities.
*The authors would like to thank Stony Brook Research Computing and Cyberinfrastructure and the Institute for Advanced Computational Science at Stony Brook University for access to the high-performance SeaWulf and Ookami computing systems. SeaWulf was made possible by $1.85M in grants from the National Science Foundation (Awards No. 1531492 and No. 2215987) and matching funds from the Empire State Development's Division of Science, Technology and Innovation (NYSTAR) program (Contract C210148), and Ookami through a $5M National Science Foundation Grant (No. 1927880). This work was funded in part by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, under Award No. DE-SC0019394, as part of the CCS Program.