md

md, a Python code which carries out a molecular dynamics simulation.

The computation involves following the paths of particles which exert a distance-dependent force on each other. The particles are not constrained by any walls; if particles meet, they simply pass through each other.

To compute the next position of each particle requires the evaluation of the right hand side of its corresponding differential equation. Since each of these calculations is independent, there is a potential speedup if the program can take advantage of parallel computing.

This version of the program is intended as a starting point; students are asked to examine the program and add the appropriate modifications to run it in parallel.

Usage:

md ( nd, np, step_num, dt )

• nd is the spatial dimension, usually 2 or 3.
• np is the number of particles in the system. A typical value might be 1000.
• step_num is the number of time steps, which might be 500 or 1,000.
• dt is the time step size, which might be 0.1.

Licensing:

The information on this web page is distributed under the MIT license.

Languages:

md is available in a C version and a C++ version and a Fortran90 version and a MATLAB version and an Octave version and a Python version.

Related Data and Programs:

hello_mpi, a Python code which prints out "Hello, world!" using the MPI parallel programming environment, under MPI4PY.

prime_mpi, a Python code which counts the number of primes between 1 and N, using MPI for parallel execution.

quad_mpi, a Python code which approximates an integral using a quadrature rule, and carries out the computation in parallel using MPI.

search_mpi, a Python code which searches integers between A and B for a value J such that F(J) = C, using MPI and MPI4PY for parallel execution.

Reference:

1. Matthew Allen, Dominic Tildesley,
   Computer Simulation of Liquids,
   Oxford University Press, 1987,
   ISBN: 0198556454,
   LC: QC145.2.
2. Peter Arbenz, Wesley Petersen,
   Introduction to Parallel Computing - A practical guide with examples in C,
   Oxford University Press,
   ISBN: 0-19-851576-6,
   LC: QA76.58.P47.
3. Rohit Chandra, Leonardo Dagum, Dave Kohr, Dror Maydan, Jeff McDonald, Ramesh Menon,
   Parallel Programming in OpenMP,
   Morgan Kaufmann, 2001,
   ISBN: 1-55860-671-8,
   LC: QA76.642.P32.
4. Barbara Chapman, Gabriele Jost, Ruud vanderPas, David Kuck,
   Using OpenMP: Portable Shared Memory Parallel Processing,
   MIT Press, 2007,
   ISBN13: 978-0262533027,
   LC: QA76.642.C49.
5. Furio Ercolessi,
   A Molecular Dynamics Primer.
6. Philipp Janert,
   Gnuplot in Action: Understanding Data with Graphs,
   Manning, 2008,
   ISBN: 1-933988-39-8.
7. Dennis Rapaport,
   An Introduction to Interactive Molecular-Dynamics Simulation,
   Computers in Physics,
   Volume 11, Number 4, July/August 1997, pages 337-347.
8. Dennis Rapaport,
   The Art of Molecular Dynamics Simulation,
   Cambridge University Press, 2004,
   ISBN: 0521825687.

Source Code:

• md.py, the source code.
• md.sh, runs all the tests.
• md.txt, the output file;