<div dir="ltr"><br><div dir="ltr">Hi again, I have an update to my post. So far I figured out, that restraining atoms is quite easy:<br>The set up equals that for fixed atoms, just in addition a restraint constant has to be defined.<br><br>&MOTION<br>  &CONSTRAINT<br>     &FIXED_ATOMS<br>        COMPONENTS_TO_FIX XYZ<br>        LIST 1 ...<br>        &RESTRAINT<br>           K 0.001<br>        &END RESTRAINT<br>     &END FIXED_ATOMS<br>  &END CONSTRAINT<br>&END MOTION<br><br>But looking at my trajectory I realized a strange movement of the restrained region compared to the flexible and the totally fixed one. I think that the restrained region is really restrained to the initial position of those atoms, while the fixed and flexible part is translating/rotating in space, which I think could be caused by the QM region, which is fit into the QM cell, thus moving the rest (not the restraint reference positions!) with it.<br><br>Is this behaviour correct or is it a bug?<br><br>Florian<br></div></div>