Astrophysics

The bar a wave of gravity barely touches

The first gravitational-wave detectors were tonnes of aluminium, hung in a vacuum and listened to for ringing. A passing wave stretches and squeezes the bar, and the bar can answer only by vibrating against its own stiffness, so the fraction of the wave it catches is set by the speed of sound over the speed of light, squared: about three parts in ten billion. That one ratio explains why the bars had to be enormous and cold, why Joseph Weber's reported detections of 1969 would have required the Galaxy to be converting dozens of suns a day into gravitational waves, and why the waves were finally found by light instead.

Assumes: The wave that stretches one way and squeezes the other · What the instrument actually hears

The wave that stretches one way and squeezes the other mentioned Joseph Weber’s detectors in passing: aluminium cylinders, with piezoelectric crystals glued round their waists, designed to ring when a gravitational wave at the right frequency went by, and replaced in the end by interferometers. The cylinders deserve more than a mention, because they are the simplest possible gravitational-wave detector and the physics of how little they catch is clean and surprising. A bar is a lump of matter with one mode of vibration that a wave can excite. Working out how strongly the wave excites it turns into a single formula, and the formula has the speed of sound in it.

A passing wave changes the distances between free particles by a fraction hh, the strain, which for the waves detected since 2015 is about 10−2110^{-21}. A bar is not made of free particles. Its atoms are held at fixed spacings by the stiffness of the metal, so the wave’s stretch has to work against that stiffness, and the bar responds by vibrating rather than by simply following the stretch. That is both the opportunity and the problem. The opportunity is resonance: a wave at the bar’s own frequency adds to the vibration cycle after cycle. The problem is that the coupling between the wave and the vibration is weak in a way that no engineering can fix.

Where a wave of gravity pulls on a bar, and how the bar can answer. Along a cylinder 3 m long lying across a passing wave: the wave's tidal force per unit length, proportional to distance from the centre (dashed), and the shape of the bar's fundamental mode of vibration, sin(πx/L) (solid), in which the two ends move in opposite directions while the middle stays still. Both are scaled to 1 at the ends. Only the part of the force that matches the mode's shape can set it ringing; the overlap makes the bar behave like two masses, each half of the bar's, 0.405 of its length apart (dots). For an aluminium bar of 2,296 kg ringing at 900 Hz, as Louisiana's ALLEGRO did, the mode has a sound speed of 5400 m/s, and the tide it answers to is a strain of the order of 10⁻²¹: a stretch of the whole bar by 3.0 × 10⁻²¹ m.
Fig. 1 Along a cylinder 3 m long lying across a passing wave: the wave’s tidal force per unit length, proportional to distance from the centre (dashed), and the bar’s fundamental mode shape sin⁡(πx/L)\sin(\pi x/L) (solid), ends moving oppositely, middle still; both scaled to 1 at the ends. Only the part of the force matching the mode sets it ringing, and the bar behaves like two masses, each half of its own, 0.405 of its length apart (dots). An aluminium bar of 2,296 kg at 900 Hz, like ALLEGRO, has a sound speed of 5,400 m/s.

The force a wave exerts on a bar

In the frame of the bar’s centre, a gravitational wave acts as a tidal force. A piece of the bar a distance xx from the centre feels a force proportional to xx and to the wave’s second time derivative, 12ρh¨x\tfrac12\rho\ddot h x per unit volume, pulling the two ends apart and then pushing them together as the wave oscillates. The first figure draws that force along the bar as a straight line through the middle.

The bar’s lowest mode of longitudinal vibration has the ends moving in opposite directions while the middle stays still, with a displacement along the bar shaped like sin⁡(πx/L)\sin(\pi x/L). It has the right symmetry to be driven: the force pulls the ends apart exactly when the mode wants them apart. Projecting the force onto the mode shape gives the effective drive, and the integral comes out simply. The bar’s fundamental mode responds as if it were two point masses, each with half the bar’s mass, held 4L/π24L/\pi^2 apart by a spring tuned to the bar’s frequency, about four-tenths of its length. The third of itself a spring carries met the same reduction for a spring with mass: a distributed object vibrating in one mode behaves like a smaller lumped one, with an effective mass and an effective length set by the shape of the mode.

The frequency of that mode is the time for sound to cross the bar and back, f0=v/2Lf_0 = v/2L. A three-metre aluminium bar, with sound at 5.4 kilometres a second, rings near 900 hertz. That is where the bars listened, because it is where the sources thought most likely in the 1960s and 1970s — the final instants of collapsing stellar cores and merging neutron stars — were expected to put most of their power.

