Communications: The Metropolis Monte Carlo finite element algorithm for electrostatic interactions.
The Metropolis Monte Carlo algorithm with the finite element method applied to compute electrostatic interaction energy between charge densities is described in this work. By using the finite element method to integrate numerically Poisson's equation, it is shown that the computing time to obtain the acceptance probability of an elementary trial move does not, in principle, depend on the number of charged particles present in the system.
Mazars
Martial
M
Laboratoire de Physique Théorique (UMR 8627), Université Paris Sud 11 and CNRS, Bâtiment 210, Orsay Cedex 91405, France. martial.mazars@th.u-psud.fr