Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

6.4: Statistical Mechanics of Independent Identical Particles

( \newcommand{\kernel}{\mathrm{null}\,}\)

6.4.1 Partition function

Now that we have an energy eigenbasis, the obvious thing to do is to calculate the canonical partition function

Z(β)= states eβE,

where for fermions and bosons, respectively, the term “state” implies the occupation number lists:

fermions (n1, n2, · · · , nM), nr = 0, 1, subject to r nr = N
bosons (n1, n2, · · · , nM), subject to r nr = N

As we have seen, it is difficult to even count these lists, much less enumerate them and perform the relevant sum! It can be done, but there is a trick that renders it unnecessary. (Don’t be ashamed if you don’t see the trick. . . neither did Einstein or Fermi. They both did it the hard, canonical way.)

The trick here, as in so many places in statistical mechanics, is to use the grand canonical ensemble. In this ensemble, the partition function is

Ξ(β,μ)= states eβE+βμN= states eβr(nrϵrμnτ)= states Mr=1eβnr(ϵrμ)

where the term “state” now implies the occupation number lists without any restriction on total particle number:

fermions (n1, n2, · · · , nM), nr = 0, 1
bosons (n1, n2, · · · , nM)

Writing out the sum over states explicitly, we have for fermions

Ξ(β,μ)=1n1=01n2=01nM=0Mr=1eβnr(ϵrμ)=[1n1=0eβn1(ϵ1μ)][1n2=0eβn2(ϵ1μ)][1n=0eβnM(ϵ1μ)].

A typical factor in the product is

[1nr=0eβnr(ϵrμ)]=1+eβ(ϵrμ),

so for fermions

Ξ(β,μ)=Mr=11+eβ(ϵrμ).

Meanwhile, for bosons, the explicit state sum is

Ξ(β,μ)=n1=0n2=0nM=0Mr=1eβnr(ϵrμ)=[n1=0eβn1(ϵ1μ)][n2=0eβn2(ϵ1μ)][nM=0eβnM(ϵ1μ)].

and a typical factor in the product is

[nr=0eβnr(ϵrμ)]=1+[eβ(ϵrμ)]+[eβ(ϵrμ)]2+[eβ(ϵrμ)]3+=11eβ(ϵrμ),

where in the last line we have summed the geometric series under the assumption that ϵr > µ. Thus for bosons

Ξ(β,μ)=Mr=111eβ(ϵrμ).

The two results are compactly written together as

Ξ(β,μ)=Mr=1[1±eβ(ϵrμ)]±1,

where the + sign refers to fermions and the − sign to bosons.

6.4.2 Mean occupation numbers

In our previous work, we have always found the partition function and worked from there. Surprisingly, however, for the situation of quantal ideal gases it is more useful to find the mean occupation numbers, such as

n5= states n5eβ(EμN) states eβ(EμN)= states n5eβr(nrϵrnrμ)Ξ(β,μ).

Note that the averages nr are functions of β and µ (as well as of r) but it is notationally clumsy to show that dependence.

How can such averages be evaluated? A slick trick would be helpful here! Consider the derivative

lnΞϵ5=1ΞΞϵ5=1Ξ states (βn5)eβr(nrϵrnrμ)=βn5.

Using the explicit expression (6.37) for Ξ (where the + sign refers to fermions and the − sign to bosons), this gives

n5=1βlnΞϵ5

=1β{ϵ5Mr=1ln[1±eβ(ϵrμ)]±1}

=1β{ϵ5ln[1±eβ(ϵSμ)]±1}

=1β{±±(β)eβ(ϵ5μ)[1±eβ(e5μ)]}

=1eβ(ϵSμ)±1,

leaving us with the final result.

nr=1eβ(ϵrμ)±1.

As before, the + sign refers to fermions and the − sign to bosons.

The mean occupation numbers play such an important role that it is easy to forget that they are only averages, that there will be fluctuations, that for a given T and µ not all states will have exactly n5 building blocks of level 5 (see problem 6.13). Keep this in mind if you ever find yourself saying “occupation number” rather than “mean occupation number”.

