Hi,<br>
<br>
I need some information about how Gromacs handles PME for free energy
calculations when sc-alpha is nonzero and the perturbation is
electrostatic. That is, how is PME handled when soft core potentials
are used to turn off electrostatics for a molecule in solution, for
example? For what it's worth, I'm using the CVS version of Gromacs 3.2
referenced here
(<a href="http://www.gromacs.org/pipermail/gmx-developers/2004-December/000982.html">http://www.gromacs.org/pipermail/gmx-developers/2004-December/000982.html</a>),
but I assume things are handled the same in Gromacs 3.3, which I'll be
switching to shortly.&nbsp; <br>
<br>
What I'd like to see is a derivation or description of a derivation for
what it means to correctly deal with PME for this case. The recent J.
Comp. Chem. paper on Gromacs (v. 26, n. 16) states that &quot;Differences in
constraints or in long-range interactions are<br>
properly handled,&quot; but doesn't explain what it means to handle these
properly in the soft core case. The papers referenced in that section
(one on Gromacs 3.0, and van Gunsteren's paper introducing the soft
core potential) are similarly unhelpful.<br>
<br>
In other words, I know the claim is that Gromacs now properly handles
electrostatics when a soft-core potential is used to turn off
electrostatics in a free energy calculation. But what does it mean to
&quot;properly handle&quot; this? How exactly are they handled? Again, I'd like
to see a derivation (or reference to one) or at the very least a
description of a derivation so I can work through it myself.<br>
<br>
Thanks!<br>
David Mobley<br>
Dill group<br>
UCSF<br>
<br>
<br>