CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace:
exercises:2017_ethz_mmm:lennard_jones_cluster

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:2017_ethz_mmm:lennard_jones_cluster [2017/02/24 09:07] – dpasseroneexercises:2017_ethz_mmm:lennard_jones_cluster [2018/04/14 14:49] – replace pdf with doi link oschuett
Line 1: Line 1:
 ====== 38 atom Lennard-Jones cluster ====== ====== 38 atom Lennard-Jones cluster ======
 +
 +{{:exercises:2017_ethz_mmm:lj38bs.jpg?400|}} (picture by Luke Abraham)
 <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 31: Line 33:
  
 <note tip> <note tip>
-All files of this exercise (**input and scripts are all commented**) can be downloaded from the wiki: {{exercise_1.1.zip|}} +All files of this exercise be downloaded from the wiki: {{exercise_1.1.zip|}} 
 </note> </note>
  
 In this exercise you will test the Lennard-Jones potential. In particular, we will focus on the system described in the following paper about the energy landscape of the 38 atom Lennard-Jones cluster: In this exercise you will test the Lennard-Jones potential. In particular, we will focus on the system described in the following paper about the energy landscape of the 38 atom Lennard-Jones cluster:
-<note tip>{{ :exercises:2017_ethz_mmm:1999_the_double-funnel_energy_landscape_of_the_38-atom_lennard-jones_cluster.pdf |}}+<note tip>[[doi>10.1063/1.478595]]
 </note> </note>
 Login to euler using your nethz credentials. Login to euler using your nethz credentials.
Line 166: Line 168:
   - extract the q4 order parameter from **fcc.xyz** and from **fcc_rand.xyz** and compare the values.<code>module load new gcc/4.8.2 python/2.7.12     - extract the q4 order parameter from **fcc.xyz** and from **fcc_rand.xyz** and compare the values.<code>module load new gcc/4.8.2 python/2.7.12  
 python stein.py file.xyz </code>. You will be asked the cutoff radius for the neighbors, it is **1.391** in sigma units. **You should input it in Angstrom**.  python stein.py file.xyz </code>. You will be asked the cutoff radius for the neighbors, it is **1.391** in sigma units. **You should input it in Angstrom**. 
-  - before running the simulation, copy the input coordinate file into in.xyz <code>cp fcc.xyz in.xyz</code> +  - before running the simulation, copy the input coordinate file into in.xyz <code>cp fcc_rand.xyz in.xyz</code> 
-  - run cp2k <code>module load cp2k</code>(this has only to be done once)<code>cp2k.opt -i geo_opt.inp -o geo_opt.out </code> +  - run cp2k <code>module load cp2k</code>(this has only to be done once)<code>cp2k.popt -i geo_opt.inp -o geo_opt.out </code> 
   - in the output file, note the final energy, **transform it in the unit of the paper (epsilon units)**   - in the output file, note the final energy, **transform it in the unit of the paper (epsilon units)**
   - load vmd module and play with the optimization trajectory <code>vmd OPT-pos-1.xyz</code> (ask the teacher)   - load vmd module and play with the optimization trajectory <code>vmd OPT-pos-1.xyz</code> (ask the teacher)
exercises/2017_ethz_mmm/lennard_jones_cluster.txt · Last modified: 2020/08/21 10:15 by 127.0.0.1