if Data.integrated_dos[lower_index][0] - Data.integrated_dos[VBM_index][0] > 0:
        energy = (Data.integrated_dos[upper_index][0] -
                  Data.integrated_dos[VBM_index][0]) * proportion