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>

<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/2af2ec1a-183a-442b-b389-2e6a26de0305n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/2af2ec1a-183a-442b-b389-2e6a26de0305n%40googlegroups.com</a>.<br />