  We introduce a new method for simulating colloidal suspensions with spherical
colloidal particles of dielectric constant different from the surrounding
medium. The method uses exact calculation of the Green function to obtain the
ion-ion interaction potential in the presence of a dielectric discontinuity at
the surface of the colloidal particle. The new method is orders of magnitude
faster than the traditional approaches based on series expansions of the
interaction potential.
