<div dir="ltr"><div>Dear CP2K community,</div><div><br></div><div>I'm performing a constrained MD simulation to use later the Lagrange Multipliers to compute the force along some reaction path. I'm constraining the distances between all atoms in a 33 atoms molecule. However, if I compute the norm of the vector formed by all Lagrange Multipliers in each step, I get the behavior shown in the plot attached. This norm monotonically increases during the simulation and reset itself when I continue the run using a restart file, and of course  you can see the same pattern replicated if I compute the force along the reaction path. I expect that the norm of the Lagrange Multipliers should show some equilibrated behavior and not have this jumps every time I continue the simulation with the restart file. Does anyone have a clue why is this happening? What could be wrong in the input file? <br></div><div><br></div><div>I'm attaching the graph showing the norm of the Lagrange Multipliers, time scale is in picoseconds, and Lagrange Multipliers units are in Hartree/Bohr.</div><div><br></div><div>I'm also attaching the input file.</div><div><br></div><div>Thanks in advance,</div><div><br></div><div>Alejandro.<br></div></div>