but the long-range nature of the Coulomb interactions in direct space poses problems for the evaluation of multicenter integrals such as those appearing in HF or MP2 correlation energy . Indeed, discrete sums of the type of Eq.7 with Φ(G) = Φ(G) are facing a singularity in G = 0 which is only properly integrable in the thermodynamic limit (Ω → ∞, ∑ G → Ω/(2π) 3 dG) and is a consequence of the finite simulation cell imposed by the numerical computation. Simply ignoring the problematic component makes the convergence of the integrals very slow and requires either many replicas of the unit cell (large supercell), or a much finer and careful sampling of the eventual k-points mesh. Therefore, schemes have been suggested in order to screen the G = 0 divergence and obtain a faster convergence .