Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

The fundamental relation

The puzzle: why temperatures, and not energies?

Drop a hot pebble into a cold lake. A minute later the pebble is at the lake’s temperature, and it holds a vanishing share of the energy — the lake holds essentially all of it. Put two identical blocks together, one hot and one cold, and they end at the same temperature and with the same energy. Put a big block and a small one together, and they end at the same temperature again, with energies that are not remotely equal.

Module 01 watched this happen quantum by quantum and wrote down where it stops. What it never said is why it stops there. Nothing microscopic forbids the flow from carrying on past equal temperatures. Quanta keep hopping in both directions for ever afterwards; each hop conserves energy exactly; a state with more energy in the pebble is just as legal as the one the pebble settles in.

This module answers it, and the answer turns out to be much larger than the question. There is a single function of a system’s state, its entropy written as S(U,V,N)S(U, V, N). Its slopes are the temperature, the pressure and the chemical potential. Its maximum under a constraint is the equilibrium state. Its curvature decides whether that equilibrium is stable. And the one property of being extensive forces two exact relations among everything else. Hand over S(U,V,N)S(U, V, N), and every property of the system follows from it.

Predict before you calculate

Commit to an answer for each of these before reading on. Write them down.

  1. A hot 2 kg2\ \mathrm{kg} block of copper is pressed against a cold 0.2 kg0.2\ \mathrm{kg} block of copper, and the pair is wrapped in insulation. When they stop changing, do they have equal energies, equal temperatures, both, or neither?

  2. Module 01’s simulation moves one quantum per step between two solids. Track their total entropy step by step as they relax. Can it ever tick down, even for a single step?

  3. A gas doubles its volume by leaking into an empty chamber through a pinhole, with no heat exchanged and no work done. Does its entropy rise, fall, or stay fixed?

  4. Suppose some substance had an entropy S(U)S(U) that bulged the wrong way — curving upwards over some range of energy. Put two identical pieces of it in contact at the same temperature. What happens?

Explore the model

The laboratory’s centrepiece is a composite: two subsystems inside one isolated box, separated by an internal wall. The totals of energy, volume and particle number are fixed. The wall decides which of them may be redistributed: a wall that conducts heat lets energy move, a wall that slides like a piston lets volume move, and a wall with holes in it lets particles move. The explorer asks one question over and over: of all the redistributions the wall now permits, which one has the largest total entropy?

Here is the simplest version, two Einstein solids — module 01’s model — with a wall that conducts heat.

Two Einstein solids share a fixed total energy; the horizontal axis is the fraction of it that solid A holds. Above: their total entropy, measured from its largest value, which peaks once. Below: the slope of each solid’s own entropy — the red one for A, the blue one for B — which is each solid’s 1/T1/T. A’s slope falls as A takes more energy and warms, and B’s rises as B is left with less and cools. The dots sweep across the partitions, and the dotted segment is the gap between the two slopes. Wherever the gap is open the total is still climbing toward the peak; directly beneath the peak, and only there, the gap closes. Solid B is ten times the size of A, and the peak sits at one-eleventh of the energy in A — not at half.

The second picture is a gas in a box divided by a wall, and it releases the constraints one at a time.

Contours of the total entropy of a divided monatomic gas, over the share of the energy and the share of the volume on side A; brighter is higher. The dot starts where the gas was prepared, off the summit. First the wall is made to conduct heat: only the energy share can move, so the dot travels horizontally, and it stops where the temperatures agree — the dashed line. Then the wall is also freed to slide: the dot climbs along that line to the summit, where the pressures agree too. Punching holes in the wall after that moves nothing at all, for a reason the derivation gives.

And the third runs module 01’s simulation again, this time keeping the entropy books.

Above: module 01’s two solids relaxing to a shared temperature, quantum by quantum. Below: the pair’s total entropy, counted exactly from the number of microstates of each solid at every step, minus its starting value. It rises while the temperatures converge and flattens exactly when they meet, just under the dashed line — the largest total entropy any partition of this energy has. Nothing in the simulation knows about entropy; it only moves quanta at random.

Run the laboratory now, then come back. Four things are worth doing before you read on:

Derive the result

System and boundary. A single system described by its energy UU, volume VV and particle number NN, and, for most of the argument, a composite of two such systems inside one isolated box, separated by an internal wall.

Independent variables and constraint. UU, VV and NN of each subsystem. Isolation fixes the three totals; the wall decides which of them may be shared out differently between the two sides.

Sign convention. The course convention holds throughout,

dU=δQ+δWon,\mathrm{d}U = \dbar Q + \dbar \Won ,

with δQ\dbar Q positive into the system and δWon\dbar \Won the work done on it.

Route. Eight steps, each using only the ones before it. First the four properties an entropy function must have, and where modules 07 and 08 already established them. Then its slopes, which turn out to be the temperature, pressure and chemical potential. Then the central result, that equilibrium is where those slopes match; then a worked fundamental relation, the Einstein solid’s, which settles the puzzle. Then two exact identities that extensivity alone forces; then why the surface must curve downwards; and finally how much entropy an irreversible process produces.

Step 1: what an entropy function has to do

Modules 07 and 08 built entropy twice, from opposite ends. Module 07 defined it by heat,

ΔS=∫abδQrevT,\Delta S = \int_a^b \frac{\dbar Q_{\text{rev}}}{T} ,

along any reversible path, and proved that the integral depends only on the endpoints. Module 08 defined it by counting, S=kBln⁡ΩS = \kB \ln \Omega, and proved that the logarithm is forced by additivity. This module stops asking where entropy comes from and asks what it does. From here on we take four of its properties as the starting point.

These are postulates in the sense that everything below is derived from them and nothing below derives them. But they are not guesses. Each has already been met in this course.

Extensivity deserves a warning. Module 08 measured that kBln⁡Ω\kB \ln \Omega is extensive only in the limit of large NN, with corrections of order ln⁡N\ln N that never vanish exactly. The postulate takes the limit as given. It is excellent for a gram of anything and wrong for a nanoparticle, and it fails outright for a star — step 7 comes back to that.

Step 2: the slopes of the surface are the intensive variables

Write the change in SS between two neighbouring equilibrium states as

dS=(∂S∂U)V,NdU+(∂S∂V)U,NdV+(∂S∂N)U,VdN.\mathrm{d}S = \left(\frac{\partial S}{\partial U}\right)_{V,N} \mathrm{d}U + \left(\frac{\partial S}{\partial V}\right)_{U,N} \mathrm{d}V + \left(\frac{\partial S}{\partial N}\right)_{U,V} \mathrm{d}N .

The three partial derivatives are properties of the state, like SS itself. Give them names.

