Field

Thermodynamics

Heat, disorder, and the one law with a direction in it.
Molecular speeds at 4 temperatures. The distribution of molecular speeds in a gas, with each curve enclosing the same area. Raising the temperature moves the peak right and lowers it: the same molecules, spread over a wider range of speeds.

The speeds in a still room

The air in a quiet room is not still. Every molecule in it is moving at hundreds of metres per second, and temperature is a single number summarising an entire distribution.

Ways to arrange 10 coins. The number of distinct arrangements giving each number of heads, for 10 coins. Every individual arrangement is equally likely; the middle wins because there are more ways to reach it.

Entropy is a count, and the arrow of time is arithmetic

Nothing in mechanics prefers a direction. Entropy is not a force pushing things toward disorder — it is the observation that some outcomes have vastly more ways of happening than others.

A Carnot cycle on pressure–volume axes. Two isothermal steps joined by two adiabatic ones, forming a closed loop. The gas expands 9.6-fold in reaching the cold reservoir at 0.50 of the hot one, and the area enclosed is the net work done over one cycle.

The ceiling on every engine, set before it was designed

There is a maximum efficiency no heat engine can exceed, and it depends on nothing but two temperatures. Not the fuel, not the working substance, not the cleverness of the engineer.

Pressure, counted as momentum arriving at a wall. Molecules with speeds drawn from the Maxwell–Boltzmann distribution and directions drawn uniformly in the plane of the figure. Those moving toward the wall will bounce off it, reversing the component perpendicular to it and delivering twice that momentum each — and pressure is nothing but the rate at which that momentum arrives. Because the directions here lie in a plane rather than in space, the perpendicular component carries half the energy rather than the third it carries in a real gas.

Pressure is a rate of arrival, and the gas law falls out of counting

Nothing in a gas is pushing on the walls. Molecules arrive, bounce, and leave, and pressure is the momentum they deliver per second — from which the ideal gas law follows with no thermodynamics in it at all.

A spike, spreading. The solution of the diffusion equation at three times, with a seeded random walk histogrammed behind it. The area under every curve is the same because nothing is lost; only the width changes, and it grows as the square root of the time.

The equation that only runs forwards, and the walk underneath it

A drop of ink spreads and never gathers. The equation describing it is one of the few in physics that is not reversible — and underneath it is nothing but a coin being tossed.

Heating 1 kg of water from -20°C to 130°C. Temperature against heat added for 1 kilogram of water taken from -20 to 130 degrees Celsius. The two flat stretches are the melting and the boiling, where 334 and 2260 kilojoules go in and the temperature does not move. Melting costs as much as warming the water by 80 degrees; boiling costs as much as warming it by 541, which is 73 per cent of the whole journey.

The heat that changes no temperature, and where it actually goes

A kettle reaches a hundred degrees in a minute and takes five more to boil dry. The heat going in during those five minutes changes nothing a thermometer can see, and it is most of the energy in the whole process.

The isothermal atmosphere against the real one. Pressure as a fraction of its sea-level value, against altitude. The curves are the isothermal barometric formula at 220, 288, 400 kelvin, whose scale heights are 6.4, 8.4, 11.7 kilometres. The points are the measured standard atmosphere. At 20 kilometres the 288 kelvin model is 71 per cent out, because the air up there is not at 288 kelvin.

Why the air thins with height, and why that is the same law as the speeds

The pressure of the atmosphere falls exponentially with altitude, and the distribution of molecular speeds falls exponentially with energy. These are not two results that happen to look alike. They are one statement read on two axes.

A path through a crowd. A point crossing a field of 90 scatterers, rebounding off each. The mean length of 4000 such segments is 0.2983 box widths, against the textbook form 1/2nr = 0.3086 — a departure of -3.3 per cent, from a measurement that knows nothing of the formula. It does not agree exactly and should not: the closed form is derived for a vanishingly dilute field and these discs cover 9.2 per cent of the plane. Two finite-density effects pull opposite ways — crowding shortens the path, and discs shadowing one another lengthen it — so which side of the formula a given field lands on is not something the formula can tell.

How far a molecule gets

A molecule of air travels about sixty-eight nanometres between collisions — some two hundred times its own size, and a ten-millionth of the width of a room. That ratio is the reason a gas can be treated as a continuous substance at all, and the reason it sometimes cannot.

Seven walks from one point. 7 random walks of 400 steps each, all starting at the same place. None of them goes anywhere in particular and none of them stays put; the typical distance reached after n steps is the square root of n, so quadrupling the time doubles the spread. The tracks are seeded, so this is a property of the figure rather than of any one run.

The jiggle that proved atoms

A pollen grain in still water never stops moving. For eighty years that was a curiosity with no explanation; then it became the measurement that settled whether matter is made of particles, by turning a microscope and a stopwatch into a count of how many molecules are in a mole.

One mass, two entropies. The entropy of a solar mass as ordinary gas, generously counted at ten Boltzmann constants per proton, against the entropy of a solar-mass horizon. The first is 1.19·10⁵⁸ k and the second 1.05·10⁷⁷ k — a factor of 8.83·10¹⁸. This is why a horizon had to be given an entropy: without one, dropping anything at all through it destroys entropy and the second law fails.

The entropy that lives on a surface

Throw a cup of tea through a horizon and the entropy of the outside world falls. Either the second law is wrong or the horizon has an entropy of its own — and the only quantity available for it turns out to be its area, in units of a length made from gravity, quantum mechanics and the speed of light together.

The staircase equipartition cannot climb. The heat capacity of hydrogen at constant volume, in units of R, against temperature on a logarithmic axis. The three translational directions contribute 3/2 at every temperature. The two rotations switch on near 85.4 K — the temperature at which kT matches the first rotational step — taking the total to 5/2, computed here by summing the rigid rotor's partition function over two hundred levels rather than by assuming the plateau. The vibration switches on near 6332 K, taking it to 7/2. At 50 K the value is 2.469; at 300 K the value is 2.502; at 5000 K the value is 3.376. Classical equipartition predicts 7/2 at every temperature, including at four kelvin, and the size of the steps is set by ħ — which is how a heat capacity measures a quantum constant.

Half a kT for every way of moving

A heat capacity ought to be a count. Every quadratic term in a system's energy carries half a kT of it, so warming a gas is a matter of enumerating the ways its molecules can move — and the fact that the count comes out wrong for hydrogen is how a thermometer measured Planck's constant.

A ratio in the exponent. The Boltzmann factor against the energy of a state measured in units of kT, with the logarithm up the axis so that the straight line is the whole content. Every kT of energy costs a factor of e, so 10, 20, 30 times kT are factors of 10^-4.3, 10^-8.7, 10^-13.0. At 300 K, kT is 25.9 meV, so a barrier of 0.35 eV is 13.5 kT and a factor of 1.3e-6. That is the sense in which a third of an electronvolt is not a small energy: it is small compared with a chemical bond and enormous compared with kT, and it is the second comparison that decides whether anything happens.

The exponential that decides everything

Maximising the number of ways a reservoir can arrange what is left after taking E out of it gives one factor, e to the minus E over kT. Its exponent is a ratio, which is why a barrier of a third of an electronvolt — nothing at all by chemical standards — is the difference between instantly and never.

Ways to arrange 12 coins. The number of distinct arrangements giving each number of heads, for 12 coins. Every individual arrangement is equally likely; the middle wins because there are more ways to reach it.

What a system actually minimises

A ball falls to the bottom of a bowl and a gas fills a room, and neither of those is the rule. A system in contact with a large reservoir minimises U − TS, and the minus sign is the reservoir's own entropy written in the system's variables — which is why a rubber band pulls harder when it is heated.

The phase boundary of water, from one equation. Pressure against temperature for water on a logarithmic pressure axis spanning 6.6 decades. The vaporisation curve is integrated from Clausius–Clapeyron between the triple point at 273.16 kelvin and 611.7 Pa and the critical point at 647.096 kelvin, with a single latent heat of 43.32 kilojoules per mole — the value the two published points on the curve imply. The measured latent heats are 45.05 at the triple point and 40.65 at the reference point, and the fitted value sits between them, because a constant latent heat is an average over the interval. The sublimation curve below the triple point is not measured but predicted, from the two latent heats adding where all three boundaries meet: 51.1 kilojoules per mole, which reaches 103.2 Pa at 253.1 kelvin against a measured 253.15. The melting curve is drawn at the slope Clapeyron gives it, -13.5 megapascals per kelvin, which is a volume ratio and nothing else: water's solid is 917 against 1000 kilograms per cubic metre for its liquid, so melting shrinks it and the line leans backwards. Across the whole 6.6 decades of this axis that line moves 5.2 kelvin, and one atmosphere shifts the melting point by 0.0075 kelvin. At 1 atmosphere the boundary is crossed at 373.1 kelvin, where water boils. The one place the curve fails is its top end: a constant latent heat reaches 37.4 MPa at the critical temperature where the measured critical pressure is 22.1 MPa, 70 per cent high, because the latent heat falls to zero at the critical point and this curve does not know that.

A boiling point is a pressure, not a temperature

Nothing about water names one hundred degrees; the air does. Where a liquid turns to vapour in its bulk is fixed by what pushes on it, which makes the familiar figure a coordinate on a curve — 71 °C on Everest, 119.5 °C in a sealed pot — and one latent heat draws the whole curve.

The ceiling, inverted. How many joules of heat a perfect machine can move per joule of work, against the outside temperature, with the inside held at 21 °C. The upper curve is heating — T_h/(T_h − T_c), which is what the Carnot argument becomes when the cycle is run backwards — and the lower one is cooling the outside, T_c/(T_h − T_c). They differ by exactly one everywhere, to 1.8e-15 across the whole range as drawn, because the work put in is delivered as heat along with whatever was moved. The dashed line at one is a resistive heater, which is 100% efficient and is the worst option on the figure. At 7 °C and −7 °C the ideal coefficients are 21.0 and 10.5; a real machine reaching 25% of the ideal gets 5.3 and 2.6, which is still several times what burning the same energy would give.

The engine that pays back more than it takes

Carnot's argument puts a ceiling on how much work a flow of heat can be made to do. Run the same cycle backwards and the ceiling inverts into a floor that is greater than one — so a machine can deliver three or four joules of heat for every joule it consumes, and a perfectly efficient electric heater is the worst way to warm a room.

The entropy of mixing, and the entropy of not mixing. The entropy gained when two ideal gases at the same temperature and pressure are allowed to mix, per particle and in units of Boltzmann's constant, against the proportion of the mixture that is the first gas. The curve has no property of either gas in it — not their masses, not their sizes, not how strongly they interact, since ideal gases do not — and it is largest at 0.500, where it reaches ln 2 = 0.6931. Below it is the same quantity for two samples of the SAME gas, which is zero at every proportion: removing the partition between two halves of a box of nitrogen changes nothing that can be measured, and putting it back recovers the original state. The two results are correct and they do not join up. Make the two gases more and more alike — two isotopes, then two nuclear spin states, then nothing at all — and the upper curve does not descend to meet the lower one; it stays exactly where it is until the two species become identical, and then jumps. What the figure is really about is that the jump is in the counting and not in the gas.

Mixing what is already mixed

Let two different gases into each other's halves of a box and the entropy rises by a fixed amount that contains nothing about either gas. Do it with the same gas on both sides and it rises by nothing. Make the gases more and more alike and the answer does not converge — it jumps.

The part of the curve no fluid follows. Van der Waals' isotherms in reduced units, at 5 temperatures either side of the critical one, so that nothing about any particular substance appears. Above the critical temperature the pressure falls monotonically as the volume grows, which is what a fluid does. Below it the curve develops a loop with a rising section in the middle, and that section says the pressure increases as the substance expands — a material with negative compressibility, which cannot exist, because any fluctuation would run away. What happens instead is drawn as the horizontal line: the substance separates into two phases at one pressure, and the volume moves along the line as the proportions change. Its height, 0.6470 of the critical pressure, is fixed by requiring the two areas the line cuts off to be equal, which is the condition that the two phases have the same Gibbs energy. It meets the curve at volumes 0.603 and 2.349, a ratio of 3.9, and those are the densities of the liquid and its vapour. The two turning points of the loop, at 0.72 and 1.53, bound the part that is not merely unobserved but impossible; between them and the construction the substance can be made to sit, superheated or supercooled, until something nucleates.

The part of the curve no fluid follows

One equation for a real gas produces isotherms with a rising middle section, which says a substance would expand as the pressure on it grows. Nothing does that. What replaces it is a horizontal line whose height is fixed by making two areas equal, and the condition is not a convenience.

The viscosity of a gas, over 8 decades of pressure. The viscosity of 3 gases at 300 K against pressure, on a logarithmic pressure axis and a linear viscosity one. The lines are flat, and that is the whole figure. Viscosity is the rate at which momentum is carried across a shear, which is the density of carriers times the distance each one carries it: ⅓ρv̄λ. Doubling the pressure doubles the density and halves the free path, and the two cancel exactly, so the same gas at a hundredth of an atmosphere is exactly as viscous as at one — which is not what anybody expects of a thinner gas and is what is measured. nitrogen comes out at 17.9 μPa·s against a measured 17.9, helium comes out at 19.3 μPa·s against a measured 19.9, argon comes out at 21.7 μPa·s against a measured 22.7. The flatness ends when the free path reaches the apparatus rather than the next layer of gas: at a vessel 10 mm across that is around 0.68 Pa for nitrogen, 1.94 Pa for helium, 0.69 Pa for argon, below which there is no gas-to-gas hand-off left to make.

The viscosity that does not care how much gas there is

Pump most of the air out of a vessel and the air that is left is exactly as viscous as it was. Maxwell derived that in 1860, did not believe it, and spent six years building an apparatus to measure it — which is a better description of how a prediction becomes knowledge than any amount of agreement would have been.

The average position of something that is only shaking. The mean displacement of an oscillator against temperature, taken as a Boltzmann average over each well rather than from any expansion of it, with the temperature measured against each well's own depth so that unlike bonds can share an axis. A symmetric well gives exactly zero at every temperature — heating a harmonic solid makes it vibrate harder and does not make it larger. The others drift outward, because the outward side is the shallower one, and the measured slopes are a pendulum 0.000, a chemical bond 0.766, a pair of atoms 0.150. Thermal expansion is not a property a spring has; it is one a spring lacks.

