[gmx-users] How to calculate work

[Sai Kumar Ramadugu]

Hi Sanku,
I was using the pullf.xvg and multiplying it with pulling rate and time.
f*v*dt = W
and getting the total work for each SMD simulation.

I'm not sure if this is the best/correct way to do it. But from original
Jarzynski's article (PRL (2007) 78(14), 2690-2693) this is what I deduced.
I asked this question but I think very few users of Gromacs do Jarzynski
Equality (JE) to get free energy differences.

If anyone can comment on this, it will be useful for the community.

PS: If you use AMBER, it gives the work profile as the output. You can use
the total work which is printed at the last line of the output and get the
exponential average of beta*W.

PPS: May be you are aware of this issue, but the JE suffers from the fact
that you do the exponential average and the smaller work values determine
everything.


Regards
Sai



On Fri, Nov 4, 2011 at 7:46 PM, Sanku M <msanku65 at yahoo.com> wrote:

> Hi,
> I am performing steered MD simulation using gromacs.
>  I was wondering how one can get the time profile of the irreversible work
> from the gromacs pull-code out put . From constant pulling-rate SMD, we get
> time profile of force pullf.xvg and pullx.xvg. I wonder does multiplying
> the the value from pullx.xvg and value from pullf.xvg will provide the work
> . Or, will it be force ( obtained from pullf.xvg) multiplied by pulling
> rate multiplied by time ?
> Sanku
>
> --
> gmx-users mailing list    gmx-users at gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20111104/4905faaf/attachment.html>