Efficient Algorithms for Electrostatic Interactions Including Dielectric Contrasts

Abstract

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.

Computing the Coulomb interaction in inhomogeneous dielectric media via a local electrostatics lattice algorithm

“…which leads to the more commonly known expression for equation (11). With these two last results, (13) and (14), two more Maxwell equations can be obtained by inserting into the constraint equation (5), namely Ampère's and Faraday's law:Ḋ…”

Lattice model of ionic liquid confined by metal electrodes

“…[58][59][60][61] Alternatively, there are electrostatic layer correction methods (ELC) or methods that explicitly sum infinite series of image charges. [62][63][64] Recently, we introduced a new simulation method which allows us to easily calculate sums over replicas. 65 The formalism is based on periodic Green's functions and was applied to a lattice model of a symmetric RTIL liquid; see Fig.…”

ESPResSo 4.0 – an extensible software package for simulating soft matter systems

“…ESPResSo has seen many additions since 2001, living up to this extensible character [48]. This includes support for the simulation of hydrodynamic interactions via an efficient GPU implementation of a lattice Boltzmann solver, rigid-body dynamics, a scheme for the study of irreversible agglomeration, and electrostatic solvers that allow for dielectric inclusions or even locally varying dielectric permittivities [48,[58][59][60][61][62].…”