Gap edge eigenpairs from density matrix purification using moments of the Dirac distribution

In this work, we propose a simple method to resolve the eigenstates located at the band gap edges of an electronic eigenspectrum using only the quasi-purified one-particle density matrix as input. The theoretical framework relies on the decomposition of the occupation number variance into a particle