The first two had better be the temperature and pressure you already know, and the definitions have to be checked against them rather than assumed. Here is the check.

Temperature. Heat the system reversibly at fixed VV and NN. No work is done, so the first law and module 07’s definition of entropy give

dU=δQrev,dS=δQrevT.\mathrm{d}U = \dbar Q_{\text{rev}} , \qquad \mathrm{d}S = \frac{\dbar Q_{\text{rev}}}{T} .

Divide one by the other: the slope ∂S/∂U\partial S/\partial U is 1/T1/T, with TT the thermodynamic temperature of module 07.

Pressure. Compress the system reversibly and adiabatically at fixed NN. No heat flows, so its entropy does not change, while the first law gives

dU=δWon=−Pmech dV,\mathrm{d}U = \dbar \Won = -P_{\text{mech}}\,\mathrm{d}V ,

with PmechP_{\text{mech}} the force per area the piston feels. Setting dS\mathrm{d}S to zero in the expansion above,

0=1T(−Pmech dV)+(∂S∂V)U,NdV,0 = \frac{1}{T}\left(-P_{\text{mech}}\,\mathrm{d}V\right) + \left(\frac{\partial S}{\partial V}\right)_{U,N} \mathrm{d}V ,

so T (∂S/∂V)T\,(\partial S/\partial V) is exactly the mechanical pressure.

Substituting the three slopes into dS\mathrm{d}S and solving for dU\mathrm{d}U gives the most used equation in thermodynamics.

A test on a real relation. For a monatomic ideal gas the fundamental relation is known in closed form.

Take its slopes. Only one term depends on UU, and only one on VV:

1T=∂S∂U=32 NkBU,PT=∂S∂V=NkBV.\frac{1}{T} = \frac{\partial S}{\partial U} = \frac{3}{2}\,\frac{N \kB}{U} , \qquad \frac{P}{T} = \frac{\partial S}{\partial V} = \frac{N \kB}{V} .

The first is module 04’s kinetic temperature, and the second is module 02’s ideal gas law:

U=32NkBT,PV=NkBT.U = \tfrac{3}{2} N \kB T , \qquad PV = N \kB T .

Neither was put in. Both came out of a single function by differentiating it, which is the whole claim of this module in miniature.

Step 3: equilibrium is where the slopes match

Now the central result. Take the composite: two systems in an isolated box, with a wall between them that conducts heat but does not move and has no holes. The total energy is fixed; how it is shared is not. Give subsystem 1 the energy U1U_1, and subsystem 2 the rest:

Stot(U1)=S1(U1)+S2(U−U1).S_{\text{tot}}(U_1) = S_1(U_1) + S_2(U - U_1) .

By the maximum postulate, the wall’s energy settles where this is largest. At a maximum the derivative vanishes, and by the chain rule the second term contributes a minus sign:

dStotdU1=∂S1∂U1−∂S2∂U2=1T1−1T2=0.\frac{\mathrm{d} S_{\text{tot}}}{\mathrm{d} U_1} = \frac{\partial S_1}{\partial U_1} - \frac{\partial S_2}{\partial U_2} = \frac{1}{T_1} - \frac{1}{T_2} = 0 .

The pressure and particle versions come the same way: release the volume at fixed total volume and differentiate with respect to V1V_1, and the slopes P/TP/T must match; with the temperatures already equal, so must the pressures. One caution the laboratory enforces: a wall that moves or lets particles through always carries energy across too — a moving piston does work, and a crossing particle takes its energy with it — so those constraints are only ever released together with the energy one.

Which way does energy flow? Before equilibrium, the total entropy is still rising, so a small transfer dU1\mathrm{d}U_1 into subsystem 1 must have

dStot=(1T1−1T2)dU1>0.\mathrm{d} S_{\text{tot}} = \left(\frac{1}{T_1} - \frac{1}{T_2}\right) \mathrm{d}U_1 > 0 .

If T1<T2T_1 < T_2, the bracket is positive, and so is dU1\mathrm{d}U_1: energy flows into the colder body. A steeper S(U)S(U) — a larger 1/T1/T — means a colder body, and giving it energy buys more entropy than the hotter body loses by giving that energy up. Heat flowing from hot to cold is not a separate law. It is the maximum postulate, applied one transfer at a time.

The puzzle, resolved. The condition is on slopes, and nothing ties a slope to a size. A big body and a small body can have the same slope with wildly different energies, and at equilibrium they do. The pebble ends at the lake’s temperature holding almost none of the energy; the lake’s S(U)S(U) is so gently curved that absorbing the pebble’s heat barely changes its slope, which is exactly what makes it a reservoir. What is maximised is the total entropy, and the only thing that condition ever equalises is ∂S/∂U\partial S/\partial U.

Is it really a maximum? A vanishing derivative could be a minimum too. The second derivative settles it:

d2StotdU12=∂2S1∂U12+∂2S2∂U22.\frac{\mathrm{d}^2 S_{\text{tot}}}{\mathrm{d} U_1^2} = \frac{\partial^2 S_1}{\partial U_1^2} + \frac{\partial^2 S_2}{\partial U_2^2} .

By the concavity postulate both terms are negative or zero, so the equal-slope point is a maximum. Step 7 shows what goes wrong without that postulate.

Step 4: a fundamental relation worked through — the Einstein solid

Module 01’s solid is nn oscillators sharing qq identical quanta of energy ε\varepsilon, so its energy is U=qεU = q\varepsilon. How many ways can qq identical quanta be shared among nn oscillators? Line the quanta up as dots and separate the oscillators with n−1n - 1 bars: every arrangement of the dots and bars in a row is one microstate. With four quanta among three oscillators, the row ●●|●|● gives the first oscillator two quanta and each of the others one. So the count is the number of ways to choose which of the q+n−1q + n - 1 positions hold dots:

Ω(q,n)=(q+n−1q).\Omega(q, n) = \binom{q + n - 1}{q} .

For four quanta among three oscillators that is 15. For macroscopic qq and nn, Stirling’s approximation from module 08 turns kBln⁡Ω\kB \ln \Omega into

S(U,n)=kB[(q+n)ln⁡(q+n)−qln⁡q−nln⁡n],q=Uε.S(U, n) = \kB \Big[ (q + n) \ln (q + n) - q \ln q - n \ln n \Big] , \qquad q = \frac{U}{\varepsilon} .

This is a fundamental relation. The solid has no volume to speak of, so SS depends on UU and on its size nn; scale both together and SS scales with them, so it is extensive. Its energy slope is the temperature:

1T=∂S∂U=kBεln⁡ ⁣(1+nεU).\frac{1}{T} = \frac{\partial S}{\partial U} = \frac{\kB}{\varepsilon} \ln\!\left(1 + \frac{n \varepsilon}{U}\right) .

