Coulomb interactions in a computer simulation of a system periodic in two directions

Abstract

The integral representation of the gamma function and the Poisson summation formula are used to calculate the interaction energy of charged particles in a 3-dimensional system periodic in two directions. A parallelogram shape simulation box is considered. Calculations are carried out for interactions described by any inverse power, and analytical continuation of the energy function leads to the final expression for the Coulomb interaction energy. Summation over the simulation box replica along one or the other side of the box base is replaced by summation in reciprocal space. Therefore there are two equivalent formulas for the potential energy that offer the possibility of avoiding slowly convergent series. The energy expressions are identical to those obtained from the Lekner method. The special case is considered where the functions defining the energy are infinite, i.e. when two charges lie on a line parallel to the simulation box side that was chosen to convert real space summation into reciprocal space.