Physics Thermal Physics 100% Free Open Access
Chapter 5 • Theory & Derivations

Maxwell’s Thermodynamic Relations

Cross-derivative Maxwell identities, Clausius-Clapeyron equation, rigorous heat capacity differences Cp - Cv, the two TdS equations, internal energy equations of state, and Joule-Kelvin throttling coefficients.

§5.1 Derivation of the Four Maxwell Relations

1. Mathematical Foundation of Maxwell's Relations

James Clerk Maxwell (1871) realized that because the four thermodynamic potentials ($U, H, F, G$) are exact state functions, their mixed second partial derivatives must commute by Schwarz's theorem:

$$\frac{\partial^2 \Phi}{\partial x \partial y} = \frac{\partial^2 \Phi}{\partial y \partial x}$$

2. Systematic Derivations from Fundamental Differentials

1. First Maxwell Relation (from Internal Energy $U(S, V)$):

$$dU = T dS - P dV$$

Since $dU$ is exact:

$$\left(\frac{\partial U}{\partial S}\right)_V = T, \quad \left(\frac{\partial U}{\partial V}\right)_S = -P$$

Equating mixed second derivatives $\frac{\partial^2 U}{\partial V \partial S} = \frac{\partial^2 U}{\partial S \partial V}$:

$$\left(\frac{\partial T}{\partial V}\right)_S = -\left(\frac{\partial P}{\partial S}\right)_V$$

2. Second Maxwell Relation (from Enthalpy $H(S, P)$):

$$dH = T dS + V dP$$
$$\left(\frac{\partial H}{\partial S}\right)_P = T, \quad \left(\frac{\partial H}{\partial P}\right)_S = V$$

Equating $\frac{\partial^2 H}{\partial P \partial S} = \frac{\partial^2 H}{\partial S \partial P}$:

$$\left(\frac{\partial T}{\partial P}\right)_S = \left(\frac{\partial V}{\partial S}\right)_P$$

3. Third Maxwell Relation (from Helmholtz Free Energy $F(T, V)$):

$$dF = -S dT - P dV$$
$$\left(\frac{\partial F}{\partial T}\right)_V = -S, \quad \left(\frac{\partial F}{\partial V}\right)_T = -P$$

Equating $\frac{\partial^2 F}{\partial V \partial T} = \frac{\partial^2 F}{\partial T \partial V}$:

$$\left(\frac{\partial S}{\partial V}\right)_T = \left(\frac{\partial P}{\partial T}\right)_V$$

4. Fourth Maxwell Relation (from Gibbs Free Energy $G(T, P)$):

$$dG = -S dT + V dP$$
$$\left(\frac{\partial G}{\partial T}\right)_P = -S, \quad \left(\frac{\partial G}{\partial P}\right)_T = V$$

Equating $\frac{\partial^2 G}{\partial P \partial T} = \frac{\partial^2 G}{\partial T \partial P}$:

$$\left(\frac{\partial S}{\partial P}\right)_T = -\left(\frac{\partial V}{\partial T}\right)_P$$

3. Max Born's Mnemonic Square

To quickly reconstruct these relations, Max Born proposed the mnemonic square:

$$\begin{array}{ccc} -V & & P \\ F & & H \\ -T & & S \end{array} \quad \text{or mnemonic: "Good Physicists Have Studied Under Very Fine Teachers"}$$

The potentials flank the edges between their natural variables. Arrow directions encode negative signs.

§5.2 The Clausius-Clapeyron Equation

1. Phase Coexistence Equilibrium

Consider two distinct macroscopic phases of a pure substance (e.g., liquid 1 and vapor 2) coexisting in dynamic phase equilibrium at temperature $T$ and pressure $P$. Equilibrium requires equality of their specific Gibbs free energies (chemical potentials):

$$g_1(T, P) = g_2(T, P)$$

If temperature shifts by an infinitesimal $dT$, the equilibrium coexistence pressure must shift along the phase line by $dP$ such that:

$$g_1(T + dT, P + dP) = g_2(T + dT, P + dP)$$

Expanding in total differentials $dg = -s dT + v dP$:

$$-s_1 dT + v_1 dP = -s_2 dT + v_2 dP$$

Rearranging terms:

