All requests for technical support from the VASP group must be addressed to:


From Vaspwiki
Jump to: navigation, search

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.
  • 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).
Mind: 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. For the calculation of phonon frequencies based on forces we recommend the method of Methfessel-Paxton (ISMEAR>0). For semiconductors and insulators the forces are correct, because partial occupancies do not vary and are either zero or one. 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 a wrong 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). 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.

Mind: Avoid using 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 ISMEAR=0 or ISMEAR=-5.

The Gaussian smearing method also leads to 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 with systematically reducing 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 SIGMA→0. Overall the Methfessel-Paxton method is easier to use for metallic systems.

  • 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.05.
  • For relaxations in metals always 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 sensible value is usually SIGMA= 0.2 (which is the default).
Mind: Avoid to use ISMEAR>0 for semiconductors and insulators, since it might cause 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


Example Calculations using this Tag

Alpha-AlF3, Alpha-SiO2, Bandgap of Si in GW, Bandgap of Si using different DFT+HF methods, Bandstructure of Si in GW (VASP2WANNIER90), Bandstructure of SrVO3 in GW, Beta-tin Si, Cd Si, Cd Si relaxation, Cd Si volume relaxation, CO, CO on Ni 111 surface, CO partial DOS, CO vibration, constrained MD using a canonical ensemble, Constrained MD using a microcanonical ensemble, Constraining the local magnetic moments, Determining the Magnetic Anisotropy, Dielectric properties of Si, Dielectric properties of Si using BSE, Dielectric properties of SiC, Equilibrium volume of Si in the RPA, Estimation of J magnetic coupling, Fcc Ni, Fcc Ni (revisited), Fcc Ni DOS, Fcc Ni DOS with hybrid functional, Fcc Si, Fcc Si bandstructure, Fcc Si DOS, Graphite TS binding energy, Graphite MBD binding energy, Graphite interlayer distance, H2O, H2O molecular dynamics, H2O vibration, Including the Spin-Orbit Coupling, liquid Si, MgO optimum mixing, Model BSE calculation on Si, Ni 100 surface bandstructure, Ni 100 surface DOS, Ni 100 surface relaxation, Ni 111 surface high precision, NiO, NiO GGA+U, NiO HSE06, NiO LSDA, NiO LSDA+U, O atom, O atom spinpolarized, O atom spinpolarized low symmetry, O dimer, Partial DOS of CO on Ni 111 surface, Si bandstructure, Si HSE bandstructure, Standard relaxation, STM of graphene, STM of graphite, Vibrational frequencies of CO on Ni 111 surface