Waves

How many ways there are to vibrate

A drum has infinitely many modes, and below any given frequency it has a finite number of them. That number turns out to depend on the drum's area and the length of its rim and — to the accuracy anybody uses — on nothing else about its shape. Almost every result in thermal physics that involves waves is an application of that count.

Assumes: Only some notes fit, and that is where discreteness comes from · The drum that has no harmonics

Only some notes fit on a string, and which ones is a question about that particular string. This essay is about a different question — how many fit below a given frequency — which turns out to have an answer that barely depends on the object at all.

How many modes a square has below a given wavenumber. The number of vibration modes of a square with wavenumber below k, counted exactly — every eigenvalue of this region is a closed form, so the staircase is the true count and not an estimate. There are 265 of them below k = 60. The smooth curves are what Weyl's law predicts. The upper one is the leading term alone, the area times k² over 4π, and it is too high by 21.5 modes at the right-hand edge; the lower one subtracts the perimeter term, the perimeter times k over 4π, and is out by 2.4. The content of the law is that the count depends on the region through its area and its perimeter and — to this order — through nothing else at all: not through its shape, not through where its corners are, not through whether it is convex. The staircase's steps are the individual modes, and they cluster where two different pairs of indices give the same wavenumber. That the count is smooth in the large while being a staircase in the small is what makes a mode count usable in thermodynamics, where it appears as a density of states and never as a list.
Fig. 1 The exact number of vibration modes of a square membrane with wavenumber below k, drawn as a staircase, against what Weyl’s law predicts. There are 265 of them below k = 60. The upper smooth curve is the area term alone and is 21.5 modes too high; the lower one subtracts the perimeter term.

Counting instead of solving

The modes of a rectangle are known exactly: with the edges held, the wavenumber of the mode with mm half-waves across and nn along is k=πm2/a2+n2/b2k = \pi\sqrt{m^2/a^2 + n^2/b^2}. Every mode is a lattice point in the quarter-plane of positive integers, and asking how many modes lie below kk is asking how many lattice points lie inside a quarter-ellipse.

That question has an obvious answer and a subtlety. The obvious answer is the area of the quarter-ellipse, which works out to Ak2/4πAk^2/4\pi with AA the area of the membrane. The subtlety is the boundary: lattice points on the axes are excluded, because a mode with m=0m=0 is a membrane at rest, and the correction is proportional to the perimeter of the region.

A string is the one-dimensional case and the one where the counting is trivial: its modes are evenly spaced in frequency, so the number below any frequency is that frequency divided by the spacing, and there is nothing to estimate. Nothing about that generalises to two dimensions except the idea of counting itself — which turns out to be the part worth keeping.

Weyl proved in 1911 that the leading term is universal — that for any reasonable region in the plane the number of modes below kk is Ak2/4πAk^2/4\pi plus something smaller, whatever the shape — and the perimeter term was conjectured at the same time and proved much later. Written out:

N(k)=Ak24πLk4π+N(k) = \frac{Ak^2}{4\pi} - \frac{Lk}{4\pi} + \cdots

with LL the length of the boundary. The minus sign is the statement that a fixed edge costs modes: the boundary is where the wave must vanish, and the more boundary a region has for its area, the fewer ways there are to fit a given wavelength into it.

6 modes of a drum, and their frequency ratios. Nodal-line diagrams for 6 modes of a circular membrane, each labelled with its frequency as a multiple of the lowest mode's. A mode (m, n) has m nodal diameters and n − 1 nodal circles, and the circles are drawn at the radii where the computed radial function J_m(j(m,n)·r/R) crosses zero — not at guessed fractions of the radius. The two tints are the two directions the head is moving in at that instant, and the lines between them are the parts of it that never move. The ratios are 1.000, 1.593, 2.136, 2.295, 2.653, 2.917: each one is a quotient of two zeros of Bessel functions, computed here from the power series and checked against their published values to 4.4e-7. Not one is a whole number, which is why a drum has no harmonic series and no pitch in the sense a string has one — and why these same ratios belong to every circular membrane ever stretched, whatever it is made of and however tightly it is pulled.
Fig. 2 Six modes of a circular membrane. Each of them is a separate solution of a boundary-value problem and each has its own frequency; the counting law is a statement about the whole infinite list of them, and it does not require any member of the list to be known.

