User Tools

Site Tools


exercises:2015_ethz_mmm:surface_cu

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
Next revisionBoth sides next revision
exercises:2015_ethz_mmm:surface_cu [2015/03/06 09:20] yakutovichexercises:2015_ethz_mmm:surface_cu [2015/03/06 10:58] yakutovich
Line 1: Line 1:
 +
 +====== Surface energies of Copper high-symmetry surfaces ======
 +
 <note warning> <note warning>
 TO USE THE FUNCTION LIBRARY (VERSION UP TO DATE) IN THE INTERACTIVE SHELL: TO USE THE FUNCTION LIBRARY (VERSION UP TO DATE) IN THE INTERACTIVE SHELL:
Line 7: Line 10:
 <note important> **REMEMBER: this is the command to load the  module for the cp2k program:** <note important> **REMEMBER: this is the command to load the  module for the cp2k program:**
  
-<code>+<code bash>
 you@eulerX ~$ module load new cp2k you@eulerX ~$ module load new cp2k
 </code> </code>
Line 13: Line 16:
 **and to submit the job:** **and to submit the job:**
  
-<code>+<code bash>
 you@eulerX ~$ bsub < jobname you@eulerX ~$ bsub < jobname
 </code>  </code> 
Line 19: Line 22:
  
  
-====== Generation of slab with high symmetry surfaces ====== 
 {{:exercises:2015_ethz_mmm:lect2_surf.001.jpg?nolink&400|}} {{:exercises:2015_ethz_mmm:lect2_surf.001.jpg?nolink&400|}}
 +
 +----
 +
 {{:exercises:2015_ethz_mmm:lect2_surf.002.jpg?400|}} {{:exercises:2015_ethz_mmm:lect2_surf.002.jpg?400|}}
 +
 +----
 +
 {{:exercises:2015_ethz_mmm:lect2_surf.003.jpg?400|}} {{:exercises:2015_ethz_mmm:lect2_surf.003.jpg?400|}}
  
  
-====== Surface energies of Copper high-symmetry surfaces ====== 
 In this exercise we will compute the surface energies of Cu using the EAM potential. In this exercise we will compute the surface energies of Cu using the EAM potential.
 As a reference, we report the table from the Gross book: As a reference, we report the table from the Gross book:
Line 31: Line 38:
 {{screen_shot_2014-03-07_at_6.55.49_am.png?direct&500 |http://example.com|External Link}} {{screen_shot_2014-03-07_at_6.55.49_am.png?direct&500 |http://example.com|External Link}}
  
- +---- 
-  * Download the 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: {{exercise_3.1.zip|exercise_3.1.zip}} (**all inputs are commented**) in your home directory and unzip it: 
-<code>+<code bash>
 you@eulerX ~$ wget http://www.cp2k.org/_media/exercises:2015_ethz_mmm:exercise_3.1.zip you@eulerX ~$ wget http://www.cp2k.org/_media/exercises:2015_ethz_mmm:exercise_3.1.zip
 you@eulerX ~$ unzip exercises:2015_ethz_mmm:exercise_3.1.zip you@eulerX ~$ unzip exercises:2015_ethz_mmm:exercise_3.1.zip
 you@eulerX ~$ cd exercise_3.1 you@eulerX ~$ cd exercise_3.1
 </code> </code>
-  * Run the optimizations 100.inp, 110.inp, 111.inp. +  * Run the optimizations 100.inp, 110.inp, 111.inp and the bulk
-<code>+<code bash>
 you@eulerX exercise_3.1$ bsub cp2k.popt -i 100.inp -o 100.out you@eulerX exercise_3.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_3.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_3.1$ bsub cp2k.popt -i 111.inp -o 111.out
 +you@eulerX exercise_3.1$ bsub cp2k.popt -i bulk.inp -o bulk.out
 </code> </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.   * 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>+<code bash>
 you@eulerX exercise_3.1$ vmd 100.xyz you@eulerX exercise_3.1$ vmd 100.xyz
 you@eulerX exercise_3.1$ vmd 110.xyz you@eulerX exercise_3.1$ vmd 110.xyz
Line 51: Line 59:
 </code> </code>
  
-  Compute the three surface energies: you need to compute the area, subtract bulk contribution, take care of the units. +  Compute the three surface energies: you need to compute the area, subtract bulk contribution, take care of the units. 
-  - In the class we will learn how to compute the Wulff crystal from these three numbers. +  * Unordered List ItemIn the class we will learn how to compute the Wulff crystal from these three numbers. 
-  Basically, you should edit the input file "DP", replacing the column beginning with **"PUT_HERE_THE_XXX_SURFACE_ENERGY"** with quantities that are proportional to the three surface free energies (see Wulff theorem). Note that for the {100} and {110} surfaces you need to edit **TWO** lines each, since the **z** direction is treated differently in this program.  +  Basically, you should edit the input file "DP", replacing  **"PUT_HERE_THE_XXX_SURFACE_ENERGY"** with corresponding surface energy. Note that for the {100} and {110} surfaces you need to edit **TWO** lines eachthe program takes care of assigning the surface energy to equivalent surfaces (like (010)) but since it was developed for island on surfaces it does not do it automatically in all directions. **z** direction is treated differently in this program.  
-  At this point, you run the program with the command **./sowos.v02.00.02.x**  +  At this point, you can run the sowos program
-  There will be many output files. Important are  +<code bash> 
-  -   * 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. If you appropriately modify it (how?) you can have a xyz file to be opened in vmd.   +you@eulerX exercise_3.1$ ./sowos.v02.00.02.x 
-  -   * the **out.plot-gnuplot.plt** file. If you launch gnuplot, then "**load "out.plot-gnuplot.plt"**  a plot of the crystal will be shown. To have it more or less to scale in all directions: **set xrange [-80:80];set yrange [-80:80]set view equal xyzreplot** +</code> 
-  -   * you can rotate with the mouse!+  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. If you appropriately modify it (how?) you can have a xyz file to be opened in vmd.   
 +  - the **out.plot-gnuplot.plt** file. You can open it with gnuplot: 
 +<code bash> 
 +you@eulerX exercise_3.1$ gnuplot 
 +gnuplot> load "out.plot-gnuplot.plt" 
 +gnuplot> set xrange [-80:80] 
 +gnuplot> set yrange [-80:80] 
 +gnuplot> set view equal xyz 
 +gnuplot> replot 
 +</code> 
 +  * you can rotate with the mouse!
  
 <note important>**Assignment: use the cluster generated with SOWOS as an input configuration for a cluster optimization with cp2k. Comment on the final geometry. BEWARE OF THE CELL! It is not a periodic system!**</note> <note important>**Assignment: use the cluster generated with SOWOS as an input configuration for a cluster optimization with cp2k. Comment on the final geometry. BEWARE OF THE CELL! It is not a periodic system!**</note>
exercises/2015_ethz_mmm/surface_cu.txt · Last modified: 2020/08/21 10:15 by 127.0.0.1