<div dir="ltr">Dear cp2k community,<br><br>I am new to cp2k and try to set up a QMMM system with the MM part gradually restraint, starting totally flexible close to the QM region and getting more and more stiff further away. So far I've just found the <br>CP2K_INPUT / MOTION / CONSTRAINT / COLLECTIVE / RESTRAINT subsection, where I can specify some harmonic force constant for some collective variable. <br>To restrain an atom I thought about the CP2K_INPUT /FORCE_EVAL / SUBSYS / COLVAR / XYZ_DIAG subsection to specify the atom and its component. With that it should be possible to assign individual restraints.<br>I thought about a setup like that:<br><br>Definition of the colvar<br><br>&FORCE_EVAL<br>  &SUBSYS<br>    &COLVAR<br>      &XYZ_DIAG<br>        ATOM 1<br>        COMPONENT X<br>      &END XYZ_DIAG<br>    &END COLVAR<br>  &END SUBSYS<br>&END FORCE_EVAL<br><br>and definition of the restraint<br><br>&MOTION<br>  &CONSTRAINT<br>     &COLLECTIVE<br>        COLVAR 1<br>        INTERMOLECULAR .TRUE.<br>        &RESTRAINT<br>           K 0.015<br>        &END RESTRAINT<br>     &END COLLECTIVE<br>  &END CONSTRAINT<br>&END MOTION<br><br>Is this the right way to go? If so, then it should be defined for each component of several hundred atoms. Is there an easier way to achieve this?<br><br>Many thanks<br><br><br>Florian<br><big class="uctt"><a href="https://manual.cp2k.org/cp2k-3_0-branch/CP2K_INPUT/MOTION/CONSTRAINT/COLLECTIVE.html"></a></big></div>