Three shapes, one area

The claim that the count depends on the shape only through two numbers is testable, and the test does not need any approximation: there are shapes whose spectra are exactly known.

A rectangle is one. A right isoceles triangle is another, because it is half of a square and its modes are the square’s modes that vanish on the diagonal — which are the ones with m>nm > n. So a square, an oblong and a triangle can all be given exactly the same area and their modes counted exactly, with no numerical eigenvalue solver anywhere.

Three shapes of the same area, counted. A square, a right triangle and an oblong four times as long as it is wide, all of exactly the same area, with the leading term of Weyl's law subtracted from each one's exact mode count. If the area were the whole story the three residuals would sit on top of one another. They do not: each falls along a straight line of its own, and the measured slopes are -0.323, -0.387, -0.384 against the −L/4π the law predicts, -0.318, -0.384, -0.398. The three shapes have perimeters 4.000, 4.828, 5.000 and that is exactly what separates them. So a spectrum reports two things about a region without being told either: its area, from how fast the count grows, and its perimeter, from how far the count falls behind. The square is the shape of least perimeter here and therefore the one with the most modes below any given k. What the spectrum does not report is the shape itself — two regions of the same area and perimeter can have different spectra, and famously two different shapes can have the same one.
Fig. 3 A square, a right triangle and a 4:1 oblong, all of area exactly one, with the area term of Weyl’s law subtracted from each one’s exact count. If the area were the whole story the three residuals would coincide. Instead each falls along a line of its own, with measured slopes matching the −L/4π the law predicts.

The three perimeters are 4, 4.83 and 5 for the same area, and those are exactly what separate the three lines. Subtract the perimeter term as well and the three collapse onto one another.

So a spectrum reports two things about a region without being asked: its area, from how fast the count grows, and its perimeter, from how far the count lags behind. Both are read off a list of frequencies by somebody who has never seen the shape. That is a striking amount of geometry to recover from a sequence of numbers — and it is also the whole of what can be recovered this way, which is the other half of the story.

Can one hear the shape of a drum? Kac asked it in that form in 1966. The answer is no: Gordon, Webb and Wolpert built in 1992 a pair of differently shaped regions with identical spectra. They have the same area and the same perimeter, necessarily, because the count says so — and they are not the same shape.

What the count is for

Counting modes would be a curiosity if it were not the step on which several of the largest results in thermal physics turn.

The mode count in one, two and three dimensions. The number of modes below a wavenumber for a string, a square membrane and a cubical box, all of unit size, on logarithmic axes. The three staircases straighten to slopes of 1.02, 2.04, 3.08, measured over the top fifth of the range and approaching one, two and three from above — the edge and surface corrections are deficits that shrink relative to the leading term, so a count is always still catching up. The count goes as the volume times k to the power of the dimension, and the constant in front is the volume of the unit ball divided by the appropriate power of 2π. That single fact is the origin of several results that look unrelated. The density of modes per unit frequency in three dimensions goes as ω², which is the factor in front of the Planck distribution and the reason the classical equipartition argument produced an ultraviolet catastrophe. The same ω² is what makes a solid's heat capacity rise as T³ at low temperature. And in one dimension the count is linear, which is why a string's modes are evenly spaced and a room's are not. Counting is the step that turns a wave equation into a thermodynamic quantity, and it is almost always done by this estimate rather than by listing anything.
Fig. 4 The mode count in one, two and three dimensions on logarithmic axes, with measured slopes of 1.02, 2.04 and 3.08 — approaching one, two and three from above, since the edge and surface corrections are deficits that shrink relative to the leading term. In three dimensions the count goes as k³, so the number of modes per unit frequency goes as ω².

That ω2\omega^2 is the factor in front of the Planck distribution. A cavity’s radiation is its electromagnetic modes; the number of them per unit frequency is fixed by the volume and by the dimension and by nothing else; and multiplying that count by the average energy of a mode gives the spectrum.

The count is what makes a blackbody spectrum the shape it is. The rise on the short-wavelength side is the mode count winning and the collapse beyond the peak is the Boltzmann factor winning, and the curve that would not come down is what happens when only the first of those is included. The catastrophe was never a failure of the counting; the counting was right, and what was missing was the factor that switches the high modes off.

