exercise_7

# exercise_7 - Molecular Dynamics Simulations 2011 Exercise 7...

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

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

View Full Document
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. It’stime....
View Full Document

{[ snackBarMessage ]}

### Page1 / 7

exercise_7 - Molecular Dynamics Simulations 2011 Exercise 7...

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

View Full Document
Ask a homework question - tutors are online