The electronic charge density ρ(r) is a fundamental observable that uniquely defines molecules and materials with non-degenerate electronic ground states . It plays a central role in density-functional approximations (DFAs), which are widely used in computational chemistry and materials science. Existing DFAs often produce reasonable approximations to the ground state energy, which makes them so useful in practice. However, the accuracy of DFAs in predicting energies has not necessarily translated into better predictions of ρ(r). This trade-off was highlighted in the seminal work by Lyssenko et al. , revealing discrepancies between DFA-generated densities and near-exact ab initio references. Foundational works by Burke and co-authors on density-corrected DFAs found that evaluating DFAs on the Hartree-Fock (HF) density without self-consistency yields more accurate energies for important prototypical examples including e.g. stretched heteronuclear bonds or torsional barriers in flexible molecules . Hence, accurate models for electron densities would benefit the developments of next-generation DFAs as well as machine-learned surrogate models for force fields and other properties of atomistic systems throughout physics, chemistry, and biology.