EFERMI: Difference between revisions

From VASP Wiki
m (Schlipf moved page Vaspwiki:EFERMI to EFERMI)
No edit summary
 
(9 intermediate revisions by 2 users not shown)
Line 1: Line 1:
{{TAGDEF|EFERMI|MIDGAP {{!}} LEGACY {{!}} [real]|MIDGAP}}
{{TAGDEF|EFERMI|MIDGAP {{!}} LEGACY {{!}} [real]|LEGACY}}


Description: Defines how the Fermi energy is calculated in VASP.
Description: Defines how the Fermi energy is calculated in VASP. It is recommended to use EFERMI = MIDGAP.


----
----
== Fermi energy in semiconductors ==


For semiconducting materials, the Fermi energy is not uniquely defined in the bandgap of the material.
For semiconducting materials, the Fermi energy is not uniquely defined in the bandgap of the material.
Any value that produces the correct number of electrons is allowed.
Any value that produces the correct number of electrons ({{TAG|NELECT}}) is allowed.
By default, VASP will put the Fermi energy in the middle of the gap because this is the most consistent with increasing the smearing {{TAG|SIGMA}}.
By default, VASP places the Fermi energy at a somewhat arbitrary value within the bandgap.
The precise value depends on values chosen for the smearing ({{TAG|ISMEAR}} and {{TAG|SIGMA}}) and the density of states ({{TAG|EMIN}}, {{TAG|EMAX}}, and {{TAG|NEDOS}}).
Typically, this places the Fermi energy towards the bottom of the bandgap.
 
You can change this behavior by setting {{TAG|EFERMI}} = MIDGAP (recommended).
VASP will then put the Fermi energy in the middle of the gap because this is the most consistent with increasing the smearing {{TAG|SIGMA}}.
This algorithm to determine the Fermi energy was introduced in VASP.6.4.
This algorithm to determine the Fermi energy was introduced in VASP.6.4.
The value of the Fermi energy should not affect the outcome of the calculation.


For prior versions, VASP places the Fermi energy at a somewhat arbitrary value within the bandgap.
== Fermi energy in metals ==
The precise value depends on values chosen for the smearing ({{TAG|ISMEAR}} and {{TAG|SIGMA}}) and the density of states ({{TAG|EMIN}}, {{TAG|EMAX}}, and {{TAG|NEDOS}}).
Typically this places the Fermi energy towards the bottom of the bandgap.
The value of the Fermi energy should not affect the outcome of the calculation but you can recover the old behavior set {{TAG|EFERMI}} = LEGACY.
 
Enforce a particular value for the Fermi energy by setting it to a numeric value.
Note that this may change the number of electrons.


== Related Tags and Sections ==
MIDGAP and LEGACY should yield the same value when your system does not have a gap.
The Fermi energy is placed precisely at the value so that underneath are enough states to accommodate {{TAG|NELECT}} electrons. The evaluation of the Fermi energy involves an integral over the 1. BZ.
Therefore, the Fermi energy in metals needs to be converged with respect to the {{FILE|KPOINTS}} mesh and smearing ({{TAG|ISMEAR}}, {{TAG|SIGMA}}).
{{NB|tip| If you are interested in the properties at the Fermi energy (e.g., for transport calculations), you should compute the Fermi energy with a very dense '''k'''-point mesh. To save computational time, you can fix the charge density ({{TAG|ICHARG}} {{=}} 11) once the Kohn-Sham states have converged with respect to the '''k'''-point density and increase the number of '''k'''-point further to converge the value of the Fermi energy.}}
== Fixed Fermi energy ==
Occasionally, you want to compute systems with fixed Fermi energy for a given charge density.
To this end, set {{TAG|EFERMI}} to a numeric value and {{TAG|ICHARG}} = 11.
A possible use case is to set {{TAG|EFERMI}} to the converged Fermi energy in a band-structure calculation.
You may use this to introduce electron doping/depletion to a system.
{{NB|warning|The Fermi energy cannot be computed based on '''k''' points along a path.}}
This is, in particular, important for band-structure calculation of metallic systems because for gaped systems the Fermi energy often still ends up at a valid value within the gap despite the inaccurate computation. Hence, to plot a band structure the band energies should be taken from the calculation with fixed charge density based on '''k'''-points along a path, but the Fermi energy should be taken from the calculation based on a '''k''' mesh (e.g., the scf calculation by which the charge density was obtained or a more precisely converged Fermi energy based on the same fixed charge density). {{TAG|EFERMI}} can fix the Fermi energy to the proper value during the band-structure calculation.
== Related tags and articles ==
{{TAG|ISMEAR}},
{{TAG|ISMEAR}},
{{TAG|SIGMA}},
{{TAG|SIGMA}},
Line 25: Line 38:
{{TAG|NEDOS}}
{{TAG|NEDOS}}


