Electrostatic interaction in the presence of dielectric interfaces and polarization-induced like-charge attraction

Abstract: Electrostatic polarization is important in many nano-/micro-scale physical systems such as colloidal suspensions, biopolymers, and nanomaterials assembly. The calculation of polarization potential requires an efficient algorithm for solving 3D Poisson's equation. We have developed a useful image charge method to rapid evaluation of the Green's function of the Poisson's equation in the presence of spherical dielectric discontinuities. This paper presents an extensive study of this method by giving an convergence analysis and developing a coarse-graining algorithm. The use of the coarse graining could reduce the number of image charges to around a dozen, by 1-2 orders of magnitude. We use the algorithm to investigate the interaction force between likely charged spheres in different dielectric environments. We find the size and charge asymmetry leads to an attraction between like charges, in agreement with existing results. Furthermore, we study three-body interaction and find in the presence of an external interface, the interaction force depends on the curvature of the interface, and behavior a non-monotonic electrostatic force.