

Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
The canonical and grand canonical distribution functions in statistical mechanics. The canonical distribution describes a system in equilibrium with a heat bath, while the grand canonical distribution describes a system in equilibrium with a heat bath and able to exchange particles. Both distributions are derived and their forms are presented. The document also includes an example of the canonical distribution for a monatomic ideal gas and suggests further reading.
Typology: Study notes
1 / 3
This page cannot be seen from the preview
Don't miss anything!


The Canonical Distribution Function L
The Canonical Distribution
ρ(q, p) = ρ 1 (q 1 , p 1 ) ρ 2 (q 2 , p 2 )
of probability distributions for the two sets of variables (q 1 , p 1 ) and (q 2 , p 2 ).
ln ρ(q, p) = ln ρ 1 (q 1 , p 1 ) + ln ρ 2 (q 2 , p 2 ) ,
which means that ln(p, q) must be an additive constant of the motion (recall that ρ is a constant of the motion by Liouville’s theorem). A known result is that the only possibility is that it be a linear combination of the ones related to spacetime and internal symmetries, (E, ~p, L~, Q, ...). Restricting our attention to energy, we get
ln ρ(q, p) = α − β H(q, p) , for some constants α and β.
The constant α is fixed by renormalization, so we get
ρ(q, p) =
Zc
e−βH(q,p)^ , with Zc :=
Γ
dq dp e−βH(q,p)^.
ρ(qb, pb, qs, ps) = constant × δ(Hb + Hs − E) , Eb Es.
Then, using the number Ω(E) of states of energy E accessible to the system,
ρ(qs, ps) =
dqb dpb ρ(qb, pb, qs, ps)
= const × Ωb(E − Hs(qs, ps)) = ∫ Ωb(E^ −^ Hs(qs, ps)) dqbdpb Ωb(E − Hs(qs, ps))
where we have again used implicitly the principle of equal a priori probabilities. Now, since Es E, we can think of Hs as a small perturbation in Ωb. However, for most systems Ωb varies extremely rapidly with its argument, while its log is a function for which a truncated power series is likely to give a better approximation, so
ρ(qs, ps) =
Ωb(E − Hs(qs, ps))/Ωb(E) ∫ dqbdqb Ωb(E − Hs(qs, ps)/Ωb(E)
= N eln[Ωb(E−Hs(qs,ps))/Ωb(E)]^ =
e−βH(qs,ps)^ ,
where
β :=
ln
Ωb(E − Hs(qs, ps)) Ωb(E)
or, if S(E − Hs) := kB ln[Ωb(E − Hs(qs, ps))/Ωb(E)], then β := (1/kB) (∂S(E + x)/∂x)x=0.