The fast Ewald methods are widely used to compute the point-charge electrostatic interactions in molecular simulations. The key step that introduces errors in the computation is the particle-mesh interpolation. In this work, the optimal interpolation basis is derived by minimizing the estimated error of the fast Ewald method. The basis can be either general or model specific, depending on whether or not the charge correlation has been taken into account. By using the TIP3P water as an example system, we demonstrate that the general optimal basis is always more accurate than the B-spline basis in the investigated parameter range, while the computational cost is at most 5% more expensive. In some cases, the optimal basis is found to be two orders of magnitude more accurate. The model specific optimal basis further improves the accuracy of the general optimal basis, but requires more computational effort in the optimization, and may not be transferable to systems with different charge correlations. Therefore, the choice between the general and model specific optimal bases is a trade-off between the generality and the accuracy.
REFERENCES
The evaluation of Kaiser-Bessel basis requires a square root and a hyperbolic sine function, which are usually much more expensive than the polynomials. It should be noted that in the productive codes, the Kaiser-Bessel basis is usually implemented by cubic interpolation of tabulated values; thus, each evaluation needs as many floating point operations as the optimal basis.
