if states_per_unit_cell < np.absolute(
                Data.integrated_dos[i][1] - Data.integrated_dos[CBM_index][1]):