PREC: Difference between revisions

From VASP Wiki
No edit summary
(Added description of SingleN)
(15 intermediate revisions by 3 users not shown)
Line 3: Line 3:


Description: {{TAG|PREC}} specifies the "precision"-mode.
Description: {{TAG|PREC}} specifies the "precision"-mode.
----
We recommend to use {{TAG|PREC}}=Normal or {{TAG|PREC}}=Accurate.
The settings {{TAG|PREC}}=Medium and {{TAG|PREC}}=High (available already in VASP.4.4)
are no longer recommended. {{TAG|PREC}}=Normal can be used for most
routine calculations, whereas {{TAG|PREC}}=Accurate reduces egg-box effects
and errors even further.  Specifically,
{{TAG|PREC}}=Accurate strictly avoids any aliasing or wrap around errors and uses an augmentation grid ({{TAG|NGXF}}×{{TAG|NGYF}}×{{TAG|NGZF}}) that is exactly twice as large as the coarser grid ({{TAG|NGX}}×{{TAG|NGY}}×{{TAG|NGZ}}) used for the representation of the pseudo orbitals. {{TAG|PREC}}=Accurate increases the memory requirements somewhat, but it should be used if very accurate forces (phonons and second derivatives) are required. The accuracy of forces can be sometimes further improved by specifying {{TAG|ADDGRID}}=.TRUE. (but user reports are somewhat contradictory, whether
this really helps).
The use of {{TAG|PREC}}=High is no longer recommend (and is  only available for compatibility). Essentially {{TAG|PREC}}=High only increases the energy cutoff
by 30 %.  The same effect can be achieved by just manually increasing {{TAG|ENCUT}}
in the INCAR file. If, additionally, very accurate forces are required, {{TAG|PREC}}=Accurate can be used in combination with an increased energy cutoff.
'''Mind''': we recommend to specify the energy cutoff ({{TAG|ENCUT}}) always manually in the {{FILE|INCAR}} file, to avoid incompatibilities between calculations.
----
----
The {{TAG|PREC}} parameter sets the defaults for four sets of parameters ({{TAG|ENCUT}}; {{TAG|NGX}}, {{TAG|NGY}}, {{TAG|NGZ}}; {{TAG|NGXF}}, {{TAG|NGYF}}, {{TAG|NGZF}}, and {{TAG|ROPT}}):
The {{TAG|PREC}} parameter sets the defaults for four sets of parameters ({{TAG|ENCUT}}; {{TAG|NGX}}, {{TAG|NGY}}, {{TAG|NGZ}}; {{TAG|NGXF}}, {{TAG|NGYF}}, {{TAG|NGZF}}, and {{TAG|ROPT}}):
Line 11: Line 25:
| Normal ||align="center"|max({{TAG|ENMAX}}) || 3/2&times;G<sub>cut</sub> ||align="center"|2&times;{{TAG|NGX}} ||-5&times;10<sup>-4</sup>
| Normal ||align="center"|max({{TAG|ENMAX}}) || 3/2&times;G<sub>cut</sub> ||align="center"|2&times;{{TAG|NGX}} ||-5&times;10<sup>-4</sup>
|-
|-
| Single ||align="center"|max({{TAG|ENMAX}}) || 3/2&times;G<sub>cut</sub> ||align="center"|{{TAG|NGX}} ||-5&times;10<sup>-4</sup>
| Single (VASP.5) ||align="center"|max({{TAG|ENMAX}}) || 3/2&times;G<sub>cut</sub> ||align="center"|{{TAG|NGX}} ||-5&times;10<sup>-4</sup>
|-
| Single (VASP.6) ||align="center"|max({{TAG|ENMAX}}) || 2 &times;G<sub>cut</sub> ||align="center"|{{TAG|NGX}} ||-5&times;10<sup>-4</sup>
|-
|SingleN (VASP.6 development)
|max({{TAG|ENMAX}})
|3/2&times;G<sub>cut</sub>
|{{TAG|NGX}}
| -5&times;10<sup>-4</sup>
|-
|-
| Accurate ||align="center"|max({{TAG|ENMAX}}) ||align="center"|2&times;G<sub>cut</sub> ||align="center"|2&times;{{TAG|NGX}}||-2.5&times;10<sup>-4</sup>
| Accurate ||align="center"|max({{TAG|ENMAX}}) ||align="center"|2&times;G<sub>cut</sub> ||align="center"|2&times;{{TAG|NGX}}||-2.5&times;10<sup>-4</sup>
|-
| '''old settings:'''
|-
|-
| Low ||align="center"|max({{TAG|ENMIN}}) || 3/2&times;G<sub>cut</sub> ||align="center"|3&times;G<sub>aug</sub> ||-1&times;10<sup>-2</sup>
| Low ||align="center"|max({{TAG|ENMIN}}) || 3/2&times;G<sub>cut</sub> ||align="center"|3&times;G<sub>aug</sub> ||-1&times;10<sup>-2</sup>
Line 25: Line 49:
with E<sub>cut</sub>={{TAG|ENCUT}} and E<sub>aug</sub>={{TAG|ENAUG}}.
with E<sub>cut</sub>={{TAG|ENCUT}} and E<sub>aug</sub>={{TAG|ENAUG}}.