Why heating a perfect spring changes nothing

A harmonic solid vibrates harder when heated and does not get any longer. Thermal expansion lives entirely in the term that the harmonic approximation throws away — and so does the fact that a solid conducts heat at a finite rate, which is the same discarded term doing a second job nobody would have connected to the first.

Two straight lines that are not the same line. Half the difference between the liquid and vapour densities, in units of the critical density, against the distance from the critical temperature — both logarithmic, over 5 decades. The van der Waals curve is solved for the coexisting pair at each temperature and its slope in the last decade is 0.500, which is the mean-field ½. Real fluids give 0.326. The two differ by 11 per cent at t = 0.1 and by a factor of 7.9 at the bottom of the axis, which is why an equation with the wrong exponent in it looked right for eighty years.

The point at which the two become one

Heat a sealed tube of carbon dioxide and the meniscus inside it does not boil away — it fades, the two densities converging until there is nothing to separate. Twenty millikelvin before that happens the fluid turns milky, and the exponent describing the last approach is a number van der Waals got wrong and could not have got right.

The work one molecule and one bit are worth. The pressure of a gas of one molecule against its volume, at 300 K, in units of the volume it starts in. The shaded area is the work the molecule does pushing a partition out isothermally, and it is measured here by integrating the drawn curve rather than written down: expanding by 1.5× yields 0.4055 kT against ln 1.5 = 0.4055; expanding by 2× yields 0.6931 kT against ln 2 = 0.6931; expanding by 4× yields 1.3863 kT against ln 4 = 1.3863; expanding by 8× yields 2.0794 kT against ln 8 = 2.0794, agreeing to 2.0e-10. The doubling is the one that matters, because a partition inserted in the middle leaves the molecule on one side or the other, and knowing which is what lets the load be attached to the right face. That single expansion delivers kT·ln2 = 2.87 zeptojoules at 300 K. It looks like work extracted from one temperature, and it is — until the engine is asked to run again, which requires forgetting which side the molecule was on.

The bit that has to be paid for

One molecule in a box, a partition, and the knowledge of which side it went — enough, between them, to extract work from a single reservoir, which the second law forbids. The engine is real and the arithmetic is right. What closes the loophole is that the cycle does not finish until the knowledge has been thrown away, and throwing away one bit costs exactly what the expansion delivered.

The ceiling, the estimate, and three power stations. Two efficiencies against the ratio of the cold reservoir's temperature to the hot one. The upper curve is Carnot's 1 − Tc/Th, which is a ceiling on the work per unit of heat and is reached only by an engine that runs infinitely slowly, because a reversible heat flow needs a vanishing temperature difference to drive it and therefore infinite time. The lower curve is 1 − √(Tc/Th), the efficiency of an engine with finite thermal contact run for the most power rather than the most work. Three measured plants are marked: West Thurrock, coal runs at 36 per cent against a ceiling of 64 and a finite-time estimate of 40; CANDU, nuclear runs at 30 per cent against a ceiling of 48 and a finite-time estimate of 28; Larderello, geothermal runs at 16 per cent against a ceiling of 33 and a finite-time estimate of 18. Every one of them is closer to the lower curve — within 4.4 points at worst, against 16.5 at best from the ceiling. The second law is not what limits a working power station. What limits it is that somebody wants the electricity this year.

The engine that has to finish

Carnot's ceiling is exact and it is reached only by an engine that takes for ever, because a reversible heat flow needs a vanishing temperature difference to drive it. Ask instead for the most power rather than the most work per joule of heat, and the answer is a different function of the same two temperatures — and three measured power stations sit on it rather than on the ceiling.

What a coordinate is worth, by the shape of its energy. The mean energy stored in one coordinate, in units of kT, against the power with which that coordinate enters the energy. The curve is 1/n and the points are the same average obtained by integrating the Boltzmann weight numerically, agreeing to 3.3e-7 per cent at worst. an energy going as |x|^1 holds 1.0000 kT; an energy going as |x|^1.5 holds 0.6667 kT; an energy going as |x|^2 holds 0.5000 kT; an energy going as |x|^3 holds 0.3333 kT; an energy going as |x|^4 holds 0.2500 kT; an energy going as |x|^6 holds 0.1667 kT. The familiar half a kT is the case n = 2 and nothing more general than that: a coordinate whose energy is linear in it — the momentum of an ultrarelativistic particle, or a field with no restoring force but a constant tension — carries a whole kT, and one confined by a very steep wall carries almost nothing. Equipartition is a theorem about quadratic terms, and calling it a theorem about degrees of freedom is the substitution that makes it fail.

The share that is not half a kT

Equipartition is quoted as half a kT for every degree of freedom, and it is nothing of the kind. It is half a kT for every *quadratic* term. A coordinate whose energy is linear in it carries a whole kT, and a gas hot enough that its particles' energy is pc rather than p²/2m therefore holds twice what the counting says — which drops its ratio of specific heats to four thirds and puts a star on the edge of being able to hold itself up.

What a thermal camera reads off a shiny surface. The temperature a camera calibrated for a black body reports, against the emissivity of the surface it is pointed at, for surfaces truly at 323 K, 373 K, 473 K in a room at 293 K. Every curve begins at the true temperature when the emissivity is one and ends at the room's temperature when it is zero, because a surface that emits nothing reflects everything and the camera is then looking at the room. polished aluminium at 0.05 reads 312 K; stainless steel, oxidised at 0.8 reads 451 K; matt black paint at 0.95 reads 468 K; human skin at 0.98 reads 471 K for a surface truly at 473 K. That is not a fault in the instrument. It follows from Kirchhoff's law — a poor emitter is a poor absorber and therefore a good reflector — so the shortfall in a shiny surface's own glow is made up almost exactly by whatever is reflected in it, and no measurement of the light leaving a surface can separate the two without knowing the emissivity in advance.

The glow that says nothing about the surface

A thermal camera pointed at a saucepan of boiling water reads a hundred degrees. Pointed at a polished aluminium block at the same temperature it reads about twenty-five, and the instrument is working perfectly. What it is measuring is emissivity as much as temperature — and a surface's emissivity is forced to equal its absorptivity, at every wavelength and every angle, by an argument with no physics of matter in it at all.

Each stage takes 25.0 per cent of what is left. Entropy against temperature for a spin-½ paramagnet at 0.25 T and 1 T, with the cooling cycle drawn between them: a vertical drop is isothermal magnetisation, a horizontal move is adiabatic demagnetisation. Starting from 1 K the treads are at 1.000 K, 0.250 K, 0.062 K, 0.016 K, 3.91 mK. Each is 0.2500 of the one before — a ratio read back off the drawn treads rather than written into them, and equal to the field ratio 0.25/1 because this refrigerant's entropy depends on the field and the temperature only through their quotient. The steps therefore shrink in proportion to what is left, and no finite number of them arrives.

The staircase that never reaches the floor

Absolute zero is unreachable, and the reason is not that the apparatus is not good enough. Every stage of cooling removes a fixed fraction of what is left rather than a fixed amount, so the steps shrink in proportion to the distance remaining — and the fixed fraction cannot be made one, because the entropy curves at two field strengths are required to meet where the axis is.

The fraction of walks that have come home, against how long they have walked. The proportion of 1,600 lattice random walks that have returned to their starting point at least once, against the number of steps taken, on a logarithmic horizontal axis, in one, two and three dimensions. In one dimension almost every walk is home almost at once and the fraction climbs toward one. In two it climbs more slowly — the return is still certain, but only logarithmically, so a two-dimensional walk that has not come back after a thousand steps is unremarkable. In three the curve flattens: it reaches 0.3481 and stops, against the exact value 0.3405, drawn as a line. That number is Watson's integral, and it is the probability that a three-dimensional walk ever comes home at all. The difference between the cases is not one of degree. In one and two dimensions the expected number of returns is infinite and a diffusing particle visits every site eventually; in three it is finite, and a molecule released in a room will, with probability two-thirds, never pass through its starting point again. The same statement runs the other way round: a reaction that needs two diffusing partners to meet is a very different problem on a membrane from what it is in a cell.

The walk that comes home

A particle wandering at random on a line returns to where it started, with certainty. On a plane it returns, with certainty. In three dimensions the probability is 0.3405 — so two out of three molecules released in a room never pass through their starting point again, and the difference between the cases is not a matter of degree.

Entropy that is still there at absolute zero, counted and measured. Four substances whose entropy does not go to zero when they are cooled as far as anybody can cool them, with the entropy counted from the arrangements they froze into beside the entropy measured by integrating their heat capacities. Ice is the famous one: every oxygen has four hydrogen bonds and the rule is that two hydrogens sit near it and two far, which leaves six of the sixteen placements legal, and Pauling's count of the whole crystal collapses to R ln(3/2) = 3.371 J per mole per kelvin against a measured 3.41 — an agreement to one per cent from an argument on one line. Carbon monoxide and nitrous oxide are the easy cases, molecules that can lie either way round in the lattice and have too little to gain by choosing. The worst of the four is off by 25 per cent, which is the honest state of this subject: the counts are crude, they ignore the correlations between neighbouring choices, and they still land within sight of a calorimeter. What the figure is really about is that the third law has an escape clause and the escape clause is measurable. A perfect crystal has one arrangement and zero entropy; a crystal that ran out of time while it had many has the logarithm of however many it stopped at, permanently.

The entropy that is still there at zero

The third law says a perfect crystal has no entropy at absolute zero. Ice has 3.41 joules per kelvin per mole left over, and the number can be recovered from one line of counting — two hydrogens near each oxygen and two far, six legal arrangements out of sixteen, R ln(3/2). The law has an escape clause and the escape clause is measurable.

The gas that leaves is not the gas inside. The distribution of molecular speeds inside a container at 300 K and in the beam that escapes through a small hole in it, each normalised to its own peak. They are not the same distribution. A molecule's chance of reaching the hole in a given time is proportional to how fast it is going, so the beam carries one more factor of speed than the gas does — v³ rather than v² times the Boltzmann factor — and the beam is therefore faster and hotter than what it came from. The mean speed inside is 476 m/s and in the beam 561 m/s, a ratio of 1.1781 against the exact 3π/8; and the mean kinetic energy is 1.500 kT inside against 2.000 kT in the beam, which are exactly 3/2 and 2. That difference is not a subtlety. A molecular beam made by effusion has a temperature, in the sense of a mean energy, a third higher than its source; a gas slowly leaking from a container leaves the remainder cooler than it would be if a fair sample had gone; and every calculation of a rate through an aperture that uses the bulk distribution is wrong by this factor.

The gas that leaves is not the gas inside

Put a small hole in a container of gas and what comes out is faster and hotter than what stays behind — its mean kinetic energy is 2kT against the 3/2 kT of the gas it came from. Nothing has heated it. A fast molecule simply reaches the hole more often than a slow one, so the sample that escapes is biased by exactly one factor of speed, and every consequence of effusion is that factor.

The barrier a new phase has to climb. The free energy of a droplet against its radius, at four supersaturations. Two terms compete: the volume term is a gain and goes as r³, the surface term is a cost and goes as r². At small radius the surface wins, so a droplet that forms by chance is more expensive than the vapour it came from and evaporates again; past a critical radius the volume wins and the droplet grows without limit. The maximum between them is the barrier. At S = 1.5 the critical radius is 2.66 nm and the barrier 533.8 kT, S = 2 the critical radius is 1.56 nm and the barrier 182.6 kT, S = 3 the critical radius is 0.98 nm and the barrier 72.7 kT, S = 5 the critical radius is 0.67 nm and the barrier 33.9 kT. The critical radius contains a few hundred molecules at low supersaturation and a handful at high, which is the first sign that a theory built on a surface tension and a bulk free energy is being applied outside its comfort. Both the critical radius and the barrier are located here by searching the drawn curve and checked against the closed forms 2γ/|Δg| and 16πγ³/3Δg², which agree to a part in a thousand.

The barrier a new phase has to climb

Water vapour three times supersaturated is thermodynamically desperate to condense and will sit there indefinitely if it is clean enough. The obstacle is that a droplet has to start small, and a small droplet is nearly all surface — so the first nanometre of every phase transition costs energy rather than releasing it, and what decides whether anything happens is the height of that cost divided by kT.

The temperature that runs off the top of the scale. On the left, the entropy of a collection of two-level systems against how many of them are in the upper level. It rises, reaches 0.69314 per spin — which is ln 2, and is where half of them are up — and then falls, because a system with every spin up is as orderly as one with every spin down. On the right, the slope of that curve, which is one over the temperature. Below the maximum it is positive and ordinary: adding energy adds entropy, and the temperature is what everybody expects. At the maximum it is zero, which means the temperature is infinite. Past it the slope is negative — adding energy now removes entropy — and the temperature is negative. Such a system is not cold. It is hotter than any positive temperature whatever: put it in contact with anything at all and energy flows out of it, because that raises the total entropy. The quantity that orders systems by which way heat flows is not the temperature but its reciprocal, which runs smoothly from large and positive through zero to negative, and it is the temperature that has the discontinuity. None of this is possible unless the energy has a ceiling, which is why it happens in a spin system and not in anything that can move: a gas has no upper bound on its kinetic energy, so its entropy never turns over and its temperature is never negative.

Hotter than any temperature there is

A system whose energy has a ceiling can be pushed past the point where adding energy adds entropy. Its temperature is then negative — and negative temperatures are not cold. They sit above every positive temperature on the only scale that decides which way heat flows, and a working laser is at one.

Throttling a gas, and the curve that says which way it goes. Curves of constant enthalpy for a van der Waals gas, in temperature and pressure both measured against the critical values. A gas pushed slowly through a plug or a valve keeps its enthalpy, so it moves along one of these curves — from right to left, since the pressure falls. Where a curve slopes upward to the right the gas cools as it expands; where it slopes downward it warms. An ideal gas would give horizontal lines and no change at all, because its enthalpy depends on the temperature alone; every curve here is bent, and the bending is the attraction between molecules and the room they take up, fighting. The dashed line through the tops of the curves is the inversion curve, and the maxima were found on the drawn points rather than put there — they lie on the closed form to 0.65 per cent. Which side of it a gas starts on decides the sign of the effect, and that is the whole of why air can be liquefied by throttling at room temperature and hydrogen cannot: hydrogen has to be pre-cooled below its own inversion temperature first, which is why Dewar needed liquid air before he could get liquid hydrogen, and why Onnes needed liquid hydrogen before he could get helium.

