elif all(item > self.fermi_energy + 30 for item in energies[index]):
                energies = np.delete(energies, index, axis=0)
                occupancy = np.delete(occupancy, index, axis=0)
                index_count = len(occupancy)