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

Optimization of Reactive Force Field Simulation: Refactor, Parallelization, and Vectorization for Interactions

Millimeter-Scale and Billion-Atom Reactive Force Field Simulation on Sunway Taihulight

Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques

Parameterization of a reactive force field using a Monte Carlo algorithm

An Initial Design-enhanced Deep Learning-based Optimization Framework to Parameterize Multicomponent ReaxFF Force Fields

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

PhyNEO: A Neural-Network-Enhanced Physics-Driven Force Field Development Workflow for Bulk Organic Molecule and Polymer Simulations.

Optimizing Force Fields with Experimental Data Using Ensemble Reweighting and Potential Contrasting

Development of a New Parameter Optimization Scheme for a Reactive Force Field (ReaxFF) Based on a Machine Learning Approach

Development of a new parameter optimization scheme for a reactive force field based on a machine learning approach

INDEEDopt: a deep learning-based ReaxFF parameterization framework

Efficient and Precise Force Field Optimization for Biomolecules Using DPA-2

JAX-ReaxFF: A Gradient-Based Framework for Fast Optimization of Reactive Force Fields.

Generating High-Precision Force Fields for Molecular Dynamics Simulations to Study Chemical Reaction Mechanisms using Molecular Configuration Transformer

General Multiobjective Force Field Optimization Framework, with Application to Reactive Force Fields for Silicon Carbide

ReaxFF Parameter Optimization with Monte-Carlo and Evolutionary Algorithms: Guidelines and Insights

A Neural-Network-Based Mapping and Optimization Framework for High-Precision Coarse-Grained Simulation

Accurate interatomic force field for molecular dynamics simulation by hybridizing classical and machine learning potentials

Global optimization of parameters in the reactive force field ReaxFF for SiOH

An Efficient Reactive Force Field without Explicit Coordination Dependence for Studying Caustic Aluminum Chemistry

Implementing Reactivity in Molecular Dynamics Simulations with Harmonic Force Fields