<div dir="ltr">Hi all,<div><br></div><div>I am posting this question on behalf of a colleague.</div><div><br></div><div>He is performing a CI-NEB calculation on copper substrates. The substrate is quite thin within the simulation cell, and has to be fixed in position during the calculation, otherwise the forces on the substrate atoms will be too large.</div><div><br></div><div>The CI-NEB calculation has 7 replicas, in which some atoms are fixed using CONSTRAINTS subsection in MOTION section. No initial velocities were set.  My colleague noticed the fixed atoms moved in replicas 2 to 7 at the first step, while the fixed atoms in replica 1 did not move (as intended). In the next steps, the atoms were fixed in all replicas.</div><div><br></div><div>It seems that this behaviour did not happen with CP2K 2.3 but is observed in versions 2.4 and 2.5.</div><div><br></div><div>Furthermore, running the calculation on a single processor gave the same behaviour as running on multiple processors, so this may not be a parallelisation issue.</div><div><br></div><div>I had a look at the source code, and it seems to me that the constraint information isn't even used for the default POT_OPT == FULL setting, and only COLVARs are used for POT_OPT == FE or ME.</div><div><br></div><div>So, the question is:</div><div><br></div><div>Can we perform a standard CI-NEB calculation with some atoms fixed in each replica?  And if so, what is the correct way of doing such calculations with CP2K?</div><div><br></div><div>We want to make sure we are doing things correctly before reporting this as a possible bug.</div><div><br></div><div>I have attached the restart file, which contains the CP2K input.</div><div><br></div><div>Many thanks and best regards,</div><div><br></div><div>Lianheng</div><div><br></div></div>