The Tkatchenko-Scheffler (TS) dispersion correction method which uses fixed neutral atoms as a reference to estimate the effective volumes of atoms-in-molecule (AIM) and to calibrate their polarizabilities and dispersion coefficients (see Sec. 6.77.3) fails to describe the structure and the energetics of ionic solids. As shown in Ref. [133,134], this problem can be solved by replacing the conventional Hirshfeld partitioning used to compute properties of interacting atoms by the iterative scheme proposed by Bultinck . In this iterative Hirshfeld algorithm (HI), the neutral reference atoms are replaced with ions with fractional charges determined together with the AIM charge densities in a single iterative procedure. The algorithm is initialized with a promolecular density defined by non-interacting neutral atoms. The iterative procedure then runs in the following steps:
where lint expresses the integer part of and uintlint.
Steps (1) to (3) are iterated until the difference in the electronic populations between two subsequent steps ( ) is less than a predefined threshold for all atoms. The converged interative Hirshfeld weights () are then used to define the AIM properties needed to evaluate dispersion energy, see Sec. 6.77.3.
The DFT-TS calculation with iterative Hirshfeld partitioning (DFT-TS/HI) is invoked by setting IVDW=21. The convergence criterion for iterative Hirshfeld partitioning (in e) can optionally be defined via parameter HITOLER (the default value is 5e-5). Other optional parameters controlling the input for the calculation are as in the conventional TS method (Sec. 6.77.3). The default value of the adjustable parameter VDW_SR is 0.95 and corresponds to the PBE functional.
The PBE-TS/HI method is described in detail in J. Chem. Theory Comput. 9, 4293 (2013)
and its performance
in optimization of
various crystalline systems is examined in
J. Chem. Phys. 141, 034114 (2014).