Parameterization of Coulomb interaction in three-dimensional periodic systems

Abstract

A simple method is proposed to calculate Coulomb interactions in three-dimensional periodic cubic systems. It is based on the parameterization of the interaction on polynomials and rational functions. The parameterized functions are compared to tabulation methods, to the Ewald calculations and cubic harmonic function fits found in the literature. Our parameterizations are computationally more efficient than, the use of tabulations at all cases and seem to be more efficient than the cubic harmonic parameterizations in the case of simultaneous potential energy and force calculations. In comparison to the Ewald method, it is feasible to use the parameterizations on small systems and on systems, where pair-wise additive short-range interactions are dominant. One also may prefer the parameterizations to the Ewald method for large systems, if limited accuracy is needed. The embedding of the method into existing molecular dynamics and Monte Carlo simulation codes is very simple. The presented investigation contains some numerical experimental data to support the correct theoretical partition of potential energy in periodic systems, as well.