physics2u
Tier
⌕ Search ⌘K
physics2u.com/vault/maxwell-relations.html
Derivation · the atom

The four Maxwell relations from exact differentials

D-016 Home PU-203 Threads energy Depends on D-004 verified

Statement

Because the thermodynamic potentials are state functions, their mixed second derivatives commute — and each commutation is a non-obvious relation between measurable quantities.

Why it matters

These relations connect quantities that are hard to measure to quantities that are easy. Entropy cannot be measured directly; (∂S/∂V)T equals (∂P/∂T)V — a pressure gauge and a thermometer. This is the derivation that makes thermodynamics useful rather than merely true.

Assumptions
Assumption. The potentials are state functions with continuous second derivatives.Fails exactly at a phase transition, where a second derivative diverges — and the failure is diagnostic, not a nuisance. It is how transitions are classified.
Assumption. The system is described by two independent variables.A simple compressible substance. Magnetic, elastic or multi-component systems have more variables and correspondingly more relations.
Assumption. Changes are quasi-static.dU = T dS − P dV holds only for reversible paths, though the resulting relations, connecting state functions, hold generally.
Derivation
1
dU = T dS − P dV
First and second laws combined. Natural variables are S and V, meaning U expressed in those variables carries complete thermodynamic information. A
2
T = (∂U/∂S)V ,   −P = (∂U/∂V)S
Read off by comparing with the generic differential. The subscript is not decoration: it specifies which variable is held fixed, and omitting it makes the expression meaningless. A
3
2U / ∂V ∂S = ∂2U / ∂S ∂V
Mixed partials of a smooth function commute — Clairaut's theorem, pure mathematics. This is the entire mechanism; the physics was all spent at step 1. B
4
(∂T/∂V)S = − (∂P/∂S)V
Apply step 3 to the two expressions in step 2. The first Maxwell relation, obtained from nothing but the existence of U. B
5
H = U + PV ,  F = U − TS ,  G = U + PV − TS
Legendre-transform to each choice of natural variables — the identical operation performed on the Lagrangian in D-004. Each transform swaps one variable for its conjugate: V ↔ P and S ↔ T. C
6
dH = T dS + V dP ,  dF = −S dT − P dV ,  dG = −S dT + V dP
Differentiate each and cancel, exactly as the dq̇ terms cancelled in D-004. Then apply step 3 to each potential in turn. C
Result
(∂T/∂V)S = − (∂P/∂S)V     (∂T/∂P)S = (∂V/∂S)P
(∂S/∂V)T = (∂P/∂T)V     (∂S/∂P)T = − (∂V/∂T)P

Reading. The last two are the useful ones: they express entropy derivatives — unmeasurable — in terms of a thermal expansion coefficient and an isothermal compressibility, both routinely measured.

Units check. Each side is a ratio of the same conjugate pair. In (∂S/∂V)T = (∂P/∂T)V, the left is (J·K−1)/m3 = J·K−1·m−3; the right is Pa/K = (J·m−3)/K = J·K−1·m−3. Both are entropy per volume per unit change — they match.

Limiting cases
  • Ideal gas: with P = N kB T / V, the third relation gives (∂S/∂V)T = N kB/V directly from the equation of state.
  • Reversibility: the differential is written for a reversible path, but the four relations connect state functions and therefore hold for any process between the same endpoints.
  • Single potential: even the bare existence of U (step 4 alone) already yields one relation; the other three add nothing physical, only new variable pairings.
Breaks when
  • At a phase transition, where a second derivative diverges and mixed partials no longer commute.
  • For systems needing more variables — a magnetic system adds B and M, and the corresponding relations must be derived afresh.
  • Out of equilibrium, where the potentials are not defined.
Failure modes
  • Omitting the subscript. (∂S/∂V) without specifying what is held constant is not a defined quantity, and dropping it produces confident nonsense.
  • Memorising the four relations rather than the two-line mechanism. The mechanism generalises to any system; the memorised set does not.
  • Sign errors from the transform conventions. Derive each d-expression rather than recalling it.
Worked number

From the third relation, for an ideal gas (∂P/∂T)V = N kB/V, hence (∂S/∂V)T = N kB/V. Integrating at fixed T gives ΔS = N kB ln(V2/V1) for isothermal expansion — recovered from an equation of state alone, with no microscopic model. Combined with D-019 it delivers the Sackur–Tetrode result of D-017.

Run the check → See the same transform in D-004 →

Discussion

The content of a Maxwell relation is almost embarrassingly thin at its source and enormous in its consequences. Step 3 of the derivation is Clairaut's theorem — a statement about smooth functions of two variables that has nothing to do with heat, work or entropy. All the physics was spent at step 1, in the single assertion that U is a state function whose differential is dU = TdSPdV. Everything after that is bookkeeping. The remarkable thing is that this bookkeeping equates a caloric quantity (a derivative of entropy or temperature) to a mechanical one (a derivative of pressure or volume). That is why the two lower relations are the workhorses: entropy has no gauge, but (∂P/∂T)V and (∂V/∂T)P are read off a pressure sensor and a dilatometer.

