Abstract:Reactive force field (ReaxFF), a recent and novel bond order potential, allows for reactive molecular dynamics (ReaxFF MD) simulations for modeling larger and more complex molecular systems involving chemical reactions when compared with computation intensive quantum mechanical methods. However, ReaxFF MD can be approximately 10-50 times slower than classical MD due to its explicit modeling of bond forming and breaking, the dynamic charge equilibration at each time-step, and its one order smaller time-step than the classical MD, all of which pose significant computational challenges in simulation capability to reach spatio-temporal scales of nanometers and nanoseconds. The very recent advances of graphics processing unit (GPU) provide not only highly favorable performance for GPU enabled MD programs compared with CPU implementations but also an opportunity to manage with the computing power and memory demanding nature imposed on computer hardware by ReaxFF MD. In this paper, we present the algorithms of GMD-Reax, the first GPU enabled ReaxFF MD program with significantly improved performance surpassing CPU implementations on desktop workstations. The performance of GMD-Reax has been benchmarked on a PC equipped with a NVIDIA C2050 GPU for coal pyrolysis simulation systems with atoms ranging from 1378 to 27,283. GMD-Reax achieved speedups as high as 12 times faster than Duin et al.'s FORTRAN codes in Lammps on 8 CPU cores and 6 times faster than the Lammps' C codes based on PuReMD in terms of the simulation time per time-step averaged over 100 steps. GMD-Reax could be used as a new and efficient computational tool for exploiting very complex molecular reactions via ReaxFF MD simulation on desktop workstations.

Several Optimized Algorithms for Ultra-soft Psedupotential DFT Molecular Dynamics Simulation

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

The Analysis of a Plane Wave Pseudopotential Density Functional Theory Code on a GPU Machine

GPU Tuning for First-Principle Electronic Structure Simulations

First Principle Calculation on Multi-GPU Machines

Large Scale Plane Wave Pseudopotential Density Functional Theory Calculations on GPU Clusters

Fast plane wave density functional theory molecular dynamics calculations on multi-GPU machines

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

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

An Efficient Multi-GPU Implementation for Linear-Response Time-Dependent Density Functional Theory

GPU_PBTE: an Efficient Solver for Three and Four Phonon Scattering Rates on Graphics Processing Units

First-principle nonlocal projector potential calculation on GPU cluster

High-speed and low-power molecular dynamics processing unit (MDPU) with ab initio accuracy

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

A GPU parallel scheme for accelerating 2D and 3D peridynamics models

DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization

GPGPU Supported Cooperative Acceleration in Molecular Dynamics

GPU Acceleration of Numerical Atomic Orbitals-Based Density Functional Theory Algorithms within the ABACUS package

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

Long Timestep Molecular Dynamics on the Graphical Processing Unit.