The gas that cools by being let go

Push a gas through a plug and its temperature changes, although no work is done on anything and no heat goes anywhere. Which way it changes depends on where it starts: inside a dome in the pressure–temperature plane it cools, outside it warms, and hydrogen at room temperature is outside — which is why liquid hydrogen needed liquid air first.

Energy, pressure and entropy of a gas nobody counted. The energy density, pressure and entropy density of blackbody radiation against temperature, on logarithmic axes, together with the pressure a monatomic gas of the same energy density would have. Every curve is a power of the temperature — the fourth for energy and pressure, the third for entropy — because the only length in the problem is the thermal wavelength and the only energy is kT. The pressure is exactly a third of the energy density, where an ordinary gas's is two thirds, a factor of 2: a photon carries momentum E/c and a slow molecule carries √(2mE), and that difference is the whole of it. Some values: at room temperature the radiation pressure is 1.86e-6 pascals, which is a ten thousand millionth of an atmosphere; at 1e+7 kelvin it is 2.52e+12, which is where radiation rather than matter holds a star up.

The gas that nobody counted

A box of gas holds however many molecules were put in it. A box of radiation holds however many photons the temperature says, because the walls make and destroy them until the free energy is least — and one dropped assumption changes every result. The pressure becomes a third of the energy density instead of two thirds, the entropy goes as the cube of the temperature, and the adiabatic index comes out at exactly four thirds.

Two experiments, two computations, one coefficient. The Seebeck coefficient of a resonant conductor against where its resonance sits relative to the chemical potential, together with the Peltier coefficient divided by the temperature. The first is obtained by applying a temperature difference and finding the voltage that stops the current; the second by applying a voltage at uniform temperature and taking the ratio of the heat flow to the current. Different driving, different measurement, different integral — and the two curves agree to 3.8e-7 of the sweep's own scale across the whole of where they pass through zero and change sign. That equality is Kelvin's relation Π = ST, guessed in 1854 from an argument its author knew was not sound and proved by Onsager in 1931 from microscopic reversibility. It is not a property of this conductor; it holds for every one.

The second experiment that cannot disagree

Heat one end of a wire and a voltage appears across it. Pass a current through the same wire at uniform temperature and it carries heat. Those are two different experiments with two different apparatus, and the coefficient in front is the same number in both — not approximately, and not for some materials. The reason is that the equations of motion underneath look the same run backwards.

A stiffness that does not go to zero — it falls off a cliff. The renormalised stiffness of a two-dimensional superfluid against temperature, obtained by integrating the flow of the bare stiffness and the vortex fugacity, with the bare stiffness drawn for comparison. Below the transition the vortices are bound in pairs, screen one another and leave a finite stiffness; above it they unbind and there is none. What is remarkable is the value at which it disappears: 0.6430, against 2/π = 0.6366. That number does not depend on the material, on the bare fugacity, on the core energy, or on anything else — every two-dimensional superfluid loses its stiffness at exactly two over π times its own transition temperature, and measurements on helium films of widely different thicknesses fall on that line. A transition at which a quantity jumps to zero from a value nobody can adjust is unlike anything the usual classification of transitions describes.

The transition with nothing to order

Every transition in this collection so far has an order parameter — a quantity that is zero on one side and not on the other. In two dimensions a continuous symmetry cannot break at any temperature above zero, so there is nothing for such a quantity to be, and by the usual reckoning there can be no transition. There is one anyway, and what changes at it is whether vortices are bound in pairs.

The year, arriving underground at four different times. Temperature against depth in soil of diffusivity 0.5 mm²/s, at four times of year, measured from the annual mean. The surface swings by ±12 K; the swing underground is smaller and later, and both are governed by one length, δ = √(2D/ω) = 2.24 m. The amplitude falls as e^(−z/δ) — the dashed envelope — and the phase lags by z/δ radians, so the curves lean over as they go down and cross the axis at different depths. At 7.0 m the lag is half a year: the ground there is at its coldest in August and its warmest in February, in antiphase with the sky, with a swing of 0.52 K left. That is a cellar, and it is also why a water main below about a metre and a half does not know that it froze last week.

The summer that reaches the cellar in December

Drive the diffusion equation at its boundary instead of releasing something into it and the solution is a decaying, lagging oscillation with a single length in it. That length governs both the shrinking and the delay, which is why the depth at which the ground is coldest in August is fixed by the same number as the depth at which the seasons stop being felt at all.

Three curves that can only meet at a point. Water's phase diagram within a kelvin of its triple point, with each boundary drawn at the slope Clausius and Clapeyron give it from measured latent heats and densities: 44.4 Pa/K for boiling, 50.3 for sublimation, and -135 bar per kelvin for melting — negative, and so steep that it is drawn vertical here, because ice is the less dense phase and pressure therefore melts it. Two things follow that no amount of measurement could adjust. The three slopes are not independent: the sublimation and vaporisation curves differ in slope by 5.92 Pa/K, which is exactly what the fusion latent heat predicts, because going solid → gas directly and solid → liquid → gas must cost the same enthalpy. And the sublimation curve is the steeper of the two, so the two cross rather than touch — which is why the triple point is a crossing and not a tangency, and why ice sublimes below it and melts above it.

Why the triple point is a point

Three phases of one substance coexist at one temperature and one pressure and nowhere else, and the reason is arithmetic rather than chemistry: count the numbers that describe the state, count the conditions equilibrium imposes, and subtract. The same subtraction says four phases of one substance are impossible, and it says so without knowing what the substance is.

A ratio in the exponent. The Boltzmann factor against the energy of a state measured in units of kT, with the logarithm up the axis so that the straight line is the whole content. Every kT of energy costs a factor of e, so 5, 12, 20 times kT are factors of 10^-2.2, 10^-5.2, 10^-8.7. At 295 K, kT is 25.4 meV, so a barrier of 0.4 eV is 15.7 kT and a factor of 1.5e-7. That is the sense in which a third of an electronvolt is not a small energy: it is small compared with a chemical bond and enormous compared with kT, and it is the second comparison that decides whether anything happens.

The temperature a molecule does not have

Temperature fixes a system's average energy and nothing more. The actual energy wanders, by an amount tied to the heat capacity, and the relative size of the wandering falls as one over the square root of the number of degrees of freedom — so a mole has a temperature and a molecule does not.

Melting curves, and the one that leans the wrong way. Melting temperature against pressure for water, benzene, naphthalene, each measured from its own melting point at one atmosphere, with pressure in bars. The slope of every coexistence line is the latent heat divided by the temperature and the change in volume, and the latent heat of melting is positive for everything — so the sign of the slope is the sign of the volume change, and nothing else. Almost everything expands on melting and its line leans forwards. Water's solid is less dense than its liquid, so its line leans backwards at 135 bars a kelvin: pressing on ice at just below zero melts it, and it takes 135 atmospheres to gain a single degree. The anomaly is not in the thermodynamics; it is in the fact that ice floats.

The melting curve that leans the wrong way

The slope of any coexistence line is the latent heat divided by the temperature and the change in volume. Latent heat is always positive, so the sign of the slope is the sign of the volume change — and for water the volume change is negative, which is the whole of why ice floats and why the melting curve leans backwards.

Effusion rate against molecular mass. The rate at which a gas escapes through a small hole, against its molar mass, normalised to hydrogen. The rate is a quarter of the number density times the mean speed times the area, and the mean speed goes as the inverse square root of the mass — so the rate does too, which is Graham's law of 1848. Hydrogen escapes four times faster than oxygen and 13.3 times faster than uranium hexafluoride. The practical consequence is isotope separation, and its difficulty is on this chart. The two uranium hexafluorides differ by 3 out of 352 in mass, so a single stage enriches by a factor of only 1.00429 — four parts in a thousand. Reaching 90 per cent from natural uranium's 0.72 per cent therefore takes about 1665 ideal stages, and a real cascade needs more because each stage is imperfect. That number is why gaseous-diffusion plants were among the largest industrial structures ever built, and why centrifuges — which separate by mass directly rather than by the square root of it — replaced them.

The gradient that drives the other thing

A concentration gradient drives a flow of matter and a temperature gradient drives a flow of heat. Each also drives the other, by coefficients that are equal — a relation nobody could have guessed and which follows from the fact that the underlying motion runs the same forwards and backwards in time.

Runs that break the second law, and how often. The work done in a process repeated many times, and the same for the process run in reverse with its work reflected, for a free-energy change of 4 kT and a dissipation of 3 kT. The average work exceeds the free-energy change, which is the second law, and individual runs do not have to: the shaded tail is the fraction of runs that do less work than the free energy — trajectories in which the entropy of the universe went down — and it is 11.03% here. The two curves cross exactly at the free-energy change, whatever the dissipation, which is what makes an irreversible measurement able to report an equilibrium quantity.

The second law, with a probability attached

Entropy increases, on average. For a small system pulled quickly, individual runs go the other way — and how often is not a matter of taste but an exact number, fixed by a relation with no adjustable constant in it and no requirement that anything be near equilibrium.

The pair potential, and the two things it does to a gas. The Lennard-Jones potential between two molecules, in units of its own depth and range, with the Mayer function it produces at 1, 3, 8 times the well depth in temperature. The virial coefficient is minus the integral of that function over volume, so the two parts of the potential contribute with opposite signs: the steep repulsive core makes the function minus one there, giving a positive contribution — molecules take up room — and the attractive well makes it positive, giving a negative one. At low temperature the attraction dominates and a gas is easier to compress than an ideal one; at high temperature the core dominates and it is harder. Between them is one temperature at which they cancel.

The first correction to the gas law

An ideal gas has no forces between its molecules. The first correction to what it does is computable from those forces alone — one integral over the pair potential — and its sign flips at a temperature where a real gas obeys the ideal law without being ideal at all.

One temperature, four pawls, and nothing gained. The net rate of a ratchet whose gas and whose pawl are at the same temperature, against the load, for notches 2, 5, 10, 20 times the thermal energy deep. Every curve passes through zero at zero load and is negative everywhere else. The device is not merely unable to lift a weight; under any load at all it turns the wrong way and lets the weight down, converting its potential energy into heat in the gas. Making the notch deeper slows everything down — an exponential in the depth — and does not change the sign anywhere. That is the second law arriving as a mechanism rather than as a prohibition. Nothing was assumed about entropy; the pawl was simply allowed to be as warm as everything else, and its own fluctuations undo exactly the rectification it was there to provide. Any rectifier small enough for thermal noise to matter has this problem, and the rectifier being clever does not help, because the same noise reaches the cleverness.

The engine a fluctuation cannot run

A ratchet lets a shaft turn one way and not the other. Put a paddle in a gas on the same shaft and molecular collisions appear to become a lifted weight — an engine running on one reservoir. It does not work, and following exactly why turns the second law from a prohibition into a mechanism: the pawl is as warm as the gas, and it lifts whenever it is asked to.

The energy that goes and comes back. A chain of 32 masses with springs a few per cent nonlinear, started with all its energy in its longest mode, with the energy of the first five modes followed against time in units of that mode's own period. The first mode gives up most of what it has — down to 9 per cent by 104 periods — and the energy appears in the second, third and fourth. Then it comes back: at 154 periods the first mode holds 98 per cent of the total again. Equipartition would put an equal share in every one of the thirty-two modes and leave it there. What happens instead is that a handful of modes trade with each other and return almost exactly to where they began, and go on doing so. The total energy is checked against its starting value throughout and holds to 9.2e-5, so nothing here is the integrator losing track of what it was given.

The energy that refuses to be shared

Put all the energy of a chain of masses into its longest mode and add a few per cent of nonlinearity, and equipartition says it should spread out among all thirty-two modes and stay there. It does not. It leaks into three or four neighbours and then comes back — almost exactly — and goes on doing so, and the calculation that found this was expected to be a demonstration that it would not happen.

The temperature a planet ought to be. Each body's measured surface temperature against the temperature at which it would radiate away exactly the sunlight it absorbs — computed from two numbers, the sunlight reaching it and the fraction it reflects, with no other property of the body used. Points on the diagonal are bodies the one-line argument gets right. the Moon: balance 270 K, surface 250 K, −20 K with no atmosphere; Mercury: balance 433 K, surface 340 K, −93 K with no atmosphere; Mars: balance 210 K, surface 210 K, +0 K; Earth: balance 255 K, surface 288 K, +33 K; Venus: balance 227 K, surface 737 K, +510 K. The airless bodies fall below the line rather than on it, and the reason is the fourth power: a surface running from noon heat to night cold radiates like its hottest parts and averages like its coldest, so a mean thermometer reading is lower than the temperature that matches the emitted flux. Mars, whose atmosphere is thin and whose surface is nearly isothermal by comparison, sits on the line. Everything with a substantial atmosphere sits above it, by tens of kelvin on Earth and hundreds on Venus, and always in the same direction. Nothing here explains why. What the figure fixes is the size and the sign of what has to be explained.

The height a planet is seen from

A body in sunlight settles where it radiates away what it absorbs, and that takes two numbers and one line of arithmetic. It gets the Moon right and Earth wrong by thirty-three kelvin. The correction is not that the atmosphere traps heat but that it moves the level space sees the planet from, and the rest is done by a lapse rate that is not a radiative quantity at all.

The entropy a substance keeps depends on how fast it was cooled. The entropy of a supercooled liquid in excess of its crystal's, against temperature, for a substance melting at 305 K with an entropy of fusion of 43 joules per kelvin per mole and a liquid heat capacity exceeding the crystal's by 62. The equilibrium curve — the one the liquid follows while it can still relax — falls steadily and would reach the crystal's entropy at 152.4 kelvin. It never gets there, because the liquid falls out of equilibrium first, at a temperature that depends on how long it is given: 0.01 s/K freezes at 196.1 K with 15.61 left, 1 s/K freezes at 188.1 K with 13.03 left, 100 s/K freezes at 182.1 K with 11.03 left. Ice's residual entropy is a count and does not move; a glass's is whatever it happened to have when it stopped being able to change, and that is a property of the experiment.

