Differences

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

--- exercises:2017_ethz_mmm:surface_cu [2017/03/17 04:20] – dpasserone
+++ exercises:2017_ethz_mmm:surface_cu [2017/03/17 14:59] – dpasserone
@@ Line 2: / Line 2: @@
 ====== Surface energies of Copper high-symmetry surfaces ======
+<note warning>
+TO USE THE FUNCTION LIBRARY (VERSION UP TO DATE) IN THE INTERACTIVE SHELL:
+you@eulerX ~$ module load courses mmm vmd
+you@eulerX ~$ mmm-init
+</note>
 <note important> **REMEMBER: this is the command to load the  module for the cp2k program:**
@@ Line 25: / Line 32: @@
 ----
-  * Download all the necessary files from  from the wiki: {{exercise_3.1.zip|exercise_3.1.zip}} (**all inputs are commented**) in your home directory and unzip it:
+  * Download all the necessary files from  from the wiki: {{e4.1.zip|e4.1.zip}} (**all inputs are commented**) in your home directory and unzip it:
 <code bash>
-you@eulerX ~$ wget http://www.cp2k.org/_media/exercises:2017_ethz_mmm:exercise_3.1.zip
+you@eulerX ~$ wget http://www.cp2k.org/_media/exercises:2017_ethz_mmm:e4.1.zip
-you@eulerX ~$ unzip exercises:2016_ethz_mmm:exercise_3.1.zip
+you@eulerX ~$ unzip exercises:2017_ethz_mmm:e4.1.zip
-you@eulerX ~$ cd exercise_3.1
+you@eulerX ~$ cd exercise_4.1
 </code>
+If you cannot download with **wget**, you can copy the zip file directly from my directory
+<code bash>
+you@eulerX ~$ cp /cluster/home/danielep/2017/e4.1.zip .
+you@eulerX ~$ unzip e4.1.zip
+you@eulerX ~$ cd exercise_4.1
+</code>
   * Run the optimizations 100.inp, 110.inp, 111.inp and the bulk.
 <code bash>
-you@eulerX exercise_3.1$ bsub cp2k.popt -i 100.inp -o 100.out
+you@eulerX exercise_4.1$ bsub cp2k.popt -i 100.inp -o 100.out
-you@eulerX exercise_3.1$ bsub cp2k.popt -i 110.inp -o 110.out
+you@eulerX exercise_4.1$ bsub cp2k.popt -i 110.inp -o 110.out
-you@eulerX exercise_3.1$ bsub cp2k.popt -i 111.inp -o 111.out
+you@eulerX exercise_4.1$ bsub cp2k.popt -i 111.inp -o 111.out
-you@eulerX exercise_3.1$ bsub cp2k.popt -i bulk.inp -o bulk.out
+you@eulerX exercise_4.1$ bsub cp2k.popt -i bulk.inp -o bulk.out
 </code>
   * While geometry optimization is running you can have a look at the corresponding initial coordinate files 100.xyz, 110.xyz, 111.xyz. In vmd it is also possible to open a console, and give the command **pbc set { a b c 90 90 90 }**  where a, b, c can be extracted from the cp2k input file. Then you can make several periodic copies for visualization.
 <code bash>
-you@eulerX exercise_3.1$ vmd 100.xyz
+you@eulerX exercise_4.1$ vmd 100.xyz
-you@eulerX exercise_3.1$ vmd 110.xyz
+you@eulerX exercise_4.1$ vmd 110.xyz
-you@eulerX exercise_3.1$ vmd 111.xyz
+you@eulerX exercise_4.1$ vmd 111.xyz
 </code>
@@ Line 55: / Line 70: @@
   * At this point, you can run the sowos program:
 <code bash>
-you@eulerX exercise_3.1$ ./sowos.v02.00.02.x
+you@eulerX exercise_4.1$ ./sowos.v02.00.02.x
 </code>
   * There will be many output files. Important are:
   - the atomistic model **out.atomistic-inside-gnuplot.xyz** which will contain many atoms if you chose a proportionality constant in the file DP which is too large. Unfortunately this file does NOT have a proper structure of an xyz file.
+  - This can be solved by using a simple m_* function that adds the atomic symbol at the beginning of each line, and adding two lines at the beginning of the resulting file (number of atoms; comment):
+<code bash>
+you@eulerX exercise_4.1$ m_addcolumn Au < out.atomistic-inside-gnuplot.xyz > cluster.xyz
+you@eulerX exercise_4.1$ nano cluster.xyz # Add two lines at the beginning.
+</code>
   - the **out.plot-gnuplot.plt** file. You can open it with gnuplot
 <code bash>
-you@eulerX exercise_3.1$ gnuplot
+you@eulerX exercise_4.1$ gnuplot
 gnuplot> load "out.plot-gnuplot.plt"
 gnuplot> set xrange [-80:80]