Coulomb singularity: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
Line 7: Line 7:
V(q)=\frac{4\pi}{q^2}
V(q)=\frac{4\pi}{q^2}
</math>
</math>
It has an (integrable) singularity at <math>q=\vert\mathbf{k}'-\mathbf{k}+\mathbf{G}\vert=0</math> that leads to a very . In order to alleviate this issue and to improve the convergence of the exact exchange with respect to the supercell size (or the k-point mesh density) different methods have been proposed: the auxiliary function {{cite|gygi:prb:86}}, probe-charge Ewald {{cite|massidda:prb:93}} ({{TAG|HFALPHA}}), and Coulomb truncation{{cite|spenceralavi:prb:08}} methods (selected with {{TAG|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, which can also be applied to screened Coulomb operators, are described below.
It has an (integrable) singularity at <math>q=\vert\mathbf{k}'-\mathbf{k}+\mathbf{G}\vert=0</math> 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 and to improve the convergence of the exact exchange with respect to the supercell size (or the k-point mesh density) different methods have been proposed: the auxiliary function {{cite|gygi:prb:86}}, probe-charge Ewald {{cite|massidda:prb:93}} ({{TAG|HFALPHA}}), and Coulomb truncation{{cite|spenceralavi:prb:08}} methods (selected with {{TAG|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, which can also be applied to screened Coulomb operators, are described below.


=== Auxiliary function ===
=== Auxiliary function ===

Revision as of 12:42, 10 May 2022

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 and to improve the convergence of the exact exchange with respect to the supercell size (or the k-point mesh density) 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, which can also be applied to screened Coulomb operators, are described below.

Auxiliary function

Probe-charge Ewald

Truncation

In this method the bare Coulomb operator is 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 Coulomb operators

and

have representations in the reciprocal space that are given by

and

respectively. Thus, 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

Related tags and articles

HFRCUT, Hybrid_functionals: formalism, Downsampling_of_the_Hartree-Fock_operator

References