$$(v_2 - v_1) dP = (s_2 - s_1) dT \implies \frac{dP}{dT} = \frac{s_2 - s_1}{v_2 - v_1} = \frac{\Delta s}{\Delta v}$$

2. The Latent Heat Formula

The entropy change during a reversible phase transformation at constant temperature $T$ is related to the latent heat $L$ per mole by $\Delta s = L / T$:

$$\frac{dP}{dT} = \frac{L}{T (v_2 - v_1)}$$

This is the Clausius-Clapeyron Equation.

3. Physical Applications

1. Boiling Point Elevation with Pressure: For vaporization, $v_{\text{vapor}} \gg v_{\text{liquid}} \implies \Delta v > 0$. Since vaporization is endothermic ($L > 0$), $\frac{dP}{dT} > 0$. Increasing pressure elevates the boiling point (utilized in pressure cookers and autoclave sterilization).

2. Depression of Freezing Point of Water: For water, ice has a lower density than liquid water at $0^\circ\text{C}$ ($v_{\text{water}} < v_{\text{ice}} \implies \Delta v = v_{\text{water}} - v_{\text{ice}} < 0$). Because latent heat of fusion $L_f > 0$:

$$\frac{dP}{dT} = \frac{L_f}{T (v_{\text{water}} - v_{\text{ice}})} < 0$$

Applying pressure lowers the melting temperature of ice (facilitating glacier regelation and ice skating).

§5.3 Specific Heat Relations (Cp - Cv) & The Two TdS Equations

1. Derivation of the First TdS Equation

Let entropy be expressed as a function of temperature and volume, $S = S(T, V)$:

$$dS = \left(\frac{\partial S}{\partial T}\right)_V dT + \left(\frac{\partial S}{\partial V}\right)_T dV$$

Multiplying by absolute temperature $T$:

$$T dS = T \left(\frac{\partial S}{\partial T}\right)_V dT + T \left(\frac{\partial S}{\partial V}\right)_T dV$$

Using the definition of isochoric heat capacity $C_v = T (\partial S/\partial T)_V$ and the Third Maxwell Relation $(\partial S/\partial V)_T = (\partial P/\partial T)_V$:

$$T dS = C_v dT + T \left(\frac{\partial P}{\partial T}\right)_V dV$$

This is the First $T dS$ Equation.

2. Derivation of the Second TdS Equation

Now let entropy be expressed as a function of temperature and pressure, $S = S(T, P)$:

$$dS = \left(\frac{\partial S}{\partial T}\right)_P dT + \left(\frac{\partial S}{\partial P}\right)_T dP$$

Multiplying by $T$:

$$T dS = T \left(\frac{\partial S}{\partial T}\right)_P dT + T \left(\frac{\partial S}{\partial P}\right)_T dP$$

Using isobaric heat capacity $C_p = T (\partial S/\partial T)_P$ and the Fourth Maxwell Relation $(\partial S/\partial P)_T = -(\partial V/\partial T)_P$:

$$T dS = C_p dT - T \left(\frac{\partial V}{\partial T}\right)_P dP$$

This is the Second $T dS$ Equation.

3. Rigorous Derivation of $C_p - C_v$

Equating the First and Second $T dS$ equations:

$$C_p dT - T \left(\frac{\partial V}{\partial T}\right)_P dP = C_v dT + T \left(\frac{\partial P}{\partial T}\right)_V dV$$
$$(C_p - C_v) dT = T \left(\frac{\partial P}{\partial T}\right)_V dV + T \left(\frac{\partial V}{\partial T}\right)_P dP$$

Dividing by $dT$ at constant pressure ($dP = 0$):

$$C_p - C_v = T \left(\frac{\partial P}{\partial T}\right)_V \left(\frac{\partial V}{\partial T}\right)_P$$

Expressing in terms of thermal expansion coefficient $\alpha = \frac{1}{V}(\partial V/\partial T)_P$ and isothermal compressibility $\kappa_T = -\frac{1}{V}(\partial V/\partial P)_T$, with the cyclic identity $(\partial P/\partial T)_V = \alpha / \kappa_T$:

$$C_p - C_v = T V \frac{\alpha^2}{\kappa_T}$$

