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>  
 
To alleviate this issue and 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{{cite|gygi:prb:86}}, probe-charge Ewald {{cite|massidda:prb:93}} ({{TAG|HFALPHA}}), and Coulomb truncation methods{{cite|spenceralavi:prb:08}} ({{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.
To alleviate this issue and improve the convergence of the exact exchange integral with respect to supercell size (or k-point mesh density) different methods have been proposed: the auxiliary function methods{{cite|gygi:prb:86}}, probe-charge Ewald {{cite|massidda:prb:93}} ({{TAG|HFALPHA}}), and Coulomb truncation methods{{cite|spenceralavi:prb:08}} ({{TAG|HFRCUT}}).
These mostly involve modifying the Coulomb Kernel in a way that yields the same result as the unmodified kernel within the limit of large supercell sizes.

Revision as of 09:12, 10 May 2022

In the unscreened HF exchange, the bare Coulomb operator

is singular in the reciprocal space at :

To alleviate this issue and 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.