Loading [MathJax]/jax/output/HTML-CSS/jax.js

Saturday, January 24, 2015

The conformationally averaged free energy


The free energy averaged over Nconf conformations of molecule X is

G=RTln(NconfIeGI/RT)  (1)

However, the way to compute an average property of the conformations x is

x=NconfIpIxI  (2)

where

pI=eGI/RTNconfJeGJ/RT (3)

But this equation cannot be used to compute the conformationally averaged free energy because

GG

Why not? First, let's see where the first equation comes from.

Obtaining equation (1)

The first equation comes from the statistical mechanical definition of the Helmholtz free energy found any any p-chem textbook

A=RTln(statesieεi/kT)=RTln(q)

The sum over microstates can be split into sums over microstates of each conformation I

A=RTln(NconfIstatesiIeεi/kT)

=RTln(NconfIqI)=RTln(NconfIeAI/RT)

To get the corresponding expression for the Gibbs free energy [Eq (1)]:

G=A+pV

=RTln(NconfIeAI/RT)RTln(epV/RT)

=RTln(NconfIeGI/RT)

The difference between equation (1) and (3) is the conformational entropy
GG=RTln(NconfJeGJ/RT)NconfIpI

+RTNconfIpIGIRT

=RTNconfIpI(GIRTln(NconfJeGJ/RT))

=RTNconfIpI(ln(eGI/RT)ln(NconfJeGJ/RT))

=RTNconfIpIln(pI)=TSconf


No comments: