> May it be just a velocity of the Langevin part of the AD_LANGEVIN thermostat slipping through the code and producing the bug? Could certainly be. I'll take a look, since I wasn't thinking about constraints when I implemented it. Noam