----
[[Category:INCAR tag]][[Category:Electronic ground-state properties]]
<!-- Activate this when you want to add it to the list of INCAR tags
[[Category:INCAR]][[Category:Electronic Minimization]][[Category:Density of States]]
-->

Latest revision as of 09:31, 17 October 2024

EFERMI = MIDGAP | LEGACY | [real]
Default: EFERMI = LEGACY 

Description: Defines how the Fermi energy is calculated in VASP. It is recommended to use EFERMI = MIDGAP.


Fermi energy in semiconductors

For semiconducting materials, the Fermi energy is not uniquely defined in the bandgap of the material. Any value that produces the correct number of electrons (NELECT) is allowed. By default, VASP places the Fermi energy at a somewhat arbitrary value within the bandgap. The precise value depends on values chosen for the smearing (ISMEAR and SIGMA) and the density of states (EMIN, EMAX, and NEDOS). Typically, this places the Fermi energy towards the bottom of the bandgap.

You can change this behavior by setting EFERMI = MIDGAP (recommended). VASP will then put the Fermi energy in the middle of the gap because this is the most consistent with increasing the smearing SIGMA. This algorithm to determine the Fermi energy was introduced in VASP.6.4. The value of the Fermi energy should not affect the outcome of the calculation.

Fermi energy in metals

MIDGAP and LEGACY should yield the same value when your system does not have a gap. The Fermi energy is placed precisely at the value so that underneath are enough states to accommodate NELECT electrons. The evaluation of the Fermi energy involves an integral over the 1. BZ. Therefore, the Fermi energy in metals needs to be converged with respect to the KPOINTS mesh and smearing (ISMEAR, SIGMA).

Tip: If you are interested in the properties at the Fermi energy (e.g., for transport calculations), you should compute the Fermi energy with a very dense k-point mesh. To save computational time, you can fix the charge density (ICHARG = 11) once the Kohn-Sham states have converged with respect to the k-point density and increase the number of k-point further to converge the value of the Fermi energy.

Fixed Fermi energy

Occasionally, you want to compute systems with fixed Fermi energy for a given charge density. To this end, set EFERMI to a numeric value and ICHARG = 11. A possible use case is to set EFERMI to the converged Fermi energy in a band-structure calculation. You may use this to introduce electron doping/depletion to a system.

Warning: The Fermi energy cannot be computed based on k points along a path.

This is, in particular, important for band-structure calculation of metallic systems because for gaped systems the Fermi energy often still ends up at a valid value within the gap despite the inaccurate computation. Hence, to plot a band structure the band energies should be taken from the calculation with fixed charge density based on k-points along a path, but the Fermi energy should be taken from the calculation based on a k mesh (e.g., the scf calculation by which the charge density was obtained or a more precisely converged Fermi energy based on the same fixed charge density). EFERMI can fix the Fermi energy to the proper value during the band-structure calculation.

Related tags and articles

ISMEAR, SIGMA, EMIN, EMAX, NEDOS