ISMEAR: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
(15 intermediate revisions by 3 users not shown)
Line 4: Line 4:
----
----
*{{TAG|ISMEAR}}=''N'' (''N''>0): method of Methfessel-Paxton order ''N''.
*{{TAG|ISMEAR}}=''N'' (''N''>0): method of Methfessel-Paxton order ''N''.
:'''Mind''': For the Methfessel-Paxton scheme the partial occupancies can be negative.
:'''Mind''': For the Methfessel-Paxton scheme the partial occupancies can be negative, as well as larger than 1. This can yield erroneous results for insulators.
*{{TAG|ISMEAR}}=0: Gaussian smearing.
*{{TAG|ISMEAR}}=0: Gaussian smearing.
*{{TAG|ISMEAR}}=−1: Fermi smearing.
*{{TAG|ISMEAR}}=−1: Fermi smearing.
*{{TAG|ISMEAR}}=−2: partial occupancies are read in from the {{FILE|WAVECAR}} or {{FILE|INCAR}} file, and kept fixed throughout run.
*{{TAG|ISMEAR}}=−2: partial occupancies are read in from the {{FILE|WAVECAR}} or {{FILE|INCAR}} file, and kept fixed throughout run.
:To set the occupancies specify
:To set the occupancies, specify
   {{TAG|FERWE}} = f1 f2 f3 ... f({{TAG|NBANDS}})
   {{TAG|FERWE}} = f(1) f(2) f(3) ... f({{TAG|NBANDS}}&times;N<sub>'''k'''</sub>)
:and for spin-polarized calculations, additionally
:and for spin-polarized calculations, additionally
   {{TAG|FERDO}} = f1 f2 f3 ... f({{TAG|NBANDS}})
   {{TAG|FERDO}} = f(1) f(2) f(3) ... f({{TAG|NBANDS}}&times;N<sub>'''k'''</sub>)
: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).  
: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.
:Note that the 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.


*{{TAG|ISMEAR}}=&minus;3: perform a loop over smearing-parameters supplied in the {{FILE|INCAR}} file.
*{{TAG|ISMEAR}}=&minus;3: perform a loop over smearing-parameters supplied in the {{FILE|INCAR}} file.
Line 21: Line 21:
: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.
: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:
*{{TAG|ISMEAR}}=&minus;4: tetrahedron method (use a [[KPOINTS|&Gamma;-centered '''k'''-mesh]]).
*{{TAG|ISMEAR}}=&minus;5:
*{{TAG|ISMEAR}}=&minus;5: tetrahedron method with Blöchl corrections (use a [[KPOINTS|&Gamma;-centered '''k'''-mesh]]).
 
----
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. Only for semiconductors and insulators the forces are correct, because the partial occupancies do not vary and are either zero or one. For the calculation of forces and phonon frequencies in metals, we recommend the method of Methfessel-Paxton ({{TAG|ISMEAR}}>0).
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 an incorrect 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).
 
'''A serious warning''': Avoid using {{TAG|ISMEAR}}>0 for semiconductors and insulators, since this often leads to incorrect results (the occupancies of some states might be smaller than 0, or larger than 1). Errors for e.g. phonons frequencies can be substantial exceeding 20 %. These errors are very hard to spot if you do not look carefully. For insulators use {{TAG|ISMEAR}}=0 or {{TAG|ISMEAR}}=-5.
 
The Gaussian smearing method leads to very 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 without systematically reducing {{TAG|SIGMA}}. In this respect, 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 somewhat easier to use for metallic systems.
 
'''Summary''':
 
*If you have no a priori knowledge of your system, for instance, if you do not know whether your system is an insulator, semiconductor or metal than always use Gaussian smearing {{TAG|ISMEAR}}=0 in combination with a small {{TAG|SIGMA}}=0.03-0.05.
This is not the default in VASP yet, so to be on the safe side, you might want to include this setting in all your {{TAG|INCAR}} files.
 
*For semiconductors or insulators, use the tetrahedron method ({{TAG|ISMEAR}}=-5), if the cell is too large (or if you use only a single or two '''k'''-points) use {{TAG|ISMEAR}}=0 in combination with a small {{TAG|SIGMA}}=0.03-0.05.
 
*For relaxations ''in metals'' use {{TAG|ISMEAR}}=1 or {{TAG|ISMEAR}}=2 and an appropriate {{TAG|SIGMA}} value (the entropy term should be less than 1 meV per atom). For metals a reasonable value is often SIGMA= 0.2 (which is the default).
:'''Mind again''': Avoid to use {{TAG|ISMEAR}}>0 for semiconductors and insulators, since it might cause severe problems.
 
