I'm using a langevin thermostat in my simulation and I noticed that my test particle was bouncing around at a higher rate than I would have expected due to the noise term. I'm am trying to get some test particles to stick to a surface (the experimentalist tell me they should stick), but looking at the video output I see the particles just bouncing off the surface.
I noticed that the formula to calculate the noise coefficient in core/thermostat.cpp is:
sqrt(24 * temperature * langevin_gamma / time_step)
sqrt(2 * temperature * langevin_gamma / time_step)
A 2 instead of a 24. I changed this, recompiled and now I can see my particles sticking to the surface. I was just wanting to confirm if the 24 in the code is correct or should that be a 2?
I'm not sure if it's relevant but I have the MASS module enabled and my particle masses are around ~1000. I know that will change the units of time in my simulation, does the noise generation function account for this change?
Thanks,
David