No abstract

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.

Coarse-grained models of soft matter are usually combined with implicit solvent models that take the electrostatic polarizability into account via a dielectric background. In biophysical or nanoscale simulations that include water, this constant can vary greatly within the system. Performing molecular dynamics or other simulations that need to compute exact electrostatic interactions between charges in those systems is computationally demanding. We review here several algorithms developed by us that perform exactly this task. For planar dielectric surfaces in partial periodic boundary conditions, the arising image charges can be either treated with the MMM2D algorithm in a very efficient and accurate way or with the electrostatic layer correction term, which enables the user to use his favorite 3D periodic Coulomb solver. Arbitrarily-shaped interfaces can be dealt with using induced surface charges with the induced charge calculation (ICC*) algorithm. Finally, the local electrostatics algorithm, MEMD (Maxwell Equations Molecular Dynamics), even allows one to employ a smoothly varying dielectric constant in the systems. We introduce the concepts of these three algorithms and an extension for the inclusion of boundaries that are to be held fixed at a constant potential (metal conditions). For each method, we present a showcase application to highlight the importance of dielectric interfaces.

Dissolved ions can alter the local permittivity of water; nevertheless most theories and simulations ignore this fact. We present a novel algorithm for treating spatial and temporal variations in the permittivity and use it to measure the equivalent conductivity of a salt-free polyelectrolyte solution. Our new approach quantitatively reproduces experimental results unlike simulations with a constant permittivity that even qualitatively fail to describe the data. We can relate this success to a change in the ion distribution close to the polymer due to the buildup of a permittivity gradient.

The ion distribution around charged colloids in solution has been investigated intensely during the last decade. However, few theoretical approaches have included the influence of variation in the dielectric permittivity within the system, let alone in the surrounding solvent. In this article, we introduce two relatively new methods that can solve the Poisson equation for systems with varying permittivity. The harmonic interpolation method (HIM) approximately solves the Green's function in terms of a spherical harmonics series, and thus provides analytical ion-ion potentials for the Hamiltonian of charged systems. The Maxwell Equations Molecular Dynamics (MEMD) algorithm features a local approach to electrostatics, allowing for arbitrary local changes of the dielectric constant. We show that the results of both methods are in very good agreement. We also found that the renormalized charge of the colloid, and with it the effective far field interaction, significantly changes if the dielectric properties within the vicinity of the colloid are changed.