CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace: • gga • defects_in_graphene • nudged_elastic_band • t_melting
exercises:2015_pitt:gga

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
exercises:2015_pitt:gga [2015/02/26 14:29] – vondeleexercises:2015_pitt:gga [2020/08/21 10:15] (current) – external edit 127.0.0.1
Line 23: Line 23:
   * Use vmd to vizualize the geometries (provided below) named ''mode1.xyz'' and ''mode2.xyz''   * Use vmd to vizualize the geometries (provided below) named ''mode1.xyz'' and ''mode2.xyz''
   * To edit the input files provided below, use an editor such as ''vi'' or ''nano''. While ''nano'' is simple to use, ''vi'' can be [[http://www.cp2k.org/tools| configured to colour-code cp2k inputs]].   * To edit the input files provided below, use an editor such as ''vi'' or ''nano''. While ''nano'' is simple to use, ''vi'' can be [[http://www.cp2k.org/tools| configured to colour-code cp2k inputs]].
-  * You will need files named ''BASIS_MOLOPT'' ''GTH_POTENTIALS'' ''dftd3.dat'' that are provided as part of CP2K in a directory ''cp2k/data'' +  * You will need files named ''BASIS_MOLOPT'' ''GTH_POTENTIALS'' ''dftd3.dat'' that are provided as part of CP2K in a directory ''cp2k/data'', unless the code has been compiled with the proper flag ('' -D__DATA_DIR'') so that these are found automatically. 
-  * Use a job script to submit jobs on the cluster, an example job submission script might look like (TODO: adjust for environment)+  * Use a job script to submit jobs on the cluster, an example job submission script might look like 
  
 <code - job> <code - job>
-#!/bin/bash --login +#PBS -N mode1 
-#SBATCH --job-name=mode1 +#PBS -j oe 
-#SBATCH --nodes=6 +#PBS -q dist_small 
-#SBATCH --time=0:10:00 +#PBS -l nodes=4:ppn=16 
-#SBATCH --account=ch5 +#PBS -l walltime=10:00 
-aprun -n 144 -N 24 -d 1 cp2k.popt -i mode1.inp -o mode1.out+#PBS -A cp2k2015 
 + 
 +cd $PBS_O_WORKDIR 
 + 
 +module purge 
 +module load cp2k/2.6 
 + 
 +prun cp2k.popt -i mode1.inp -o mode1.out
 </code> </code>
  
Line 89: Line 96:
 ===== 3. Task: relative stabilities  ===== ===== 3. Task: relative stabilities  =====
  
-TODO+In order to compute the relative stability of mode1 and mode2, both configurations need to be geometry optimized. 
  
 +To do so, turn off the generation of cubes (''&E_DENSITY_CUBE OFF'') in mode1.inp, change to ''RUN_TYPE GEO_OPT'' and adjust the project name. Create and run a similar input file for mode2.
 +
 +input topics:
 +  * BFGS vs LBFGS
 +  * EPS_SCF, CUTOFF, MAX_DR, ..
 +
 +output topics:
 +  * ''Informations at step''
 +  * Trajectory ''MODE1_GEO-pos-1.xyz''
 +
 +Compare the final energies (''ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):''), and determine which mode is most stable. Does this agree with the values in table 1 of the manuscript cited ?
 ===== 4. Task: ab initio molecular dynamics  ===== ===== 4. Task: ab initio molecular dynamics  =====
  
-TODO+<note>This task is optional, and can be done near the end if time is available. Ab initio MD will be briefly discussed in a next exercise.</note> 
 + 
 +Perform a short ab initio molecular dynamics simulation of the system (~1000 steps, ~0.5ps) by changing to ''RUN_TYPE MD''. After a couple of hours the job should be finished. Now analyze the OH distance in VMD. A possible outcome is : 
 + 
 +{{exercises:2015_pitt:acetic_acid_mode1_md.png?600 |}} 
 + 
 +What can you say about the hydrogen bond to the surface, relative acidity of the two oxygens ?  
 +Note that, in order to be statistically relevant, longer trajectories should be employed, and surface slab thickness will play an important role. Also compare to Fig. 7 of the paper referenced. 
 +====== Required Files ====== 
 + 
 +(right) click on the filename to download to your local machine.
  
-===== Required Files ===== 
 <code - mode1.inp> <code - mode1.inp>
 &GLOBAL &GLOBAL
Line 118: Line 145:
     &MGRID     &MGRID
        ! PW cutoff ... depends on the element (basis) too small cutoffs lead to the eggbox effect.        ! PW cutoff ... depends on the element (basis) too small cutoffs lead to the eggbox effect.
-       ! certain calculations (e.g. geometry optimization, vibrational frequencies, NPT and cell optimizations, need higher cutoffs)+       ! certain calculations (e.g. geometry optimization, vibrational frequencies, 
 +       ! NPT and cell optimizations, need higher cutoffs)
        CUTOFF [Ry] 400         CUTOFF [Ry] 400 
     &END     &END
  
     &QS     &QS
-       METHOD GPW  ! use the GPW method (i.e. pseudopotential based calculations with the Gaussian and Plane Waves scheme). +       ! use the GPW method (i.e. pseudopotential based calculations with the Gaussian and Plane Waves scheme). 
-       EPS_DEFAULT 1.0E-10 ! default threshold for numerics ~ roughly numerical accuracy of the total energy per electron, sets reasonable values for all other thresholds +       METHOD GPW  
-       EXTRAPOLATION ASPC ! used for MD, the method used to generate the initial guess. +       ! default threshold for numerics ~ roughly numerical accuracy of the total energy per electron, 
 +       ! sets reasonable values for all other thresholds. 
 +       EPS_DEFAULT 1.0E-10  
 +       ! used for MD, the method used to generate the initial guess. 
 +       EXTRAPOLATION ASPC 
     &END     &END
  
Line 182: Line 214:
     &END CELL     &END CELL
  
-    ! atom coordinates can be in the &COORD section, or provided as an external file.+    ! atom coordinates can be in the &COORD section, 
 +    ! or provided as an external file.
     &TOPOLOGY     &TOPOLOGY
       COORD_FILE_NAME mode1.xyz       COORD_FILE_NAME mode1.xyz
Line 188: Line 221:
     &END     &END
  
-    ! MOLOPT basis sets are fairly costly, but in the 'DZVP-MOLOPT-SR-GTH' available for all elements +    ! MOLOPT basis sets are fairly costly, 
-    ! their contracted nature makes them suitable for condensed and gas phase systems alike.+    ! but in the 'DZVP-MOLOPT-SR-GTH' available for all elements 
 +    ! their contracted nature makes them suitable 
 +    ! for condensed and gas phase systems alike.
     &KIND H                                   &KIND H                              
       BASIS_SET DZVP-MOLOPT-SR-GTH                BASIS_SET DZVP-MOLOPT-SR-GTH         
Line 222: Line 257:
    TEMPERATURE [K] 300    TEMPERATURE [K] 300
    TIMESTEP [fs] 0.5    TIMESTEP [fs] 0.5
-   STEPS 100+   STEPS 1000
    # GLE thermostat as generated at http://epfl-cosmo.github.io/gle4md     # GLE thermostat as generated at http://epfl-cosmo.github.io/gle4md 
    # GLE provides an effective NVT sampling.    # GLE provides an effective NVT sampling.
exercises/2015_pitt/gga.1424960964.txt.gz · Last modified: (external edit)