LREAL: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
Line 27: Line 27:
This expression can be evaluated in reciprocal or real space: In reciprocal space (second line) the number of operations scales with the size of the basis set (i.e. number of plane-waves). In real space (first line) the projection-operators are confined to spheres around each atom. Therefore the number of operations necessary to evaluate one C<sub>in'''k'''</sub> does not increase with the system size (usually the number of grid points within the cut-off-sphere is between 500 and 2000). One of the major obstacles of the method working in real space is that the projection operators must be optimized, i.e. all high Fourier components must be removed from the projection operators. If this is not done "aliasing" can happen (i.e. the high Fourier components of the projection operators are downfolded to low Fourier components and a random noise is introduced).
This expression can be evaluated in reciprocal or real space: In reciprocal space (second line) the number of operations scales with the size of the basis set (i.e. number of plane-waves). In real space (first line) the projection-operators are confined to spheres around each atom. Therefore the number of operations necessary to evaluate one C<sub>in'''k'''</sub> does not increase with the system size (usually the number of grid points within the cut-off-sphere is between 500 and 2000). One of the major obstacles of the method working in real space is that the projection operators must be optimized, i.e. all high Fourier components must be removed from the projection operators. If this is not done "aliasing" can happen (i.e. the high Fourier components of the projection operators are downfolded to low Fourier components and a random noise is introduced).


Currently VASP supports three different schemes to remove the high Fourier components from the projectors. {{TAG|LREAL}}=.TRUE. is the simplest one. For {{TAG|LREAL}}=.TRUE., the real space projectors that have been generated by the pseudopotential generation code are used. This requires no user interference but is potentially very inaccurate. For the outdated {{TAG|LREAL}}=On, the real space projectors are optimized by VASP using an algorithm proposed by King-Smith et al.<ref name="kingsmith:prb:91"/> For the recommended {{TAG|LREAL}}=Auto, a new scheme<ref name="kresse:tobepublished"/> is used which is considerably better (resulting in more localized) projector functions than the King-Smith et al. method. To fine-tune the optimization procedure the flag {{TAG|ROPT}} can and should be used, if {{TAG|LREAL}}=Auto (or {{TAG|LREAL}}=On) is used. Specifically, perform first reference calculations using {{TAG|LREAL}}=.False. and decrease {{TAG|ROPT}} until an acceptable accuracy of say 1 meV/atom is attained. Please also check carefully the documentation for {{TAG|ROPT}}.
Currently VASP supports three different schemes to remove the high Fourier components from the projectors. {{TAG|LREAL}}=.TRUE. is the simplest one. For {{TAG|LREAL}}=.TRUE., the real space projectors that have been generated by the pseudopotential generation code are used. This requires no user interference but is potentially very inaccurate. For the outdated {{TAG|LREAL}}=On, the real space projectors are optimized by VASP using an algorithm proposed by King-Smith et al.<ref name="kingsmith:prb:91"/> For the recommended {{TAG|LREAL}}=Auto, an unpublished scheme<ref name="kresse:tobepublished"/> is used which  
results in more accurate and localized projector functions than for the King-Smith et al. method. To fine-tune the optimization procedure the flag {{TAG|ROPT}} can and should be used, if {{TAG|LREAL}}=Auto (or {{TAG|LREAL}}=On) is used. Specifically, perform first reference calculations using {{TAG|LREAL}}=.False. and decrease {{TAG|ROPT}} until an acceptable accuracy of say 1 meV/atom is attained. Please also check carefully the documentation for {{TAG|ROPT}}.


We recommend to use the real-space projection scheme for systems containing more than about 30 atoms. We also strongly recommend to use only {{TAG|LREAL}}=Auto.
We recommend to use the real-space projection scheme for systems containing more than about 30 atoms. We also strongly recommend to use only {{TAG|LREAL}}=Auto.

Revision as of 08:32, 11 November 2020

LREAL = .FALSE. | Auto (or A) | On (or O) | .TRUE.
Default: LREAL = .FALSE. 

Description: LREAL determines whether the projection operators are evaluated in real-space or in reciprocal space.


