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

RevelsMD: Reduced Variance Estimators of the Local Structure in Molecular Dynamics

Use the force! Reduced variance estimators for densities, radial distribution functions and local mobilities in molecular simulations

Learning Effective Molecular Models from Experimental Observables.

drMD: Molecular Dynamics for Experimentalists

MDAnalysis: A toolkit for the analysis of molecular dynamics simulations

Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques

Learning the Quantum Centroid Force Correction in Molecular Systems: A Localized Approach

MDRefine: a Python package for refining Molecular Dynamics trajectories with experimental data

Learning reduced kinetic Monte Carlo models of complex chemistry from molecular dynamics

Rigid Base Biasing in Molecular Dynamics Enables Enhanced Sampling of DNA Conformations

Robust Automated Equilibration Detection for Molecular Simulations

MDSCAN: RMSD-based HDBSCAN clustering of long molecular dynamics

Reaction analysis and visualization of ReaxFF molecular dynamics simulations

Statistical mechanical refinement of protein structure prediction schemes. II. Mayer cluster expansion approach

RELION: Implementation of a Bayesian approach to cryo-EM structure determination

DiffMD: A Geometric Diffusion Model for Molecular Dynamics Simulations

MD-TASK: a software suite for analyzing molecular dynamics trajectories

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

A local resampling trick for focused molecular dynamics

Enhanced sampling of robust molecular datasets with uncertainty-based collective variables

MolSieve: A Progressive Visual Analytics System for Molecular Dynamics Simulations