Dear Steven, <div><br /></div><div>There should be no problem in applying constraints in a QM/MM simulation,</div><div>but if you are using the interface between Gromacs and CP2K, the dynamics is driven by Gromacs</div><div>not by CP2K. </div><div><br /></div><div>Regards</div><div>Marcella</div><div><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, June 24, 2024 at 1:21:47 PM UTC+2 steve...@gmail.com wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi all,<div><br></div><div>I am currently trying to run a QM/MM simulation using CP2K with Gromacs. Running the simulation alone works fine, but I would like to add a distance restraint or constraint to pull two QM modelled atoms together. Below is the relevant input parameters i have used, but no pulling effect is observed (Note that the numbers have been set very high just for testing purposes). I'm unsure if it is possible to constrain a system when performing a FORCE_EVAL?</div><div><br></div><div>Any ideas or assistance from the community would be most appreciated.</div><div><br></div><div> &COLVAR<br> &DISTANCE<br> ATOMS 33293 3157<br> &END DISTANCE<br> &PRINT<br> &END<br> &END COLVAR<br> &END SUBSYS<br>&END FORCE_EVAL<br>&MOTION<br> &CONSTRAINT<br> &CONSTRAINT_INFO<br> COMMON_ITERATION_LEVELS 3<br> FILENAME COLVAR_CONSTR.OUT<br> &EACH<br> MD 1<br> &END<br> &END<br> &COLLECTIVE<br> COLVAR 1<br> INTERMOLECULAR T<br> TARGET [angstrom] 3.9<br> TARGET_GROWTH [angstrom*fs^-1] -50.00<br> TARGET_LIMIT [angstrom] 2.0<br> &END<br> &END<br>&END<br></div><div><br></div><div><div>Many thanks,</div><div>Steven</div></div></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/8bbd63f1-6908-439d-840d-f3440eafb1acn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/8bbd63f1-6908-439d-840d-f3440eafb1acn%40googlegroups.com</a>.<br />