AMBER Archive (2008)

Subject: AMBER: energy leakage?

From: Sally Pias (sallypias_at_gmail.com)
Date: Thu Oct 23 2008 - 02:25:59 CDT

I have an explicitly solvated protein system for which I am attempting
to look at tens of nanosecond molecular dynamics using Sander.
Although the system appears to equilibrate well at first, I have been
unable to maintain constant temperature (at 300 K) upon switching from
ntt=3 to ntt=0. In previous simulations, the temperature increased
gradually from 298 K (after an initial dip) to 310 K over 30
nanoseconds. When I examined the energy profiles, I found that both
kinetic and potential energy are increasing throughout the ntt=0
simulations. I understand this to be an indicator of poor
equilibration, but I am unclear as to the cause of this sort of
problem in terms of the way the force field works.

Sally
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
      to majordomo_at_scripps.edu