CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace:
exercises:2018_ethz_mmm:h2o_md

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
Last revisionBoth sides next revision
exercises:2018_ethz_mmm:h2o_md [2018/03/02 11:13] – dpasseroneexercises:2018_ethz_mmm:h2o_md [2018/03/02 15:36] – dpasserone
Line 9: Line 9:
 Download the 2.1 exercise into your **EXERCISES** folder and unzip it.  Download the 2.1 exercise into your **EXERCISES** folder and unzip it. 
  
-<code> +
-max@qmobile:~$ cd ; cd EXERCISES +
-max@qmobile:~$ wget https://polybox.ethz.ch/index.php/s/ULaq6rjhg4Pu63A +
-max@qmobile:~$ unzip exercise_2.1.zip +
-max@qmobile:~$ cd exercise_2.1 +
-</code>+
  
 This exercise is mostly taken by a [[https://www.cp2k.org/exercises:2015_uzh_molsim:index|previous lab session]] by Marcella Iannuzzi, UZH, who should be credited and acknowledged here. This exercise is mostly taken by a [[https://www.cp2k.org/exercises:2015_uzh_molsim:index|previous lab session]] by Marcella Iannuzzi, UZH, who should be credited and acknowledged here.
Line 247: Line 242:
 <code>./vel_autocorr.py T300-vel-1.xyz 3 10 100  <code>./vel_autocorr.py T300-vel-1.xyz 3 10 100 
 ./vel_autocorr_plot.py vel_autocorr.out ./vel_autocorr_plot.py vel_autocorr.out
-./mean_square_displ.py T300.xyz 3 10 100+./mean_square_disp.py T300.xyz 3 10 100
 ./msd_fit_and_plot.py mean_square_disp.out ./msd_fit_and_plot.py mean_square_disp.out
 </code> </code>
 +Remember: to display the **png** files, the command is <code>display file.png</code>
 <note important>Assignments <note important>Assignments
-  - look at the file T300-1.ener at the end of the simulation. It contains several quantities. Check the header+  - look at the file T300-1.ener at the end of the simulation. It contains several quantities. Check the header!
   - Using gnuplot, check the stability of the MD, you can also use <code> ./simpleplot.py Col1 Col2 </code>   - Using gnuplot, check the stability of the MD, you can also use <code> ./simpleplot.py Col1 Col2 </code>
 +  - How do you interpret the oscillation of the potential energy? And of temperature?
   - Compare the values of the diffusion coefficient with the two methods. Are they compatible?    - Compare the values of the diffusion coefficient with the two methods. Are they compatible? 
   - Look at the trajectories with vmd. Why are the particles "exiting" the box with time? What is the importance of this for the algorithm?   - Look at the trajectories with vmd. Why are the particles "exiting" the box with time? What is the importance of this for the algorithm?
Line 258: Line 255:
   - Copy md300.in into md400.in and change to 400 in the two places. Run the md. What is the final temperature? Why?   - Copy md300.in into md400.in and change to 400 in the two places. Run the md. What is the final temperature? Why?
   - Check the temperature dependence of the diffusion coefficient. Plot the result.   - Check the temperature dependence of the diffusion coefficient. Plot the result.
-</note> 
 </note> </note>
  
  
exercises/2018_ethz_mmm/h2o_md.txt · Last modified: 2020/08/21 10:15 by 127.0.0.1