A burst, and the ringing it leaves

A burst sets the bar ringing long after it has passed. A burst of gravitational wave at the bar's own frequency, 900 Hz, lasting a few milliseconds with peak strain 10⁻¹⁸ (thin, scaled to its peak), and the bar's response, integrated from its equation of motion (thick, scaled to its largest swing, 6.1·10⁻¹⁸ m). The bar keeps every cycle the wave puts in, adding them up while the burst lasts and then ringing at constant amplitude, because at a quality factor of 10⁶ its ringing takes 354 s to fall by a factor of e. The swing it is left with is of order the strain times the bar's length, a few thousandths of the diameter of a proton for this strain; a real burst near 10⁻²¹ leaves a thousand times less, and the ringing has to be read against the bar's own thermal jiggling.
Fig. 2 A burst at the bar’s own frequency, 900 Hz, lasting a few milliseconds with peak strain 10−1810^{-18} (thin, scaled to its peak), and the bar’s response, integrated from its equation of motion (thick, scaled to its largest swing, 6.1×10−186.1\times10^{-18} m). The bar adds up every cycle the wave puts in, then rings at constant amplitude; at a quality factor of 10610^6 the ringing takes 354 s to fall by a factor of e.

The second figure drives the mode with a short burst of wave at its own frequency and integrates the result. While the burst lasts, the bar’s swing grows, each cycle of the wave adding to it in step; when the burst has gone, the swing stays. At a quality factor of a million, typical of a well-made aluminium bar suspended in a vacuum, the ringing would take six minutes to die away by a factor of e, so the bar holds the burst’s energy long enough for an amplifier to read it.

The swing itself is of order the strain times the bar’s length. For the burst drawn, with a strain a thousand times larger than anything detected since, it is 6×10−186\times10^{-18} metres, a few thousandths of the diameter of a proton. A burst of the size actually detected would leave a swing a thousand times smaller still. The whole art of the bars was in reading motions of that size against everything else that makes the bar move.

How much of a wave a bar can catch

The energy the burst leaves in the bar can be compared with the energy the wave carried. A gravitational wave carries energy at a rate set by its strain and frequency, c3/16πGc^3/16\pi G times the mean square of h˙\dot h, per square metre. That coefficient is enormous — about 103410^{34} in SI units — which is the sense in which spacetime is stiff: a strain of 10−2110^{-21} at a kilohertz carries about a tenth of a watt per square metre, and a burst of the strengths Weber was looking for would have carried megajoules. The bar’s integrated cross-section, the ratio of the energy it absorbs to the energy per unit area per unit frequency that crossed it, works out to

Σ0=8π GMc(vc)2.\Sigma_0 = \frac{8}{\pi}\,\frac{GM}{c}\left(\frac{v}{c}\right)^2 .

The figures obtain it by integrating the driven mode and dividing by the fluence, and check the two against each other.

How much of a wave of gravity a bar can catch. A resonant bar's cross-section for gravitational waves, integrated over its resonance, (8/π)(GM/c)(v/c)², against its mass, on logarithmic axes, for sound speeds of 3.0 km/s, 5.4 km/s, 10.0 km/s (aluminium is about 5.4). It is proportional to the mass, which is why the bars weighed tonnes, and to the square of the speed of sound over the speed of light — the bar answers a strain only by being stretched against its own stiffness, and (v/c)² is about 3 × 10⁻¹⁰. Weber's first bar caught 2.3·10⁻²⁵ m² Hz: to absorb a single joule, a burst must carry 4.4·10²⁴ J per square metre per hertz at the bar's frequency. Dots: Weber's bar, and the two cryogenic aluminium bars that ran longest, ALLEGRO in Louisiana and NAUTILUS near Rome.
Fig. 3 A resonant bar’s integrated cross-section, (8/π)(GM/c)(v/c)2(8/\pi)(GM/c)(v/c)^2, against its mass on logarithmic axes, for sound speeds of 3.0, 5.4 and 10 km/s. It grows with the mass and with the square of the speed of sound over the speed of light, about 3×10−103\times10^{-10} for aluminium. Weber’s first bar caught 2.3×10−25 m2 Hz2.3\times10^{-25}\ \mathrm{m^2\,Hz}: to absorb a single joule, a burst must carry 4.4×10244.4\times10^{24} J per square metre per hertz at the bar’s frequency. Dots: Weber’s bar, ALLEGRO and NAUTILUS.