*For the calculations of the DOS and very accurate ''total energy'' calculations (no relaxation in metals) use the tetrahedron method ({{TAG|ISMEAR}}=-5).


== Related Tags and Sections ==
== Related Tags and Sections ==
{{TAG|SIGMA}}
{{TAG|SIGMA}},
{{TAG|FERWE}},
{{TAG|FERDO}},
{{TAG|SMEARINGS}}
 
{{sc|ISMEAR|Examples|Examples that use this tag}}
----
----
[[The_VASP_Manual|Contents]]
[[The_VASP_Manual|Contents]]


[[Category:INCAR]]
[[Category:INCAR]][[Category:Electronic Minimization]][[Category:Electronic Minimization Methods]][[Category:Density of States]]

Revision as of 06:56, 30 April 2020

ISMEAR = -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.


  • ISMEAR=N (N>0): method of Methfessel-Paxton order N.
Mind: For the Methfessel-Paxton scheme the partial occupancies can be negative, as well as larger than 1. This can yield erroneous results for insulators.
  • ISMEAR=0: Gaussian smearing.
  • ISMEAR=−1: Fermi smearing.
  • ISMEAR=−2: partial occupancies are read in from the WAVECAR or INCAR file, and kept fixed throughout run.
To set the occupancies, specify
 FERWE = f(1) f(2) f(3) ... f(NBANDS×Nk)
and for spin-polarized calculations, additionally
 FERDO = f(1) f(2) f(3) ... f(NBANDS×Nk)
in the 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).
Note that the partial occupancies are also written to the OUTCAR file, but in this case they are multiplied by 2, i.e. they are between 0 and 2.
  • ISMEAR=−3: perform a loop over smearing-parameters supplied in the INCAR file.
In this case a tag
 SMEARINGS= ismear1 sigma1  ismear2 sigma2  ...
must be present in the INCAR file, supplying different smearing parameters. IBRION has to be set to -1 and NSW to the number of supplied pairs ismeari/sigmai. The first loop is done using the tetrahedron method with Blöchl corrections.

For the calculation of the total energy in bulk materials we recommend the tetrahedron method with Blöchl corrections (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. Only for semiconductors and insulators the forces are correct, because the partial occupancies do not vary and are either zero or one. For the calculation of forces and phonon frequencies in metals, we recommend the method of Methfessel-Paxton (ISMEAR>0). The method of Methfessel-Paxton (ISMEAR>0) also results in a very accurate description of the total energy, nevertheless the width of the smearing (SIGMA) must be chosen carefully. Too large smearing-parameters might result in an incorrect total energy, small smearing parameters require a dense mesh of k-points. SIGMA should be as large as possible keeping the difference between the free energy and the total energy (i.e. the term entropy T*S) in the OUTCAR file negligible (1 meV/atom).

A serious warning: Avoid using ISMEAR>0 for semiconductors and insulators, since this often leads to incorrect results (the occupancies of some states might be smaller than 0, or larger than 1). Errors for e.g. phonons frequencies can be substantial exceeding 20 %. These errors are very hard to spot if you do not look carefully. For insulators use ISMEAR=0 or ISMEAR=-5.

The Gaussian smearing method leads to very reasonable results in most cases. Within this method it is necessary to extrapolate from finite SIGMA results to SIGMA=0 results. You can find an extra line in the OUTCAR file: energy( SIGMA→0 ), giving the extrapolated results. Large 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 without systematically reducing SIGMA. In this respect, 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 SIGMA→0. Overall the Methfessel-Paxton method is somewhat easier to use for metallic systems.

Summary:

  • If you have no a priori knowledge of your system, for instance, if you do not know whether your system is an insulator, semiconductor or metal than always use Gaussian smearing ISMEAR=0 in combination with a small SIGMA=0.03-0.05.

This is not the default in VASP yet, so to be on the safe side, you might want to include this setting in all your INCAR files.

  • For semiconductors or insulators, use the tetrahedron method (ISMEAR=-5), if the cell is too large (or if you use only a single or two k-points) use ISMEAR=0 in combination with a small SIGMA=0.03-0.05.
  • For relaxations in metals use ISMEAR=1 or ISMEAR=2 and an appropriate SIGMA value (the entropy term should be less than 1 meV per atom). For metals a reasonable value is often SIGMA= 0.2 (which is the default).
Mind again: Avoid to use ISMEAR>0 for semiconductors and insulators, since it might cause severe problems.
  • For the calculations of the DOS and very accurate total energy calculations (no relaxation in metals) use the tetrahedron method (ISMEAR=-5).

Related Tags and Sections

SIGMA, FERWE, FERDO, SMEARINGS

Examples that use this tag


Contents