Processing math: 100%
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

3.2: Calculating Chemical Potentials

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

Now let us discuss properties of ideal gases of free, indistinguishable particles in more detail, paying special attention to the chemical potential μ – which, for some readers, may still be a somewhat mysterious aspect of the Fermi and Bose distributions. Note again that particle indistinguishability requires the absence of thermal excitations of their internal degrees of freedom, so that in the balance of this chapter such excitations will be ignored, and the particle’s energy εk will be associated with its “external” energy alone: for a free particle in an ideal gas, with its kinetic energy (3.1.3).

Let us start from the classical gas, and recall the conclusion of thermodynamics that μ is just the Gibbs potential per unit particle – see Equation (1.5.7). Hence we can calculate μ=G/N from Eqs. (1.4.26) and (3.1.17). The result,

μ=TlnVN+f(T)+Tln[NgV(2π2mT)3/2],

which may be rewritten as

exp{μT}=NgV(2π2mT)3/2,

gives us some information about μ not only for a classical gas but for quantum (Fermi and Bose) gases as well. Indeed, we already know that for indistinguishable particles, the Boltzmann distribution (2.8.1) is valid only if Nk<<1. Comparing this condition with the quantum statistics (2.8.5) and (2.8.8), we see again that the condition of the gas behaving classically may be expressed as

exp{μεkT}<<1

for all εk. Since the lowest value of εk given by Equation (3.1.3) is zero, Equation (???) may be satisfied only if exp{μ/T}<<1. This means that the chemical potential of a classical gas has to be not just negative, but also “strongly negative” in the sense

μ>>T.

According to Equation (???-???), this important condition may be represented as

T>>T0,

with T0 defined as

Quantum scale of temperature:

T02m(NgV)2/32m(ng)2/32g2/3mr2ave ,

where rave is the average distance between the gas particles:

rave1n1/3=(VN)1/3.

In this form, the condition (???-???) is very transparent physically: disregarding the factor g2/3 (which is typically of the order of 1), it means that the average thermal energy of a particle, which is always of the order of T, has to be much larger than the energy of quantization of particle’s motion at the length rave. An alternative form of the same condition is19

rave>>g1/3rc, where rc(mT)1/2.

For a typical gas (say, N2, with m14mp2.3×1026 kg) at the standard room temperature (T=kB×300 K 4.1×1021 J), the correlation length rc is close to 1011 m, i.e. is significantly smaller than the physical size a3×1010 m of the molecule. This estimate shows that at room temperature, as soon as any practical gas is rare enough to be ideal (rave>>a), it is classical, i.e. the only way to observe quantum effects in the translational motion of molecules is very deep refrigeration. According to Equation (???), for the same nitrogen molecule, taking rave102a108 m (to ensure that direct interaction effects are negligible), the temperature should be well below 1 mK.

In order to analyze quantitatively what happens with gases when T is reduced to such low values, we need to calculate μ for an arbitrary ideal gas of indistinguishable particles. Let us use the lucky fact that the Fermi-Dirac and the Bose-Einstein statistics may be represented with one formula:

N(ε)=1e(εμ)/T±1,

where (and everywhere in the balance of this section) the top sign stands for fermions and the lower one for bosons, to discuss fermionic and bosonic ideal gases in one shot.

If we deal with a member of the grand canonical ensemble (Figure 2.7.1), in which not only T but also μ is externally fixed, we may use Equation (???) to calculate the average number N of particles in volume V. If the volume is so large that N>>1, we may use the general state counting rule (3.1.13) to get

N=gV(2π)3N(ε)d3k=gV(2πh)3d3pe[ε(p)μ]/T±1=gV(2πh)304πp2dpe[ε(p)μ]/T±1.

In most practical cases, however, the number N of gas particles is fixed by particle confinement (i.e. the gas portion under study is a member of a canonical ensemble – see Figure 2.4.1), and hence μ rather than N should be calculated. Let us use the trick already mentioned in Sec. 2.8: if N is very large, the relative fluctuation of the particle number, at fixed μ, is negligibly small (δN/N1/N<<1), and the relation between the average values of N and μ should not depend on which of these variables is exactly fixed.

Hence, Equation (???), with μ having the sense of the average chemical potential, should be valid even if N is exactly fixed, so that the small fluctuations of N are replaced with (equally small) fluctuations of μ. Physically, in this case the role of the μ-fixing environment for any sub-portion of the gas is played by the rest of it, and Equation (???) expresses the condition of self-consistency of such chemical equilibrium.

So, at N>>1, Equation (???) may be used for calculating the average μ as a function of two independent parameters: N (i.e. the gas density n=N/V) and temperature T. For carrying out this calculation, it is convenient to convert the right-hand side of Equation (???) to an integral over the particle’s energy ε(p)=p2/2m, so that p=(2mε)1/2, and dp=(m/2ε)1/2dε, getting

