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

Abstract We present an efficient and highly accurate GP-GPU implementation of our community code, PuReMD, for reactive molecular dynamics simulations using the ReaxFF force field. PuReMD and its incorporation into LAMMPS (Reax/C) is used by a large number of research groups worldwide for simulating diverse systems ranging from biomembranes to explosives (RDX) at atomistic level of detail. The sub-femtosecond time-steps associated with ReaxFF strongly motivate significant improvements to per-timestep simulation time through effective use of GPUs. This paper presents, in detail, the design and implementation of PuReMD-GPU, which enables ReaxFF simulations on GPUs, as well as various performance optimization techniques we developed to obtain high performance on state-of-the-art hardware. Comprehensive experiments on model systems (bulk water and amorphous silica) are presented to quantify the performance improvements achieved by PuReMD-GPU and to verify its accuracy. In particular, our experiments show up to 16× improvement in runtime compared to our highly optimized CPU-only single-core ReaxFF implementation. PuReMD-GPU is a unique production code, and is currently available on request from the authors.

[1]  Markus J. Buehler,et al.  Hierarchical chemo-nanomechanics of proteins: entropic elasticity, protein unfolding and molecular fracture , 2007 .

[2]  Emppu Salonen,et al.  Polarizable force fields. , 2013, Methods in molecular biology.

[3]  S. Stuart,et al.  A reactive potential for hydrocarbons with intermolecular interactions , 2000 .

[4]  John L. Klepeis,et al.  Anton, a special-purpose machine for molecular dynamics simulation , 2007, ISCA '07.

[5]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[6]  Ananth Grama,et al.  Reactive Molecular Dynamics: Numerical Methods and Algorithmic Techniques , 2012, SIAM J. Sci. Comput..

[7]  A. V. van Duin,et al.  Simulations on the thermal decomposition of a poly(dimethylsiloxane) polymer using the ReaxFF reactive force field. , 2005, Journal of the American Chemical Society.

[8]  D. Brenner,et al.  Empirical potential for hydrocarbons for use in simulating the chemical vapor deposition of diamond films. , 1990, Physical review. B, Condensed matter.

[9]  Ananth Grama,et al.  Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques , 2012, Parallel Comput..

[10]  A. Thompson,et al.  Atomistic simulation of orientation dependence in shock-induced initiation of pentaerythritol tetranitrate. , 2013, The journal of physical chemistry. B.

[11]  Hasan Metin Aktulga PRISM Seminar Series - Spring 2010 A FAST ILU PRECONDITIONING-BASED SOLVER FOR THE CHARGE EQUILIBRATION PROBLEM , 2010 .

[12]  Joseph C. Fogarty,et al.  A reactive molecular dynamics simulation of the silica-water interface. , 2010, The Journal of chemical physics.

[13]  Ananth Grama,et al.  Algorithmic and numerical techniques for atomistic modeling , 2010 .

[14]  A. V. van Duin,et al.  Development of the ReaxFF reactive force field for describing transition metal catalyzed reactions, with application to the initial stages of the catalytic formation of carbon nanotubes. , 2005, The journal of physical chemistry. A.

[15]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[16]  Ashish Sharma,et al.  De Novo Ultrascale Atomistic Simulations On High-End Parallel Supercomputers , 2008, Int. J. High Perform. Comput. Appl..

[17]  Alex Fit-Florea,et al.  Precision and Performance: Floating Point and IEEE 754 Compliance for NVIDIA GPUs , 2011 .

[18]  Li Guo,et al.  Algorithms of GPU-enabled reactive force field (ReaxFF) molecular dynamics. , 2013, Journal of molecular graphics & modelling.

[20]  Michael Garland,et al.  Efficient Sparse Matrix-Vector Multiplication on CUDA , 2008 .

[21]  Donald W. Brenner,et al.  A second-generation reactive empirical bond order (REBO) potential energy expression for hydrocarbons , 2002 .

[22]  Aiichiro Nakano,et al.  Parallel multilevel preconditioned conjugate-gradient approach to variable-charge molecular dynamics , 1997 .

[23]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[24]  A. V. Duin,et al.  ReaxFFSiO Reactive Force Field for Silicon and Silicon Oxide Systems , 2003 .

[25]  A. V. Duin,et al.  ReaxFF: A Reactive Force Field for Hydrocarbons , 2001 .

[26]  Ananth Grama,et al.  Strain relaxation in Si/Ge/Si nanoscale bars from molecular dynamics simulations , 2009 .

[27]  Sandeep Patel,et al.  Revised charge equilibration potential for liquid alkanes. , 2008, The journal of physical chemistry. B.