Local Simulation Algorithms for Coulomb Gases with Dynamical Dielectric   Effects
  We discuss the application of the local lattice technique of Maggs and Rossetto to problems that involve the motion of objects with different dielectric constants than the background. In these systems the simulation method produces a spurious interaction force which causes the particles to move in an unphysical manner. We show that this term can be removed using a variant of a method known from high-energy physics simulations, the multiboson method, and demonstrate the effectiveness of this corrective method on a system of neutral particles. We then apply our method to a one-component plasma to show the effect of the spurious interaction term on a charged system. 
A. Duncan, R.D. Sedgewick, R.D. Coalson