Long Range Interactions

where are the lattice vectors . This sum is conditionally convergent; the result depends on the order the terms are summed in. The Ewald sum [127,115] is commonly used to model the effect of long range forces. It was originally developed for the study of ionic crystals. The Ewald sum decomposes the sum in (3.46) into two rapidly convergent sums

Physically the Ewald sum works by surrounding each point charge in the system by a charge distribution of equal magnitude and opposite sign. This is commonly taken to be a Gaussian distribution, although this choice is arbitrary. The counter-charge screens the original potential making it short ranged. This is then summed in real space. Then a second imaginary charge distribution of opposite sign to the first (and of the same sign as the point charges) is added to cancel out the screening charge. As this screening distribution is a smooth function, its Fourier transform is rapidly convergent. Thus this second is summed in reciprocal space. Two further terms are present in the Ewald sum. The first is the

The Ewald sum is possibly the most common method for evaluating long-range interactions in simulations. However, it can be expensive for large systems due to its scaling. Thus several, cheaper, alternative methods have been proposed. One such technique is the Particle-Particle Particle-Mesh method (PPPM) [128] which scales as . In this interactions at short range are directly calculated. Interactions at long range are evaluated by interpolating the charges onto a mesh. The Poisson equation can then be solved on this mesh to calculate the potential and forces. There are several variants of this technique that use different methods for discretizing the charges onto a mesh [129]. Another set of alternatives to the Ewald sum are the Fast Multipole Methods (FMM). These work by breaking the system into cells. The potential and force on each molecule are evaluated by considering the interaction between the molecule and each of these cells, where the distribution of charges in each of these cells are approximated by a multipole expansion. For a large enough system size this has been shown to scale as , although due to a large prefactor, these methods only become favourable for systems of about 100000 particles [130]. Another alternative to the Ewald sum is the reaction field method [131,115]. Here contributions from molecules within a cavity of radius are calculated explicitly. Molecules outside this are considered to form a dielectric continuum producing a reaction field within the cavity.