The formula is the argument of this essay in one line. The cross-section grows with the mass, which is why the bars weighed tonnes. It contains Newton’s constant, because the coupling of a wave of gravity to matter is gravitational. And it contains the square of v/cv/c, the speed of sound in the bar over the speed of light. That last factor is about three parts in ten billion for aluminium, and it is the reason a bar catches so little. A bar responds to a strain only by deforming against its own elasticity, and the elastic forces that resist deformation are, compared with anything relativistic, very soft. A body that could transmit stress at nearly the speed of light would couple ten billion times better; no such material exists, because sound speed is bounded by how stiff matter can be, and the stiffest solids are diamond and sapphire, two or three times aluminium.

The third figure draws the cross-section against mass. Weber’s first bar, 1.4 tonnes, caught 2.3×10−252.3\times10^{-25} square metres for each hertz of bandwidth the wave covered — compared with the 6.7×10−296.7\times10^{-29} square metres with which a single electron scatters light, which the cross-section that forgets the colour found to be the smallest cross-section a charge can offer, the bar is not small. It is the energy per unit area that must be enormous. To absorb a single joule, a wave would have to bring more than 102410^{24} joules per square metre per hertz to the bar.

Why a sharper resonance does not catch more

It would be natural to think that a bar with a higher quality factor catches more of a wave, since it rings longer and more sharply. It does not, and the reason is the relation the width that is a lifetime established for every resonance. A resonance’s response, plotted against frequency, has a height proportional to its quality factor and a width inversely proportional to it, and the cross-section integrated across the resonance — the quantity in the formula — is the area under that curve, which does not depend on the quality factor at all. A bar with ten times the QQ absorbs from a narrower slice of the wave’s spectrum, ten times more strongly; for a burst whose spectrum is broad and flat across the resonance, the two effects cancel exactly.

What the quality factor does change is how long the bar remembers. The energy a burst leaves stays for Q/πf0Q/\pi f_0 seconds, and the thermal motion it has to be read against changes on the same slow timescale. That is the lever the cryogenic bars pulled: a high QQ does not enlarge the signal, but it slows the noise down, so a measurement short compared with the ringing time compares the bar’s state just before and just after a burst and sees only the small part of the thermal energy that could have changed in between.

A planet as a bar

The formula says nothing about what the bar is made of except its mass and its sound speed, and the largest elastic bodies available are planets. Freeman Dyson worked out in 1969 how the solid Earth responds to a passing gravitational wave, and the Earth’s own quadrupole oscillation, the slowest way it can ring, is excited by a wave of the right frequency exactly as a bar’s mode is. The tone the whole Earth rings at met that oscillation as the free vibration set going by great earthquakes, with a period of about 54 minutes. Its mass is 6×10246\times10^{24} kilograms, twenty-one orders of magnitude more than a bar’s, and the seismic wave speed is comparable with aluminium’s, so the Earth’s cross-section is vast.

It is also deafened by its own noise, which is why no one has heard a gravitational wave in it. Earthquakes, the atmosphere and the oceans keep the Earth’s free oscillations permanently excited far above any plausible gravitational signal at millihertz frequencies. The Moon is quieter, and seismometers left by the Apollo astronauts were used, much later, to set upper limits on a background of gravitational waves; proposals to place more sensitive instruments on the lunar surface, turning the whole Moon into a resonant detector for the band between the interferometers on the ground and in space, are under study.

Standing out from the bar’s own motion

The burst it takes to stand out from the bar's own jiggling. The peak strain of a one-millisecond burst at 900 Hz that leaves, in a 2.3-tonne aluminium bar, as much energy as its fundamental mode already holds from heat or from quantum mechanics, on a logarithmic axis: room temperature, 300 K, 1.6·10⁻¹⁶; liquid helium, 4.2 K, 1.8·10⁻¹⁷; cooled to 0.1 K, 2.8·10⁻¹⁸; effective noise, 1 mK, 2.8·10⁻¹⁹; one quantum, ħω₀, 1.9·10⁻²¹. A bar at room temperature, as Weber's were, needs a strain near 10⁻¹⁶ to rise above one kT of thermal energy. Cooling helps only as the square root of the temperature, and the best cryogenic bars, with amplifiers that kept their effective noise near a millikelvin, reached a few times 10⁻¹⁹. The bursts that interferometers detected in 2015 and since have strains near 10⁻²¹, below every bar's reach; and they sweep through a bar's narrow resonance in a fraction of the time it would need to build up.
Fig. 4 The peak strain of a one-millisecond burst at 900 Hz that leaves, in a 2.3-tonne aluminium bar, as much energy as its fundamental mode already holds: 1.6×10−161.6\times10^{-16} for kT at room temperature, 1.8×10−171.8\times10^{-17} at 4.2 K, 2.8×10−182.8\times10^{-18} at 0.1 K, 2.8×10−192.8\times10^{-19} for an effective noise of 1 mK, and 1.9×10−211.9\times10^{-21} for one quantum ℏω0\hbar\omega_0. The dashed line is the strain of the bursts interferometers detected from 2015.

