On Fri, Jan 21, 2005, Arantxa Sanz wrote:
>
> I have been running the following Sander script to run MD dynamics at
> constant T with Andersen thermostat
> Discarding the first 200ps, it seems that average T oscillates
> around 310K. That is, 10K above the given target T.
Try setting vrand to a lower value than its default value of 1000. Also,
run some tests with a small value of ntpr and look closely at the temperature
profile....maybe that will give some hints.
Also: be sure that you really want/need this sort of thermostat. It is
implemented in only a very naive way.
....dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sat Jan 22 2005 - 00:53:00 PST
