Calculating the Fermi energy in a metal
Posted: Fri Apr 26, 2024 2:57 pm
VASP team,
I would like bit more clarification as to how VASP calculates the Fermi energy for a metal. I understand that it follows the following algorithm as described in the fermi_energy.F file.
If possible, would you be able to explain this method in a bit more detail?
Would the Fermi energy as a result of this algorithm get you the same Fermi energy found by integrating the DOS*occupations over energy constrained by the number of electrons?
Thanks for your time!
Brian
I would like bit more clarification as to how VASP calculates the Fermi energy for a metal. I understand that it follows the following algorithm as described in the fermi_energy.F file.
Code: Select all
The basic algorithm to determine the Fermi energy works as follows:
* We consider a single k-point and evaluate the number of occupied states
per spin channel there. Here, we need to account for the fact that the
states may be split over multiple CPUs. We use a bisection algorithm with
step function smearing to do so.
* For metals or when the smearing lead to an incorrect number of electrons,
we determine the Fermi energy with a bisection method between the band
extrema. We compute the number of electrons for a Fermi energy in the
middle of the interval. When the result is smaller or larger than the
expected number of electrons, we replace the lower or upper interval
boundary, respectively.
Would the Fermi energy as a result of this algorithm get you the same Fermi energy found by integrating the DOS*occupations over energy constrained by the number of electrons?
Thanks for your time!
Brian