ADDGRID: Difference between revisions
No edit summary |
No edit summary |
||
| Line 3: | Line 3: | ||
Description: {{TAG|ADDGRID}} determines whether an additional support grid is used for the evaluation of the augmentation charges. | Description: {{TAG|ADDGRID}} determines whether an additional support grid is used for the evaluation of the augmentation charges. | ||
---- | ---- | ||
When {{TAG|ADDGRID}}=.TRUE. VASP uses an additional support grid for the evaluation of the augmentation charges. This grid contains 8 times more points than the standard "fine" grid ({{TAG|NGXF}}×{{TAG|NGYF}}×{{TAG|NGZF}}). Whenever terms involving augmentation charges are evaluated, this additional grid is used. For instance: The augmentation charge is evaluated first in real space on this additional grid, FFT-transformed to reciprocal space, and then added to the total charge density on the standard "fine" grid ({{TAG|NGXF}}×{{TAG|NGYF}}×{{TAG|NGZF}}). The additional grid often helps to reduce the noise in the forces. In some cases, it even allows to perform calculations with {{TAG|NGXF}}={{TAG|NGX}}. | When {{TAG|ADDGRID}}=.TRUE. VASP uses an additional support grid for the evaluation of the augmentation charges. This grid contains 8 times more points than the standard "fine" grid ({{TAG|NGXF}}×{{TAG|NGYF}}×{{TAG|NGZF}}). Whenever terms involving augmentation charges are evaluated, this additional grid is used. For instance: The augmentation charge is evaluated first in real space on this additional grid, FFT-transformed to reciprocal space, and then added to the total charge density on the standard "fine" grid ({{TAG|NGXF}}×{{TAG|NGYF}}×{{TAG|NGZF}}). The additional grid often helps to reduce the noise in the forces. In some cases, it even allows to perform calculations with {{TAG|NGXF}}{{=}}{{TAG|NGX}}. | ||
{{NB|important|If there is any contribution in the density or potential at the highest Fourier component <math>G</math> of | |||
If there is any contribution in the density or potential at the highest Fourier component <math>G</math> of | |||
the conventional fine grid (given by {{TAG|NGXF}}×{{TAG|NGYF}}×{{TAG|NGZF}}), then Fourier interpolation to twice the grid density leads to | the conventional fine grid (given by {{TAG|NGXF}}×{{TAG|NGYF}}×{{TAG|NGZF}}), then Fourier interpolation to twice the grid density leads to | ||
oscillations in real space. These oscillations correspond to the largest wave vector <math> G_{cut} </math> i.e. <math>e^{i G_{cut} r}</math>. | oscillations in real space. These oscillations correspond to the largest wave vector <math> G_{cut} </math> i.e. <math>e^{i G_{cut} r}</math>. | ||
| Line 15: | Line 13: | ||
The upshot is that Fourier-interpolation can lead to termination wiggles with oscillations <math>e^{i G_{cut} r}</math> in the interpolated potential (where <math> G_{cut}</math> corresponds to | The upshot is that Fourier-interpolation can lead to termination wiggles with oscillations <math>e^{i G_{cut} r}</math> in the interpolated potential (where <math> G_{cut}</math> corresponds to | ||
the largest Fourier components on the fine grid). Fourier smoothing, which is in essence used for the augmentation densities, is generally less problematic, but it can also result in negative density in real space. | the largest Fourier components on the fine grid). Fourier smoothing, which is in essence used for the augmentation densities, is generally less problematic, but it can also result in negative density in real space. | ||
Therefore, we recommend performing careful tests, on whether {{TAG|ADDGRID}} works as desired; please do not use this tag as default in all your calculations! | Therefore, we recommend performing careful tests, on whether {{TAG|ADDGRID}} works as desired; please do not use this tag as default in all your calculations!}} | ||
== Related tags and articles == | == Related tags and articles == | ||
Revision as of 11:32, 6 May 2026
ADDGRID = .TRUE. | .FALSE.
Default: ADDGRID = .FALSE.
Description: ADDGRID determines whether an additional support grid is used for the evaluation of the augmentation charges.
When ADDGRID=.TRUE. VASP uses an additional support grid for the evaluation of the augmentation charges. This grid contains 8 times more points than the standard "fine" grid (NGXF×NGYF×NGZF). Whenever terms involving augmentation charges are evaluated, this additional grid is used. For instance: The augmentation charge is evaluated first in real space on this additional grid, FFT-transformed to reciprocal space, and then added to the total charge density on the standard "fine" grid (NGXF×NGYF×NGZF). The additional grid often helps to reduce the noise in the forces. In some cases, it even allows to perform calculations with NGXF=NGX.
| Important: If there is any contribution in the density or potential at the highest Fourier component [math]\displaystyle{ G }[/math] of
the conventional fine grid (given by NGXF×NGYF×NGZF), then Fourier interpolation to twice the grid density leads to oscillations in real space. These oscillations correspond to the largest wave vector [math]\displaystyle{ G_{cut} }[/math] i.e. [math]\displaystyle{ e^{i G_{cut} r} }[/math]. In real space, the charge density or potential will therefore alternate between positive and negative values on the ultra-fine grid, in particular, in regions where the density or potential are small. The terminus techniques is "termination wiggles". Although this is a somewhat oversimplified presentation, it is fairly straightforward to derive more rigorous results in 1D. The upshot is that Fourier-interpolation can lead to termination wiggles with oscillations [math]\displaystyle{ e^{i G_{cut} r} }[/math] in the interpolated potential (where [math]\displaystyle{ G_{cut} }[/math] corresponds to the largest Fourier components on the fine grid). Fourier smoothing, which is in essence used for the augmentation densities, is generally less problematic, but it can also result in negative density in real space. Therefore, we recommend performing careful tests, on whether ADDGRID works as desired; please do not use this tag as default in all your calculations! |
Related tags and articles
PREC, NGX, NGY, NGZ, NGXF, NGYF, NGZF, ENCUT, ENAUG, ENMAX, PRECFOCK