Physical meaning of #tex2html_wrap_inline770# next up previous
Next: Example of a dissociation Up: Free energy profiles Previous: Free energy profiles

Physical meaning of tex2html_wrap_inline758

Consider the free energy difference between two values of the reaction coordinate tex2html_wrap_inline772 and tex2html_wrap_inline774 , which can be written as

displaymath115

The integrand can be written as

displaymath124

Now, we introduce a coordinate transformation from Cartesian coordinates to a set of generalized coordinates that contains q:

displaymath131

where n=3N-1. In addition, a corresponding transformation is made on the conjugate momenta according to:

displaymath133

so that, in the measure, no overall Jacobian appears:

displaymath137

Thus, we have

displaymath139

Next, the derivative with respect to tex2html_wrap_inline780 is changed to a derivative with respect to q:

displaymath146

Then, an integration by parts is performed, which yields

eqnarray153

where the final average is an ensemble average conditional upon the restriction that the reaction coordinate q is equal to tex2html_wrap_inline780 . Generally,

displaymath168

Thus, the free energy difference becomes:

displaymath173

Note that the average in the above expression can also be performed with respect to Cartesian positions and momenta, in which case the derivative can be carried out via the chain rule. If

displaymath180

then

eqnarray184

The quantity tex2html_wrap_inline788 is the generalized force on the generalized coordinate q. Thus, let the conditional average of this force be tex2html_wrap_inline792 . Then,

displaymath203

and

displaymath210

Thus, the free energy difference can be seen to be equal to the work performed against the averaged generalized force in bringing the system from tex2html_wrap_inline774 to tex2html_wrap_inline772 , irrespective of the values of the other degrees of freedom. Such an integration is called thermodynamic integration.

Another useful expression for the free energy derivative can be obtained by integrating out the momenta before performing a transformation. We begin with

displaymath124

Now, noting that the tex2html_wrap_inline746 -function condition is independent of momenta, we can integrate out the Cartesian momenta to yield:

displaymath224

Next, the coordinate transformation to generalized coordinates is performed:

displaymath131

associated with which there is a Jacobian given by

displaymath233

Introducing the coordinate transformation, we obtain

eqnarray236

or

displaymath264


next up previous
Next: Example of a dissociation Up: Free energy profiles Previous: Free energy profiles

Mark Tuckerman
Thu Mar 4 11:03:16 EST 1999