Re: [AMBER] Huge discrepancy in el part of "solvation energy" using Thermodynamic Integration

From: case <case.biomaps.rutgers.edu>
Date: Sat, 5 May 2012 09:45:09 -0400

On Wed, May 02, 2012, Marek Maly wrote:

> I repeated calculation of dG_bulk_elec of phenol using input files from
> tut. A9 and have obtained value -38.49138 kcal/mol which is almost
> identical with the previously reported value ( -38.64 see below) which
> is really different from the value reported in below mentioned article (
> -6.03 kcal/mol).

If dG_bulk_elec refers to phenol in bulk water, then this number is way wrong
(as you point out). It was kind of hard to follow the scripts to supplied,
but I don't see where you are computing the electrostatic change in vacuum:
this needs to be subtracted from the charging free energy in water to get the
electrostatic contribution to the solvation free energy.

[The "softcore" code makes this vaccum step unnecessary, but you don't seem
to be using that option(?)]

...hope this helps....dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sat May 05 2012 - 07:00:03 PDT
Custom Search