local density matrices must be stored with charge density for subsequent nscf calculations