exercises:2016_ethz_mmm:t_melting
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision | ||
exercises:2016_ethz_mmm:t_melting [2016/04/08 00:02] – yakutovich | exercises:2016_ethz_mmm:t_melting [2016/04/08 10:27] – dpasserone | ||
---|---|---|---|
Line 29: | Line 29: | ||
* Download the 5.1 exercise into your $HOME folder and unzip it: | * Download the 5.1 exercise into your $HOME folder and unzip it: | ||
<code bash> | <code bash> | ||
- | you@eulerX ~$ wget http:// | + | you@eulerX ~$ wget http:// |
- | you@eulerX ~$ unzip exercises: | + | you@eulerX ~$ tar -zxvf exercises: |
you@eulerX ~$ cd exercise_5.1 | you@eulerX ~$ cd exercise_5.1 | ||
</ | </ | ||
Line 62: | Line 62: | ||
you@eulerX exercise_5.1$ cp2k.popt -i cell.inp -o cell.out | you@eulerX exercise_5.1$ cp2k.popt -i cell.inp -o cell.out | ||
</ | </ | ||
+ | |||
+ | <note tip> | ||
+ | * Repeat the simulation with a larger cutoff. How does the cohesive energy per atom change? Why? How does the equilibrium lattice parameter change? Why? | ||
+ | </ | ||
<note important> | <note important> | ||
- | Pay attention at the folowing | + | Pay attention at the following |
MULTIPLE_UNIT_CELL 4 2 2 | MULTIPLE_UNIT_CELL 4 2 2 | ||
Line 88: | Line 92: | ||
- | It is a 5000 step molecular dynamics. While it is running you can complete the folowing assignments: | + | It is a 10000 step molecular dynamics. While it is running you can complete the folowing assignments: |
<note tip> | <note tip> | ||
- Take a look at the file ELEMENT_opt_unit.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 ELEMENT_opt_unit.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. | ||
Line 115: | Line 119: | ||
- | Now, starting from the restart of this simulation, we equilibrate the system in nve, and we move all particles. But again some parameters needs to be specified in the file nve.inp. Use unit cell from the file half-1.restart. Also file with coordinates (X_init_nve.xyz) needs to be created. Use the coordinates from the last frame of the previous run. | + | Now, starting from the restart of this simulation, we equilibrate the system in nve, and we move all particles. But again some parameters needs to be specified in the file nve.inp. Use the simulation |
And run the simulations: | And run the simulations: | ||
Line 126: | Line 130: | ||
The resulting configuration (check) will be an equilibrated system (which profile?). | The resulting configuration (check) will be an equilibrated system (which profile?). | ||
- | Now we have a file called "1400nve-1.restart" | + | Now we have a file called "nve-1.restart" |
Line 132: | Line 136: | ||
===== SIMULATIONS AT DIFFERENT TOTAL ENERGIES FOR DETERMINING THE MELTING TEMPERATURE ===== | ===== SIMULATIONS AT DIFFERENT TOTAL ENERGIES FOR DETERMINING THE MELTING TEMPERATURE ===== | ||
- | As explained in the class, we will run NPE (that is, constant energies but variable cell) simulations at energies which are above and below the supposed " | + | As explained in the class, we should |
+ | in this case. | ||
**THE TEMPERATURE WILL NOT BE CONTROLLED DURING THE RUN** | **THE TEMPERATURE WILL NOT BE CONTROLLED DURING THE RUN** | ||
Line 139: | Line 144: | ||
- | * 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 TEMPnve.init.inp and TEMPnve.inp into 100nve.init.inp and 100nve.inp (for T=100) and then edit them in the appropriate points: PROJECT name, INITIAL temperature and RESTART filename |
- | * Run the first simulation: bsub cp2k.popt -i 1300npe.init.inp | + | * Run the first simulation: bsub cp2k.popt -i 100nve.init.inp |
- | * Run the second simulation: bsub cp2k.popt -i 1300npe.inp > 1300npe.out | + | * Run the second simulation: bsub cp2k.popt -i 100nve.inp > 100nve.out |
* 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? | ||
Line 147: | Line 152: | ||
<note tip> | <note tip> | ||
- | * What is the melting temperature of copper | + | * What is the melting temperature of the noble gas that you have chosen? |
</ | </ |
exercises/2016_ethz_mmm/t_melting.txt · Last modified: 2020/08/21 10:15 by 127.0.0.1