The first term is calculated on a global grid using FFT techniques and the remaining two contributions are computed numerically using atomic grids. Due to the non-local character of the Coulomb operator, the decomposition is however more complex for the electrostatic energy, requiring to introduce atom-dependent screening densities n 0 A ,