0% found this document useful (0 votes)
79 views5 pages

V.B The Cluster Expansion

1) The document discusses methods for organizing perturbation series expansions for short-range interacting systems, including cluster expansions in powers of density and virial expansions of the equation of state. 2) Diagrammatic representations are introduced to classify terms according to linked clusters, with the contribution of each graph given by the product of its cluster contributions. 3) The grand partition function is reformulated exactly as the exponential of the sum over only connected linked clusters, providing a simple geometrical interpretation of the relationship between moments and cumulants.

Uploaded by

Requale Morter
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
79 views5 pages

V.B The Cluster Expansion

1) The document discusses methods for organizing perturbation series expansions for short-range interacting systems, including cluster expansions in powers of density and virial expansions of the equation of state. 2) Diagrammatic representations are introduced to classify terms according to linked clusters, with the contribution of each graph given by the product of its cluster contributions. 3) The grand partition function is reformulated exactly as the exponential of the sum over only connected linked clusters, providing a simple geometrical interpretation of the relationship between moments and cumulants.

Uploaded by

Requale Morter
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

V.

B The Cluster Expansion


For short range interactions, specially with a hard core, it is much better to replace
the expansion parameter V(~q ) by f (~q ) = exp (V(~q ))1, which is obtained by summing
over all possible number of bonds between two points on a cumulant graph. The resulting
series is organized in powers of the density N/V , and is most suitable for obtaining a virial
expansion, which expresses the deviations from the ideal gas equation of state in a power
series "  2 #
P N N N
= 1 + B2 (T ) + B3 (T ) + . (V.14)
kB T V V V

The temperature dependent parameters, Bi (T ), are known as the virial coefficients and
originate from the inter-particle interactions. Our initial goal is to compute these coeffi-
cients from first principles.
To illustrate a different method of expansion, we shall perform computations in the
grand canonical ensemble. With a macro-state M (T, , V ), the grand partition function
is given by
N
1 e
X X 
N
Q(, T, V ) = e Z(N, T, V ) = SN , (V.15)
N ! 3
N=0 N=0

where
N
Z Y Y
SN = d3 ~qi (1 + fij ), (V.16)
i=1 i<j

and fij = f (~qi ~qj ).


The 2N(N1)/2 terms in S N can now be ordered in powers of fij as

Z Y N X X
SN = d3 ~qi 1 + fij + fij fkl + . (V.17)
i=1 i<j i<j,k<l

An efficient method for organizing the perturbation series is to represent the various con-
tributions diagrammatically. In particular we shall apply the following conventions:
(a) Draw N dots labelled by i = 1, , N to represent the coordinates ~q1 through ~qN ,

.
1 2 N
(b) Each term in eq.(V.17) corresponds to a product of fij , represented by drawing lines
connecting i and j for each fij . For example, the graph,

,
1 2 3 4 5 6 N
102
represents the integral
Z  Z  Z  Z 
3 3 3 3 3 3 3
d ~q1 d ~q2 d ~q3 f23 d ~q4 d ~q5 d ~q6 f45 f56 d ~qN .

As the above example indicates, the value of each graph is the product of the contri-
butions from its linked clusters. Since these clusters are more fundamental, we reformulate
the sum in terms of them by defining a quantity b , equal to the sum over all -particle
linked clusters (one-particle irreducible or not). For example
Z
b1 = = d3 ~q = V, (V.18)

and Z
b2 = = d3 ~q1 d3 ~q2 f (~q1 ~q2 ). (V.19)

There are four diagrams contributing to b3 , leading to


Z
b3 = d3 ~q1 d3 ~q2 d3 ~q3 f (~q1 ~q2 )f (~q2 ~q3 ) + f (~q2 ~q3 )f (~q3 ~q1 ) + f (~q3 ~q1 )f (~q1 ~q2 )



+ f (~q1 ~q2 )f (~q2 ~q3 )f (~q3 ~q1 ) .
(V.20)
A given N -particle graph can be decomposed to n1 1-clusters, n2 2-clusters, , n -
clusters, etc. Hence,
X Y
SN = bn W ({n }), (V.21)
{n }

where the restricted sum is over all distinct divisions of N points into a set of clusters {n },
P
such that n = N . The coefficients W ({n }) are the number of ways of assigning N
particle labels to groups of n -clusters. For example, the divisions of 3 particles into a
1-cluster and a 2-cluster are


, , and .
1 2 3 2 1 3 3 2 1

All above graphs have n1 = 1 and n2 = 1, and contribute a factor of b1 b2 to S3 ; thus


W (1, 1) = 3.
In general, W ({n }) is the number of distinct ways of grouping the labels 1, . . . , N
into bins of n -clusters. It can be obtained from the total number of permutations, N !,
after dividing by the number of equivalent assignments. Within each bin of n particles,
equivalent assignments are obtained by: (i) permuting the labels in each subgroup in !

