🤖 AI Summary
This work addresses the longstanding bottleneck in orbital-free density functional theory (OF-DFT): the absence of energy functionals that simultaneously ensure numerical stability and chemical accuracy. We propose a machine-learned, SE(3)-equivariant graph neural network (GNN) functional that directly maps electron density to total molecular energy—bypassing Kohn–Sham equation solving entirely. To overcome scarcity of physically meaningful training data, we introduce a novel active augmentation strategy based on perturbed external potentials, enabling generation of diverse, physically consistent density configurations. Evaluated on the full QM9 dataset, our method achieves a mean absolute error of 1.2 kcal/mol—meeting the chemical accuracy threshold—while attaining >98% electron density convergence. Crucially, computational cost is reduced by approximately two orders of magnitude relative to standard KS-DFT. This represents the first demonstration of OF-DFT achieving unified high accuracy, robustness, and efficiency for large-scale molecular systems.
📝 Abstract
Hohenberg and Kohn have proven that the electronic energy and the one-particle electron density can, in principle, be obtained by minimizing an energy functional with respect to the density. Given that decades of theoretical work have so far failed to produce this elusive exact energy functional promising great computational savings, it is reasonable to try and learn it empirically. Using rotationally equivariant atomistic machine learning, we obtain for the first time a density functional that, when applied to the organic molecules in QM9, yields energies with chemical accuracy while also converging to meaningful electron densities. Augmenting the training data with densities obtained from perturbed potentials proved key to these advances. Altogether, we are now closer than ever to fulfilling Hohenberg and Kohn's promise, paving the way for more efficient calculations in large molecular systems.