[gmx-users] Energy too high for free energy

[Ing. Vojtěch Spiwok]

Hi gromacs users / developers

I'm trying to estimate a free energy from MD of the
complex of large protein (> 1000 AA) and small ligand
in order to compare free energies of two binding modes
which are different in conformation but same in topology.
This system is huge (> 80 000 atoms) and the potential
energy in absolute magnitude is high (- 10^6 kJmol-1).
When trying g_energy i recieve "inf" instead of free
energy estimate. I would like to ask wether this is due
to the fact that E/kT is too low so e^(E/kT) is close to
zero (out of precision of calculation) ? Than it would be
possible to calculate it from the output file. Is it correct
to set potential energy (# 9 in g_energy menu) for such
free-energy estimation ?
Thank you very much for previous answers.

Vojtech Spiwok
ICT Prague

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20030529/fc1eb6c5/attachment.html>