I have run the Self-Consistent Field calculations and found the highest unoccupied and lowest unoccupied energy level for Bismuth (hR2) crystal structure. Now, bismuth is a semi metal with Fermi level around 0eV. But after taking the mean/average of the highest occupied and unoccupied energy level, my Fermi energy is coming out to be around 9eV. Even in the SCF.out file, the Fermi energy is around 9eV.
I don't understand how to determine the correct Fermi energy, and why is this value coming around 9 when it should be around 0. The aim of my work is to determine the direct band gap at some specific k point (L and T). But without knowing the fermi energy, I can't find the valence and conduction band and hence I can't find the direct band gap at these k points. I have also run the band structure calculation to obtain the band structure for the higher symmetry points. For this, I used seek k-path to find the correct k-path for band structure calculation. But it's useless unless I know how to find the band gap. I don't even have the correct fermi energy I think. Please help me understand where am I going wrong.