Thermo Dynamics Script
Thermo Dynamics Script
Ph.W. Courteille
Universidade de São Paulo
Instituto de Fı́sica de São Carlos
12/06/2024
2
.
3
.
4
Preface
The script was developed for the course Physical Chemistry and Thermodynamics
o Solids (SFI5769) offered by the Institute of Physics of São Carlos (IFSC) of the
University of São Paulo (USP).
The course is intended for masters and PhD students in physics. The script
is a preliminary version continually being subject to corrections and modifications.
Error notifications and suggestions for improvement are always welcome. The script
incorporates exercises the solutions of which can be obtained from the author.
I Thermodynamics 1
1 Foundations and mathematical formalism 3
1.1 Phenomenological thermodynamics at equilibrium . . . . . . . . . . . 3
1.1.1 Temperature . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.1.2 Kinetic theory and microscopic interpretation of temperature . 7
1.1.3 Heat and work . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.1.4 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.2 Canonical formulation of thermodynamics . . . . . . . . . . . . . . . . 20
1.2.1 Tackling thermodynamic systems . . . . . . . . . . . . . . . . . 21
1.2.2 The laws of thermodynamics . . . . . . . . . . . . . . . . . . . 24
1.2.3 Thermodynamic potentials . . . . . . . . . . . . . . . . . . . . 28
1.2.4 Strategy for deriving thermodynamic relations . . . . . . . . . 33
1.2.5 Ideal gases . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
1.2.6 Cyclic processes . . . . . . . . . . . . . . . . . . . . . . . . . . 41
1.2.7 Real gases, liquids and solids . . . . . . . . . . . . . . . . . . . 44
1.2.8 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
1.3 Thermodynamic equilibrium . . . . . . . . . . . . . . . . . . . . . . . . 58
1.3.1 Conditions for equilibrium . . . . . . . . . . . . . . . . . . . . . 59
1.3.2 Entropy maximization in two-phase systems, chemical potential 60
1.3.3 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
1.4 Thermodynamic ensembles . . . . . . . . . . . . . . . . . . . . . . . . 62
1.4.1 Coupling of thermodynamic ensembles to reservoirs . . . . . . 63
1.4.2 Thermodynamic potentials associated to specific ensembles . . 64
1.4.3 Microcanonical ensemble . . . . . . . . . . . . . . . . . . . . . . 65
1.4.4 Canonical ensemble . . . . . . . . . . . . . . . . . . . . . . . . 69
1.4.5 Grand canonical ensemble . . . . . . . . . . . . . . . . . . . . . 71
1.4.6 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
1.5 Further reading . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
5
6 CONTENT
2.3 Multi-component, heterogeneous, non-reacting systems . . . . . . . . . 96
2.3.1 Equilibrium conditions . . . . . . . . . . . . . . . . . . . . . . . 96
2.3.2 Structure of phase diagrams . . . . . . . . . . . . . . . . . . . . 98
2.3.3 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
2.4 Continuous, non-uniform systems exposed to external forces . . . . . . 99
2.4.1 Thermodynamic densities . . . . . . . . . . . . . . . . . . . . . 99
2.4.2 Equilibrium conditions . . . . . . . . . . . . . . . . . . . . . . . 100
2.4.3 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
2.5 Reacting systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
2.5.1 Univariant chemical reactions in the gas phase . . . . . . . . . 106
2.5.2 Multi-variant chemical reactions in the gas phase . . . . . . . . 108
2.5.3 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
2.6 Classification of thermodynamic phase transitions . . . . . . . . . . . . 111
2.6.1 Solid-liquid-vapor . . . . . . . . . . . . . . . . . . . . . . . . . . 111
2.6.2 Bose-Einstein condensation . . . . . . . . . . . . . . . . . . . . 112
2.7 Materials . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
2.7.1 Electrons in solids . . . . . . . . . . . . . . . . . . . . . . . . . 112
2.7.2 Plasmas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
2.7.3 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
2.8 Further reading . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
Thermodynamics
1
Chapter 1
Foundations and
mathematical formalism
Thermodynamics is a central branch of modern science, and its general laws govern
the physical and chemical processes which occur in our world. An important early
application of thermodynamics dealt with steam (or heat) engines, in which heat is
converted to mechanical energy. Phenomenological thermodynamics was developed in
the nineteenth and in the beginning of the twentieth century by Watt, Carnot, Clau-
sius, Joule, von Helmholtz, Lord Kelvin, Nernst, Boltzmann, and Gibbs, culminating
in the discovery of the Laws of Thermodynamics. These laws set general limits for
the conversion of one form of energy, for example heat or chemical energy, to another
one, for example mechanical work.
The generic question addressed by thermodynamics is, how a given system re-
sponds to environmental changes. Indeed, the sudden modification of an environment
will force the system to seek a new state of equilibrium, as illustrated in Fig. 1.1. On
the other hand, thermodynamics is limited to describing equilibrium. It does not tell
us step by step, how the new equilibrium is reached, only what the final state will be.
Nevertheless, this is sufficient to establish complete phase diagrams, which are maps
of equilibrium states. Fig. 1.2 shows as an example the phase diagram of water.
3
4 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
Figure 1.3: (a) Cylinder with a piston containing a gas. (b) Interaction between two systems
through a wall.
1.1.1 Temperature
Temperature is, in general, measured by observing some quantity which is sensitive
to its variation. Defining a temperature scale on this quantity in a particular sys-
tem, we construct a thermometer. An example is the mercury thermometer, where a
certain volume of liquid mercury is placed in a capillary glass tube and the thermal
expansion of mercury is observed as a function of temperature. The length L of the
mercury column (usually calibrated in degrees Celsius) varies approximately linearly
6 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
with temperature T , that is, T = aL + b, where a and b are two constants that depend
on chosen reference temperatures. Conventionally, the melting temperatures of ice
(triple point) at 0 ◦ C and the boiling temperature of water at 100 ◦ C are used, so that,
100 100Lg
a= and b= . (1.1)
Lv − Lg Lv − Lg
The gas thermometer illustrated in Fig. 1.4 is another possible realization, in which
the volume of a gas is used as thermometric quantity.
P V = N kB T . (1.2)
where T is the absolute temperature (in K), N is the number of molecules contained
in the volume V , and kB is the Boltzmann constant. This is the most famous equation
of state, and it is based on experimental observations of Boyle, Mariotte, and Gay-
Lussac. Eq. (1.3) can also be written in terms of the number of moles, which is a
quantity defined by n = N/NA , where NA is Avogadro’s number. In this case,
P V = nRg T , (1.3)
1.1. PHENOMENOLOGICAL THERMODYNAMICS AT EQUILIBRIUM 7
Since half of the particles move to the left, the number of collisions in the time interval
∆t is given by,
1 Ni
#= Avxi ∆t . (1.5)
2 V
The total change of momentum is,
2
Ni mvxi A∆t
∆I = ∆pxi # = , (1.6)
V
and the pressure is,
2
∆I Ni mvxi
Pi =
= . (1.7)
A∆t V
Hence, as the movement is isotropic,
X mX 2 m N 2
P = Pi = Ni vxi ≡ vx2 = mv 2 = N Ēkin . (1.8)
i
V i V 3V 3
Figure 1.5: Box with N molecules colliding with the wall of a container.
Using the ideal gas state equation (1.3), Eq. (1.8) yields,
3 mv 2
kB T = , (1.9)
2 2
from which we conclude that temperature is associated with translational energy of
ideal gas molecules. This expression is also known as the theorem of energy equiparti-
tion theorem. Generally speaking, we assign to each degree of freedom of a system the
term 12 kB T . In the example above, we have 3 degrees of freedom, which correspond
to translations in the x, y and z directions. If energy can be stored in vibrations or
rotations of a molecule, we also have to assign a term 12 kB T to each of these degrees
of freedom.
Example 1 (Johnson noise): As an example demonstrating the usefulness of
the energy equipartition theorem, let us consider a resistor R subject to a certain
temperature T . If we associate 21 kB T with the average power P̄J dissipated by
the resistor within a time ∆t we have,
r
ŪJ2 1 kB T R
P̄J ∆t = ∆t = kB T =⇒ ŪJ = , (1.10)
R 2 2∆t
that is, a small voltage ŪJ appears at the resistor terminals, which is known as
Johnson noise. For a 1 Ω resistor at ambient temperature we find ŪJ ≈ 4.0 nV
averaged over 1 second. This voltage is small, but must be taken into account
in high precision measurements.
The surface expansion of a body is then, with A = L1 L2 = L10 L20 + 2αL10 L20 ∆T +
L10 L20 α2 ∆T 2 and assuming α to be very small,
∆A ≃ A0 2α∆T . (1.13)
∆V ≃ V0 3α∆T . (1.14)
Figure 1.6: (a) Interaction energy between two atoms. (b) temperature-dependence of the
density of water.
Liquids and gases also undergo volume variations with temperature. In this case,
it is quite common to work with the fluid density instead of volume:
m m m
ρ= =⇒ ∆ρ = − ∆V = − 2 γV ∆T , (1.15)
V V2 V
where γ is the volumetric expansion coefficient. Therefore, ∆ρ = −γρ0 ∆. In general,
γ is positive and the density of the fluid decreases with temperature. An exception to
10 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
this rule is the case of water, as shown in Fig. 1.6(b), which below 4 ◦ C it has γ < 0
and thus, between 4 ◦ C and 0 ◦ C the density increases with the temperature. This
explains why during winter lakes freeze starting from the surface. Do the Excs. 1.1.4.9
to 1.1.4.12.
(i) Type of bar material. – There are materials that conduct heat better than
others, e.g. copper conducts heat better than steel. A measure is provided by
the thermal conductivity K.
(ii) The cross section A. – The larger it is, the greater the thermal current, as more
atoms are participating in the conduction process.
(iii) The temperature gradient. – The thermal current depends on the difference in
temperature between adjacent layers of atoms (left and right of the plane at
position x).
Based on the above considerations, we can write the following expression for the
thermal current H:
δQ dT
H= = −KA . (1.21)
dt dx
If the bar is thermally insulated, as is the case in Fig. 1.7, the current is conserved,
that is, all the heat entering one end of the bar will come out at the other, as there are
no losses. In this situation, H is independent of x and, consequently, dT dx is constant.
Hence,
dT T2 − T1
= , (1.22)
dx L
and consequently,
T2 − T1
H = KA . (1.23)
L
In this case, the temperature distribution is a straight line, as shown in Fig. 1.7(b).
On the other hand, if the side surface of the bar is not insulated, there will be heat
12 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
losses by convection and the thermal current decreases as x increases. In this case,
dT
dx also decreases and, as a consequence, we have a temperature distribution like the
one illustrated by the dotted line in Fig. 1.7(b).
Figure 1.7: (a) Heat conduction through a laterally insulated bar. (b) Temperature distri-
bution along an insulated bar (a) and not isolated (b) laterally.
We will now analyze two examples where the equation of conductivity (1.23) ap-
plies.
T1 − T T − T2
H = K1 A = K2 A , (1.24)
L1 L2
yielding,
K2 L1 T2 + K1 L2 T1
T = . (1.25)
K2 L1 + K1 L2
Substituting this in the expression for H we get,
K1 K2 A
H= (T1 − T2 ) . (1.26)
K2 L1 + K1 L2
In the particular case in which K1 = K2 we recover the result (1.23) derived for
a single isolated bar. The temperature distribution along the bars depends on
the ratio between K1 and K2 . If K1 > K2 , we have the temperature distribution
shown in Fig. 1.8(b).
dT dT
H = −KA = −K2πrL . (1.27)
dr dr
1.1. PHENOMENOLOGICAL THERMODYNAMICS AT EQUILIBRIUM 13
Figure 1.8: (a) Bars of different materials placed in series. (b) Temperature distribution
along two bars of different materials placed in series.
As H is constant, since there are no losses, we can integrate this equality between
a ≤ r ≤ b assuming T1 ≥ T ≥ T2 . We find the result,
2πLK
H= (T1 − T2 ) . (1.28)
ln(b/a)
Let us now briefly cover the other two types of heat transmission mentioned in the
beginning. In the case of convection, we are typically interested in the following type
of problem: given a body at a temperature T surrounded by atmospheric air colder
by an amount ∆T , how much heat does it lose per unit of time? The thermal current
from the body to the air is described by a similar formula as thermal conduction,
H = hA∆T , (1.29)
where A is the area through which heat is being lost and h is a number that depends
on ∆T (in general h ∝ ∆T 1/4 ), the geometry of the body and its orientation in space
(since convection is due to the fact that hot air rises). Therefore, a lying plate has a
different h as compared to a standing plate.
Heat transport by radiation is proportional to T 4 , where T is the absolute tem-
perature (in Kelvin):
R = eσT 4 . (1.30)
Here, R is the thermal current emitted per unit area, e is the emissivity of the body
(0 ≤ e ≤ 1), and σ is the Stefan-Boltzmann constant. Do the Exc. 1.1.4.15 to 1.1.4.18.
1.1.3.2 Work
Heat can be injected into a system in various ways, e.g. via thermal contact with a
reservoir, via electrical dissipation or discharge, or by carrying out mechanical work.
The amount of energy (translational, rotational and vibrational) contained in a system
is called internal energy. The law of energy conservation demands that the heat
δQ supplied to a system is used to change its internal energy E and/or to perform
mechanical work δW . This principle is known as the 1st law of thermodynamics and
can be expressed mathematically as:
dE = δQ + δW . (1.31)
14 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
Figure 1.9: (a) P V diagram of a gas showing work carried out in two different processes
leading from an initial state (1) to a final state (2). (b) Work performed in a complete gas
cycle.
E = 32 N kB T , (1.35)
1.1. PHENOMENOLOGICAL THERMODYNAMICS AT EQUILIBRIUM 15
that is, for processes maintaining the temperature constant, ∆E = 0. This relation-
ship was experimentally verified by Joule, who found that when a gas adiabatically
expands (without heat exchange nor work executed) its temperature stays constant.
After the warm up provided by the previous sections, let us now start setting up
the framework of phenomenological thermodynamics in the following sections.
1.1.4 Exercises
1.1.4.1 Ex: Gas thermometer
A gas thermometer filled with an ideal gas and working at a constant volume is cali-
brated on the one hand in dry ice (carbon dioxide in its solid state at a temperature
of −80.0 ◦ C) and on the other hand in boiling alcohol (78.0 ◦ C). At these respective
temperatures, the pressure in the gas thermometer is 0.900 bar or 1.635 bar. At ab-
solute zero, the gas in the thermometer is still gaseous, but the pressure has dropped
to 0.000 bar.
a. At what ◦ C is the absolute zero?
b. What is the pressure at the freezing point of water and what is it at the boiling
point?
g dm gρ dV
dP = = = −gρ dh . (1.36)
A A
Here dV = −A dh is the change in volume of the air column (located above the base).
The air should be treated approximately as a substance with a uniform molar mass
M̃ .
a. Show with the help of the ideal gas equation that Eq. (1.36) can be cast into the
form dP = −kP dh under the given conditions with the constant k. Which is the
expression for k?
b. What is the integral relationship P = P (h)? At what height h is the air pressure
at T = 273 K and M̃ = 29 g/mol only half the size of P (0)?
c. What is the air pressure on the Mont Blanc (4794 m) and the Mount Everest
(8848 m) at T = 273 K as compared to P (0) at sea level? How big is the pressure
difference ∆P compared to normal zero on the Tübingen market place (h = 341 m)?
b. Calculate the volume of the bubbles as a function of water depth. How big is the
volume just below the water surface? Why is it important for the diver to exhale
continuously as he ascends?
Numerical values: Water depth h0 = 40 m, V0 = 1 cm3 , T1 = 20 ◦ C, α = 0.2 ◦ C/m,
ρ = 1 kg/l, and P1 = 1013 hPa.
A closed box with end face A and side length L is divided into two equal halves by a
movable plate (see figure). Both halves contain one mole of helium under a pressure
of P0 . The movable plate is now shifted to the right by the distance x. The shift
takes place at constant temperature T = 20 ◦ C.
a. Give the volume of the right or left sub-box as a function of x. Give the pressure
in the right or left sub-box as a function of x.
b. Calculate the work W that needs to be done to move the plate from x = 0 to
x = L/4. Specify W in Joules.
18 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
264 Termologia e termodinâmica
x
10- Encontre o gradiente de temperatura e a corrente térmica numa barra de
condutividade K, comprimento L e secção transversal irregular, como
mostra a Fig. 13.16.
A
L
L/2 L/2
L
T1>T2 FigureA1.12: Ideal2A T2
gas in a box.
L/2
1.1.4.9 Ex: Bi-metal Fig. 13.16
Two bars12-of different
Duas barras materials, withdiferentes,
de materiais lengths, Young’s modules and
com comprimentos, thermal
módulos de expansion
coefficients given
Youngrespectively
e coeficientesby de L 1 , L2 , Y
dilatação 1 , Y2 , dados
térmica α1 , and α2 , are pinched
respectivamente por L1between
, two
walls, as shown
L2, Yin
1 , Ythe
2 , αfigure.
1 e α2 , Calculate
estão presas the
entre distance
duas traveled
paredes como by
mostrathe
a junction
Fig. point
of the bars when
13.14.the system
Calcule is heated
a distância by an amount
percorrida pelo ponto ∆Tde. junção
What das
is the tension on the
barras
bars?
quando o sistema é aquecido de ∆T. Qual é a tensão nas barras?
α1 , Y1 α2 , Y2
L1 L2
Fig. 13.17
Figure 1.13:
T: d
L
T+∆T :
20
(
CHAPTER
A
)
F = Y ∆L está presa entre duas paredes, conforme mostra a Fig. 13.15.
L 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
Calcule a tensão na barra quando a temperatura é acrescida de ∆T.
Fig. 13.15
Figure 1.15:
7-
264 Qual a quantidade de calor necessária para transformar
Termologia1g de gelo a –10
e termodinâmica
0
C (cgelo = 0.55 cal/g C, Lf = 80 cal/g) em vapor a 100 0C (LV = 540
0
1.1.4.17 Ex: Thermal conduction
10- Encontre
cal/g)? o gradiente de temperatura e a corrente térmica numa barra de
Find the temperature
Coloca-segradient
8- condutividade umaK,barraand thermal
comprimento
de metal current
(CL =e0,2 in
C) aa bar
cal/gotransversal
secção ofirregular,
100oC conductivity
sobre grande K, length
um como
L and irregular cross
mostra
bloco section,
adeFig.
gelo o
13.16. as Qual
a 0 C. shown é a in thedafigure.
massa barra se quando o sistema atingir
o equilíbrio térmico 500 g de gelo se derreteram?
L
9- Coloca-se um bloco de gelo a –20 0C dentro de um recipiente
hermeticamente fechado com 200g de vapor de água a 100 0C. Se a massa
T1>T2 A 2A T2
do gelo é 500 g, qual será a temperatura final do sistema?
L/2
Fig. 13.16
Figure 1.16:
12- Duas barras de materiais diferentes, com comprimentos, módulos de
S. C. Zilio e V. S. Bagnato Mecânica, calor e ondas
Young e coeficientes de dilatação térmica dados respectivamente por L1,
L2, Y1, Y2, α1 e α2, estão presas entre duas paredes como mostra a Fig.
1.1.4.18 Ex: Tramway
13.14. Calcule a distância percorrida pelo ponto de junção das barras
A tram with quando
mass m B = 12500
o sistema kg brakes
é aquecido de ∆T. from
Qual éaa tensão
speednas
v= 57.6 km/h to standstill.
barras?
What is the temperature of the eight cast iron brake blocks when the mass of each
block is 9.0 kg and 60% of the kinetic energy flows into the heating of the blocks?
α1 , Y1 α2 , Y2
In the second step, exposed in the subsequent chapter 4, is called statistical thermo-
dynamics. It aims at explaining the laws found in phenomenological thermodynamics
by deriving them from features of the microscopic subsystems. In particular, some
subsystems are small (or cold) enough to behave according to the rules of quantum
mechanics, and this can have an important impact on the macroscopic behavior. This
is studied in the field of quantum statistics exposed in Sec. 4.2.
The methods introduced in the first two chapters will be illustrated with simple
examples (mainly ideal gases and solids). In Chp. 2 we apply them to complex real
systems, e.g. multi-component, heterogeneous, or chemically interacting systems.
system. The most common ones are temperature, pressure, entropy, volume, particle
number, chemical potential, and internal energy,
T, S, P, V, N, µ . (1.37)
The pressure and volume may be replaced by other mechanical variables. Their
equilibrium dynamics is confined to trajectories obeying so-called state equations,
that is, functional relationships such as,
f (T, S, P, ...) = 0 , (1.38)
or their illustration in phase diagrams, such as the one shown in Fig. 1.2. Such
state equations are obviously extremely useful, as they allow to abstract from the
physical process which led to a particular state of a system, and which might be very
complicated. In fact, their discovery represents one of the major achievements of
thermodynamics.
In contrast to state functions, the values of process variables depend upon the
path followed by the process. Consequently, they only have a meaning for changing
systems, i.e. for systems traversing a sequence of different states. The process variables
fall into two category, mechanical work and exchanged heat,
δQ, δW . (1.39)
The concept of work is developed in classical mechanics
R as the path integral over a
force acting on a body along a given path, W = S F(s) · ds, and it can be brought
into the context of thermodynamics considering, e.g. a piston working against the
pressure of a confined gas, as illustrated in Fig. 1.3. Work depends on changes ds
and thus cannot be associated with stationary systems. Any of the forces known in
physics, inertial forces in accelerated or rotating systems, electromagnetic forces, or
molecular forces, can work.
Work is associated with a displacement of macroscopic matter, e.g. the movement
of a piston. Systems may, however, exchange energy without net displacement of
masses via the exchange of heat. By itself the notion of a ’quantity of heat’ in a
system is meaningless. Only the amount of heat exchanged with another system in a
given process can be quantified.
An example for a coefficient relation emanating from the real gas equation T =
T (P, V ) = N 1kB (P + P ∗ )(V − V ∗ ) is,
V −V∗ P + P∗
∂T ∂T
dT = AdP + BdV = dP + dV = dP + dV .
∂P V ∂V P N kB N kB
T1 = T2 ∧ T2 = T3 =⇒ T1 = T3 . (1.40)
1. The first law states that the total energy is always conserved,
dE = δQ + δW . (1.41)
2. The second law states that the entropy of any closed system goes always in-
creasing,
dS ≥ 0 . (1.42)
3. The third law states that for T → 0, the entropy difference between systems
connected by a reversible process vanishes,
lim S = 0 . (1.43)
T →0
Figure 1.17: Any change of a system’s internal energy is due to either work done on it or
heat transferred through its borders.
change the overall entropy balances, the total entropy production of the universe,
consisting of the system plus its environment, must remain positive 2 .
For reversible processes the process variables are readily calculated, since each
intermediate state (red dashed line in Fig. 1.18) is described by just a few state
functions, e.g. temperature or pressure. In contrast, this is very complicated for
irreversible processes carrying the system to a state whose state functions depend on
the trajectory on which the state was reached.
2 We may formulate a continuity-type equation for entropy density and entropy flow.
1.2. CANONICAL FORMULATION OF THERMODYNAMICS 27
On the other hand, changes of state functions are easy to calculate also for ir-
reversible processes, since these only depend on the initial and final states. Thus,
changes in state functions for a given complex irreversible process may be calculated
by substituting the process by a reversible one. E.g. in Fig. 1.18 instead of following
solid red lines, we substitute them by dashed red lines connecting the same initial and
final states. This procedure illustrates the fundamental role of state functions and
reversible processes in thermodynamics.
Let us now analyze entropy transfer for reversible processes. Let δQrev be the
heat absorbed in an infinitesimal step and T the temperature of the system. Now,
although δQrev is a process variable, δQrev /T is the differential of a state function,
that is, the state function for an infinitesimal process. To prove this, it is sufficient
to show that for a general cyclic process, such as the Carnot cycle, the path integral
I
δQrev
=0 (1.45)
T
vanishes. The state function δQrev /T is defined to be the entropy,
δQrev = T dS . (1.46)
An analogous treatment can also be done for work. We express reversible work
by a combination of state functions,
δWrev = −P dV . (1.50)
Combining the first and second law of thermodynamics Eqs. (1.41), (1.46), and (1.50),
we find for reversible processes,
dE = T dS − P dV . (1.51)
28 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
Q, W (1.53)
α, κ, CP , CV . (1.55)
These material properties may also depend on temperature or pressure. For a given
substance the state equations can be computed from those properties.
The general procedure to tackle a problem is the following:
1.2. CANONICAL FORMULATION OF THERMODYNAMICS 29
1. Identify the properties of the system about which information is available. These
are the independent variables, for example, T and P .
3. Such functions will necessarily contain material properties, which will have to
be looked up from data bases.
The crucial step in this procedure is evidently the second one, which consists in finding
the appropriate state function.
E, H, F, G . (1.56)
Which one of the defined potentials is used as a state function in a particular problem
is, in principle, arbitrary. However, some processes are easier to describe in terms
of particular potentials. Before we define them below, let us present a useful mathe-
matical framework facilitating the conversion between state functions called Legendre
transform.
The partial derivatives A, B, C are evaluated assuming that the derived quantities do
not depend on other variables. Since the order of the derivatives can be inverted, we
know,
! !
∂A ∂ ∂W ∂ ∂W ∂B
= = = (1.58)
∂Y X,Z ∂Y ∂X Y,Z ∂X ∂Y X,Z ∂X Y,Z
X,Z Y,Z
! !
∂A ∂ ∂W ∂ ∂W ∂C
= = = ,
∂Z X,Y ∂Z ∂X Y,Z ∂X ∂Z X,Y ∂X Y,Z
X,Y Y,Z
and analogously for all other second derivatives. Applied to thermodynamic systems
these expression are called Maxwell relations.
30 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
Now, we want to predict how the system will evolve when a different set of varia-
tions is applied, for instance (dA, dY, dZ). To solve the problem we first define a new
state function V = V (A, Y, Z) via,
V ≡ W − AX . (1.59)
and each of these three expressions has a physical meaning. The thermodynamic
potentials introduced in the next section will exemplify the procedure. Do the
Exc. 1.2.8.1.
Useful mathematical identities when working with partial derivatives are,
∂X ∂Z ∂Z ∂X ∂Y
=1 , = −1 . (1.61)
∂Z Y ∂X Y ∂X Y ∂Y Z ∂Z X
1.2.3.2 Enthalpy
The energy function called enthalpy is defined as,
H ≡ E + PV . (1.62)
dH = dE + P dV + V dP = T dS + V dP . (1.63)
It has the same level of generality as the combined first and second laws.
Historically, the enthalpy was introduced to simplify the description of heat engines
taken through cycles at atmospheric pressure, dP = 0. For isobaric processes in
simple systems, dHP = T dSP = δQrev,P , the enthalpy provides a direct measure of
the reversible heat exchange of the engine with the environment.
F ≡ E − TS . (1.64)
It has the same level of generality as the combined first and second laws.
This function was devised to simplify the description of processes occurring at
a fixed (if necessary stabilized) temperature, dT = 0. For isothermal processes in
simple systems, dFT = −P dVT = δWrev,T , the Helmholtz free energy reports the
total reversible work done on the system.
1.2. CANONICAL FORMULATION OF THERMODYNAMICS 31
G ≡ E + PV − TS . (1.66)
It has the same level of generality as the combined first and second laws.
This function was introduced to simplify the description of processes occurring
at both temperature, dT = 0, and constant pressure, dP = 0. For isobaric and
isothermal processes in simple systems, dGT,P = 0. But in systems undergoing phase
transformations or chemical reactions, the Gibbs free energy yields the total work
′
other than mechanical work, dGT,P = δWT,P .
δQrev = T dS
δWrev = −P dV
E = Qrev + Wrev =⇒ dE = T dS − P dV
(1.68)
H = E + PV =⇒ dH = T dS + V dP
F = E − TS =⇒ dF = −SdT − P dV
G = E + PV − TS =⇒ dG = −SdT + V dP
Thermal properties are grasped by the concept of heat capacity. This quantity
is measured via the temperature rise of a substance due to reversible absorption
of a defined quantity of heat. Since heat is a process variable, the heat capacity
measurement will depend on the circumstances, i.e. whether the pressure is kept
constant during the measurement or the volume. The heat capacity will be for the
respective cases,
δQrev δQrev
CP ≡ or CV ≡ . (1.72)
dT P dT V
V = V (T, P ) =⇒ dV = αV dT − κV dP
CP
S = S(T, P ) =⇒ dS = T dT − αV dP
E = E(T, P ) =⇒ dE = (CP − αP V )dT + (κP − αT )V dP
(1.80)
H = H(T, P ) =⇒ dH = CP dT + (1 − αT )V dP
F = F (T, P ) =⇒ dF = −(S + αP V )dT + κP V dP
G = G(T, P ) =⇒ dG = −SdT + V dP
P V = N kB T , (1.88)
where we consider for now a fixed number of particles N = const. The thermal expan-
sion coefficient and the compressibility are readily calculated from their definitions
(1.69),
1 ∂V 1 ∂N kB T /P 1
αid = = = (1.89)
V ∂T P V ∂T P T
1 ∂V 1 ∂N kB T /P 1
κid = − =− =
V ∂P T V ∂P T P
1 ∂P 1 ∂N kB T /V 1
βid = = = .
P ∂T V P ∂T V T
The relationship between the heat capacities at constant volume and pressure is di-
rectly obtained from (1.86),
α2 T V
∂S
CV = T = CP − id = CP − N kB . (1.90)
∂T V κid
Hence, for an ideal gas, the heat capacities do not depend on T not P , but only on
the number of atoms N and their configuration in each gas molecule,
CV = 32 N kB . (1.91)
36 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
We will see later, that for a molecular gas with f degrees of freedom (accessible at
the temperature), the results must be generalized to,
CV = f2 N kB . (1.92)
Thus, E only depends on T , and since E is a state function, this holds for any process,
whether reversible or irreversible. For the other thermodynamic potentials we get,
exploiting the relationships (1.68),
E = CV T , F = (CV − S)T
(1.94)
H = CP T , G = (CP − S)T
δQadiab = 0 . (1.95)
δQadiab,rev = T dS = 0 , (1.96)
we finally find,
N kB /CV N kB /CP CP /CV
T2 V1 P2 P2 V1
= = or = . (1.100)
T1 V2 P1 P1 V2
1.2. CANONICAL FORMULATION OF THERMODYNAMICS 37
The ratio
CP
γ≡ , (1.101)
CV
is called the adiabaticity coefficient. For the ideal gas studied in the previous section
we have, γ = 53 , and for a molecular gas with f degrees of freedom γ = 1 + f2 .
With the adiabaticity coefficient the state functions for adiabatic reversible pro-
cesses can be written,
where the second relation follows from the ideal gas equations (1.88). Such a process
is shown in the P V -diagram of Fig. 1.21.
Since the specific heat can be obtained from the heat capacity simply by dividing
by the mass, c ≡ C/m, the ratio of the specific heat at constant pressure and to the
specific heat at constant volume is also equal to the constant γ, that is, γ = cP /cV .
Do the Excs. 1.2.8.13 to 1.2.8.22.
Since no entropy can flow out of the system, any change in entropy must come
from local production arising from irreversible processes. The expansion of the gas is
a complicated process occurring far from equilibrium. Nevertheless, since entropy is
a state function, the gain in entropy is same as for a fictive reversible process leading
to the same final state. Since the initial state of the gas is known (V1 , Ti ) as well as
the final volume, Vf = V1 + V2 , the final temperature Tf is determined. Now, since
the internal energy does not change during expansion, dE = δQ + δW = 0, and since
for an ideal gas the internal energy is proportional to temperature, the temperature
doesn’t change,
P V = N kB T = const . (1.103)
38 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
Such a process is called isothermal and is shown in the P V -diagram of Fig. 1.21.
On an isotherm of an ideal gas, we have P dV + V dP = N kB T = 0 and therefore
δW = P dV = −V dP . Nonetheless,
N kB T N kB T dV
P = ⇒ dP = − dV ⇒ δW = N kB T , (1.104)
V V2 V
and so, the work done when the gas goes from V1 to V2 is,
Z V2
dV V2
∆W1−2 = N kB T = N kB T ln . (1.105)
V1 V V1
Similar considerations can be made for isobaric and isochoric processes. The results
are summarized in Tab. 1.2.
V2 P2 P2 V2γ T2 V2γ−1
∆S1→2 = CP ln + CV ln = CV ln = C V ln . (1.109)
V1 P1 P1 V1γ T1 V1γ−1
(a) (b)
4 3 δW = 0
3
2
S
T
dE = 0
1
1 δQ = 0
2 2
0 0 dE = 0
0 1 0 1
1 1
2 0 V 2 0 V
P T
Figure 1.20: (code) (a) Qualitative behavior of the ideal gas law in T P V space. (b) Qual-
itative behavior of the Sackur-Tetrode formula in T SV space with planes corresponding to
δQ = 0, δW = 0, and δQ + δW .
Z S2 Z S2 Z V2
process ∆S1→2 = dS = ∆Q1→2 = T dS = ∆W1→2 = − P dV =
S1 S1 V1
Z Z Z
κ
dST → S2 − S1 T dST → T (S2 − S1 ) − α
P dST → −∆Q1→2
Z Z Z
2 2
isotherm − αV dPT → −N kB ln P P1
− αV dPT → −N kB T ln P P1
κP V dPT → −∆Q1→2
Z Z Z
α α
κ
dVT → N kB ln VV12 κ
dVT → N kB T ln VV12 − P dVT → −∆Q1→2
Z
P
0 0 ( κC
αT
− αV )P dTS → CV (T2 − T1 )
Z
2 P2 V2 −P1 V1
isentropic 0 0 ( αCTPV − κ)P V dPS → γ−1
Z
P2 V2 −P1 V1
0 0 − P dVS → γ−1
Z Z Z
CP T2
T
dTP → CP ln T1
CP dTP → CP (T2 − T1 ) − αP V dTP → −N kB (T2 − T1 )
Z Z Z Z
αT V
isobar dSP → S2 − S1 T dSP → CP T1 (e∆SP /CP − 1) − CP
dSP → − NCkPB T dSP
Z Z Z
CP V2 CP CP P
αT V
dVP → CP ln V1 αV
dVP → N kB
(V2 − V1 ) − P dVP → −P (V2 − V1 )
Z Z
CP α2 V T2 α2 T V
T
− κ
dTV → CV ln T1
(CP − κ
)dTV → CV (T2 − T1 ) 0
Z Z
isochor dSV → S2 − S1 T dSV → CV T1 (e∆SV /CV − 1) 0
Z Z
P2 P CV V
dPV → CV ln P1
( κC
αT
− αV )dPV → N kB
(P2 − P1 ) 0
Table 1.2: Summary of entropy, heat, and work balances upon various state changes.
CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
1.2. CANONICAL FORMULATION OF THERMODYNAMICS 41
We define the efficiency of a cyclic process as the ratio between net work performed
and heat absorbed (not delivered),
P
− j ∆Wj→j+1
η≡ P >0 . (1.112)
j ∆Qj→j+1
Figure 1.21: P V -diagram (left) and T S-diagram (right) for the Clément-Desormes cyclic
process.
change in the volume of the gas is ∆V = yA, so that there is a pressure variation
∆P accompanying ∆V . This increase in pressure produces a restoring force
F = ∆P A. Assuming that the process is nearly static and adiabatic, we have,
P V γ = const ⇒ γP V γ−1 ∆V + V γ ∆P = 0 . (1.115)
Using ∆P = F/A and ∆V = yA, we obtain,
γP A2
γP V γ−1 yA + V γ F/A = 0 ⇒ F =− y, (1.116)
V
1.2. CANONICAL FORMULATION OF THERMODYNAMICS 43
d2 y γP A2
m + y=0, (1.117)
dt2 V
which is the differential equation of a simple harmonic motion of frequency:
γP A2
ω02 = . (1.118)
mV
In this way, knowing P , V , m and A, we can measure τ and obtain γ as,
r
mV
τ = 2π . (1.119)
γP A2
∆Qhot
∆W = (Thot − Tcold )∆S = (Thot − Tcold ) , (1.120)
Thot
where Qhot is heat transferred from the hot reservoir to the system per cycle.
The heat and work balance of the Carnot cycle is summarized in the following
table (see Fig. 1.23):
Figure 1.23: P V -diagram (left) and T S-diagram (right) for the Carnot cycle.
Figure 1.24: T S-diagram for an arbitrary cyclic process (left) and for the Carnot cycle
(right).
The ideal gas law (1.3) is only valid for non-interacting particles. In reality, inter-
particle interactions increase the effective pressure and the finite size of the molecules
reduces the effective volume. Indeed, even at T = 0 the volume V of a real gas cannot
be zero, because the molecules have their own volume V ∗ . And as molecules interact
attract each other, the pressure is zero even before T = 0. In the van der Waals model
the ideal gas equation is generalized to,
2
a N N
(P +P ∗ )(V −V ∗ ) = N kB T with P ∗ = and V ∗ = b (1.123)
V2 NA NA
where a and b are empirical constants specifically depending on the gas. Real gases
are studied in Excs. 1.2.8.34 to 1.2.8.36. Using molar functions denoted by a tilde (.̃),
N → NA , NA kB → Rg , V NNA → Ṽ , etc.,
a
P+ (Ṽ − b) = Rg T . (1.124)
Ṽ 2
1.5
1
1
P/Pc
T /Tc
0
0.5
-1 1.2
1
0 1
-0.5 0.5 -0.5
0 0 0.8 T /Tc
0.5 1 0 P/Pc 0.5 1
lg (V /Vc ) lg (V /Vc )
The Joule-Thomson effect describes the temperature change of a real gas or liquid
(as differentiated from an ideal gas), when it is forced through a valve or porous plug
while keeping it insulated so that no heat is exchanged with the environment (see
Exc. 1.2.8.37).
The relationships (1.80)(ii) and (iii) express energy and entropy as a function of
temperature and pressure, Ẽ = Ẽ(T, P ) and S̃ = S̃(T, P ). Using the procedure
outlined in Sec. 1.2.4 to express energy and entropy as a function of temperature and
46 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
Hence, ! !
∂ Ẽ ∂ S̃
=T −P . (1.126)
∂ Ṽ T
∂ Ṽ T
From the Maxwell relation derived from, dF̃ = −S̃dT − P dṼ , we find,
! ! ! ! !
∂P ∂ ∂ F̃ ∂ ∂ F̃ ∂ S̃
=− =− = . (1.127)
∂T V ∂T ∂ Ṽ ∂V ∂T ∂ Ṽ T
T V V T
Hence, !
∂ Ẽ ∂P
=T −P . (1.128)
∂ Ṽ ∂T V
T
Now, we consider a dense real van der Waals gas according to Eq. (1.124) but
neglecting the volume parameter, b = 0. The pressure P then behaves as a function
of temperature T and molar volume Ṽ according to the following state equation,
Rg T a
P = − , (1.129)
Ṽ Ṽ 2
where a is a positive constant and Rg is the universal gas constant. From the equation
of state we obtain,
∂P Rg
= , (1.130)
∂T V Ṽ
which, replaced in the expression (1.128), gives,
!
∂ Ẽ Rg a
=T −P = . (1.131)
∂ Ṽ T Ṽ Ṽ 2
Inserting the integral of (1.133), K(T ) = C̃V T + K0 , into (1.132), we end up with,
a
Ẽ = − + C̃V T + K0 . (1.134)
Ṽ
In a process of free expansion, Ẽ remains invariant and Ṽ grows. Thus, resolving
(1.134) by temperature and deriving by volume,
!
∂T ∂ Ẽ + a/Ṽ − K0 −a
= = . (1.135)
∂ Ṽ E ∂V CV CV Ṽ 2
E
For an ideal gas, a = 0, we expect no temperature change. For a gas with positive
(negative) a the variation (1.135) will be negative (positive).
Example 11 (Microscopic interpretation of Joule-Thomson cooling ): In
the discussion of the the process illustrated in Fig. 1.19 we stated that an ideal
gas does not changes its temperature when it expands flowing through a nozzle
from one volume into another. This is in contrast to the behavior of a real gas.
Indeed, in a compressed gas the molecules are closer to each other and thus
feel attractive (or repulsive) van der Waals forces. When the gas expands, the
molecules must overcome these forces at the cost (gain) of kinetic energy.
1.2.8 Exercises
1.2.8.1 Ex: Material parameters
The compressibility κ, the thermal expansion coefficient α, and the stress coefficient β
are important material parameters. In the case of a unary system (single substance),
they are defined by,
1 ∂V 1 ∂V 1 ∂P
κ=− , α= , β= .
V ∂P ν,T V ∂T ν,P P ∂T ν,V
Here ν is the number of moles.
a. Show that these relationships can be rewritten, using the molar volume Ṽ , to,
! !
1 ∂ Ṽ 1 ∂ Ṽ 1 ∂P
κ=− , α= , β= .
Ṽ ∂P T Ṽ ∂T P P ∂T Ṽ
a. Explain why the specific heat at constant volume is less than the specific heat at
constant pressure.
CP
b. Show that for a diatomic gas γ = C V
= 75 .
a. Two gas containers are brought into thermal contact. They contain gases with
the temperatures T1 and T2 , as well as the heat capacities C1 and C2 . The thermal
capacity of the containers is negligible. What is the temperature of the gases after an
equilibrium has been reached?
b. Now consider the temperature equilibrium of three containers, each with 100 g of
the gas H2 at the temperature TH2 = 10 ◦ C, 50 g of the gas He at the temperature
THe = 15 ◦ C, and 200 g of the gas N2 at temperature TN2 = 20 ◦ C. What is the final
temperature?
An oxygen bottle with the volume V2 = 40 l contains a filling ex works that would
have the volume V1 = 6 m3 at atmospheric pressure P1 = 101 kPa. The bottle, which
has been emptied to atmospheric pressure, is refilled at a constant temperature of
T1 = 18 ◦ C. What mechanical work W must be added to the gas to compress it
isothermally from P1 to the filling pressure?
1.2. CANONICAL FORMULATION OF THERMODYNAMICS 49
1 kmol of nitrogen under normal conditions (P0 = 1.01 × 105 Pa, T = 0 ◦ C) adiabat-
ically expands from V1 to V2 = 5V1 . Calculate the change in the internal energy of
the gas and the amount of work the gas does as it expands.
a. During the adiabatic expansion of a gas the pressure P and the volume V of the
gas satisfy the relationship P V γ = α, where α is a constant, and γ is the factor
of the gas that gives the ratio between the specific heats at constant pressure and
volume, i.e. γ = cP /cV . A gas was placed in a cylinder with a movable (frictionless)
plunger completely insulated from the external environment. The assembly makes it
possible to measure the volume and pressure of the gas during its expansion and the
experimental values obtained are given in the table below.
a. From the values in the table below, and using the least squares method, deter-
mine the gas factor and the constant α. (Hint: To obtain a linear relationship, take
x = lg V and y = lg P .).
b. Determine, through the method of least squares, the uncertainties in the values
obtained for γ and α.
c. Using a log×log paper, prepare a graph P × V and determine the values of γ and
α. Compare with the results obtained by the least squares method.
Notes: When displaying the values of γ and α, be sure to indicate the units in which
they are expressed.
Display the values of Sx , Sy , Sx2 and Sxy used in the least squares method calcula-
tions.
V (l) P (atm)
40 1.20
41 1.16
43 1.10
44 1.05
46 0.98
47 0.96
49 0.90
50 0.87
In a thermally insulated container B there are n mol of an ideal gas and a body K
with the heat capacity C. Specify the relationship between pressure P and volume
V , whereby the change in V is carried out so slowly that the following always applies
to body and gas: TK = TG .
Note: Body and gas exchange heat. Assume that the container itself does not take
heat from the gas or the body.
50 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
dV
B K
Figure 1.26: Potato.
1.2.8.24 Ex: Heat power engine and heat pump based on Carnot cycle
Consider a heat power engine and a heat pump based on the Carnot cycle. Calculate
the efficiency and the generated, respectively, consumed power as a function of the
temperature difference between the hot and the cold bath.
1.2. CANONICAL FORMULATION OF THERMODYNAMICS 53
3
1 0
2
Figure 1.27: Principle of a heat pump.
Figure 1.28: Otto cyclic process. The black branches are adiabatic processes.
Figure 1.29: P V -diagram (left) and T S-diagram (right) for the Diesel cycle. The black
branches are adiabatic processes.
54 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
Figure 1.30: P V -diagram (left) and T S-diagram (right) for the Brayton cycle.
An engine, whose four work steps consist of two isothermal and two isochoric pro-
cesses, runs at a speed of N = 500 min-1 . There is ν = 0.5 mol of an ideal, mono-
atomic gas in the volume of the engine. The parameters for the individually labeled
working steps are T1 = 50 ◦ C, P1 = 2 bar, and P2 = 5 bar.
a. Determine the volumes V1 = V2 and V3 = V4 , the pressure P4 , and the temperature
T2 .
b. What efficiency η does the cycle have?
c. How high is the net power P of the engine?
A thermodynamic system is brought from an initial state (1) to another state (2), then
to (3), and finally back to (1), as illustrated in the diagram in Fig. 1.34. Calculate
the work and heat balance for the entire cycle.
56 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
conditions,
dP d2 P
= =0?
dṼ dṼ 2
b. Show that the temperature of the critical point Tc and the volume of the system at
the critical point Ṽc are linked to the material constants a and b by the relationships:
8a
Tc = and Ṽc = 3b .
27Rg b
P eα/(Rg T Ṽ ) (Ṽ − β) = Rg T ,
increasing or decreasing their volumes. Finally, it is assumed that the gas is adiabat-
ically isolated from the environment and therefore only exchanges with it energy in
the form of work.
a. Show that the enthalpy H remains constant in both recipients during this process.
b. Show that,
∂T 1 ∂V
= T −V ,
∂P H CP ∂T P
and calculate (∂T /∂P )H explicitly for an ideal gas.
c. For a real gas, the so-called inversion curve P (T ) defined by (∂T /∂P )H is obtained
in the P T -plane. Physically interpret the areas above and below this curve. Calculate
the inversion curve for the van der Waals gas using the thermal equation of state for
real gases.
d. Discuss the behavior of entropy in this process.
Figure 1.37: (a) Stable and unstable equilibrium positions in a mechanical system. (b) A
thermally insulated copper rod conducting heat from a hot source to a cold sink develops a
stationary temperature profile across its length (solid line in c). However, when disconnected
from its surroundings (source and sink) it changes a temperature profile (dashed line in c).
3. Substitute dXk in Eq. (1.138) and collect the coefficients Cl of the remaining
terms dXl|l̸=k ,
dW = C1 dXl1 + C2 dXl2 + ... . (1.140)
transiently with some pump or sink. If it stays stationary once the pump or sink is removed, it should
be considered at equilibrium.
60 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
(a) (b)
1 1
Z(X, Y )
Z(X, Y )
0.5 0.5
0 Y (X) 0 Z2 (X, Y )
1 1 1 1
0.5 0.5 0.5 0.5
0 0 0 0
Y X Y X
Figure 1.38: (code) Illustration of state function and constraints. The surface in (a) visualizes
a state function depending on two variables; the black curve illustrates a possible constraint
Y = Y (x) between the variables. Fig. (b) shows the same surface as in (a), now constrained
by and additional surface Z2 = Z2 (X, Y ).
AE + BG + C = 0 = AF + BH + D , (1.143)
S1 + S2 = Stot , (1.144)
between the phases, as the state functions depend on the particle number, we need
to consider this dependency,
Ej = Ej (Sj , Vj , Nj ) , (1.145)
and so on for the other state functions. This leads to a generalization of the differential
form (1.68),
dEj = Tj dSj − Pj dVj + µj dNj , (1.146)
1 1 P1 P2 µ1 µ2
dStot = dS1 + dS2 = dE1 + dE2 + dV1 + dV2 − dN1 − dN2 . (1.148)
T1 T2 T1 T2 T1 T2
The fact that the system is closed implies, dE1 + dE2 = 0, dV1 + dV2 = 0, and
dN1 + dN2 = 0. Using these constraints in (1.148) we get,
1 1 P1 P2 µ1 µ2
dStot = − dE1 + − dV1 − − dN1 . (1.149)
T1 T2 T1 T2 T1 T2
Hence, the systems comes to equilibrium once the following conditions are fulfilled,
T1 = T2 thermal equilibrium
P1 = P2 mechanical equilibrium . (1.150)
µ1 = µ2 chemical equilibrium
Application of any of these criteria leads to the same set of conditions for equilibrium.
62 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
1.3.3 Exercises
1.3.3.1 Ex: Gibbs free energy
δQrev = T dS
δWrev = −P dV
R
E = (T dS − P dV + µdN ) =⇒ dE = T dS − P dV + µdN
Ψ = E − µN =⇒ dΨ = T dS − P dV − N dµ
H = E + PV =⇒ dH = T dS + V dP + µdN
(1.151)
Φ = H − µN =⇒ dΦ = T dS + V dP − N dµ
F = E − TS =⇒ dF = −SdT − P dV + µdN
Ω = F − µN =⇒ dΩ = −SdT − P dV − N dµ
G = H − TS =⇒ dG = −SdT + V dP + µdN
O = G − µN = 0 =⇒ dO = −SdT + V dP − N dµ
4 Note that [13] claims on page 65 that this relationship holds for the enthalpy H, which is wrong.
1.4. THERMODYNAMIC ENSEMBLES 63
CP σ
dS = T
dT − αV dP + µ
dµ
dV = αV dT − κV dP + ηV dµ . (1.154)
σ
dN = µ
dT − ηV dP + κN dµ
The question is now, what are the additional material constants γ, κ, and η?
Do the Exc. 1.4.6.2.
Figure 1.39: Insulated or heat conducting, closed or open, rigid or working systems.
In the following we will qualitatively discuss the most common ones, which are
the microcanonical, the canonical, and the grand canonical ensembles.
5 For this reason, the microcanonical ensemble is also called N V E-ensemble.
1.4. THERMODYNAMIC ENSEMBLES 65
Pmc = W −1 . (1.155)
1
|ψk ⟩f ( E−ε
P
ρ̂mc = ∆E )⟨ψk |
k
Ξmc k
. (1.156)
E−εk
f (x) = θ( 21 − |x|)
P
Ξmc = k f ( ∆E ) with
1.4.3.1 Applicability
Because of its connection with the elementary assumptions of equilibrium statistical
mechanics (particularly the postulate of a priori equal probabilities), the microcanoni-
cal ensemble is an important conceptual building block in the theory and is sometimes
considered to be the fundamental distribution of equilibrium statistical mechanics. It
66 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
dv
Ssur = kB ln = SB − kB ln ∆E . (1.159)
dE
P ∂S µ ∂S
= , =− . (1.161)
T ∂V T ∂N
1.4. THERMODYNAMIC ENSEMBLES 67
The volume entropy Svol and associated temperature Tvol form a close analogy to
thermodynamic entropy and temperature. It is possible to show exactly that,
dE = Tvol dSvol − ⟨P ⟩dV , (1.162)
where ⟨P ⟩ is the ensemble average pressure, as expected for the first law of thermo-
dynamics. A similar equation can be found for the surface entropy and its associated
temperature Tsur , however the ’pressure’ in this equation is a complicated quantity
unrelated to the average pressure.
The microcanonical Tvol and Tsur are not entirely satisfactory in their analogy to
temperature. Outside of the thermodynamic limit, a number of artefacts occur.
• Nontrivial result of combining two systems: Two systems, each described by
an independent microcanonical ensemble, can be brought into thermal contact
and be allowed to equilibriate into a combined system also described by a mi-
crocanonical ensemble. Unfortunately, the energy flow between the two systems
cannot be predicted based on the initial T ’s. Even when the initial T ’s are
equal, there may be energy transferred. Moreover, the T of the combination is
different from the initial values. This contradicts the intuition that temperature
should be an intensive quantity, and that two equal-temperature systems should
be unaffected by being brought into thermal contact.
• Strange behavior for few-particle systems: Many results, such as the micro-
canonical equipartition theorem acquire a one- or two-degree of freedom offset
when written in terms of Tsur . For a small systems this offset is significant, and
so if we make Ssur the analogue of entropy, several exceptions need to be made
for systems with only one or two degrees of freedom.
• Spurious negative temperatures: A negative Tsur occurs whenever the density
of states decreases with energy. In some systems the density of states is not
monotonic in energy, and so Tsur can change sign multiple times as the energy
is increased. The preferred solution to these problems is to avoid using the
microcanonical ensemble. In many realistic cases a system is thermostatted to
a heat bath so that the energy is not precisely known. Then, a more accurate
description is the canonical ensemble or grand canonical ensemble, both of which
have complete correspondence to thermodynamics.
68 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
1 ∂S 3 N kB
≡ = .
T ∂E 2 E
which agrees with analogous result from the kinetic theory of gases. Calculating
the pressure gives the ideal gas law,
P ∂S N kB
≡ = .
T ∂V V
Finally, the chemical potential µ is,
" 3/2 #
∂S V 4πmE
µ ≡ −T = kB T ln .
∂N N 3N
mgz 5N
5N mg 2
−2
ρ(z) = −1 1− .
2 E E
1.4. THERMODYNAMIC ENSEMBLES 69
Similarly, the distribution of the velocity magnitude |v| (averaged over all heights)
is,
5(N2−1)
m3/2 |v|2 m|v|2
Γ(5N/2)
f (|v|) = × 1/2 3/2 × 1 − .
Γ(3/2)Γ(5N/2 − 3/2) 2 E 2E
The analogues of these equations in the canonical ensemble are the barometric
formula and the Maxwell-Boltzmann distribution, respectively. In the limit N →
∞, the microcanonical and canonical expressions coincide; however, they differ
for finite N . In particular, in the microcanonical ensemble, the positions and
velocities are not statistically independent. As a result, the kinetic temperature,
defined as the average kinetic energy in a given volume Adz, is nonuniform
throughout the container,
3E mgz
Tkinetic = 1− .
5N − 2 E
where E is the total energy of the microstate and Ξcn the canonical partition function.
In quantum mechanics the density operator and partition function are,
e−β Ĥ
eβ(F −Ĥ) = 1
|ψk ⟩eβ(F −Ek ) ⟨ψk |
P
ρ̂cn = Ξcn = Ξcn k
. (1.164)
Tr e−β Ĥ = e−βF = −βEk
P
Ξcn = ke
The Helmholtz free energy F is constant for the ensemble. However, the prob-
abilities and F will vary if different N, V, T are selected. The free energy F serves
two roles: first, it provides a normalization factor for the probability distribution (the
probabilities, over the complete set of microstates, must add up to one); second, many
important ensemble averages can be directly calculated from the function F (N, V, T ).
70 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
1.4.4.1 Applicability
The canonical ensemble is the ensemble that describes the possible states of a system
that is in thermal equilibrium with a heat bath. It applies to systems of any size;
while it is necessary to assume that the heat bath is very large (i.e. take a macroscopic
limit), the system itself may be small or large.
The condition that the system is mechanically isolated is necessary in order to
ensure it does not exchange energy with any external object besides the heat bath. In
general, it is desirable to apply the canonical ensemble to systems that are in direct
contact with the heat bath, since it is that contact that ensures the equilibrium.
In practical situations, the use of the canonical ensemble is usually justified either
(1) by assuming that the contact is mechanically weak, or (2) by incorporating a
suitable part of the heat bath connection into the system under analysis, so that the
connection’s mechanical influence on the system is modeled within the system.
When the total energy is fixed but the internal state of the system is otherwise
unknown, the appropriate description is not the canonical ensemble but the micro-
canonical ensemble. For systems where the particle number is variable (due to contact
with a particle reservoir), the correct description is the grand canonical ensemble.
∂F
P =− , (1.165)
∂V
the entropy is,
∂F
S = −kB ⟨ln Pcn ⟩ = − , (1.166)
∂T
the partial derivative ∂F/∂N is approximately related to chemical potential, although
the concept of chemical equilibrium does not exactly apply to canonical ensembles of
small systems, and the average energy is,
From the above expressions, it can be seen that, for a given N , the function F (V, T )
has the exact differential,
dF = −S dT − P dV . (1.168)
Substituting the above relationship for ⟨E⟩ into the exact differential of F , an equation
similar to the first law of thermodynamics is found,
dE = T dS − P dV . (1.169)
The energy in the system is uncertain and fluctuates in the canonical ensemble. The
variance of the energy is,
∂E
⟨Ĥ 2 ⟩ − ⟨Ĥ⟩2 = kB T 2 . (1.170)
∂T
Example 17 (Ising model of strongly interacting systems): In a system of
strongly interacting particles, it is usually not possible to find a way to separate
the system into independent subsystems as done in the Boltzmann distribution.
In these systems it is necessary to resort to using the full expression of the
canonical ensemble in order to describe the thermodynamics of the system when
it is thermostatted to a heat bath. The Ising model, which is a widely discussed
toy model for the phenomena of ferromagnetism, is one of the simplest models
showing a phase transition.
where N is the number of particles in the microstate and E is the total energy of the
microstate.
6 In the case where more than one kind of particle is allowed to vary in number, the probability
expression generalizes to Pgc = eβ(Ω+µ1 N1 +µ2 N2 +...−E) , where µj is the chemical potential for the
j-th kind of particles, Nj the number of that kind of particle in the microstate.
72 CHAPTER 1. FOUNDATIONS AND MATHEMATICAL FORMALISM
The quantum mechanics the density operator and partition function are,
e−β(Ĥ−µN̂ )
eβ(Ω+µN̂ −Ĥ) = 1
|ψk ⟩eβ(Ω+µnk −Ek ) ⟨ψk |
P
ρ̂gc = Ξgc = Ξgc k
. (1.172)
−β(Ĥ−µN̂ ) −βΩ β(µnk −Ek )
P
Ξgc = Tr e =e = k e
The grand potential Ω is constant for the ensemble. However, the probabilities
and Ω will vary if different µ, V, T are selected. The grand potential Ω serves two roles:
to provide a normalization factor for the probability distribution (the probabilities,
over the complete set of microstates, must add up to one); second, many important
ensemble averages can be directly calculated from the function Ω(µ, V, T ).
1.4.5.1 Applicability
The grand canonical ensemble is the ensemble that describes the possible states of
an isolated system that is in thermal and chemical equilibrium with a reservoir. The
grand canonical ensemble applies to systems of any size, small or large; it is only
necessary to assume that the reservoir with which it is in contact is much larger
(i.e. to take the macroscopic limit).
The condition that the system is isolated is necessary in order to ensure it has well-
defined thermodynamic quantities and evolution. In practice, however, it is desirable
to apply the grand canonical ensemble to describe systems that are in direct contact
with the reservoir, since it is that contact that ensures the equilibrium. The use of
the grand canonical ensemble in these cases is usually justified either (1) by assuming
that the contact is weak, or (2) by incorporating a part of the reservoir connection
into the system under analysis, so that the connection’s influence on the region of
interest is correctly modeled. Alternatively, theoretical approaches can be used to
model the influence of the connection, yielding an open statistical ensemble.
Another case in which the grand canonical ensemble appears is when considering
a system that is large and thermodynamic (a system that is ’in equilibrium with
itself’). Even if the exact conditions of the system do not actually allow for variations
in energy or particle number, the grand canonical ensemble can be used to simplify
calculations of some thermodynamic properties. The reason for this is that various
thermodynamic ensembles (microcanonical, canonical) become equivalent in some
aspects to the grand canonical ensemble, once the system is very large. Of course,
for small systems, the different ensembles are no longer equivalent even in the mean.
As a result, the grand canonical ensemble can be highly inaccurate when applied to
small systems of fixed particle number, such as atomic nuclei 7 .
Grand ensembles are apt for use when describing systems such as electrons in a
conductor or photons in a cavity, where the shape is fixed but the energy and number
of particles can easily fluctuate due to contact with a reservoir (e.g. an electrical
ground or a dark surface, in these cases). The grand canonical ensemble provides a
natural setting for an exact derivation of the Fermi-Dirac statistics or Bose-Einstein
statistics for a system of non-interacting quantum particles.
7 Note that even in the thermodynamic limit, in the presence of long range interactions, the
1.4.6 Exercises
1.4.6.1 Ex: Thermodynamic potential
What is the most suitable thermodynamic potential to describe a compressible system
(compressible in the sense that the systems always adjusts its pressure to that of
a large environment) with fixed particle number held at constant temperature via
thermal contact with a large reservoir. How is the corresponding ensemble called?
Thermodynamics applied to
fluids and solids
In the present chapter we will apply the notions in thermodynamics acquired in Chp. 1
to various physical systems including multi-component, heterogeneous, chemically
reacting gases and solids.
75
76 CHAPTER 2. THERMODYNAMICS APPLIED TO FLUIDS AND SOLIDS
Let us define intensive variables as extensive ones per particle (or per mole),
∂X
X̃ ≡ , (2.1)
∂N T,P
from now on we designate intensive properties which have been converted from ex-
tensive ones by a tilde (.̃). With the Gibbs free energy G = E + P V − T S defined in
(1.66) we find the differential form (1.152),
S V 1 G̃
= − dT + dP + µ dN − dN = −S̃dT + Ṽ dP ,
N N N N
or
dµ = −S̃dT + Ṽ dP . (2.5)
For heterogeneous systems,
µα = µα (T α , P α ) and µβ = µβ (T β , P β ) . (2.6)
Tα = Tβ ∧ Pα = Pβ =⇒ µα (T α , P α ) = µβ (T β , P β ) , (2.7)
Now, we first hold the temperature fixed, dT = 0, and exploiting Eq. (2.5)
integrate the chemical potential over pressure,
Z P
P
µα (T0 , P ) = µα (T0 , P0 ) + Ṽ α dP ′ = µα (T0 , P0 ) + Rg T ln . (2.10)
P0 P0
Then we hold the pressure fixed, dP = 0, and exploiting the same Eq. (2.5)
integrate the chemical potential over temperature,
Z T
µα (T, P ) = µα (T0 , P ) − S̃ α (T )dT (2.11)
T0
If the temperature-dependent heat capacities of both phases, CPα and CPβ , are
known, as well as the entropy change, ∆S α→β (Tm ) ≡ S β (Tm ) − S α (Tm ), upon phase
transition at a specific temperature Tm , then the entropy difference between two states
in different phases and at different temperatures is simply,
Z Tm α ′ Z T β ′
CP (T ) ′ CP (T ) ′
S̃ β (T ) = S̃ α (T0 ) + ′
dT + ∆ S̃ α→β
(Tm ) + dT . (2.12)
T0 T Tm T′
Such a path is illustrated in Fig. 2.1.
In order to determine the change in chemical potentials between states localized
at different phases and at different temperatures T0 and T , we need to integrate
µα (T, P0 ) from T0 to the phase boundary at Tm and continue integrating µβ (T, P0 )
from Tm to the end point Tm , remembering that within a single phase γ = α, β,
Z CPγ (T )
dCPγ
Z T γ ′
γ γ CP (T ) ′
S̃ (T ) = S̃ (T0 ) + = S̃ γ (T0 ) + dT , (2.13)
γ
CP (T0 ) T T0 T′
78 CHAPTER 2. THERMODYNAMICS APPLIED TO FLUIDS AND SOLIDS
×104
0 (T, P0 )
μ (J)
-1 (Tm , P0 )
(T0 , P0 )
-2
0
100 1
0.5
200 0
T (◦ C) P (bar)
Figure 2.1: (code) Intersecting chemical potential surfaces. The points (Tj , P0 ) can be
related to each other (see text).
and
µγ (T ) T ′′
!
T
CPγ (T ′ ) ′
Z Z Z
γ γ
µ (T ) = µ (T0 ) + γ
dµ = µ (T0 ) − γ γ
S̃ (T0 ) + dT dT ′′ .
µγ (T0 ) T0 T0 T′
(2.14)
We calculate,
Z µα (Tm ) Z µβ (T )
µβ (T ) = µα (T0 ) + dµα + dµβ (2.15)
µα (T0 ) µβ (Tm )
T ′′ T ′′
! !
Tm
CPα (T ′ ) ′ T
CPβ (T ′ ) ′
Z Z Z Z
′′
α
= µ (T0 ) − α
S̃ (T0 ) + dT dT − β
S̃ (T0 ) + dT dT ′′
T0 T0 T′ Tm Tm T′
= µα (T0 ) − S̃ α (T0 )(Tm − T0 ) − S̃ α (Tm )(T − Tm ) − ∆S̃ α→β (Tm )(T − Tm )
Z Tm Z T ′′ α ′ Z T Z T ′′ β ′
CP (T ) ′ ′′ CP (T ) ′ ′′
− ′
dT dT − dT dT
T0 T0 T Tm Tm T′
Z Tm α ′′
CP (T ) ′′
= µα (T0 ) − S̃ α (T0 )(T − T0 ) − ∆S̃ α→β (Tm ) + dT (T − Tm )
T0 T ′′
Z Tm Z T ′′ α ′ Z T Z T ′′ β ′
CP (T ) ′ ′′ CP (T ) ′ ′′
− dT dT − dT dT .
T0 T0 T′ Tm Tm T′
Hence, knowing entropy S̃ α (T0 ) at a specific temperature in one of the phases, the
entropy change between the phases at a specific temperature, ∆S̃ α→β (Tm ), and the
temperature-dependent heat capacities on both phases, CPα (T ) and CPβ (T ), we can
relate the chemical potentials at any temperature of both phases.
2.1. UNARY HETEROGENEOUS SYSTEMS 79
dP ∆S̃ α→β
= , (2.18)
dT ∆Ṽ α→β
which is one form of the Clausius-Clapeyron equation. It states that state changes
along a phase coexistence curve (red line in Fig. 2.1) are ruled by the ratio between
the molar entropy change and the molar volume change. Knowing this ratio, we can
integrate the P (T ) dependence.
In practice, ∆S̃ α→β is not measured in experiments, but rather the heat produced
or absorbed under the transformation (e.g. condensation or evaporation). Since, trans-
formation occurs isobarically under reversible conditions,
δQα→β = ∆H α→β . (2.19)
Recalling that G = H − T S and Gα = Gβ , we have,
∆H̃ α→β = T ∆S̃ α→β , (2.20)
so that
dP ∆H̃ α→β
= , (2.21)
dT T ∆Ṽ α→β
which is the most frequently used form of the Clausius-Clapeyron equation.
In practice, it turns out that the prefactor of the pressure differential is zero for pres-
sure changes below 100000 bar, so that the enthalpy change is very well approximated
by,
d∆H̃ α→β (T, P ) = ∆CPα→β dT . (2.25)
Using an empirical fit equation to describe the heat capacities in both phases γ = α, β,
CPγ = aγ + bγ T + cγ T −2 + dγ T 2 , (2.26)
we find,
Z T
∆H̃ α→β (T, P ) = (∆a + ∆b T + ∆c T −2 + ∆d T 2 )dT (2.27)
T0
= ∆a T + ∆b T 2 − ∆c T −1 + 13 ∆d T −3 + ∆H̃ (0) .
∆H̃ α→β
P 1 1
ln =− − . (2.33)
P0 Rg T T0
and represents a good estimation for phase transitions between gases and fluids or
solids.
2.1. UNARY HETEROGENEOUS SYSTEMS 81
Figure 2.2: (a) P T -phase diagram. (b) Stable and metastable equilibrium lines near a triple
point.
and the same holds for S̃ and H̃, e.g. the enthalpy change upon fusion and vaporization
sum up to the enthalpy change upon sublimation,
The sublimation curve [red line in Fig. 2.2(a)] and the vaporization curve [blue
line in Fig. 2.2(b)] both share a boundary to the gaseous phase, to that we may use
the solution (2.33) of the Clausius-Clapeyron relation for both,
vap
/Rg )(1/Tvap −1/T )
Pvap (T ) = P0 e(∆H̃ (2.36)
(∆H̃ sub /Rg )(1/Tt −1/T )
Psub (T ) = Pt e .
If the boiling temperature Tvap is known, setting P0 = 1 bar obviously fully determines
the vapor curve Pvap = Pvap (T ). In contrast, the sublimation curve Psub = Psub (T ) is
not fully determined as long as the triple point is unknown. However, approximating,
Tt = Tf us (2.37)
and relating the enthalpy change via (2.35), we are a good step further. We just need
to find the triple pressure Pt using the information, that both curves share the triple
point,
Pvap (Tt ) = Pt = Psub (Tt ) = Psub (Tf us ) . (2.38)
82 CHAPTER 2. THERMODYNAMICS APPLIED TO FLUIDS AND SOLIDS
10 1010
(bar)
Tf us
10−2
(Pt , Tt ) 8
10
P
10−3
10−4 106
600 800 1000 1200 1400 300 500 700 900
T (◦ C) T (◦ C)
Figure 2.3: (code) (a) Vapor pressure of strontium obtained with ∆H̃ f us = 144 kJ/mol,
∆H̃ vap = 8.3 kJ/mol, Tf us = 1050 K, and Tvap = 1650 K. The solid lines corresponds to a
path along the phase transition [(see arrow in Fig. 2.2(a)]. The sublimation curve (marked
in red) and the vaporization curve (blue) correspond to those emphasized in Fig. 2.2(a) by
the same colors. The dash-dotted lines show extensions of the phase boundaries beyond the
sublimation curve, respectively, vaporization curve helping to emphasize the discontinuity
at the triple point. (b) Density of the strontium vapor corresponding to the partial pressure
according to the ideal gas law: n = P/kB T .
2.1.3 Exercises
2.1.3.1 Ex: Chemical potential surface
Compute and plot the chemical potential surface: (T, P ) for a monatomic ideal gas in
the range 5 K < T < 1000 K and 10−5 bar < P < 10 bar. Suppose the gas is helium
with S298 = 126.04 J/(mol K) and CP = 5193.2 J/(kg K).
2.1. UNARY HETEROGENEOUS SYSTEMS 83
dP mL
= ,
dT T (Vgas − Vf l )
where L is the latent heat and m the mass of the component undergoing a phase
transition from fluid to gas, can also be derived via a cycle analogous to Carnot’s
cycle. The working fluid is an evaporating liquid, and the efficiency of this fictitious
machine is dT /T , because the temperature difference between theH two isotherms is
dT . From the heat Q used upon evaporation and the work done P dV , which are
related on one hand to the latent heat L and on the other to the volume difference of
Vf l and Vgas , results a rise in vapor pressure dP/dT [see Eq. (1.120)].
How much heat is needed to transform 1 g of ice at −10 ◦ C (cice = 0.55 cal/g/K,
Lf us = 80 cal/g) in a vapor at 100 ◦ C (caq = 1 cal/g/K, Lvap = 540 cal/g)?
A metal bar with specific heat capacity cmt = 0.2 cal/g/K at 100 ◦ C is placed on a
large block of ice at 0 ◦ C. What is the mass of the bar if, when the system reaches
thermal equilibrium, maq = 500 g of ice have melted?
An ice block with the mass mice = 500 g and the temperature −20 ◦ C is put in an
airtight container together with mvap = 200 g of water vapor at 100 ◦ C. What will be
the final temperature of the system?
84 CHAPTER 2. THERMODYNAMICS APPLIED TO FLUIDS AND SOLIDS
This important result is termed the Gibbs-Duhem equation. This equation shows that
intensive properties are not independent but related. When pressure and temperature
are variable, for a system with J components only J −1 components have independent
values of chemical potential 1 .
1 The Gibbs-Duhem equation cannot be used for small thermodynamic systems due to the influence
we must conclude, X
−SdT + V dP = Nj dµj .
j
by holding pressure, temperature, and the number of moles of all other components
fixed. Note that, in contrast to the previous decoration (.̃), which referred to quantities
per total number of particles, the new decoration (.̄) refers to quantities per number
of particles of that species. Then, the total differential form is,
∂X ∂X X
dX = dT + dP + X̄j dNj . (2.44)
∂T P,{Nj } ∂P T,{Nj } j
with the definitions (1.69) of the thermal expansion coefficient and the com-
pressibility,
1 ∂V 1 ∂V
α≡ , κ≡− . (2.46)
V ∂T P,{Nj } V ∂P T,{Nj }
The quantity V̄j is called partial molal volume, and analogous procedures can
be followed for Ēj , S̄j , V̄j , H̄j , F̄j , Ḡj .
86 CHAPTER 2. THERMODYNAMICS APPLIED TO FLUIDS AND SOLIDS
For the particular extensive property X → G and X̄j → µj we derived the Gibbs-
Duhem equation (2.42), but the result holds for any extensive property satisfying
(2.43). E.g. holding temperature and pressure constant, (2.44) becomes,
X
dXT,P = X̄j dNj . (2.47)
j
That is, changes of the partial molal properties of the components add up to a total
change of the system. The value of the extensive state function XT,P is obtained
by integrating (2.47). Fortunately, as an intensive property, X̄j can only depend on
other intensive properties, that is, it cannot depend on Nj . Furthermore, changes of
state functions are path-independent. Hence, the total state XT,P reached by adding
all the components is,
XZ Nj X Z Nj X
XT,P = X̄j dNj = X̄j dNj = X̄j Nj . (2.48)
j 0 j 0 j
which only coincides with Eq. (2.47), if the second term vanishes. Hence,
X
Nj dX̄j = 0 . (2.50)
j
Substituting X̄j → µj , we recover the Gibbs-Duhem equation (2.42) for the case
T, P = cnst.
The important message of the Gibbs-Duhem equation is, that the partial molal
properties are not all independent. Its integration provides a recipe, how to cal-
culate values of partial molal properties of one component from those of the other
components, as we will see later.
Now, X X
∆X mix = X̄jsol − X̄j0 Nj = ∆X̄jmix Nj .
(2.53)
j j
∆X̄jmix measures the change per particle of type j that the state function suffers from
being put into the surrounding composed by all other particle types, and ∆X mix is
the weighted sum of all these changes. Differentiating Eq. (2.53),
X X
d∆X mix = ∆X̄jmix dNj + Nj d∆X̄jmix . (2.54)
j j
Here, the second term is zero because X̄j0 are properties of the reference state, which
is fixed by definition. The summation over the first term is zero by the Gibbs-Duhem
equation (2.50). Hence, X
Nj d∆X̄jmix = 0 , (2.56)
j
which is the Gibbs-Duhem equation applied to the mixing process. With this we
deduce from (2.54),
X
d∆X mix = ∆X̄jmix dNj . (2.57)
j
Let us consider a binary system and rewrite Eq. (2.53) for normalized quantities,
d∆X̃ mix = ∆X̄1mix dη1 + ∆X̄2mix dη2 with dη1 = −dη2 , (2.59)
which gives,
d∆X̃ mix
= ∆X̄1mix − ∆X̄2mix . (2.60)
dη2
Isolating the terms X̄1mix and X̄2mix from the system of equations (2.58) and (2.60),
d∆X̃ mix
∆X̄jmix = ∆X̃ mix + (1 − ηj ) with j = 1, 2 . (2.61)
dηj
Evaluating,
d∆H̄ mix
= a(1 − 2ηj ) ,
dηj
we get,
d∆H̄ mix
∆H̄1mix = ∆H mix + (1 − η1 ) = aη22 ,
dη1
and finally,
∆H̄jmix = aηi̸2=j .
A consistency check yields,
2.2.2.3 Partial molal properties of one component from those of the oth-
ers
The partial molal version of the Gibbs-Duhem relation for the mixing process (2.55)
is, X
ηj d∆X̄jmix = 0 , (2.62)
j
Now, limη2 →0 X̄2sol = X̄20 , since there is nothing to mix. Hence, limη2 →0 ∆X̄2mix = 0.
Integrating the last expression,
∆X̄2mix η2
η2 d∆X̄2mix
Z Z
η2
∆X̄1mix =− d∆X̄2mix = − dη2 . (2.64)
0 η1 0 η1 dη2
Do the Exc. 2.2.5.2.
∆H2mix = aη12 .
to any total property X yields the corresponding molal property. Applying this
operator to the definitions, laws, coefficient relations, and Maxwell equations, we
obtain the corresponding molal expressions. For example, the counterpart of the
relation H = E + P V is,
∂H ∂E ∂V
H̄j ≡ = +P ≡ Ēj + P V̄j , (2.66)
∂Nj T,P,Nk̸=j ∂Nj T,P,Nk̸=j ∂Nj T,P,Nk̸=j
and analogously for all other equations. For a solution with various components, we
can simply substitute any extensive variable X for each component j,
X −→ X̄j . (2.67)
with
∂E
µj ≡ . (2.69)
∂Nj S,V,Nk̸=j
Comparing this expression with the first law, we identify,
J
X
δArev = µj dNj (2.70)
j=1
X
P j ≡ ηj P such that P = Pj . (2.77)
j
dT =0
dµj = dḠj = −S̄j dT + V̄j dP −→ V̄j dP , (2.78)
since temperature does not change. The partial molal volume can be evaluated from,
∂V
V̄j = (2.79)
∂Nj T,P,Nk̸=j
! !
Nk kB T
∂ N kPB T
P
ideal ∂ k P kB T
−→ = = ,
∂Nj ∂Nj P
T,P,Nk̸=j T,P,Nk̸=j
assuming ideal gases. With this we can now proceed to integrate the chemical poten-
tial,
Z Pj
∆Ḡj = ∆µj = µj − µ0j = V̄j dPT (2.80)
P0
Pj
= kB T ln = kB T ln ηj = kB T ln aj .
P0
We see that, for ideal gas solutions the abundance and the activity are identical, which
is to say the activity coefficient is γj = 1. Based on the expressions (2.73), we can
92 CHAPTER 2. THERMODYNAMICS APPLIED TO FLUIDS AND SOLIDS
also calculate the dependence of other partial molal properties on the abundances,
6 0
(J/mol)
(J/mol)
(a)
T 2000 (c) T
all T
(J/mol K)
-2000
4 0
-4000
ΔGmix , ΔF mix
ΔGmix , ΔF mix
-2000
-6000
2
ΔS mix
-4000
-8000
-6000
(b)
0 -10000
0 0.5 1 0 0.5 1 0 0.5 1
η2 η2 η2
Figure 2.5: (code) Dependence of state functions on composition and temperature in the
range between 300 K and 1700 K. (a) The entropy for an ideal gas and a regular solution
are identical. (b) Gibbs free and Helmholtz free energy for an ideal gas solution and (c) for
a regular solution.
For real solutions for which the ideal gas equation does not hold, the partial molal vol-
ume cannot be determined from Eq. (2.79), but must be experimentally determined.
It is often interesting to define a ’deviation volume’,
kB T
Vj ≡ V̄j − . (2.82)
P
2.2. MULTI-COMPONENT, HOMOGENEOUS, NON-REACTING SYSTEMS 93
As the deviation from ideal gas behavior, measured by Vj diminishes, the fugacity
of component j approaches its partial pressure. Measurement of the fugacity of one
component over a range of temperature, pressure, and composition is sufficient to
describe the behavior of real gas mixtures completely in that range.
fj = P0 aj . (2.85)
From (2.80),
The activity coefficient γj = γj (T, P, ηj ) quantifies the deviation from ideal solution
behavior; if γj > 1 then the component j ’acts’ as if its abundance were more than
expected from a supposed ideal gas behavior.
Using the expressions for the thermodynamic potentials of the first column of
Tab. (2.81), we can also express the other partial molal properties (PMP) in terms of
activity or activity coefficient,
Because, the excess mixing entropy vanishes ∆S̃ mix,xs ≡ ∆S̃ mix,rl −∆S̃ mix,id = 0,
the excess Gibbs free energy G̃ = H̃ −T S̃ becomes equal to the excess mixing enthalpy,
∆G̃mix,xs = (∆H̃ mix,rl − T ∆S̃ mix,rl ) − (∆H̃ mix,id − T ∆S̃ mix,id ) (2.90)
mix,xs
= ∆H̃ .
Since S̃ = (∂ G̃/∂T )P,Nj , we conclude that ∆G̃mix,xs = ∆H̃ mix,xs are temperature-
independent. For the components, we have from Eq. (2.80),
mix,xs
∆Ḡmix,xs = kB T ln γj =⇒ γj = e∆H̄j /kB T
. (2.91)
Since, all properties of a solution can be calculated from a known activity coefficient,
the regular solution model with focuses upon the heat of mixing as a function of
composition.
= kB T (η1 ln η1 + η2 ln η2 ) + a0 η1 η2 .
A solution is called ’non-regular’ when the coefficients of the solution model ad-
ditionally depend on temperature, for example,
The presence of collisions different particles motivates the ansatz for the mixing en-
tropy made in example 23,
∆H mix = aη1 η2 . (2.94)
2.2.5 Exercises
2.2.5.1 Ex: Partial pressures
A closed cylindrical reservoir with the base area S = 10 cm2 is kept at a constant
temperature T = 27 ◦ C. It is divided in two volumes by an airtight mobile disk with
the mass m = 10 kg. The upper volume VO2 contains ηO2 = 1 mol of oxygen, the
lower volume VN2 contains the same amount of nitrogen. Due to its weight the disc
finds an equilibrium position when the lower volume is VN2 = 10 l.
a. What are the masses mO2 and mN2 of the gases?
b. What are the pressures PO2 and PN2 ?
c. What is the upper volume VO2 ?
d. What are the densities nO2 and nN2 ?
e. Now the disc has a hole, so that the gases can mix and the disc falls to the bottom
of the reservoir. What is the final pressure of the mixture?
∆Ḡ2 = Rg T ln η2 .
Use the Gibbs-Duhem integration (2.64) to derive the corresponding relation for com-
ponent 1.
Now, we assume that each component exists in Γ different phases. Each phase viewed
as a system exchanges heat, work, and matter with the other phases and with the
reservoir,
J
X
dE γ = T γ dS γ − P γ dV γ + µγj dNjγ , (2.96)
j=1
for all phases γ = α, β, .... For the extensive properties V, S, E, H, F, G the value for
the property of the system is the sum of the parts,
Γ
X
X syst = Xγ . (2.97)
γ=1
When the system is taken through an arbitrary change of state, then the change of
X syst is simply the sum of the changes that each phase experiences,
Γ
! Γ
X X
syst γ
dX =d X = dX γ . (2.98)
γ=1 γ=1
2.3. MULTI-COMPONENT, HETEROGENEOUS, NON-REACTING SYSTEMS97
Figure 2.6: Heat work, and particle exchange in a micro structure with 3 components and
two phases.
To find the equilibrium conditions, we need to express the entropy. We solve the
Eq. (2.96) for dS γ and sum up the entropies for all phases,
Γ Γ γ J
X X 1 P 1 X
dS syst = dS γ = dE γ +
γ γ
dV γ − γ µγj dNjγ . (2.100)
γ=1 γ=1
T T T j=1
The condition for equilibrium is S syst = maximal. Eq. (2.100) contains Γ×(2+J)
variables, but their number is reduced if the system is isolated from the reservoir,
dE syst = dV syst = dN syst = 0. Let us suppose there are only two phases γ = α, β.
Then the isolation assumption reads,
J
Pα Pβ
syst,iso 1 1 α α
X 1 α 1 β
dS = α
− β dE + − dV − µ − β µj dNjα
α j
T T Tα Tβ j=1
T T
!
=0. (2.102)
Tα = Tβ thermal equilibrium
α β
P =P mechanical equilibrium (2.103)
µα
j = µβj , ∀j ≤ J chemical equilibrium
98 CHAPTER 2. THERMODYNAMICS APPLIED TO FLUIDS AND SOLIDS
This rule can be easily generalized to more than two phases by pairwise comparison
of the phases. We immediately see that,
These equations form the basis for the construction and calculation of multi-component
multi-phase phase diagrams.
2.3.3 Exercises
2.3.3.1 Ex: Volume change in a multi-phase multi-component system
Use Eq. (2.98) to write out a general expression for the change of volume of a three-
phase two-component system including all 12 terms.
in the whole system. Changes of this property along a thermodynamic process are
expressed as usual, Z
dXV = dx(r)δV , (2.109)
V
dδSV X dδNjV X
=T + µj = T ds(r) + µj dnj (r) ,
δV j
δV j
where s is the entropy density and nj the molar density of the component j. That is,
we can derive local versions for all thermodynamic relationships.
If during the thermodynamic process the system is isolated against energy and particle
exchange,
Z Z
0 = dE = de(r)δV and 0 = dNj = dnj (r)δV . (2.112)
V V
Equilibrium is reached when the entropy is maximum, and a necessary condition for
that is, X
dS − kB αj dNj − kB βdE = 0 (2.113)
j
for any value of the Lagrange multipliers αj and β. Substituting (2.111) and (2.112),
Z Z Z
1 1 X X
0= de(r) − µj dnj (r) δV − kB β
de(r)δV − kB αj dnj (r)δV ,
V T T j V j V
(2.114)
which implies,
1 X µj
0= − kB β de(r) − + kB αj dnj (r) . (2.115)
T j
T
2.4. CONTINUOUS, NON-UNIFORM SYSTEMS EXPOSED TO EXTERNAL FORCES101
We conclude,
1 µj
β= , αj = − . (2.116)
kB T kB T
Since αj and β are constants, this means that temperature and chemical potential
cannot depend on position,
∇T = 0 = ∇µj . (2.117)
To derive the condition for mechanical equilibrium, we apply the operator (2.43)
on (2.42), that is,
∂ X
applied to Nj dµj = −SdT + V dP − Nk dµk , (2.118)
∂Nj P,T,Nk̸=j
k̸=j
Multiplying both sides with ηj and summing over all components, and exploiting the
fact
P that summation over all partial molar volumes produces the total molar volume,
j ηj V̄j = V ,
X X ∂∇µk
V ∇P = ηj Nk . (2.122)
j
∂Nj P,T,Nk̸=j
k̸=j
The right-hand side can be shown to vanish using the Gibbs-Duhem equation (2.50).
Hence,
∇P = 0 . (2.123)
The results (2.117) and (2.121) were expected and only tell us that there is nothing
wrong with the density formalism.
The work executed by this force, δWext = dWext = F · dr, then does not depend on
the path taken between start and end point,
Z r2 Z r2 Z Φ(r2 )
Wext = F · dr = − ∇Φ(r) · dr = − δΦ(r) = Φ(r1 ) − Φ(r2 ) . (2.125)
r1 r1 Φ(r1 )
The total potential energy is the sum of the potential energies to which the mass
elements δm are subject,
Z Z Z
δΦ(r) δΦ(r) δΦ(r) X
Epot = δV = δm(r) = Mj nj δV . (2.127)
V δV V δm V δm j
which implies,
X µj
1 δΦ(r)
0= − kB β de(r) − + kB β Mj + kB αj dnj (r) . (2.130)
T j
T δm
We conclude,
1 µj δΦ(r)
T = , + kB β Mj + kB αj = 0 . (2.131)
kB β T δm
To derive the condition for mechanical equilibrium, we use again the expres-
sion (2.120), but now inserting the equilibrium conditions (2.132), we obtain an equa-
tion generalizing (2.121),
X
δΦ(r) ∂∇µk
V̄j ∇P = −∇ Mj − Nk . (2.133)
δm ∂Nj P,T,Nk̸=j
k̸=j
Again, multiplying both sides with ηj and summing over all components,
X δΦ(r)
V ∇P = −∇ η j Mj , (2.134)
j
δm
or,
M δ∇Φ(r)
∇P = − . (2.135)
V δm
In particular, for an ideal gas,
M ideal M̃ g
dP = − gdz −→ − P dz . (2.142)
V Rg T
This is a first-order differential equation in P with the solution,
2.4.3 Exercises
2.4.3.1 Ex: Pressure in a harmonically trapped ideal gas
Calculate the local pressure in a harmonically trapped ideal gas in thermodynamic
equilibrium.
2 H2 + O2 −−→ 2 H2 O . (2.158)
106 CHAPTER 2. THERMODYNAMICS APPLIED TO FLUIDS AND SOLIDS
2 H2 + O2 −−→ 2 H2 O (2.160)
H2 + O2 −−→ H2 O2 .
2 H2 O + O2 −−→ H2 O2 (2.161)
are linear combinations of (2.160). Systems with R > 1 are called bivariant, respec-
tively, multivariant.
In the next section we will first focus on the gas phase, where the components
unambiguously exist as molecules. In solids, which are characterized by the existence
of multiple bonds between molecules, this is more complicated. Also, we will first
treat univariant systems before progressively generalizing to multivariant systems.
known, then the chemical potential of the components and hence the affinity can be
computed.
The affinity can be positive or negative, depending on the chemical potentials of
the reactants and products. The third law of thermodynamics, however, requests
that for any process dS iso > 0, so that in Eq. (2.164) the affinity and dNH2 O must
have opposite signs. In other words, the affinity decides in which direction a chemical
reaction will take place. A state of equilibrium is reached, when dS iso > 0, which
implies A = 0. The chemical reaction comes to a halt, when the sum of the chemical
potentials of the reactants and products are equal, and this condition may require an
excess of reactants or products, as illustrated in Fig. 2.7.
Figure 2.7: (a) Direction of a chemical reaction as a function of affinity. (b) The equilibrium
of a chemical reaction may be on the left or right side.
where G̃0j is the molar Gibbs free energy of component j, when it is in its reference
state. Let us now consider an arbitrary chemical reaction,
p P + q Q −−→ x X + y Y , (2.166)
for which,
where the abbreviation ∆G0 comprises the four Gibbs free energy terms and describes
the change in Gibbs free energy upon complete conversion of p moles of P and q moles
of Q in their standard states into x moles of X and y moles of Y in their standard
states. The quotient in the logarithm is called the ratio of activities for the reaction,
aX aY
Q≡ . (2.168)
aP aQ
As already mentioned, a reaction will go on until the system reaches equilibrium,
A gas mixture may be considered an ideal solution, so that the activities of their
components are given by their molar fractions,
a2H2 O ηH2
2O
Q= 2
= 2
≈ 3.1 × 105 .
a2H2 aO2 η2H2 ηO2
Thus, Q/Qeq ≈ 7.2 × 10−19 , which means that the equilibrium is far on the
H2 O side. In the reaction process, the 0.01 mol of hydrogen will be completely
used up. However, there is an excess of oxygen in the system, as only 0.005 mol
are used. Hence, the final composition will be ηH2 ≈ 0, ηO2 = 0.0025 mol,
and ηH2 O = 0.97 mol. Note also that the total number of moles is reduced to
ηH2 + ηO2 + ηH2 O = 0.995 mol.
1. Identify the types j of molecules involved and their numbers Njin (or moles or
molar fractions) before a chemical reaction is initiated.
3. Calculate the total number of atoms of each atomic element in the system from,
J
X
Ak = Akj Njin ∀ k = 1, ..., K . (2.170)
j=1
4. The numbers of atoms per species Akj and the chemical composition of the
mixture Nj change when the atoms reorganize into different molecules through
chemical reactions (labeled r = 1, ..., R). Thus, in parallel, we need to identify
all R possible chemical reactions transforming a group of molecules {jr } into a
group {jr′ }, and vice versa, with jr , jr′ = 1, ..., J labeling the types of molecules.
independent reactions in this system, and each reaction has its equilibrium
condition,
Q eq
0 eq j ′ aj ′
∆Gr = −Rg T ln Qr = −NA kB T ln Q r eqr (2.171)
jr ajr
Q eq Q eq
ideal jr′ ηjr′ jr′ Njr′
−→ −NA kB T ln Q eq = −NA kB T ln Q eq ∀ r = 1, ..., R ,
jr ηjr jr Njr
for an ideal solution. Thus, each of the equilibrium constants Qeq r ruling a
chemical reaction r is related to the composition of the system through a non-
linear equation.
6. The number of atoms Ak of a given atomic element is conserved and does not
change with time during the reaction process. Hence, the K equations (2.170)
also hold after chemical equilibrium is reached. Together with R equations
(2.171) we obtain a sufficient set of equations,
PJ
j=1 Akj Njeq = Ak ∀ k = 1, ..., K
(2.172)
∆G0
ln Njeq eq
P P
jr′ ′ − jr ln Njr = − Rg Tr ∀ r = 1, ..., R
r
While the first equations are linear in the J unknown values Njeq , the second
equations are non-linear.
Note, that the equilibrium composition of the mixture depends on the initial
composition that the system had while it was still isolated, but only via the numbers
of atoms of all involved atomic elements, no matter how these were distributed among
the molecular components before these were mixed. The task is to express the final
composition Njeq as a function of the Ak and the Qeq r .
2 H2 + O2 ←−→ 2 H2 O (i)
2 CO + O2 ←−→ 2 CO2 (ii) . (2.178)
CH4 + 2 O2 ←−→ 2 H2 O + CO2 (iii)
The standard free energy for these reactions at the specified temperature are,
With (2.177) and (2.180) there are thus seven equations for the six unknown
equilibrium molar fractions and the ratio of total numbers of molecules before
in eq
and after the reactions Ntot /Ntot , that can be solved numerically. A MATLAB
code which is available here, yields,
2.5.3 Exercises
2.5.3.1 Ex: Final composition of an ideal gas mixture
An ideal gas at 1000 K has the following composition,
2.6.1 Solid-liquid-vapor
In the case of the liquid-vapor transition the two phases are only quantitatively dis-
tinct, but have the same symmetry. Therefore, a discontinuity of the thermodynamic
potentials is required to reveal the phase transition.
In the case of the solid-liquid transition the two phases are qualitatively distinct
due to different symmetries. We do not need a discontinuity to distinguish the phases.
Landau’s theory holds for this class of transitions. It establishes a relationship between
symmetry considerations and physical characteristics by introducing the notion of the
order parameter and free energy.
112 CHAPTER 2. THERMODYNAMICS APPLIED TO FLUIDS AND SOLIDS
2.7 Materials
2.7.1 Electrons in solids
2.7.1.1 Types of solids
In contrast to a gas, which in most cases consists of isolated particles, the interparti-
cle interaction plays a dominant role in crystals. Solids, or more specifically crystals,
are classified according to the predominant type of binding. 1. Molecular binding
is responsible for the solidification of binary gases like O2 . Here, fluctuating dipole
moments inducing dipole moments in neighboring molecules lead to van der Waals
attractive forces on the order of Ebind ≃ 10−2 eV going like r−7 . 2. Ionic binding
gives rise to periodic structures alternating positive and negative charges, as in NaCl.
3. Covalent binding is directional. This directionality determines the crystalline struc-
ture, such as in graphite and diamonds. In those three binding types there are no
free electrons and hence no conductance. However, covalently bound crystals can
sometimes be semiconductors or transparent. 4. Metallic binding is a limit of cova-
lent binding, in which the valence electrons, which shared by all atoms, overrule the
repulsion between the ions. The ionic lattice is immersed in a gas of free electrons.
The ions have filled shells and are spherically symmetric. The electrons can easily
absorb light, which makes the crystal opaque. The type of binding is studied via
X-ray diffraction, via the dielectric properties, etc..
The number of orbitals in the isolated atoms forming the crystal gives the number
of states available to the free electron gas. The exchange interaction of the fermionic
electrons lifts the degenerescence (generalization of the H2 molecule) and gives rise to
a band structure. The electronic localization determines the width of the band: very
delocalized electrons move in large bands. The interatomic distance also influences
the band width. The closer the atoms the stronger the interaction, the larger the
bands.
Bands connecting to different orbital may finally overlap. Note that the ml degen-
eracy is lifted because spherical symmetry is broken by the crystal [E(3s) ̸= E(3p)].
2.7. MATERIALS 113
j eE λ
= vd = × , (2.182)
ne m v̄
where j is the current density and ne the charge density. The fact that v̄ and λ do
not depends on the electric field is known as Ohm’s law. The mobility µ ≡ vd /E
allows to write the electrical conductance as
ρ−1 = n− e− µ− + n+ e+ µ+ . (2.183)
The value and sign of the Hall coefficient 1/ne can be measured by the Hall effect. It
is positive if the conductance occurs primarily through holes and negative if it occurs
through electrons.
2.7.1.4 Semiconductors
There is an intrinsic temperature-dependent conductivity (for Si ρ(600 K)/ρ(300 K) ≈
109 ). Extrinsic conductivity can be induced by photoexcitation or doping. E.g. Ar/Ga
in a Ge crystal has one weakly bound electron more/less than required to fit into the
lattice. This generates discrete energy levels slightly below the conducting/above the
valence band, min(Ec ) − En , Ep − max(Ev ) ≃ 0.01 eV.
The Fermi energy EF is the energy, where half of the electrons are below that.
In an isolator is between max(Ev ) and min(Ec ). In the presence of doping EF is
shifted by the additional amount of electrons/holes toward En , Ep . If n and p-doped
materials are combined, electrons drift from the n to the p region, such as to minimize
energy and obtain a uniform EF across the hole crystal.
Thermally excited electrons may drift and recombine with holes. The junction
is maintained by a steady flux in a dynamic equilibrium. An external voltage can
higher/lower the barrier, because the potential drops mostly near the junction, where
the resistance is highest. In this case the thermal current is not equilibrated, the
diode either blocks or opens. The electrons move to try to rectify EF .
A transistor is a series of junctions in npn of pnp configuration. The base-emitter
current can be used to switch a collector-emitter current by injecting electrons. A
tunnel diode acts like a normal diode except that when the bands come closer together
within the junction (at low voltages in conduction polarization), electrons may pass
114 CHAPTER 2. THERMODYNAMICS APPLIED TO FLUIDS AND SOLIDS
by tunneling from the conducting into the valence band. This flow gradually stops
when EF is leveled (for zero voltage). Tunneling currents react much faster than
thermal drift currents.
2.7.2 Plasmas
2.7.2.1 Debye length
Consider a mixture of charges + and −, that is, a plasma. Energy seeks to be
minimized by local compensation of charge imbalance. However, thermal motion
spoils perfect homogeneity. That is, if on the one hand, looking at large scales, the
environment seems neutral and homogeneous, at small scales there may be charge
imbalances producing potential sites with exponentially decreasing ranges,
ne2
1 1 1 1 1
= + = + . (2.184)
λD λD− λD+ ε0 kB T+ kB T−
The Debye length naturally enters the thermodynamic description of large sys-
tems of mobile charges. We consider a system of 2 different species of charges q±
and n± (r) at locations r. According to the so-called primitive model, these charges
are distributed in a continuous medium characterized only by its relative static per-
mittivity, εr . This distribution of charges through the medium generates an electric
potential Φ(r) that satisfies the Poisson equation:
where ε ≡ εr ε0 , ε0 is the dielectric constant, and ρE is the charge density outside the
medium (logically, not spatially).
The mobile charges do not only generate Φ(r), but also are moved according to the
associated Coulomb force, −q± ∇Φ(r). Assuming the system to be in thermodynamic
equilibrium with a heat reservoir at an absolute temperature T , the concentrations
of discrete charges, n± (r), can be considered as thermodynamic averages (ensem-
ble average) and the associated electrical potential as a thermodynamic mean field.
With these assumptions, the concentration of species j is described by the Boltzmann
distribution,
n± (r) = n0± e−q± Φ(r)/kB T , (2.186)
where n0j is the mean field concentration of the charge species j.
Identifying the instantaneous concentrations and the potentials in the Poisson
equation with their mean-field counterparts in the Boltzmann distribution, we obtain
the Poisson-Boltzmann equation:
ε∇2 Φ(r) = −q+ n0+ e−q+ Φ(r)/kB T − q− n0− e−q− Φ(r)/kB T − ρE (r) . (2.187)
Solutions of this nonlinear equation are known for simple systems. Solutions for more
general systems can be obtained in the high-temperature (or low-coupling) limit,
qj Φ(r) ≪ kB T , by Taylor expansion of the exponential,
q± Φ(r)
e−q± Φ(r)/kB T ≈ 1 − . (2.188)
kB T
2.7. MATERIALS 115
n0+ q+
2
n0 q 2
2
ε∇ Φ(r) = + − − Φ(r) − n0+ q+ − n0− q− − ρe (r) (2.189)
kB T kB T
also known as Debye-Hückel equation. The second term on the right side disappears
for electrically neutral systems. The term in parentheses divided by ε, has the unit
1/m2 . By a dimensional analysis, it leads to a definition of a characteristic length
scale,
1/2
εkB T
λD = (2.190)
n0+ q+
2 + n0 q 2
− −
usually called Debye-Hückel length. Being the only characteristic length scale of the
Debye-Hückel equation, λD defines the scale of variations in the potential and the
concentrations of the charged species. All charged species contribute to the Debye-
Hückel length in the same manner regardless of the charge signal. For an electrically
neutral system, the Poisson equation is,
ρe (r)
∇2 Φ(r) = λ−2
D Φ(r) − . (2.191)
ε
To illustrate the Debye shielding, the potential produced by an external point-like
charge ρe = Qδ(r) is,
Q −r/λD
Φ(r) = e . (2.192)
4πεr
The bare Coulomb potential is exponentially shielded by the medium over a distance
corresponding to the Debye length.
The length of Debye-Hückel can be expressed in terms of the length of Bjerrum
λB as,
−1/2
XN
λD = 4πλB n0j zj2 , (2.193)
j=1
where zj = qj /e.
where T± are the temperatures of the electrons and ions, n− is the density of the
electrons and n+ that of the atomic species j, with positive ionic charge z+ qe . The
ion term is often neglected, giving,
s
ε0 kB Te
λD = , (2.195)
ne qe2
although this is valid only, when the mobility of ions is negligible on the time scale of
the process.
2.7.3 Exercises
Appendices to
’Thermodynamics’
3.1 Quantities and formulas in thermodynamics
3.1.1 Statistical formulas
Stirling’s formula is,
3.1.2 Polylogarithm
The polylogarithm (or Joncquière’s function) is a function defined as,
∞ ∞
Zt xη−1 dx
Z
X 1
Liη (Z) ≡ = . (3.2)
t=1
tη Γ(η) 0 Z −1 ex − 1
It serves to express the Bose and Fermi functions used in quantum statistics,
The upper sign holds for bosons, the lower for fermions.
117
118 CHAPTER 3. APPENDICES TO ’THERMODYNAMICS’
where the second equation represents an expansion. The derivative satisfies a useful
relationship,
∂gξ± (Z) ∞
X ∂ (±Z)ℓ
∞
X ±(±Z)ℓ−1
∞
1 X (±Z)ℓ
±
gξ−1 (Z)
= ξ
= ξ−1
= ξ−1
= . (3.6)
∂Z ∂Z ℓ ℓ Z ℓ Z
ℓ=0 ℓ=1 ℓ=1
d −Z ex
= . (3.7)
dx Z −1 ex + 1 (ex /Z ∓ 1)2
We calculate,
∂gξ± (Z) ∞ ∂
−xξ−1 ∂Z (Z −1 ex ∓ 1) ∞
ex
Z Z
1 1
= −1
dx = 2 xξ−1 dx
∂Z Γ(ξ) 0
x
(Z e ∓ 1) 2 Z Γ(ξ) 0 (Z −1 ex ∓ 1)2
(3.8)
Z ∞ ∞
xξ−1 ex 1 −Z
= x − 2 (ξ − 1)xξ−2 −1 x dx
(e /Z ∓ 1)20 Z Γ(ξ) 0 Z e ∓1
Z ∞ ±
1 xξ−2 gξ−1 (Z)
=0+ −1 x
dx = .
ZΓ(ξ − 1) 0 Z e ∓ 1 Z
Another useful relationship is the Sommerfeld expansion, which holds for Fermi
functions,
∞ y ∞ x
η(y + x)ξ−1 dx η(y − x)ξ−1 dx
Z Z Z Z
η(x)dx
= η(x)dx + − (3.9)
0 ex−y + 1 0 0 ex + 1 0 ex + 1
Z y 2
π ′
≈ η(x)dx + η (x) + ...
0 6
2.5
L (J/kg)
2
fusion
vaporization
1.5
-50 0 50 100 150 200
T (◦ C)
Figure 3.1: (code) Specific latent heat of fusions and condensation of water.
121
Chapter 4
Statistical thermodynamics
All thermodynamic quantities studied in Chp. 1 (extensive or intensive) are quasi-
continuous, i.e. macroscopic. The laws of thermodynamics found to rule the behavior
of large systems were discovered empirically via experimental observations. The na-
ture of the laws is thus phenomenological, i.e. not derived from first principles. Until
now we totally neglected the fact that matter (gases, fluid, or solids) is composed of
microscopic elementary particles (atoms or molecules). Nevertheless, it already be-
came clear that the behavior of a system is somehow related to the properties of the
particles that compose it. E.g. the degrees of freedom of a molecule that can be ex-
cited have an influence on the heat capacity of a gas composed of these molecules; the
Joule-Thomson effect is due to intermolecular forces; and what we experience as heat,
is actually an outward manifestation of molecular and atomic motion, as we already
pointed out in Sec. 1.1.2. Tracing back macroscopic properties and phenomena to mi-
croscopic models bears a formidable potential of deepening our level of understanding
thermodynamic systems. It may even provide insight into the physical meaning of
mysterious or elusive phenomenological concepts such as entropy production.
An atomistic description acknowledges the fact that matter is quantized into small
portions called molecules 1 . Each molecule is understood as a (not necessarily rigid)
body characterized by its center-of-mass coordinates, but also its rotations or internal
vibrations. With typically 1023 atoms in just one liter of air the task of describing
the microstate by all its coordinates is hopeless. The mathematical discipline that
provides the tools capable of handling such big numbers is statistics, and the primary
tool supplied for the purpose is the concept of the distribution function. The idea is
to lump atoms having similar properties together to classes, e.g. energy levels. The
distribution function then simply reports the number of particles in each class, which
dramatically reduces the amount of information. The task of statistical thermodynam-
ics is now the description of a thermodynamic state in terms of a distribution function
called macrostate. The formulation of statistical thermodynamics by Boltzmann and
Gibbs provided a solid microscopic foundation of phenomenological thermodynamics.
We will begin this chapter with a calculation of the Boltzmann distribution of
microstates over the macrostates in Sec. 4.1 and introduce the concept of partition
function, from which all macroscopic state functions may be computed. As appli-
cations of this algorithm we will revisit the ideal gas and the Einstein model of a
1 The ’quantization’ of matter is not to be understood in the quantum mechanical sense. Nev-
ertheless, the particles themselves are generally microscopic and, under certain circumstances, may
behave following rules dictated by quantum mechanics. This can lead to macroscopically observable
phenomena studied in the area of quantum statistics, as we will learn in Chp. 4.2.
123
124 CHAPTER 4. STATISTICAL THERMODYNAMICS
crystalline solid.
Of all possible macrostates, there will be one containing the largest number of
microstates, and the probability to encounter the system in this macrostate is highest.
Examination of P{nj } for a variety of macrostates {nj } reveals that the probability
distribution is sharply peaked, and that macrostates deviating only slightly are already
very unlikely. The most probable state is now interpreted as the state of equilibrium,
and this hypothesis forms the basis for connecting phenomenological thermodynamics
to an atomistic statistical description.
4.1. MICROSTATES, MACROSTATES, AND ENTROPY 125
Figure 4.1: Illustration of micro- and macrostates with identical indistinguishable: (a-c) Dis-
tribution of 12 particles over 2 boxes. (d-f) Distribution of 13 particles over 4 energy levels.
All schemes show different microstates, but only the schemes (a) and (b), respectively, (d)
and (e) correspond to same macrostates.
S = kB ln W . (4.3)
This expression allows to compute the entropy of any macrostate of the system.
To find the equilibrium macrostate {n1 , n2 , .., nj , .., nr }eq in the atomistic descrip-
tion, we have to maximize the entropy (4.5). That is, we have to evaluate the to-
tal differential of entropy in the direction of changes {dn1 , dn2 , .., dnj , .., dnr } of the
126 CHAPTER 4. STATISTICAL THERMODYNAMICS
Pr
macrostate under the constraint N = j=1 nj ,
r r r
X ∂S ∂S X nj X nj
dS = dnj + dN = −kB 1 + ln dnj +kB dN , (4.6)
j=1
∂nj ∂N j=1
N j=1
N
yielding,
r
X nj
dS = −kB ln dnj . (4.7)
j=1
N
The maximum of the entropy function (4.5) under the constraints (4.9) can be
found using the technique of Lagrange multipliers, which consists in solving the equa-
tion
r
X nj
0 = dS − αkB dN − βkB dE = kB − ln − α − βεj dnj (4.10)
j=1
N
leaving us with,
nj e−βεj
= Pr −βεj
, (4.13)
N j=1 e
2 Note that dε = 0, if the energy levels do not vary along a thermodynamic process, only their
j
population with particles.
4.1. MICROSTATES, MACROSTATES, AND ENTROPY 127
1 P µ
dS = dE + dV − dN , (4.16)
T T T
where µ is the chemical potential per atom and dV = 0, since we assumed in this
derivation, that every atom has access to the whole volume of the system. A compar-
ison of the expressions (4.15) and (4.16) then yields,
1
β= and α = −βµ = ln Ξcn . (4.17)
kB T
r
nj 1 −εj /kB T X
= e with Ξcn = e−εj /kB T . (4.18)
N Ξcn j=1
mann distribution holds for any ensemble of classical particles, we can use it to derive the distribution
function for canonical ensembles.
128 CHAPTER 4. STATISTICAL THERMODYNAMICS
∂ 2 ln Ξcn
∂E ∂ ln Ξcn
CV = = 2kB T + kB T 2 . (4.23)
∂T V ∂T V ∂T 2 V
r
X
Ξcn = e−βεj = 1 + e−βε , (4.24)
j=1
n1 1 1 n2 e−βε e−βε
= = , = = . (4.25)
N Ξcn 1 + e−βε N Ξcn 1 + e−βε
In particular, the ratio between populations of consecutive levels is, n2 /n1 = e−βε .
At low temperature, kB T ≪ ε, the excited state population is negligibly small, while
at high temperature, kB T ≫ ε, both energy levels have almost the same population.
Do the Exc. 4.1.7.8.
With the partition function it is easy to evaluate the potentials,
εj = (j + 32 )ℏω . (4.28)
The normal-mode frequency ω is related to the spring constant of the atomic bond
and the atomic mass. The spectrum (4.27) completely defines the model.
r r ∞
X X X e−3βℏω/2
Ξcn = e−βεj = e−3βℏω/2 e−βℏωj ≃ e−3βℏω/2 e−βℏωj = . (4.29)
j=1 j=1 j=0
1 − e−βℏω
The discrete energies nℏω are identified with quasi-particles called phonons. The
quantum nature of atoms does not matter, they just provide the medium supporting
the phonons.
130 CHAPTER 4. STATISTICAL THERMODYNAMICS
where v is the velocity of sound propagation, the formula is totally equivalent to the
density-of-states for photons in a cavity. Assuming thatR there is an upper bound νm
ν
for the vibrational frequencies, we normalize as 3N0 = 0 m ρ(ν)dν. The energy now
is 4 ,
Z νm
T 4 θ/T x3 dx
Z
ℏω 4πV
E= νdν = 9N kB 3 . (4.35)
0 eℏω/kB T − 1 v 3 θ 0 ex − 1
4 The fact that the electron gas also has a heat capacity is neglected.
4.1. MICROSTATES, MACROSTATES, AND ENTROPY 131
The Debye temperature θ = hνm /kB is characteristic for the metal. The derivative is
then, " Z #
3 θ/T
T x3 dx θ 1
CV = 9N kB 4 − . (4.36)
θ 0 ex − 1 T eθ/T − 1
At low temperatures this formula reproduces the Debye law,
" Z #
3 ∞
T x3 dx θ 12π 4 3
CV ≃ 9N kB 4 x
− θ/T
= N kB (T /θ) . (4.37)
θ 0 e − 1 Te 5
n(ε) 1 −βmv2 /2
= e , (4.40)
N Ξcn
4.1.7 Exercises
4.1.7.1 Ex: Probabilities
In a game, 5 ideal dice are rolled.
a. What is the probability that exactly two of these dice show the number one?
b. What is the probability that at least one die shows the number one?
is a free spin.
a. Calculate the probability for one free spin per game.
b. What is the probability of getting exactly 3 free spins in 10 games?
c. What is the probability of winning at least once in 10 games?
This gives the average kinetic energy for each molecule of ⟨Ekin ⟩ = 21 kB T . According
to the equipartition theorem, Maxwell’s velocity distribution of a three-dimensional
gas is given by f (vx )dvx f (vy )dvy f (vz )dvz .
a. Write down the velocity distribution explicitly and determine the average kinetic
energy of a molecule in the three-dimensional gas at temperature T .
Determine the average absolute velocity ⟨v⟩ = ⟨|v|⟩ and compare ⟨v⟩2 with ⟨v 2 ⟩ for
the three-dimensional case.
c. What is the number of particles F (v)dv with an absolute velocity v = |v| in the
range v and v + dv.
134 CHAPTER 4. STATISTICAL THERMODYNAMICS
d. Consider a gas made of rubidium atoms (m = 87u) and sketch F (v) for tempera-
tures between 100 K and 300 K.
e. Consider the rubidium gas at room temperature (T = 300 K). What is the propor-
tion of molecules whose average velocity ⟨v⟩ is greater than 1000 m/s?
ε1 + ε4 = ε2 + ε3 . (4.47)
nj = f (εj ) . (4.48)
To satisfy Eqs. (4.46) and (4.47) f must have the functional form,
1
f (εj ) = , (4.49)
Ceβεj −s
Figure 4.3: (a) Detailed balance entails thermal equilibrium. (b) Subdivision of energy levels
j in subboxes gj . Red circles are fermions and green particles bosons.
4.2. QUANTUM STATISTICS 137
In analogy to the procedure outlined in Sec. 4.1.2 we derive the Boltzmann distri-
bution by first taking the logarithm from (4.51) and then simplifying it using Stirling’s
formula (3.1),
X X gj
ln W = ln N ! + [nj ln gj − ln nj !] ≃ ln N ! + [nj ln + nj ] , (4.51)
j j
n j
4.2.2.2 Bosons
Boltzmann’s fundamental equation (4.3) relates the thermodynamic entropy S to
the logarithm of the number of microstates W{nj } . It was pointed out by Gibbs
however, that the above expression (4.51) does not yield an extensive entropy, and is
therefore faulty 5 This problem is known as the Gibbs paradox. The problem is that
the particles considered by the above equation are not indistinguishable. In other
words, for two particles (i and j) in two energy sublevels the population represented
by [i, j] is considered distinct from the population [j, i], while for indistinguishable
particles, they are not. Indeed, bosons have anti-symmetric wavefunctions, Fermions
have symmetric ones. Boltzons have all wavefunctions as eigenfunctions. In the limit
of high temperatures all particles behave like boltzons. Discretize the one-particle
energies in small cells labeled j of constant energy εj . Let nj be their number and gj
their degeneracy [see Fig. 4.3(b)].
Figure 4.4: Distribution of nj bosons (green) and gfermions (red) over gj boxes the number
j −1
of possibilities being, respectively, nj +g
nj
and njj .
For bosons, each level gj can hold arbitrarily many of the nj particles. If we carry
out the argument for indistinguishable particles, we are led to the expression for the
partition function for bosons 6 ,
r
Y nj + gj − 1
W{nj } = . (4.58)
j=1
nj
can be seen as follows: Consider two identical systems, r′ = r and gj′ ′ = gj , with atom
5 This
numbers N = j nj and N ′ = j n′j . The partition function for boltzons is not multiplicative,
P P
r g j n r g j n r g j n′
j j j
Y Y Y
′ ′
(N + N )! ̸= N ! ×N !
j=1
nj ! j=1
nj ! j=1
n′j !
while for fermions it is. To see this we set n′j ′ ≡ nr+j for j ′ = r + j and j = 1, ..., r. Then,
2r r r
Y gj Y gj Y gj
= × .
j=1
nj j=1
nj j=1
n′j
The same argument holds for bosons. A critical discussion of the above statements can be read in
[33].
6 Note that this partition function converges toward the one for boltzons for g ≫ n ≫ 1, which
j j
can be seen by simplifying it using the Stirling formula.
4.2. QUANTUM STATISTICS 139
the differential,
X ∂ ln W X nj + gj − 1
d ln W = dnj = ln dnj , (4.60)
j
∂nj j
nj
with the same Lagrange multipliers. This yields the Bose-Einstein distribution,
gj − 1
nj = . (4.62)
eβ(εj −µ) − 1
The Boltzmann distribution follows from this Bose-Einstein distribution for temper-
atures well above absolute zero, implying that gj ≫ 1. The Boltzmann distribution
also requires low density, implying that gj ≫ nj . Under these conditions, we may use
Stirling’s approximation (3.1) for the factorial: N ! ≈ N N e−N .
4.2.2.3 Fermions
For fermions, each level gj can hold at most one of the nj particles, which implies
that necessarily gj > nj [see Fig. 4.3(b)]. Let us consider a single energy level j. The
first of the nj particles has the choice between gj boxes. Since no box can be filled
with more than one particle, the second particle has only gj − 1 boxes at its disposal,
and so on until all particles have been assigned. This corresponds to gj !/nj ! possible
choices. However, we still need to respect the indistinguishability requirement. The
overcounting can be removed by dividing by (gj −nj )!. The procedure is now repeated
with all energy levels j, which leads to the partition function for fermions,
r
Y gj
W{nj } = . (4.63)
j=1
n j
the differential,
X ∂ ln W X gj − nj
d ln W = dnj = ln dnj . (4.65)
j
∂nj j
nj
with the same Lagrange multipliers. This yields the Fermi-Dirac distribution for
gj ≫ 1,
gj
nj = β(ε −µ) . (4.67)
e j +1
Do the Exc. 4.2.6.1.
With this, knowing the energy spectrum εi and the distribution of states gj of
the system, we are able to calculate all thermodynamic potentials. E.g. the entropy
reads,
X
gj
nj
S = kB ln W{nj } = kB nj ln s + + sgj ln 1 + s (4.70)
j
nj gj
X sgj β(εj − µ)
β(εj −µ)
= kB + sgj ln e − s .
j
1 − se−β(εj −µ)
The Bose-Einstein and the Fermi-Dirac distribution both have many applications
in quantum mechanics, e.g. for the explanation of the blackbody radiation, the heat
capacity of metals, the laser, the Bose-Einstein condensation, and much more. In fact,
these distributions must be used whenever quantum statistical effects are important.
Prominent examples of systems where a quantum statistical treatment is crucial are
electrons in metals and ultracold quantum gases. We will discuss the latter in Secs. 4.3
and 4.4.
4.2. QUANTUM STATISTICS 141
(a) 2 (b) 2
kB T = 0.3μ E = 0.8μ
0.05μ μ
1.2μ
1.5 1.5
εj /μ
1 1
nj
0.5 0.5
0 0
0 1 2 0 1 2
nj kB T /μ
Figure 4.5: (code) Quantum statistical weight (4.49) for fermions (red dash-dotted line),
bosons (green dashed line), and boltzons (black solid line). (a) Weight nj as a function of
level energy εj for two different temperatures (solid and dash-dotted lines). (b) Weight nj
as a function of temperature for various level energy εj (solid, dash-dotted, and dashed).
ℏ2 2
Ĥ = − ∇ + U (r) . (4.71)
2m
As the wavefunction is localized, the spectrum of possible energies organizes into
discrete levels, and the atoms are allocated in populations of these levels. Such mul-
tidimensional systems are often degenerate, which means that the same total energy
can be realized with different sets of quantum numbers 7 . The way an atomic cloud
accommodates itself inside a trapping potential is governed by the density of available
states. We now introduce the density-of-states η(ϵ) for an arbitrary potential via,
(2m)3/2
Z Z Z Z
1 p
η(ε)dε ≡ d3 rd3 k = d3 r dε ε − U (r) , (4.72)
(2π)3 (2π)2 ℏ3
q
with the substitution k = 2m
ℏ2 [ε − U (r)].
As an example, let us consider a box potential of volume V . In this case, the
7 This can be checked easily with separable potentials, such as the rectangular 3D box potential or
the 3D harmonic oscillator, where the same energy E = Ex + Ey + Ez can be reached with different
combinations of Ex , Ey , and Ez .
142 CHAPTER 4. STATISTICAL THERMODYNAMICS
(2m)3/2 √ (2m)3/2 √
Z
η(ε) = d3 r ε = V ε (box potential) . (4.73)
(2π)2 ℏ3 V (2π)2 ℏ3
In the following we derive the density-of-states for the case of an harmonic oscillator
potential. More general potentials are discussed in the Excs. 4.2.6.2 and 4.2.6.3.
Figure 4.6: Artists’s view of phase space cells in a trapping potential in two dimensions.
where the coefficients nj with j = x, y, z are integer numbers. For the cylindrical
harmonic trap defined in (4.73), we find with a little help from Dr. Bronstein
[6],
(2m)3/2
Z r
3 m
η(ε) = d r ε − ωr2 ρ2 (4.78)
(2π)2 ℏ3 2
Z 1 √
Z 1−x̃2 √
Z 1−x̃2 −ỹ2
1 8ε2 p
= 2 3
dx̃ √ dỹ √ dz̃ 1 − x̃2 − ỹ 2 − z̃ 2 .
(2π) (ℏω̄) −1 − 1−x̃2 − 1−x̃2 −ỹ 2
ε2
η(ε) = (harmonic potential) . (4.79)
2(ℏω̄)3
4.2. QUANTUM STATISTICS 143
(a) (b)
100 100
(MHz)
50 50
0 2 0 2
U
2 0 1 0
0 -2 -2 0 -1 -2 x (mm)
y (mm) x (mm)
y (mm)
(c) 40 (d) 2
(MHz)
(108 )
30
20 1
η(ε)h̄ω
10
U
0
-1 0 1 0 50 100
x (mm) ε/kB (μK)
Figure 4.7: (code) (a) The figure shows two dimensions of a Ioffe-Pritchard type magnetic
trapping potential (characterized by being approximately linear at large distances from the
center and harmonic near the center). (b) Harmonic approximation (most experimentally
feasible potentials are approximately harmonic near the center). (c) One-dimensional cut
through the potential of (a,b). (d) Density-of-states for a harmonic (dotted line) and a
Ioffe-Pritchard type potential (solid line).
In the following sections we will calculate all system variables based on the expressions
(4.90) in the thermodynamic limit.
144 CHAPTER 4. STATISTICAL THERMODYNAMICS
p̂2m
ĥm = + Vtrap (r̂m ) , (4.83)
2m
with m = 1, . . . , N , whose spectrum is εm = ⟨ψm |ĥm |ψm ⟩. That is, every single
particle is completely characterized by the quantum number m 9 . A microstate |Ψk ⟩
is now identified as an eigenstate of the total many-particle Hamiltonian,
N
X N
Y
Ĥ = ĥm with |Ψk ⟩ = |ψm ⟩k . (4.84)
m=1 m=1
The request that the particles do not interact makes the system separable. The
density operator and the grand canonical partition function are [10],
e−β(Ĥ−µN̂ )
ρ̂ = and Ξgc = e−βΩ = Tr e−β(Ĥ−µN̂ ) , (4.85)
Ξgc
obviously satisfying Tr ρ̂ = 1. For the grand canonical ensemble the basis states of
the total Hamiltonian Ĥ are all microstates composed of many particles, and the
operators N̂ and ρ̂ can be expressed in the same basis.
8 In case of multi-species ensembles, the potentials add up like µ1 N1 + µ2 N2 + ....
9 In practice, a set of several quantum numbers may be required.
4.2. QUANTUM STATISTICS 145
We now migrate from the single-particle product state basis {|Ψk ⟩} to a Fock state
basis assigning a given number of particles nj to every possible energy level εj , where
j = 1, . . . , ∞, as illustrated in Fig. 4.8,
|Ψk ⟩ −→ |n1 , n2 , . . . , nj , . . .⟩ . (4.86)
I.e. we replace the distribution of microstates by a distribution of populations {nj }
among the energy levels. Since the energy and particle numbers are separately con-
served, the corresponding operators commute,
[Ĥ, N̂ ] = 0 , (4.87)
and therefore it is possible to find a complete basis of simultaneous eigenstates,
Ĥ| . . . nj . . .⟩ = E| . . . nj . . .⟩ with N̂ | . . . nj . . .⟩ = N | . . . nj . . .⟩ (4.88)
with,
∞
X ∞
X
E= εj nj and N= nj . (4.89)
j=0 j=0
This means that the number of particles is a conserved quantity and that Ĥ and N̂
can be simultaneously diagonalized.
We can now evaluate the partition function (4.85),
X
Ξgc = ⟨Ψk |e−β(Ĥ−µN̂ ) |Ψk ⟩ (4.90)
k∈{microstates}
X X
= ⟨. . . nj . . . |e−β(Ĥ−µN̂ ) | . . . nj . . .⟩ = e−β(E−µN ) .
{nj } {nj }
Figure 4.8: (a) Ensemble of N particles with different positions and velocities. (b) Distri-
bution of the particles over the spectrum of allowed energies.
e−β(Ĥ−µN̂ )
P
1 −β {nj } (εj −µ)n̂j
ρ̂ = = e (4.94)
Ξgc Ξgc
∞
1 −β(ε1 −µ)n̂1 −β(ε2 −µ)n̂2 Y
= e e × ... = ρ̂j .
Ξgc j=1
P
Note, that breaking down the exponential of a sum of operators, e− n̂j , into a prod-
Q −n̂
uct of exponentials of that operators, e j , is only possible because the operators
commute, [n̂k , n̂j ] = 0. In the last step we defined,
The problem with this expression is, that the global wavefunction |Ψ⟩ has not
yet been (anti-)symmetrized according the particles’ bosonic or fermionic nature.
For bosons, nj may be any non-negative integer and each value of nj counts as
one microstate due to the indistinguishability of particles. For fermions, the Pauli
exclusion principle allows only two microstates for the orbital (occupation of 0 or 1),
giving a two-term series 10 ,
P∞ −β(nj εj −nj µ) 1
nj =0 e =
1−e−β(εj −µ)
for bosons
Ξj = (4.96)
P1 −β(nj εj −nj µ) −β(εj −µ)
nj =0 e = 1 + e for fermions
Hence,
∞
Y
Ξgc = (1 − se−β(εj −µ) )−s , (4.97)
j=1
automatically introducing field operators satisfying bosonic or fermionic commutation rules. Indeed,
we can rewrite the Hamiltonian and the number operator of any non-interacting system like [49],
εj â†j âj
X X †
Ĥ = and N̂ = âj âj ,
{nj } {nj }
where â†j
and âj are the particle creation and annihilation operators introduced in the occupation
number representation.
4.2. QUANTUM STATISTICS 147
Considering again the entire system, the total Landau grand potential is found by
adding up the Ωj for all orbitals,
∞
X
Ω= Ωj . (4.99)
j=1
11
In any case the value
∂Ωj 1
nj = − = β(ε −µ) ≡ wT,µ (εj ) (4.100)
∂µ e j −s
gives the thermodynamic average number of particles on the orbital: the Fermi-Dirac
distribution for fermions, and the Bose-Einstein distribution for bosons.
The problem is completely analogous to Planck’s treatment of blackbody radia-
tion, where the Bose-Einstein distribution function followed as a corollary from the
Boltzmann statistics in thermal equilibrium and Planck’s quantization hypothesis,
E = N εj .
1 = Tr ρ̂
∂Ω
N = ⟨N̂ ⟩ = Tr ρ̂N̂ =−
∂µ
T,V
∂Ω
S = Tr ρ̂ ln ρ̂ = − . (4.101)
∂T
µ,V
∂Ω
P =−
∂V T,µ
E = ⟨Ĥ⟩ = Tr ρ̂Ĥ = T S + µN + Ω
11 Note the absence of the degeneracy factor g in comparison to the formula (4.68), which is
j
simply due to the fact that here we only consider a potential with non-degenerate eigenstates. The
degeneracy factor gj can, however, simply added ad hoc.
148 CHAPTER 4. STATISTICAL THERMODYNAMICS
∂Ξgc
= βTr N̂ e−β(Ĥ−µN̂ ) = βΞgc ⟨N̂ ⟩ (4.103)
∂µ
∂Ω Tr N̂ e−β(Ĥ−µN̂ )
−β = = β⟨N̂ ⟩
∂µ Ξgc
∂2Ω
∂
Ξgc ∂µ Tr N̂ e−β(Ĥ−µN̂ ) − Tr N̂ e−β(Ĥ−µN̂ ) ∂µ
∂
Ξgc
−β 2
= 2
= β(⟨N̂ 2 ⟩ − ⟨N̂ ⟩2 ) .
∂µ Ξgc
The key behind second quantization is to remove the restriction that the number of
particles is fixed. Instead, the theory is built around the idea of Fock space, where the
number of particles is not fixed. This is highly advantageous when dealing with many-
body systems. This same idea, when extended to finite temperatures, is what we call
the grand canonical ensemble. What we want is to consider some finite temperature
density matrix ρ̂ ∼ e−β Ĥ , where the number of particles is not fixed, but can fluctuate
[39].
However, we cannot let it fluctuate arbitrarily since that would make no physical
sense. Instead, the basic idea of the grand canonical ensemble is to impose that the
number of particles in the system is only fixed on average. That is, we impose that,
⟨N̂ ⟩ = N . (4.104)
In some systems, the number of particles does indeed fluctuate. This happens, for
instance, in chemical solutions: if we look at a certain region of a liquid, the number
of molecules there is constantly fluctuating due to molecules moving in and out from
other regions. Of course, in many other systems, the number of particles is fixed.
However, even in these cases, pretending it can fluctuate may still give good answers
for large N (thermodynamic limit). The reason is that, as we have seen above, the
variance of N̂ scales as,
√
∆N̂ ∝ N , (4.105)
which is small. Hence, when N is large, the grand canonical ensemble will give
accurate answers, even if the number of particles is not actually allowed to fluctuate.
This is the idea behind ensemble equivalence: we are allowed to use an ensemble
where the number of particles fluctuates, even though it actually doesn’t, because in
the thermodynamic limit the fluctuations are small.
sector,
E(N1 , 1)
E(N1 , 2)
..
.
Ĥ =
E(N2 , 1) .
(4.106)
E(N2 , 2)
..
.
..
.
Suppose now that the system is in thermal equilibrium with exactly N particles,
which corresponds to a canonical ensemble. As resumed in Tab. 1.4, the conditions
for equilibrium are then obtained minimizing the Helmholtz free energy, dF = 0,
and the corresponding canonical density operator and partition function are (see also
Sec. 1.4.4),
e−β Ĥ X
ρ̂cn = , Ξcn (N ) = e−βE(N,m) , F = −kB T ln Ξcn (N ) . (4.107)
Ξcn (N ) m∈sector
This is a constrained sum, since we are only summing over that sector that has
exactly N particles. This constraint makes it notoriously difficult to compute the
sum in practice solving a Schrödinger equation with Ĥ.
Instead, in the grand canonical ensemble we allow the number of particles to
fluctuate but only fix them on average (4.104). To accomplish this we had to intro-
duce a new parameter µ, called the chemical potential, so that the grand canonical
equilibrium state is transformed to (see also Sec. 1.4.5),
e−β(Ĥ−µN̂ )
ρgc = , Ξgc = Tr e−β(Ĥ−µN̂ ) , Ω = −kB T ln Ξgc . (4.108)
Ξgc
Apparently, the chemical potential enters by shifting the Hamiltonian,
Ĥ → Ĥ − µN̂ . (4.109)
As resumed in Tab. 1.4, in grand canonical ensembles the conditions for equilib-
rium are obtained minimizing the Landau energy, dΩ = d(F − µN ) = 0. To obtain
the energy spectrum in the case of fluctuating particle numbers, we need to solve
a many-body Schrödinger equation (such as the Gross-Pitaevski equation) with the
Hamiltonian substituted by Ω̂ = Ĥ − µN̂ [17].
The logic behind µ is twofold. When the number of particles is allowed to fluctuate,
the value of µ is fixed externally (like the temperature). As a consequence the number
of particles ⟨N̂ ⟩ = N (µ, T ) is interpreted as a function of µ and T . Conversely, if the
number of particles N is fixed, then µ = µ(N, T ) is to be interpreted as a function of
N and T , which is to be determined as the solution of the implicit equation,
Tr N̂ e−β(Ĥ−µN̂ )
⟨N̂ ⟩ = =N . (4.110)
Tr e−β(Ĥ−µN̂ )
• In high energy particle physics, ordinary particles can be spawned out of pure
energy, if a corresponding antiparticle is created. Then, neither the number of
particles nor antiparticles are conserved, only their difference.
• The grand canonical ensemble is particularly useful for developing the thermo-
dynamics of large ideal trapped quantum gases. While the phenomenon of BEC
can be derived in any ensemble (in Sec. 4.2.2 we derived the bosonic partition
function from the detailed balanced assumption using combinatorial arguments),
when the dynamics of a condensate is the subject under study, it is often useful
to consider it as a separate system being in thermal and chemical equilibrium
with a reservoir. The role of a reservoir is played by the thermal cloud, which
always coexists with the condensate and which exchanges particles and energy
with it.
4.2.4.3 Fluctuations
Fluctuations if the system can be readily calculated as well [39],
∂⟨Ĥ⟩ ∂⟨Ĥ⟩
(∆Ĥ)2 = ⟨Ĥ 2 ⟩ − ⟨Ĥ⟩2 = kB T 2 + kB T µ (4.111)
∂T ∂µ
∂⟨N̂ ⟩
(∆N̂ )2 = ⟨N̂ 2 ⟩ − ⟨N̂ ⟩2 = kB T .
∂µ
12 Note that photons in a highly reflective cavity can be conserved and caused to have a non-zero
chemical potential µ.
4.2. QUANTUM STATISTICS 151
∂⟨N̂2 ⟩ ∂⟨N̂1 ⟩
⟨N1 N2 ⟩ − ⟨N1 ⟩⟨N2 ⟩ = kB T = kB T (4.112)
∂µ1 ∂µ2
∂⟨Ĥ⟩
⟨N̂1 Ĥ⟩ − ⟨N̂1 ⟩⟨Ĥ⟩ = kB T .
∂µ1
From the above expressions, it can be seen that the function Ω has the exact
differential,
dΩ = −S dT − ⟨N̂ ⟩dµ − P dV . (4.113)
Substituting the relationship (4.101)(v) for E into the exact differential of Ω, an equa-
tion similar to the first law of thermodynamics is found, except that some quantities
only appear as averages,
which, for confined ensembles, can even be simplified using the concept of density-of-
states η(ε). As long as we are deep in the thermodynamic limit, N → ∞, we expect to
obtain reliable results. Let us now do this exercise for an ideal quantum gas confined
in a box potential of volume V , whose density-of-states is given by (4.73).
We begin with the request that the chemical potential satisfies the normalization
condition, √ √
V 2m ∞
Z Z
εdε
N = wT,µ (ε)η(ε)dε = . (4.116)
(2π)2 ℏ3 0 eβ(ε−µ) ∓ 1
Introducing the thermal de Broglie wavelength,
s
2πℏ2
λth ≡ , (4.117)
mkB T
At this point, to simplify the notation, we introduce the Bose function and its
integral representation,
∞ Z ∞
X Zt 1 xξ−1 dx
gξ+ (Z) = = ≡ gξ (Z) , (4.120)
t=1
tξ Γ(ξ) 0 Z −1 ex − 1
where Γ(η) denotes the Gamma function. Analogically, we can define the Fermi
function via 13 ,
∞ Z ∞
−
X (−Z)t 1 xξ−1 dx
gξ (Z) = − ξ = ≡ fξ (Z) . (4.121)
t=1
t Γ(ξ) 0 Z −1 ex + 1
(a) (b) 6
cη
g3 5
2
g3/2
f3 4
gξ (Z)
ζ(ξ)
f3/2
1 3
0 1
0 1 2 1 2 3 4
Z ξ
Figure 4.9: (code) (a) Bose and Fermi functions for box potentials (g3/2 and f3/2 ) and for
harmonic potentials (g3 and f3 ). Also shown is the Boltzmann limit (4.122). (b) Riemann
function.
With all these definitions we can now rewrite the expression (4.119),
V (s)
N= g (Z) , (4.124)
λ3th 3/2
13 When the context is clear, we will use the shorter notations g and f for Bose and Fermi
ξ ξ
functions, respectively.
4.2. QUANTUM STATISTICS 153
4.2.6 Exercises
4.2.6.1 Ex: Quantum statistics
n particles are distributed over g > n different cells with the same probability. Cal-
culate the probabilities
a. that there is exactly one particle in each one of the first n cells;
b. that there is no cell with more than one particle.
Use the three different assumptions that:
i. the particles are boltzon, i.e. they are identifiable and arbitrarily many particles
can be assigned to each cell;
ii. the particles are bosons, i.e. they are NOT identifiable and arbitrarily many par-
ticles can be assigned to each cell;
iii. the particles are fermions, i.e. they are NOT identifiable and only a single particle
may be assigned to each cell.
V (2m3 )1/2 √
ρ(ε)dε = εdε ,
π 2 ℏ3
1
n(ε)ρ(ε)dε = (ε−ε )/k T ρ(ε)dε .
e F B +1
Calculate the maximum energy at T = 0.
1 N →∞
Z= −→ 1 , (4.127)
1 + 1/N
we find that the fugacity approaches unity. Thus, Z = 1 is the condition for a
macroscopic ground state population.
Let us now calculate the ground state population at finite temperatures. For
a free gas with energy spectrum, ε = p2 /2m, we derived the density-of-states η(ε)
in (4.73) 15 . Using the occupation number wT,µ (ε) for the Bose-Einstein distribu-
tion (4.100) in the thermodynamic limit, we express the total number of atoms as we
already did in Eq. (4.124),
Z ∞
V
N= wT,µ (ε)η(ε)dε = g3/2 (Z) . (4.128)
0 λ3th
14 See script on Quantum mechanics (2023).
15 We must, however, keep in mind that the state density approach is an approximation not valid
for experiments with a limited number of atoms.
4.3. CONDENSATION OF AN IDEAL BOSE GAS 155
The problem with the expression (4.128) now is, that the thermal de Broglie wave-
length diverges for T → 0, while the phase space density g3/2 (Z) is bounded be-
tween g3/2 (0) = 0 and g3/2 (1) ≈ 2.612, which we realize after a quick inspection of
Fig. 4.9(a). Hence, according to this formula, even taking the largest possible value
T →0
of the fugacity, Z −→ 1, the number of atoms in the lowest energy state tends to 0,
3/2
V mkB T T →0
N = 3 g3/2 (Z) < V g3/2 (1) −→ 0 . (4.129)
λth 2πℏ2
This is obviously in contrast to the expectation of a large ground state population for
T → 0.
The reason is, that in the process of converting the sum to an integral (4.115),
the density-of-states disappears as we approach the ground state, thus removing the
ground state from the spectrum of energies that can be occupied. Einstein’s idea to
resolve the problem, was to explicitly maintain a discrete term accounting for the
ground state population Nc and to add it to the expression (4.128),
V
N = Nc + g3/2 (Z) . (4.130)
λ3th
V V
N= g3/2 (Z) = 3 g3/2 (1) for T ≥ Tc , (4.131)
λ3th λc
with g3/2 (1) = 2.612. The first part of Eq. (4.130) holds for T ≥ Tc and provides a
mean of determining Z from temperature and total atom number. The second part
of Eq. (4.130) holds at T = Tc . Resolving it by Tc we obtain,
2/3
2πℏ2
N
kB Tc = . (4.132)
m V g3/2 (1)
Below the critical temperature we need to add an additional term Nc . Resolving the
full expression (4.130) by the fraction Nc /N of atoms condensed in the ground state
and substituting N from (4.131), we obtain,
λ3
(
Nc V λ3c g3/2 (Z) 1 − λ3c for T ≤ Tc
=1− g3/2 (Z) = 1 − 3 = th (4.133)
N N λ3th λth g3/2 (1) 0 for T ≥ Tc
156 CHAPTER 4. STATISTICAL THERMODYNAMICS
16
The superscript (3/2) denotes the box potential shape of the trapping potential .
In summary we have,
!3/2 2/3
(3/2)
2πℏ2
Nc min(T, Tc ) N
=1− (3/2)
with kB Tc(3/2) = .
N Tc m V g3/2 (1)
(4.134)
The abrupt occurrence of a finite occupation in a single quantum state at temperature
(3/2)
below Tc indicates a spontaneous change in the system and a thermodynamic
phase transition. Solve Exc. 4.3.4.1.
1
Nc /N
0.5
0
0 0.5 1 1.5
T /Tc
Figure 4.10: (code) Condensed fraction for an ideal Bose gas as a function of reduced
temperature for a (blue) in a box potential and (green) in a harmonic trap. Red circles
denote experimentally measured data points [28]. The red dashed line is a fit to the data.
The cyan dash-dotted line is a theoretical curve taking into account finite size effects and
interatomic interactions.
e−β(ĤN −µN̂ ) X
and ρ̂ ≡ and Tr . . . ≡ ⟨ψj | . . . |ψj ⟩ . (4.135)
Ξgc j
The parameters µ, V, T are held fixed. As we have seen, for large systems in the
thermodynamic limit, the sum can be replaced by an integral, which, in turn, may
be expressed by the Riemann zeta-function (see Secs. 4.2.5 and 3.1.2). The thermo-
dynamic potentials and their expressions are summarized in the following table 17 .
Table 4.1: Thermodynamic potentials for an ideal Bose gas (upper signs) or Fermi
gas (lower signs) trapped in a box potential.
P∞
Tr ρ̂ ln Ξgc j=1 limN →∞
−βεj
ln(1∓Ze )
− β1 ln Ξgc
P
Ω j ±β
1
µ β
ln Z
1 Tr ρ̂
∂Ωj
nj − ∂µ
wT,µ − β1 ∂
∂εj
ln Ξgc βεj
1
T,V e /Z∓1
n o
±
− ∂Ω ∂ V 1
P
N ∂µ
Tr N̂ ρ̂ Z ∂Z ln Ξgc j nj λ3
g3/2 + 1/Z−1
T,V th
βε
∂Ω
P nj e j 5V ±
S/kB − kB ∂T Tr ρ̂ ln ρ̂ ln Ξgc ± j ln Z 2λ3
g5/2 −{ln Z}
µ,V th
∂Ω 1 V ±
P − ∂V T,µ βV
ln Ξgc λ3
g5/2 −{N ln(1 − Z)}
th
∂ 3kB T V ±
≃ 3P V
P
E T S + µN + Ω Tr Ĥ ρ̂ − ∂β ln Ξgc j nj εj 2λ3
g5/2
th 2±
± 9N g3/2
∂E 15V
CV ∂T N,V 4λ3
g5/2 − ±
th 4g1/2
In the thermal Bose-gas phase, T ≥ Tc , we get from (4.137) the state equation,
+
PV g5/2 (Z) T →∞
= + −→ 1 . (4.138)
N kB T g3/2 (Z)
158 CHAPTER 4. STATISTICAL THERMODYNAMICS
In the classical limit, obtained by noticing gξ0 (Z) = Z, follows the well-known classical
ideal gas equation. In the Bose-condensate phase, T ≤ Tc , using the definition of the
critical temperature, we recover from (4.137) the equation of state (4.134).
2
4 (f)
optical density
(b) (d)
1
1 2
0 0 0
0 120 240 360 0 120 240 360 0 120 240 360
position (μm) position (μm) position (μm)
Figure 4.11: Ultracold 87 Rb gas at various temperatures (a,b) T > Tc , (c,d) T ≃ Tc , and
(e,f) T < Tc measured in experiment [28]. The figures (a,c,e) are two-dimensional false color
images of the momentum distribution. The figures (b,d,f) are cuts through the images.
In the same way as for a potential well we find for a harmonic potential,
3 3
kB T T
Nth = g3 (1) = N (3)
, (4.140)
ℏω̄ Tc
with g3 (1) = 1.202. Since Nc + Nth = N , the number of particles in the ground state
is,
!3 1/3
(3)
Nc min(T, Tc ) N
=1− (3)
with kB Tc(3) = ℏω̄ . (4.141)
N Tc g3 (1)
(a) (b) 2
(μK)
4 0
(105 )
-2
μ/kB
2
Nc
-4
0
0 0.5 1 1.5 0 0.5 1 1.5
T /Tc T /Tc
(c) 10 (d)
10
(μK)
C/NkB
5
5
E/kB
0 0
0 0.5 1 1.5 0 0.5 1 1.5
T /Tc T /Tc
Fig. 4.10 traces the condensed fraction Nc /N measured as a function of the reduced
(3)
temperature T /Tc . Experiments [28, 19] confirm Bose’s ideal gas theory in the
thermodynamic limit. A particularity of inhomogeneous trapping potentials is, that
the condensed and the normal phase separate in momentum space, simply because the
condensed atoms occupy only the ground state, whose spatial extend is small, while
thermal atoms are distributed over all energy levels. Fig. 4.11 shows a measurement
of velocity distributions of a cloud of atoms close to the critical temperature.
We note that smaller trapping volumes (or tighter potentials) increase the critical
temperature Tc , thus allowing for quantum degeneracy at higher temperatures, which
can be advantageous in experimentation. Also, at a given temperature, a strongly
confining potential reduces the total minimum number of atoms required to reach
condensation.
160 CHAPTER 4. STATISTICAL THERMODYNAMICS
2.0
1.5
E / NkBTo
0.0
1.0
-0.2
∆
0.5 -0.4
0.0
0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8
T / To(N)
The heat capacity quantifies the system’s ability to secure its energy. In conven-
tional systems, the heat capacity is typically either specified at constant volume or
at constant pressure. With this specification heat capacities are extensive state vari-
ables. When crossing a phase transition, the temperature-dependent heat capacity
measures the degree of change in the system above and below the critical temperature
and provides valuable information about the general type of phase transition.
Using (4.81), the total energy per particle is given by,
For a confined gas, volume and temperature are interdependent, and the concept of
pressure is somewhat vague. In this case, we can not refer to the heat capacity at
constant volume or pressure. However, one can define the heat capacity for a fixed
number of particles,
∂E(T )
C(T ) = . (4.143)
∂T N
Fig. 4.12 shows the temperature dependence of some thermodynamic potentials for a
harmonically trapped ultracold Bose gas. The discontinuity of the heat capacity at
the critical temperature is known as λ-point.
Calculating the second moments of the distributions obtained for the same density
by time-of-flight of absorption images, we obtain the kinetic energy,
Z 2
p
Ekin = n(p)d3 p . (4.144)
2m
4.3. CONDENSATION OF AN IDEAL BOSE GAS 161
For confined ideal gases, the virial theorem ensures Ekin + Epot = 2Ekin . For real
gases, the repulsive energy of the mean field adds to this energy, E = Ekin + Epot +
Eself . The sudden extinction of the trapping potential before time-of-flight takes away
the potential energy Epot non-adiabatically. The kinetic energy and the self-energy of
the condensate are fully converted into kinetic energy during ballistic expansion. It
is this energy, p2 /2m = Ekin + Eself , which is sometimes called release energy, which
is measured after ballistic expansion 18 . Fig. 4.13(right) shows a measurement of the
release energy. Solve the Exc. 4.3.4.3.
The question which ensemble is the correct assumption depends on the experimental
situation. The question is particularly interesting in the context of Bose-Einstein
condensation: Here it is related to the question which state better describes a BEC:
A Fock state characterized by a fixed atom number or a Glauber state, where the
atom number is fluctuating.
The condensates experimentally produced in alkali gases consisted of relatively
small atom numbers between 1000 to 107 , so that the validity of the thermodynamic
approximation and the use of the density-of-states approach has been questioned [25].
Also, the decision whether to use the grand canonical, the canonical or the micro-
canonical ensemble for calculating the thermodynamic quantities noticeably influences
the results. Herzog and Olshanii [29] have shown that for small atom numbers on
the order of 100 the canonical and grand canonical statistics lead to predictions on
the condensed fraction that differ by up to 10% [see Fig. 4.10(right)]. On the other
hand, they give the same results if the particle numbers are large. Which canonical
statistics is more appropriate is not a trivial question and depends on the experimen-
tal setup and in particular on the time scale of the measurements. If we look at the
sample for short times, the number of condensed atoms will be fixed, and we can
assume a canonical ensemble. For longer times, however, the atom number may be
an equilibrium parameter depending on the contact of the sample with a reservoir,
and the grand canonical statistics is better suited.
but we note that the calculations are analogous to the calculations for fermionic gases
presented in Sec. 4.4.4.
For an ideal Bose gas the density and momentum distributions are expressed by
Bose functions g3/2 (Z) [9]. For example, as will be derived in Exc. 4.3.4.4, the density
and momentum distributions are,
1
n(x) = g3/2 (e−β[U (x)−µ] )
λ3th
(bosonic distribution functions) (4.145)
a6ho 2
n(k) = g3/2 (eβ(µ−p /2m) )
λ3th
In the classical limit, we can calibrate the chemical potential by Eq. (4.136) for a
box potential or by (4.139) for a harmonic potential,
( N 3
βµ βµ βµ V λth (for a box potential)
g3/2 (e ) → c3/2 (e )=e = 3
c3 (eβµ ) =N ℏω̄
kB T (for a harmonic potential)
(4.146)
Hence, we obtain for the classical density distribution,
1 eβµ
n(x) = 3 c3/2 (e−β[U (x)−µ] ) = 3 e−βU (x) (4.147)
λ λth
th
N
V
x∈V
(for a box potential)
= q
mω̄ 2
3 2 2
−βmω̄ x /2
N
2πkB T e (for a harmonic potential)
where
p we used the spatial extend of the ground state of the harmonic oscillator aho =
ℏ/mω. We see that we recover the Maxwell-Boltzmann velocity distribution, as
seen in Fig. 4.14,
3
m3
r
m 2
n(v) = n(k) =N e−βmv /2 . (4.149)
ℏ3 2πkB T
4.3. CONDENSATION OF AN IDEAL BOSE GAS 163
(a) 1 (b) 3
0 0
-50 0 50 -50 0 50
x (aho ) p (h̄/aho )
Figure 4.14: (code) (a) Density and (b) momentum distribution of a Bose gas (red) and a
Boltzmann gas (green) at T = 1.1Tc (solid line) and at T = 2Tc (dotted line).
(4.150)
where we defined, r
kB T
rrms ≡ tT oF . (4.151)
m
This distribution does not directly depend on the potential U (r), that is, the expansion
is isotropic. In Exc. 4.3.4.4(b) we determine the time-of-flight density distribution of
an ultracold Bose gas. For very long flight times (usually several 10 ms) the density
resembles a Gaussian distribution [9]. Note however, that in interacting non-ideal
gases the chemical potential does depend on the potential.
In a time-of-flight experiment, any deviation observed between the results (4.150)
and (4.151) points towards an impact of quantum statistics. However, absorption
images only record column densities, i.e. projections of the time-of-flight distribution
on a plane, which tends to smear out the non-Gaussian features.
Example 33 (Heat capacity measurement): For an ideal Bose gas trapped
in a harmonic potential the temperature dependence of the heat capacity at
the threshold to condensation can easily be obtained as follows. The condensed
fraction determines the chemical potential through,
3
kB T
N = N0 + g3 (Z) , (4.152)
ℏω
where Z(T ) = eµ/kB T for a grand canonical ensemble and gn denotes the
Riemann zeta function. The condensed fraction vanishes above the critical
164 CHAPTER 4. STATISTICAL THERMODYNAMICS
For an interacting Bose-gas we expect that the Eqs. (4.152) and (4.153) are
not scrupulously obeyed. Indeed, the abrupt discontinuous change in the heat
capacity at the phase transition to BEC, expected for ideal gases, is smeared
out by atomic collisions [1].
Figure 4.15: Population variation during a slow adiabatic compression followed by a sudden
non-adiabatic decompression.
4.3.4 Exercises
4.3.4.1 Ex: Monoatomic gas as a canonical ensemble
Consider a classical monoatomic gas made up of N non-interacting atoms of mass m
confined in a container of volume V , at temperature T . The Hamiltonian correspond-
ing to an atom is given by Ĥ = (p̂2x + p̂2y + p̂2z )/2m.
a. Show that the atomic canonical partition function is ξ = V /λ3th , where λth is the
thermal de Broglie wavelength defined in Eq. (4.113).
b. Using ξ of the previous item, obtain the system’s partition function Ξcn and the
Helmholtz free energy F . Also obtain the free energy per atom f = F/N in the
thermodynamic limit N −→ ∞, V −→ ∞, such that v = N/V fixed.
c. Obtain internal energy E and the gas pressure p.
d. Calculate the chemical potential and entropy per atom in the thermodynamic limit,
thus deriving the so-called Sackur-Tetrode formula.
4.3. CONDENSATION OF AN IDEAL BOSE GAS 165
(bosonic potentials)
ξ
Nc min(T, Tc )
= 1−
N T
c ξ
E gξ+1 (Z) min(T, Tc )
= ξ
N kB T gξ (Z) Tc
S gξ+1 (Z) 2µ
= 4 −
N kB gξ (Z) kB T .
ξ
C gξ+1 (Z) min(T, Tc ) gξ (Z) max(T − Tc , 0)
= ξ(ξ + 1) − ξ2
N kB gξ (Z) Tc gη−1 (Z) T − Tc
CT >Tc gξ+1 (Z) 2 gξ (Z) CT <Tc gξ+1 (1)
= ξ(ξ + 1) −ξ , = ξ(η + 1)
N kB gξ (Z) gξ−1 (Z) N kB gξ (1)
∆CTc CTc− − CTc+ gξ (1)
= = ξ2
N kB N kB gξ−1 (1)
For low temperatures, βµ ≫ 1, we can use the Sommerfeld expansion of the Fermi
function, which in first order gives fξ (ex ) ≃ xξ /Γ(ξ + 1), where x is a placeholder for
βµ, Γ is the Γ-function, and ξ = 3 for a harmonic potential. From this we immediately
obtain the chemical potential at zero temperature defined as the Fermi energy,
and from that the momentum of free particles and the Fermi radius,
r s s
2mEF 2EF 2EF
KF ≡ and rF ≡ , zF = . (4.156)
ℏ2 mωr2 mωz2
19 We note that meanwhile ultracold two-components Fermi gas have been demonstrated to form
bosonic Cooper-pairs, similarly to the phenomena known as superconductivity in some metals and
as superfluidity of the fermionic 3 He.
4.4. QUANTUM DEGENERACY OF AN IDEAL FERMI GAS 167
xξ π 2 ξ(ξ − 1)
x
fξ (e ) ≃ 1+ + ... , (4.157)
Γ(ξ + 1) 6x2
For highly excited atoms, ε−µ ≫ kB T , the Fermi function approaches the identity,
Z→0
fξ (Z) −→ Z (see Fig. 4.9), so that,
3 3
kB T kB T
N= eβµ = (1 + βµ + ...) , (4.159)
ℏω̄ ℏω̄
" 3 # 3
ℏω̄ 1 EF
µ = kB T ln Z ≃ kB T ln N = kB T ln ,
kB T 6 kB T
where in the last step we substituted the definition of the Fermi energy. This means
that highly excited fermions behave like a Boltzmann gas, which satisfies an ideal gas
equation similar to that of classical particles in a box potential,
3
kB T
N= . (Boltzmann) . (4.160)
ℏω̄
Fig. 4.16(a) shows calculations of the chemical potential for an ideal Fermi gas
along with the chemical potentials of a Boltzmann gas and a Bose gas.
4.4.2 Energy
Using (4.81), the total energy per particle, E/N ≡ N −1 εwT µ d3 xd3 k, is given by,
R
−1
εη(ε) eβ(ε−µ) + 1
R
εwT,µ (x, k)d3 xd3 k
R
E dε f4 (Z)
= R 3 3
= R −1 = 3kB T , (4.161)
N wT,µ (x, k)d xd k η(ε) e β(ε−µ) +1 dε f3 (Z)
in analogy to the expression (4.142) holding for a Bose gas. Again using the Sommer-
feld approximation, we see that for low temperatures, T → 0, the energy is limited
by [see Fig. 4.16(b)],
3µ4 2π 2
3 βµ T →0 3
E= f 4 (e ) = 1+ + ... −→ EF . (Fermi) (4.162)
β(βℏω̄)3 4EF3 (βµ)2 4
Hence, the total energy per fermion does not vanish for T → 0. The reason is that
the atoms are forced to adopt states in the outermost regions of the harmonic trap.
168 CHAPTER 4. STATISTICAL THERMODYNAMICS
In the limit of high temperatures, T → ∞, a classical gas has the energy per
particle,
3 3
E= 3
f4 f3−1 (βE6F ) ≃ 3N kB T . (Boltzmann) (4.163)
β(βℏω̄)
Z→0
which is seen by taking the high temperature limit fη (Z) −→ Z and extrapolating
T →∞
to all Z. This implies, E1 /EF −→ 3kB T /EF .
Hence, the total energy per boson decreases very rapidly for T → 0. The reason
is that the atoms are bosonically encouraged to pile up in the inner region of the
harmonic trap.
and becomes,
f4 (Z) µ 4E1 µ
− =
S1 = 4kB − . (4.166)
f3 (Z) T 3T T
fη−1 (Z),
f4 (Z) 3µ f4 (Z)f2 (Z) E1 3µ f4 (Z)f2 (Z)
C1 = 3kB − 1− = − 1 − . (4.167)
f3 (Z) T f3 (Z)2 T T f3 (Z)2
For fermions well below the Fermi temperature, T → 0, using the Sommerfeld
approximation, we calculate,
T →0 3π 2 kB T
C1 −→ . (Fermi) (4.168)
2 TF
(a) 1 (b) 3
(μK)
(μK)
0 2
EF
E/kB
μ/kB -1 1
3
4 EF
-2 0
0 0.5 1 0 0.5 1
T /TF T /TF
×106
(c) (d)
10 2
Tc
C/NkB
C/kB
5 1
0 0
0 0.5 1 0 0.5 1
T /TF T /TF
Figure 4.16: (code) Calculation of thermodynamic potentials for Bose (red), Fermi (green),
and Boltzmann gases as a function of temperature for a given harmonic trapping potential.
The gases are assumed to have same mass, same atom number N = 200000, and same trap
frequencies ωho /2π = 200 Hz. (a) Chemical potential, (b) energy, (c) heat capacity per
particle, and (d) total heat capacity. The dotted magenta line in (a) shows the chemical
potential calculated from the Sommerfeld approximation.
4πρ2 dρ
Z
1
= 3 β[ε+mω 2 2
(2π) e r ρ /2−µ] + 1
3/2 Z √ 3/2
1 2 tdt 1 2
= = Γ(3/2)f3/2 (eβ(µ−ε) ) ,
(2π)2 βmωr2 eβ[ε+t−µ] + 1 (2π)2 βmωr2
such that,
ñ(k) = λ−3 6
th aho f3/2 (e
β(µ−ε)
) (Fermi) . (4.177)
At low temperatures, T → 0,
3/2
1 2 Γ(3/2)
ñ(k) ≈ 2 2
(β [µ − ε])3/2 (4.178)
(2π) βmωr Γ(5/2)
3/2 3/2 3/2
ℏ2 k 2 k2
1 2 2 8 N
≈ EF − = 2 3 1− 2 .
(2π)2 mωr2 3 2m π KF KF
This can easily be integrated by dimensions,
Z ∞Z ∞ 3/2
k2
Z Z
8 N
ñT →0 (kz ) = ñcl (k)dkx dky = 2 3 1− 2 dkx dky
−∞ −∞ π KF |k|≤KF KF
(4.179)
√ 2 2 !3/2
KF −kz 5/2
2π Z
kz2 kρ2 kz2
Z
8 N 16 N
= 1− − 2 kρ dkρ dϕ = 1− 2 .
π 2 KF3 0 0 KF2 KF 5π KF KF
4.4. QUANTUM DEGENERACY OF AN IDEAL FERMI GAS 171
R∞
It is easy to check −∞ ñT →0 dkz = N , with Maple.
At high temperatures, T → ∞, we should recover the Boltzmann gas situation,
3/2
ℏ2 ω̄ 2
ñ(k) ≈ N e−βε (Boltzmann) . (4.180)
2πmωr2
√
Since ε is the kinetic energy, the rms-radius k 2 of this distribution is βℏ2 ⟨k 2 ⟩ = m.
In comparison,
h i
β (µ−p2 /2m)
ñ(k) = λ−3 6
th aho g3/2 e (Bose gas above Tc ) . (4.181)
1/3
where ω̄ ≡ (ωx ωy ωz )1/3 and ω̌ ≡ (ω̌x ω̌y ω̌z ) .
1 ω̌ 3
βµ− 12 βmΣj ω̌j2 x2j
nT oF (x, t) = f 3/2 e . (4.186)
λ3th ω̄ 3
×1010
5
4
(cm−3 )
3
nT oF
0
0 100 200 300
r (μm)
At low temperatures,
3/2
m 3 N 8 (mx/ℏt)2
nT oF (x, t) = 1− (4.188)
ℏt KF3 π 2 KF2
!2 3/2
m 3 R 3 R mx/ℏt
F F
= 1−
ℏt 6π 2 λ (48N λ)
1/3
4.4. QUANTUM DEGENERACY OF AN IDEAL FERMI GAS 173
At high temperatures,
1 1 2 2
nT oF (x, t) = f3/2 (eβ(µ−mx /2t ) ) (4.189)
λ3th ω̄ 2 t2
1 1 2 2
≈ 3 2 2 eβ(µ−mx /2t )
λth ω̄ t
3/2 3 3/2
mkB T 1 ℏω̄ −βmx2 /2t2 ω̄ m 2 2
≈ N e ≈ N e−βmx /2t .
2πℏ2 ω̄ 2 t2 kB T t2 2πkB T
A rms-width is,
Z
rT2 oF = r2 nT oF (x, t)d3 x (4.190)
1 ω̌ 3
Z
1 2 2
= 3 3 r2 f3/2 eβµ− 2 βmΣj ω̌j xj d3 x
λth ω̄
Z
2 εg(ε)dε kB T g4 (Z)
= = .
mω̌r2 N eβ(ε−µ) + 1 mω̌r2 g3 (Z)
This shows that the width of the flight-of-time distribution
√ can simply be obtained
from the spatial distribution by substituting ω → ω/ 1 + ω 2 t2 . Of course this does
not hold for condensed gases Bose.
This confirms the equipartition theorem for confined particles, which postulates,
Or we determine the temperature Tg where the Boltzmann gas turns into a Fermi gas
3µ/4 = 3kB Tg ,
3
32 kB Tg
N= . (4.194)
3 ℏω̄
174 CHAPTER 4. STATISTICAL THERMODYNAMICS
(2m)3/2
Z p
η(ε) = ε − U (r, z)rdrdz . (4.195)
2πℏ3
⟨k 2 ⟩
Z
E1 1 εη(ε)dε
2 = = β(ε−µ)
. (4.197)
kF EF N EF e +1
⟨k 2 ⟩ 3T f4 (Z)
2 = , (4.199)
kF TF f3 (Z)
⟨k 2 ⟩ n+1
2 = . (4.200)
kF βEF
Must for a single dimension the value be divided by three? ℏ2 ⟨kj2 ⟩ = 2mkB T f4 (Z)/f3 (Z)
setting ε = ℏ2 k 2 /m.
For a harmonic potential η(ε) ∝ ε2 and for a linear potential η(ε) ∝ ε7/2 . Inter-
mediate values are possible for non isotropic traps, which are linear in some directions
and harmonic in others, e.g. for a radially quadrupolar and axially harmonic trap, we
expect η(ε) ∝ ε3 and thus E = 4N kB T . In general, we may have more complicated
4.4. QUANTUM DEGENERACY OF AN IDEAL FERMI GAS 175
situations, where the trap becomes non-harmonic beyond a certain distance from the
origin. In those cases, the density-of-states may be approximated by series,
η(ε) ∝ ε2 + κε3 , (4.201)
where η is a small parameter, so that,
1 (ε3 + κε4 )(eβ(ε−µ) + 1)−1 dε
R
⟨k 2 ⟩ T 3f4 (Z) + 12κf5 (Z)
= = , (4.202)
kF2
R
EF (ε2 + κε3 )(eβ(ε−µ) + 1)−1 dε TF f3 (Z) + 3κf4 (Z)
which in the classical limit gives rise to energies E = 3..4N kB T depending on the
value of κ.
Such effects must be considered when the time-of-flight method is used for temper-
atures measurements. For example, if we underestimate η(ε) by assuming a harmonic
potential at all ε, although the potential is quadrupolar at large ε ≫ kB T , we get a
wrong estimate for the temperature Twrng = E/3N kB instead of Tcorr = E/4N kB .
3/2 3 s 3
ℏ2
1 ℏωho −βε 2 2
ñT →∞ (k) = 2 N e =N e−ℏ k /2mkB T .
2πβmωho kB T 2πmkB T
(4.207)
176 CHAPTER 4. STATISTICAL THERMODYNAMICS
(a) ×10−15
1 (b) 0.04 total atom number N = 20000
temperature T = 0.5 μ K
Fermi temperature T F = 2.446
0.03
n(k)
n(k)
0.5 0.02
0.01
0 0
0 1 2 -2 0 2
k/kF k/kF
Figure 4.18: (code) (a) Radial momentum distribution and (b) distribution of momentum
classes in the direction of kz for a Fermi gas at T /TF = 0.2 µK (red solid), a classical gas
(black), and a Fermi gas at T = 0 (red dash-dotted).
is independent on the atomic mass. This means that a Li and a Rb cloud in the same
harmonic trap at the same temperature have the same radius. This ensures good
overlap. E.g. at T = 10 µK assuming the Rb secular frequencies ωr ≃ 2π × 300 Hz
and ωz ≃ 2π × 30 Hz, we expect rrms = 16 µm and zrms = 160 µm. However below
the temperature 0.5TF , which is TF ≃ 1 µK for NF = 104 , the quantum pressure
stops the reduction of the fermion cloud while cooling. This evtl. reduces the overlap
with the boson cloud, disconnects the two clouds and stops the evaporative cooling.
On the other hand, the interaction energy of the boson cloud also increases its size,
when the Rb cloud approaches the critical temperature Tc ≃ 0.6 µK for NB = 106 .
The Pauli blocking of sympathetic cooling is a signature for the advent of quantum
statistics [15, 24, 45]. It is due to a reduced mobility (or better reduced available phase
space at collisions) of the atoms and not to be confused with the prohibition of s-
wave collisions due to the Pauli exlusion principle. Furthermore, elastic collisions are
suppressed [14], because atoms cannot be scattered into occupied trap levels [31, 48,
26, 27].
The fermions inside the bosonic cloud can be regarded as impurities. If they travel too
slow, v < c, and if the condensed fraction is too large, the motion will be frictionless
and thermalization stops. If they travel fast, quasiparticles are excited, which can
be√ removed by evaporation. With the typical velocity of sound in the BEC c =
ℏ 16πna/2mB ≈ 2 mm/s, or m 2
2 c ≈ kB × 20 nK, we see that this is no real danger.
178 CHAPTER 4. STATISTICAL THERMODYNAMICS
Γ ns + 1 Γ nt + 1
n/s+n/t
N =g 2 (kB T ) fn/s+n/t (z) . (4.215)
n n/s n/t n
(2ℏ) a b Γ 2 + 1
1 2
This gives for a harmonic trap where ε = 2m p +m 2 2
2 ω r and with the spin degeneracy
factor g = 1, n
kB T
N= fn (z) . (4.216)
ℏω
The Fermi energy again follows from Sommerfeld’s expansion,
In 1D,
1.02
N= (kB T )3/4 f3/4 (z)
ℏb
EF ≈ 0.87(N ℏb)4/3 .
vanishes, i.e. the electron gas does not contribute to the heat capacity of a metal. It
is only at very low temperatures of a few K, when the heat capacity of the atomic
lattice drops due to the underlying bosonic statistics, that the electrons contribute.
Now, make the metallic box potential having a finite depth. An electron can then
leave the metal, if it surmounts the exit work W = −Vmin − EF ≃ 10 eV, which is the
difference between the potential depth and the Fermi energy. At high temperatures,
the tail of the Fermi-Dirac distribution can leak into the unbound regime, which gives
rise to thermoionic emission. This feature explains the existence of contact potentials:
Metals with different W and EF brought into contact exchange charges until their
Fermi level is at same height.
4.4.8 Exercises
4.4.8.1 Ex: Li Fermi gas
Programs on Li Fermi gases.
Quantum theory of
non-relativistic scalar particle
fields
In the preceding sections we mostly restricted to single particles except in the con-
text of the electron shell of atoms, where we presented the Pauli principle in Sec. ??.
In the context of optical lattices [see Eq. (??)] we introduced particle creation and
annihilation operators, which we only used to simplify the notation without consid-
ering the possible presence of several atoms per lattice site. On the other hand, in
the context of the harmonic oscillator (Sec. ??), we allowed for large population of a
quantum states by quasi-particles identified as phonon or photons (see Sec. ??). The
commutation rule that the field operators of the quasi-particles had to obey was the
expression of their bosonic nature.
The question now arises how to treat the presence of several massive particles. In
quantum mechanics (differently from classical mechanics) there is no way of differen-
tiating identical particles, and this feature has far-reaching consequences in the way
of counting the number of possible microstates that a system may adopt. The field
of physics preoccupied with this counting is quantum statistics.
The quantum mechanical state of a system is described by a global wavefunc-
tion, which must somehow be composed from the wavefunctions of all the particles
of the system in such a way that the indistinguishability of particles is respected.
In mathematical terms this means that the global wavefunction must satisfy some
symmetry requirements upon particle exchange. We will show in Sec. 5.1.2 that the
system’s Hamiltonian must commute with the operator describing particle exchange,
and consequently the eigenstates of the Hamiltonian must also be eigenstates of the
particle exchange operator. We will find that only two eigenstates are possible: a
symmetric and an antisymmetric one. Which one of the symmetries the system must
exhibit depends on the nature of the particles: Particles behaving symmetrically under
particle exchange are called bosons, particles behaving antisymmetrically are called
fermions. Obviously, the symmetry requirement drastically limits the ways how the
single-particle wavefunctions may be joined to a global wavefunction. We will see that
bosons have increased probability to occupy the same quantum state, while fermions
cannot occupy the same quantum state at all.
The quantum statistical nature of particles has far-reaching consequences for the
description of many-body problem; when the particles do not interact with each other
181
182 CHAPTER 5. QUANTUM THEORY OF NON-RELATIVISTIC SCALAR PARTICLE FIELD
and even more when interactions between the atoms couple the single-atom states.
The formalism that rigorously enforces the exchange symmetry of the many-body
wavefunction is known as the second quantization method, which will be introduced
in Sec. 5.2, although a less confusing name would be occupation number (or Fock)
representation. In Sec. 5.2.2 we introduce construction operators for the creation or
annihilation of particles in properly symmetrized many-body states and show that the
proper quantum statistics is enforced by an algebra defined through the commutation
relations between the construction operators.
All this being said it is often not so obvious what is actually meant by particle.
In most cases, the particles under study are composed from other particles bound
together through strong forces. E.g. protons and neutrons bound together by nuclear
forces form a nucleus, electrons bound to a nucleus by Coulomb forces form an atom,
and atoms bound together by van der Waals forces form a molecule. The bosonic
or fermionic nature of composite particles somehow depends on the nature of the
components and on the way they interact. We consider a composite system as a
particle as long it has a well-defined internal structure. If all particles of a system
have the same internal structure they are called identical. When their structure brakes
down, for instance, induced by strong inter-particle forces, new quantum numbers
appear (related to the symmetries of the new structure), but the statistical nature
under exchange of complete particles (e.g. all components of an atom) is conserved
(see Sec. 5.1.2) [49].
ℏ2 2
− ∇ ψk (r) = εk ψk (r) , (5.1)
2m
ℏ2 k 2
ψk (r) = √1
V
eık·r and εk = . (5.2)
2m
The wavefunctions ψk (r) represent plane wave solutions, normalized to the volume of
the box, with k the wave vector of the atom, k = |k| = 2π/λdB the wave number,
and λdB the de Broglie wavelength. The periodic boundary conditions give rise to a
discrete set of wavenumbers, kα = (2π/L)nα with nα = 0, ±1, ±2, ... and α = x, y, z.
5.1. QUANTIZING SCALAR FIELDS 183
The Hamiltonian for the motion of two atoms with interatomic interaction V(r12 )
and confined by the cubic box potential U(r) defined above is given by,
X ℏ2
H= − ∇2i + U(ri ) + V(r12 ) . (5.3)
i=1,2
2mi
When the interatomic interaction can be neglected, the Schrödinger equation takes
the form 1 ,
ℏ2 2 ℏ2 2
− ∇1 − ∇2 Ψk1 ,k2 (r1 , r2 ) = Ek1 ,k2 Ψk1 ,k2 (r1 , r2 ) . (5.4)
2m1 2m2
In this limit, we have complete separation of variables so that the pair solution can
be written in the form of a product wavefunction
ℏ2 k12 ℏ2 k22
Ek1 ,k2 = + and P = ℏk1 + ℏk2 . (5.6)
2m1 2m2
Hence, the total momentum is not affected by particle exchange, but the energy,
Importantly, only for pairs of unlike atoms the product wavefunctions (5.5) represent
uniquely defined quantum mechanical eigenstates for the eigenvalues Ek1 ,k2 . By unlike
we mean that the atoms may be distinguished from each other because they are of
different species. For identical atoms, i.e. atoms of the same species such that m1 =
m2 , the situation is fundamentally different. In this case the product wavefunctions
(5.5) are degenerate with pair wavefunctions in which the atoms are exchanged,
This is called exchange degeneracy (two states, one energy). This phenomenon shows
that we cannot determine which atom is in which state by measuring the energy of
the pair. Actually, it is fundamentally impossible to distinguish the atoms by any
measurement because (by definition of being identical) the expectation values of all
observables of the pair are invariant under exchange of the atoms. For identical atoms
in the same state (k1 = k2 ) the exchange degeneracy is manifestly absent because there
is nothing to be distinguished.
The exchange degeneracy implies that any linear combination of the type,
1 1
Ψk1 ,k2 (r1 , r2 ) = p (c1 eık1 ·r1 eık2 ·r2 + c2 eık2 ·r2 eık1 ·r1 ) (5.9)
V |c1 | + |c2 |2
2
1 We use the letter Ψ to denote single-particle orbitals and the letter Ψ for many-particle orbitals.
184 CHAPTER 5. QUANTUM THEORY OF NON-RELATIVISTIC SCALAR PARTICLE FIELD
represents a properly normalized energy eigenstate of the pair. For example, we can
choose c1 = 1 and c2 = ±1, yielding,
1 1 ık1 ·r1 ık2 ·r2
Ψ±
k1 ,k2 (r1 , r2 ) = √ (e e ± eık2 ·r2 eık1 ·r1 ) . (5.10)
V 2
These are called the symmetric and antisymmetric eigenstates of the pair. In the
absence of exchange degeneracy (i.e. for k1 = k2 ) Eq. (5.5) represents the proper
solution. This state is symmetric and unit normalized to start with, so there is no
need for explicit symmetrization. Do the Exc. 5.1.3.1.
Ψ(r1 , σ1 ; r2 , σ2 ) , (5.11)
where r1 and r2 are the position coordinates and σ1 and σ2 the spin coordinates
of the particles 1 and 2, respectively. The squared modulus of the wavefunction,
|Ψ(r1 , σ1 ; r2 , σ2 )|2 , corresponds to the probability of observing particle 1 at position
r1 in spin state σ1 with particle 2 at position r2 in spin state σ2 . As before, two
particles are called identical, if there is no physical way to establish whether or not
they particles have been exchanged, a condition which is satisfied for particles with
identical internal structure.
To describe the particle exchange, we introduce the exchange operator P. For two
identical particles in an arbitrary pair state Ψ(r1 , σ1 ; r2 , σ2 ) the operator P is defined
by
PΨ(r1 , σ1 ; r2 , σ2 ) ≡ Ψ(r2 , σ2 ; r1 , σ1 ) , (5.12)
The effect of this operator is to exchange the particle labels. As tracking is impossible
in quantum mechanics, we have no physical means to determine whether or not two
identical particles have been exchanged, i.e. |Ψ(r1 , σ1 ; r2 , σ2 )|2 = |Ψ(r2 , σ2 ; r1 , σ1 )|2 .
Because P is a norm-conserving operator, ⟨Ψ|P † P|Ψ⟩ = 1, we have
P †P = 1 . (5.13)
Furthermore, exchanging the particles twice must leave the pair state unchanged.
Therefore, we have,
P2 = 1 . (5.14)
and writing P = P † PP = P, we see that P † is Hermitian, i.e. the eigenvalues are
real and have to ±1 for the norm to be conserved.
5.1. QUANTIZING SCALAR FIELDS 185
Any pair state Ψ(r1 , σ1 ; r2 , σ2 ) can be written as the sum of a symmetric (+) and
an antisymmetric (−) part, as we will demonstrate in Exc. 5.1.3.2. Therefore, the
eigenstates of P span the full Hilbert space of the pair and P is not only Hermitian
but also an observable. We do not deepen the discussion about the relation between
spin and statistics and only mention that the bosons turn out to have integral total
spin and the fermions half-integral total spin, so that the eigenvalues of the exchange
operator are given by,
where s is the (semi-)integral spin of the particle. Apparently, for identical particles
the pair wavefunction has to be an eigenfunction of the exchange operator, i.e. the
exchange symmetry of the wavefunction is conserved in time and the pair Hamiltonian
Ĥ is invariant under exchange of the two particles,
For fermions in the same spin state (σ1 = σ2 ≡ σ) and at the same position
(r1 = r2 ≡ r) this condition can only be satisfied, if Ψ(r, σ; r, σ) = −Ψ(r, σ; r, σ).
Hence, two fermions in the same spin state have zero probability to be found at the
same position. Therefore, the fermions show correlated motion. Importantly, these
kinematic correlations occur irrespective of the presence or absence of forces between
the particles.
We first look at the symmetry of the spin states. A Clebsch-Gordan decomposition
(see Sec. ?? and ??),
X
|S, MS ⟩ = |s1 , ms1 ; s2 , ms2 ⟩⟨s1 , ms1 ; s2 , ms2 |S, MS ⟩ (5.18)
mS1 ,mS2
1 1
yields for the 2 × 2 two fermion spin system,
|1, +1⟩ = q
| ↑↑⟩
1 (S = 1)
|1, 0⟩ = 2 (| ↑↓⟩ + | ↓↑⟩)
(5.19)
|1, −1⟩ = q
| ↓↓⟩
Hence, assuming the fermion to be in one of the symmetric spin states, S = 1, the
orbital wavefunction must be antisymmetric,
q
1 ık1 ·r1 ık2 ·r2
ΨAk1 ,k2 (r1 , r 2 )|1, MS ⟩ = 1
V 2 (e e − eık1 ·r2 eık2 ·r1 )|1, MS ⟩ . (5.21)
If the fermions are in different orbital states (k1 ̸= k2 ) this gives rise to the above
mentioned kinematic correlations. For two fermions in the same orbital state (k1 =
k2 = k) the total wavefunction Eq. (5.20) vanishes and also its norm |Ψ− k,k (r1 , r2 )|
2
is zero. Apparently two (identical) fermions cannot occupy the same state; such
a coincidence is entirely destroyed by interference. This is the essence of Pauli ’s
exclusion principle which holds for all fermions. Starting from an antisymmetric spin
state the orbital part should symmetric,
1 1 ık1 ·r1 ık2 ·r2
ΨA
k1 ,k2 (r1 , r2 )|0, MS ⟩ = V 2 (e e + eık1 ·r2 eık2 ·r1 )(| ↑↓⟩ − | ↓↑⟩) . (5.22)
In this case, no restriction is found for the relative positions of the fermions.
We found that the quantum mechanical indistinguishability of identical particles
affects the distribution of particles over the single-particle states. Remarkably, these
kinematic correlations happen in the complete absence of forces between the particles:
it is a purely quantum statistical effect.
Example 37 (Kinematic correlation in ortho and para isomers of the
hydrogen and deuterium): A consequence of such correlations can be observed
in the rotational properties of homonuclear diatomic molecules: depending on
the total spin of the molecules either the even or the odd rotational levels are
observed. This is illustrated in Tab. 5.1 for the ortho and para isomers of
hydrogen (bosonic atoms) and deuterium (fermionic atoms).
Table 5.1: In the ortho and para isomers of the hydrogen and deuterium molecule the
distribution over the rotational levels is affected by quantum statistics (only even or odd
levels) [49].
symmetry
1
species J I Ψnucl Ψrot Ψvib Σ+
g Ψtot
ortho-H2 1, 3, 5, .. 1 S A S A S
ortho-D2 0, 2, 4, .. 1 A S S A A
para-H2 0, 2, 4, .. 0 A S S A S
para-D2 1, 3, 5, .. 0, 2 S A S A A
For k1 ̸= k2 this form is appropriate because it is symmetric and also has the proper
normalization of unity, ⟨k1 , k2 |k1 , k2 ⟩ = 1 (see Exc. 5.1.3.2). For k1 = k2 = k the
situation is different. Eq. (5.23) yields norm 2 rather than the physically required
value unity. In this case the properly symmetrized and normalized wavefunction is
the product wavefunction,
Ψk,k (r1 , r2 ) = 1
V eık1 ·r1 eık2 ·r2 , (5.24)
With this, Eq. (5.21) can be written in the form of a determinant. Taking, for instance,
the spin state |S, MS ⟩ = | ↑, ↑⟩, we write, introducing the abbreviation (i) ≡ (ri , σi ),
ϕk (r1 ) = ⟨r1 |k⟩ = ⟨r|k⟩1 and χs,ms (σ1 ) = ⟨σ1 |s, ms ⟩ = ⟨σ|s, ms ⟩1
denote the spatial and the spin wavefunction of particle 1 in state |k, s, ms ⟩.
Hence,
Hence,
Θki ,si ,msi ;kj ,sj ,msj (r1 , σ1 ; r2 , σ2 ) = ⟨r1 |⟨r2 |⟨σ1 |⟨σ2 |ki , si , msi ; kj , sj , msj ⟩ .
188 CHAPTER 5. QUANTUM THEORY OF NON-RELATIVISTIC SCALAR PARTICLE FIELD
Similarly, the symmetric spin state |1, 0⟩ in combination with the antisymmetric
orbital state Ψ−
k1 ,k2 (r1 , r2 ) can be written as the sum of two determinants,
where ψkj ↕ (i) = ψkj (ri )χ↕ (σi ) with i = 1, 2. The two-body state (5.22) consisting of
an antisymmetric spin state and a symmetric orbital state (k1 = k2 ) takes the form
Indeed, the property of determinants to vanish when two columns or two rows are
equal assures that the wavefunction vanishes when two fermions are in the same state
α or share the same (position and spin) coordinates (i) = (ri , σi ), while exchanging
two rows or two columns yields the minus sign required for antisymmetric wavefunc-
tions. One can easily show that any two-body fermion state can be expressed as a
linear combination of determinantal spin-orbital states.
In this form the determinant is called a Slater determinant. It is the simplest gener-
alization of the product wavefunction with the proper symmetry under interchange
of any two fermions and consistent the Pauli principle. In Dirac notation the anti-
symmetrized form of N fermions in states α1 , ..., αN is given by
X
|α1 , ..., αN ⟩ ≡ N1 ! (−1)p P |α1 , ...αN ) , (5.30)
P
where |α1 , ..., αN ) ≡ |α1 ⟩1 ⊗ |α2 ⟩2 ⊗ ...|αN ⟩N is the N -body product state of the
single-particle states |ακ ⟩i , where κ = 1, .., N is the state index and i = 1, .., N is
the particle index. The sum runs over all permutations P of the particles, p being
the parity (number of transpositions; i.e., binary interchanges) required to realize the
permutation starting from an initial ordering fixed by convention. As the sum runs
over all permutations, it makes no difference whether we permute all particles or
permute all states of the particles. We choose the permutation operator P to act on
the state index (κ) and not on the particle index (i). With this choice, the interchange
of the states of particles 1 and 2 is written as,
To assure a uniquely defined sign of the Slater determinants we shall adopt the stan-
dard ordering convention of atomic configurations (see below). The state labeling
α1 , ...αN represents both the orbital and the spin quantum numbers, e.g. ψακ (ri , σi ) =
ψk1 ↑ (ri , σi ) and is called spin orbital.
5.2. OCCUPATION NUMBER REPRESENTATION 189
5.1.3 Exercises
5.1.3.1 Ex: Norm of symmetrized wavefunctions
Show that ψk1 ,k2 (r1 , r2 ) has the norm N = 1 for k1 ̸= k2 , but N = 1 for k1 = k2 .
where γ = {1, 1, ..., 2, 2, ..., nl } with l ≲ N . In this way the states take the shape of
number states, which are the basis states of the occupation number representation
(see next section). For the case of N bosons in the same state |ks ⟩ the number state
is given by |ns ⟩ ≡ |ks , ..., ks ⟩; for a single particle in state |ks ⟩ we have |1s ⟩ ≡ |ks ⟩.
Note that the Bose symmetrization procedure puts no restriction on the value or
order of the occupations n1 , ..., nl as long as they add up to the total number of
particles, n1 + n2 + ... + nl = N . For fermions the notation is the same but because
the wavefunction changes sign under permutation the order in which the occupations
are listed becomes subject to convention (for instance in order of growing energy of
the states). Up to this point and in view of Eqs. (??) and (??) the number states
(5.32) have normalization,
⟨n′1 , n′2 , ...|n1 , n2 , ...⟩ = δn′1 ,n1 δn′2 ,n2 ... (5.33)
and closure
′
X
|n1 , n2 , ...⟩⟨n1 , n2 , ...| = I , (5.34)
n1 ,n2 ,...
where the prime indicates that the sum over all occupations equals the total number
of particles, n1 + n2 + ... + nl = N . This is called closure within HN .
where the r|s⟩ = ψs (r) are the wavefunctions of an arbitrary single-particle basis
{|s⟩}. With this transformation we introduced two operator densities,
X X
ψ̂(r) ≡ ψ(r)âs and ψ̂ † (r) ≡ ψ ∗ (r)â†s , (5.36)
s s
which are called field operators in view of their dependence on position. In terms of
these field operators the total number operator takes the form,
Z Z
†
N̂ = d rψ̂ (r)ψ̂(r) = d3 rn̂(r) ,
3
(5.37)
where we defined the density operator n̂(r) as the diagonal part of the density matrix
operator. The field operators are construction operators that create or annihilate
particles at a given position. Let us demonstrate this for ψ̂(r): This field operator is
a creation operator because it is defined in terms of creation operators,
X X
ψ̂ † (r)|0⟩ = ψs∗ (r)â†s |0⟩ = |s⟩⟨s|r⟩ = |r⟩ . (5.38)
s s
5.2. OCCUPATION NUMBER REPRESENTATION 191
Using the closure relation, we found that ψ̂ † (r) creates from the vacuum a particle
in state |r⟩; i.e., a particle at position r. Similarly we can show that ψ̂(r) is the
corresponding annihilation operator,
The field operators are important quantities because (at least in principle) the po-
sitions of the particles can be measured to arbitrary accuracy in any many-body
system, also when the concept of stationary single-particle states has lost meaning
due to coupling by the interactions. Do the Exc. 5.2.4.2.
This expression follows immediately from Eq. (5.36). As only a single single-
particle state is occupied, only a single term contributes. If in contrast, |ΨN ⟩ =
|ñγ ⟩ is a pure number state of the single-particle representation {|s⟩}, with
many single-particle states are occupied, the expression is replaced by the linear
combination, X√
⟨ΨN −1 |ψ̂(r)|ΨN ⟩ = ns ψs (r) .
s
[ψ̂(r), ψ̂ † (r′ )]± = δ(r − r′ ) and [ψ̂(r), ψ̂(r′ )]± = 0 = [ψ̂ † (r), ψ̂ † (r′ )]± . (5.40)
Like previously, the commutators (−) refer to to case of bosons and the anti-commutators
(+) to the case of fermions. Further we have for both bosons and fermions,
[n̂(r), ψ̂ † (r′ )]± = ψ̂ † (r)δ(r − r′ ) and [n̂(r), ψ̂(r′ )]± = −ψ̂(r)δ(r − r′ ) . (5.41)
Using the expression (5.37) for the total number operator the latter commutation
relation leads to,
[ψ̂ † (r′ ), N̂ ]± = ψ̂(r) . (5.42)
We separated the operators into a sum of two contributions: (i) a part diagonal in the
number representation of {|s⟩}: For pure number states, |ΨN ⟩ = |ñγ ⟩, we calculate,
X X
⟨ΨN |n̂(r)|ΨN ⟩ = |φs (r)|2 ⟨ñγ |n̂s |ñγ ⟩ = |φs (r)|2 |n̂s = n̄(r) . (5.44)
s s
For pure number states the cross term contribution vanishes. For linear combinations
of number states the cross terms correspond to density fluctuations about the average.
5.2.3 Correlations
5.2.4 Exercises
5.2.4.1 Ex: Quantum statistics
Show that the Fermi-Dirac statistics and the Bose-Einstein statistics directly follow
from the fermionic, respectively, bosonic commutation rules.
for bosons.
3. Give the definition of the annihilation operator âs acting on the number state
|ns , nt , ..., nl ⟩ for bosons.
4. Evaluate â†s |0q , 0s , ...⟩ for fermions.
5. Evaluate â†s |1q , 0s , ...⟩ for fermions.
6. Evaluate â†s |1q , 1s , ...⟩ for fermions.
7. Show that also the vacuum state |0⟩ is normalized.
8. Calculate the norm ⟨2s |2s ⟩ for bosons and fermions without using the fermion rule
â†s |1q , 1s , ...⟩ = 0.
9. Derive [â†q , â†s ] = 0 for bosons.
10. Derive [â†q , â†s ]+ = 0 for fermions.
11. Derive [âq , â†s ]+ = δq,s for fermions.
12. Show that [n̂q , â†s ]± = +â†s δq,s for both bosons and fermions.
where
p2i
H(i) = + U(ri )
2m
is the free particle contribution of particle i and H(i,j) = V(ri , rj ) represents the
potential energy of interaction between the particles i and j. At sufficiently low
temperatures this potential may be approximated by the expression
V(ri , rj ) = gδ(ri − rj ) ,
where g = (4πℏ2 /m)a and a is called the s-wave scattering length. In the occupation
number representation for a system of identical particles the Hamiltonian takes the
form
Ĥ = Ĥ (1) + Ĥ (2) ,
where
â†s′ ⟨s′ |H(1) |s⟩âs ,
X
Ĥ (1) =
s,s′
and
â†s′ â†t′ ⟨s′ , t′ |H(1,2) |s, t⟩ât âs ,
XX
Ĥ (2) = 1
2
t,t′ s,s′
Here |s⟩, |s′ ⟩ ∈ {|s⟩}, and |t⟩, |t′ ⟩ ∈ {|t⟩} represent single particle eigenstates of the
Hamiltonians H(1) and H(2) , respectively. Rewriting the Hamiltonian in terms of field
operators we obtain Z
Ĥ (1) = drψ̂ † (r)H0 ψ̂(r)
194 CHAPTER 5. QUANTUM THEORY OF NON-RELATIVISTIC SCALAR PARTICLE FIELD
Here r is the position coordinate of a single particle and r, r′ are the position coordi-
nates of a pair of particles. The total number operator is defined as the integral over
the density operator Z
N̂ = drn̂(r) .
Three commutation relations for boson field operators are crucial for the understand-
ing of the ground state of an interacting bosonic superfluid:
ϕ(r, t) ≡ µt − Φ(r, t) .
The phase Φ(r, t) is called the fluctuating phase and represents the deviation from the
dynamical phase evolution µt/ℏ of the stationary state Ψ0 (r, t). The current density
is defined as
ıℏ
j(r, t) = (Ψ∇Ψ∗ − Ψ∗ ∇Ψ) .
2m
a. Derive the following expression for the current density,
ℏ
j(r, t) = n0 (r, t)∇Φ(r, t) .
m
b. Show that the time dependence of the phase of the order parameter can be expressed
in the form,
∂ϕ ℏ2 ℏ2 1 2
−ℏ = (∇Φ)2 + U(r) + g|Ψ|2 − ∇ |Ψ| .
∂t 2m 2m |Ψ|
Hint: Use the difference of time-dependent GP equation and its complex conjugate.
5.3. QUANTUM STATISTICS 195
5.3.1 Exercises
196 CHAPTER 5. QUANTUM THEORY OF NON-RELATIVISTIC SCALAR PARTICLE FIELD
Bibliography
[1] V. S. Bagnato, G. P. Lafyatis, A. G. Martin, E. L. Raab, R. N. Ahmad-Bitar,
and D. E. Pritchard, Continuous stopping and trapping of neutral atoms, Phys.
Rev. Lett. 58 (1987), 2194.
[7] D. A. Butts and D. S. Rokhsar, Trapped Fermi gases, Phys. Rev. A 55 (1997),
4346, DOI.
[8] E. A. Cornell, Very cold indeed: The nanokelvin physics of Bose-Einstein con-
densation, J. Res. Natl. Inst. Stand. Tech. 101 (1996), 419.
[11] J. de Boer and G. E. Uhlenbeck, Construction operator formalism for many par-
ticle systems, In Studies in Statistical Mechanics, J. de Boer and G.E. Uhlenbeck
(Eds.), (North Holland, 1965)] III (1996), 212.
197
198 BIBLIOGRAPHY
[16] P. A. M. Dirac, The quantum theory of the emission and absorption of radiation,
Proc. Roy. Soc. A 114 (1927), 243.
[22] S. Giovanazzi, A. Görlitz, and T. Pfau, Tuning the dipolar interaction in quantum
gases, Phys. Rev. Lett. 89 (2002), 130401, .
[23] B. Goss-Levi, Ultra-cold fermionic atoms team up as molecules: Can they form
cooper pairs as well, Phys. Today (2004), .
[29] C. Herzog and M. Olshanii, Trapped bose gas: The canonical versus grand canon-
ical statistics, Phys. Rev. A 55 (1997), 3254, .
[32] K. Huang, Statistical mechanics, John Wiley and Sons, 1987, ISBN.
[35] P. Jordan and E. Wigner, über das paulische äquivalenzverbot, Zeitschrift für
Physik 47 (1928), 631.
[36] W. Ketterle and N. J. Van Druten, Evaporative cooling of trapped atoms, Adv.
At. Mol. Opt. Phys. 37 (1996), 181, .
[40] Mingzhe Li, Zijun Yan, Jincan Chen, Lixuan Chen, and Chuanhong Chen, Ther-
modynamic properties of an ideal Fermi gas in an external potential with U = brt
in any dimensional space, Phys. Rev. A 58 (1998), 1445, .
[46] D. S. Petrov, Three-boson problem near a narrow Feshbach resonance, Phys. Rev.
Lett. 93 (2004), 143201, .
201
202 INDEX
ideal, 92
real, 92
solutions
regular, 94
Sommerfeld expansion, 118
specific heat, 10
spin orbital, 187, 188
spin-charge separation, 179
spin-statistics theorem, 184
state equation, 21
state functions, 21
state index, 188
Stirling
James, 117
Stirling’s formula, 117, 125
stress coefficient, 31
sublimation pump, 95
temperature
absolute zero, 27
thermal de Broglie wavelength, 151
thermal expansion coefficient, 31
thermodynamic ensembles, 63
thermodynamic limit, 143
thermodynamic potential, 31
thermodynamic potentials, 29
thermoionic emission, 180
thermometer, 5
Thomson
William, 45
triple point, 81
Watt
James, 4