In principle, the PAW method can be used in the same manner as the US-PP method. Only special PAW POTCAR files are required. In principle, also no additional user interference is required. However there are a few flags that control the behavior of the PAW implementation. The first one is LMAXPAW:
LMAXPAW = LThis flag sets the maximum -quantum number for the evaluation of the on-site terms on the radial support grids in the PAW method. The default for LMAXPAW is , where is the maximum angular quantum number of the partial waves. Useful settings for LMAXPAW are for instance:
LMAXPAW = 0In this case, only spherical terms are evaluated on the radial grid. This does not mean that a-spherical terms are totally neglected, because the compensation charges are always expanded up to on the plane wave grid.
Finally, LMAXPAW=-1 has a special meaning. For LMAXPAW=-1, no on-site correction terms are evaluated on the radial support grid, which effectively means that the behavior of US-PP's is recovered with PAW input datasets. Usually this allows very efficient and fast calculations, and this switch might be of interest for relaxations and molecular dynamics runs. Energies should be evaluated with the default setting for LMAXPAW.
An additional flag controls up to which l quantum number the onsite PAW charge densities are passed through the charge density mixer and written to the CHGCAR file:
LMAXMIX = lThe default is LMAXMIX=2. Higher l-quantum numbers are usually not handled by the mixer, i.e. a straight mixing is applied for them (the PAW on-site charge density for higher l quantum numbers is reset precisely to the value corresponding to the present orbitals). Usually, it is not required to increase LMAXMIX, but the following two cases are exceptions:
The second switch, that is useful in the context of the PAW method (and US-PP) is ADDGRID. The default is ADDGRID=.FALSE. If
ADDGRID = .TRUE.is written in the INCAR file, an additional (third) support grid is used for the evaluation of the augmentation charges. This third grid contains 8 times more points than the fine grid NGXF, NGYF, NGZF. Whenever terms involving augmentation charges are evaluated, this third grid is used. For instance: The augmentation charge is evaluated first in real space on this fine grid, FFT-transformed to reciprocal space and then added to the total charge density on the grid NGXF, NGYF, NGZF. The additional grid helps to reduce the noise in the forces significantly. In many cases, it even allows to perform calculations in which NGXF=NGX etc. This can be achieved by setting
ENAUG = 1 ; ADDGRID = .TRUE.in the INCAR file. The flag can also be used for US-PPs, in particular, to reduce the noise in the forces.