The entropy that depends on how fast it was cooled

Ice's residual entropy is a count, and it comes out the same whoever measures it. A glass's does not. A glass keeps whatever entropy it happened to have when its own relaxation time crossed the experiment's, so cooling ten times more slowly leaves less behind — and extrapolating the equilibrium liquid below that point takes its entropy under the crystal's at a finite temperature, which cannot happen and does not, for a reason that is still argued about.

The entropy each degree of freedom has not yet given up. The entropy carried by each kind of degree of freedom, drawn against the temperature at which it orders and hands that entropy over. Lattice vibrations freeze out around room temperature; electron spins in a paramagnetic salt order in the millikelvin range, which is what makes adiabatic demagnetisation work; and nuclear spins hold R ln(2I+1) — 11.5 joules per kelvin per mole for copper — down to some tens of nanokelvin, where their own dipolar interactions finally sort them out. A copper sample at a microkelvin therefore has a large entropy and violates nothing: its nuclear spin system has not reached its ground state, and the third law is a statement about ground states rather than about thermometers.

A law about spectra, not about heat

The third law is usually met as a statement about cooling. Its statistical form is a statement about a spectrum: the entropy of a system in its ground state is k ln g, and it vanishes only when the ground state is unique. Every apparent exception is a degeneracy or a system that never reached its ground state — and copper nuclei carry eleven joules per kelvin per mole down to a hundred nanokelvin without violating anything.

The rectangle every cycle is equal to. An ideal Otto cycle for air on a temperature–entropy diagram: compression ratio 9, intake at 300 K, 1400 kJ/kg added at constant volume. Compression takes the charge to 722 K, combustion to 2672 K, expansion back to 1110 K, and the exhaust cools at constant volume. The shaded loop is the work; the region under the lower curve is the heat rejected. Heat enters over a range of temperatures, and its entropy-weighted mean — heat divided by the entropy it brings — is 1491 K; the heat leaves at a mean of 619 K. The dashed rectangle between those two temperatures has exactly the loop's area, and 1 − 619/1491 = 58.5%, which is the Otto efficiency, checked to rounding. A Carnot engine between the coldest and hottest points of the same cycle would reach 88.8%.

The temperature an engine really takes its heat at

Carnot's ceiling is set by two temperatures, and no engine that burns fuel takes its heat in at one temperature or gives it out at another. It takes heat over a range, from the moment combustion starts to the moment it ends. For any reversible cycle there is an exact replacement for Carnot's two numbers: the average temperature at which heat arrives and the average at which it leaves, each weighted by the entropy the heat carries. The gap between a real cycle and Carnot is a gap between those averages and the extremes.

Two bodies an engine draws together. Two equal bodies of 4.186 kJ/K — a kilogram of water each — one at 90.0 °C and one at 10.0 °C, against the heat drawn from the hot one. The solid curves are the best possible engine running between them, a reversible one, which leaves the product of the two temperatures unchanged and brings both to the geometric mean, 320.7 K (47.5 °C). It draws 177.8 kJ from the hot body and delivers 20.8 kJ of work, C(√T₁ − √T₂)², checked against the heat balance. The dashed lines are the same bodies simply touching: they meet at the arithmetic mean, 323.1 K, having exchanged 167.4 kJ and delivered nothing. The 2.5 K between the two endpoints is the work, left behind as heat.

The work left in two buckets of water

Carnot's ceiling assumes reservoirs so large that taking heat from one and giving it to the other changes neither temperature. Two buckets of water are not reservoirs. Run the best possible engine between a hot one and a cold one and both temperatures move, the efficiency available shrinks as they do, and the engine stops when they meet — at the geometric mean of the starting temperatures, not the ordinary one. The work it delivered is exactly the difference between those two meeting points, and it is far less than the starting temperatures promise.

Two entropies that agree until half filling. The entropy per unit of 100 two-level units against the fraction excited, by Boltzmann's definition, the logarithm of the number of arrangements at that energy, and by Gibbs's, the logarithm of the number at or below it. Below half filling they nearly coincide: at a quarter excited they are 0.538 and 0.542 per unit, and both approach the dashed large-system curve. Boltzmann's entropy then turns over and falls back to zero when every unit is excited; at three-quarters it is 0.538. Gibbs's cannot fall, because a running total cannot, and it levels off at ln 2 = 0.693, reaching 0.693 at three-quarters. The slope of each is one over its temperature.

The count that decides which entropy is right

There are two ways to count the states of an isolated system: the states at its energy, which is Boltzmann's entropy, and the states at or below it, which is Gibbs's. For large systems in ordinary conditions they agree to the last measurable digit. For a system whose energy has a ceiling, past the halfway point, one gives negative temperatures and the other forbids them. Definitions cannot settle which is right, but a temperature is for something — saying which way heat will flow — and putting two such systems in contact lets the count of states answer.

Nitrogen that runs the wrong way. The nitrogen mole fraction in each of two bulbs joined by a capillary, as in Duncan and Toor's experiment: one bulb starts with 0.50086 nitrogen and the rest carbon dioxide, the other with 0.49879 nitrogen and the rest hydrogen, at 35 °C. Solid curves: the capillary solved at each instant from the Maxwell–Stefan equations with the three pairs' diffusivities, 83.8, 68.0 and 16.8 mm²/s. Dashed: Fick's law for nitrogen alone, which can only let the two start values relax together. At the start the nitrogen gradient is 0.00207, yet nitrogen flows at 309 times the rate that gradient would drive. From 0.1 to 6.5 hours it flows from the bulb with less nitrogen into the bulb with more, opening a difference of 0.1468 at 6.5 hours; at 6.6 hours its flux passes through zero with a difference of 0.1468 still in place. Each gas is conserved to a part in 10⁹. The carbon dioxide moving out of the first bulb drags nitrogen with it, because the nitrogen–carbon dioxide pair has by far the smallest diffusivity and so the strongest friction.

The gas that flows towards more of itself

Fick's law says a substance diffuses from where there is more of it to where there is less. In a mixture of three gases, nitrogen can do the opposite for hours — flowing into the bulb that already holds more nitrogen, and then stopping while a difference remains — and in a welded bar of steel, carbon crosses into the side that is already richer. Nothing is wrong with the second law. Diffusion flattens chemical potential, and with more than two components, or a second element changing it, that is not the same as flattening concentration.

The condition three modes never quite satisfy. By how much three modes of a thirty-two mass chain fail to be in resonance — the sum of two mode frequencies minus the frequency of their sum, on a logarithmic scale, against the second of the two modes for four choices of the first. The leading nonlinear term couples modes in threes and the exchange accumulates only where this quantity is zero. It never is: the chain's dispersion is a sine, a sine is concave, and the sum of two of its values always exceeds the value at their sum. The smallest mismatch anywhere on the chain is at the two lowest modes and equals 2.156e-4 — which the scan finds and which is the cube of pi over four times the cube of one more than the mode count, checked here on chains from eight masses to two hundred and fifty-six. That closed form is the whole of why this is a finite-chain problem: the mismatch falls as the cube of the length, so a long enough chain is arbitrarily close to resonant and the continuum limit is exactly resonant, which is where the solitary waves come from.

The condition three modes never meet

Whether two modes of a chain can hand energy to a third is arithmetic on the dispersion relation, and for a chain of masses the answer is never: a sine is concave, so the sum of two frequencies always exceeds the frequency of their sum. The smallest shortfall anywhere on a chain of N masses is π³/4(N+1)³ — never zero, and never far from it — and a shortfall turns a transfer into a beat.

The modes of a wire, counted. Nyquist's argument of 1928, with its count performed. Two resistors joined by a lossless line are in equilibrium, and the line is a one-dimensional cavity whose standing waves are spaced c/2L apart in frequency. Each mode has an electric and a magnetic energy, both quadratic, so equipartition gives it kT — the same half a kT per quadratic term that a heat capacity counts. The upper points are how many modes fall in a band of a hundred and thirty-seven megahertz, for five line lengths, from 17 on a thirteen-metre line to 8,558 on one of six kilometres. The lower points are the power each end therefore receives, as a fraction of kTΔf, and they converge on one: a longer line has proportionally more modes and takes proportionally longer to deliver them, so the length cancels. The longest line lands within 0.005 per cent and the shortest is 4.5 per cent low, because thirteen metres holds only seventeen whole modes in the band and the remainder is a real granularity rather than an error. What is left in the limit is kTΔf — a noise power with no resistance in it at all, and none of the line's properties either.

Half a kT in a piece of wire

Count the quadratic terms in a molecule's energy and equipartition gives a heat capacity. Count them on a transmission line instead and the same theorem gives a resistor's noise voltage — 4kTRΔf, with nothing in it about what the resistor is made of. A fifty-ohm input at room temperature says 0.91 nanovolts in every root hertz, and no design removes it.

Twice the kinetic energy, and what it equals. Twice the time-averaged kinetic energy of a bound orbit, divided by its time-averaged potential energy, against the power with which that potential depends on separation. Each point is a measurement: an eccentric orbit integrated for more than a hundred radial periods, with the two averages accumulated along it, and the radius checked to vary by at least a fifth so that the orbit is not trivially circular. The line is the exponent itself, and the points miss it by at most 6.6e-4. Two cases carry everything. At n = 2, a harmonic well, the two energies are equal — which is the ordinary equipartition statement, half a kT to the kinetic term and half a kT to the potential one. At n = −1, which is gravity and the Coulomb force, twice the kinetic energy equals minus the potential energy, so the total energy of a bound system is minus its kinetic energy. Nothing about temperature entered, and nothing about equilibrium: the relation holds for one orbit averaged over time as well as for a crowd averaged over members.

Weighing what cannot be put on a scale

Summed over every coordinate of a bound system, equipartition stops being a statement about temperature and becomes a relation between two averages: twice the kinetic energy equals n times the potential energy for a potential going as the nth power. For gravity that fixes a bound system's total energy from how fast its parts move — so a Doppler shift and an angular size return a mass, and for the Coma cluster the mass they return is fifty times the mass that shines.

A fridge with no work going into it, and its ceiling. How much heat a three-reservoir machine can lift out of a cold space per unit of heat supplied to drive it, against the temperature of the driving heat, for 3 cold temperatures and an ambient of 300 kelvin. No work enters or leaves: the machine takes heat in at the top, takes heat in at the bottom, and rejects the sum at ambient. That it can do anything at all is the surprise — the second law allows heat to be moved up a gradient provided a larger flow is moved down one, and the accounting is a single inequality in the three entropy flows. The ceiling is the product of two familiar expressions, and at 450 kelvin driving a 253-kelvin space it is 1.79. The marks are what real machines achieve, which is a fifth to a third of it — absorption refrigeration is not efficient and is chosen when the heat is free and the silence and the absence of moving parts are worth something.

A fridge with no work going into it

Every engine here so far turns heat into work or work into a heat flow. A machine exchanging heat with three reservoirs and doing no work at all can still move heat from cold to hot, and the ceiling on how much is the product of two Carnot expressions — an engine's efficiency times a fridge's coefficient of performance. A gas flame makes ice, and the accounting is one inequality in three entropy flows.

One dimensionless group between an engine and Carnot. The efficiency of a thermoelectric couple against the temperature of its hot side, with the cold side at 300 kelvin, for 4 values of the figure of merit, and the Carnot ceiling drawn above them. The expression has exactly one material quantity in it — the dimensionless group formed from the Seebeck coefficient squared, the electrical conductivity, the temperature and the thermal conductivity — and everything else is the two temperatures. At 600 kelvin, a figure of merit of one gives 10.8 per cent against a Carnot ceiling of 50.0, and a figure of merit of four gives 22.6. The approach to the ceiling is slow: every doubling of the group buys less than the last, so the difference between a good material and a perfect one is smaller than the difference between a poor material and a good one.

An engine with one number in it

A thermoelectric couple has no moving part and no working fluid, and its efficiency is the Carnot value multiplied by a factor containing exactly one dimensionless group of material properties. Sixty years of effort have moved that group from about one to about two, and the reason it is hard is that its three ingredients are not independent: raising the conductivity ruins the coefficient it is squared against, and the only lever that is really free is the heat the lattice carries.

Every reaction's free energy has its lowest point inside. An ideal reaction A ⇌ B at 298 K. Across: how far it has gone, from pure A towards pure B. Up: the Gibbs energy of the mixture per mole, relative to pure A. Left, for standard reaction Gibbs energies ΔG° of −4, 0, +4 kJ/mol: the dashed straight lines are what the energy would be if A and B did not mix, and the solid curves add the entropy of mixing them. For ΔG° = −4 kJ/mol the lowest point is at 83.4 per cent B; for ΔG° = 0 kJ/mol the lowest point is at 50.0 per cent B; for ΔG° = +4 kJ/mol the lowest point is at 16.6 per cent B. Right, magnified near pure A, a reaction with ΔG° = +10 kJ/mol, whose straight line climbs from the start and which looks as if it should not proceed at all: its curve first falls, to a minimum of −43 J/mol at 1.74 per cent B, because the mixing term falls infinitely steeply away from a pure end. Each minimum was found by search and sits where the ratio of B to A equals exp(−ΔG°/RT).

The reaction that cannot go all the way

Chemistry speaks of reactions that go to completion and reactions that do not happen, and at equilibrium there are neither. The reason is a logarithm. The free energy of a half-finished reaction contains the entropy of mixing, whose slope is infinite at both pure ends, so every reaction's lowest point lies strictly inside — and the slope of that free energy, the chemical potential, is to particles what temperature is to heat.