When the solid is hot, UU is much larger than nεn\varepsilon, so the logarithm’s argument is only slightly above one and the logarithm can be expanded in powers of the small excess. Keeping the first two terms,

T≈UnkB+ε2kB,T \approx \frac{U}{n \kB} + \frac{\varepsilon}{2 \kB} ,

which is module 01’s map from quanta to temperature, plus an offset of half a quantum that matters less and less as the solid warms. Module 01 assumed that map by equipartition; here it has come out of counting. Inverted, the exact slope gives

U=nεeε/kBT−1,U = \frac{n \varepsilon}{e^{\varepsilon / \kB T} - 1} ,

Einstein’s formula of 1907, whose falling heat capacity at low temperature module 16 takes up.

Two unequal solids. Put solid A, with nAn_A oscillators, against solid B, with nBn_B. Equal slopes means equal logarithms, and so

nAqA=nBqB⟹UAUB=nAnB.\frac{n_A}{q_A} = \frac{n_B}{q_B} \qquad\Longrightarrow\qquad \frac{U_A}{U_B} = \frac{n_A}{n_B} .

Each oscillator ends with the same average energy, so each solid’s share of the energy is its share of the oscillators. With nB=10 nAn_B = 10\,n_A, solid A ends with one part in eleven. That is prediction 1 answered: equal temperatures, and energies in the ratio of the sizes.

Step 5: extensivity forces the Euler relation

The first postulate said that scaling a system scales its entropy. As an equation, for every λ>0\lambda > 0,

S(λU,λV,λN)=λ S(U,V,N).S(\lambda U, \lambda V, \lambda N) = \lambda\, S(U, V, N) .

A function with this property is called homogeneous of degree one, and there is a short theorem about such functions. Differentiate both sides with respect to λ\lambda. On the left the chain rule gives one term per argument; on the right the result is just SS:

U∂S∂(λU)+V∂S∂(λV)+N∂S∂(λN)=S(U,V,N).U \frac{\partial S}{\partial (\lambda U)} + V \frac{\partial S}{\partial (\lambda V)} + N \frac{\partial S}{\partial (\lambda N)} = S(U, V, N) .

This holds for every λ\lambda, so set λ=1\lambda = 1 and replace each slope by its name from step 2:

S=UT+PVT−μNT.S = \frac{U}{T} + \frac{P V}{T} - \frac{\mu N}{T} .

Multiply by TT and rearrange.

There is a second way to see it, which shows where extensivity enters. Build the system up from nothing by adding more and more of it at a fixed temperature, pressure and chemical potential, as a crystal grows in its own melt. Along that one path TT, PP and μ\mu never change, so the differential form of step 2 integrates term by term, and it gives exactly the Euler relation. What makes that path available is extensivity: only if a bigger piece of the same stuff has the same intensive state can the system be grown without changing TT, PP and μ\mu. Along any other path they change, and the naive integration is wrong.

For one mole of argon at 298.15 K298.15\ \mathrm{K} and 1 bar1\ \mathrm{bar}, Sackur–Tetrode gives these values, and the Euler relation can be checked with a calculator:

TS=46 167 J,PV=2479 J,μN=−39 970 J,TS−PV+μN=3718 J=32RT.\begin{gathered} T S = 46\,167\ \mathrm{J} , \qquad P V = 2479\ \mathrm{J} , \qquad \mu N = -39\,970\ \mathrm{J} , \\ T S - P V + \mu N = 3718\ \mathrm{J} = \tfrac{3}{2} R T . \end{gathered}

Step 6: the Gibbs–Duhem relation

Now there are two expressions for UU in hand. Differentiate the Euler relation, using the product rule on each term:

dU=T dS+S dT−P dV−V dP+μ dN+N dμ.\mathrm{d}U = T\,\mathrm{d}S + S\,\mathrm{d}T - P\,\mathrm{d}V - V\,\mathrm{d}P + \mu\,\mathrm{d}N + N\,\mathrm{d}\mu .

Subtract the differential form from step 2. The three terms of that equation cancel, and what is left must vanish.

This is why, in the second animation, punching holes in the wall moved nothing. Once the temperatures and the pressures on the two sides agreed, the Gibbs–Duhem relation left the chemical potentials no freedom to differ: the same gas at the same TT and PP has the same μ\mu. The three equilibrium conditions of step 3 are not three independent demands.

A useful special case: hold the temperature fixed. Then the change in μ\mu is V/NV/N times the change in pressure, and for an ideal gas V/NV/N is kBT/P\kB T/P, so

μ(T,P2)−μ(T,P1)=kBTln⁡P2P1.\mu(T, P_2) - \mu(T, P_1) = \kB T \ln\frac{P_2}{P_1} .

Module 13 builds on this line.

Step 7: concavity is stability

The postulates asked for a concave SS. Here is why a system that broke the rule could not exist in equilibrium with anything.

Differentiate the temperature slope once more. With 1/T1/T as the first derivative,

∂2S∂U2=∂∂U ⁣(1T)=−1T2 ∂T∂U=−1T2CV,\frac{\partial^2 S}{\partial U^2} = \frac{\partial}{\partial U}\!\left(\frac{1}{T}\right) = -\frac{1}{T^2}\,\frac{\partial T}{\partial U} = -\frac{1}{T^2 C_V} ,

where the heat capacity CVC_V is the derivative (∂U/∂T)V,N(\partial U/\partial T)_{V,N}. So the sign of the curvature is the sign of the heat capacity, reversed.

The runaway. Take two identical pieces of a substance, each with energy UU, in contact. Move a small amount of energy δ\delta from one to the other. The total entropy goes from 2 S(U)2\,S(U) to

S(U+δ)+S(U−δ).S(U + \delta) + S(U - \delta) .

If SS curves downwards, the chord between those two points lies below the curve, the new total is smaller, and the transfer is undone: the fluctuation dies. If SS curves upwards, the new total is larger, the maximum postulate rewards the transfer, and the next transfer in the same direction is rewarded more. The energy segregates. In the language of temperature: with a negative heat capacity, the piece that loses energy gets hotter, so heat flows back out of it faster, and the difference grows instead of shrinking.

Two plots of entropy against energy. Left, a curve bending downwards, with a chord between two points lying below the curve at the midpoint. Right, a curve bending upwards, with the chord lying above it.

The runaway in one picture. Two identical pieces start at the middle energy UU. Sharing out a fluctuation ±δ\pm\delta moves them to the two ends of the chord, and the chord’s midpoint is their new average entropy. (a) A concave SS: the chord’s midpoint is below the curve, so the fluctuation lowers the total entropy and is undone. (b) A convex SS: the midpoint is above the curve, so the fluctuation raises the total entropy and grows.

