CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace: • simple_stm • hybrid
exercises:2018_uzh_acpc2:mol_sol

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:2018_uzh_acpc2:mol_sol [2018/05/04 14:09] – jglanexercises:2018_uzh_acpc2:mol_sol [2020/08/21 10:15] (current) – external edit 127.0.0.1
Line 4: Line 4:
 [[http://www1.lsbu.ac.uk/water/water_models.html|Water molecular models]] are computational techniques that have been developed in order to help discover the structure of water. In this section, you will be asked to calculate some physical properties based on classical molecular dynamics simulation. The TIP3/Fw model will be usded in the simulations.  [[http://www1.lsbu.ac.uk/water/water_models.html|Water molecular models]] are computational techniques that have been developed in order to help discover the structure of water. In this section, you will be asked to calculate some physical properties based on classical molecular dynamics simulation. The TIP3/Fw model will be usded in the simulations. 
  
-We have prepared a CP2K input file ''water.inp'' for running a MD simulation of liquid water using the force field from the first exercise (parametrized by [[https://aip.scitation.org/doi/pdf/10.1063/1.1884609|Praprotnik et al.]]). Download {{ water.zip |water.zip}}+We have prepared a CP2K input file ''water.inp'' for running a MD simulation of liquid water using the force field from the first exercise (parametrized by [[https://aip.scitation.org/doi/pdf/10.1063/1.1884609|Praprotnik et al.]]). 
 +<note important>Download {{ water.zip |water.zip}} and extract it.</note> 
  
  
Line 32: Line 33:
   * Plot RMSD for the water at 300K and calculate corresponding diffusion coefficient, are they expected?   * Plot RMSD for the water at 300K and calculate corresponding diffusion coefficient, are they expected?
 </note> </note>
 +<note important>The diffusion coefficient is calculated using MSD but NOT RMSD.</note>
  
 We will compute the vibrational spectrum, and dielectric constant of water based on molecular dynamics. The spectra for water are available in this paper [[https://aip.scitation.org/doi/pdf/10.1063/1.1884609|https://doi.org/10.1063/1.1884609]]. The provided program computes the correlation function of the (derivative of) the dipole moment and performs the Fourier transform.  We will compute the vibrational spectrum, and dielectric constant of water based on molecular dynamics. The spectra for water are available in this paper [[https://aip.scitation.org/doi/pdf/10.1063/1.1884609|https://doi.org/10.1063/1.1884609]]. The provided program computes the correlation function of the (derivative of) the dipole moment and performs the Fourier transform. 
Line 65: Line 67:
   * Does IR or dielectric constant match the experimenal value? If not, why?   * Does IR or dielectric constant match the experimenal value? If not, why?
 </note> </note>
- 
  
 ===== Ramachandran plot ===== ===== Ramachandran plot =====
Line 84: Line 85:
 Visualize the structure ''glyala.pdb'' with VMD and determine the atomic indices of the atoms defining the dihedral angles. Visualize the structure ''glyala.pdb'' with VMD and determine the atomic indices of the atoms defining the dihedral angles.
 </note> </note>
 +
 +<note important>//Note:// While VMD starts counting atoms from 0, CP2K starts counting from 1, i.e. the VMD indices need to be increased by 1.</note>
 +
  
 With this knowledge at hand, With this knowledge at hand,
Line 90: Line 94:
 <note>**TASK 2** <note>**TASK 2**
  
-  - The atomic indices defining the dihedral indices in the input file ''geo.in'' are missing. Replace ''I1'' to ''I4'' by the atomic indices determined previously. //Note:// While VMD starts counting atoms from 0, CP2K starts counting from 1, i.e. the VMD indices need to be increased by 1.+  - The atomic indices defining the dihedral indices in the input file ''geo.in'' are missing. Replace ''I1'' to ''I4'' by the atomic indices determined previously. 
   - Use ''perform-gopt.sh'' to perform the grid of geometry optimizations.   - Use ''perform-gopt.sh'' to perform the grid of geometry optimizations.
   - Use gnuplot to plot the potential energy surface (we have provided a script ''epot.gp''). Which are the two most favoured conformations? <code> $ gnuplot</code><code> gnuplot > load "epot.gp"</code>   - Use gnuplot to plot the potential energy surface (we have provided a script ''epot.gp''). Which are the two most favoured conformations? <code> $ gnuplot</code><code> gnuplot > load "epot.gp"</code>
 </note> </note>
- 
 ===== Glyala in water ===== ===== Glyala in water =====
 Now, we will move to a more realistic system - Glyala in water. We will preformed a MD of glyala in water and save the trajectory. Now, we will move to a more realistic system - Glyala in water. We will preformed a MD of glyala in water and save the trajectory.
exercises/2018_uzh_acpc2/mol_sol.1525442964.txt.gz · Last modified: (external edit)