The light of a diode at room temperature is as bright as a surface thousands of kelvin hot. How many photons occupy each mode of the light, on a logarithmic scale, against photon energy. The lowest curve is the thermal glow of a 1.42 eV semiconductor at 300 K with no voltage across it. The solid curve above it is the same device with 1.3 V across it, emitting only above its gap. The dashed curve is a blackbody at 2571 K, the temperature whose light has the same occupation as the diode's at 1.472 eV, just above the gap. They cross there and nowhere else: the diode's occupation falls a factor e every 25.9 meV, as its lattice's temperature requires, and the blackbody's every 222 meV. No single temperature describes the diode's light. At each photon energy it has a brightness temperature, and that temperature is 300 K multiplied by ε/(ε − qV).

The glow that carries a voltage

Thermal radiation has no chemical potential, because walls make and destroy photons freely. A light-emitting diode is a body that glows at room temperature with a voltage written into its light — Planck's law with the voltage as the photons' chemical potential — which is why its light can be as bright as a surface thousands of kelvin hot, why at low voltage it can put out more light than the power it draws and cool itself doing so, and why a reverse voltage makes a surface look colder than it is.

How full a surface is, against the pressure of the gas above it. The fraction of adsorption sites occupied on a surface that binds a molecule with 0.35 eV, against the gas pressure on a logarithmic axis, at 250, 300, 350 K. Each curve is θ = p/(p + p½) — Langmuir's isotherm — and each is the same S-shape moved along the axis: on a logarithmic pressure scale the coverage is a function of ln p alone, shifted by the binding energy over kT. The half-full pressure is 5.58 Pa at 250 K, 132 Pa at 300 K, 1341 Pa at 350 K. Warming the surface does not change the shape; it slides the curve to higher pressure by a factor that is a Boltzmann factor, e^(ε/kT), which is how a binding energy is read off a family of isotherms.

The site that fills like an electron level

A patch of catalyst holding a gas molecule, a haem in a muscle cell holding oxygen, an energy level in a metal holding an electron: each is a place that can hold one particle or none, in contact with a reservoir it trades particles with. Each fills according to the same function of one variable — how far the reservoir's chemical potential lies above the energy of the site. Langmuir's isotherm and the Fermi function are not analogous. They are the same law, and haemoglobin works because its sites break it.

Electrons and holes in silicon against the donors added. The concentrations of free electrons and holes in silicon at 300 K against the concentration of donor atoms added, on logarithmic axes, solved from charge neutrality and the product law n·p = nᵢ². Below the intrinsic concentration of 10¹⁰ per cubic centimetre the donors change almost nothing. Above it the electrons follow the donors one for one and the holes fall in exact proportion: a hundred times more donors, a hundred times fewer holes. At 10¹⁶ donors the electrons outnumber the holes by 10¹² to one. The dashed line is √(n·p), which does not move at all across ten decades of doping — checked to a part in a billion.

The product doping cannot move

Add one phosphorus atom for every five million silicon atoms and the crystal's free electrons multiply a millionfold — and its holes fall a millionfold, although nothing was done to them. Add a drop of acid to water and the hydroxide ions vanish in exactly the proportion the hydrogen ions grow. Both are one law: when a reaction makes two species together, equilibrium fixes the sum of their chemical potentials, and so the product of their concentrations. Every semiconductor device is built on the minority carriers that law leaves behind.

Fifty-nine millivolts for every factor of ten. The voltage that holds a singly charged ion at one concentration on one side of a boundary and another concentration on the other, in equilibrium, against the ratio of the concentrations on a logarithmic axis: (kT/e) ln(c₁/c₂), from requiring the electrochemical potential to be the same on both sides. At 0 °C, 54.20 mV per decade; at 25 °C, 59.16 mV per decade; at 37 °C, 61.54 mV per decade; at 100 °C, 74.04 mV per decade. The voltage depends on the ratio and the temperature and on nothing about the ion, the solvent or the boundary: it is the entropy of dilution, expressed per unit charge.

The voltage that is a chemical potential

A battery, a nerve cell and a pH meter all turn a difference in chemical potential into a voltage, and the conversion is the same in each: fifty-nine millivolts for every factor of ten in concentration, at room temperature, whatever the ion. The voltage a battery holds as it discharges is the chemical potential of the atoms stored in its electrodes — and a battery whose voltage barely changes is one whose electrode is going through a phase transition.

The price of precision never falls below two. For a clock ticking forwards and backwards with a free-energy drop a per step, the product of its squared relative uncertainty and the entropy it has produced, in units of k: exactly a·coth(a/2), whatever the time elapsed. It approaches the bound of 2 only as the driving vanishes and grows as a once the driving is strong; the three clocks compared elsewhere sit at 2.04 (0.5 kT), 2.63 (2.0 kT), 8.01 (8.0 kT). No steady-state process of any design lies below the line at 2 — the thermodynamic uncertainty relation.

The precision a clock pays for in heat

Anything that counts — a molecular motor stepping along a filament, a chemical oscillator keeping time in a bacterium, a clock ticking on a wall — has to be driven, because a count that is not driven runs backwards as often as forwards. The heat it gives off buys it regularity, and there is an exact floor under the price: the squared relative uncertainty of the count, times the entropy produced, is never less than two in units of Boltzmann's constant.

Hydrogen comes apart far below its binding temperature. The fraction of hydrogen that is ionised in equilibrium, against temperature on a logarithmic axis, for densities of hydrogen nuclei of 10²³, 10¹⁸, 10¹², 10⁶ per cubic metre, from Saha's equation. Hydrogen's binding energy, 13.6 eV, corresponds to 157,803 K. It is half ionised at 13,980 K at 10²³ per cubic metre, 7,235 K at 10¹⁸ per cubic metre, 4,519 K at 10¹² per cubic metre, 3,270 K at 10⁶ per cubic metre — at a temperature where kT is between an eleventh and a fiftieth of the energy needed to free the electron. The thinner the gas, the cooler the temperature at which it comes apart, because a free electron in a thin gas has vastly more room to be in than a bound one.

The electron that leaves because there is room

It takes 13.6 electronvolts to pull the electron off a hydrogen atom, the energy of a particle at a temperature of 158,000 kelvin. Hydrogen is half ionised at the surface of a star at 14,000 kelvin, in a thin nebula at 7,000, and in the early universe at under 4,000 — a fortieth of its binding energy. The electron leaves not because it has enough energy but because a free electron has so much more room to be in than a bound one, and at equilibrium room counts as much as energy. Saha's equation is that trade written down, and it turned out to be the key to what stars are made of.

The melting curve that turns back, then lies flat. The pressure at which liquid and solid helium-3 coexist, against temperature from 0.9 millikelvin to 1 kelvin on a logarithmic axis, from the international provisional low-temperature scale PLTS-2000, which is defined by this curve. Above 315.2 mK the curve behaves like any other: a higher temperature needs a higher pressure to freeze the liquid. Below it the curve turns back — at its minimum, 2.93113 MPa, it takes a higher pressure to freeze colder helium — and it rises to 3.4391 MPa at 1 mK, where it lies almost flat, as the third law requires of every coexistence line. The kinks that mark the liquid becoming superfluid and the solid's spins ordering sit on that flat end, at 2.444, 1.896 and 0.902 mK.

The melting curve that has to arrive flat

The slope of every line on which two phases coexist is their difference in entropy divided by their difference in volume. The third law says that entropy differences vanish at absolute zero, so every such line — melting, boiling, a superconductor's critical field — must reach the bottom of the temperature scale lying flat. Helium-3 reaches it in the strangest way available: below a third of a kelvin its solid is more disordered than its liquid, the melting curve runs backwards, and squeezing the liquid into solid absorbs heat. That backwards curve was the refrigerator in which superfluid helium-3 was found, and it is now the definition of temperature below one kelvin.

The same rate for every size of molecule. The diffusion-limited rate constant for two equal spherical molecules meeting in water at 25 °C, against their radius from 0.1 to 100 nm: 8kT/3η = 7.43·10⁹ per molar per second, flat. The dashed line is what a larger target would give if its partner kept the mobility of a 0.3 nm molecule — 7.4·10¹¹ at 30 nm — the intuition that a bigger molecule is easier to hit. Stokes and Einstein's law takes it away again: a molecule twice as big is twice as easy to hit and diffuses half as fast, and the two cancel exactly. Every reaction whose partners react at the first touch runs at close to this one number, whether the partners are ions, small organic molecules or proteins.

The speed limit that does not care how big the molecules are

Two molecules in water cannot react faster than they can find each other, and how fast they find each other is set by diffusion. A bigger molecule is an easier target, so it seems a reaction between big molecules should be faster. It is not: a molecule twice as big is twice as easy to hit and wanders half as fast, and the two cancel exactly. Every reaction that happens at the first touch runs at the same rate, about seven billion per molar per second in water, whether the partners are ions or proteins — a speed limit written by the viscosity of water and the temperature, and nothing else.

The same thermal motion, spread three ways. How the thermal jiggling of a spring-mounted mass is distributed over frequency, for quality factors 3, 30, 300, with frequency in units of the resonance and the spectrum in units that make the area under each curve its mean-square displacement. All three curves enclose exactly the same area, kT/k — 1.000, 1.000, 1.000 of it in the units drawn — because equipartition fixes the energy in the spring and not how fast it changes. What the damping chooses is where that motion goes: a high Q concentrates it in a peak 100 times taller than the low-Q one, and leaves the rest of the spectrum 100 times quieter at 10 times the resonance frequency.

The noise a high Q moves out of the way

Anything held by a spring at a temperature jiggles, and equipartition says exactly how much: the mean square of its displacement is kT divided by the spring's stiffness, whatever the spring is made of and however it is damped. That looks like a floor no measurement can go below. It is a floor on the total, and the total can be moved. An oscillator with little damping puts almost all of its thermal motion into a narrow peak at its own frequency and leaves every other frequency quiet — which is why the mirrors of gravitational-wave detectors hang from fibres of glass that ring for minutes.

Hydrogen ionised by heat at one end and by crowding at the other. The ionised fraction of hydrogen against the density of nuclei, at 5,000 K, 10,000 K, 20,000 K, from Saha's equation with the binding energy lowered by the neighbouring charges (solid) and without the lowering (dashed). Thin gas is ionised by the room a free electron gains; compression takes the room away and the gas recombines, as the dashed curves continue to do for ever. The solid curves turn round: at 5,000 K the gas is least ionised, 2.5·10⁻⁴ per cent, near 10²⁶ per cubic metre; at 10,000 K the gas is least ionised, 0.41 per cent, near 10²⁷ per cubic metre; at 20,000 K the gas is least ionised, 12 per cent, near 10²⁸ per cubic metre, and by 6·10²⁸ per cubic metre — 0.10 g/cm³ — every curve is fully ionised, because the ground state no longer exists. That density does not depend on the temperature at all.

The pressure that ionises

Squeeze a hot gas of hydrogen and Saha's equation says it recombines: a free electron has less room to be in, and room was what paid for its freedom. Squeeze it far enough and it comes apart again, completely, at any temperature — not because the room has grown but because the atom has nowhere left to be. Crowded to within a few atomic radii of each other, hydrogen atoms have no bound state, and their electrons, freed by pressure alone, form a metal. Deep inside Jupiter most of the planet's hydrogen is in that state.

The proton fraction a cold star settles at. The fraction of nucleons that are protons in cold matter in beta equilibrium, where the neutron's chemical potential equals the proton's plus the electron's, against density on logarithmic axes, for free neutrons, protons and electrons with no nuclear forces. Below 1.22·10⁷ g/cm³ the electrons cannot pay the 1.29 MeV difference between a neutron and a proton and the matter is all protons. Above it the fraction collapses, to 0.014 per cent at its lowest near 7.9·10¹¹ g/cm³, and at nuclear density, 2.66·10¹⁴ g/cm³, it is 0.49 per cent. At higher densities it climbs slowly towards one ninth, 11 per cent, which it never reaches. The dashed part lies where real matter is made of nuclei rather than free nucleons.

The protons a star cannot afford

A free neutron decays in fifteen minutes into a proton, an electron and an antineutrino. Inside a neutron star it cannot, and the star is made of neutrons because of it. Beta decay is a reaction like any other, and at equilibrium the neutron's chemical potential must equal the proton's plus the electron's. In a crowded star an electron can only be added at the top of a sea tens of megaelectronvolts deep. So the matter settles where protons are rare: half a per cent of the nucleons at nuclear density, if nucleons were free, and never more than one in nine at any density.

The gas that has to stop expanding. The thermal expansion coefficient at constant pressure times the temperature, αT, against temperature over the Fermi temperature, on logarithmic axes, for an ideal classical gas and an ideal gas of spin-half fermions at the same density. For the classical gas αT is exactly one at every temperature, so α = 1/T grows without limit as the gas is cooled, which the third law forbids: by a Maxwell relation α is minus the pressure derivative of the entropy per unit volume, and that must vanish at absolute zero. The Fermi gas agrees with the classical one when hot and departs from it near the Fermi temperature; at a tenth of it αT is 0.048, and at a hundredth 0.002, falling as π²/2 times the square of the temperature over the Fermi temperature. Its expansion vanishes as the law requires.

The two laws absolute zero will not allow

The ideal-gas law and Curie's law are the first two equations of state most people meet. The third law of thermodynamics says both must be wrong at low temperature, and it says so without any model of what replaces them. A Maxwell relation turns every derivative of the entropy into something measurable with a ruler or a magnetometer: the thermal expansion, the slope of the magnetisation with temperature. If the entropy stops depending on pressure and field at absolute zero, those must vanish. The classical gas expands as 1/T and Curie's magnetisation climbs as 1/T, and both slopes grow without limit instead.

What one trajectory says, and what the crowd says. For the walk with α = 0.5: the time-averaged mean-square displacement of each of 12 single trajectories, each watched for 10000 time units — the average of the squared displacement over every window of a given lag along that one path — against the lag, on logarithmic axes (thin lines), and the mean-square displacement of an ensemble of 800 walkers at the same times (thick). The ensemble spreads as t^0.5. Every single trajectory's time average grows nearly in proportion to the lag — mean slope 0.82, where the ensemble's is 0.5 — as if it were ordinary diffusion, and the trajectories disagree with one another about the diffusion coefficient by more than a factor of ten. A time average along one path and an average over many paths, which for ordinary diffusion are the same thing, here give different laws. That is weak ergodicity breaking, and it is what a single-molecule tracking experiment has to be read against.