The volume direction. With two variables in play, “curves downwards in every direction” has a precise meaning. Write SUUS_{UU}, SVVS_{VV} and SUVS_{UV} for the second partial derivatives of S(U,V)S(U, V) at fixed NN, the last being the mixed one. A small step dU\mathrm{d}U, dV\mathrm{d}V in any direction must change SS, at second order, by a negative amount:

SUU dU2+2SUV dU dV+SVV dV2≤0.S_{UU}\,\mathrm{d}U^2 + 2 S_{UV}\,\mathrm{d}U\,\mathrm{d}V + S_{VV}\,\mathrm{d}V^2 \le 0 .

A quadratic expression like this is never positive, whatever the step, exactly when

SUU≤0,SVV≤0,SUUSVV−SUV2≥0.S_{UU} \le 0 , \qquad S_{VV} \le 0 , \qquad S_{UU} S_{VV} - S_{UV}^2 \ge 0 .

Now squeeze the system at fixed temperature. Its first slope SUS_U is 1/T1/T, so keeping TT fixed means an energy change must accompany the volume change and cancel its effect on that slope:

SUU dU+SUV dV=0.S_{UU}\,\mathrm{d}U + S_{UV}\,\mathrm{d}V = 0 .

The volume slope SVS_V is P/TP/T, and with TT fixed its change is the change in pressure divided by TT. Substituting the energy change from the line above,

dPT=SUV dU+SVV dV=(SVV−SUV2SUU)dV,\frac{\mathrm{d}P}{T} = S_{UV}\,\mathrm{d}U + S_{VV}\,\mathrm{d}V = \left( S_{VV} - \frac{S_{UV}^2}{S_{UU}} \right) \mathrm{d}V ,

and so

(∂P∂V)T=T (SUUSVV−SUV2)SUU≤0.\left(\frac{\partial P}{\partial V}\right)_T = \frac{T\,\big(S_{UU} S_{VV} - S_{UV}^2\big)}{S_{UU}} \le 0 .

The numerator is not negative and the denominator is negative, so compressing at fixed temperature can only raise the pressure. That is the statement that the isothermal compressibility,

κT=−1V(∂V∂P)T,\kappa_T = -\frac{1}{V}\left(\frac{\partial V}{\partial P}\right)_T ,

is never negative. On the Sackur–Tetrode surface the mixed derivative vanishes, and the formula returns exactly one over the pressure — the ideal gas’s compressibility, read off the curvature of its entropy. The runaway argument has a mechanical twin here: if a substance’s pressure fell when it was squeezed, a movable wall between two identical samples would be pushed further by any small displacement, and the pair would separate into a dense region and a dilute one. That is exactly what the van der Waals loop below its critical temperature would do, and it is why a real gas condenses there instead.

The energy runaway, for its part, is not hypothetical. A star is a ball of gas held together by its own gravity, and the virial theorem says its total energy is minus its kinetic energy:

U=−32NkBT,C=−32NkB.U = -\tfrac{3}{2} N \kB T , \qquad C = -\tfrac{3}{2} N \kB .

A star that radiates energy away heats up. Two stars could never settle at a common temperature, and a star cannot be put in equilibrium with a heat bath. It escapes the argument by breaking the first postulate: gravity is long-ranged, a star twice the mass is not two stars, and its entropy is not extensive. The laboratory includes this relation, and the Euler relation of step 5 fails on it by exactly twice its energy.

Step 8: entropy produced, computed

Module 07 defined the entropy produced by a process as the total change of the system and of everything it exchanged heat with, and proved it is never negative. This step computes it for two processes that module 07 could only bound.

Two bodies at different temperatures. Put body A, heat capacity CAC_A, at TA,0T_{A,0} against body B, CBC_B at TB,0T_{B,0}, and insulate the pair. Energy conservation fixes where they end, as module 01 found:

Teq=CATA,0+CBTB,0CA+CB.T_{\text{eq}} = \frac{C_A T_{A,0} + C_B T_{B,0}}{C_A + C_B} .

Each body’s entropy change is module 07’s constant-heat-capacity result, and nothing crosses the pair’s outer boundary, so their sum is all produced:

ΔStot=CAln⁡TeqTA,0+CBln⁡TeqTB,0.\Delta S_{\text{tot}} = C_A \ln\frac{T_{\text{eq}}}{T_{A,0}} + C_B \ln\frac{T_{\text{eq}}}{T_{B,0}} .

Why is this positive? Divide by CA+CBC_A + C_B and write wAw_A and wBw_B for the two bodies’ shares of the total heat capacity. The final temperature is the weighted average of the starting ones, and the entropy change is the logarithm of that average minus the weighted average of the logarithms. The logarithm is concave, so the logarithm of an average always beats the average of the logarithms — the same chord picture as step 7:

ln⁡ ⁣(wATA,0+wBTB,0)≥wAln⁡TA,0+wBln⁡TB,0,\ln\!\big(w_A T_{A,0} + w_B T_{B,0}\big) \ge w_A \ln T_{A,0} + w_B \ln T_{B,0} ,

with equality only when the two starting temperatures are equal.

A free expansion. A gas leaks from V1V_1 into an evacuated chamber until it fills V2V_2. No heat, no work, so UU is unchanged, and so is NN. The Sackur–Tetrode equation depends on VV through one term, NkBln⁡VN \kB \ln V, so

ΔS=NkBln⁡V2V1.\Delta S = N \kB \ln\frac{V_2}{V_1} .

That makes three independent routes to one number. Module 07 took a reversible isotherm between the same two end states and integrated δQrev/T\dbar Q_{\text{rev}}/T. Module 08 counted: each atom has V2/V1V_2/V_1 times the room, so Ω\Omega grows by that factor to the power NN. And here it is a difference of two values of one function. Prediction 3 is answered three times over: the entropy rises, although nothing flowed in. Nothing in between was an equilibrium state, and it did not need to be — SS belongs to the endpoints, and both endpoints are equilibrium states.

Flow and production are different things. It is tempting to think entropy is produced wherever it moves. It is not.

A gas expanding reversibly along an isotherm at TT, in contact with a bath at the same TT, takes in heat QQ. Entropy Q/TQ/T flows into the gas and the same Q/TQ/T flows out of the bath, and nothing is produced: a great deal of flow, zero production. Now let the same heat QQ cross from a body at TAT_A into one at TBT_B with TA>TBT_A > T_B. It leaves A carrying Q/TAQ/T_A and arrives in B bringing Q/TBQ/T_B, which is more; the difference,