§5.4 The Energy Equations & Internal Pressure

1. The First Energy Equation (Volume Dependence of U)

From the First Law, $dU = T dS - P dV$. Dividing by $dV$ at constant temperature $T$:

$$\left(\frac{\partial U}{\partial V}\right)_T = T \left(\frac{\partial S}{\partial V}\right)_T - P$$

Substituting the Third Maxwell Relation $(\partial S/\partial V)_T = (\partial P/\partial T)_V$:

$$\left(\frac{\partial U}{\partial V}\right)_T = T \left(\frac{\partial P}{\partial T}\right)_V - P$$

This is the First Energy Equation. The term $(\partial U/\partial V)_T$ is called the internal pressure ($P_i$), representing the internal force per unit area exerted by intermolecular attractive forces.

Examples:

1. Ideal Gas: $P = nRT / V \implies (\partial P/\partial T)_V = nR/V$.

$$\left(\frac{\partial U}{\partial V}\right)_T = T \left(\frac{n R}{V}\right) - P = P - P = 0$$

proving Joule's law theoretically without approximations.

2. Van der Waals Gas: $P = \frac{n R T}{V - n b} - \frac{n^2 a}{V^2} \implies (\partial P/\partial T)_V = \frac{n R}{V - n b}$.

$$\left(\frac{\partial U}{\partial V}\right)_T = T \left(\frac{n R}{V - n b}\right) - \left(\frac{n R T}{V - n b} - \frac{n^2 a}{V^2}\right) = \frac{n^2 a}{V^2}$$

Internal energy increases during isothermal expansion as intermolecular bonds stretch.

2. The Second Energy Equation (Pressure Dependence of H)

From enthalpy $dH = T dS + V dP$, dividing by $dP$ at constant temperature $T$:

$$\left(\frac{\partial H}{\partial P}\right)_T = T \left(\frac{\partial S}{\partial P}\right)_T + V$$

Substituting the Fourth Maxwell Relation $(\partial S/\partial P)_T = -(\partial V/\partial T)_P$:

$$\left(\frac{\partial H}{\partial P}\right)_T = V - T \left(\frac{\partial V}{\partial T}\right)_P$$

This is the Second Energy Equation. For an ideal gas, $(\partial V/\partial T)_P = V/T$, so $(\partial H/\partial P)_T = 0$.

§5.5 Joule-Kelvin Throttling & Inversion Temperature

1. The Porous Plug Throttling Experiment

In the Joule-Thomson (Joule-Kelvin) experiment (1852), a continuous gas stream under steady high pressure $P_1$ is forced through an insulated porous plug or throttling valve to a lower pressure $P_2$.

  • Upstream work done on the gas by piston: $W_1 = -P_1 V_1$.
  • Downstream work done by the gas on piston: $W_2 = +P_2 V_2$.
  • Net work done by gas: $W = P_2 V_2 - P_1 V_1$.
  • Because the tube is thermally insulated, $Q = 0$.

By the First Law:

$$\Delta U = Q - W \implies U_2 - U_1 = -(P_2 V_2 - P_1 V_1) \implies U_1 + P_1 V_1 = U_2 + P_2 V_2$$
$$H_1 = H_2$$

A throttling process is strictly isenthalpic ($\Delta H = 0$).

2. Derivation of the Joule-Thomson Coefficient

The Joule-Thomson Coefficient $\mu_{\text{JT}}$ is the rate of temperature change with pressure under constant enthalpy:

$$\mu_{\text{JT}} \equiv \left(\frac{\partial T}{\partial P}\right)_H$$

By the cyclic permutation identity for $(T, P, H)$:

$$\left(\frac{\partial T}{\partial P}\right)_H = -\frac{(\partial H/\partial P)_T}{(\partial H/\partial T)_P} = -\frac{1}{C_p} \left(\frac{\partial H}{\partial P}\right)_T$$

Substituting the Second Energy Equation $(\partial H/\partial P)_T = V - T(\partial V/\partial T)_P$:

$$\mu_{\text{JT}} = \frac{1}{C_p} \left[ T \left(\frac{\partial V}{\partial T}\right)_P - V \right]$$

3. Joule-Thomson Cooling for Real Gases