The diffusion that slows the longer it is watched

A random walker whose steps come at irregular times still diffuses normally, provided the waits have a finite average. Let the waits have no average — mostly short, occasionally enormous, as they are for a molecule that keeps getting stuck — and three things change at once. The spread grows more slowly than time. A diffusion coefficient measured by following one particle depends on how long it was followed. And following one particle for a long time no longer tells what many particles do: averages along a path and averages over a crowd give different laws. Tracking single molecules inside living cells runs straight into all three.

The melting line of water, through five ices. The pressure at which water freezes, in megapascals on a logarithmic axis, against temperature, along each of the five solid phases that meet the liquid, from the international standard equations. Along ordinary ice, Ih, the melting pressure rises as the temperature falls, from 611 Pa at 273.16 K to 209 MPa at 251.165 K — the backwards slope of the familiar melting line. There the liquid meets ice III, and from then on every ice melts at a higher pressure as the temperature rises: ice V, ice VI, and ice VII, which reaches 20.6 GPa at 715 K. The dots are the four triple points where two ices and the liquid coexist. The lowest temperature at which liquid water can exist in equilibrium, 251.165 K, is at the first of them.

The melting line that turns round

Squeeze water and it freezes at a lower temperature — the backwards slope of its melting line, and the reason ice floats. Keep squeezing and the slope reverses. Past two thousand atmospheres ordinary ice gives way to a denser ice, then another, and from there on every ice melts at a higher temperature under more pressure, exactly as a normal substance's does. At twenty thousand atmospheres water freezes into an ice that is still solid at a hundred degrees Celsius; at two hundred thousand, one still solid at four hundred. The famous anomaly belongs to one ice out of five, and the same relation that explained it explains its end.

The boiling curve, which turns back. Heat flux from a horizontal surface into water boiling at one atmosphere, on a logarithmic axis, against how far the surface is above the boiling point, also logarithmic, from standard correlations. At a few degrees heat leaves by convection; from about 5 K bubbles nucleate and the flux climbs as the cube of the superheat, to the critical heat flux of 1.11 MW/m² at 20 K. Beyond it the flux FALLS as the surface gets hotter — the transition regime, drawn dashed as the interpolation it is — down to 19 kW/m² at 67 K, the Leidenfrost point, where a continuous film of vapour insulates the surface. Past that, film boiling and radiation let the flux rise again, slowly: 310 kW/m² at 1000 K above boiling. Between the crisis and the Leidenfrost point, a surface 3 times hotter passes 58 times less heat.

The boiling curve that turns back

Heat a surface under water and more heat flows from it the hotter it gets — up to a point. Past about twenty degrees above boiling the bubbles crowd together, merge, and begin to blanket the surface with vapour, and from then on a hotter surface passes less heat, not more, until by a few hundred degrees above boiling it is wrapped in an insulating film and passes fifty times less than it did at its best. Shigeo Nukiyama found the turn in 1934 by melting wires. It is why a water drop skates for a minute on a very hot pan but vanishes in two seconds on a merely hot one, and why every boiler, steam generator and water-cooled reactor is designed around a number called the critical heat flux.

The race to equilibrium that the hotter start wins. How far each system is from equilibrium with the bath — the relative entropy of its distribution from the bath's, in nats, on a logarithmic axis — against time since it was placed in the bath at 0.30, also logarithmic, for a start at 0.50 (warm, dashed) and one at 0.901 (hot, solid). The hot system starts further away: 1.243 against 0.254. Both close in quickly as their left wells equilibrate internally. Then the warm system must wait for probability to leak between the wells, and stalls, while the hot one, whose wells already had the right shares, carries on falling. The hot curve overtakes the warm one at a time of about 6.9, and from then on the hotter start is always nearer equilibrium.

The hotter start that cools first

Put two identical systems, one warm and one hot, into the same cold bath, and common sense says the warm one reaches the bath's temperature first, since the hot one has further to go. The claim that hot water sometimes freezes before warm water — the Mpemba effect — has been argued over for sixty years and does not survive careful measurement. But a cleaner version of it is real and has been observed: in systems whose relaxation has a slow route and a fast one, a hotter start can miss the slow route entirely and overtake a cooler start, exponentially. Whether a system cools quickly depends not on how far it is from equilibrium but on how much of the slowest way back it needs.

The temperature through growing lake ice. Temperature against depth below the ice surface in a lake under air 10 K below freezing, after 1, 4 and 16 days, from Stefan's law with the ice's top held at the air's temperature. The ice is 11.1, 22.3, 44.6 cm thick; within it the temperature rises in a straight line from −10 °C at the top to 0 °C at the bottom, and below it the water is at 0 °C throughout. Ice forms only at the bottom, where water meets ice at the freezing point, and the latent heat it releases there has to be carried up through all the ice above by conduction. The straight lines get shallower as the ice thickens: the same temperature difference drives heat across a longer path, so less heat escapes and the ice grows more slowly. Sixteen times as long gives four times the thickness.

The ice that grows more slowly the thicker it gets

A lake under air ten degrees below freezing grows eleven centimetres of ice on the first day and only another eighteen in the following six. The ice forms at its underside, and every gram of it releases latent heat that must be carried up through the ice already there before more can form: the ice is its own insulation. The thickness therefore grows as the square root of time, and the rule, Stefan's, depends on temperature and time only through their product — the freezing degree-days that lake-ice forecasters count. A few centimetres of snow on top slows it more than any amount of cold speeds it up.

How well a diamond-like crystal conducts heat, from 1 K to 1000 K. The thermal conductivity of a model crystal with diamond's sound speed and Debye temperature, from kinetic theory applied to its phonons, against temperature on logarithmic axes: natural carbon in a sample 1 mm across, isotopically enriched carbon, and natural carbon in a sample 0.1 mm across; copper's room-temperature 400 W/m K is marked for comparison. The model is calibrated to natural diamond's 2,200 W/m K at room temperature. At low temperature the conductivity rises as the cube of the temperature, limited by the sample's walls; at high temperature it falls as umklapp collisions shorten the phonons' paths. Between, it peaks: 3.7·10⁴ W/m K at 100 K for the natural crystal; 7.8·10⁴ W/m K at 94 K for the enriched crystal; 1.2·10⁴ W/m K at 126 K for the thin natural crystal. The model gets the shape — the rise as the cube of the temperature, a peak near 100 K, the fall — but overestimates the heights, because it leaves out the phonon collisions that conserve momentum and every defect but isotopes; measured enriched diamond peaks near 4 × 10⁴ W/m K. An insulator with no free electrons still conducts heat five times better than copper at room temperature.

The gas of sound that carries a diamond's heat

Diamond has no free electrons and conducts heat five times better than copper. What carries the heat is a gas — not of molecules but of quantised sound waves, phonons — and the kinetic theory that gives the conductivity of air gives the conductivity of a crystal with the same formula: a third of heat capacity times speed times mean free path. The gas is strange in two ways. At low temperature its particles cross the whole crystal without meeting anything, so the conductivity depends on how big the sample is. And at room temperature one atom in a hundred being a neutron heavier is enough to scatter its fastest-oscillating members.

How a solute moves both ends of the liquid range. The molar chemical potential of water against temperature near its freezing point (left) and its boiling point (right): liquid (blue), ice (green), vapour (red), each a line whose slope is minus its entropy, drawn relative to the liquid at the transition. A solute at mole fraction 0.05 — about 2.9 mol/kg — lowers the liquid's line by RTx (dashed) and leaves the ice and the vapour unchanged. The crossings move: the freezing point down by RTx divided by the entropy of melting, 5.16 K, the boiling point up by RTx divided by the entropy of boiling, 1.42 K. The same vertical shift moves the freezing crossing 3.6 times as far, because ice's line is nearly parallel to the liquid's — melting changes the entropy by 22.0 J/mol K — while the vapour's line crosses steeply, with 109 J/mol K between them. Lowering the liquid makes it the winner over a wider range of temperature at both ends: the solution is liquid from below 0 °C to above 100 °C.

The degree a spoonful of solute buys

Dissolve anything in water and the water freezes lower and boils higher. The two shifts come from one cause — the solute lowers the liquid water's chemical potential, so the liquid wins against ice and against vapour over a wider range — but they are not the same size. For every mole of particles per kilogram the freezing point falls 1.86 degrees and the boiling point rises only 0.51, three and a half times less, and the reason is in the entropies: ice and water are so alike that a small push moves their crossing a long way. The same arithmetic says salting pasta water does nothing a cook could notice, weighed molecules for a century of chemistry, and in 1887 showed that salts in water are made of ions.

How supersaturated a cloud is, to an ice crystal. The supersaturation over ice of air that is just saturated over liquid water, eᵥᵥ/eᵢ − 1, in per cent, against temperature; dashed, the same computed from the chemical potentials, exp(LΔT/RT·Tₘ) − 1, with L the latent heat of melting per mole, ΔT the supercooling and Tₘ the melting point. A supercooled cloud is in balance with its droplets, so to any ice crystal inside it the air is 5.1 per cent supersaturated at −5 °C, 10.5 at −10, 22.0 at −20 and 41.2 at −35. Over the droplets themselves it is exactly saturated, so they neither grow nor shrink by vapour. The excess is the chemical-potential difference between supercooled water and ice, 217 J per mole at −10 °C, expressed as a pressure ratio.

The ice that grows by stealing from the droplets

Most clouds colder than freezing are not made of ice. They are made of droplets of liquid water, some tens of degrees below their freezing point, waiting for something to start them freezing. When a few ice crystals do appear among them, the crystals grow and the droplets shrink, without touching, because supercooled water has a higher chemical potential than ice and so a higher vapour pressure: air balanced with the droplets is ten per cent supersaturated, to an ice crystal, at minus ten degrees. A few crystals among a hundred thousand droplets each take everything, grow heavy enough to fall, and that is how most of the rain outside the tropics begins — as snow.

The highest note the air can carry, with height. The frequency at which a sound wave in the air loses its amplitude by a factor e in a single wavelength, in hertz on a logarithmic scale, against height, with the band of human hearing shaded. Because viscosity does not depend on pressure, the collision time grows as the pressure falls, and the ceiling falls with it. It is about 212 MHz at sea level, 0.17 MHz at 50 km and 3 kHz at 80 km. It passes the top of human hearing, 20 kHz, at about 68 km, and the bottom, 20 Hz, at about 110 km; at 100 km it is 92 Hz. Above about a hundred kilometres the air cannot carry anything a person could hear, at any loudness: the molecules are too far apart to pass a compression on before it has dispersed.

The note too high for thin air

Sound in a gas is a compression handed from molecule to molecule by collisions, and it can only be handed on if the molecules collide many times in each period of the wave. At sea level they collide seven thousand million times a second, and the limit is far above anything audible. But a gas's viscosity does not depend on its pressure, so as the air thins with height the time between collisions grows in proportion, and the highest note the air can carry falls with it: two hundred megahertz at the ground, three kilohertz at eighty kilometres, below twenty hertz above a hundred and ten. The upper atmosphere is not silent for want of loudness. It is silent because its molecules are too far apart to pass a sound along.

Why silt heaves and sand and clay do not. The relative rate at which water is drawn to a growing ice lens, against grain diameter on a logarithmic scale, for three loads on the soil above it, in a capillary model: the rate is the soil's permeability, which grows as the grain diameter squared, times the suction the lens can sustain less the load. The suction is capped both by the supercooling at the lens, here 2 K, and by the Gibbs–Thomson pressure at which ice would break into the pores. Coarse soil has open pores that ice simply invades, so the suction never exceeds the load: under 10 kPa, half a metre of soil, nothing heaves above 58 μm. Fine clay can sustain the full suction but passes water too slowly to feed it. Between lies a fastest grain size, 29 μm under 10 kPa, 10 μm under 30 kPa and 2.9 μm under 100 kPa: silt, the soil engineers have long called frost-susceptible. Real heave pressures exceed this model's, and the fit is to the shape, not the scale.

The frost that lifts by drinking

Frozen ground rises in winter, sometimes by tens of centimetres, and the usual explanation — water expands by nine per cent when it freezes — is wrong: soil soaked in benzene, which shrinks as it freezes, heaves too. The lift comes from layers of clear ice that grow by pulling liquid water up through the soil from below. Ice and water touching each other need not be at the same pressure, and when the water is free to flow, ice a degree below its melting point can push with 1.12 megapascals while it drinks — twelve times less than freezing water in a sealed pipe, and enough to lift a house at a tenth of a degree.

The distance at which the air takes over. Evaporation of water at 20 °C into dry air, in kilograms per square metre per second, against the thickness of the layer the vapour must cross before it is carried away, both on logarithmic scales, with the kinetic and diffusive resistances added in series: 1/u + δ/D, where u = √(RT/2πM) is the speed at which the kinetic flux carries vapour off and D the diffusion coefficient. Dashed, the kinetic ceiling. At one atmosphere the two resistances are equal at δ = D/u = 0.17 μm, about two and a half mean free paths of the air; a layer much thinner than that leaves the surface at its ceiling, and a layer much thicker throttles it in proportion to its thickness. At a tenth of an atmosphere diffusion is ten times faster and the crossover moves out to 1.7 μm. Only in vacuum, or across gaps shorter than a mean free path, does the liquid itself set the rate.

The puddle that could dry in a second

A liquid in equilibrium with its vapour loses molecules exactly as fast as the vapour delivers them, and kinetic theory says how fast that is. For water at room temperature it is two and a half kilograms per square metre every second — a puddle two and a half millimetres deep gone in a second. Puddles take hours because almost every molecule that leaves is knocked back by the air within a fraction of a micrometre. The ceiling binds only where the vapour can escape in a single free flight: in vacuum, where water freezes itself, in the smallest droplets, and on the filament of a lamp.

Latent heat of boiling against boiling point. The heat needed to boil one mole of a liquid at its normal boiling point, against that boiling point, both on logarithmic axes, for twenty-one substances from helium to sodium. The line is ΔH = 87.2 J/(mol K) × Tb, the mean entropy of boiling of the ordinary liquids boiling between 250 and 700 K, which scatter about it by 3.3 per cent; the ordinary liquids as a whole (blue), from neon at 27 K to sodium at 1,156 K, by 10.0 per cent, the coldest lying low. Water, ammonia and the alcohols (red) need more heat than the line predicts, acetic acid (green) less; the quantum liquids helium and hydrogen (violet) fall far below.

