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
Competing routes to spontaneous flow in confined active nematics
Rahil N. Valani, Vedad Dzanic, Sumesh P. Thampi et al.
Scaling and Condensation of Dry Active Matter Around Circular Obstacles
Felipe P. S. Júnior, F. Q. Potiguar, Jorge L. C. Domingos et al.
Active Hydrodynamics Couples Polymer Organization, Shape Fluctuations, and Motility in Deformable Droplets
Ritu Raj, P. B. Sunil Kumar
Inferring interactions between active particles using harmonic traps
Arnaud Compagnie, Joscha Mecke, Ivo Buttinoni et al.
Spontaneous filament formation and network self-assembly via active phase separation
Elena Lucas, Varun Venkatesh, Amin Doostmohammadi
Sensitivity of Nucleation Thermodynamics and Kinetics to the Treatment of Long-Range Interactions
Fernanda Sulantay Vargas, Kimia Sinaeian, Amir Haji-Akbari