Hi,
>
> On an aside I believe a number of people now consider 2fs with shake to
> really be on the bleeding edge of what you can safely use. This is with
> regards to energy conservation etc when running in the NVE ensemble. For
> safety you should probably use a 1.5fs timestep even with shake.
>
...Just my 0.02$: Long ago I played around with droplet dynamics.
Though that is bad (except as a toy), I observed *no* discernible
drift of total energy over 7 nanoseconds with NVE and a 2 fs timestep.
The statistics can be seen here:
http://www.tbi.univie.ac.at/~svrci/sample.ps
(-> a page from my thesis)
So if there's a drift with 2 fs, it might not come from SHAKE
(or home-brewn RATTLE, in that particular case, which should be the same).
By the way: I did a RNA tetraloop (~400 atoms) + ~1800 TIP3P waters at
300K.
good luck
Andreas
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Jun 04 2006 - 06:07:07 PDT