103
ways, for a total of (!)n permutations; and (ii) the n ! rearrangements of the n subgroups.
Hence,
N!
W ({n }) = Q n
. (V.22)
n !(!)

(We can indeed check that W (1, 1) = 3!/(1!)(2!) = 3 as obtained above.)


Using the above value of W , the expression for S N in eq.(V.21) can be evaluated.
P
However, the restriction of the sum to configurations such that n = N complicates
the evaluation. Fortunately, this restriction disappears in the expression for the grand
partition function in eq.(V.16),

N X
1 e

X N! Y
Q= bn . (V.23)
N ! 3 n
Q
n !(!)
N=0 {n }

P P
The restriction in the second sum is now removed by noting that N=0 {n } P n ,N =
P

{n } . Therefore,

P
X  e  n Y bn X Y 1  e b n
Q= =
3 n !(!)n n ! 3 !
{n } {n }
#n
Y X 1  e  b Y  #
" "
e b
= 3
= exp 3
(V.24)
n ! ! !
{n }

X e  b
"  #
= exp .
3 !
=1

The above result has the simple geometrical interpretation that the sum over all graphs,
connected or not, equals the exponential of the sum over connected graphs. This is a
quite general result that is also related to the graphical connection between moments and
cumulants discussed in sec.II.B.
The grand potential is now obtained from


e

PV X b
ln Q = G = = . (V.25)
kT 3 l!
=1

In eq.(V.25), the extensivity condition is used to get G = E T S N = P V . Thus


the terms on the right hand side of the above equation must also be proportional to the
volume V . This can be explicitly verified by noting that in evaluating each b there is an

104
integral over the center of mass coordinate that explores the whole volume. For example,
b2 = d3 ~q1 d3 ~q2 f (~q1 ~q2 ) = V d3 ~q12 f (~q12 ). Quite generally, we can set
R R

lim b = V b , (V.26)
V

and the pressure is now obtained from


e

P X b
= . (V.27)
kT 3 !
=1

The linked cluster theorem ensures G V , since if any non-linked cluster had appeared in
ln Q, it would have contributed a higher power of V .
Although an expansion for the gas pressure, eq.(V.27) is quite different from eq.(V.14)
in that it involves powers of e rather than the density n = N/V . This difference can be
removed by solving for the density in terms of the chemical potential, using


e

ln Q X V b
N= = . (V.28)
() 3 !
=1

The equation of state can be obtained by eliminating the fugacity x = e /3 , between


the equations

X x P X x
n= b , and = b , (V.29)
( 1)! kT !
=1 =1

using the following steps:


d3 ~q/V = 1)
R
(a) Solve for x(n) from (b1 =

b3 3
x = n b2 x2 x . (V.30)
2

The perturbative solution at each order is obtained by substituting the solution at the
previous order in eq.(V.30),

x1 = n + O(n2 )
x2 = n b2 n2 + O(n3 ) (V.31)
b3 b3
x3 = n b2 (n b2 n)2 n3 + O(n4 ) = n b2 n2 + (2b22 )n3 + O(n4 ).
2 2
105
(b) Substitute the perturbative result for x(n) into eq.(V.29), yielding

b2 2 b3 3
P = x + x + x +
2 6
b3 b2 b3
= n b2 n2 + (2b22 )n3 + n2 b22 n3 + n3 + (V.32)
2 2 6
b2 2 b
3
= n n + (b22 )n3 + O(n4 ).
2 3
The final result is in the form of the virial expansion of eq.(V.14),

X
P = n + B (T )n .
=2

The first term in the series reproduces the ideal gas result. The next two corrections are
b2 1
Z  
B2 = = d3 ~q e V (~q ) 1 , (V.33)
2 2
and
b3
B3 = b22
3
Z 2
V (~

3 q)
= d ~q e 1
 Z  (V.34)
1
Z
3 3 3 3
3 d ~q12 d ~q13 f (~q12 )f (~q13 ) + d ~q12 d ~q13 f (~q12 )f (~q13 )f (~q12 ~q13 )
3
1
Z
= d3 ~q12 d3 ~q13 f (~q12 )f (~q13 )f (~q12 ~q13 ).
3
The above example demonstrates the cancellation of the one particle reducible cluster
that appears in b3 . While all clusters (reducible or not) appear in the sum for b , as
demonstrated in the previous section, only the one particle irreducible ones can appear in
an expansion in powers of density. The final expression for the th virial coefficient is

( 1)
B (T ) = d , (V.35)
!
where d is defined as the sum over all oneparticleirreducible clusters of points. Note
that in terms of d , the partition function can be organized as

X n
ln Z = ln Z0 + V d , (V.36)
!
=2

reproducing the above virial expansion from P = ln Z/V .

106

You might also like