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.