Accelerated multiboson algorithm for Coulomb gases with dynamical dielectric effects
Abstract
A recent reformulation [Phys. Rev. E 73, 016705 (2006)] of the problem of Coulomb gases in the presence of a dynamical dielectric medium showed that finite-temperature simulations of such systems can be accomplished on the basis of completely local Hamiltonians on a spatial lattice by including additional bosonic fields. For large systems, the Monte Carlo algorithm proposed in the above-cited work becomes inefficient due to a low acceptance rate for particle moves in a fixed background multiboson field. We show here how this problem can be circumvented by use of a coupled particle-multiboson update procedure that improves acceptance rates on large lattices by orders of magnitude. The method is tested on a one-component plasma with neutral dielectric particles for a variety of system sizes.
- Publication:
-
Physical Review E
- Pub Date:
- June 2006
- DOI:
- 10.1103/PhysRevE.73.066711
- arXiv:
- arXiv:cond-mat/0602363
- Bibcode:
- 2006PhRvE..73f6711D
- Keywords:
-
- 05.10.-a;
- 61.20.Ja;
- 61.20.Qg;
- 05.50.+q;
- Computational methods in statistical physics and nonlinear dynamics;
- Computer simulation of liquid structure;
- Structure of associated liquids: electrolytes molten salts etc.;
- Lattice theory and statistics;
- Condensed Matter - Soft Condensed Matter;
- High Energy Physics - Lattice
- E-Print:
- 13 pages, 2 figures, fixed typos, added references