1. Ideal Gas: $(\partial V/\partial T)_P = V/T \implies \mu_{\text{JT}} = 0$ (no temperature change).

2. Van der Waals Gas: Using $(P + a/V_m^2)(V_m - b) \approx RT$, we find:

$$\mu_{\text{JT}} \approx \frac{1}{C_p} \left[ \frac{2a}{R T} - b \right]$$
  • Cooling Zone ($\mu_{\text{JT}} > 0$): Since $dP < 0$ across the throttle, $dT = \mu_{\text{JT}} dP < 0$. Gas cools if $T < T_i$.
  • Heating Zone ($\mu_{\text{JT}} < 0$): Gas heats up upon throttling if $T > T_i$.
  • Temperature of Inversion $T_i$: The boundary where $\mu_{\text{JT}} = 0$:
$$T_i = \frac{2a}{R b}$$

For air, $T_i \approx 600\text{ K}$ (cools at room temperature). For hydrogen ($T_i \approx 200\text{ K}$) and helium ($T_i \approx 40\text{ K}$), pre-cooling below their inversion temperatures is mandatory before throttling can produce liquefaction.

Standard University Exam Examination Problems

Rigorous Analytical & Numerical Solved Problems

Comprehensive step-by-step mathematical proofs, dimensional evaluations, and calculations matching B.Sc. Honors university examinations.

Standard Example 5.1: Clausius-Clapeyron Vapor Pressure & Boiling Point Elevation

