Hardware accelerator for molecular dynamics: MDGRAPE-2

Abstract We developed MDGRAPE-2, a hardware accelerator that calculates forces at high speed in molecular dynamics (MD) simulations. MDGRAPE-2 is connected to a PC or a workstation as an extension board. The sustained performance of one MDGRAPE-2 board is 15 Gflops, roughly equivalent to the peak performance of the fastest supercomputer processing element. One board is able to calculate all forces between 10 000 particles in 0.28 s (i.e. 310000 time steps per day). If 16 boards are connected to one computer and operated in parallel, this calculation speed becomes ∼10 times faster. In addition to MD, MDGRAPE-2 can be applied to gravitational N -body simulations, the vortex method and smoothed particle hydrodynamics in computational fluid dynamics.

[1]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[2]  Junichiro Makino.  Yet Another Fast Multipole Method without Multipoles-Pseudoparticle Multipole Method , 1999 .

[3]  B. Quentrec,et al.  New method for searching for neighbors in molecular dynamics computations , 1973 .

[4]  T. Narumi,et al.  46 TFLOPS special-purpose computer for molecular dynamics simulations: WINE-2 , 2000, WCC 2000 - ICSP 2000. 2000 5th International Conference on Signal Processing Proceedings. 16th World Computer Congress 2000.

[5]  Weitao Yang,et al.  A simple O(N log N) algorithm for the rapid evaluation of particle-particle interactions , 1995 .

[6]  Toshikazu Ebisuzaki,et al.  A Highly Parallelized Special-Purpose Computer for Many-Body Simulations with an Arbitrary Central Force: MD-GRAPE , 1996 .

[7]  Toshikazu Ebisuzaki,et al.  A special-purpose computer for gravitational many-body problems , 1990, Nature.

[8]  Toshiyuki Fukushige.  WINE-1--Special-Purpose Computer for N-Body Simulations with a Periodic Boundary Condition (GRAPE) , 1993 .

[9]  Toshikazu Ebisuzaki,et al.  A highly vectorised "link-cell" FORTRAN code for the DL_POLY molecular dynamics simulation package , 2000 .

[10]  Minoru Saito,et al.  Molecular dynamics simulations of proteins in water without the truncation of long-range Coulomb interactions , 1992 .

[11]  Toshikazu Ebisuzaki,et al.  A Special-Purpose Computer for N-Body Simulations:GRAPE-2A (GRAPE) , 1993 .

[12]  B. Berne,et al.  A new molecular dynamics method combining the reference system propagator algorithm with a fast multipole method for simulating proteins and other complex systems , 1995 .

[13]  Masayuki Hata,et al.  Hydrolysis Mechanism of the Phenylalanine−Proline Peptide Bond Specific to HIV-1 Protease: Investigation by the ab Initio Molecular Orbital Method , 1999 .

[14]  L. Lucy.  A numerical approach to the testing of the fission hypothesis. , 1977 .

[15]  Makoto Taiji,et al.  Fast and accurate molecular dynamics simulation of a protein using a special‐purpose computer , 1997 .

[16]  T Ito,et al.  A special‐purpose computer for molecular dynamics: GRAPE‐2A , 1994, Proteins.

[17]  Toshikazu Ebisuzaki,et al.  Smoothed particle hydrodynamics with GRAPE-1A , 1993 .

[18]  Atsushi Kawai,et al.  Pseudoparticle Multipole Method: A Simple Method to Implement a High-Accuracy Tree Code , 2000, astro-ph/0012041.

[19]  Makoto Haraguchi,et al.  Parallel molecular dynamics simulation of a protein , 2001, Parallel Comput..

[20]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[21]  Makino,et al.  Astrophysics on the GRAPE family of special-purpose computers , 1999, Science.

[22]  T Fukushige,et al.  A high performance system for molecular dynamics simulation of biomolecules using a special-purpose computer. , 1996, Pacific Symposium on Biocomputing. Pacific Symposium on Biocomputing.

[23]  R. W. Hockney,et al.  A 10000 particle molecular dynamics model with long range forces , 1973 .

[24]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[25]  Takashi Amisaki,et al.  Development of MD Engine: High‐speed accelerator with parallel processor design for molecular dynamics simulations , 1999 .

[26]  Masuhiro Mikami,et al.  Computationally efficient canonical molecular dynamics simulations by using a multiple time‐step integrator algorithm combined with the particle mesh Ewald method and with the fast multipole method , 2000 .

[27]  J. Monaghan,et al.  Smoothed particle hydrodynamics: Theory and application to non-spherical stars , 1977 .

[28]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .