Coulomb singularity: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
Line 20: Line 20:
V(0)=2\pi R_{\text{c}}^{2}
V(0)=2\pi R_{\text{c}}^{2}
</math>
</math>
The screened potentials
The screened Coulomb operators
:<math>
:<math>
V(\vert\mathbf{r}-\mathbf{r}'\vert)=\frac{e^{-\lambda\left\vert\mathbf{r}-\mathbf{r}'\right\vert}}{\left\vert\mathbf{r}-\mathbf{r}'\right\vert}
V(\vert\mathbf{r}-\mathbf{r}'\vert)=\frac{e^{-\lambda\left\vert\mathbf{r}-\mathbf{r}'\right\vert}}{\left\vert\mathbf{r}-\mathbf{r}'\right\vert}
Line 36: Line 36:
V(q)=\frac{4\pi}{q^{2}}\left(1-e^{-q^{2}/\left(4\lambda^2\right)}\right)
V(q)=\frac{4\pi}{q^{2}}\left(1-e^{-q^{2}/\left(4\lambda^2\right)}\right)
</math>
</math>
respectively. Thus, the screened potentials have no singularity at <math>q=0</math>.
respectively. Thus, the screened potentials have no singularity at <math>q=0</math>. Nevertheless, it is still beneficial for accelerating convergence with respect to the number of k-points to multiply these screened operators by <math>\theta(R_{\text{c}}-\left\vert\mathbf{r}-\mathbf{r}'\right\vert)</math>.




=== Auxiliary function methods ===
=== Auxiliary function methods ===

Revision as of 11:47, 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 a singularity at , and 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 methods[1], probe-charge Ewald [2] (HFALPHA), and Coulomb truncation methods[3] (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 are described below.

Truncation methods

The potential 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

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 convergence with respect to the number of k-points to multiply these screened operators by .


Auxiliary function methods