The classical argument gave every mode kTkT, which is the equipartition result and was not the mistake. The mistake was arithmetic that had nothing wrong with it: an ω2\omega^2 count times a constant energy per mode integrates to infinity. The ultraviolet catastrophe is a statement about how many modes there are, and Planck’s constant enters as a bound on what the high-frequency ones can hold rather than on how many of them exist.

The same ω2\omega^2 count applied to sound waves in a crystal instead of light waves in a cavity gives Debye’s model of a solid’s heat capacity, with its T3T^3 at low temperature and a cut-off at the shortest wavelength the lattice allows. That cut-off is the one physical difference between the two problems: a cavity can hold arbitrarily short waves and a lattice cannot, because below the atomic spacing there is nothing left to wave.

Debye’s calculation of a solid’s heat capacity is the identical argument with sound in place of light and a highest mode in place of an unbounded spectrum. The T3T^3 law that comes out of it is one of the most accurately confirmed results in low-temperature physics, and every step of it after “count the modes” is bookkeeping.

How the universal part is proved

The reason the leading term cannot depend on the shape is worth seeing, because the argument is elementary and explains what “universal” is doing here.

Cut the region into small squares. Two comparison problems can then be built. In the first, every small square is given a fixed edge of its own, which is a more restrictive boundary condition than the region had — the wave is now forced to vanish on lines where it was previously free — and adding constraints can only raise frequencies, so this count is a lower bound. In the second, every internal edge is left free, which is less restrictive, and that count is an upper bound. Both are sums over squares of the count for a single square, which is known exactly.

As the squares are made smaller the two bounds close on each other, and what survives is the area divided by the area of one square, times the count for one square. Nothing about the region’s outline survives that limit except how much area it enclosed. The boundary contributes only through the squares it cuts, whose number is proportional to the perimeter — which is the second term, arriving as the error of the first.

That is Weyl’s argument in outline, and it is the reason the law’s leading term holds for regions with corners, with holes, and with boundaries too rough to have a tangent anywhere.

The frequency at which a room stops having modes

The staircase and the smooth curve describe the same object and are useful at opposite ends, and the frequency where the description changes over can be worked out.

Modes have widths, because a real room absorbs, and the width in frequency is set by how fast the sound decays. Below the frequency at which neighbouring modes start to overlap, a listener hears individual resonances and moving a metre changes the level of a tone by a great deal; above it, the modes are a continuum and the count is what matters. That crossover — the Schroeder frequency — is around 200 Hz in a domestic room and around 20 Hz in a concert hall, because the mode density grows as the volume.

It is why the same building acoustics that works well for a hall says nothing useful about a control room, why bass is the difficult part of every small listening space, and why the modes of a box are the right language at low frequency and a hopeless one at high. The same crossover appears in every system with a spectrum: individual states at the bottom, a density of states above.

The density of states

The count differentiated is the density of states — how many modes there are per unit frequency, or per unit energy — and in that form it appears in nearly every calculation involving a large number of quantum levels.

The same law turns up again under another name. The allowed states of a particle in a box are standing waves subject to the same boundary condition as a drum’s, so the counting law for them is Weyl’s law with different symbols — and the electron count in a metal is that law evaluated at the Fermi level.

A metal’s electrons fill states from the bottom up, and the energy at which they run out — the Fermi energy — is set by requiring the count of states below it to equal the number of electrons. That is Weyl’s law used as an equation to be solved rather than as a description.

Which energies exist at all is the mode count; how many of them are filled is the statistics. The two are separate questions, answered separately and then multiplied, and keeping them apart is what makes a metal’s electronic heat capacity tractable: the density of states comes from the geometry of the box and the occupancy from the temperature, and neither knows about the other.

The same quantity decides reaction rates, emission rates and conductivities, always in the same shape: a rate is a matrix element times the number of states available to go to. Fermi’s golden rule is that sentence written as an equation, and the second factor is a count.

It is also what makes the exclusion principle produce a pressure: the electrons must occupy distinct states, the number of states below an energy is fixed by the volume, and squeezing the volume raises the energy of the topmost occupied one. Every step of that argument except the exclusion is mode counting.

