<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>