Dear Amber users,
This is a follow up question to a thread in 2007 regarding NVE simulations (
http://archive.ambermd.org/200708/0140.html). I had
the same problem with energy drifting. Thus, I followed the recommendations
from Case and Cheatham; lower shake tolerance, increase the accuracy of the
PME by reducing DSUM_TOL by an order of magnitude, and dt=0.001. The
simulation is much "smoother" now in terms of energy and temperature. The
shift has now been reduced to 5 kcal/mol (-75489 to -75484, first to last
step) during a 10ns simulation. The energy does not yet seem to converge.
Should I be concerned about this energy drift, or is it small enough such
that I can accept it? My apologies if this has been answered earlier.
Best regards,
Lars Skjærven
Department of Biomedicine
University of Bergen
Norway
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Mar 09 2008 - 06:07:31 PDT
