Hi Jingagang,<div><br></div><div>I believe the issue is that although the constraint (shake/rattle) is in principle implemented for RESPA (src/motion/integrator.F) it's only applied to the inner loop of the algorithm which likely uses the high-low force_eval force difference as far as I can see from the code. </div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir<br> </div><div class="gmail_quote"><div dir="auto" class="gmail_attr">вторник, 26 октября 2021 г. в 22:40:46 UTC+2, Jinggang Lan: <br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">When the blue-moon ensemble combines with RESPA, the Lagrangian multipliers are incorrect. The SHAKE Lagrangian Multipliers from RESPA are much smaller than the normal constrained MD.<br></blockquote></div>

<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 on the web visit <a href="https://groups.google.com/d/msgid/cp2k/c46037e7-da3a-4f7b-9e8d-da8527a6ac6cn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/c46037e7-da3a-4f7b-9e8d-da8527a6ac6cn%40googlegroups.com</a>.<br />