Sgen=QTB−QTA>0,S_{\text{gen}} = \frac{Q}{T_B} - \frac{Q}{T_A} > 0 ,

was produced in the crossing. That is the whole of the contact formula above, one sliver of heat at a time.

Verify computationally

The derivation is a chain of exact statements. The laboratory checks that each link was assembled correctly, on the two relations of this module and on a third, the star, that deliberately breaks the rules. Each check is also a test in the project’s suite.

1. The slopes are the equations of state. Numerical derivatives of the Sackur–Tetrode function, with nothing else supplied, return the temperature and pressure the state was built with to about one part in 108, and a chemical potential that matches kBTln⁡(Nλ3/V)\kB T \ln(N \lambda^3 / V), with λ\lambda the thermal wavelength, to one part in 106.

2. The Einstein slope is the Einstein temperature. At every energy tried, from half a quantum per oscillator to two hundred, the numerical slope of the Stirling-form S(U,n)S(U, n) agrees with the closed form of step 4 to one part in 107. As the solid heats, the temperature approaches module 01’s map with an offset that settles at exactly half a quantum.

3. The maximiser finds equal slopes. Released through a heat-conducting wall, a solid of 40 oscillators and one of 400 end at the same temperature to one part in 106, with the energy split exactly 1:101 : 10. The divided gas, with its wall freed to slide, ends with equal temperatures and equal pressures; perforating the wall afterwards changes the total entropy by less than one part in 1012.

4. Euler and Gibbs–Duhem close — and can fail. For both extensive relations, the Euler relation balances to about one part in 108 of UU, and the Gibbs–Duhem sum along a small displacement vanishes to one part in 105 of its terms, shrinking fourfold every time the displacement is halved. The star’s Euler residual is not small. It is 2.000, in units of ∣U∣|U|: the relation misses exactly twice the energy, as its non-extensive entropy predicts.

5. The curvature is the heat capacity. The second derivative of Sackur–Tetrode gives a heat capacity of 32NkB\tfrac{3}{2} N \kB to one part in 105. The star gives −32NkB-\tfrac{3}{2} N \kB, and the stability test rejects it. The full curvature in energy and volume, fed through the formula of step 7, gives the gas’s isothermal compressibility as one over its pressure, again to one part in 105.

6. The ledger of module 01’s simulation. Two solids of 300 and 100 oscillators, started at 500 K500\ \mathrm{K} and 250 K250\ \mathrm{K} with quanta of ε=(5 K) kB\varepsilon = (5\ \mathrm{K})\,\kB, relax as in module 01. The ledger counts each solid’s microstates exactly at every step.

The module-01 hop rule moves quanta as if they were labelled, so its own long-run distribution is not exactly the Einstein count: it peaks at almost the same partition — here within two tenths of a percent — but its jitter about the peak is narrower than a true Einstein solid’s. The rise and the plateau, which are what the ledger tests, are unaffected.

7. Three entropies of one contact. The same relaxation can be priced three ways. The exact count says 15.27 kB15.27\,\kB. The Stirling-form surface of step 4 says 15.69 kB15.69\,\kB. The contact formula of step 8, with module 01’s heat capacities C=nkBC = n \kB, says 15.90 kB15.90\,\kB.

8. Free expansion, three ways. Doubling the volume of 10 00010\,000 atoms gives NkBln⁡2N \kB \ln 2 from the Sackur–Tetrode difference, from module 07’s ideal-gas formula, and from the heat absorbed along a reversible isotherm divided by its temperature — to nine significant figures in every case.

Transfer the idea

Looking back:

Looking ahead:

And in the laboratory:

Worked examples

Try each one before opening its solution.

Worked example 1 — sharing out the quanta. Solid A has 100 oscillators and solid B has 300. Together they hold 2000 quanta of size ε\varepsilon. At equilibrium, how many quanta does each hold, and what is their common temperature?

Solution

Equal slopes means equal energy per oscillator:

qA100=qB300,qA+qB=2000⟹qA=500,qB=1500.\frac{q_A}{100} = \frac{q_B}{300} , \qquad q_A + q_B = 2000 \qquad\Longrightarrow\qquad q_A = 500 , \quad q_B = 1500 .

Each oscillator holds 5 quanta on average, so the temperature from step 4 is

T=εkBln⁡(1+1/5)=ε0.1823 kB=5.485 εkB.T = \frac{\varepsilon}{\kB \ln(1 + 1/5)} = \frac{\varepsilon}{0.1823\,\kB} = 5.485\,\frac{\varepsilon}{\kB} .

Module 01’s map would say 5 ε/kB5\,\varepsilon/\kB; the difference is almost exactly the half-quantum offset. B holds three times A’s energy at the same temperature.

Worked example 2 — mixing water. Half a kilogram of water at 360 K360\ \mathrm{K} is poured into 1.5 kg1.5\ \mathrm{kg} at 280 K280\ \mathrm{K} in an insulated flask. Take the specific heat as 4184 J kg−1 K−14184\ \mathrm{J\,kg^{-1}\,K^{-1}}. Find the final temperature, each part’s entropy change, and the entropy produced.

Solution

The heat capacities are 2092 J K−12092\ \mathrm{J\,K^{-1}} and 6276 J K−16276\ \mathrm{J\,K^{-1}}, and the final temperature is their weighted average:

Teq=0.5×360+1.5×2802.0=300.0 K.T_{\text{eq}} = \frac{0.5 \times 360 + 1.5 \times 280}{2.0} = 300.0\ \mathrm{K} .

Each part’s entropy change is Cln⁡(Teq/T0)C \ln(T_{\text{eq}}/T_0):

ΔShot=2092ln⁡300360=−381.4 J K−1,ΔScold=6276ln⁡300280=+433.0 J K−1,ΔStot=51.6 J K−1.\begin{aligned} \Delta S_{\text{hot}} &= 2092 \ln\frac{300}{360} = -381.4\ \mathrm{J\,K^{-1}} , \\ \Delta S_{\text{cold}} &= 6276 \ln\frac{300}{280} = +433.0\ \mathrm{J\,K^{-1}} , \\ \Delta S_{\text{tot}} &= 51.6\ \mathrm{J\,K^{-1}} . \end{aligned}

The hot water’s entropy fell, and nothing is wrong: the second law constrains the total. The same 125.5 kJ125.5\ \mathrm{kJ} of heat left the hot water and entered the cold, but it left from hotter water, carrying less entropy per joule than it delivered. The 51.6 J K−151.6\ \mathrm{J\,K^{-1}} was produced in the crossing.

