self.energies = np.array(
            [Data.energies[band, k] for k in kpoint_indices])  # in units eV