A random batch Ewald method for particle systems with Coulomb interactions
arXiv:2010.01559
Abstract
We develop a random batch Ewald (RBE) method for molecular dynamics simulations of particle systems with long-range Coulomb interactions, which achieves an complexity in each step of simulating the -body systems. The RBE method is based on the Ewald splitting for the Coulomb kernel with a random "mini-batch" type technique introduced to speed up the summation of the Fourier series for the long-range part of the splitting. Importance sampling is employed to reduce the induced force variance by taking advantage of the fast decay property of the Fourier coefficients. The stochastic approximation is unbiased with controlled variance. Analysis for bounded force fields gives some theoretic support of the method. Simulations of two typical problems of charged systems are presented to illustrate the accuracy and efficiency of the RBE method in comparison to the results from the Debye-Hückel theory and the classical Ewald summation, demonstrating that the proposed method has the attractiveness of being easy to implement with the linear scaling and is promising for many practical applications.
arXiv comment23 pages, 7 figures
References in corpus (8)
- Electrostatic correlations: from Plasma to Biology
- How to mesh up Ewald sums (II): An accurate error estimate for the P3M algorithm
- Direct observation of charge inversion by multivalent ions as a universal electrostatic phenomenon
- On the Random Batch Method for second order interacting particle systems
- Charge inversion at minute electrolyte concentrations
- Simulation of charge reversal in salty environments: Giant overcharging?
- A stochastic version of Stein Variational Gradient Descent for efficient sampling
- A fast spectral method for electrostatics in doubly-periodic slit channels