A bar at temperature TT is never still. Half a kT for every way of moving gives each mode of vibration an average thermal energy of kTkT, half kinetic and half potential, and the fundamental mode is no exception. A burst is detectable only if it changes the mode’s energy by something comparable. The fourth figure turns that into a strain: the burst that leaves one kTkT in the bar.

At room temperature, which is where Weber’s bars ran, it takes a strain of 1.6×10−161.6\times10^{-16} — five orders of magnitude larger than the waves that have since been detected. Cooling helps less than it seems to, because the energy a burst deposits goes as the square of its strain, so the strain needed falls only as the square root of the temperature. Liquid helium gains a factor of eight; a dilution refrigerator at a tenth of a kelvin, another six. The cryogenic bars of the 1990s and 2000s — ALLEGRO at Louisiana State University, EXPLORER at CERN and NAUTILUS at Frascati, AURIGA near Padua — did better than their temperatures alone allow by a trick the high quality factor makes possible. The noise a high Q moves out of the way found that a resonator with a high quality factor concentrates its thermal motion into a slow random walk of amplitude and phase, so a measurement short compared with the ringing time sees only a small part of the thermal energy. With superconducting amplifiers the bars reached effective noise temperatures of a few millikelvin and strains of a few times 10−1910^{-19}.

Below that is the floor no cooling removes. The mode is a quantum oscillator, and its energy comes in units of ℏω0\hbar\omega_0; a burst that deposits less than one quantum cannot be told from no burst by a measurement of energy. For a 900-hertz bar that quantum corresponds to a strain of about 2×10−212\times10^{-21}. The noise every amplifier must add found the related limit on any amplifier that reads both quadratures of an oscillation: it must add half a quantum of noise of its own. The bars never reached it. Proposals to evade it by measuring only one quadrature — the first ideas of quantum non-demolition measurement, worked out by Vladimir Braginsky and Kip Thorne largely for the bars — were developed in the end for the interferometers, which use squeezed light to do it.

What Weber’s events would have cost

What Weber's events would have cost their source. If a burst deposited one kT of room-temperature energy in Weber's 1.4-tonne bar at 1,661 Hz — the size of the coincident events he reported from 1969 — it carried 1.8·10⁴ J per square metre per hertz at the resonance, about 1.8·10⁷ J per square metre over a millisecond burst's kilohertz of bandwidth. Emitted equally in every direction, that is the rest energy of the number of suns drawn here against the source's distance: 82 solar masses per event from the centre of the Galaxy, 8.2 kpc away. At the roughly one event a day he reported, the Galaxy would have been converting 3·10⁴ solar masses a year into gravitational waves — enough to dissolve it in far less than its age. Beaming the waves into the Galaxy's plane, or tuning the source's spectrum to the bars' frequency, lowered the bill by large factors and never far enough; the events were not reproduced by other bars, and none had a source.
Fig. 5 If a burst deposited one room-temperature kT in Weber’s 1.4-tonne bar at 1,661 Hz, it carried 1.8×1041.8\times10^4 J per square metre per hertz at the resonance, about 1.8×1071.8\times10^7 J per square metre over a millisecond burst’s kilohertz of bandwidth. Emitted equally in every direction, that is the rest energy of the number of suns drawn against the source’s distance: 82 per event from the Galactic centre, 8.2 kpc away.

In 1969 Weber reported coincident excitations of bars at Maryland and at Argonne, a thousand kilometres apart, at a rate of about one a day, and later claimed that they came preferentially from the direction of the Galactic centre. The fifth figure puts his claim through the cross-section. An event large enough to stand out in a room-temperature bar had to deposit roughly kTkT, and to deposit that, the wave had to bring some eighteen thousand joules per square metre per hertz to the bar, about twenty million joules per square metre over the bandwidth of a millisecond burst. If its source was at the centre of the Galaxy and radiated equally in every direction, it emitted the rest energy of about eighty suns in each event — at one a day, tens of thousands of solar masses a year.

