1 / 12

Running molecular dynamics with constraints included

Running molecular dynamics with constraints included. Restraint or Constraint Dynamics?. Constraint: Must be satisfied (eg fixed bond length) Restraint: penalties are imposed for a derivation from the reference value (eg increased force constant) (Chapter 9.2).

huy
Download Presentation

Running molecular dynamics with constraints included

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Running molecular dynamics with constraints included

  2. Restraint or Constraint Dynamics? Constraint: Must be satisfied (eg fixed bond length) Restraint: penalties are imposed for a derivation from the reference value (eg increased force constant) (Chapter 9.2)

  3. What is Constraints Dynamics? One of more internal coordinate is kept fixed during the computation. Why Constraints Dynamics? The time step is dictated by the highest frequency. Therefore, using constraint can reduce the time step and make the computation faster. It reduces the numbers of degrees of freedom to 3N - k .

  4. Homolytic Constraints ƒ(q1,q2,q3,…,t)=0 eg. r2 - a2 = 0 Non-Homolytic Constraints ƒ(q1,q2,q3,…,t)>0 eg. r2 - a2 > 0 The SHAKE Algorithm uses Homolytic Constraints.

  5. dij 1 2

  6. The Verlet Algorithm gives (from last week): eqn 7.8 Applying the expressions for the forces in the Verlet Algorithm... Combining both...

  7. For the case with two atoms: What else do we know? This is a problem with 3 equations and 3 variables ( r1(t+δt), r2(t+δt), λ12 ) !!!

  8. Let After painful manipulation of variables... This can be solved algebraically! However, this is not always the case...

  9. For the case with 3 atoms, we get the following equations: This will become more complex very rapidly!

  10. One way of solving this system of equation is to drop the λ2 terms, to form linear equations with respect to λ. We then get a k x k matrix, were k is the number of constraints. The SHAKE algorithm doesn’t do that... It satisfies one constraint at one time and reiterate around the constraints until a convergence is obtained. The tolerance needs to be tight enough that the fluctuations in this algorithm don’t affect the remainder of the simulation.

  11. So, What’s the catch? You need to make sure that the constraints don’t prevent a full exploration of all the possibilities (eg by preventing a rotation or a torsion). In other words, the constrained degrees of freedom must be weakly coupled to the remaining degrees of freedom.

  12. But I don’t want to use the Verlet Algorithm... This method has also been applied to the other algorithms, ie: Velocity Verlet (renamed to RATTLE) leap-frog predictor-corrector What about angles and torsions? Angles can be constrained with an additional distance constraint (eg SPC water). There’s also a method by Tobias and Brooks that enable constraints to be applied to arbitrary internal coordinates.

More Related