View source for 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{{cite|tkatchenko:prl:12}}. In this method, termed TS+SCS, the frequency-dependent screened polarizabilities <math>\alpha^{SCS}(\omega)</math> are obtained by solving the self-consistent screening equation:

:<math>\alpha_{i}^{SCS}(\omega) = \alpha_{i}(\omega) - \alpha_{i}(\omega) 
\sum_{i \neq j} \tau_{ij} \alpha_{j}^{SCS}(\omega)</math>

where <math>\tau_{ij}</math> is the dipole-dipole interaction tensor and <math>\alpha_{i}(\omega)</math> is the effective frequency-dependent polarizability, approximated by

:<math>\alpha_{i}(\omega) = \frac{\alpha_{i}}{1+\left (
\omega / \omega_i \right )^2}</math>

with the characteristic mean excitation frequency <math>\omega_i = \frac{4}{3} \frac{C_{6ii}}{(\alpha_{i})^2}</math>. The dispersion coefficients are computed from the Casimir-Polder integral:

:<math>C_{6ii} = \frac{3}{\pi} \int_0^{\infty} \alpha_{i}^{SCS}(\omega) 
\alpha_{i}^{SCS}(\omega) \,d\omega.</math>

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

:<math>R_{0i}^{SCS} = \left ( \frac{\alpha_{i}^{SCS}}{\alpha_{i}} \right )^{1/3} R_{0i}. </math>

The dispersion energy is computed using the same equation as in the original {{TAG|Tkatchenko-Scheffler method}} but with corrected parameters <math>C_{6ii}^{SCS}</math>, <math>\alpha_{i}^{SCS}</math>, and <math>R_{0i}^{SCS}</math>. The TS+SCS method is invoked by setting {{TAG|IVDW}}=2|20 and {{TAG|LVDWSCS}}=''.TRUE.''. In addition to parameters controlling the {{TAG|Tkatchenko-Scheffler method}}, the following optional parameters can set by the user:

*{{TAG|VDW_SR}}=0.97 scaling factor <math>s_{R}</math>
*{{TAG|SCSRAD}}=120 cutoff radius (in <math>\AA</math>) used in the calculation of <math>\tau_{ij}</math>
*{{TAG|LSCSGRAD}}=.TRUE. decides whether to compute SCS contribution to gradients ({{TAG|LSCSGRAD}}=''.TRUE.'') or not
*{{TAG|LSCALER0}}=.TRUE. decides whether to use the equation above for <math>R_{0i}^{SCS}</math> to re-scale the parameter <math>R_{0}</math> ({{TAG|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 {{cite|bucko:prb:13}}.

{{NB|mind|This method requires the use of {{TAG|POTCAR}} files from the PAW dataset version 52 or later.}}
{{NB|mind|This method is incompatible with the setting {{TAG|ADDGRID}}{{=}}''.TRUE.''.}}
{{NB|mind|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 {{TAG|LSCSGRAD}}{{=}}''.FALSE.''.}}
{{NB|mind|The default value for the parameter {{TAG|VDW_SR}} (which is, in general, different from that used in the unscreened {{TAG|Tkatchenko-Scheffler method}} method) is available only for the PBE functional. If a functional other than PBE is used, the value for {{TAG|VDW_SR}} must be specified in the {{TAG|INCAR}} file.}}

== Related Tags and Sections ==
{{TAG|LVDWSCS}}
{{TAG|VDW_SR}}
{{TAG|SCSRAD}}
{{TAG|LSCSGRAD}}
{{TAG|LSCALER0}}
{{TAG|IVDW}},
{{TAG|DFT-D2}},
{{TAG|DFT-D3}},
{{TAG|Tkatchenko-Scheffler method}},
{{TAG|Tkatchenko-Scheffler method with iterative Hirshfeld partitioning}},
{{TAG|Many-body dispersion energy}},
{{TAG|dDsC dispersion correction}}

== References ==
<references/>

----
[[Category:Exchange-correlation functionals]][[Category:van der Waals functionals]][[Category:Theory]]