Worked example 3 — checking Euler on argon. For one mole of argon at 298.15 K298.15\ \mathrm{K} and 1 bar1\ \mathrm{bar}, Sackur–Tetrode gives S=154.85 J K−1S = 154.85\ \mathrm{J\,K^{-1}}. Find UU, PVPV and μN\mu N, and check the Euler relation.

Solution

The slopes of step 2 give the energy and the product of pressure and volume:

U=1.5×8.3145×298.15=3718 J,PV=2479 J.U = 1.5 \times 8.3145 \times 298.15 = 3718\ \mathrm{J} , \qquad P V = 2479\ \mathrm{J} .

For μ\mu, differentiate Sackur–Tetrode with respect to NN. Every term of S/NS/N depends on NN only through V/NV/N and U/NU/N, which together contribute −52kB-\tfrac{5}{2}\kB per particle, so

μ=−T(SN−52kB),μN=−298.15×(154.85−20.79)=−39 970 J.\mu = -T\left(\frac{S}{N} - \frac{5}{2}\kB\right) , \qquad \mu N = -298.15 \times (154.85 - 20.79) = -39\,970\ \mathrm{J} .

Then

TS−PV+μN=46 167−2479−39 970=3718 J,T S - P V + \mu N = 46\,167 - 2479 - 39\,970 = 3718\ \mathrm{J} ,

which is UU. Note what μ\mu is doing: it is large and negative, because adding an atom to a dilute gas at fixed energy and volume raises the entropy by a lot.

Worked example 4 — Gibbs–Duhem at work. Argon at 298.15 K298.15\ \mathrm{K} is compressed isothermally from 1 bar1\ \mathrm{bar} to 2 bar2\ \mathrm{bar}. By how much does its molar chemical potential change?

Solution

At fixed temperature Gibbs–Duhem leaves

N dμ=V dP,N\,\mathrm{d}\mu = V\,\mathrm{d}P ,

and for an ideal gas V/NV/N is kBT/P\kB T / P. Per mole,

Δμmolar=RTln⁡P2P1=8.3145×298.15×ln⁡2=1718 J mol−1.\Delta \mu_{\text{molar}} = R T \ln\frac{P_2}{P_1} = 8.3145 \times 298.15 \times \ln 2 = 1718\ \mathrm{J\,mol^{-1}} .

The chemical potential rises with pressure: squeezed atoms have less room, so each one added at fixed energy and volume buys less entropy. Module 13 will show that particles move from higher μ\mu to lower, which is why a gas leaks from high pressure to low.

Worked example 5 — the price of a free expansion. One mole of an ideal gas at 300 K300\ \mathrm{K} doubles its volume by leaking into a vacuum. How much entropy is produced, and how much work was lost compared with the best possible route to the same final state?

Solution

The entropy change is

ΔS=Rln⁡2=5.763 J K−1,\Delta S = R \ln 2 = 5.763\ \mathrm{J\,K^{-1}} ,

and since nothing flowed in, all of it was produced. A reversible isothermal expansion against a piston, drawing heat from surroundings at 300 K300\ \mathrm{K}, would have delivered

RTln⁡2=8.3145×300×0.6931=1729 JR T \ln 2 = 8.3145 \times 300 \times 0.6931 = 1729\ \mathrm{J}

of work — which is exactly the temperature times the entropy produced, the lost-work rule of module 07.

Worked example 6 — a star. A star of NN particles radiates away energy at a rate LL. Use its heat capacity to find how its temperature changes.

Solution

With the star’s energy from step 7, losing energy at the rate LL means

dTdt=1C dUdt=−L−32NkB=+2L3NkB.\frac{\mathrm{d}T}{\mathrm{d}t} = \frac{1}{C}\,\frac{\mathrm{d}U}{\mathrm{d}t} = \frac{-L}{-\tfrac{3}{2} N \kB} = + \frac{2 L}{3 N \kB} .

It heats up. The star contracts, its gravitational energy falls by twice what it radiates, and half of that goes into making its gas hotter. Put it in contact with anything cooler and the exchange runs away rather than settling — the convex picture of step 7, in the sky.

