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

PuReMD-GPU: A reactive molecular dynamics simulation package for GPUs

Algorithms of GPU-enabled reactive force field (ReaxFF) molecular dynamics

Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques

Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

Implementation of the moving particle semi-implicit method for free-surface flows on GPU clusters.

86 PFLOPS Deep Potential Molecular Dynamics simulation of 100 million atoms with ab initio accuracy

Graphics Processing Unit Acceleration and Parallelization of GENESIS for Large-Scale Molecular Dynamics Simulations

Molecular dynamics simulation of macromolecules using graphics processing unit

RBMD: A molecular dynamics package enabling to simulate 10 million all-atom particles in a single graphics processing unit

Molecular dynamics simulations with many-body potentials on multiple GPUs - The implementation, package and performance.

Robust and effective ab initio molecular dynamics simulations on the GPU cloud infrastructure using the Schrödinger Materials Science Suite

Optimizing the Performance of Reactive Molecular Dynamics Simulations for Many-Core Architectures

GPGPU Supported Cooperative Acceleration in Molecular Dynamics

A GPU accelerated mixed-precision Smoothed Particle Hydrodynamics framework with cell-based relative coordinates

GPU acceleration of ab initio simulations of large-scale identical particles based on path integral molecular dynamics

Efficient molecular dynamics simulations with many-body potentials on graphics processing units

Langevin dynamics simulations of biomolecules on graphics processors

Multi-GPU RI-HF Energies and Analytic Gradients $-$ Towards High Throughput Ab Initio Molecular Dynamics

Quantum Chemistry for Solvated Molecules on Graphical Processing Units (GPUs)using Polarizable Continuum Models