Accelerating molecular dynamics simulations using fast Ewald summation with prolates.
basic_science · Level V
Where this comes from
- Record sourced from PubMed, PMID 42168174.
- Also identified by DOI 10.1038/s41467-026-73232-8.
- No licence information is recorded for this record.
- Because redistribution is not established, this page shows the abstract only. Follow the links below for the full text.
Abstract
The evaluation of long-range Coulomb interactions is a significant cost in molecular dynamics (MD), even when using Particle Mesh Ewald (PME) or Particle-Particle-Particle-Mesh (PPPM) methods, which rely on Ewald splitting and the fast Fourier transform to achieve near-linear scaling. We introduce ESP-Ewald summation with prolate spheroidal wave functions (PSWFs)-which leads to a more efficient Fourier representation and a reduction in the required grid size, global communication, and particle-grid operations, without loss of accuracy. We have integrated the ESP method into two widely-used open-source MD packages, LAMMPS and GROMACS, enabling rapid comparison and adoption. Relative to PME/PPPM baselines at error tolerances 10<sup>-3</sup> to 10<sup>-4</sup>, ESP gives roughly a 3-fold acceleration of electrostatic interactions, and a 2.5-fold speed-up in the MD simulation when using about 10<sup>3</sup> compute cores. At high accuracy (10<sup>-5</sup>), these increase to 10-fold for the far-field electrostatics and 5-fold for MD simulation. Furthermore, we show that the accelerated codes have improved strong scaling with core count, and validate them in realistic long-time biological and material simulations. ESP thus offers a practical, drop-in path to reduce the time-to-solution and energy footprint of MD workflows.