'''Mind''': the settings {{TAG|PREC}}=Normal and {{TAG|PREC}}=Accurate are only available in VASP.4.5 and newer versions. The setting {{TAG|PREC}}=Single is only available as of VASP.5.1.
'''Mind''': the settiLow ||align="center"|max({{TAG|ENMIN}}) || 3/2&times;G<sub>cut</sub> ||align="center"|3&times;G<sub>aug</sub> ||-1&times;10<sup>-2</sup>ngs {{TAG|PREC}}=Normal and {{TAG|PREC}}=Accurate are only available in VASP.4.5 and newer versions. The setting {{TAG|PREC}}=Single is only available as of VASP.5.1, and has been updated in VASP.6.


*Plane wave kinetic energy cutoff ({{TAG|ENCUT}}):
*Plane wave kinetic energy cutoff ({{TAG|ENCUT}}):
:The {{TAG|PREC}}-tag determines the energy cutoff {{TAG|ENCUT}}, if {{TAG|ENCUT}} is not specified in the {{FIL|INCAR}} file. For {{TAG|PREC}}=Low, {{TAG|ENCUT}} will be set to the maximal {{TAG|ENMIN}} value found on the {{FILE|POTCAR}} file. For {{TAG|PREC}}=Medium and {{TAG|PREC}}=Accurate, {{TAG|ENCUT}} will be set to maximal {{TAG|ENMAX}} value found on the {{FILE|POTCAR}} file. Finally for {{TAG|PREC}}=High, {{TAG|ENCUT}} is set to the maximal {{TAG|ENMAX}} value on the {{FILE|POTCAR}} file plus 30%. {{TAG|PREC}}=High guarantees that the absolute energies are converged to a few meV, and it ensures that the stress tensor is converged within a few kBar. In general, an increased energy cutoff is only required for accurate evaluation of quantities related to the stress tensor (e.g. elastic properties).
:The {{TAG|PREC}}-tag determines the energy cutoff {{TAG|ENCUT}}, if {{TAG|ENCUT}} is not specified in the {{FILE|INCAR}} file. For {{TAG|PREC}}=Low, {{TAG|ENCUT}} will be set to the maximal {{TAG|ENMIN}} value found on the {{FILE|POTCAR}} file. For {{TAG|PREC}}=Medium and {{TAG|PREC}}=Accurate, {{TAG|ENCUT}} will be set to maximal {{TAG|ENMAX}} value found on the {{FILE|POTCAR}} file. Finally for {{TAG|PREC}}=High, {{TAG|ENCUT}} is set to the maximal {{TAG|ENMAX}} value on the {{FILE|POTCAR}} file plus 30%. {{TAG|PREC}}=High guarantees that the absolute energies are converged to a few meV, and it ensures that the stress tensor is converged within a few kBar. In general, an increased energy cutoff is only required for accurate evaluation of quantities related to the stress tensor (e.g. elastic properties).


*FFT grids ({{TAG|NGX}},{{TAG|NGY}},{{TAG|NGZ}},{{TAG|NGXF}},{{TAG|NGYF}},{{TAG|NGZF}}):
*FFT grids ({{TAG|NGX}},{{TAG|NGY}},{{TAG|NGZ}},{{TAG|NGXF}},{{TAG|NGYF}},{{TAG|NGZF}}):
:For {{TAG|PREC}}=High and {{TAG|PREC}}=Accurate, wrap around errors are avoided, all ''G''-vectors that are twice as large as the vectors included in the basis set are taken into account in the FFT's. For {{TAG|PREC}}=Low, {{TAG|PREC}}=Medium or {{TAG|PREC}}=Normal, the FFT grids are reduced, and 3/4 of the required values are used. Usually {{TAG|PREC}}=Medium and {{TAG|PREC}}=Normal, are sufficiently accurate with errors less than 1 meV/atom.
:For {{TAG|PREC}}=High and {{TAG|PREC}}=Accurate, wrap around errors are avoided, all ''G''-vectors that are twice as large as the vectors included in the basis set are taken into account in the FFT's. For {{TAG|PREC}}=Low, {{TAG|PREC}}=Medium or {{TAG|PREC}}=Normal, the FFT grids are reduced, and 3/4 of the required values are used. Usually {{TAG|PREC}}=Medium and {{TAG|PREC}}=Normal, are sufficiently accurate with errors less than 1 meV/atom.


