Heat capacity and thermal transport diagnostics
Distinguish thermodynamic state counting from heat-current transport and test the assumptions behind low-temperature fits.
On this page
The question for today
What you’ll understand
- Extract electronic and phonon coefficients from a conditional low-temperature heat-capacity model.
- Convert a measured thermal conductance into a conductivity using the correct geometry.
- Apply a Wiedemann–Franz benchmark without assigning a unique carrier mechanism from one ratio.
Heat capacity asks how much additional energy a sample stores when its temperature rises. Thermal conductivity asks how efficiently excitations carry heat through it. A large density of low-energy states can increase the first while strong scattering suppresses the second. Comparing these observables with electrical transport is therefore informative precisely because they weight states and scattering differently. Neither a fitted linear heat-capacity term nor a residual thermal signal uniquely identifies one exotic excitation without additional controls and a justified extrapolation.
A low-temperature fit is a physical hypothesis
For a conventional three-dimensional metal at temperatures well below relevant electronic and Debye scales, a useful molar model combines a linear electronic heat capacity with a cubic acoustic-phonon term. Plotting heat capacity divided by temperature against temperature squared separates intercept and slope. The electronic coefficient reflects a thermodynamic density of states, including appropriate quasiparticle renormalization in a Fermi-liquid interpretation. It is not the same weighted quantity as a Hall number or one photoemission intensity.
Magnetic excitations, nuclear or impurity Schottky anomalies, low-dimensional modes, glassy defects, and additional transitions can invalidate that two-term fit. Addenda from the platform and grease must be subtracted, the sample amount must be known, and thermal equilibration must be demonstrated. Experiments often measure heat capacity at constant pressure; comparison with a constant-volume model requires that the thermal-expansion correction be negligible or included. Changing the fit window is a physical stability test, not merely a way to improve a regression coefficient.
A thermal current needs a temperature profile
In a steady one-dimensional diffusive bar, heat current density is minus thermal conductivity times temperature gradient. A heater power and a measured temperature difference give conductance, while conductivity additionally requires length and cross section. The power entering the sample can differ from electrical heater power because of heat leaks through wires, radiation, and the mounting structure. A contact temperature drop can also make the apparent sample gradient too large and the inferred conductivity too small.
Use thermometers calibrated in the same field and temperature regime as the measurement, check linearity with heater power, and allow a stationary profile to develop. At very low temperatures the electronic and lattice subsystems or the thermometers can decouple, so the reported temperature may not represent every heat-carrying subsystem. Boundary-limited or ballistic transport can require a description beyond a bulk local conductivity. In a magnetic field, transverse heat flow and thermoelectric corrections can require a tensor treatment rather than a scalar bar formula.
Charge and heat need not be carried in the same way
For degenerate conventional charged quasiparticles dominated by elastic scattering, a low-temperature electronic Wiedemann–Franz benchmark relates thermal and electrical conductivity through the Sommerfeld Lorenz number. The thermal conductivity is the open-circuit electronic value; phonons and neutral excitations add heat without carrying the same electrical current. Inelastic scattering, multiband structure, finite-temperature thermoelectric effects, and nonlocal transport can alter a naive comparison. The benchmark is not an unrestricted identity for all quantum matter.
A residual intercept of thermal conductivity divided by temperature depends on how phonons and other contributions are extrapolated. A fitted intercept can change with boundary scattering or a mistaken temperature calibration, while a genuine electronic intercept need not determine a unique gap structure. Combine temperature and field trends with calorimetric entropy and independent gap probes. Shared geometry errors can partly cancel in a Lorenz ratio when exactly the same segment carries both currents, but mismatched contacts and different heat paths can destroy that cancellation.
Work it through
Two calorimetric points determine a conditional fit
A hypothetical sample follows the two-term molar model. Its is at 1 K and 7 in the same units at 2 K. Find gamma, beta, and the predicted heat capacity at 3 K.
The horizontal coordinates T squared are 1 and 4 square kelvin. The slope is , which is beta.
Subtract the phonon contribution at 1 K: . Two points determine the model coefficients but cannot test whether extra terms are absent.
At 3 K the electronic contribution is and the phonon contribution is . Their sum is .
The conditional coefficients are gamma three and beta one in the stated units, predicting 36 millijoules per mole kelvin at 3 kelvin.
The living notebook
Explore the model
Interactive lab · heat
Diffusion depends on scale as well as diffusivity
For independent runs Reset and set spatial mode n and diffusivity to and . Press Step twice for t=1 second. Predict amplitudes exp(-0.5), about 0.607, and exp(-2), about 0.135. This is one exact periodic heat-equation mode on a dimensionless spatial interval; it does not simulate calorimetry, thermometers, contacts, or a microscopic thermal-conductivity fit.
- Current amplitude
- 1
- Spatial L² norm
- 1.772
- Decay rate · κn²
- 3.6 /s
The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forA diffusion decay rate depends on spatial scale as well as the transport coefficient.
Pause & reason
Check your understanding
Your turn
Try it yourself
For an illustrative conventional elastic metal at 10 K, take and rounded . If total thermal conductivity is , find the benchmark electronic contribution and residual.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Use kappa_e=L0 T/rho, then subtract it only within the stipulated benchmark assumptions.
Compare with the solution
The benchmark gives , leaving . That residual is conditional on the electronic benchmark, matched geometry, and calibration. It could include phonons or other contributions and is not by itself a unique measurement of neutral quasiparticles.
Keep the idea
Take it with you
- A linear-plus-cubic heat-capacity fit assumes a specific low-energy spectrum and reliable addenda subtraction.
- Thermal conductivity depends on heat paths, equilibration, and geometry as well as excitations.
- A Lorenz comparison is a conditional electronic benchmark, not a unique mechanism detector.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.