bd giving incorrect results, such as a temperature 2x as high as it should be in release-4-5-patches
Crashes with bd -- when I reduce the time step, I can get it to run, but the temperature is generally about 2x what it should be. This appears to have been the case since at least January.
fixed the reported temperature with Brownian dynamics
Reduced the (meaningless) masses with BD by a factor 2.
Since v with BD is defined as dx/dt, the reported Ekin and T are
now correct (and if there are too high, dt is too large).
#1 Updated by Berk Hess over 9 years ago
It should crash with such a timestep/friction ratio.
The reported temperature always seems to have been a factor 2 too high.
This is because the velocity is set to the displacement divided by the time step
and the mass is set to the friction coefficient times the time step.
As both velocity and mass are meaningless, there are two solutions:
1) set the velocity to 0.5*displacement/dt
2) set the mass to 0.5*fric*dt
The velocity has no meaning, but I think it makes more sense to set that equal to the displacement and therefore use 2). I can do this when we agree.
Note that all this has no effect on the simulations, only on the reported ekin and T.