DFT-ulg

In the DFT-ulg method of Kim et al., the correction term takes the form:

$E_{\mathrm{disp}} = -\frac{1}{2} s_{lg}\sum_{i=1}^{N_{at}} \sum_{j=1}^{N_{at}} \sum_{\mathbf{L}} {}^{\prime} \frac{2D_{0}^{ij}(R_{0}^{ij})^{6}}{r_{ij,L}^{6}+b_{lg}(R_{0}^{ij})^{6}}$

where the first two summations are over all $N_{at}$ atoms in the unit cell and the third summation is over all translations of the unit cell ${\mathbf{L}}=(l_1,l_2,l_3)$ where the prime indicates that $i\not=j$ for ${\mathbf{L}}=0$.

• VDW_RADIUS=50.0 : cutoff radius (in $\AA$) for pair interactions
• VDW_S6=0.7012 : global scaling factor $s_{lg}$ (available in VASP.5.3.5 and later)
• VDW_D=0.6966 : damping parameter $b_{lg}$
• VDW_C6=[real array] : $C_6$ parameters ($\mathrm{Jnm}^{6}\mathrm{mol}^{-1}$) for each species defined in the POSCAR file
• VDW_R0=[real array] : $R_0$ parameters ($\AA$) for each species defined in the POSCAR file
• LVDW_EWALD=.FALSE. : the lattice summation in $E_{\mathrm{disp}}$ expression is computed by means of Ewald's summation (.TRUE. ) or via a real space summation over all atomic pairs within cutoff radius VDW_RADIUS (.FALSE.). (available in VASP.5.3.5 and later)

Related tags and articles

VDW_RADIUS, VDW_S6, VDW_D, VDW_C6, VDW_R0, LVDW_EWALD, IVDW

---