Where the smooth estimate stops describing anything is a periodic medium. A band gap is a range of frequencies with no modes at all, so the counting function is flat across it, and Weyl’s smooth estimate is an average through the structure rather than a description of it. The gap a repeat opens is where the missing modes went: they are not destroyed but pushed to the band edges, and the total count is conserved — which is why the smooth law is still right about the total while being wrong everywhere in detail.

The difference between two infinities

The most startling use of mode counting takes two counts that are each infinite and subtracts them, and gets a force that has been measured.

Every mode of the electromagnetic field has a zero-point energy of half a quantum, even with no photons in it. Summing that over all modes in any region gives infinity, and the infinity is usually discarded as an unmeasurable additive constant. What cannot be discarded is a difference: put two parallel conducting plates a distance apart and the modes between them are restricted — only wavelengths that fit are allowed — while the modes outside are not.

So the zero-point energy between the plates depends on their separation. Both the constrained and the free counts are infinite; their difference is finite, and differentiating it with respect to the separation gives a force. The plates attract, with a pressure

FA=π2c240d4,\frac{F}{A} = \frac{\pi^2\hbar c}{240\,d^4},

which contains no property of the plates at all — not their material, not their charge, not their mass. Only the geometry and two constants.

The fourth power makes it a short-range effect and a violent one. At a hundred nanometres the pressure is thirteen pascals, which is small and measurable; at ten nanometres it is over an atmosphere. That is why the force is irrelevant to anything at ordinary separations and is a dominant nuisance in micro-electromechanical devices, whose moving parts sit micrometres or less from their neighbours and occasionally snap together and stay.

Measuring it took until 1997 to do convincingly, and the difficulty is that everything else at that separation is also large — electrostatic patch potentials, residual gas, surface roughness. What makes the measurement possible is the exponent: nothing else in the vicinity falls as the fourth power of the separation, so sweeping the distance separates the Casimir term from its competitors by its shape rather than by its size.

The argument’s structure is worth carrying beyond its subject. Two divergent quantities whose difference is finite and observable is a pattern that appears throughout physics — in the self-energy of a charge, in the binding energy of a quantum field, in the running of a coupling — and the discipline in every case is the same: never compute either count, compute the difference, and check that the answer is independent of how the counting was cut off.

Changing the count changes the rate

The density of states enters every rate, and the useful consequence is that a rate can be changed by changing the states rather than by touching the thing that is doing the emitting.

An excited atom emits a photon at a rate proportional to how many modes are available for the photon to go into at the transition’s frequency. In free space that number is the ω2\omega^2 count of this essay, and the resulting lifetime is what atomic physics tabulates. Put the same atom inside a cavity and the count is a different number entirely.

Two things can happen. A cavity tuned to the transition concentrates the available modes into a narrow band around it, so the density of states at the atom’s frequency is much larger than in free space and the atom emits faster — by a factor that can be thousands. A cavity whose smallest dimension is below half the transition’s wavelength has no mode at that frequency at all, so the atom cannot emit, and its excited state’s lifetime is extended indefinitely.

Both have been done. The enhancement is used in light-emitting devices and in single-photon sources, where a fast, directional emission into one mode is exactly what is wanted; the suppression was demonstrated with atoms passed between closely spaced mirrors and with atoms in a waveguide below cutoff.

What makes it worth putting in this essay is the reassignment of blame. A spontaneous emission rate looks like a property of an atom, and it is not: it is a property of an atom and of the modes surrounding it, and the second factor is the mode count. An excited state has no lifetime of its own, only a lifetime given what it is allowed to decay into.

The same reasoning applies to a photonic crystal, whose band gap is a range of frequencies with no modes at all. An atom inside one, emitting into the gap, is in the strongest form of the suppressed case — and the prospect of building materials in which chosen transitions simply cannot happen is what the field was started for.

Where the smooth count is not enough

Weyl’s law is asymptotic, and there are two quite different ways of caring about what it misses.

The first is that the staircase is a staircase. Between one mode and the next the count does not move at all, and near the bottom of the spectrum the smooth curve is a poor description of anything: a room at 40 Hz has a handful of modes and the estimate is useless, which is why small-room acoustics is a different subject from concert-hall acoustics. The crossover is the frequency at which the modes overlap, and it is where the language changes from listing modes to counting them.

