Abstract:Coulomb interaction, following an inverse-square force-law, quantifies the amount of force between two stationary and electrically charged particles. The long-range nature of Coulomb interactions poses a major challenge to molecular dynamics simulations, which are major tools for problems at the nano-/micro-scale. Various algorithms are developed to calculate the pairwise Coulomb interactions to a linear scale, but poor scalability limits the size of simulated systems. Here, we use an efficient molecular dynamics algorithm with the random batch Ewald method on all-atom systems where the complete Fourier components in the Coulomb interaction are replaced by randomly selected mini-batches. By simulating the N-body systems up to 108 particles using 10 000 central processing unit cores, we show that this algorithm furnishes O(N) complexity, almost perfect scalability, and an order of magnitude faster computational speed when compared to the existing state-of-the-art algorithms. Further examinations of our algorithm on distinct systems, including pure water, a micro-phase separated electrolyte, and a protein solution, demonstrate that the spatiotemporal information on all time and length scales investigated and thermodynamic quantities derived from our algorithm are in perfect agreement with those obtained from the existing algorithms. Therefore, our algorithm provides a promising solution on scalability of computing the Coulomb interaction. It is particularly useful and cost-effective to simulate ultra-large systems, which is either impossible or very costly to conduct using existing algorithms, and thus will be beneficial to a broad range of problems at nano-/micro-scales.

Ergodicity and long-time behavior of the Random Batch Method for interacting particle systems

On Random Batch Methods (RBM) for interacting particle systems driven by Lévy processes

Random Batch Methods for classical and quantum interacting particle systems and statistical samplings

On the Random Batch Method for second order interacting particle systems

Mean field error estimate of the random batch method for large interacting particle system

Convergence of Random Batch Method with replacement for interacting particle systems

Random Batch Method with Momentum Correction

Uniform Error Estimates for the Random Batch Method to the First‐order Consensus Models with Antisymmetric Interaction Kernels

Convergence of Random Batch Method for Interacting Particles with Disparate Species and Weights

Convergence of Random Batch Method for interacting particles with disparate species and weights

Reduced variance random batch methods for nonlocal PDEs

Energy Stable Scheme for Random Batch Molecular Dynamics

On the Mean Field Limit of the Random Batch Method for Interacting Particle Systems

A Random-Batch Monte Carlo Method for Many-Body Systems with Singular Kernels.

Some remarks on the effect of the Random Batch Method on phase transition

Ergodicity of Particle Systems

A RANDOM BATCH EWALD METHOD FOR PARTICLE SYSTEMS WITH COULOMB INTERACTIONS

A collision-oriented interacting particle system for Landau-type equations and the molecular chaos

Uniform-in-time Error Estimate of the Random Batch Method for the Cucker–Smale Model

LONG-TERM BEHAVIOUR OF LARGE MECHANICAL SYSTEMS WITH RANDOM INITIAL DATA