Abstract
Long ranged electrostatic interactions are time consuming to calculate in molecular dynamics and Monte Carlo simulations. We introduce an algorithmic framework for simulating charged particles which modifies the dynamics so as to allow equilibration using a local Hamiltonian. The method introduces an auxiliary field with constrained dynamics so that the equilibrium distribution is determined by the Coulomb interaction. We demonstrate the efficiency of the method by simulating a simple, charged lattice gas.
- Received 1 November 2001
DOI:https://doi.org/10.1103/PhysRevLett.88.196402
©2002 American Physical Society