In practice, these results from the grand canonical ensemble are used as follows: One uses these results to find quantities of interest as functions of temperature, volume, and chemical potential, such as the pressure p(T, V, µ). But most experiments are done with a fixed number of particles N, so at the very end of your calculation you will want to find µ(T, V, N) in order to express your final answer as p(T, V, N). You can find µ(T, V, N) by demanding that

N=Mr=1nr=Mr=11eβ(erμ)±1.

In other words, the quantity µ serves as a parameter to insure normalization, very much as the quantity Z serves to insure normalization in the canonical ensemble through

1=neβEnZ.

You might wonder, in fact, about the relation between the canonical probability

eβEnZ,

which we have seen many times before, and the recently derived occupancy probability

1N1eβ(ϵrμ)±1.

The first result applies to both interacting and non-interacting systems, both classical and quantal. The second applies only to non-interacting quantal systems. Why do we need a new probability? What was wrong with our derivation (in section 4.1) of the canonical probability that requires us to replace it with an occupancy probability? The answer is that nothing was wrong and that the occupancy probability doesn’t replace the canonical probability. The canonical probability and the occupancy probability answer different questions. The first finds the probability that the entire system is in the many-body state n. The second finds the probability that the one-body level r is used as a building block in constructing the many-body state. Indeed, although we derived the occupancy probability result through a grand canonical argument, it is also possible to derive the occupancy probabilities from strict canonically arguments, proof that these two probabilities can coexist peacefully.

6.4.3 The Boltzmann limit

This is the limit where particles are far enough apart that overlap of wavefunction is minimal, so we needn’t worry about symmetrization or antisymmetrization. Equivalently, it is the limit where nr for all r.

6.4.4 Problems

6.10 Evaluation of the grand canonical partition function

Can you find a simple expression for Ξ(β, µ) for non-interacting particles in a one-dimensional harmonic well? For non-interacting particles in a one-dimensional infinite square well? For any other potential? Can you do anything valuable with such an expression once you’ve found it?

6.11 Entropy of quantal ideal gases

This problem derives an expression for the entropy of a quantal ideal gas in terms of the mean occupation numbers hnri. (Compare problem 4.3.) Throughout the problem, in the symbols ± and ∓, the top sign refers to fermions and the bottom sign refers to bosons.

a. Use the connection between thermodynamics and statistical mechanics to show that, for any system,

S(T,V,μ)kB=lnΞβlnΞβ.

b. Show that for the quantal ideal gas,

lnΞ(T,V,μ)=rln(1nr).

c. The mean occupation numbers nr are functions of T, V, and µ (although it is notationally clumsy to show this dependence). Show that

βnrβ)V,μ=β(ϵrμ)nr(1nr)=[ln(1nr)lnnr|nr(1nr).

d. Finally, show that

S(T,V,μ)=kBr[nrlnnr±(1nr)ln(1nr)].

e. Find a good approximation for this expression in the Boltzmann limit, nr1.

f. (Optional.) Find an expression for CV in terms of the quantities nr.

6.12 Isothermal compressibility of quantal ideal gases

a. Show that in a quantal ideal gas, the isothermal compressibility is

κT=1ρkBT[1r(nr)2rnr],

where as usual the top sign refers to fermions and the bottom sign to bosons. (Clue: Choose the most appropriate expression for κT from those uncovered in problem 3.33.)

b. Compare this expression to that for a classical (“Maxwell-Boltzmann”) ideal gas.

c. The negative sign in the expression for fermions opens the possibility that κT could be negative. Prove that this potential horror never happens. d. Do the relative sizes of the three compressibilites (fermion, classical, boson) adhere to your qualitative expectations? (Compare problem 6.29.)

6.13 Dispersion in occupation number

Find an expression analogous to (6.45) giving the dispersion in the occupation numbers. (Clue: A slick trick would be helpful here.) Answer:

Δnr=1eβ(ϵrμ)/2±eβ(ϵrμ)/2=nr(1nr)


This page titled 6.4: Statistical Mechanics of Independent Identical Particles is shared under a CC BY-SA license and was authored, remixed, and/or curated by Daniel F. Styer.

Support Center

How can we help?