Re: AMBER: 1 fs timestep for GB?

From: Andreas Svrcek-Seiler <>
Date: Tue, 25 Nov 2003 19:32:15 +0100 (CET)

> I noticed that people tend to use 1 fs time-step for integration in GB/SA
> simulations. I understand that the scheme implemented does not involve
> friction with the environment and this can give rise to faster
> motions in some degrees of freedom of the solute. I am not sure how much
> faster bonds fluctuate in this case, but a 2 fs integration step
> did not produce anything strange for my nucleic acid simulation:
....I believe a 2 fs timestep won't hurt (at least not more than in other
simulations). The fast motions frozen
by SHAKE or RATTLE are bound to be similar in GB/SA and
other types of simulations.

> no drift in total energy over 10 ns - although the simulation was
> run at constant temperature (300 K) with tautp=5.0. check whether a large timestep/finite cutoff/error in the gradient
leads to drifts, tautp= *say* 1.e9
(or -better- no heat bath at all) would be more appropriate.
The thermostat masks drifts quite efficiently.

Hope it helps,
good luck

           ( O O )
              o Wolfgang Andreas Svrcek-Seiler
              o (godzilla)
      .oooO Tel.:01-4277-52733
      ( ) Oooo.
-------\ (----( )--------------------------------------------------------
        \_) ) /

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Tue Nov 25 2003 - 18:53:01 PST
Custom Search