proportion = (states_per_unit_cell -
                  np.absolute(Data.integrated_dos[lower_index][1] -
                              Data.integrated_dos[VBM_index][1])
                  ) / np.absolute(Data.integrated_dos[upper_index][1] -
                                  (Data.integrated_dos[lower_index][1]))