[ESPResSo-devel] Langevin thermostat (cool/heat)

[Jakub Krajniak]

Dear all,

I have a question on Langevin thermostat, mainly the heat/cool down part.

Looking into the code:
https://github.com/espressomd/espresso/blob/master/src/core/thermostat.cpp#L142

void thermo_heat_up()
{
if(thermo_switch & THERMO_LANGEVIN) {
langevin_pref2_buffer = langevin_pref2;
langevin_pref2_rotation_buffer = langevin_pref2_rotation;
langevin_pref2 *= sqrt(3);

Why pref2 is multiply by sqrt(3), is it some rule of thumb or is there 
any theory after that?

Could this heat/cool part cause some instability if the integrator is 
restarted too often?

Best regards,
Jakub

Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm