Local Simulation Algorithms for Coulombic Interactions

Abstract

We consider dynamically constrained Monte-Carlo dynamics and show that this leads to the generation of long ranged effective interactions. This allows us to construct a local algorithm for the simulation of charged systems without ever having to evaluate pair potentials or solve the Poisson equation. We discuss a simple implementation of a charged lattice gas as well as more elaborate off-lattice versions of the algorithm. There are analogies between our formulation of electrostatics and the bosonic Hubbard model in the phase approximation. Cluster methods developed for this model further improve the efficiency of the electrostatics algorithm.