<div dir="ltr"><div>Hello,</div><div><br></div><div>I know this is an old post, but I have the same issue and I haven't found any solution in this group. I manage to notice that all Lagrange multipliers are printed zero when I use Fixed Atoms as a constraint, but are printed normal when I constraint some Collective Variable. So, I'm starting to think that this is a bug of the program, because somewhere have to the information of the force acting on the fixed atoms (If I print the forces, are also zero). I will appreciate any help regarding this topic.<br></div><div><br></div><div>Best regards,</div><div><br></div><div>Alejandro.<br></div><div><br></div><div><br></div><br>El miércoles, 27 de mayo de 2015, 10:21:38 (UTC-4), brad...@gmail.com escribió:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hello!<div><br></div><div>I'm trying to calculate a PMF by integrating forces over a reaction coordinate defined by one interatomic distance (defined in the colvar/collective sections).  For the attached input, I've turned off all other constraints for testing purposes.  The Lagrange multiplier files print and, as expected, I get one number each for SHAKE and RATTLE, but they're both zero.</div><div><br></div><div>What am I doing wrong?</div><div><br></div><div>Thanks!</div><div><br></div><div>-Brad</div></div></blockquote></div>