Long range forces in a performance portable Molecular Dynamics framework
This work addresses the need for efficient and portable MD simulations for physicists and chemists, but it is incremental as it extends an existing framework with specific algorithms.
The authors tackled the challenge of incorporating long-range electrostatic interactions into a performance-portable Molecular Dynamics framework, achieving competitive scaling and performance compared to existing codes through the implementation of an Ewald summation algorithm and OpenMP optimizations.
Molecular Dynamics (MD) codes predict the fundamental properties of matter by following the trajectories of a collection of interacting model particles. To exploit diverse modern manycore hardware, efficient codes must use all available parallelism. At the same time they need to be portable and easily extendible by the domain specialist (physicist/chemist) without detailed knowledge of this hardware. To address this challenge, we recently described a new Domain Specific Language (DSL) for the development of performance portable MD codes based on a "Separation of Concerns": a Python framework automatically generates efficient parallel code for a range of target architectures. Electrostatic interactions between charged particles are important in many physical systems and often dominate the runtime. Here we discuss the inclusion of long-range interaction algorithms in our code generation framework. These algorithms require global communications and careful consideration has to be given to any impact on parallel scalability. We implemented an Ewald summation algorithm for electrostatic forces, present scaling comparisons for different system sizes and compare to the performance of existing codes. We also report on further performance optimisations delivered with OpenMP shared memory parallelism.