The Galaxy cannot afford that. Its whole stellar mass is about a hundred billion suns, and at that rate it would have radiated itself away in a few million years. The argument was made at the time, by Dennis Sciama and others, in versions that allowed the waves to be beamed into the plane of the Galaxy or concentrated near the bars’ frequency, which reduced the bill by large factors and never far enough to be credible. Other groups built bars and found nothing, and by the mid-1970s the community had concluded that the events were not gravitational waves. Weber never accepted the conclusion. He is nonetheless remembered as the founder of the field, because the instruments, the coincidence technique and the very idea that the waves could be looked for were his.

Why light won

The interferometers that detected the first waves in 2015 escape the bar’s limitation in a specific way. A bar measures a strain through the elastic deformation of a solid, and pays (v/c)2(v/c)^2 for doing so. An interferometer’s mirrors are free masses, hung on pendulums whose own resonance is far below the frequencies being measured, so they move with the wave as free particles would; and the distance between them is measured with light, which travels at cc. There is no factor of (v/c)2(v/c)^2 to pay, and the arms can be made kilometres long, so the change in distance to be read is 10−1810^{-18} metres rather than 10−2110^{-21}. What the instrument actually hears took apart how that reading is made and where it has blind spots of its own.

The interferometers are also broadband. A bar responds strongly only within a small fraction of a hertz of its resonance, extended by clever readout to tens of hertz and, in the best of the bars, about a hundred, and the waves from merging black holes and neutron stars are chirps that sweep through hundreds of hertz in a fraction of a second. The orbit that has to shrink found why: the frequency of the waves rises as the orbit decays, faster and faster. A chirp spends only milliseconds crossing a bar’s band, less time than the bar needs to ring up, and the bar sees a fraction of it.

Where the cross-section formula holds

The formula is for a thin cylinder lying broadside to the wave with the wave’s polarisation along its axis — the most favourable orientation. Any other orientation reduces it, by the fourth power of the sine of the angle between the bar and the direction the wave comes from, and by the cosine of twice the polarisation angle. It treats the bar as a single mode; higher modes couple too, at their own frequencies, but more weakly. It assumes the wave’s wavelength, about 300 kilometres at a kilohertz, is far longer than the bar, which it is by five orders of magnitude.

Spheres were proposed as the successors to bars because a sphere has five quadrupole modes at nearly the same frequency, one for each independent pattern of tidal stretch, so it hears waves from every direction and can tell their polarisation. Two small spherical detectors, MiniGRAIL in Leiden and Mario Schenberg in São Paulo, were built and cooled in the 2000s. They run into the same factor of (v/c)2(v/c)^2 as bars, and none approached the interferometers’ sensitivity.

A ringing drawn without its noise

The figures draw the mode’s response to a single burst in a noiseless bar; the ringing a real bar holds is that response added to a thermal motion whose amplitude, at room temperature, is a hundred times larger, and the detection is a statistical decision made on the output of an amplifier, not a glance at a trace. Nor do they show how the bar’s vibration is converted into a signal: Weber’s piezoelectric crystals, the resonant transducers and superconducting quantum interference devices of the cryogenic bars, each with its own noise and its own coupling losses. And the energy budget of the fifth figure assumes a source radiating equally in all directions with a spectrum about a kilohertz wide; every serious attempt to rescue Weber’s events altered those assumptions, and the figure’s numbers would move with them.

Still open: whether a resonant detector can still find something an interferometer misses

Bars and spheres are not entirely retired as ideas. A resonant detector listens at a single frequency for a long time, and that suits sources that emit continuously at a known frequency, or high-frequency waves above the interferometers’ bands, where several speculative sources — from the early universe, or from hypothetical light dark-matter particles clustered round black holes — have been proposed. Small quartz and sapphire resonators cooled to millikelvin temperatures are being run at megahertz frequencies for exactly that reason, and two unexplained events seen in one in 2021 are still being argued over. Whether any gravitational wave at megahertz frequencies exists, and whether a resonant detector of any size can reach the strains needed to see it, is not known.

The habit worth carrying away is to find the soft link in a chain of coupling. A bar catches a gravitational wave with a cross-section (8/π)(GM/c)(v/c)2(8/\pi)(GM/c)(v/c)^2, about 2×10−252\times10^{-25} square metre-hertz for Weber’s 1.4 tonnes, because it can answer a strain only by deforming against stiffness that transmits stress at the speed of sound; so a burst large enough to stand out in a room-temperature bar would have cost a source at the Galactic centre about eighty solar masses, and the waves were found instead by free mirrors read with light.

Part 9 of 9

This essay is one argument about Gravitational 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.

The objects named here

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

Cross-sectionGravitational wavesNormal modesQuality factorResonanceStrainThermal noiseTidal force