ISMEAR: Difference between revisions

From VASP Wiki
No edit summary
 
(35 intermediate revisions by 8 users not shown)
Line 1: Line 1:
{{TAGDEF|ISMEAR|-5 {{!}} -4 {{!}} -3 {{!}} -2 {{!}} -1 {{!}} 0 {{!}} [integer]>0 |1}}
{{TAGDEF|ISMEAR|-15 {{!}} -14 {{!}} -5 {{!}} -4 {{!}} -3 {{!}} -2 {{!}} -1 {{!}} 0 {{!}} [integer]>0 |1}}


Description: {{TAG|ISMEAR}} determines how the partial occupancies ''f''<sub>n'''k'''</sub> are set for each orbital. {{TAG|SIGMA}} determines the width of the smearing in eV.
Description: {{TAG|ISMEAR}} determines how the partial occupancies ''f''<sub>n'''k'''</sub> are set for each orbital. {{TAG|SIGMA}} determines the width of the smearing in eV.
----
----
*{{TAG|ISMEAR}}=''N'' (''N''>0): method of Methfessel-Paxton order ''N''.
:'''Mind''': For the Methfessel-Paxton scheme the partial occupancies can be negative.
*{{TAG|ISMEAR}}=0: Gaussian smearing.
*{{TAG|ISMEAR}}=&minus;1: Fermi smearing.
*{{TAG|ISMEAR}}=&minus;2: partial occupancies are read in from the {{FILE|WAVECAR}} or {{FILE|INCAR}} file, and kept fixed throughout run.
:To set the occupancies specify
  {{TAG|FERWE}} = f1 f2 f3 ... f({{TAG|NBANDS}})
:and for spin-polarized calculations, additionally
  {{TAG|FERDO}} = f1 f2 f3 ... f({{TAG|NBANDS}})
:in the {{TAG|INCAR}} file. The (partial) occupancies must be specified for all bands and k-points. The band-index runs fastest. The occupancies must be between 0 and 1 (for spin-polarized and non-spin-polarized calculations).


:'''Mind''': Partial occupancies are also written to the {{FILE|OUTCAR}} file, but in this case they are multiplied by 2, i.e. they are between 0 and 2.
Please consider how-to guide to choose the optimal [[smearing technique]].


*{{TAG|ISMEAR}}=&minus;3: perform a loop over smearing-parameters supplied in the {{FILE|INCAR}} file.
== Tag options ==
:In this case a tag
  {{TAG|SMEARINGS}}= ismear<sub>1</sub> sigma<sub>1</sub>  ismear<sub>2</sub> sigma<sub>2</sub>  ...
:must be present in the {{FILE|INCAR}} file, supplying different smearing parameters. {{TAG|IBRION}} has to be set to -1 and {{TAG|NSW}} to the number of supplied pairs ismear<sub>i</sub>/sigma<sub>i</sub>. The first loop is done using the tetrahedron method with Blöchl corrections.


*{{TAG|ISMEAR}}=&minus;4: tetrahedron method (use a [[KPOINTS|&Gamma;-centered '''k'''-mesh]]).
*{{TAGO|ISMEAR|0|op=>}}: method of Methfessel-Paxton order {{TAG|ISMEAR}} with width {{TAG|SIGMA}}.
*{{TAG|ISMEAR}}=&minus;5: tetrahedron method with Blöchl corrections (use a [[KPOINTS|&Gamma;-centered '''k'''-mesh]]).
{{NB|mind|Methfessel-Paxton can yield erroneous results for insulators because the partial occupancies can be unphysical.|:}}


----
*{{TAGO|ISMEAR|0}}: Gaussian smearing with width {{TAG|SIGMA}}.
For the calculation of the ''total energy'' in bulk materials we recommend the tetrahedron method with Blöchl corrections ({{TAG|ISMEAR}}=-5). This method also gives a good account for the electronic density of states (DOS). The only drawback is that the method is not variational with respect to the partial occupancies. Therefore the calculated forces and the stress tensor can be wrong by up to 5 to 10% for metals.
 
For the calculation of phonon frequencies based on forces we recommend the method of Methfessel-Paxton ({{TAG|ISMEAR}}>0). For semiconductors and insulators the forces are correct, because partial occupancies do not vary and are either zero or one.
*{{TAGO|ISMEAR|-1}}: Fermi smearing with width {{TAG|SIGMA}}.
The method of Methfessel-Paxton ({{TAG|ISMEAR}}>0) also results in a very accurate description of the total energy, nevertheless the width of the smearing ({{TAG|SIGMA}}) must be chosen carefully. Too large smearing-parameters might result in a wrong total energy, small smearing parameters require a dense mesh of '''k'''-points. {{TAG|SIGMA}} should be as large as possible keeping the difference between the free energy and the total energy (i.e. the term <tt>entropy T*S</tt>) in the {{FILE|OUTCAR}} file negligible (1 meV/atom). In most cases  and  leads to very similar results. The method of Methfessel-Paxton is also the method of choice for large supercells, since the tetrahedron method is not applicable, if less than three '''k'''-points are used.
 
