exercise_7 - Molecular Dynamics Simulations 2011 Exercise...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Molecular Dynamics Simulations 2011 Exercise 7: Suggested solutions Eero Holmstr om November 2, 2011 1. The pressure control was implemented in main.c/f90 , and the virial calcula- tion is done in forces.c/f90 . So, first we run the simulation with the usual atoms.in to get the equilibrium box size, which is found to be 18.833 A. Then we take the last frame of atoms.out , and scale all the coordinates by 1.02. This can be done with the following command (on a single line): tail -502 atoms.out | awk BEGIN{s=1.02} {if($1=="Cu"){print $1, $2*s, $3*s, $4*s, $5} else {print $0} } > atoms.in.2 Also, we must set the correct box size in mdmorse.in , which is 1 . 02 18 . 833. Then we run the simulation using the file atoms.in.2 as the input atom coordinates. And similarly for the case where we scale all coordinates by 0.98. Results are plotted in Figs 1 and 2. In Fig. 1 we have the box size in the x-direction (the same as in y and z due to symmetry) as a function of time. We see that for both cases the same equilibrium box size is found. Itstime....
View Full Document

Page1 / 7

exercise_7 - Molecular Dynamics Simulations 2011 Exercise...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online