<p>I am trying to run a CP2K molecular dynamics simulation where specific bond lengths should remain fixed. Based on the documentation, SHAKE should be used, but I keep encountering errors when defining it under &CONSTRAINT or &MD.</p><p>Here is my attempt:</p><div><div>```</div><div><div><div><span aria-label="Copy">Copy</span>Edit</div></div></div><div dir="ltr">&MOTION
<br />         &CONSTRAINT
<br />                &SHAKE<br />                      ATOMS 1 2<br />                      TYPE BOND
<br />                &END SHAKE
<br />         &END CONSTRAINT
<br />&END MOTION<br />```</div></div><p>However, CP2K returns an error:</p><div><div><br /></div><div dir="ltr">``` unknown subsection SHAKE of section CONSTRAINT```</div></div><p>I also tried defining SHAKE within &MD, but CP2K does not recognize that either.</p><p>What is the correct way to fix bond lengths using SHAKE in CP2K? If SHAKE is not supported in the latest versions, what alternative approach should I use?</p><p>Thank you for any suggestion and help</p>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/256bce8e-ceee-4151-83be-e729b198abban%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/256bce8e-ceee-4151-83be-e729b198abban%40googlegroups.com</a>.<br />