CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace:
exercises:2015_ethz_mmm:c2h2_bond_energy

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:2015_ethz_mmm:c2h2_bond_energy [2015/02/27 07:39] – dpasseroneexercises:2015_ethz_mmm:c2h2_bond_energy [2015/06/02 11:31] – dpasserone
Line 1: Line 1:
 ====== C2H2 and C2H4 bond energy ====== ====== C2H2 and C2H4 bond energy ======
 <note warning> <note warning>
-Before you start it is **ABSOLUTELY NEСESSARY** to update the function library. So from you $HOME directory please run: +TO USE THE FUNCTION LIBRARY (VERSION UP TO DATE) IN THE INTERACTIVE SHELL: 
-<code> + 
-you@eulerX ~$ ./update_functionlibrary +you@eulerX ~$ module load courses mmm vmd 
-</code> +  
-** Pay attention to the output of this program. It may give you some helpful information**+you@eulerX ~$ mmm-init 
 +</note>
  
-If you don't have the file "update_functionlibrary" please download it from here {{update_functionlibrary.zip|}} and unzip into your $HOME directory 
-<code> 
-you@eulerX ~$ unzip update_functionlibrary.zip 
-you@eulerX ~$ chmod 755 update_functionlibrary 
-</code> 
-Or if you are working on the remote computer you can also do the following: 
-<code> 
-you@eulerX ~$ wget http://www.cp2k.org/_media/exercises:2015_ethz_mmm:update_functionlibrary.zip 
-you@eulerX ~$ unzip exercises:2015_ethz_mmm:update_functionlibrary.zip 
-you@eulerX ~$ chmod 755 update_functionlibrary 
-</code> 
  
-</note> 
  
-<note important> **REMEMBER** 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>
Line 27: Line 16:
 </code> </code>
  
-**and to submit the job chain:**+**and to submit the job:**
  
 <code> <code>
-you@eulerX ~$ bsub < c2h2.a.chain+you@eulerX ~$ bsub < jobname
 </code>  </code> 
 </note>  </note> 
Line 126: Line 115:
 # #
 # loads the functions # loads the functions
-# this script uses the functions: change replace getcolumn add list+# this script uses the functions: m_change m_replace m_getcolumn m_sum m_list
 # #
 module load new cp2k module load new cp2k
-. ~/Scripts/myfunctions.bash+. /cluster/apps/courses/mmm/m_functions.bash
  
  
Line 141: Line 130:
 # #
  
-for n in `list -5 5`+for n in `m_list -5 5`
 do do
-  dist_cc=`change $dist_cc0 $n $delta` +  dist_cc=`m_change $dist_cc0 $n $delta` 
-  c2=`add $dist_hc $dist_cc ` +  c2=`m_sum $dist_hc $dist_cc` 
-  h2=`add $c2 $dist_hc `+  h2=`m_sum $c2 $dist_hc `
   #   #
   # replaces coordinates in the input   # replaces coordinates in the input
   #   #
-  replace _C1_ $dist_hc < c2h2.amber.inp.templ | replace _C2_ $c2  | replace _H2_ $h2 > c2h2.$dist_cc.inp+  m_replace _C1_ $dist_hc < c2h2.inp.templ | m_replace _C2_ $c2  | m_replace _H2_ $h2 > c2h2.$dist_cc.inp
   #   #
   # Runs cp2k   # Runs cp2k
Line 159: Line 148:
   # gets the energy from the output (in a.u.)   # gets the energy from the output (in a.u.)
   #   #
-  energy=` getcolumn 'ENERGY|' 9 < c2h2.$dist_cc.out `+  energy=` m_getcolumn 'ENERGY|' 9 < c2h2.$dist_cc.out ` 
   echo $dist_cc $energy  >> curve.amber.c2h2   echo $dist_cc $energy  >> curve.amber.c2h2
 # #
Line 172: Line 161:
 mv c2h2.*.inp c2h2.*.out Logs mv c2h2.*.inp c2h2.*.out Logs
 rm RUN* *restart* rm RUN* *restart*
 +
 +
  
 </code> </code>
Line 179: Line 170:
  
 <code> <code>
-you@brutusX c2h2$ module load new cp2k +you@eulerX c2h2$ module load new cp2k 
-you@brutusX c2h2$ bsub < c2h2.chain+you@eulerX c2h2$ bsub < c2h2.chain
 </code> </code>
  
Line 188: Line 179:
  
 <code> <code>
-you@brutusX c2h2$ gnuplot fit.gnu+you@eulerX c2h2$ gnuplot fit.gnu
 </code> </code>
  
Line 262: Line 253:
 # #
 # loads the functions # loads the functions
-# this script uses the functions: change replace getcolumn add list+# this script uses the functions: m_change m_replace m_getcolumn m_sum m_list
 # #
 module load new cp2k module load new cp2k
-. ~/Scripts/myfunctions.bash+#. ~/Scripts/myfunctions.bash 
 +. /cluster/apps/courses/mmm/m_functions.bash
  
  
Line 277: Line 269:
 # #
  
-for n in `list -5 5`+for n in `m_list -5 5`
 do do
-  dist_cc=`change $dist_cc0 $n $delta` +  dist_cc=`m_change $dist_cc0 $n $delta` 
-  c2=`add $dist_hcx $dist_cc` +  c2=`m_sum $dist_hcx $dist_cc` 
-  h3=`add $c2 $dist_hcx `+  h3=`m_sum $c2 $dist_hcx `
   #   #
   # replaces coordinates in the input   # replaces coordinates in the input
   #   #
-  replace _C1_ $dist_hcx < c2h4.amber.inp.templ | replace _C2_ $c2  | replace _H3_ $h3 > c2h4.$dist_cc.inp+  m_replace _C1_ $dist_hcx < c2h4.inp.templ | m_replace _C2_ $c2  | m_replace _H3_ $h3 > c2h4.$dist_cc.inp
   #   #
   # Runs cp2k   # Runs cp2k
Line 294: Line 286:
   # gets the energy from the output (in a.u.)   # gets the energy from the output (in a.u.)
   #   #
-  energy=` getcolumn 'ENERGY|' 9 < c2h4.$dist_cc.out ` +  energy=` m_getcolumn 'ENERGY|' 9 < c2h4.$dist_cc.out `
   echo $dist_cc $energy  >> curve.amber.c2h4   echo $dist_cc $energy  >> curve.amber.c2h4
 # #
Line 307: Line 299:
 mv c2h4.*.inp c2h4.*.out Logs mv c2h4.*.inp c2h4.*.out Logs
 rm RUN* *restart* rm RUN* *restart*
-                                                                                                                                                                        + 
 +                                                                                                                                                  
 </code> </code>
  
exercises/2015_ethz_mmm/c2h2_bond_energy.txt · Last modified: 2020/08/21 10:15 by 127.0.0.1