The entropy every liquid gains on boiling

Divide the heat it takes to boil a liquid by the temperature at which it boils, and for benzene, chloroform, acetone, hexane, ether and carbon tetrachloride the answer comes out within a few per cent of the same number — about ten and a half times the gas constant. Frederick Trouton noticed it in 1884. It is a statement about entropy, and once it is read that way the liquids it fails for become the interesting ones: water and the alcohols, whose hydrogen bonds hold the liquid in order; acetic acid, whose vapour is made of pairs; and helium, whose liquid is already as disordered as quantum mechanics allows.

A cluster grown by random walkers. 3000 particles, each released far away to wander at random on a square grid until it first touches the cluster and sticks there, starting from one seed at the centre. Colour marks the order of arrival: the first third (violet), the second (blue), the last (red). The cluster is 220 grid steps across, branched and open, and almost every late arrival is at the outer tips: walkers wander into the outer branches long before they can find their way into the gaps between them. The gaps never fill.

The cluster that grows where the walkers arrive

Let particles wander in one at a time from far away, each sticking where it first touches a growing clump. The rule has nothing in it about shape, and the clump that results is a lacy, branched thing whose mass grows more slowly than its area, so it gets emptier the bigger it gets. The reason is that a random walker reaches the outermost tips long before it can find its way into a gap between branches. The same arithmetic — the probability that a walker arrives is a solution of Laplace's equation — shapes electrodeposits, lightning-like discharges in plastic, fingers of water pushed into oil, and colonies of bacteria short of food.

A heat capacity that goes up and comes down. The heat capacity per system, in units of Boltzmann's constant, against temperature in units of the level spacing over k: for a system with two levels (red), Schottky's hump x²eˣ/(1 + eˣ)², x = Δ/kT, and for a harmonic oscillator with the same spacing (blue), which climbs to one k and stays there. Both are frozen out when kT is well below the gap. The oscillator then fills its unending stack of levels and reaches equipartition's value. The two-level system peaks at 0.439k when kT = 0.417Δ, and then, with both levels nearly equally occupied and nothing higher to go to, it stops absorbing heat: its capacity falls back towards zero as 1/T².

The heat capacity that goes up and comes down

Equipartition promises every way a system can move a fixed share of heat, and for a harmonic oscillator the promise is kept: once the temperature is high enough, it holds one k per oscillator and goes on holding it. A system with only two levels cannot keep it. It absorbs heat while its upper level is filling, then stops — the two levels equally full and nowhere higher to go — so its heat capacity rises to a peak and falls back to zero. The hump, named after Walter Schottky, is the signature of a degree of freedom with a ceiling, and the area under it counts how many states were there.

Where skin and the thing it touches meet. Temperature across the contact between skin at 33 °C (left) and a block at 20 °C (right), against distance from the contact in millimetres, at 0.1, 1 and 10 s after touching, for steel (blue) and pine wood (red). Each pair of curves is two error functions meeting at the contact, and the meeting temperature does not change with time: 21.0 °C against steel and 30.3 °C against wood. The nerve endings a fraction of a millimetre under the skin report something close to that interface temperature, so the steel feels like an object at 21 °C and the wood like one at 30 — although both are at 20. The profiles spread as the square root of time, much further into the steel than into the wood.

The metal that feels colder than the wood

A steel railing and a wooden bench in the same shade are at the same temperature, and the railing feels much colder. The usual explanation, that metal conducts heat better, is half the answer, and it gets water wrong. When two bodies touch, the surface between them jumps at once to a temperature that then stays fixed, a weighted average of the two in which each body's weight is the square root of its conductivity times its heat capacity per volume. Skin at 33 °C touching steel at 20 °C meets it at 21; touching wood, at 30. The nerves report the meeting temperature, not the object's.

How far below the air a surface facing the sky falls. How much colder than the air a surface facing the night sky becomes, when insulated from below and with air at 5 °C, against how strongly the air can warm it, measured by the convection coefficient (3 to 5 W/m²K in still air, 20 or more in a breeze), for sky emissivities of 0.72, 0.82, 0.97 — clear and dry, clear and humid, and overcast. Under a clear dry sky in still air the surface sits 12.3 °C below the air; a breeze of 20 W/m²K halves that and more; under cloud the drop is 1.3 °C. The surface is radiating to something colder than the air: the upper atmosphere and, through the window near ten micrometres, space.

The cold reservoir overhead

On a clear, still night a car roof can grow frost while the air around it is several degrees above freezing. Nothing in the air is cold enough to freeze water. The roof is radiating its heat upward to a sky whose radiation back is that of something twenty degrees colder than the air — part of it the cold upper atmosphere, part of it space itself, seen through a window in the infrared where the air is transparent. The clear sky is a cold reservoir available everywhere, and it can be used: to make frost, to cool buildings below the air temperature even in sunshine, and in principle to run an engine at night.

How much of supercooled water freezes in an instant. The fraction of a body of supercooled water that freezes at once when it is nucleated, against how far below 0 °C it had been cooled: the ice that forms releases latent heat until the whole of it is back at 0 °C, and freezing stops there until heat is taken away. Solid: with water's measured heat capacity, which rises steeply below about −20 °C; dashed: with the heat capacity it has at 0 °C. At −15 °C, 19.1 per cent freezes at once; at −20 °C, 25.7 per cent; at −38 °C, the coldest liquid water reaches before it freezes by itself, 51 per cent. The rest stays liquid at 0 °C. Water would have to be cooled about 80 degrees for all of it to freeze at once, and no liquid water survives that long.

The supercooled water that only partly freezes

Water can be cooled well below 0 °C and stay liquid, and when a crystal finally starts, the whole bottle seems to turn to ice at once. It does not. Freezing releases heat, and with nowhere for that heat to go in the instant of freezing, it warms everything back to 0 °C — at which point freezing must stop. Water cooled to −15 °C freezes about a fifth at once; to −38 °C, the coldest liquid water gets, about half. The rest has to wait for the surroundings to take the heat away, and a hand warmer is the same arithmetic run above body temperature.

Insulating a thin wire makes it lose more heat. The heat lost per metre from a wire 2 mm across held 40 K above the air, against the outer radius of the PVC insulation round it (k = 0.19 W/m K), with the surface losing 10 W/m²K to still air by convection and radiation, on a logarithmic radius axis; dashed, the bare wire. The first layers of insulation increase the loss, from 2.51 to a maximum of 12.11 W/m at an outer radius of 19 mm — the critical radius k/h — because they add surface faster than they add resistance. Beyond the critical radius the loss falls again, but so slowly — the resistance of the insulation grows only as the logarithm of its radius — that it does not get back down to the bare wire's until the PVC is 178 km thick. Any sheath a cable could carry makes it lose more heat, not less.

The insulation that makes a wire lose more heat

Wrap a hot water pipe in glass wool and it loses less heat. Wrap a thin electrical wire in PVC and it loses more — four times more, at the best thickness — and runs cooler for it. The insulation adds a resistance to heat flowing through it, but it also adds surface for the heat to leave from, and around a thin enough cylinder the surface wins. The crossover is a radius, k/h, set by the insulation and the air; below it, more insulation means more heat lost, and a small enough ball cannot be insulated by anything wrapped round it.

Each foil divides the heat by one more. The radiant heat flowing from a wall at 293 K to one at 77 K — room temperature to liquid nitrogen — in watts per square metre on a logarithmic axis, against the number of thin shields placed between them, for walls and shields all black-painted, dull metal and aluminised film. With nothing between them, black walls exchange 340 W/m²; aluminised walls 6.3. Each shield adds one more equal resistance in series, so N shields divide the flux by N + 1: thirty aluminised foils bring it to 204 milliwatts per square metre. Between black walls a single aluminised foil passes 6.2 W/m², less than thirty black shields do, 11.0 W/m², because a shield's resistance is 2/ε − 1 and a shiny one is worth 54 black.

The thin foils that stop radiant heat

In a vacuum, heat can cross a gap only as radiation, and a single sheet of foil a few hundredths of a millimetre thick, hung in the gap and touching nothing, halves it. Two sheets cut it to a third, thirty to a thirty-first. The foil does not need to be thick or cold or clever; it needs only to be shiny, because a sheet that absorbs little also emits little, and each one becomes a radiator in its own right, floating at a temperature between its neighbours. The golden blankets on spacecraft and the silvering inside a vacuum flask are the same arithmetic of resistances in series.

The fold beyond which no heap is steady. The centre's steady temperature rise, in units of RTₐ²/E, against the Frank-Kamenetskii parameter δ, for a slab, a long cylinder and a sphere, found by integrating the steady heat equation outward from the centre: each shape has two steady states for every δ below a critical value — a cool, stable one and a hot, unstable one — and none above it. The fold is at δ = 0.878, 2.000, 3.322 for the slab, cylinder and sphere: a sphere tolerates nearly four times a slab's δ at the same half-width, because its surface is larger in proportion to its volume and heat from its centre has more ways out. A real heap, which is neither, has a critical δ between, and fixing it for a given shape is the first step in calculating the size at which a store becomes unsafe.

The heap that sets itself alight

A haystack, a coal stockpile, a bin of oily rags or a silo of damp grain can catch fire with no spark at all. A slow reaction inside makes a little heat; the heap holds it in; the warmer middle reacts faster and makes more. Whether that ends in a warm heap or a burning one is decided by a single number built from the heap's size, the air's temperature and the reaction's speed — and past one value of it there is no steady temperature at all. The critical size shrinks steeply as the weather warms, and the heap gives almost no warning before it goes.

The period of a ball bouncing on a gas. The period of a 16 mm steel ball (16.8 g) oscillating in a vertical tube that fits it on top of a flask, against the flask's volume on a logarithmic axis, for helium, air and carbon dioxide at room temperature, if every compression is adiabatic: T = 2π√(mV/γPA²). On a 10 litre flask the ball swings with a period of 0.983 s on helium, 1.072 s on air and 1.117 s on carbon dioxide. The gas is the spring, and its stiffness is γ times the pressure: the gas whose molecules have fewest ways to store energy is the stiffest, because less of the work done in compressing it goes into motions that do not push on the walls.

The ball that weighs how a gas stores heat

Drop a steel ball into a glass tube that fits it closely, on top of a flask of gas, and it does not fall to the bottom. It bounces, up and down, on the gas trapped beneath it, with a period of about a second. The gas is a spring, and its stiffness is the pressure multiplied by a number that counts how many ways the gas's molecules can hold energy. So a stopwatch on a bouncing ball measures something that otherwise needs a calorimeter: whether the molecules can rotate, whether they can vibrate, and how far a vibration has woken up at the temperature of the room.

A race between the flame and the end gas. The Livengood–Wu integral of the end gas — the fraction of its ignition delay used up — against time from the spark, for 95-octane fuel at compression ratios 9, 11 and 13; the dotted line marks when nine-tenths of the charge has burned. If the integral reaches one before then, the end gas ignites by itself while a sizeable part of the charge is still unburned, and knocks. At a ratio of 9 the flame wins easily: the integral ends at 0.90. At 11 it reaches one when 92 per cent of the charge has burned, just short of knock; at 13, when only 78 per cent has, and the rest detonates. The integral climbs fastest at the end, where the end gas is hottest and most compressed, which is why knock is heard late in the burn.

The squeeze a petrol engine is not allowed

Every textbook gives a petrol engine's ideal efficiency as one minus the compression ratio raised to a small negative power, and every reader can see that the efficiency keeps climbing as the ratio rises. Real petrol engines stop near ten or eleven. Nothing in the thermodynamics stops them. The fuel does: the charge being compressed is already fuel and air, and the burning gas behind the flame squeezes the rest until it ignites on its own, in a hammer-blow that wrecks pistons. An octane number buys nothing but time against that, and a diesel escapes it by compressing air and adding the fuel only when it wants ignition.

A ring of balls that forgets, and then remembers. Kac's ring: 400 balls on 400 sites round a circle, all white at the start, with 52 of the edges between sites marked; each tick every ball moves one site and changes colour if it crosses a marked edge. The solid line is the excess of white over black, as a fraction of all the balls, against ticks; the dashed line is Boltzmann's prediction (1 − 2μ)^t with μ = 0.1, the fraction of edges marked. For the first few dozen ticks the ring follows Boltzmann's decay exactly as an irreversible process would: the colours mix and the excess falls to zero, then wanders by a few hundredths. At tick 400 every ball has gone once round, through all 52 markers — an even number of colour changes — and the ring is all white again. A rule with no randomness in it, run forward, looked irreversible for 400 ticks and then undid itself.

The ring that forgets and then remembers

Put a few hundred balls on a ring, all white, and mark a handful of the gaps between them. Each tick, every ball moves one place round and changes colour if it crosses a mark. Nothing is random and everything can be run backwards. Yet the colours mix, half and half, exactly as Boltzmann's equation for a gas says they should, and the entropy rises to its maximum and stays there — until, after one full turn, every ball is white again. Mark Kac built the model in 1956 to show where the arrow of time in Boltzmann's equation comes from, and what it quietly assumes.

Bands of precipitate laid down in a tube of gel. A tube of gel holding one salt, B, with a second salt, A, diffusing in from its left end, computed at the end of the run: the shading shows where solid precipitate has formed, darker for more. The product of their reaction precipitates only where its concentration first passes a nucleation threshold, 0.5 in the units used, and then wherever solid is already present. The solid is not spread smoothly. It forms in 23 sharp bands, close together near the entrance and further and further apart down the tube — the last two are 26.0 units apart, against 0.8 between the first two. Each band holds more precipitate than the last, in proportion to its distance. Nothing in the gel was patterned to begin with.

The rings a reaction lays down in a gel

