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 msank...@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 listgmx-users@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-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
--
gmx-users mailing listgmx-users@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-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
