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.