Extending the fast multipole method for charges inside a dielectric sphere in an ionic solvent: High-order image approximations for reaction fields

As a sequel to our previous paper on extending the Fast Multipole Method (FMM) for charges inside a dielectric sphere [J. Comput. Phys. 223 (2007) 846–864], this paper further extends the FMM to the electrostatic calculation for charges inside a dielectric sphere immersed in an ionic solvent, a scenery with more relevance in biological applications. The key findings include two fourth-order multiple discrete image approximations in terms of u = λa to the reaction field induced by the ionic solvent, provided that u = λa < 1 where λ is the inverse Debye screening length of the ionic solvent and a is the radius of the dielectric sphere. A 10−4 relative accuracy in the reaction field of a source charge within the sphere can be achieved with only 3–4 point image charges. Together with the image charges, the FMM can be used to speed up the calculation of electrostatic interactions of charges in a dielectric sphere immersed in an ionic solvent.