LREAL=.FALSE. projection done in reciprocal space
LREAL=Auto or A projection done in real space, fully automatic optimization of projection operators (little to no user interference required)
LREAL=On or O projection done in real space, projection operators are re-optimized (not recommended)
LREAL=.TRUE. projection done in real space, use projectors on file (not recommended)


The non local part of the pseudopotential requires the evaluation of an expression:


.

where the "projected wavefunction character" is defined as:

This expression can be evaluated in reciprocal or real space: In reciprocal space (second line) the number of operations scales with the size of the basis set (i.e. number of plane-waves). In real space (first line) the projection-operators are confined to spheres around each atom. Therefore the number of operations necessary to evaluate one Cink does not increase with the system size (usually the number of grid points within the cut-off-sphere is between 500 and 2000). One of the major obstacles of the method working in real space is that the projection operators must be optimized, i.e. all high Fourier components must be removed from the projection operators. If this is not done "aliasing" can happen (i.e. the high Fourier components of the projection operators are downfolded to low Fourier components and a random noise is introduced).

Currently VASP supports three different schemes to remove the high Fourier components from the projectors. LREAL=.TRUE. is the simplest one. For LREAL=.TRUE., the real space projectors that have been generated by the pseudopotential generation code are used. This requires no user interference but is potentially very inaccurate. For the outdated LREAL=On, the real space projectors are optimized by VASP using an algorithm proposed by King-Smith et al.[1] For the recommended LREAL=Auto, an unpublished scheme[2] is used which

results in more accurate and localized projector functions than for the King-Smith et al. method. To fine-tune the optimization procedure the flag ROPT can and should be used, if LREAL=Auto (or LREAL=On) is used. Specifically, perform first reference calculations using LREAL=.False. and decrease ROPT until an acceptable accuracy of say 1 meV/atom is attained. Please also check carefully the documentation for ROPT.

We recommend to use the real-space projection scheme for systems containing more than about 30 atoms. We also strongly recommend to use only LREAL=Auto.

For LREAL=A (and LREAL=O) the projection operators are optimized by VASP on the fly (i.e. on startup). Several flags influence the optimization:

  • ENCUT (i.e. the energy cutoff), components beyond the energy cutoff are 'removed' from the projection operators.
  • PREC tag specifies how precise the real space projectors should be and sets the variables ROPT accordingly to the following values:
For LREAL=Auto
PREC=Fast accuracy 2 ·  10-3 (ROPT=0.002)
PREC=Normal accuracy 5 · 10-4 (ROPT=5E-4)
PREC=Accurate accuracy 2.5 · 10-4 (ROPT=2.5E-4)
PREC=Low accuracy 1 · 10-2 (ROPT=0.01)
PREC=Med accuracy 2 · 10-3 (ROPT=0.002)
PREC=High accuracy 2 · 10-4 (ROPT=2E-4)
For LREAL=On
PREC=Low 700 points in the real space sphere (ROPT=0.67)
PREC=Med 1000 points in the real space sphere (ROPT=1.0)
PREC=High 1500 points in the real space sphere (ROPT=1.5)


These defaults can be superseded by specifying the ROPT-tag in the INCAR file.


Mind: Real space optimization (LREAL=.TRUE., On, or Auto) always results in a small (not necessarily negligible) error (the error is usually a constant energy shift for each atom). If you are interested in energy differences of a few meV use only calculations with the same setup (i.e. same ENCUT, PREC, LREAL and ROPT setting) for all calculations. For example, if you want to calculate surface energies recalculate the bulk groundstate energy with exactly the same setting you are going to use for the surface. Another possibility is to relax the structure using real space projection and to perform one final total energy calculation using LREAL=.FALSE. to get "exact" energies. Anyway, for PREC=Normal, the errors introduced by the real space projection are usually of the same order magnitude as those introduced by the wrap around errors. For PREC=Accurate errors are usually less than meV. PREC=Fast should be used only for say molecular dyanamics, if compute resources are really an issue.

Related Tags and Sections

ROPT, PREC

Examples that use this tag

References