A continuum, O(N) Monte Carlo algorithm for charged particles
Abstract
We introduce a Monte Carlo algorithm for the simulation of charged particles moving in the continuum. Electrostatic interactions are not instantaneous as in conventional approaches, but are mediated by a constrained, diffusing electric field on an interpolating lattice. We discuss the theoretical justifications of the algorithm and show that it efficiently equilibrates model electrolytes and polar fluids. In order to reduce lattice artifacts that arise from the interpolation of charges to the grid we implement a local, dynamic subtraction algorithm. This dynamic scheme is completely general and can also be used with other Coulomb codes, such as multigrid based methods.
- Publication:
-
Journal of Chemical Physics
- Pub Date:
- February 2004
- DOI:
- 10.1063/1.1642590
- arXiv:
- arXiv:cond-mat/0308441
- Bibcode:
- 2004JChPh.120.3119R
- Keywords:
-
- 61.20.Gy;
- 82.45.Gj;
- 02.70.Tt;
- 02.30.Hq;
- Theory and models of liquid structure;
- Electrolytes;
- Justifications or modifications of Monte Carlo methods;
- Ordinary differential equations;
- Condensed Matter - Soft Condensed Matter;
- Condensed Matter - Materials Science;
- Physics - Chemical Physics;
- Physics - Computational Physics
- E-Print:
- doi:10.1063/1.1642590