The partials of a drum, against the harmonics of a string. The lowest 6 modes of a circular membrane and the first 6 harmonics of a string, on one axis of frequency, both in units of their own lowest note. The string's are the whole numbers, because a one-dimensional boundary admits whole numbers of half-wavelengths. The drum's are 1.000, 1.593, 2.136, 2.295, 2.653, 2.917 — the zeros of Bessel functions divided by the first zero of J₀ — and they interleave with the harmonics without ever landing on one. The closest approach on this drawing is mode (1,2) at 2.917, which is 2.8% below the 3rd harmonic, and the loudest partial, (2,1) at 2.136, misses the octave by 6.8%. A quarter of a semitone is about 1%, so these are not small errors that an ear rounds off: they are intervals with no name. The scale is set separately, and only the scale depends on the drum: at R = 0.325 m, T = 1900 N/m and σ = 0.262 kg/m² the lowest mode is 100.3 Hz and the rest follow it at these ratios.
Fig. 5 The lowest partials of a circular membrane against a string’s harmonics. At the bottom of a spectrum the individual modes are what matter and no smooth count says anything; a drum sounds as it does because of the first six numbers here, not because of how many there are below some k.

The second is that the fluctuation of the true count about the smooth one is itself a subject. Its statistics depend on whether the corresponding classical motion is regular or chaotic — a rectangle and a stadium have quite different fluctuations with identical Weyl terms — and that observation is the foundation of quantum chaos. Order out of the random has a counterpart here: the smooth part of a spectrum knows the geometry, and the fluctuating part knows the dynamics.

The habit generalises past modes. Counting configurations rather than frequencies is the same manoeuvre thermal physics uses everywhere: a quantity hopeless to enumerate one at a time is easy to count in bulk, and the bulk count is what the physics actually needs. Weyl’s law is that habit applied to a wave equation.

The count that is not of modes

One more use is worth naming because it looks like a different subject and is the same arithmetic.

The entropy of a gas is the logarithm of how many microscopic states it can be in, and counting those states means counting cells in a phase space of positions and momenta — which is the same lattice count in six dimensions rather than two. The volume factor that appears in every expression for the entropy of an ideal gas is that count, and the arbitrary cell size that classical statistical mechanics had to introduce to make it finite is Planck’s constant, arriving from the same direction it arrived from in the blackbody problem.

So the pattern is consistent across the subject: a quantity that is impossible to enumerate is easy to count in bulk, the count is a volume divided by a cell, and the physics is what multiplies it.

What the pictures cannot show

Every count here is for a fixed edge. A free edge changes the sign of the perimeter term, and a mixed boundary gives something between. The area term is untouched, which is the sense in which it is universal, and the perimeter term is where the boundary condition shows up.

The shapes were chosen because their spectra are exact. A general region needs a numerical eigenvalue solver, and the agreement with Weyl’s law is then a test of two things at once. Using a square, an oblong and a triangle keeps the test on the law.

Nothing here is about degeneracy. The square has many modes sharing a frequency, and a small change of shape splits them without changing the count. That splitting is what a real drum’s timbre is sensitive to and the count is blind to.

And the third term is not a perimeter. After the area and the perimeter comes a constant that depends on the corners and on the curvature of the boundary — a disc and a square of the same area and perimeter differ in it — and beyond that the expansion stops being a series in any ordinary sense.

The ladder from here

Later rungs on this anchor: the density of states as the object rather than the count; the same law in three dimensions with the surface term, which is what a cavity’s spectrum needs; the isospectral drums and what their construction actually does; and the spectral fluctuations, which are where the counting argument hands over to quantum chaos.

The neighbouring ladders are the modes of a string, where counting is trivial and the idea is visible, the drum’s own frequencies, which are what is being counted, and the blackbody spectrum, which is this count multiplied by an energy.

Part 4 of 7

This essay is one argument about Standing waves. The others:

What links here

Essays that reach for this one mid-argument — the half of a link its own author cannot write down.

What this makes readable

Essays that declare this one a prerequisite.

The objects named here

The third axis, after the field and the reading path: the things themselves, and every essay that touches each one.

BlackbodyBoundary conditionDebye modelDensity of statesEigenvalueMode countingNormal modesStanding wave