AMBER: How do you evaluate forces on a fixed atom?

From: Bill K <>
Date: Mon, 1 Dec 2008 13:45:42 -0500 (EST)

Hi Everybody,
     Has anyone tried to extract the systematic forces on a fixed
atom, or may know how to do it? I'm using amber7 and, if you have an
atom fixed in the system, the forces on it are canceled before they
can be evaluated. That is to say, output prints values of 0.0 0.0 0.0
for the force on a fixed atom, even though it experienced forces from
the surrounding atoms. Looking at the code, this occurs because
ibelly sets these values to zero. It looks like this task is not
possible without changing the source code. Or, perhaps this is
something that has been worked around in the more recent amber
releases. Any input would be greatly appreciated, as I've looked
through the archives on here and haven't seen this discussed before.
Thank you very much for your time and consideration.

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
Received on Fri Dec 05 2008 - 18:23:42 PST
Custom Search