Ewald summation

Content sourced from Wikipedia, licensed under CC BY-SA 3.0.

Ewald summation is a practical way to calculate long-range interactions, like electrostatic forces, in systems that repeat periodically (such as crystals or simulations that use periodic boundary conditions). The core idea is to split the interaction into two parts that are easy to handle separately: a short-range piece that converges quickly in real space, and a long-range piece that converges quickly in Fourier (reciprocal) space.

In more detail, the long-range interaction is screened by a smooth function (usually a Gaussian). This turns the original problem into two sums: a real-space sum of the short-range part, and a reciprocal-space sum for the long-range part. The system is treated as infinitely periodic by tiling the central unit cell with images, and the total charge in the unit cell must be neutral to avoid mathematical divergences.

The total electrostatic energy is written as the sum of two contributions: E_sr, the short-range real-space part, and E_lr, the long-range reciprocal-space part. The short-range term is computed directly in real space where it decays quickly. The long-range term is handled in Fourier space using the charge density of the unit cell and the Fourier transform of the long-range kernel. When written in this way, the reciprocal-space sum can be made to converge very rapidly.

A widely used version is the particle-mesh Ewald (PME) method. Here, particle charges are spread onto a grid (a mesh), the Fourier transform is computed with fast Fourier transforms (FFTs), and the long-range energy is obtained from the grid-based quantities. PME scales as roughly O(N log N), making it especially suitable for large molecular dynamics simulations with periodic boundaries. The short-range part is still summed in real space, and both parts can be truncated with little loss of accuracy.

Some caveats: the unit cell must be charge-neutral to avoid infinite sums, and the way the long-range contributions are ordered can affect the result for certain polar systems (this is a feature of conditional convergence in infinite periodic models). Boundary conditions effectively introduce a surface charge that is accounted for in the formulation.

Ewald summation has a long history, originally developed in physics in 1921 to study ionic crystals, and it has since become a standard tool for calculating long-range forces in simulations of plasmas, galaxies, and molecules. The PME variant is a common, efficient implementation used in many modern simulations.