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).

Pouring SIRAH on NAMD

The SIRAH force field 2.0: Altius, Fortius, Citius

Scalable molecular dynamics on CPU and GPU architectures with NAMD

Scalable molecular dynamics with NAMD

Development of the Quantum Inspired SIBFA Many-Body Polarizable Force Field: Enabling Condensed Phase Molecular Dynamics Simulations

New ways to boost molecular dynamics simulations

Quantum Fragment Based ab Initio Molecular Dynamics for Proteins.

Breaking the Million-Electron and 1 EFLOP/s Barriers: Biomolecular-Scale Ab Initio Molecular Dynamics Using MP2 Potentials

Breaking the mold: overcoming the time constraints of molecular dynamics on general-purpose hardware

Improved Adaptive Resolution Molecular Dynamics Simulation

Speeding-up Hybrid Functional based Ab Initio Molecular Dynamics using Multiple Time-stepping and Resonance Free Thermostat

Advancing Molecular Dynamics: Toward Standardization, Integration, and Data Accessibility in Structural Biology

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Density Fluctuations, Solvation Thermodynamics and Coexistence Curves in Grand Canonical Molecular Dynamics Simulations

Parametrization of a Reactive Force Field (ReaxFF) for Molecular Dynamics Simulations of Si Nanoparticles

Achieving an Order of Magnitude Speed-up in Hybrid Functional and Plane Wave based Ab Initio Molecular Dynamics: Applications to Proton Transfer Reactions in Enzymes and in Solution

VMD as a Platform for Interactive Small Molecule Preparation and Visualization in Quantum and Classical Simulations

Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques

An Event-Driven Hybrid Molecular Dynamics and Direct Simulation Monte Carlo Algorithm

Improved Parameterization of Amine-Carboxyate and Amine-Phosphate Interactions for Molecular Dynamics Simulations Using the CHARMM and AMBER Force Fields

Towards Electronic Structure-Based Ab-Initio Molecular Dynamics Simulations with Hundreds of Millions of Atoms