Abstract:Coulomb interaction, following an inverse-square force-law, quantifies the amount of force between two stationary, 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-/microscale. Various algorithms aim to speed up the pairwise Coulomb interactions to a linear scaling but the poor scalability limits the size of simulated systems. Here, we conduct 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 100 million particles using 10 thousand CPU 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, micro-phase-separated electrolyte and 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 breakthrough solution on scalability of computing the Coulomb interaction. It is particularly useful and cost-effective to simulate ultra-large systems, which was either impossible or very costing to conduct using existing algorithms, thus would benefit the broad community of sciences.

Achieving linear scaling for the electronic quantum Coulomb problem

Linear scaling density functional calculations with Gaussian orbitals

Achieving linear scaling in exchange-correlation density functional quadratures

Linear scaling electronic structure methods in chemistry and physics

Efficient Linear-Scaling Quantum Transport Calculations on Graphics Processing Units and Applications on Electron Transport in Graphene

Linear-scaling density-functional theory with Gaussian orbitals and periodic boundary conditions: Efficient evaluation of energy and forces via the fast multipole method

Linear scaling computation of the Fock matrix VII. Periodic Density Functional Theory at the $Γ$-point

Elongation Method: Towards Linear Scaling For Electronic Structure Of Random Polymers And Other Quasilinear Materials

Linear Scaling Algorithm for Tight-Binding Molecular Dynamics Simulations.

Linear scaling conjugate gradient density matrix search as an alternative to diagonalization for first principles electronic structure calculations

Electronic structure calculations and molecular dynamics simulations with linear system-size scaling

Linear scaling coupled cluster and perturbation theories in the atomic orbital basis

Linear Scaling Electronic Structure Calculations With Numerical Atomic Basis Set

Accurate ionic forces and geometry optimisation in linear scaling density-functional theory with local orbitals

A Sub-Linear Scaling Algorithm for Computing the Electronic Structure of Materials

O(N) methods in electronic structure calculations

Linear scaling DFT calculations for large Tungsten systems using an optimized local basis

A linear scaling method for Hartree–Fock exchange calculations of large molecules

Self-consistent field theory based molecular dynamics with linear system-size scaling

Super-Scalable Molecular Dynamics Algorithm

Efficient and Accurate Linear Algebraic Methods for Large-scale Electronic Structure Calculations with Non-orthogonal Atomic Orbitals