The optimal P3M algorithm for computing electrostatic energies in periodic systems

We optimize Hockney and Eastwood's particle-particle particle-mesh algorithm to achieve maximal accuracy in the electrostatic energies (instead of forces) in three-dimensional periodic charged systems. To this end we construct an optimal influence function that minimizes the root-mean-square (rms) errors of the energies. As a by-product we derive a new real-space cutoff correction term, give a transparent derivation of the systematic errors in terms of Madelung energies, and provide an accurate analytical estimate for the rms error of the energies. This error estimate is a useful indicator of the accuracy of the computed energies and allows an easy and precise determination of the optimal values of the various parameters in the algorithm (Ewald splitting parameter, mesh size, and charge assignment order).