:In addition, the {{TAG|PREC}} tag determines the spacing for the grids representing the augmentation charges, charge densities and potentials ({{TAG|NGXF}}, {{TAG|NGYF}}, {{TAG|NGZF}}). For {{TAG|PREC}}=Accurate and {{TAG|PREC}}=Normal, the support grid contains twice as many points in each direction as the grids for the orbitals ({{TAG|NGXF}}=2&times;{{TAG|NGX}}, {{TAG|NGYF}}=2&times;{{TAG|NGY}}, {{TAG|NGZF}}=2&times;{{TAG|NGZ}}). {{TAG|PREC}}=Single is identical to {{TAG|PREC}}=Normal, execpt that the double grid technique is not applied. This is convenient of you need to cut down on storage demands, or want to reduce the size of the {{FILE|CHG}} and {{FILE|CHGCAR}} files (for scanning tunneling microscopy simulation, it is recommended to use {{TAG|PREC}}=Single). In all other cases, {{TAG|NGXF}}, {{TAG|NGYF}}, and {{TAG|NGZF}} are determined by some rather heuristic formula from {{TAG|ENAUG}}.
:In addition, the {{TAG|PREC}} tag determines the spacing for the grids representing the augmentation charges, charge densities and potentials ({{TAG|NGXF}}, {{TAG|NGYF}}, {{TAG|NGZF}}). For {{TAG|PREC}}=Accurate and {{TAG|PREC}}=Normal, the support grid contains twice as many points in each direction as the grids for the orbitals ({{TAG|NGXF}}=2&times;{{TAG|NGX}}, {{TAG|NGYF}}=2&times;{{TAG|NGY}}, {{TAG|NGZF}}=2&times;{{TAG|NGZ}}). In all other cases, {{TAG|NGXF}}, {{TAG|NGYF}}, and {{TAG|NGZF}} are determined by some heuristic formula from {{TAG|ENAUG}} (these settings should be avoided).
 
:{{TAG|PREC}}=Single is uses the same FFT grids as  {{TAG|PREC}}=Normal in VASP.5 and the same FFT grids as {{TAG|PREC}}=Accurate in  VASP.6. The double grid technique is not used. This is convenient, if one needs to cut down on storage demands, or if one wants to reduce the size of the files {{FILE|CHG}} and {{FILE|CHGCAR}}. Furthermore, {{TAG|PREC}}=Single avoids high frequency oscillations caused the double grid technique and the resulting Fourier interpolation. It is often expedient for scanning tunneling simulations or other calculations where high frequency wiggles of the charge density in the vacuum region are undesirable.


*Real space projectors ({{TAG|ROPT}}):
*Real space projectors ({{TAG|ROPT}}):
Line 47: Line 73:
|-
|-
| Normal || 1.0 || || Normal || -5&times;10<sup>-4</sup>
| Normal || 1.0 || || Normal || -5&times;10<sup>-4</sup>
|-
| Single || 1.0 || || Single || -5&times;10<sup>-4</sup>
|-
|-
| Accurate || 1.0 || || Accurate || -2.5&times;10<sup>-4</sup>
| Accurate || 1.0 || || Accurate || -2.5&times;10<sup>-4</sup>
Line 62: Line 90:
{{TAG|ENCUT}},
{{TAG|ENCUT}},
{{TAG|ENAUG}},
{{TAG|ENAUG}},
{{TAG|ENMAX}}
{{TAG|ENMAX}},
{{TAG|PRECFOCK}}


{{sc|PREC|Examples|Examples that use this tag}} 
----
----
[[The_VASP_Manual|Contents]]


[[Category:INCAR]]
[[Category:INCAR]][[Category:Electronic Minimization]][[Category:Electronic Minimization Methods]]

Revision as of 06:50, 26 March 2021

PREC = Low | Medium | High | Normal | Single | Accurate 

Default: PREC = Medium for VASP.4.X
= Normal for VASP.5.X

Description: PREC specifies the "precision"-mode.


We recommend to use PREC=Normal or PREC=Accurate. The settings PREC=Medium and PREC=High (available already in VASP.4.4) are no longer recommended. PREC=Normal can be used for most routine calculations, whereas PREC=Accurate reduces egg-box effects and errors even further. Specifically, PREC=Accurate strictly avoids any aliasing or wrap around errors and uses an augmentation grid (NGXF×NGYF×NGZF) that is exactly twice as large as the coarser grid (NGX×NGY×NGZ) used for the representation of the pseudo orbitals. PREC=Accurate increases the memory requirements somewhat, but it should be used if very accurate forces (phonons and second derivatives) are required. The accuracy of forces can be sometimes further improved by specifying ADDGRID=.TRUE. (but user reports are somewhat contradictory, whether this really helps).

