CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace:
exercises:2017_uzh_acpc2:prot_fol

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_uzh_acpc2:prot_fol [2017/05/17 12:37] – vrybkinexercises:2017_uzh_acpc2:prot_fol [2017/05/17 12:48] – [Background] vrybkin
Line 3: Line 3:
  
 ===== Background ===== ===== Background =====
-A model protein you will have to deal with is the alanine decapeptide. The folding/unfolding will be achieved by fixing the distance between the end carbon atoms in the chain: atoms 7 and 98. This distance is called a collective variable. At each distance one runs the MD simulation (constrained MD) to extract the time-averaged forces acting on the collective variable, $F(x)$. Then, a free energy difference can be calculated via thermodynamic integration (TI):+A model protein you will have to deal with is the alanine decapeptide. The folding/unfolding will be achieved by stretching/compressing the chain and fixing the distance between the end carbon atoms in it: atoms 7 and 98. This distance is called a collective variable. At each distance one runs the MD simulation (constrained MD) to extract the time-averaged forces acting on the collective variable, $F(x)$. Then, a free energy difference can be calculated via thermodynamic integration (TI):
  
 \begin{equation} \begin{equation}
Line 54: Line 54:
 </code> </code>
  
-===== Task 3: Evaluate the free energy difference and plot unfolding energy profile =====+===== Task 3: Evaluate the free energy difference =====
 ⇒ Each constrained MD will produce a ''.LagrangeMultLog''-files, which look like this: ⇒ Each constrained MD will produce a ''.LagrangeMultLog''-files, which look like this:
 <code> <code>
Line 62: Line 62:
 Rattle Lagrangian Multipliers:            -0.318145579 Rattle Lagrangian Multipliers:            -0.318145579
 </code> </code>
 +
 +<note warning>
 +Make sure that you get the units right. The Largange multipliers are written in atomic units (Hartree/bohr), while the distances are in Angstrom.
 +</note>
  
   * From these files you can calculate the average Lagrange multiplier of the Shake-algorithm like this:   * From these files you can calculate the average Lagrange multiplier of the Shake-algorithm like this:
Line 71: Line 75:
   * From these forces the free energy difference can be obtained via TI (see **Background**)   * From these forces the free energy difference can be obtained via TI (see **Background**)
  
-<note warning> + 
-Make sure that you get the units right. The Largange multipliers are written in atomic units (Hartree/bohr), while the distances are in Angstrom.+<note tip> 
 +  * Calculate $\Delta A$ numerically using the trapezoidal rule (or equivalent) with EXCEL, ORIGIN or any scripting language.
 </note> </note>
 +
  
exercises/2017_uzh_acpc2/prot_fol.txt · Last modified: 2020/08/21 10:15 by 127.0.0.1