FYSM350, MD PART

General info

Lecture notes: L1
L2
On L3 we'll discuss some material used in the corresponding course (Kuronen and Nordlund) in HY: from the link see Lectures 5,8,9
L4
YOUR MD PROJECT IS DESCRIBED IN LECTURE 5: L5

A sample fortran MD code from openmp.org web site is here. We will discuss 28.3. a structure of a "good" MD program using this code as an example. One can take this code as platform for your project (if you are using fortran) but it will need several modifications and additions. (Do not care about the OpenMP directives - they would make the calculation run parallel if you would try this code in a machine where OpenMP works. This is off the topic of our course. Anyway this code works and can be run in serial CPU, the directives are then ignored as comments.)

Just for fun - visualisation (avi-file) of the dynamcis of 26 Coulomb particles in a harmonic trap. The temperature is increased during the run and the system "melts " in an interesting way. Matlab - simulation by Timo Santa-Nokki 2004. The visualisation works at least with QuickTime viewer.

Bits and pieces, sample f90 routines are discussed here.

Simple tests for a MD code are found here. We will discuss these tests during the lecture Tuesday 4.4.

Coordinate files for:

32-atom fcc lattice in internal LJ units, rho=1, corresponding to the lattice parameter a = sqrt(2)*r0, r0=2**(1/6)sigma (r0=position of the bottom of the LJ well) in LJ units (sigma = 1) for MD test # 1,
the same converted to xyz format for VMD (sigma=3.4 A for Argon).

256-atom fcc lattice with the same lattice parameter
the same converted to xyz format for VMD .

Back to FYSM350 homepage.