Molecular dynamics with the massively parallel APE computers

In this paper we discuss the general question of the portability of Molecular Dynamics codes for diffusive systems (liquids) on parallel computers of the APE family. The intrinsic single precision arithmetics of the today available APE platforms does not seem to affect the numerical accuracy of the simulations, while the absence of integer addressing from the CPU to individual nodes puts strong constraints on the possible programming strategies. As a test case, we report the results of the simulation of the dynamics of 512 molecules of liquid butane (C4H10) ta room temperature. After 30 to 50 ps of equilibrium, the system was followed along four long trajectories, each one more than 1.3 ns. The effective CPU time corresponding to the simulation of a trajectory of 1 ns on the Torre (512 nodes ≈ 25 Gigaflops) was ≈ 50 hours. The CPU time can be substantially reduced (by almost a factor 3), if the APE-assembler micro-code of the most time-consuming part of the program is carefully optimized.

[1]  B. Berne,et al.  A Multiple-Time-Step Molecular Dynamics Algorithm for Macromolecules , 1994 .

[2]  R. C. Weast.  Handbook of chemistry and physics , 1973 .

[3]  Anomalous fluctuations in low-temperature molecular dynamics simulations. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[4]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[5]  Masakatsu Watanabe,et al.  Dynamics of molecules with internal degrees of freedom by multiple time-step methods , 1993 .

[6]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[7]  M. Karplus,et al.  SIMULATIONS OF MACROMOLECULES BY MULTIPLE TIME-STEP METHODS , 1995 .

[8]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

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

[10]  William G. Hoover,et al.  Nonequilibrium molecular dynamics via Gauss's principle of least constraint , 1983 .

[11]  D. Osguthorpe,et al.  Structure and energetics of ligand binding to proteins: Escherichia coli dihydrofolate reductase‐trimethoprim, a drug‐receptor system , 1988, Proteins.

[12]  Piero Procacci,et al.  Computer simulation of solid C60 using multiple time‐step algorithms , 1994 .

[13]  William L. Jorgensen,et al.  Optimized intermolecular potential functions for liquid hydrocarbons , 1984 .

[14]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[15]  G. Benettin,et al.  Kolmogorov Entropy and Numerical Experiments , 1976 .

[16]  Pettini,et al.  Strong stochasticity threshold in nonlinear large Hamiltonian systems: Effect on mixing times. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[17]  Nicolai Petkov,et al.  Systolic Parallel Processing , 1992 .

[18]  Luciano Barone,et al.  Molecular dynamics on APE100 , 1995 .

[19]  W. C. Swope,et al.  A computer simulation method for the calculation of equilibrium constants for the formation of physi , 1981 .

[20]  M. Pettini,et al.  Relaxation properties and ergodicity breaking in nonlinear Hamiltonian dynamics. , 1990, Physical review. A, Atomic, molecular, and optical physics.