Typically, electron densities from quantum calcula-tions are constructed from the wave functions of the occupied states. This so-called one-particle density matrix, however, grows as the square of the system size. To keep the size of the density representation linear, we chose to express the density in terms of an atom-centered, auxiliary, "density fitting" basis: