exercises:2017_ethz_mmm:t_melting_2017
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
exercises:2017_ethz_mmm:t_melting_2017 [2017/03/24 09:50] – dpasserone | exercises:2017_ethz_mmm:t_melting_2017 [2020/08/21 10:15] (current) – external edit 127.0.0.1 | ||
---|---|---|---|
Line 39: | Line 39: | ||
<code bash> | <code bash> | ||
- | you@hypatia ~$ module load new cp2k | + | [you@hypatia ~]$ mmm-init |
+ | [you@hypatia ~]$ cd exercise_5.1 | ||
</ | </ | ||
+ | is then the only thing you need to do in order to initialize the m_* scripts. This time the exercise is already in your home directory, and the cp2k module is already loaded! | ||
- | you@eulerX ~$ module load courses mmm vmd | ||
- | |||
- | you@eulerX ~$ mmm-init | ||
</ | </ | ||
- | <note important> | + | <note important> |
<code bash> | <code bash> | ||
- | you@eulerX | + | [you@hypatia |
</ | </ | ||
Line 59: | Line 58: | ||
<code bash> | <code bash> | ||
- | you@eulerX ~$ bsub < jobname | + | [stb@hypatia exercise_5.1]$ qsub -v INP=file run |
+ | [stb@hypatia exercise_5.1]$ qstat -u stb # checks if it is running | ||
</ | </ | ||
+ | where **file** is the name of the input file without suffix. | ||
</ | </ | ||
In this exercise, we will use a slab geometry (without vacuum region, so without a surface) with full periodic boundary conditions to study the melting behavior of copper. | In this exercise, we will use a slab geometry (without vacuum region, so without a surface) with full periodic boundary conditions to study the melting behavior of copper. | ||
+ | All files of this exercise (**all inputs and scripts are commented**) can be also downloaded from the wiki: {{exercise5.1.zip|exercise5.1.zip}} | ||
- | * Download | + | * First of all, we will test the NOSE-HOOVER thermostat |
- | <code bash> | + | |
- | you@eulerX ~$ wget http:// | + | |
- | you@eulerX ~$ unzip exercises: | + | |
- | you@eulerX ~$ cd exercise_5.1 | + | |
- | </ | + | |
- | + | ||
<note tip> | <note tip> | ||
- | All files of this exercise (**all inputs | + | - Take a look at the file 111.xyz with vmd. |
+ | - We will apply a thermostat to half of the cell. | ||
+ | - copy the half_TIMECON.inp into half_20.inp | ||
+ | - Edit the half_20.inp | ||
+ | - Run the job (interactively): > cp2k.popt -i half_20.inp -o half_20.out | ||
+ | - Plot the temperature behavior with gnuplot, file half_20-1.ener | ||
+ | - Take a look at the half_20.inp file. How is the temperature controlled? Are all particles moving? Why? Which are the relevant sections for MD? Which kind of MD is it? | ||
+ | - Repeat by copying half_TIMECON.inp into half_300.inp (TIMECON 300) | ||
+ | - Plot the -growing- half*ener file with gnuplot together with the previous ones. Comment the differences | ||
</ | </ | ||
- | * Now, run the first simulation, that should melt your system: | + | * Now, run the simulation that should melt your system: |
<code bash> | <code bash> | ||
- | you@eulerX | + | you@hypatia |
</ | </ | ||
Line 89: | Line 92: | ||
<note tip> | <note tip> | ||
- Take a look at the file 111.xyz with vmd. Visualize it on the screen, and try to reproduce the figure similar to the one on the last slide of the lectures of today. Include the pbc box, create a representation with vdw, periodic images, rotate the sample, etc. Produce a snapshot and include the file in your assignment. | - Take a look at the file 111.xyz with vmd. Visualize it on the screen, and try to reproduce the figure similar to the one on the last slide of the lectures of today. Include the pbc box, create a representation with vdw, periodic images, rotate the sample, etc. Produce a snapshot and include the file in your assignment. | ||
- | - Take a look at the half.inp file. How is the temperature controlled? Are all particles moving? Why? Which are the relevant sections for MD? Which kind of MD is it? | + | |
- | | + | |
</ | </ | ||
Line 96: | Line 98: | ||
<code bash> | <code bash> | ||
- | you@eulerX | + | you@hypatia |
</ | </ | ||
Line 104: | Line 106: | ||
<code bash> | <code bash> | ||
- | you@eulerX | + | you@hypatia |
</ | </ | ||
or download the file to your local machine and open in your internet browser. It will run the animation. | or download the file to your local machine and open in your internet browser. It will run the animation. | ||
Line 116: | Line 118: | ||
<code bash> | <code bash> | ||
- | you@eulerX | + | you@hypatia |
</ | </ | ||
Line 135: | Line 137: | ||
* Copy the files TEMPnpe.init.inp and TEMPnpe.inp into 1300npe.init.inp and 1300npe.inp (for T=1300) and then edit them in the appropriate points: PROJECT name, INITIAL temperature and RESTART filename. | * Copy the files TEMPnpe.init.inp and TEMPnpe.inp into 1300npe.init.inp and 1300npe.inp (for T=1300) and then edit them in the appropriate points: PROJECT name, INITIAL temperature and RESTART filename. | ||
- | * Run the first simulation: | + | * Run the first simulation: |
- | * Run the second simulation: | + | * Run the second simulation: |
* Observe the temperature and the z profile. Can you find the melting temperature? | * Observe the temperature and the z profile. Can you find the melting temperature? | ||
exercises/2017_ethz_mmm/t_melting_2017.1490349021.txt.gz · Last modified: 2020/08/21 10:15 (external edit)