Rigid body constraints realized in massively-parallel molecular dynamics on graphics processing units

Abstract Molecular dynamics (MD) methods compute the trajectory of a system of point particles in response to a potential function by numerically integrating Newtonʼs equations of motion. Extending these basic methods with rigid body constraints enables composite particles with complex shapes such as anisotropic nanoparticles, grains, molecules, and rigid proteins to be modeled. Rigid body constraints are added to the GPU-accelerated MD package, HOOMD-blue, version 0.10.0. The software can now simulate systems of particles, rigid bodies, or mixed systems in microcanonical (NVE), canonical (NVT), and isothermal-isobaric (NPT) ensembles. It can also apply the FIRE energy minimization technique to these systems. In this paper, we detail the massively parallel scheme that implements these algorithms and discuss how our design is tuned for the maximum possible performance. Two different case studies are included to demonstrate the performance attained, patchy spheres and tethered nanorods. In typical cases, HOOMD-blue on a single GTX 480 executes 2.5–3.6 times faster than LAMMPS executing the same simulation on any number of CPU cores in parallel. Simulations with rigid bodies may now be run with larger systems and for longer time scales on a single workstation than was previously even possible on large clusters.

[1]  Klaus Schulten,et al.  Accelerating Molecular Modeling Applications with GPU Computing , 2009 .

[2]  Vijay S Pande,et al.  CCMA: A Robust, Parallelizable Constraint Method for Molecular Simulations. , 2010, Journal of chemical theory and computation.

[3]  T. Mexia,et al.  Author ' s personal copy , 2009 .

[4]  Sharon C Glotzer,et al.  Switchable helical structures formed by the hierarchical self-assembly of laterally tethered nanorods. , 2009, Small.

[5]  V. Pande,et al.  Absolute comparison of simulated and experimental protein-folding dynamics , 2002, Nature.

[6]  L. Broadbelt,et al.  Distance and angular holonomic constraints in molecular simulations. , 2010, The Journal of chemical physics.

[7]  John E. Stone,et al.  An efficient library for parallel ray tracing and animation , 1998 .

[8]  J. Elliott,et al.  A dissipative particle dynamics method for modeling the geometrical packing of filler particles in polymer composites , 2000 .

[9]  Thomas F. Miller,et al.  Symplectic quaternion scheme for biophysical molecular dynamics , 2002 .

[10]  S. Glotzer,et al.  Self-Assembly of Patchy Particles. , 2004, Nano letters.

[11]  O. Urakawa,et al.  Small - , 2007 .

[12]  R. Rosenfeld.  Nature , 2009, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

[13]  M. P. Neal,et al.  Time reversible and symplectic integrators for molecular dynamics simulations of rigid molecules. , 2005, The Journal of chemical physics.

[14]  R. Mountain.  Solvation Structure of Ions in Water , 2007 .

[15]  Joshua A. Anderson,et al.  General purpose molecular dynamics simulations fully implemented on graphics processing units , 2008, J. Comput. Phys..

[16]  P. Meijer,et al.  Simulation of C60 through the plastic transition temperatures , 1995 .

[17]  H. R. Noori.  on "Understanding molecular simulations" , 2013 .

[18]  Sharon C Glotzer,et al.  Molecular simulation study of self-assembly of tethered V-shaped nanoparticles. , 2008, The Journal of chemical physics.

[19]  Matt K. Petersen,et al.  Effect of particle shape and charge on bulk rheology of nanoparticle suspensions , 2010, 1004.2411.

[20]  William H. Press,et al.  Numerical recipes in C , 2002 .

[21]  Klaus Schulten,et al.  Challenges in protein-folding simulations , 2010 .

[22]  A. Arnold,et al.  Harvesting graphics power for MD simulations , 2007, 0709.3225.

[23]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[24]  F E Cohen,et al.  Protein conformational landscapes: Energy minimization and clustering of a long molecular dynamics trajectory , 1995, Proteins.

[25]  Pekka Koskinen,et al.  Structural relaxation made simple. , 2006, Physical review letters.

[26]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

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

[28]  S. Glotzer,et al.  Stability of the double gyroid phase to nanoparticle polydispersity in polymer-tethered nanosphere systems , 2009 .

[29]  S. Glotzer,et al.  Reconfigurable assemblies of shape-changing nanorods. , 2010, ACS nano.

[30]  P. Kollman,et al.  Pathways to a protein folding intermediate observed in a 1-microsecond simulation in aqueous solution. , 1998, Science.

[31]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[32]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[33]  Sohail Murad,et al.  Singularity free algorithm for molecular dynamics simulation of rigid polyatomics , 1977 .

[34]  S. Glotzer,et al.  Self-assembly of end-tethered nanorods in a neat system and role of block fractions and aspect ratio , 2010 .

[35]  A. Travesset,et al.  Design of polymer nanocomposites in solution by polymer functionalization. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.