<br><br><div><span class="gmail_quote">2007/12/19, Konrad Hinsen &lt;<a href="mailto:hinsen@cnrs-orleans.fr">hinsen@cnrs-orleans.fr</a>&gt;:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
On Dec 19, 2007, at 14:27, Andreas Kring wrote:<br><br>&gt; Having a look at this paper by Kneller and Mülders it is seen (if I<br>&gt; understand it correctly) that more than one set of equations of<br>&gt; motions<br>&gt; are derived:
<br>&gt;<br>&gt; 1)<br>&gt; A set of equation derived by Lagrangians, projector techniques and<br>&gt; Dirac&#39;s theory of constrained Hamitonian dynamics.<br>&gt;<br>&gt; 2)<br>&gt; A set of equations derived as generalized Euler equations for linked
<br>&gt; rigid bodies.<br>&gt;<br>&gt; My question is: which of these sets of equations are used by MMTK to<br>&gt; perform MD?<br><br>The first one. In the absence of distance constraints, it is<br>equivalent to the Andersen barostat. In the absence of a barostat, it
<br>is equivalent to the standard SHAKE/RATTLE algorithm.<br></blockquote></div><br>Ok, thanks a lot. I&#39;m not that familiar with the SHAKE/RATTLE algorithm, so I will have a closer look at this. It&#39;s nice to know how things work :-)
<br><br>Kind regards<br>Andreas<br>