The four relations are structured by the four potentials, and the potentials are structured by which pair of natural variables you choose. This is not four independent facts — it is one fact (the exactness of a differential) applied to the four Legendre transforms of a single function. The pattern of signs is dictated entirely by the transform: a variable appears with a plus sign when it is the “slope” conjugate to what you differentiate, and the minus signs migrate exactly as they do when you pass from U to F = UTS. If you understand why dF = −SdTPdV has the signs it does, you can regenerate the whole set on demand and never memorise a thing.

The deeper structure is geometric. A thermodynamic state is a point on a two-dimensional surface embedded in the higher-dimensional space of all the variables, and the potentials are generating functions for that surface in the sense of contact geometry. The Maxwell relations are precisely the integrability (closedness) conditions d(dU) = 0 for the one-forms involved — the statement that the forms TdSPdV etc. are exact rather than merely closed. Seen this way they are the thermodynamic analogue of the equality of mixed partials that gives you conservative force fields (∇×F = 0) in mechanics, and of the Onsager reciprocal relations in linear irreversible thermodynamics, which are their non-equilibrium descendants. The same skeleton reappears whenever a system is governed by a potential: it is exactly the mechanism that made D-004 produce symmetric Hessians in Lagrangian mechanics.

Common misconceptions. (i) Students often read (∂S/∂V)T = (∂P/∂T)V as a coincidence to be verified case by case; it is an identity forced on every substance by the existence of F, and if a proposed equation of state violated it, that equation of state would be thermodynamically impossible. (ii) The subscript is not optional shorthand: (∂S/∂V)T and (∂S/∂V)P are genuinely different numbers, and the relation holds only for the specific pairing shown. (iii) The relations are statements about equilibrium state functions, so “quasi-static” in the assumptions governs the derivation route (the TdS form), not the validity of the result — once derived, they constrain the state functions of the material regardless of how any particular process is run.

Worked examples

Example 1 — Internal pressure of a real gas. How much does the internal energy of one mole of CO₂ change when it expands isothermally from 1.0 L to 10 L? An ideal gas would give zero; the Maxwell relations tell us exactly what the intermolecular attraction contributes.

1
( ∂U/∂V )T = T( ∂S/∂V )TP
Divide dU = TdSPdV through by dV at fixed T.
2
( ∂U/∂V )T = T( ∂P/∂T )VP
Substitute the third Maxwell relation (∂S/∂V)T = (∂P/∂T)V. This is the thermodynamic equation of state — internal pressure from the equation of state alone.
3
P = nRT/(V − nb) − an²/V²  ⇒  ( ∂P/∂T )V = nR/(V − nb)
Van der Waals equation of state; only the first term carries T.
4
( ∂U/∂V )T = T·nR/(V−nb) − [ nRT/(V−nb) − an²/V² ] = an²/V²
The repulsive (b) term cancels exactly; the internal pressure is purely the attractive term.
5
ΔU = ∫VV an²/V² dV = an²( 1/V₁ − 1/V₂ )
Integrate at fixed T. Now insert numbers: n = 1 mol, a = 0.364 Pa·m⁶·mol⁻², V₁ = 1.0×10⁻³ m³, V₂ = 1.0×10⁻² m³.
6
ΔU = 0.364 × (1000 − 100) m⁻³ · Pa·m⁶ = 0.364 × 900 J
Arithmetic. Units: Pa·m³ = J.
ΔU ≈ +3.3 × 10² J

Answer. The internal energy rises by about 328 J per mole. It is positive because pulling attracting molecules apart stores energy — energy the ideal-gas model misses entirely. The number came from the equation of state alone, with no microscopic potential assumed, exactly as the Maxwell machinery promises.

Example 2 — Adiabatic heating of water under compression. A sealed, thermally insulated sample of liquid water at 298 K is compressed by 100 bar. By how much does its temperature rise? There is no heat flow, yet the temperature changes — and a Maxwell relation gives the coefficient.

1
dS = ( ∂S/∂T )P dT + ( ∂S/∂P )T dP = 0
Reversible adiabatic = isentropic; write S as a function of T,P and set dS = 0.
2
( ∂S/∂T )P = CP/T ,   ( ∂S/∂P )T = −( ∂V/∂T )P = −Vα
The first from the definition of CP; the second is the fourth Maxwell relation, then the definition of the isobaric expansivity α ≡ (1/V)(∂V/∂T)P.
3
( ∂T/∂P )S = TVα / CP
Solve step 1 for dT/dP. Symbolic result — a pure combination of measurable coefficients.
4
( ∂T/∂P )S = (298 K)(18.0×10⁻⁶ m³/mol)(2.57×10⁻⁴ K⁻¹) / (75.4 J mol⁻¹ K⁻¹)
Molar values for water at 298 K, kept molar throughout so V and CP are consistent.
5
( ∂T/∂P )S = 1.38×10⁻⁶ / 75.4 ≈ 1.83×10⁻⁸ K/Pa
Evaluate. Units: (K · m³ · K⁻¹) / (J K⁻¹) = m³·K/J = K/Pa.
6
ΔT ≈ ( ∂T/∂P )S · ΔP = (1.83×10⁻⁸ K/Pa)(1.0×10⁵ Pa)
100 bar = 1.0×10⁵ Pa. Coefficient is nearly constant over this range, so a linear estimate is fine.
ΔT ≈ +0.18 K

