I read the paper cited in the mmfp.doc and they have used Langevin dynamics to control the temperature of the system.
I will try rerunning my simulations using LD and update.

Thanks
Vidhya