[CP2K-user] [CP2K:21261] How to Fix Bond Lengths Using SHAKE in CP2K?

[Planck]

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.

Here is my attempt:
```
CopyEdit
&MOTION 
         &CONSTRAINT 
                &SHAKE
                      ATOMS 1 2
                      TYPE BOND 
                &END SHAKE 
         &END CONSTRAINT 
&END MOTION
```

However, CP2K returns an error:

``` unknown subsection SHAKE of section CONSTRAINT```

I also tried defining SHAKE within &MD, but CP2K does not recognize that 
either.

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?

Thank you for any suggestion and help

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/256bce8e-ceee-4151-83be-e729b198abban%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20250306/c7b8e4ee/attachment.htm>