Re: AMBER: Solvation net forces

From: David A. Case <>
Date: Tue, 18 Sep 2007 17:14:55 -0700

On Mon, Sep 17, 2007, Tyler Luchko wrote:
> I am implementing an implicit solvation model in Sander called 3D-
> RISM. The solvation forces for each solute atom are integrated off
> of a grid. By looking at several systems we have confirmed that we
> are getting the expected forces except for what appears to be a
> random drift force. I.e. there is a non-zero net solvation force on
> the solute. These forces scale linearly with the number of solute
> atoms and are roughly 0.1 kcal/mol/Angstrom in each direction. I
> would like to simply apply an equal and opposite force to the center-
> of-mass to remove the error.

This looks like a pretty large force error to me. Can you get rid of it
by decreasing the grid size? (Even if this is expensive, it seems like a good
test in order to make sure that there are not other bugs lying around.)

I don't have a particular opinion about subtracting it out, expect that it
*sounds* dangerous....


The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Wed Sep 19 2007 - 06:07:43 PDT
Custom Search