Extending electrostatics of dielectric spheres to arbitrary charge distributions with applications to biosystems

Arbitrary-Shape Dielectric Particles Interacting in the Linearized Poisson-Boltzmann Framework: An Analytical Treatment

This work considers the interaction of two dielectric particles of arbitrary shape immersed into a solvent containing a dissociated salt and assuming that the linearized Poisson-Boltzmann equation holds. We establish a new general spherical re-expansion result which relies neither on the conventional condition that particle radii are small with respect to the characteristic separating distance between particles nor on any symmetry assumption. This is instrumental in calculating suitable expansion coefficients for the electrostatic potential inside and outside the objects and in constructing small-parameter asymptotic expansions for the potential, the total electrostatic energy, and forces in ascending order of Debye screening. This generalizes a recent result for the case of dielectric spheres to particles of arbitrary shape and builds for the first time a rigorous (exact at the Debye-Hückel level) analytical theory of electrostatic interactions of such particles at arbitrary distances. Numerical tests confirm that the proposed theory may also become especially useful in developing a new class of grid-free, fast, highly scalable solvers.

Rigorous treatment of pairwise and many-body electrostatic interactions among dielectric spheres at the Debye-Hückel level

Electrostatic interactions among colloidal particles are often described using the venerable (two-particle) Derjaguin-Landau-Verwey-Overbeek (DLVO) approximation and its various modifications. However, until the recent development of a many-body theory exact at the Debye-Hückel level (Yu in Phys Rev E 102:052404, 2020), it was difficult to assess the errors of such approximations and impossible to assess the role of many-body effects. By applying the exact Debye-Hückel level theory, we quantify the errors inherent to DLVO and the additional errors associated with replacing many-particle interactions by the sum of pairwise interactions (even when the latter are calculated exactly). In particular, we show that: (1) the DLVO approximation does not provide sufficient accuracy at shorter distances, especially when there is an asymmetry in charges and/or sizes of interacting dielectric spheres; (2) the pairwise approximation leads to significant errors at shorter distances and at large and moderate Debye lengths and also gets worse with increasing asymmetry in the size of the spheres or magnitude or placement of the charges. We also demonstrate that asymmetric dielectric screening, i.e., the enhanced repulsion between charged dielectric bodies immersed in media with high dielectric constant, is preserved in the presence of free ions in the medium.

Charged dielectric spheres interacting in electrolytic solution: A linearized Poisson-Boltzmann equation model

We present an analytical theory of electrostatic interactions of two spherical dielectric particles of arbitrary radii and dielectric constants, immersed into a polarizable ionic solvent (assuming that the linearized Poisson-Boltzmann framework holds) and bearing arbitrary charge distributions expanded in multipolar terms. The presented development entails a novel two-center re-expansion analytical theory that expands upon and improves the existing ones, bypassing the conventional expansions in modified Bessel functions. On this basis, we develop a specific matrix formalism that facilitates the construction of asymptotic expansions in ascending order of Debye screening terms of potential coefficients, which are then employed to find exact closed-form expressions for the total electrostatic energy. In particular, this work allows us to explicitly and precisely quantify the k-screened terms of the potential coefficients and mutual interaction energy. Specific cases of monopolar and dipolar distributions are described in particular detail. Comprehensive numerical examples and tests of series convergence and the relative balance of leading and higher-order terms of the mutual interaction energy are presented depending on the inter-particle distance and particles' radii. The results of this work find application in soft matter modeling and, in particular, in computational biophysics and colloid science, where the availability of increasingly larger experimental structures at the atomic-level resolution makes numerical treatment challenging and calls for more efficient expressions and an increased range of validity.

Electrostatics of charged dielectric spheres with application to biological systems. III. Rigorous ionic screening at the Debye-Hückel level

The unequivocal role of electrostatic forces in biological (and colloidal) systems underscores the importance of attaining accurate and rapid calculations of electrostatic forces if one wishes to faithfully simulate the electrostatic aspect of a biological system. This paper makes significant progress toward this aspect as it rigorously incorporates ionic screening at the Debye-Hückel level for an electrolyte system containing dielectric spheres of finite radii. We investigated earlier this system without mobile ions via a surface charge method. However, the need for computing a large number of Wigner rotation matrix elements per configuration can significantly slow down the numerical calculations. This difficulty was recently overcome by our Wigner-matrix-free formalism. Unfortunately, in that method ions can only be included individually, making it impractical to investigate, for example, ionic screening in a system modeled by charged dielectric spheres immersed in a solution of mobile ions. Here, we overcome this difficulty by extending the surface charge method to treat ions implicitly. Previous treatments of charged dielectric spheres in a solution of mobile ions did not emphasize the energy reciprocity of electrostatics and are largely limited to a few spheres and/or special symmetries. Our new formalism respects reciprocity and accommodates arbitrarily many dielectric spheres of different dielectric constants and sizes while being rigorous at the Debye-Hückel level. The differences, and the relationship, between our new implicit ion treatment and our previous ion-free (or explicit ion) approach are described. A closed form for the electrostatic energy with implicit ions is also provided. This new formalism speeds up the computation of the electrostatic energy in the presence of ions, and accommodates permanent and induced multipoles that are very important when the polarization effect needs to be correctly included. We also mention how the proposed method can be transformed to a numerical method for use with arbitrary nonspherical surfaces.

Electrostatics of charged dielectric spheres with application to biological systems. II. A formalism bypassing Wigner rotation matrices

Given the crucial role of electrostatic forces in biological systems, accurate and rapid calculations of electrostatic forces are imperative in faithfully simulating biological systems. More than a decade ago, we proposed a surface charge method, applied it to a system of an arbitrary number of charged dielectric spheres, and obtained an exact solution for arbitrary configuration of the spheres. The precision depends only on the number of terms kept in a series expansion, and can therefore be controlled at will. However, the numerical implementation can be significantly slowed down by the need to compute a large number of Wigner rotation matrix elements each time the electrostatic energy is calculated. In this paper, we provide the proof of a formula introduced in 1992 and apply it to arrive at a useful closed form for the electrostatic interaction energy without computing any Wigner rotation matrix elements, hence significantly improving the efficiency for numeric implementation of the rigorous surface charge method. This new Wigner-matrix-free formalism may also be used to speed up the computation of the electrostatic energy in the presence of permanent and induced multipoles, which can be very important for atomic modeling with polarization effect included.

Surface charge method for molecular surfaces with curved areal elements I. Spherical triangles

Parametrizing a curved surface with flat triangles in electrostatics problems creates a diverging electric field. One way to avoid this is to have curved areal elements. However, charge density integration over curved patches appears difficult. This paper, dealing with spherical triangles, is the first in a series aiming to solve this problem. Here, we lay the ground work for employing curved patches for applying the surface charge method to electrostatics. We show analytically how one may control the accuracy by expanding in powers of the the arc length (multiplied by the curvature). To accommodate not extremely small curved areal elements, we have provided enough details to include higher order corrections that are needed for better accuracy when slightly larger surface elements are used.