Basic equation for μ:

N=gVm3/22π230ε1/2dεe(εμ)/T±1.

This key result may be represented in two other, more convenient forms. First, Equation (???), derived for our current (3D, isotropic and parabolic-dispersion) approximation (3.1.3), is just a particular case of the following self-evident state-counting relation

N=0g(ε)N(ε)dε,

where

g(ε)dNstates/dε

is the temperature-independent density of all quantum states of a particle – regardless of whether they are occupied or not. Indeed, according to the general Equation (3.1.4), for our simple model (3.1.3),

g(ε)=g3(ε)dNstatesdε=ddε(gV(2π)34π3p3)=gVm3/22π23ε1/2,

so that we return to Equation (???).

On the other hand, for some calculations, it is convenient to introduce the following dimensionless energy variable: ξε/T, to express Equation (???) via a dimensionless integral:

N=gV(mT)3/22π230ξ1/2dξeξμ/T±1.

As a sanity check, in the classical limit (???-???), the exponent in the denominator of the fraction under the integral is much larger than 1, and Equation (???) reduces to

N=gV(mT)3/22π230ξ1/2dξeξμ/TgV(mT)3/22π23exp{μT}0ξ1/2eξdξ, at μT.

By the definition of the gamma function Γ(ξ),20 the last integral is just Γ(3/2)=π1/2/2, and we get

exp{μT}=N2π23gV(mT)3/22π=(2πT0T)3/2,

which is exactly the same result as given by Equation (???-???), obtained earlier in a rather different way – from the Boltzmann distribution and thermodynamic identities.

Unfortunately, in the general case of arbitrary μ, the integral in Equation (???) cannot be worked out analytically.21 The best we can do is to use T0, defined by Equation (???), to rewrite Equation (???) in the following convenient, fully dimensionless form:

TT0=[12π20ξ1/2dξeξμ/T±1]2/3,

and then use this relation to calculate the ratios T/T0 and μ/T0(μ/T)×(T/T0), as functions of μ/T numerically. After that, we may plot the results versus each other, now considering the first ratio as the argument. Figure 3.2.1 below shows the resulting plots, for both particle types. They show that at high temperatures, T>>T0, the chemical potential is negative and approaches the classical behavior given by Equation (???) for both fermions and bosons – just as we could expect. However, at temperatures TT0 the type of statistics becomes crucial. For fermions, the reduction of temperature leads to μ changing its sign from negative to positive, and then approaching a constant positive value called the Fermi energy, εF7.595 T0 at T0. On the contrary, the chemical potential of a bosonic gas stays negative, and then turns into zero at a certain critical temperature Tc3.313 T0. Both these limits, which are very important for applications, may (and will be :-) explored analytically, separately for each statistics.

fig-ch01_patchfile_01.jpg
Figure 3.2.1: The chemical potential of an ideal gas of N>>1 indistinguishable quantum particles, as a function of temperature at a fixed gas density nN/V (i.e. fixed T0n2/3), for two different particle types. The dashed line shows the classical approximation (???), valid only at T>>T0.

Before carrying out such studies (in the next two sections), let me show that, rather surprisingly, for any non-relativistic, ideal quantum gas, the relation between the product PV and the energy,

Ideal gas: PV vs. E

PV=23E,

is exactly the same as follows from Eqs. (3.1.19) and (3.1.20) for the classical gas, and hence does not depend on the particle statistics. To prove this, it is sufficient to use Eqs. (2.8.4) and (2.8.7) for the grand thermodynamic potential of each quantum state, which may be conveniently represented by a single formula,

Ωk=Tln(1±e(μεk)/T),

and sum them over all states k, using the general summation formula (3.1.13). The result for the total grand potential of a 3D gas with the dispersion law (3.1.3) is

Ω=TgV(2π)30ln(1±e(μp2/2m)/T)4πp2dp=TgVm3/22π230ln(1±e(με)/T)ε1/2dε.

Working out this integral by parts, exactly as we did it with the one in Equation (2.6.10), we get

Ω=23gVm3/22π230ε3/2dεe(εμ)/T±1=230εg3(ε)N(ε)dε.

But the last integral is just the total energy E of the gas:

Ideal gas: energy

E=gV(2π)30p22m4πp2dpe[ε(p)μ]/T±1=gVm3/22π230ε3/2dεe(εμ)/T±1=0εg3(ε)N(ε)dε,

so that for any temperature and any particle type, Ω=(2/3)E. But since, from thermodynamics, Ω=PV, we have Equation (???) proved. This universal relation22 will be repeatedly used below.


This page titled 3.2: Calculating Chemical Potentials is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Konstantin K. Likharev via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?