The issue if you do simulations at non-default temperature is that the particles diffusion coefficients scale with T according to the equations of Brownian However, The BrownianDynamics class relies on whatever diffusion coefficients were set by the Diffusion decorator, regardless of the simulation temperature. Therefore, to simulate correctly at non-default temperatures, one must do two things:
(1) Call BrownianDynamics::set_temperature()
(2) Update the diffusion coefficients for all particles to scale with temperature - e.g. D is 10% higher at 330 K than at 300 K.
I will update the IMP documentation to explain this subtlety. Let me know if you think we should change the code behavior. Perhaps we can add code to automatically update D, but it's also risky and could lead to unexpected behavior.
Barak