*{{TAGO|ISMEAR|-2}}: Partial occupancies are read in from the {{FILE|WAVECAR}} and kept fixed throughout run. Alternatively, you can also choose occupancies in the {{FILE|INCAR}} file with the tag {{TAG|FERWE}} (and {{TAG|FERDO}} for {{TAGO|ISPIN|2}} calculations).
 
*{{TAGO|ISMEAR|-3}}: perform a loop over {{TAG|SMEARINGS}} parameters supplied in the {{FILE|INCAR}} file.
 
*{{TAGO|ISMEAR|-4}}: Tetrahedron method without smearing.
 
*{{TAGO|ISMEAR|-5}}: Tetrahedron method with Blöchl corrections without smearing.
 
*{{TAGO|ISMEAR|-14}}: Tetrahedron method with Fermi-Dirac smearing {{TAG|SIGMA}}.
 
*{{TAGO|ISMEAR|-15}}: Tetrahedron method with Blöchl corrections with Fermi-Dirac smearing {{TAG|SIGMA}}.
{{NB|mind|Use a [[KPOINTS|&Gamma;-centered '''k'''-mesh]] for the tetrahedron methods.|:}}


'''Mind''': Avoid using {{TAG|ISMEAR}}>0 for semiconductors and insulators, since this often leads to incorrect results (the occupancies of some states might be larger or smaller than 1). For insulators use {{TAG|ISMEAR}}=0 or {{TAG|ISMEAR}}=-5.
== Related tags and articles ==
{{TAG|SIGMA}},
{{TAG|EFERMI}},
{{TAG|FERWE}},
{{TAG|FERDO}},
{{TAG|SMEARINGS}},
[[Smearing technique]],
[[K-point integration]]


The Gaussian smearing method also leads to reasonable results in most cases. Within this method it is necessary to extrapolate from finite {{TAG|SIGMA}} results to {{TAG|SIGMA}}=0 results. You can find an extra line in the {{FILE|OUTCAR}} file: <tt>energy( SIGMA&rarr;0 )</tt>, giving the extrapolated results. Large {{TAG|SIGMA}} values lead to a similar error as the Methfessel-Paxton scheme, but in contrast to the Methfessel-Paxton scheme one can not determine how large the error due to the smearing is with systematically reducing {{TAG|SIGMA}}. Therefore the method of Methfessel-Paxton is more convenient than the Gaussian smearing method. In addition, in the Gaussian smearing method forces and the stress tensor are consistent with the free energy and not the energy for {{TAG|SIGMA}}&rarr;0. Overall the Methfessel-Paxton method is easier to use for metallic systems.
{{sc|ISMEAR|Examples|Examples that use this tag}}


== Related Tags and Sections ==
{{TAG|SIGMA}}
----
[[The_VASP_Manual|Contents]]


[[Category:INCAR]]
[[Category:INCAR tag]][[Category:Electronic occupancy]][[Category:Electronic minimization]][[Category:Density of states]]

Latest revision as of 10:34, 27 February 2025

ISMEAR = -15 | -14 | -5 | -4 | -3 | -2 | -1 | 0 | [integer]>0
Default: ISMEAR = 1 

Description: ISMEAR determines how the partial occupancies fnk are set for each orbital. SIGMA determines the width of the smearing in eV.


Please consider how-to guide to choose the optimal smearing technique.

Tag options

Mind: Methfessel-Paxton can yield erroneous results for insulators because the partial occupancies can be unphysical.
  • ISMEAR = -2: Partial occupancies are read in from the WAVECAR and kept fixed throughout run. Alternatively, you can also choose occupancies in the INCAR file with the tag FERWE (and FERDO for ISPIN = 2 calculations).
  • ISMEAR = -4: Tetrahedron method without smearing.
  • ISMEAR = -5: Tetrahedron method with Blöchl corrections without smearing.
  • ISMEAR = -14: Tetrahedron method with Fermi-Dirac smearing SIGMA.
  • ISMEAR = -15: Tetrahedron method with Blöchl corrections with Fermi-Dirac smearing SIGMA.
Mind: Use a Γ-centered k-mesh for the tetrahedron methods.

Related tags and articles

SIGMA, EFERMI, FERWE, FERDO, SMEARINGS, Smearing technique, K-point integration

Examples that use this tag