We set up and solve the problem of optimally partitioning the Coulomb
operator 1/r into a sum of two functions f(1)(r) and f(2)(r) such that
both f(1) and the Fourier transform of f(2) decay as quickly as possi
ble. The rigorous solution involves a Hermite function, but we find th
at the conventional Ewald-KWIK partition appears to be only slightly i
nferior.