<div dir="ltr">Dear all,<div><br></div><div>I am trying to run Quickstep MD with some fixed atoms positions (set in MOTION / CONSTRAINT / FIXED_ATOMS), but the positions do change over time.</div><div>Curiously enough, constrained atoms have zero velocity at frame 0 and zero force at every timestep, which would suggest that the constraint is somehow applied but ineffective.</div><div><br></div><div>This behaviour seems to be machine and version independent:</div><div>tested cp2k versions: 2.7 (Development Version) git: 0826abf, 2.6 (Development Version) svn:14226, 2.5.1 svn:13632 </div><div>tested architectures/compilers: gfortran 4.8 on Darwin for cp2k 2.7, xlf 14.1 on BlueGene/Q for cp2k 2.5 and 2.6</div><div><br></div><div>The attachments are a minimal example of the problem, showing 5 timesteps of dynamics of 3 atoms of Si of which atom 1 is (or rather, should be) fixed.</div><div><br></div><div>Is this problem known, and is there any workaround? I know there is for GEO_OPT runs as in https://groups.google.com/forum/#!topic/cp2k/KAMhME81IbM, but I haven't found anything applicable to MD</div><div><br></div><div>Kindest regards,</div><div><br></div><div>Marco</div><div><br></div></div>