For several reasons (drift during equilibration, drift as a result of force truncation and integration errors, heating due to external or frictional forces), it is necessary to control the temperature of the system.
The Berendsen algorithm(195) mimics weak coupling with first-order kinetics to an external heat bath with given temperature T0. The effect of this algorithm is that a deviation of the system temperature from T0 is slowly corrected according to