Check your understanding

  1. This module takes four properties of the entropy as its starting point. Which of the following is one that module 08 had already PROVED from counting, rather than one that this module simply assumes?

    A. Entropies of independent subsystems add, because their multiplicities multiply and the logarithm turns a product into a sum. B. The entropy surface is concave in every direction. C. Entropy is exactly extensive at every system size. D. The equilibrium temperature is the heat-capacity-weighted average of the starting temperatures.

    Check your answer

    Correct answer: A — Entropies of independent subsystems add, because their multiplicities multiply and the logarithm turns a product into a sum.

    • A. Right. Given S = k_B ln Omega and independent subsystems, additivity is a theorem: Omega_AB = Omega_A Omega_B, so S_AB = S_A + S_B. The postulates list it anyway, because the thermodynamic argument has to stand on its own, whichever construction of S you started from.

    • B. Concavity is the one postulate neither module 07 nor module 08 established. It is assumed here, and its content is stability: a concave S is a positive heat capacity. Stars break it, so it cannot be a theorem about all matter.

    • C. Module 08 measured the opposite: k_B ln Omega is extensive only in the limit of large N, with ln N corrections that never vanish exactly. The postulate takes the limit as given, which is an assumption, not something 08 proved.

    • D. That is module 01’s result from energy conservation, not a property of the entropy function at all. This module explains why the flow stops there; it does not take it as a postulate.

  2. A hot 2 kg block of copper and a cold 0.2 kg block of copper are pressed together and insulated from everything else. What is true once they stop changing?

    A. They have equal temperatures, and the big block holds about ten times the small block’s energy. B. They have equal energies, because equilibrium shares the energy out evenly. C. They have equal temperatures and equal energies, since each follows from the other. D. Neither is equal; the flow stops when the entropy of each block separately stops rising.

    Check your answer

    Correct answer: A — They have equal temperatures, and the big block holds about ten times the small block’s energy.

    • A. Right. Maximising the total entropy sets the slopes dS/dU of the two blocks equal, and that slope is 1/T. Nothing ties a slope to a size, so the energies end in the ratio of the heat capacities, here 10 : 1.

    • B. This is the equal-energy misconception. What equilibrium maximises is the total entropy, and at the maximum the two slopes of S(U) agree, not the two energies. A pebble dropped in a lake ends at the lake’s temperature holding almost none of the energy.

    • C. Equal temperatures imply equal energies only for two identical bodies. For bodies of different sizes the same temperature means different energies, in proportion to their heat capacities.

    • D. The hot block’s entropy falls the whole time, so waiting for each block’s entropy to stop rising gives nothing. What stops rising is the TOTAL, and it peaks exactly where the temperatures agree.

  3. A table of a system’s entropy at fixed volume and particle number gives S = 2.3400 J/K at U = 99.0 J and S = 2.3466 J/K at U = 101.0 J. Estimate the system’s temperature at U = 100 J, in kelvin.

    Check your answer

    Correct answer: 303.0 (± 3.0)

    The slope is dS/dU = (2.3466 - 2.3400) / (101.0 - 99.0) = 0.0033 K^-1, and that slope is 1/T, so T = 1/0.0033 = 303 K. A central difference like this one is accurate to second order in the step. If you got 0.0033 you stopped at the slope; if you got about 7 you divided S by U, which is not a derivative.

  4. Two systems share a heat-conducting wall. At the moment, the slope dS_1/dU_1 of system 1 is larger than the slope dS_2/dU_2 of system 2. Which way does energy flow?

    A. Into system 1, which is the colder of the two. B. Into system 2, because a steeper slope means a hotter body. C. Nowhere; energy flows only when the energies differ, not the slopes. D. It cannot be decided without knowing the two heat capacities.

    Check your answer

    Correct answer: A — Into system 1, which is the colder of the two.

    • A. Right. Moving dU into system 1 changes the total entropy by (dS_1/dU_1 - dS_2/dU_2) dU, which is positive, so the transfer is favoured. A larger slope is a larger 1/T, a lower temperature: heat flows to the colder body because that raises the total entropy.

    • B. The slope is 1/T, so steeper means colder, not hotter. Check with the sign of the entropy change: moving energy into system 2 lowers the total entropy here.

    • C. Energy keeps flowing until the slopes match, whatever the energies are. Two bodies of very different energy can be in equilibrium, and two of equal energy need not be.

    • D. The heat capacities decide how far each temperature moves, and so where the flow stops. The direction is fixed by the slopes alone.

  5. For some one-component system, T = 300 K, S = 150 J/K, P = 1.00e5 Pa, V = 0.0250 m^3, and the product mu N = -40000 J. Use the Euler relation to find the internal energy U, in joules.

    Check your answer

    Correct answer: 2500 (± 20)

    U = T S - P V + mu N = 45000 - 2500 - 40000 = 2500 J. Watch the two signs: the PV term enters with a minus, and mu N is itself negative here. The Euler relation needs nothing but extensivity, so these numbers are consistent with any extensive system.

  6. A one-component system’s temperature is changed by dT and its pressure by dP. Can its chemical potential change by any amount you like?

    A. No. The Gibbs-Duhem relation fixes it at dmu = (V dP - S dT) / N. B. Yes. T, P and mu are three independent state variables. C. Only if the system is an ideal gas; otherwise it is free. D. Yes, but only if N is also allowed to change.

    Check your answer

    Correct answer: A — No. The Gibbs-Duhem relation fixes it at dmu = (V dP - S dT) / N.

    • A. Right. S dT - V dP + N dmu = 0 leaves no freedom: of the three intensive variables, only two can be chosen. This is why, once two subsystems of the same substance agree on T and P, they automatically agree on mu.

    • B. They look independent, but extensivity ties them together through Gibbs-Duhem. A single phase has only two independent intensive variables.

    • C. Gibbs-Duhem uses nothing about the ideal gas. It follows from the Euler relation, which needs only extensivity, so it holds for every substance.

    • D. The relation involves N itself, not dN; changing the amount of substance at fixed T and P changes nothing intensive. mu is fixed by T and P either way.

  7. Suppose a substance’s entropy S(U) curved UPWARDS over some range of energy. Two identical pieces of it, in that range and at the same temperature, are put in contact. What happens?

    A. Any small transfer of energy between them raises the total entropy and grows, so the pieces run away to a hot one and a cold one. B. Nothing observable; they stay in equilibrium, because their temperatures are equal. C. They stay at equal temperatures but their energies oscillate back and forth forever. D. The substance would have to be at negative absolute temperature.

    Check your answer

    Correct answer: A — Any small transfer of energy between them raises the total entropy and grows, so the pieces run away to a hot one and a cold one.

    • A. Right. On a convex curve the chord lies above the curve, so S(U + delta) + S(U - delta) exceeds 2 S(U): every fluctuation is rewarded. Equivalently the heat capacity is negative, and the piece that loses energy gets hotter. A star is such a system.

    • B. Equal temperatures make the equal split a stationary point of S_total, but on a convex curve it is a minimum, not a maximum. Any fluctuation moves the pair away, and nothing moves it back.

    • C. Oscillation would need a restoring tendency, and a convex entropy supplies the opposite: every departure from the equal split raises the total entropy further.

    • D. The slope can be positive while the curvature is positive, as it is for a star. What a convex S gives is a negative heat capacity, not a negative temperature.

  8. One kilogram of water at 350 K is mixed with one kilogram of water at 290 K in an insulated flask. Take the specific heat as 4184 J/(kg K). Compute the total entropy produced, in J/K.

    Check your answer

    Correct answer: 36.94 (± 0.5)

    T_eq = 320 K. Delta S_total = 4184 ln(320/350) + 4184 ln(320/290) = -374.9 + 411.9 = 36.9 J/K. The hot water’s entropy fell and the cold water’s rose by more. If you got zero you assumed the heat Q carried equal entropy out and in, but it left at a higher temperature than it arrived: Q/T_hot is less than Q/T_cold.

  9. An ideal gas expands slowly along an isotherm at 300 K, in contact with a bath at exactly 300 K, absorbing 600 J of heat. What happens to entropy?

    A. 2 J/K flows from the bath into the gas, and none is produced. B. 2 J/K is produced, because entropy moved into the gas. C. 4 J/K is produced, 2 J/K in the gas and 2 J/K in the bath. D. No entropy moves at all, because the temperature stayed constant.

    Check your answer

    Correct answer: A — 2 J/K flows from the bath into the gas, and none is produced.

    • A. Right. The gas gains Q/T = 2 J/K and the bath loses exactly 2 J/K: entropy has moved, and the total has not changed. Production needs an irreversibility, such as heat crossing a temperature difference; flow on its own is not one.

    • B. This is the flow-is-production misconception. The bath’s entropy fell by exactly what the gas’s rose, so the total change, which is what production means, is zero.

    • C. The bath gave up heat, so its entropy fell, it did not rise. The two changes cancel exactly for a reversible exchange.

    • D. Constant temperature does not mean constant entropy. The gas took in 600 J of heat at 300 K, and that carried 2 J/K of entropy in with it.

  10. A gas at 300 K rushes freely from 1 L into an evacuated 2 L box. None of the states in between is an equilibrium state. What can be said about the gas’s entropy change?

