Re: [gmx-users] (dH/dl) calculation

2006-11-09 Thread David Mobley
Lei, I am also confused about the details of dgdl (or dHdl) calculation in GROMACS. For example, if the lambda is fixed at certain value, how did the program determine dHdl? See below. My initial understanding was kind of similiar to that of Mauricio that dHdl equals the difference in total

Re: [gmx-users] (dH/dl) calculation

2006-11-09 Thread Lei Zhou
Dear David and Mauricio,   I am also confused about the details of dgdl (or dHdl) calculation in GROMACS. For example, if the lambda is fixed at certain value, how did the program determine dHdl? My initial understanding was kind of similiar to that of Mauricio that dHdl equals the difference in to

Re: [gmx-users] (dH/dl) calculation

2006-11-06 Thread David Mobley
Mauricio, I'm somewhat confused by your question and notation. However, I think the basic answer is something like this: In molecular dynamics, you know the Hamiltonian from which you are sampling; call it H(x,p, l), where x denotes all of the positions, p the momentums, and l lambda. This, of co

[gmx-users] (dH/dl) calculation

2006-11-05 Thread Mauricio Sica
Dear experts I am doing FEP (thermodynamic integration method) simulations. I have a questions about calculation in GROMACS. Take in mind equation 3.77 from the GROMACS 3.3 manual. There, dA/dl is calculated as dA/dl = SS{ (dH/dl) exp()dp dq } / SS{ exp()dp dq = NVT;l } where SS are doble integ