At standard atmospheric pressure $P_0 = 1.013\times 10^5\text{ Pa}$, water boils at $T_0 = 373.15\text{ K}$ ($100.0^\circ\text{C}$) with a specific latent heat of vaporization of $L_v = 2.257\times 10^6\text{ J/kg}$. The specific volume of steam is $v_{\text{steam}} = 1.673\text{ m}^3\text{/kg}$, while the specific volume of liquid water is $v_{\text{water}} = 1.043\times 10^{-3\text{ m}^3\text{/kg}$. (a) Calculate the rate of change of boiling point with pressure $dT/dP$ at $100^\circ\text{C}$. (b) If atmospheric pressure at a high-altitude mountain station is $P = 70.0\text{ kPa}$, estimate the boiling temperature of water at that station using the integrated Clausius-Clapeyron equation.

Step 1: Calculate $dT/dP$ at $100^\circ\text{C}$

From the Clausius-Clapeyron equation:

$$\frac{dP}{dT} = \frac{L_v}{T (v_{\text{steam}} - v_{\text{water}})}$$
$$\Delta v = 1.673 - 0.001043 = 1.67196\text{ m}^3\text{/kg}$$
$$\frac{dP}{dT} = \frac{2.257 \times 10^6\text{ J/kg}}{(373.15\text{ K})(1.67196\text{ m}^3\text{/kg})} = \frac{2.257 \times 10^6}{623.89} = 3617.6\text{ Pa/K}$$

Inverting to obtain $dT/dP$:

$$\frac{dT}{dP} = \frac{1}{3617.6} = 2.764 \times 10^{-4}\text{ K/Pa} = 0.280\text{ K/torr} = 27.6\text{ mK/kPa}$$

Every $1\text{ kPa}$ drop in pressure lowers the boiling point of water by $0.0276^\circ\text{C}$.

Step 2: Integrated Clausius-Clapeyron Equation

Approximating steam as an ideal gas ($v_{\text{steam}} \approx R T / P M$) and neglecting $v_{\text{liquid}}$:

$$\frac{d \ln P}{dT} = \frac{L_{v, m}}{R T^2} \implies \ln\left(\frac{P}{P_0}\right) = -\frac{L_{v, m}}{R} \left( \frac{1}{T} - \frac{1}{T_0} \right)$$

Molar latent heat: $L_{v, m} = L_v \times M = (2.257 \times 10^6)(0.018015\text{ kg/mol}) = 40660\text{ J/mol}$.

$$\ln\left(\frac{70.0}{101.3}\right) = \ln(0.6910) = -0.3696$$
$$\frac{L_{v, m}}{R} = \frac{40660}{8.314} = 4890.5\text{ K}$$
$$-0.3696 = -4890.5 \left( \frac{1}{T} - \frac{1}{373.15} \right)$$
$$\frac{1}{T} - 0.0026799 = \frac{0.3696}{4890.5} = 7.5575 \times 10^{-5}$$
$$\frac{1}{T} = 0.0026799 + 0.00007558 = 0.0027555\text{ K}^{-1}$$
$$T = \frac{1}{0.0027555} = 362.91\text{ K} = 89.76^\circ\text{C}$$
Final Answer & Physical Insight

(a) $dT/dP = 2.76\times 10^{-4}\text{ K/Pa} = 27.6\text{ mK/kPa}$. (b) At $70.0\text{ kPa}$, water boils at $T = 362.9\text{ K}$ ($89.8^\circ\text{C}$), a drop of over $10^\circ\text{C}$ below standard sea-level boiling.

Standard Example 5.2: Adiabatic Compression Temperature Rise of Liquid Water

Liquid water at initial temperature $T_1 = 293.15\text{ K}$ ($20.0^\circ\text{C}$) and atmospheric pressure $P_1 = 1.0\text{ bar}$ ($1.0\times 10^5\text{ Pa}$) is compressed reversibly and adiabatically to a final pressure of $P_2 = 1000.0\text{ bar}$ ($1.0\times 10^8\text{ Pa}$). Over this pressure range, the average density is $\rho = 1010\text{ kg/m}^3$, isobaric specific heat capacity is $c_p = 4180\text{ J/(kg}\cdot\text{K)}$, and thermal expansion coefficient is $\alpha = 2.10\times 10^{-4}\text{ K}^{-1}$. Derive the adiabatic heating formula $\left(\frac{\partial T}{\partial P}\right)_S = \frac{T v \alpha}{c_p}$ and calculate the final water temperature $T_2$.

Step 1: Derive the Adiabatic Temperature-Pressure Gradient

From the Second $T dS$ equation:

$$T dS = c_p dT - T \left(\frac{\partial v}{\partial T}\right)_P dP$$

For a reversible adiabatic process, $dS = 0$:

$$c_p dT = T \left(\frac{\partial v}{\partial T}\right)_P dP$$

Recalling the definition of thermal expansion coefficient $\alpha = \frac{1}{v} \left(\frac{\partial v}{\partial T}\right)_P \implies \left(\frac{\partial v}{\partial T}\right)_P = v \alpha$:

$$\left(\frac{\partial T}{\partial P}\right)_S = \frac{T v \alpha}{c_p} = \frac{T \alpha}{\rho c_p}$$

This confirms that compressing any substance with positive thermal expansion ($\alpha > 0$) adiabatically raises its temperature.

Step 2: Integrate to Find the Temperature Rise $\Delta T$

Assuming $T \approx T_1$ in the derivative because $\Delta T$ is small relative to $T_1$:

$$\Delta T = \int_{P_1}^{P_2} \frac{T_1 \alpha}{\rho c_p} dP = \frac{T_1 \alpha}{\rho c_p} (P_2 - P_1)$$

Given values:

$$T_1 = 293.15\text{ K}, \quad \alpha = 2.10 \times 10^{-4}\text{ K}^{-1}$$
$$\rho = 1010\text{ kg/m}^3, \quad c_p = 4180\text{ J/(kg}\cdot\text{K)}$$
$$\Delta P = 1.0 \times 10^8 - 1.0 \times 10^5 \approx 1.0 \times 10^8\text{ Pa}$$
$$\rho c_p = (1010)(4180) = 4.2218 \times 10^6\text{ J/(m}^3\cdot\text{K)}$$
$$\Delta T = \frac{(293.15)(2.10 \times 10^{-4})}{4.2218 \times 10^6} \times 1.0 \times 10^8 = \frac{0.06156}{4.2218 \times 10^6} \times 1.0 \times 10^8 = (1.458 \times 10^{-8})(10^8) = +1.46\text{ K}$$
Step 3: Final Temperature

Final temperature:

$$T_2 = T_1 + \Delta T = 293.15 + 1.46 = 294.61\text{ K} = 21.46^\circ\text{C}$$

Despite an enormous pressure increase of $1000\text{ bar}$, water heats by only $1.46^\circ\text{C}$ due to its high density, large heat capacity, and low thermal expansion.

Final Answer & Physical Insight

Formula derived: $\left(\frac{\partial T}{\partial P}\right)_S = \frac{T v \alpha}{c_p}$. Temperature increase is $\Delta T = +1.46\text{ K}$, yielding final temperature $T_2 = 294.61\text{ K}$ ($21.46^\circ\text{C}$).

Standard Example 5.3: Joule-Thomson Inversion Curve & Cooling Calculation for Nitrogen

For gaseous nitrogen ($\text{N}_2$), the Van der Waals constants are $a = 0.137\text{ J}\cdot\text{m}^3\text{/mol}^2$ and $b = 3.87\times 10^{-5}\text{ m}^3\text{/mol}$, and the molar heat capacity at constant pressure is $C_p = 29.12\text{ J/(mol}\cdot\text{K)}$. (a) Calculate the maximum inversion temperature $T_i$ for nitrogen. (b) Calculate the Joule-Thomson coefficient $\mu_{\text{JT}}$ at $T = 300.0\text{ K}$. (c) If nitrogen at $300.0\text{ K}$ is throttled from $P_1 = 150.0\text{ bar}$ to $P_2 = 1.0\text{ bar}$, calculate the temperature drop $\Delta T$ across the porous plug.

Step 1: Calculate Maximum Inversion Temperature $T_i$

For a Van der Waals gas, the inversion temperature at zero pressure is:

$$T_i = \frac{2a}{R b} = \frac{2(0.137\text{ J}\cdot\text{m}^3\text{/mol}^2)}{(8.314\text{ J/(mol}\cdot\text{K)})(3.87 \times 10^{-5}\text{ m}^3\text{/mol})}$$
$$T_i = \frac{0.274}{3.2175 \times 10^{-4}} = 851.6\text{ K} = 578.4^\circ\text{C}$$

Because room temperature ($300\text{ K}$) is well below $851.6\text{ K}$, nitrogen will cool upon throttling.

Step 2: Calculate Joule-Thomson Coefficient $\mu_{\text{JT}}$ at $300\text{ K}$

The theoretical Joule-Thomson coefficient for a Van der Waals gas is:

$$\mu_{\text{JT}} = \frac{1}{C_p} \left[ \frac{2a}{R T} - b \right]$$
$$\frac{2a}{R T} = \frac{2(0.137)}{(8.314)(300.0)} = \frac{0.274}{2494.2} = 1.0985 \times 10^{-4}\text{ m}^3\text{/mol}$$
$$\frac{2a}{R T} - b = 1.0985 \times 10^{-4} - 0.387 \times 10^{-4} = 7.115 \times 10^{-5}\text{ m}^3\text{/mol}$$
$$\mu_{\text{JT}} = \frac{7.115 \times 10^{-5}\text{ m}^3\text{/mol}}{29.12\text{ J/(mol}\cdot\text{K)}} = 2.443 \times 10^{-6}\text{ K/Pa} = 0.2443\text{ K/bar}$$

Each bar of pressure reduction across the throttle drops the nitrogen temperature by $0.244^\circ\text{C}$.

Step 3: Calculate Throttling Temperature Drop $\Delta T$

The pressure drop across the porous plug is:

$$\Delta P = P_2 - P_1 = 1.0 - 150.0 = -149.0\text{ bar}$$
$$\Delta T = \mu_{\text{JT}} \Delta P = (0.2443\text{ K/bar})(-149.0\text{ bar}) = -36.40\text{ K}$$

Final temperature downstream of throttle:

$$T_2 = T_1 + \Delta T = 300.0 - 36.4 = 263.6\text{ K} = -9.55^\circ\text{C}$$

This significant cooling illustrates the physical basis for industrial regenerative air liquefaction (the Linde-Hampson process).

Final Answer & Physical Insight

(a) Maximum inversion temperature $T_i = 851.6\text{ K}$ ($578.4^\circ\text{C}$). (b) At $300\text{ K}$, $\mu_{\text{JT}} = 0.244\text{ K/bar}$. (c) Throttling across $149\text{ bar}$ produces a cooling of $\Delta T = -36.4\text{ K}$, reaching $263.6\text{ K}$ ($-9.6^\circ\text{C}$).