A. It is N k_B ln 2, found by any reversible path between the same two equilibrium end states. B. It cannot be computed, because entropy is defined only at equilibrium and the process passed through none. C. It is zero, because no heat was exchanged. D. It is negative, because the gas cooled as it expanded.

Check your answer

Correct answer: A — It is N k_B ln 2, found by any reversible path between the same two equilibrium end states.

  • A. Right. S is a function of the equilibrium state, and both end states are equilibrium states. A reversible isotherm from 1 L to 2 L connects them, and integrating delta-Q_rev/T along it gives N k_B ln 2. The same number comes from the Sackur-Tetrode function and from counting positions.

  • B. The intermediate states do not need to be equilibrium states, only the endpoints do. Entropy belongs to states, not to the route, so any reversible route between the same endpoints is a legitimate instrument for computing it.

  • C. dS = delta-Q/T holds only for a reversible exchange. In a free expansion the entropy rises by N k_B ln 2 with no heat at all: it is all produced, none received.

  • D. An ideal gas expanding into vacuum does no work and receives no heat, so its energy and temperature are unchanged. And even a gas that did cool could not lower its entropy in an isolated process.

  1. Before running the ledger in the laboratory: module 01’s two solids relax toward a common temperature, one random quantum hop per step. You track their exact total entropy after every step. Predict (a) whether it ever goes down on a single step, (b) where it ends, and (c) what its overall rise does and does not prove.

Check your answer

(a) Yes, often: on roughly a fifth of all steps a quantum hops the “wrong” way and the total entropy dips, by about 0.01 k_B at most, while the whole rise is about 15 k_B. (b) It flattens just under the maximum total entropy any partition of that energy has, at exactly the moment the temperatures meet. (c) It illustrates the maximum postulate on one run, and the fact that it is monotone only within noise is the statistical character of the second law made visible. It proves nothing: a finite run cannot establish a law, and the proof is module 08’s counting, which says the macrostates near the peak have overwhelmingly more microstates. If you predicted that it could never dip, notice that nothing in the dynamics knows about entropy at all.

  1. Einstein solid A has 100 oscillators and solid B has 300. Together they hold 2000 energy quanta. At equilibrium, how are the quanta shared?

A. 500 in A and 1500 in B, the same average per oscillator. B. 1000 in each, so that the energies are equal. C. 1500 in A and 500 in B, since the smaller solid gets hotter more easily. D. It depends on how the quanta were shared at the start.

Check your answer

Correct answer: A — 500 in A and 1500 in B, the same average per oscillator.

  • A. Right. Equal slopes of S(U, n) = k_B [(q + n) ln(q + n) - q ln q - n ln n] mean equal ln(1 + n/q), so q_A/n_A = q_B/n_B: 5 quanta per oscillator in each.

  • B. Equal energies would leave A with 10 quanta per oscillator and B with 3.3: A would be much hotter, and energy would still be flowing from A to B.

  • C. A smaller solid does heat up more for the same energy, which is exactly why it ends with LESS of the energy at a shared temperature, not more.

  • D. The starting split decides how much entropy is produced on the way, not where the pair ends. The maximum of the total entropy depends only on the total energy and the two sizes.

  1. The Sackur-Tetrode entropy of a monatomic gas depends on the volume only through a term N k_B ln V. What does the slope dS/dV at fixed U and N give?

A. P/T = N k_B / V, which is the ideal gas law. B. The pressure P itself, P = N k_B / V. C. Nothing physical; only the energy slope of S has a meaning. D. The heat capacity at constant volume.

Check your answer

Correct answer: A — P/T = N k_B / V, which is the ideal gas law.

  • A. Right. The volume slope of S is P/T, so the one ln V term alone yields P V = N k_B T. The equation of state was never put in; it was read off the fundamental relation.

  • B. The slope dS/dV is P/T, not P. Check the units: an entropy per volume is J/(K m^3), which is a pressure divided by a temperature.

  • C. All three slopes of S(U, V, N) have a meaning: 1/T, P/T and -mu/T. That the volume slope is P/T was checked on this page with a reversible adiabatic compression.

  • D. Heat capacities come from SECOND derivatives of S with respect to U: the curvature d^2S/dU^2 = -1/(T^2 C_V). A first derivative with respect to V gives P/T.

Exam-style problems for this module — the ones worth doing with a pen — are collected in the problem set.

Explain it in your own words

Answer in a few sentences each.

  1. Why does a huge body make an ideal “reservoir”? Explain it using the shape of its S(U)S(U) near the energy it holds, not the word “large”.

  2. A student says: “Entropy is produced whenever it flows.” Give a process in which a great deal of entropy flows and none is produced, and repair the student’s claim.

  3. A free expansion passes through no equilibrium states at all. Why can its entropy change nevertheless be computed, and why does the answer not depend on how it was computed?

  4. Two bodies in contact end at equal temperatures. Explain, without formulas, why “equal energies” was never a candidate for what they approach.

  5. The Euler relation and the differential form of step 2 look similar. Explain what each one assumes, and why integrating the second gives the first along only one special kind of path.

Advanced: the fixed intensive variable, and the energy picture

One intensive variable is never free. A one-component system has three intensive variables, and Gibbs–Duhem ties them together: μ\mu is a function of TT and PP, and nothing else. So the intensive state of a single phase has two degrees of freedom. Now let two phases of the same substance coexist, liquid and vapour. Equilibrium between them demands equal TT, equal PP and equal μ\mu. With the first two built in, that is one equation — the second phase’s μ(T,P)\mu(T, P) must equal the first’s — in two unknowns, leaving one free: coexistence happens along a line in the TT–PP plane, the boiling curve, and not over a region. Three phases leave none, which is why a triple point is a point. This counting is Gibbs’s phase rule, and module 14 derives it.

The energy-minimum picture. Invert the fundamental relation: since SS increases with UU, the surface can equally be written as U(S,V,N)U(S, V, N). The maximum postulate then has a twin.

The argument runs by contradiction. Suppose a composite were in equilibrium at fixed total entropy with its energy not at the minimum. Then some rearrangement at the same entropy would have lower energy. Carry that rearrangement out reversibly, drawing the excess energy off as work and returning it as heat. The entropy rises, so the original state did not have the maximum entropy for its energy, contradicting the maximum postulate. The two statements pick out the same states, in the way a circle is both the shape with the largest area for a given perimeter and the shape with the smallest perimeter for a given area. On the explorer’s surface, minimising U1+U2U_1 + U_2 along a line of constant total entropy lands on the same equal-temperature partition as maximising S1+S2S_1 + S_2 at constant total energy. The energy form is the one engineers use, and it is the one module 10 transforms.