Corrected Article: “An error-controlled fast multipole method” [J. Chem. Phys.131, 244102 (2009)]

Anisotropic reaction field correction for long-range electrostatic interactions in molecular dynamics simulations

Reaction-field (RF) methods have been extensively used in molecular dynamics simulations to efficiently compute long-range electrostatic interactions. They assume a continuous dielectric medium outside a certain cutoff, which has shown to be a reasonable approximation in many cases. However, lattice sum or fast multipole methods are nowadays often used instead, which treat long-range interactions explicitly but may introduce different artefacts. In the following work, the major issue of RFs is addressed, i.e., their inability to account for inhomogeneity even in heterogenous environments (e.g., membranes or protein binding pockets). By using a first-order Laplace series expansion of the dielectric permittivity on the cutoff sphere, local anisotropic effects can be described in a simple form. It is shown that the resulting boundary-value problem cannot be solved analytically, but instead a well-behaved approximative anisotropic reaction field (ARF) is introduced, which preserves coordinate invariance and approaches the standard RF solution for homogeneous systems. The comparison of RF to the state-of-the-art particle-particle particle-mesh (P₃M) method shows a difference in the orientation of molecules close to the interface between two different dielectrics (water/chloroform). It was found that the ARF leads to a correction of the orientational distribution toward the P₃M reference for planar and spherical interfaces.

The grid-based fast multipole method – a massively parallel numerical scheme for calculating two-electron interaction energies

A grid-based fast multipole method has been developed for calculating two-electron interaction energies for non-overlapping charge densities.

Comparison of scalable fast methods for long-range interactions

Based on a parallel scalable library for Coulomb interactions in particle systems, a comparison between the fast multipole method (FMM), multigrid-based methods, fast Fourier transform (FFT)-based methods, and a Maxwell solver is provided for the case of three-dimensional periodic boundary conditions. These methods are directly compared with respect to complexity, scalability, performance, and accuracy. To ensure comparable conditions for all methods and to cover typical applications, we tested all methods on the same set of computers using identical benchmark systems. Our findings suggest that, depending on system size and desired accuracy, the FMM- and FFT-based methods are most efficient in performance and stability.

Strongly coupled plasmas via Rydberg blockade of cold atoms

We propose and analyze a new scheme to produce ultracold neutral plasmas deep in the strongly coupled regime. The method exploits the interaction blockade between cold atoms excited to high-lying Rydberg states and therefore does not require substantial extensions of current ultracold plasma experiments. Extensive simulations reveal a universal behavior of the resulting Coulomb coupling parameter, providing a direct connection between the physics of strongly correlated Rydberg gases and ultracold plasmas. The approach is shown to reduce currently accessible temperatures by more than an order of magnitude, which opens up a new regime for ultracold plasma research and cold ion-beam applications with readily available experimental techniques.