Answer. The water warms by about 0.18 K. It is small because water is nearly incompressible and its expansivity is modest, but it is real, sign-definite (compression of a normally-expanding material heats it), and computed entirely from α, CP and V — no calorimetry of the compression itself required. Note that near 4 °C water has α < 0, and there the same formula predicts adiabatic cooling under compression.

Problems
  1. (Warm-up.) Starting from the Gibbs free energy G = U + PVTS, write down its differential, identify its natural variables, and state the Maxwell relation it yields.
    Solution

    dG = dU + PdV + VdPTdSSdT. Substitute dU = TdSPdV; the TdS and PdV terms cancel, leaving dG = −SdT + VdP. Natural variables are T and P, so −S = (∂G/∂T)P and V = (∂G/∂P)T. Equality of mixed partials gives (∂S/∂P)T = −(∂V/∂T)P — the fourth relation.

  2. (Ideal gas.) Show that the internal energy of an ideal gas depends only on temperature, i.e. (∂U/∂V)T = 0.
    Solution

    From the thermodynamic equation of state (∂U/∂V)T = T(∂P/∂T)VP, derived by dividing dU = TdSPdV at fixed T and inserting the third Maxwell relation. For an ideal gas P = NkBT/V, so (∂P/∂T)V = NkB/V = P/T. Hence T(∂P/∂T)VP = T·P/TP = 0. So U = U(T) only — Joule's law, obtained without any kinetic-theory model.

  3. (Heat-capacity gap.) Prove the general relation CPCV = TVα²/κT, where α is the isobaric expansivity and κT the isothermal compressibility, and evaluate it numerically for water at 298 K (α = 2.57×10⁻⁴ K⁻¹, κT = 4.52×10⁻¹⁰ Pa⁻¹, Vm = 18.0×10⁻⁶ m³/mol).
    Solution

    Standard identity: CPCV = T(∂P/∂T)V(∂V/∂T)P. By the triple product rule (∂P/∂T)V = −(∂V/∂T)P/(∂V/∂P)T = Vα/(VκT) = α/κT, and (∂V/∂T)P = Vα. Multiplying: CPCV = T(α/κT)(Vα) = TVα²/κT. (The Maxwell relations underlie the entropy identities that make this reduction possible; the result is manifestly ≥ 0 since V, T, κT > 0.) Numerically: TVα²/κT = 298 × 18.0×10⁻⁶ × (2.57×10⁻⁴)² / 4.52×10⁻¹⁰ = 298 × 18.0×10⁻⁶ × 6.60×10⁻⁸ / 4.52×10⁻¹⁰ ≈ 0.78 J mol⁻¹ K⁻¹. Far below the ideal-gas value R = 8.31, reflecting how little water expands.

  4. (Joule–Thomson.) Derive the Joule–Thomson coefficient μJT = (∂T/∂P)H in terms of measurable quantities, and show it vanishes for an ideal gas. What is the sign condition for cooling on expansion?
    Solution

    At constant H, μJT = (∂T/∂P)H = −(∂H/∂P)T/(∂H/∂T)P. The denominator is CP. For the numerator use dH = TdS + VdP: (∂H/∂P)T = T(∂S/∂P)T + V. Apply the fourth Maxwell relation (∂S/∂P)T = −(∂V/∂T)P, giving (∂H/∂P)T = VT(∂V/∂T)P. Therefore μJT = [T(∂V/∂T)PV]/CP = (V/CP)(Tα − 1). For an ideal gas V = NkBT/P, so T(∂V/∂T)P = V and μJT = 0 — no throttling temperature change. A real gas cools on expansion (μJT > 0) when Tα > 1, i.e. below its inversion temperature; above it, throttling warms the gas.

  5. (Second derivatives.) Using a Maxwell relation, show that (∂CV/∂V)T = T(∂²P/∂T²)V. Hence determine whether CV of a van der Waals gas depends on volume.
    Solution

    Write CV = T(∂S/∂T)V. Then (∂CV/∂V)T = T ∂/∂V[(∂S/∂T)V]T = T ∂/∂T[(∂S/∂V)T]V, where the mixed partials of S have been swapped (Clairaut — the same move as step 3 of the derivation). Insert the third Maxwell relation (∂S/∂V)T = (∂P/∂T)V, giving (∂CV/∂V)T = T(∂²P/∂T²)V. For a van der Waals gas P = nRT/(Vnb) − an²/V² is linear in T at fixed V, so (∂²P/∂T²)V = 0 and CV is independent of volume — just as for an ideal gas. Volume-dependence of CV would require an equation of state with a nonlinear temperature dependence at fixed volume; the whole result follows from one commuting pair of second derivatives.