Let one salt soak into a tube of jelly that holds another, and where they meet a solid should form, spreading smoothly behind the advancing front. It does not. It forms in sharp, separate bands, close together near the entrance and further and further apart, each one a fixed multiple of the distance of the last. Raphael Liesegang saw the rings in 1896 and they have been found since in rocks, in agate, in gallstones and in microchips. Nothing in the jelly is patterned to begin with. The pattern comes from a solid that will not start until the solution is overloaded, and that then empties the solution around it.

Salt that takes in water from the air at three-quarters humidity. The water held by a crystal of common salt, in grams per gram of salt, against the relative humidity of the air around it at 25 °C, with the humidity rising (solid) and falling (dashed). Rising, the crystal stays dry until 75.2 per cent and then dissolves at once into a drop holding 2.79 g of water per gram — the saturated solution — which goes on drinking as the air gets damper: 4.24 g at 85 per cent, 11.5 g at 95. Falling, the drop does not crystallise at 75 per cent. It sheds water along the same curve, becomes supersaturated, and stays liquid down to about 45 per cent before a crystal finally nucleates and it dries out at once.

The salt that dissolves in the air's own water

Leave a dish of salt in a bathroom and it does not just clump; on a humid enough day it turns into a puddle, though no water was poured on it. Above 75 per cent humidity, common salt pulls water out of the air until it has dissolved itself. The number is not a property of salt's attraction for water in any vague sense. It is the vapour pressure of a saturated salt solution, read as a fraction of pure water's: wherever the air holds more water than that, the solution is the state with the lower chemical potential. Bring the humidity back down, and the puddle does not turn back into salt until the air is far drier — about 45 per cent.

The boiling point in a column of water, and the depth heating reaches it first. The temperature at which water boils against depth in a geyser's conduit 2240 m above sea level, where the air's pressure is 77.6 kPa (solid): 92.7 °C at the surface, 103.6 °C at 4 m, 115.5 °C at 10 m and 131.8 °C at 22 m — rising steeply near the top and more gently below, because each extra metre adds the same pressure to water whose boiling point responds less as it gets hotter. Dashed: the water as it heats from beneath, cooled to 80 °C at the top by the pool and below that warming by only 0.8 °C per metre, mixed by convection; 6 °C and 3 °C short of the curve at 4 m, and then touching it. It touches at 4 m, not at the bottom: there the water is still 13.7 °C below boiling.

The column that boils from the middle

Water at the bottom of a geyser's tube is hotter than water anywhere else in it, and it is not the water that boils first. The weight of the column raises the boiling point with depth faster than the heat from below raises the temperature, so the column reaches its boiling curve partway down, under a cooler surface layer. A little boiling there lifts the water above. Everything beneath is then at a lower pressure than the one it was heated under, so it is superheated and flashes, and the eruption feeds itself: past about fourteen centimetres of water thrown out, each slug of steam makes more room than it was given.

How a sphere of gas rings. The acoustic pressure against distance from the centre of a spherical cavity, for its three lowest radial modes, in which the gas breathes in and out in spherical shells: sin(zr/a)/(zr/a), with z = 4.4934, 7.7253, 10.9041 the roots of tan z = z that make the gas's motion vanish at the wall. Filled with argon at the triple point of water, 273.16 K, a sphere 62 mm in radius rings in these modes at 3551, 6104, 8616 Hz. Radial modes are the ones precision work uses: they have no motion along the wall, so the wall's roughness and the gas's viscosity there barely touch them, and their frequencies are set by the volume rather than the shape — a slightly imperfect sphere shifts them only at second order.

The note that fixed the kelvin

A sphere of argon, struck, rings at a pitch set by the speed of sound in it, and the speed of sound in argon contains only one quantity nobody could compute: the constant that converts a temperature into an energy. Argon's atoms can hold energy in exactly three ways, so its heat capacity is a count rather than a measurement, and its speed of sound is (5/3)kT/m with nothing left over. Ring the sphere at the temperature of melting ice, measure its size with radio waves in the same cavity, and Boltzmann's constant comes out to a part in a million. In 2019 that number stopped being measured and became the definition.

The extra work in letting the gas expand further. Pressure against volume, the volume in units of the combustion chamber and the pressure on a logarithmic axis, for two ideal engines that compress the charge by the same ratio of 10 and add the same heat: an Otto cycle that expands by 10 too and throws away gas still at 3.8 bar (dashed), and an over-expanded cycle that lets the burned gas expand 26.8 times, until it is back at the 1 bar it came in at (solid). The shaded region between the expansion line and the intake pressure, beyond the Otto engine's stroke, is where the extra work is done (on this logarithmic axis its area is not to scale): 136 kJ per kilogram of charge, raising the efficiency from 55.3 to 63.6 per cent. It is taken from gas that an Otto engine releases through its exhaust valve still pushing.

The engine that expands further than it compresses

A petrol engine's efficiency is usually quoted as a function of its compression ratio, and its compression ratio is held near ten because a higher one makes the fuel detonate. But the efficiency does not come from the compression. It comes from the expansion — from how far the burned gas is allowed to push the piston before it is thrown away — and an ordinary engine expands its gas exactly as far as it compressed it, simply because the same cylinder does both. Let the gas expand further than it was compressed and most of what a higher compression ratio would have given comes back with no knock at all. The price is paid in power, and it is why the most efficient petrol engines on the road are built this way.

The wait after an eruption, against the eruption. For 400 simulated eruptions of the model geyser: the interval to the next eruption against the duration of the eruption just finished (dots), and the straight line through them, interval = 33.8 min + 10.5 × duration, with a correlation of 0.94. The eruptions fall in two clusters, short ones of about 2.2 minutes followed by waits near an hour and long ones of about 5.0 minutes followed by waits near an hour and a half, and the line runs through both. Its slope is the ratio of the discharge rate to the refill rate, 10.7: each minute of eruption throws out water that takes 10.7 minutes to replace. Plotted against the eruption that follows instead, the same intervals show a correlation of 0.05 — nothing. The interval is set by the past and tells nothing about the future.

The geyser that remembers one eruption

Old Faithful is famous for being predictable, and its predictability has an odd shape. Its eruptions come in two kinds, short and long, and the wait after each is short or long to match — about an hour after a short eruption, about an hour and a half after a long one. So the next eruption can be forecast from the last, to within about ten minutes, and from nothing else. The reason is a bookkeeping of water. An eruption throws out a volume of water that the geyser must take back from below before it can erupt again, and a long eruption throws out more. The geyser carries a memory of exactly one eruption: how much it has to refill.

An attraction found by counting. How often two large discs are found at each gap between their surfaces, relative to how often they are found far apart, in a direct simulation of two discs among 270 small particles that pass through each other freely and do not enter the discs (dots), and the prediction from counting the small particles' available area (line); small particles a quarter of the large radius. Nothing in the simulation attracts anything: every move that does not put a small particle inside a large one is accepted. Yet the discs are found touching 2.7 times as often as apart (2.6 predicted), because touching returns the overlap of their exclusion zones to the small particles, and the probability of a separation is proportional to the number of arrangements the small particles have at it. Beyond a gap of twice the small radius the zones no longer overlap and the effect vanishes. Written as a potential, −kT times the logarithm of that ratio, it is an attraction about 1.0 kT deep.

The attraction that is the small ones' room

Put large particles into a crowd of small ones that cannot get inside them, and the large ones begin to stick together — though nothing in the mixture attracts anything. Each large particle keeps the small ones' centres out of a thin shell around it. When two large particles touch, their shells overlap, and the overlap is handed back to the small particles as room to move in. More room is more arrangements, more arrangements is more entropy, and the large particles are pushed together by the pressure of the small ones on their outer sides. Fumio Oosawa and Sho Asakura worked it out in 1954. The force is now used to purify DNA, sort colloids onto patterned surfaces and explain why the inside of a cell behaves as though its proteins were stickier than they are.

Why dissolving beats evaporating near absolute zero. Cooling power against temperature on logarithmic axes, for pumping on liquid helium-4, pumping on liquid helium-3 — both with a pump of 100 litres a second at room temperature — and a dilution refrigerator circulating 100 micromoles of helium-3 a second (the T² law, dashed above 0.1 K where it no longer holds). Evaporative cooling is proportional to the vapour pressure, which collapses as e^(−L/T): helium-4 gives 7·10⁻¹¹ W at 0.3 K, helium-3 10⁻⁴ W at 0.3 K and 2·10⁻²⁴ W at 50 mK. The dilution refrigerator's "vapour" is the dilute phase, whose density of helium-3 stays at 6.6% however cold it gets, so its cooling falls only as T²: 2·10⁻⁵ W at 50 mK. Below about 0.3 K nothing that evaporates can compete, and most experiments that must stay continuously below a tenth of a kelvin — superconducting quantum computers, bolometers for the cosmic microwave background, searches for dark matter — run on dilution.

The mixture that stays mixed at absolute zero

The third law says a mixture's entropy must vanish at absolute zero, so every mixture must either separate completely or order into a compound as it is cooled. Liquid helium-3 dissolved in helium-4 does neither. Below 0.87 kelvin it splits into two liquids, but the one that is mostly helium-4 keeps 6.6 per cent helium-3 however cold it gets — and the third law permits it, because the dissolved atoms are a degenerate Fermi gas with a single ground state and no entropy of mixing to lose. Pushing helium-3 across that boundary absorbs heat the way evaporation does, from a 'vapour' that never thins out. That is the dilution refrigerator, the machine that holds every superconducting quantum computer at ten millikelvin.

Two metals mixing at unequal speeds. The fraction of the fast-diffusing species A against position across a welded couple, A on the left and B on the right, at three times in the ratio 1 : 4 : 16, when A diffuses 3 times as fast as B, from Darken's equations, in units in which B's diffusivity and the middle time are both one; each profile's marker — an inert wire laid at the original weld — as a dot on it. The profiles spread as the square root of time, and they are lopsided: A pours into B's side faster than B into A's, so the mixed zone extends further into B. The markers do not stay at the weld. They move into A's side — by 0.39, 0.80 and 1.63 units at the three times — because the lattice planes of A's side are being emptied faster than they are filled. That drift, first seen in copper and brass in 1947, is the Kirkendall effect.

The markers that move because atoms are unequal

In 1947 Ernest Kirkendall and Alice Smigelskas wound fine molybdenum wires around a bar of brass, plated it thickly with copper, and heated it for fifty-six days. The wires, which do not dissolve in either metal, ended closer together: the brass inside them had shrunk. Zinc had diffused out of the brass faster than copper diffused in, and the lattice planes of the brass had been taken away along with the zinc. If atoms in a metal moved by swapping places with each other, the two flows would have had to be equal. They were not, and the only way the books balance is if atoms move by jumping into vacancies — empty lattice sites — that flow the other way. The wires were the first direct evidence that diffusion in metals runs on holes.

Newcomen's dilemma. Two consequences of how cold a Newcomen engine's injection spray leaves its cylinder, against that temperature: the pressure difference the atmosphere then pushes the piston down with, one atmosphere less the vapour pressure of the water left in the cylinder; and the steam the next stroke wastes reheating the cylinder's iron skin, as a fraction of the most it wastes at 30 °C (8.2 cylinder-fulls). A cylinder cooled to 40 °C gives a push of 0.93 atmospheres and wastes 7.0 cylinder-fulls of steam a stroke; one left at 70 °C wastes half as much but pushes with only 0.69. Every choice was bad, because the one vessel had to be both hot, to receive the steam, and cold, to make the vacuum. Newcomen's engineers ran it in between.

The cylinder that must not be cooled

The first steam engines that did useful work, Thomas Newcomen's pumps of 1712, made a vacuum by spraying cold water into a cylinder full of steam. It worked, and on every stroke it cooled the cylinder, so that the next charge of steam condensed on the cold iron until the walls were hot again — several cylinders' worth of steam spent heating metal that the next spray would cool. James Watt saw in 1765 that the vessel was being asked to be hot and cold at once, and moved the cold to a separate condenser. The improvement was a factor of three or four in coal. A century later the reason had a name: heat passing from hot to cold across a large difference, doing nothing on the way, is work thrown away.

The pipe that out-conducts any solid. The temperature difference needed to carry 30 W along a rod 300 mm long and 6 mm across, on a logarithmic scale: for solid aluminium 1343 K, for copper 796 K, for diamond, the best conductor of heat at room temperature, 159 K; and for a water heat pipe of the same size 0.8 K, almost all of it across the thin wet wicks where the heat goes in and comes out. Solid metal could not do it at all: the hot end of the copper rod would have to be hundreds of degrees above the cold. The heat pipe carries the heat as latent heat in moving vapour, and the vapour needs almost no temperature difference to flow. As an effective conductivity that is 1000 times copper's — which is why every laptop and many satellites move their heat in pipes rather than through metal.

The pipe that carries heat as vapour

A copper rod six millimetres across and thirty centimetres long, asked to carry thirty watts from a laptop's processor to its fan, would need its hot end eight hundred degrees above its cold end. The flattened copper tube that actually does the job manages it with less than a degree. It is a heat pipe: sealed, emptied of air, and holding a few drops of water and their vapour. The water boils at the hot end, the vapour rushes to the cold end and condenses, and a wick of sintered copper powder drags the liquid back by capillarity. The heat goes as latent heat, and vapour needs almost no temperature difference to flow — which makes the pipe a thousand times better than copper, as long as its wick can keep up.

A cooling planet, read from the top. Temperature against depth in a body cooling by conduction from a uniform start at 3870 °C, with a surface held at nought, at four ages: the exact error-function solution for a half-space. Cooling reaches downward as the square root of time, so after a hundred million years the upper few hundred kilometres have cooled and the interior knows nothing about it. The gradient at the surface is what a mine measures, and it falls as the inverse root of time. A thermal gradient therefore gives an age, and this is what Kelvin saw in 1862. He assumed the Earth was a solid conductor from the start, on a good basis, as everyone believed then, and these curves are the mathematical consequence of that belief.

The age a cooling Earth gives

In 1862 William Thomson, later Lord Kelvin, took a thermometer reading from the bottom of a mine and a piece of mathematics from Fourier and computed the age of the Earth: about a hundred million years. Geologists, needing far longer, had no physics with which to reply. The arithmetic was right. The usual story blames radioactivity for overturning it, but half the heat is not radioactive, and the error lay in assuming the inside of the planet stays still.

All essays