The bands_above_fermi_level argument to the DiagonalizationSolver gives you the option to include only a limited number of unoccupied bands when diagonalizing the Hamiltonian. This can give both speed-up and reduce memory consumption. However, if you include too few bands (i.e. exclude bands that are not completely empty), you will not correctly diagonalize the Hamiltonian, and your calculations may be inaccurate, hard to converge, or both.