Improved random batch Ewald method in molecular dynamics simulations
arXiv:2204.13595 · doi:10.1021/acs.jpca.2c01918
Abstract
The random batch Ewald (RBE) is an efficient and accurate method for molecular dynamics (MD) simulations of physical systems at the nano-/micro- scale. The method shows great potential to solve the computational bottleneck of long-range interactions, motivating a necessity to accelerating short-range components of the non-bonded interactions for a further speedup of MD simulations. In this work, we present an improved RBE method for the non-bonding interactions by introducing the random batch idea to constructing neighbor lists for the treatment of both the short-range part of the Ewald splitting and the Lennard-Jones potential. The efficiency of the novel neighbor list algorithm owes to the stochastic minibatch strategy which can significantly reduce the total number of neighbors. We obtan the error estimate and convergence by theoretical analysis and implement the improved RBE method in the LAMMPS package. Benchmark simulations are performed to demonstrate the accuracy and stability of the algorithm. Numerical tests on computer performance by conducting large-scaled MD simulations for systems including up to 0.1 billion water molecules, run on massive cluster with up to 50 thousand CPU cores, demonstrating the attractive features such as the high parallel scalability and memory-saving of the method in comparison to the existing methods.
arXiv comment18 pages, 12 figures, to be published in J. Phys. Chem. A
References in corpus (13)
- A flexible algorithm for calculating pair interactions on SIMD architectures
- How to mesh up Ewald sums (II): An accurate error estimate for the P3M algorithm
- 86 PFLOPS Deep Potential Molecular Dynamics simulation of 100 million atoms with ab initio accuracy
- Local Simulation Algorithms for Coulomb Interaction
- Improved neighbor list algorithm in molecular simulations using cell decomposition and data sorting method
- On the Random Batch Method for second order interacting particle systems
- Onset of DNA Aggregation in Presence of Monovalent and Multivalent Counterions
- The -series: A separable decomposition for electrostatics computation with improved accuracy
- Efficiency of linked cell algorithms
- A random-batch Monte Carlo method for many-body systems with singular kernels
- A random batch Ewald method for charged particles in the isothermal-isobaric ensemble
- The Symmetry-Preserving Mean Field Condition for Electrostatic Correlations in Bulk
- HSMA: An O(N) electrostatics package implemented in LAMMPS