From chemistry-request@ccl.net Tue Jun 17 12:34:13 2003
Received: from erouter0.it-datacntr.louisville.edu (erouter0.it-datacntr.louisville.edu [136.165.5.194])
	by server.ccl.net (8.12.8/8.12.8) with ESMTP id h5HGYD5u006481
	for <chemistry=at=ccl.net>; Tue, 17 Jun 2003 12:34:13 -0400
Received: from savgw-out.louisville.edu (savgw-out.louisville.edu [136.165.2.38])
	by erouter0.it-datacntr.louisville.edu (Postfix) with SMTP id 57FE24D812
	for <chemistry=at=ccl.net>; Tue, 17 Jun 2003 12:34:13 -0400 (EDT)
Received: from gwise.louisville.edu ([136.165.1.106])
 by savgw-out.louisville.edu (SAVSMTP 3.1.1.32) with SMTP id M2003061712330331703
 for <chemistry=at=ccl.net>; Tue, 17 Jun 2003 12:33:03 -0400
Received: from gwgate-MTA by gwise.louisville.edu
	with Novell_GroupWise; Tue, 17 Jun 2003 12:34:03 -0400
Message-Id: <seef0abb.026=at=gwise.louisville.edu>
X-Mailer: Novell GroupWise Internet Agent 6.0.2 Beta
Date: Tue, 17 Jun 2003 12:33:39 -0400
From: "Daquan Gao" <d0gao001=at=gwise.louisville.edu>
To: <chemistry=at=ccl.net>
Cc: <amber=at=heimdal.compchem.ucsf.edu>
Subject: CCL: Amber free energy equation
Mime-Version: 1.0
Content-Type: text/plain; charset=US-ASCII
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

Dear CCLers or AMBER users:

 Having a question about the equation with free energy perpurbation:

The equation is , among other variations:

  Delta G = -kT ln <exp(-(H(lambda + d lambda) - H (lambda))/kT) >

Now in gibbs, I used nstlim = 1, nrun = 100, almdel =0.01, isldyn = -3.

The reason I used nstlim =1 is I wanted to test the equation.

The isldym = -3 means I have a starting state at lambda = 1.0, and a
final state of lambda =0.0 which is when the ligand disappears.

Now my questions:

After the first nrun, there is still only one set of energies, (either
Etot or EPtot or something other terms.). There is already a free energy
for "reverse" given. How was this achieved? There has not been any
energy differences yey.

At nrun = 2, I have lambda = 0.99, the previous state is lambda = 1.00.
Now I have two sets of energies, I would be able to get a delta H =
H(0.99)- H(1.00). This is the ansemble  average, since there is only one
such delta H! But the result is not that number. I don't understand the
number given in gibbs in that stage.

What exactly is the equation that is being used in gibbs? 

I have looked at the references, all says the above equation. Only
Amber7 mentioned a reference by Brooks in JCP. It looked a little
different, but only added DS and DE terms. Don't you get DG direcly from
the above eqution?

If any one knows the answer, please help. It may help me to improve my
results of binding free energy calculations.

I have posted once today, it did not appear, so again, sorry if you get
this twice.
TIA and have a good day. 
Daquan Gao