General formulation next up previous
Next: Derivation of the Van Up: Distribution functions and perturbation Previous: Distribution functions and perturbation

General formulation

Recall the expression for the configurational partition function:

displaymath38

Suppose that the potential U can be written as a sum of two contributions

displaymath41

where tex2html_wrap_inline553 is, in some sense, small compared to tex2html_wrap_inline555 . An extra bonus can be had if the partition function for tex2html_wrap_inline555 can be evaluated analytically.

Let

displaymath43

Then, we may express tex2html_wrap_inline559 as

eqnarray46

where tex2html_wrap_inline561 means average with respect to tex2html_wrap_inline555 only. If tex2html_wrap_inline553 is small, then the average can be expanded in powers of tex2html_wrap_inline553 :

eqnarray52

The free energy is given by

displaymath60

Separating A into two contributions, we have

displaymath68

where tex2html_wrap_inline571 is independent of tex2html_wrap_inline553 and is given by

displaymath70

and

eqnarray74

We wish to develop an expansion for tex2html_wrap_inline575 of the general form

displaymath82

where tex2html_wrap_inline577 are a set of expansion coefficients that are determined by the condition that such an expansion be consistent with tex2html_wrap_inline579 .

Using the fact that

displaymath89

we have that

eqnarray95

Equating this expansion to the proposed expansion for tex2html_wrap_inline575 , we obtain

displaymath113

This must be solved for each of the undetermined parameters tex2html_wrap_inline577 , which can be done by equating like powers of tex2html_wrap_inline585 on both sides of the equation. Thus, from the tex2html_wrap_inline587 term, we find, from the right side:

displaymath125

and from the left side, the l=1 and k=1 term contributes:

displaymath129

from which it can be easily seen that

displaymath133

Likewise, from the tex2html_wrap_inline593 term,

displaymath135

and from the left side, we see that the l=1,k=2 and l=2,k=1 terms contribute:

displaymath139

Thus,

displaymath143

For tex2html_wrap_inline599 , the right sides gives:

displaymath145

the left side contributes the l=1,k=3, k=2,l=2 and l=3,k=1 terms:

Thus,

displaymath156

Now, the free energy, up to the third order term is given by

eqnarray158

In order to evaluate tex2html_wrap_inline607 , suppose that tex2html_wrap_inline553 is given by a pair potential

displaymath166

Then,

eqnarray170

The free energy is therefore given by

displaymath182


next up previous
Next: Derivation of the Van Up: Distribution functions and perturbation Previous: Distribution functions and perturbation

Mark Tuckerman
Wed Mar 3 18:15:06 EST 1999