**From:** Rathinavelan, Thenmalarchelvi (*thenmcr_at_ku.edu*)

**Date:** Sat Sep 01 2007 - 17:53:55 CDT

**Next message:**Rima Chaudhuri: "calculate interaction energies for a system with amber forcefield"**Previous message:**Christopher Harrison: "Re: About QM-MM"**Next in thread:**JC Gumbart: "Re: Total external work calculation in SMD-PMF calculation"**Reply:**JC Gumbart: "Re: Total external work calculation in SMD-PMF calculation"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

Dear users,

I would greatly appreciate if anyone help me to understand the total work calculation of the SMD-PMF reconstruction that is given in Jensen et al., PNAS,2002 & Amaro et al., PNAS, 2003. It has been given in the paper that in each time slice the average total external work over all the SMD generated trajectories is calculated by,

<W(j)avg>traj = 1/N Sum N{integral over time slice(j) dt/deltaT * W(t)-

0.5k(<z(j)avg> -z0 –vt(j)avg)**2}

(Sorry for the bad representation of the mathematical symbols)

‘Sum N’ – summation over N trajectories

‘< >traj’ - average over trajectory

‘(j)avg’ - average over the time slice

Here is my question:

I like to know how the total work W(t) is calculated. According to my understanding from the papers, it is calculated for each time slice, for instance, in the case of the time slice between 200-400ps (deltaT=200ps), the total work is the summation of work carried out over this time slice (200-400ps) at a given interval (for example, 2ps). Can anyone please tell me is this correct?

Thanks,

Thenmalar

**Next message:**Rima Chaudhuri: "calculate interaction energies for a system with amber forcefield"**Previous message:**Christopher Harrison: "Re: About QM-MM"**Next in thread:**JC Gumbart: "Re: Total external work calculation in SMD-PMF calculation"**Reply:**JC Gumbart: "Re: Total external work calculation in SMD-PMF calculation"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

*
This archive was generated by hypermail 2.1.6
: Wed Feb 29 2012 - 15:45:11 CST
*