A Continuum,O(N) Monte-Carlo algorithm for charged particles
Joerg Rottler, A. C. Maggs
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 polyelectrolytes 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.
Create a lesson
Related papers
Charged Three-Dimensional Solitons in a Ferroelectric Nematic Liquid Crystal
Netra Prasad Dhakal, Zakaria Siddiquee, Hiroya Nishikawa et al.
Frank Elastic Constants of Associating Rigid-Rodlike Systems
Akihiko Matsuyama
Medial Model of Elastic Curve Contacts, from Clasps to Knits
Kathryn Chain, Haoye Sun, Michael S. Dimitriyev
How symmetry and complexity shape three-dimensional programmable self-assembly
Prarthana Agrawal, Joakim Bohlin, Ard A. Louis
Designing tunable hyperuniformity in random organization through nonlocal noise
Yusheng Lei, Yuhan Peng, Ran Ni
Long-time Brownian dynamics of hard and sticky colloids by geometric projection
Rémi Martin, Sandrine Duluard, Céline Merlet