The use of PREC=High is no longer recommend (and is only available for compatibility). Essentially PREC=High only increases the energy cutoff by 30 %. The same effect can be achieved by just manually increasing ENCUT in the INCAR file. If, additionally, very accurate forces are required, PREC=Accurate can be used in combination with an increased energy cutoff.

Mind: we recommend to specify the energy cutoff (ENCUT) always manually in the INCAR file, to avoid incompatibilities between calculations.


The PREC parameter sets the defaults for four sets of parameters (ENCUT; NGX, NGY, NGZ; NGXF, NGYF, NGZF, and ROPT):

PREC ENCUT NGX,Y,Z NGXF,YF,ZF ROPT
Normal max(ENMAX) 3/2×Gcut NGX -5×10-4
Single (VASP.5) max(ENMAX) 3/2×Gcut NGX -5×10-4
Single (VASP.6) max(ENMAX) 2 ×Gcut NGX -5×10-4
SingleN (VASP.6 development) max(ENMAX) 3/2×Gcut NGX -5×10-4
Accurate max(ENMAX) 2×Gcut NGX -2.5×10-4
old settings:
Low max(ENMIN) 3/2×Gcut 3×Gaug -1×10-2
Medium max(ENMAX) 3/2×Gcut 4×Gaug -2×10-3
High 1.3×max(ENMAX) 2×Gcut 16/3×Gaug -4×10-4

where max(ENMAX) and max(ENMIN) are the maximum ENMAX and ENMIN found on the POTCAR file, and

with Ecut=ENCUT and Eaug=ENAUG.

Mind: the settiLow ||align="center"|max(ENMIN) || 3/2×Gcut ||align="center"|3×Gaug ||-1×10-2ngs PREC=Normal and PREC=Accurate are only available in VASP.4.5 and newer versions. The setting PREC=Single is only available as of VASP.5.1, and has been updated in VASP.6.

  • Plane wave kinetic energy cutoff (ENCUT):
The PREC-tag determines the energy cutoff ENCUT, if ENCUT is not specified in the INCAR file. For PREC=Low, ENCUT will be set to the maximal ENMIN value found on the POTCAR file. For PREC=Medium and PREC=Accurate, ENCUT will be set to maximal ENMAX value found on the POTCAR file. Finally for PREC=High, ENCUT is set to the maximal ENMAX value on the POTCAR file plus 30%. PREC=High guarantees that the absolute energies are converged to a few meV, and it ensures that the stress tensor is converged within a few kBar. In general, an increased energy cutoff is only required for accurate evaluation of quantities related to the stress tensor (e.g. elastic properties).
For PREC=High and PREC=Accurate, wrap around errors are avoided, all G-vectors that are twice as large as the vectors included in the basis set are taken into account in the FFT's. For PREC=Low, PREC=Medium or PREC=Normal, the FFT grids are reduced, and 3/4 of the required values are used. Usually PREC=Medium and PREC=Normal, are sufficiently accurate with errors less than 1 meV/atom.
In addition, the PREC tag determines the spacing for the grids representing the augmentation charges, charge densities and potentials (NGXF, NGYF, NGZF). For PREC=Accurate and PREC=Normal, the support grid contains twice as many points in each direction as the grids for the orbitals (NGXF=2×NGX, NGYF=2×NGY, NGZF=2×NGZ). In all other cases, NGXF, NGYF, and NGZF are determined by some heuristic formula from ENAUG (these settings should be avoided).
PREC=Single is uses the same FFT grids as PREC=Normal in VASP.5 and the same FFT grids as PREC=Accurate in VASP.6. The double grid technique is not used. This is convenient, if one needs to cut down on storage demands, or if one wants to reduce the size of the files CHG and CHGCAR. Furthermore, PREC=Single avoids high frequency oscillations caused the double grid technique and the resulting Fourier interpolation. It is often expedient for scanning tunneling simulations or other calculations where high frequency wiggles of the charge density in the vacuum region are undesirable.
  • Real space projectors (ROPT):
LREAL=O LREAL=A
PREC ROPT PREC ROPT
Low 0.67 Low -1×10-2
Med 1.0 Med -2×10-3
Normal 1.0 Normal -5×10-4
Single 1.0 Single -5×10-4
Accurate 1.0 Accurate -2.5×10-4
High 1.5 High -4×10-4

Related Tags and Sections

NGX, NGY, NGZ, NGXF, NGYF, NGZF, ENCUT, ENAUG, ENMAX, PRECFOCK

Examples that use this tag