Abstract:Molecular dynamics modeling has provided a powerful tool for simulating and understanding diverse systems – ranging from materials processes to biophysical phenomena. Parallel formulations of these methods have been shown to be among the most scalable scientific computing applications. Many instances of this class of methods rely on a static bond structure for molecules, rendering them infeasible for reactive systems. Recent work on reactive force fields has resulted in the development of ReaxFF, a novel bond order potential that bridges quantum-scale and classical MD approaches by explicitly modeling bond activity (reactions) and charge equilibration. These aspects of ReaxFF pose significant challenges from a computational standpoint, both in sequential and parallel contexts. Evolving bond structure requires efficient dynamic data structures. Minimizing electrostatic energy through charge equilibration requires the solution of a large sparse linear system with a shielded electrostatic kernel at each sub-femtosecond long time-step. In this context, reaching spatio-temporal scales of tens of nanometers and nanoseconds, where phenomena of interest can be observed, poses significant challenges.In this paper, we present the design and implementation details of the Purdue Reactive Molecular Dynamics code, PuReMD. PuReMD has been demonstrated to be highly efficient (in terms of processor performance) and scalable. It extends current spatio-temporal simulation capability for reactive atomistic systems by over an order of magnitude. It incorporates efficient dynamic data structures, algorithmic optimizations, and effective solvers to deliver low per-time-step simulation time, with a small memory footprint. PuReMD is comprehensively validated for performance and accuracy on up to 3375 cores on a commodity cluster (Hera at LLNL-OCF). Potential performance bottlenecks to scalability beyond our experiments have also been analyzed. PuReMD is available over the public domain and has been used to model diverse systems, ranging from strain relaxation in Si–Ge nanobars, water–silica surface interaction, and oxidative stress in lipid bilayers (bio-membranes).

Reactive Programming of Simulations in Physics

Simulation of Quantum Mechanics Using Reactive Programming

Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques

Reactive multi-particle collision dynamics with reactive boundary conditions

Synthesizing Functional Reactive Programs

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

Multi-physics simulations using a hierarchical interchangeable software interface

Tackling the Awkward Squad for Reactive Programming: The Actor-Reactor Model

PyRID: A Brownian dynamics simulator for reacting and interacting particles written in Python

Differentiable programming for particle physics simulations

Differentiable Programming of Reaction-Diffusion Patterns

VisualPDE: rapid interactive simulations of partial differential equations

Differentiable Physics: A Position Piece

Accelerating and parallelizing Lagrangian simulations of mixing-limited reactive transport

Automated Verification of Reactive and Concurrent Programs by Calculation

Green's Function Reaction Dynamics: a new approach to simulate biochemical networks at the particle level and in time and space

Continuous-space event-driven simulations of reaction-diffusion processes in three dimensions

Smoldyn on Graphics Processing Units: Massively Parallel Brownian Dynamics Simulations

Simulation of reactive flows using particle methods

Diva: A Declarative and Reactive Language for In-Situ Visualization