Re: AMBER: igb and ibelly

From: David A. Case <case.scripps.edu>
Date: Fri, 30 May 2008 13:20:52 -0700

On Fri, May 30, 2008, Holly Freedman wrote:

> Can someone please explain to me why it is that the ibelly option is not
> allowed with igb>0 (Generalized Born)? I don't understand how the two
> options are connected.

Since GB is not a pairwise theory, it is not a simple matter to simply try to
ignore certain atoms. A discussion of what would be involved, and the
approximations required, is given here:

%A O. Guvench
%A J. Weiser
%A P. Shenkin
%A I. Kolossvary
%A W.C. Still
%T Application of the Frozen Atom Approximation to the GB/SA Continuum Model
for Solvation Free Energy
%J J. Comput. Chem.
%V 23
%P 214-221
%D 2002

These ideas have never been implemented in Amber.

You can set ntr=1 and restrain atoms to be (very) near to their starting
points.

[Aside: I've long been in favor of removing ibelly for any purpose in Amber,
but it hasn't been a high priority, and people seem to have some use for that
facility, so it still sticks around.]

...dac

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
      to majordomo.scripps.edu
Received on Sun Jun 01 2008 - 06:08:04 PDT
Custom Search