Coulomb singularity

From VASP Wiki
Revision as of 16:13, 11 May 2022 by Ftran (talk | contribs)

The bare Coulomb operator

in the unscreened HF exchange has a representation in the reciprocal space that is given by

It has an (integrable) singularity at that leads to a very slow convergence of the results with respect to the cell size or number of k-points. In order to alleviate this issue different methods have been proposed: the auxiliary function [1], probe-charge Ewald [2] (HFALPHA), and Coulomb truncation[3] methods (selected with HFRCUT). These mostly involve modifying the Coulomb Kernel in a way that yields the same result as the unmodified kernel in the limit of large supercell sizes. These methods can also be applied to the Thomas-Fermi and error function screened Coulomb operators given by

and

respectively, whose representations in the reciprocal space are given by

and

respectively.

Auxiliary function

In this approach an auxiliary periodic function with the same divergence as the Coulomb potential in reciprocal space is subtracted in the k points used to integrate the Hartree-Fock energy, thus regularizing the integral[1]. This function is chosen such that it has a closed analytical expression for its integral[1] or the integral is evaluated numerically[4]. This approach is currently not implemented in VASP, instead, the probe-charge Ewald method is used.

Probe-charge Ewald

A similar approach to the auxiliary function method described above is the probe-charge Ewald method [2]. In this case, the auxiliary function is chosen to have the form of the Coulomb kernel times a Gaussian function with a width (HFALPHA) comparable to the Brillouin zone diameter. This function is used to regularize the Coulomb integral that is evaluated in the regular k point grid with the divergent part being evaluated by analytical integration of the Coulomb kernel (see eq. 29 in ref. [2]). The value of the integral of the bare Coulomb potential is (see eq. 31 in ref. [2])

for the Thomas-Fermi and error function screened Coulomb kernels we have

and

respectively.

Spherical truncation

In this method[3] the bare Coulomb operator is spherically truncated by multiplying it by the step function , and in the reciprocal this leads to

whose value at is finite and is given by , where the truncation radius is chosen as with being the number of -points in the full Brillouin zone.

The screened potentials have no singularity at . Nevertheless, it is still beneficial for accelerating the convergence with respect to the number of k-points to multiply these screened operators by , which in the reciprocal space gives

and

respectively, with the following values at :

and

Note that the spherical truncation method described above works very well in the case of 3D systems. However, it is not recommended for systems with a lower dimensionality[5]. For such systems, the approach proposed in ref. [5] (not implemented in VASP) is more adapted since the truncation is done according to the Wigner-Seitz cell and therefore more general.

Related tags and articles

HFRCUT, FOCKCORR, Hybrid_functionals: formalism, Downsampling_of_the_Hartree-Fock_operator

References