[gmx-developers] [gmx-users] Wrong Nose-Hoover integrator
hessb at mpip-mainz.mpg.de
Tue Mar 29 12:12:02 CEST 2005
I have implemented the correct Nose-Hoover leap-frog scheme as written in:
But I have not been able to reproduce the incorrect 3.2.1 thermostat work
as reported by Michel.
For a box of 216 waters with delta_t=0.002 and 0.9 nm RF I get correct work
at least up to tau_t=0.1 ps.
Michel, how many waters, what cut-off scheme and what time step did you use?
The correct leap-frog algorithm (or probably any reversible integration
poses another problem.
With dt=2 fs SPC water with a full step tempature of 300 K (determined
has a half step temperature of 297.8 K. With dt=4 fs the difference will
be much larger.
The full-step temperature is the more accurate expression, but can not
in any simple integration scheme. This means that the temperature for water
simulations with reversible Nose-Hoover coupling will be off by at least
I am afraid this problem can not be solved.
More information about the gromacs.org_gmx-developers