Local simulation algorithms for Coulomb gases with dynamical dielectric effects

Abstract

We discuss the application of the local lattice technique of Maggs and Rossetto [Phys. Rev. Lett. 88, 196402 (2002)] 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.