Requests for technical support from the VASP group should be posted in the VASP-forum.

# Self-consistent screening in Tkatchenko-Scheffler method

A computationally efficient way to account for electrodynamic response effects, in particular the interaction of atoms with the dynamic electric field due to the surrounding polarizable atoms was proposed by Tkatchenko et al^{[1]}. In this method, termed TS+SCS, the frequency-dependent screened polarizabilities are obtained by solving the self-consistent screening equation:

where is the dipole-dipole interaction tensor and is the effective frequency-dependent polarizability, approximated by

with the characteristic mean excitation frequency . The dispersion coefficients are computed from the Casimir-Polder integral:

The van der Waals radii of atoms are obtained by rescaling the radii computed using DFT-TS:

The dispersion energy is computed using the same equation as in the original Tkatchenko-Scheffler method but with corrected parameters , , and . The TS+SCS method is invoked by setting IVDW=2|20 and LVDWSCS=*.TRUE.*. In addition to parameters controlling the Tkatchenko-Scheffler method, the following optional parameters can set by the user:

- VDW_SR=0.97 scaling factor
- SCSRAD=120 cutoff radius (in ) used in the calculation of
- LSCSGRAD=.TRUE. decides whether to compute SCS contribution to gradients (LSCSGRAD=
*.TRUE.*) or not - LSCALER0=.TRUE. decider whether to use the equation above for to re-scale the parameter (LSCALER0=
*.TRUE.*) or not

Details of the implementation of the TS+SCS method in VASP and the performance tests made on various crystalline systems are presented in reference ^{[1]}.

## IMPORTANT NOTES

- This method requires the use of POTCAR files from the PAW dataset version 52 or later.
- This method is incompatible with the setting ADDGRID=
*.TRUE.*. - This type of calculation may be time-consuming for large systems. Note that the SCS contribution to gradients and stress tensor is only modest (but non-negligible) in many cases. In the initial stages of relaxation of large systems, or if only energy is of interest, the calculation can be accelerated by setting LSCSGRAD=
*.FALSE.*. - The default value for the parameter VDW_SR (which is, in general, different from that used in the unscreened Tkatchenko-Scheffler method method) is available only for the PBE functional. If a functional other than PBE is used, the value for VDW_SR must be specified in the INCAR file.

## Related Tags and Sections

IVDW, DFT-D2, DFT-D3, Tkatchenko-Scheffler method, Tkatchenko-Scheffler method with iterative Hirshfeld partitioning, Many-body dispersion energy, dDsC dispersion correction