Chirp Mass and the Chirp Signal

Two compact bodies bound in a slowly decaying orbit are the loudest sustained gravitational-wave sources in the Universe, and their signal has a signature so distinctive it named the phenomenon: a chirp. As the orbit radiates its energy away it shrinks, the bodies speed up, and — because the dominant gravitational-wave frequency is exactly twice the orbital frequency — the emitted tone climbs in pitch while swelling in amplitude, right up to merger. The astonishing content of this article is that the entire shape of that frequency sweep, at leading (Newtonian-quadrupole) post-Newtonian order, is fixed by a single combination of the two masses: the chirp mass M=(m1m2)3/5/(m1+m2)1/5\mathcal{M}=(m_1m_2)^{3/5}/(m_1+m_2)^{1/5}. The result is not a coincidence of algebra but a consequence of energy balance. The orbital binding energy E(f)E(f) lost by the system must equal the luminosity L(f)\mathcal{L}(f) carried off by the waves. Both depend on the component masses only through M\mathcal{M} at this order, so the balance law dE/dt=−L\mathrm{d}E/\mathrm{d}t=-\mathcal{L} collapses to one ordinary differential equation for f(t)f(t) with M\mathcal{M} as its lone parameter. Integrating it gives the celebrated power law f∝(tc−t)−3/8f\propto(t_c-t)^{-3/8} that diverges at coalescence. This one parameter is also the best-measured. The wave's phase, accumulated over hundreds or thousands of cycles, is what matched filtering pins down most sharply, and it depends on M\mathcal{M} alone at leading order — which is why M\mathcal{M} is routinely measured to a fraction of a percent while the individual masses stay fuzzy. The amplitude carries the complement: once M\mathcal{M} and ff are known, the observed strain reads off the luminosity distance directly, making every inspiral a self-calibrating standard siren. This article derives the chirp equation, its integral, the strain, and the phase — the machinery behind modern gravitational-wave astronomy.

Key takeaways
  • The chirp mass M=(m1m2)3/5/(m1+m2)1/5=Mη3/5\mathcal{M}=(m_1m_2)^{3/5}/(m_1+m_2)^{1/5}=M\eta^{3/5} is the single mass combination that fixes, at leading post-Newtonian order, both the rate at which an inspiral's frequency rises and its strain amplitude.
  • Energy balance between the orbital binding energy and the quadrupole luminosity gives the chirp equation f˙=965π8/3(GM/c3)5/3f11/3\dot f=\tfrac{96}{5}\pi^{8/3}(G\mathcal{M}/c^3)^{5/3}f^{11/3}, whose solution sweeps the frequency upward and formally diverges at coalescence.
  • Integrating the chirp gives f∝(tc−t)−3/8f\propto(t_c-t)^{-3/8} and a time-to-merger τ=5256(GM/c3)−5/3(πf)−8/3\tau=\tfrac{5}{256}(G\mathcal{M}/c^3)^{-5/3}(\pi f)^{-8/3} — the countdown clock of an inspiral.
  • The strain scales as h∝(GM/c2)5/3f2/3/Dh\propto(G\mathcal{M}/c^2)^{5/3}f^{2/3}/D, so measuring f(t)f(t) and f˙\dot f yields M\mathcal{M}, and combining amplitude with M\mathcal{M} yields the luminosity distance DD — the standard-siren idea.
  • Because the phase depends on the masses only through M\mathcal{M}, it is by far the best-measured parameter; the mass ratio and spins enter only at higher post-Newtonian order.
  • For GW150914 the detector-frame chirp mass was about 30 M⊙30\,M_\odot, sweeping the LIGO band from roughly 3535 Hz to 250250 Hz in about 0.20.2 s — the audible chirp that gave the class of signals its name.
  • Key idea: at leading post-Newtonian order, energy balance between orbital binding energy and quadrupole luminosity makes the frequency sweep f˙\dot f and the strain amplitude hh depend on the two masses only through the single combination M=(m1m2)3/5/M1/5\mathcal{M}=(m_1m_2)^{3/5}/M^{1/5}.
The core insight

At leading post-Newtonian order, energy balance between orbital binding energy and quadrupole luminosity makes the frequency sweep f˙\dot f and the strain amplitude hh depend on the two component masses only through the single combination M=(m1m2)3/5/M1/5\mathcal{M}=(m_1m_2)^{3/5}/M^{1/5} — which is therefore the dominant observable of every inspiral.

Where this lives today

The chirp mass is the working currency of gravitational-wave astronomy. In every LIGO–Virgo–KAGRA detection the low-latency pipelines estimate M\mathcal{M} within seconds of a trigger, because a rough value sorts the event into a class — around 1 M⊙1\,M_\odot flags a binary neutron star, tens of M⊙M_\odot a binary black hole — and predicts how long remains before merger, setting whether telescopes can be turned to catch an electromagnetic counterpart. The spectacular payoff came with GW170817: the chirp mass fixed the intrinsic loudness, the strain fixed the luminosity distance as a standard siren, and an optical counterpart supplied the host-galaxy redshift, yielding a Hubble constant independent of the traditional distance ladder — a template for resolving the current tension between local and cosmological values of H0H_0. Looking ahead, the space-based detector LISA will track the chirp of supermassive black-hole binaries for months, measuring M\mathcal{M} so precisely they become distance markers across the observable Universe. In every case the first number out of the data, and the most trustworthy, is M\mathcal{M}.

The physical picture

Picture two masses spiralling inward. Gravity bleeds orbital energy away as radiation; as the orbit tightens the pair speeds up, so the pitch rises. The emitted waveform is literally a chirp — a tone climbing in frequency and swelling in loudness up to merger, the gravitational analogue of a bird call. Two facts make it quantitative: the wave frequency is twice the orbital frequency (f=2forbf=2f_{\mathrm{orb}}, because the mass quadrupole returns to itself twice per orbit), and the strain grows as f2/3f^{2/3} as the frequency climbs. The first figure is that waveform.

Figure
The gravitational-wave chirp h(t) of a compact binary. As the orbit decays the frequency climbs as (t_c − t)^−3/8 and the strain swells as f^2/3, so a rising, loudening tone runs up to merger and then rings down. The gravitational-wave frequency is twice the orbital frequency, f = 2 f_orb
The chirp: rising frequency, growing strain — f sweeps up as (t_c − t)^−3/8, h grows as f^2/3 time merger t_c inspiral — low pitch, faint chirp — pitch and strain climb ringdown

The magic of the chirp mass is that the rate of pitch-rise encodes the masses in one specific way. Heavy binaries chirp fast and loud, light ones slowly and softly — but only through M\mathcal{M}. Two systems with the same M\mathcal{M} but different mass ratios produce nearly identical early waveforms; the ratio only reveals itself as the merger nears. Think of M\mathcal{M} as the tempo knob and volume knob fused into one: both how fast the signal accelerates and how strong it is climb together as M5/3\mathcal{M}^{5/3}, so measuring the sweep alone pins it down with extraordinary precision. The second figure is what an observatory actually reconstructs — the frequency track of a real event, whose slope is the measurement of M\mathcal{M}.

Figure
Reconstructed gravitational-wave frequency track of a GW150914-like binary with detector-frame chirp mass about 30 solar masses. The (t_c − t)^−3/8 law sweeps the frequency from 35 Hz to 250 Hz in the final 0.16 s, and the steepness of this sweep alone fixes the chirp mass ℳ
The punchline: the frequency sweep 35 → 250 Hz in 0.16 s — its slope measures ℳ ≈ 30 M⊙ 0 100 200 300 0 50 100 150 time before merger (ms) GW frequency f (Hz) 35 Hz — enters band 100 Hz 250 Hz — merger f ∝ (t_c − t)^−3/8 : nearly flat, then the final plunge

The core equations

Everything begins with the definition. Of the mass combinations that could appear, it is the chirp mass M\mathcal{M}, not the total mass MM, that controls the radiation. It is a weighted geometric mean of the components, biased toward the lighter one, and can equivalently be written through the symmetric mass ratio η\eta.

M=(m1m2)3/5(m1+m2)1/5=M η3/5(1)\mathcal{M}=\frac{(m_1 m_2)^{3/5}}{(m_1+m_2)^{1/5}}=M\,\eta^{3/5} \qquad (1)

where: M\mathcal{M} — chirp mass (source frame); m1,m2m_1,m_2 — component masses; M=m1+m2M=m_1+m_2 — total mass; η=m1m2/M2\eta=m_1m_2/M^2 — symmetric mass ratio, bounded 0<η≤140<\eta\le\tfrac14. Reading it: for fixed total mass, M\mathcal{M} is largest when the masses are equal (η=14\eta=\tfrac14) and shrinks as the ratio grows extreme. The dynamical heart of the subject is the chirp equation — the rate at which the gravitational-wave frequency rises. It is the balance law dE/dt=−L\mathrm{d}E/\mathrm{d}t=-\mathcal{L} rewritten as an equation for f˙\dot f, in which M\mathcal{M} is the only surviving mass combination.

f˙=965 π8/3(GMc3)5/3f11/3(2)\dot f=\frac{96}{5}\,\pi^{8/3}\left(\frac{G\mathcal{M}}{c^3}\right)^{5/3} f^{11/3} \qquad (2)

where: ff — gravitational-wave frequency, equal to 2forb2f_{\mathrm{orb}}; f˙=df/dt\dot f=\mathrm{d}f/\mathrm{d}t — the sweep rate; G,cG,c — Newton's constant and the speed of light. Reading it: the sweep accelerates steeply (as f11/3f^{11/3}), and its overall scale is set by M5/3\mathcal{M}^{5/3} — a heavier chirp mass means a faster chirp. Integrating (2) gives the explicit frequency-versus-time law. It is a pure power law in the time remaining before coalescence, the fingerprint of the leading-order inspiral.

f(t)=1π(5256 1tc−t)3/8(GMc3)−5/8(3)f(t)=\frac{1}{\pi}\left(\frac{5}{256}\,\frac{1}{t_c-t}\right)^{3/8}\left(\frac{G\mathcal{M}}{c^3}\right)^{-5/8} \qquad (3)

where: tct_c — the coalescence (merger) time; tc−tt_c-t — time remaining before merger. Reading it: ff diverges as (tc−t)−3/8(t_c-t)^{-3/8}, so the sweep is gentle for most of the inspiral and violent in the last cycles. Inverting the same integral gives the countdown clock: how long a binary at frequency ff has left before it merges.

τ=tc−t=5256(GMc3)−5/3(πf)−8/3(4)\tau=t_c-t=\frac{5}{256}\left(\frac{G\mathcal{M}}{c^3}\right)^{-5/3}(\pi f)^{-8/3} \qquad (4)

where: τ\tau — time to coalescence from frequency ff. Reading it: light chirp masses linger (a neutron-star binary lasts tens of seconds), heavy ones flash through in a fraction of a second. The amplitude carries the second, independent piece of information. For a face-on (optimally oriented) circular binary the strain of a circular binary depends on M\mathcal{M}, the instantaneous frequency, and the luminosity distance DD.

h(t)=4D(GMc2)5/3(πf(t)c)2/3cos⁡Φ(t)(5)h(t)=\frac{4}{D}\left(\frac{G\mathcal{M}}{c^2}\right)^{5/3}\left(\frac{\pi f(t)}{c}\right)^{2/3}\cos\Phi(t) \qquad (5)

where: hh — strain for a face-on (optimally oriented) binary; DD — luminosity distance; Φ(t)\Phi(t) — accumulated gravitational-wave phase. Reading it: once M\mathcal{M} and ff are known, the amplitude is set by DD together with the binary's orientation — a generic binary radiates h+=h0 1+cos⁡2ι2cos⁡Φh_+ = h_0\,\tfrac{1+\cos^2\iota}{2}\cos\Phi and h×=h0cos⁡ι sin⁡Φh_\times = h_0\cos\iota\,\sin\Phi with h0h_0 the coefficient above, and the detector adds its antenna responses, so distance and inclination are measured together and are partially degenerate. The phase is the time integral of the frequency, and it is the quantity matched filtering measures most precisely. In the time domain it follows directly from integrating (3).

Φ(t)=Φc−2(c3(tc−t)5GM)5/8(6)\Phi(t)=\Phi_c-2\left(\frac{c^3(t_c-t)}{5G\mathcal{M}}\right)^{5/8} \qquad (6)

where: Φc\Phi_c — phase at coalescence. Reading it: the phase accumulates over hundreds to thousands of cycles and depends on M\mathcal{M} alone at this order — the origin of the superb M\mathcal{M} measurement. In matched filtering the equivalent frequency-domain phase carries a dominant term ∝M−5/3\propto\mathcal{M}^{-5/3}, so a tiny error in M\mathcal{M} misaligns the template by many radians and destroys the match. Finally, a gravitational wave redshifts exactly like light, so time and frequency at the detector are dilated by (1+z)(1+z). The directly measured quantity is therefore the redshifted, detector-frame chirp mass.

Mz=(1+z) M(7)\mathcal{M}_z=(1+z)\,\mathcal{M} \qquad (7)

where: zz — cosmological redshift of the source; Mz\mathcal{M}_z — detector-frame chirp mass. Reading it: the waveform alone yields Mz\mathcal{M}_z; recovering the intrinsic M\mathcal{M} needs an independent redshift.

Derivation AThe chirp equation from energy balance

The goal is to show that a quasi-circular Newtonian binary losing energy to quadrupole radiation obeys f˙∝M5/3f11/3\dot f\propto\mathcal{M}^{5/3}f^{11/3}, and that no other combination of the two masses survives. We assemble three ingredients the earlier articles own — the orbital binding energy, Kepler's third law, and the quadrupole luminosity — express each through the orbital frequency, and impose energy balance dE/dt=−L\mathrm{d}E/\mathrm{d}t=-\mathcal{L}. The single combination μM2/3=M5/3\mu M^{2/3}=\mathcal{M}^{5/3} then falls out and everything else cancels.

Step 1. For a circular orbit of separation rr the binding energy is E=−GμM/(2r)E=-G\mu M/(2r) with reduced mass μ=m1m2/M\mu=m_1m_2/M, and Kepler's third law gives ωorb2=GM/r3\omega_{\mathrm{orb}}^2=GM/r^3. Eliminating rr leaves E=−12μ(GM)2/3ωorb2/3E=-\tfrac12\mu(GM)^{2/3}\omega_{\mathrm{orb}}^{2/3}. Why: this trades the unobservable separation rr for the orbital frequency, which is directly tied to the wave frequency.

Step 2. Because the mass quadrupole repeats twice per orbit, the wave frequency is f=2forb=ωorb/πf=2f_{\mathrm{orb}}=\omega_{\mathrm{orb}}/\pi, so ωorb=πf\omega_{\mathrm{orb}}=\pi f and E=−12μ(GMπf)2/3E=-\tfrac12\mu(GM\pi f)^{2/3}. Why: this rewrites the energy in the one variable the detector records, the gravitational-wave frequency.

Step 3. The quadrupole formula applied to a circular binary gives the luminosity L=325(G/c5)μ2r4ωorb6\mathcal{L}=\tfrac{32}{5}(G/c^5)\mu^2 r^4\omega_{\mathrm{orb}}^6; substituting r=(GM/ωorb2)1/3r=(GM/\omega_{\mathrm{orb}}^2)^{1/3} and ωorb=πf\omega_{\mathrm{orb}}=\pi f turns it into L=325(G/c5)μ2(GM)4/3(πf)10/3\mathcal{L}=\tfrac{32}{5}(G/c^5)\mu^2(GM)^{4/3}(\pi f)^{10/3}. Why: this expresses the energy loss rate in the same variable as the energy, so the two can be balanced directly.

Step 4. Impose dE/dt=−L\mathrm{d}E/\mathrm{d}t=-\mathcal{L}, and write the left side as (dE/df)f˙(\mathrm{d}E/\mathrm{d}f)\dot f with dE/df=−13μ(GMπ)2/3f−1/3\mathrm{d}E/\mathrm{d}f=-\tfrac13\mu(GM\pi)^{2/3}f^{-1/3}. Solving for f˙\dot f cancels one power of μ\mu and combines the mass factors into μM2/3\mu M^{2/3}. Why: the balance law is the physical statement that the orbit shrinks exactly fast enough to supply the radiated power.

Step 5. Recognize that μM2/3=(m1m2/M) M2/3=m1m2 M−1/3\mu M^{2/3}=(m_1m_2/M)\,M^{2/3}=m_1m_2\,M^{-1/3} is exactly M5/3\mathcal{M}^{5/3}, giving f˙=965π8/3(GM/c3)5/3f11/3\dot f=\tfrac{96}{5}\pi^{8/3}(G\mathcal{M}/c^3)^{5/3}f^{11/3}. Why: this is the punchline — the two masses enter only through M=μ3/5M2/5=(m1m2)3/5/M1/5\mathcal{M}=\mu^{3/5}M^{2/5}=(m_1m_2)^{3/5}/M^{1/5}, so the frequency sweep measures M\mathcal{M} and nothing else at this order.

Go deeper · Derivation B · 1 minIntegrating the chirp — frequency sweep, coalescence time, and phase
Derivation BIntegrating the chirp — frequency sweep, coalescence time, and phase

The goal is to turn the differential chirp law into the three quantities an observation actually uses: the explicit frequency track f(t)f(t), the time-to-merger τ(f)\tau(f), and the accumulated phase Φ(t)\Phi(t) whose precision underlies the whole measurement. The chirp equation is separable, so we separate variables, integrate to the merger, invert for the frequency track, and integrate once more for the phase.

Step 1. Separate variables in equation (2): f−11/3 df=965π8/3(GM/c3)5/3 dtf^{-11/3}\,\mathrm{d}f=\tfrac{96}{5}\pi^{8/3}(G\mathcal{M}/c^3)^{5/3}\,\mathrm{d}t, with both sides now elementary. Why: the chirp equation couples ff and tt only through a power law, so the variables cleanly separate.

Step 2. Integrate from a time tt (frequency ff) up to coalescence tct_c (where f→∞f\to\infty). The left side gives 38f−8/3\tfrac38 f^{-8/3}, yielding 38f−8/3=965π8/3(GM/c3)5/3(tc−t)\tfrac38 f^{-8/3}=\tfrac{96}{5}\pi^{8/3}(G\mathcal{M}/c^3)^{5/3}(t_c-t), which rearranges to τ=5256(GM/c3)−5/3(πf)−8/3\tau=\tfrac{5}{256}(G\mathcal{M}/c^3)^{-5/3}(\pi f)^{-8/3}. Why: this is the time-to-merger — the clock that tells an observatory how long a binary at a given frequency has left.

Step 3. Solve the same relation for ff: f(t)=1π(52561tc−t)3/8(GM/c3)−5/8f(t)=\tfrac1\pi\big(\tfrac{5}{256}\tfrac{1}{t_c-t}\big)^{3/8}(G\mathcal{M}/c^3)^{-5/8}, the characteristic (tc−t)−3/8(t_c-t)^{-3/8} sweep. Why: this explicit track is what the second figure plots; its steepness is a direct read-out of M\mathcal{M}.

Step 4. Integrate the phase, Φ=2π∫f dt\Phi=2\pi\int f\,\mathrm{d}t, using ∫(tc−t)−3/8 dt=−85(tc−t)5/8\int(t_c-t)^{-3/8}\,\mathrm{d}t=-\tfrac85(t_c-t)^{5/8} and the algebraic identity 85(5256)3/8=(15)5/8\tfrac85\big(\tfrac{5}{256}\big)^{3/8}=\big(\tfrac15\big)^{5/8}, giving Φ(t)=Φc−2(c3(tc−t)/5GM)5/8\Phi(t)=\Phi_c-2\big(c^3(t_c-t)/5G\mathcal{M}\big)^{5/8}. Why: the phase, accrued over thousands of cycles and depending on M\mathcal{M} alone, is the quantity matched filtering measures best — hence M\mathcal{M}'s sub-percent precision. *The coefficient identity 85(5/256)3/8=(1/5)5/8\tfrac85(5/256)^{3/8}=(1/5)^{5/8} used in Step 4 can be confirmed by direct symbolic evaluation.*

Go deeper · Derivation C · 1 minThe strain amplitude and the standard siren
Derivation CThe strain amplitude and the standard siren

The goal is to obtain the strain amplitude of a circular binary and to show that it makes the source self-calibrating for distance. Starting from the far-field quadrupole strain, we insert the circular-orbit quadrupole moment and eliminate the separation with Kepler's law; the mass dependence again collapses to M5/3\mathcal{M}^{5/3}, after which the amplitude contains only the known M\mathcal{M}, the measured ff, and the distance DD — so the strain reads off DD.

Step 1. The far-field transverse-traceless strain is hijTT=(2G/c4D)Q¨ijTTh^{TT}_{ij}=(2G/c^4D)\ddot Q^{TT}_{ij}, and for a circular binary the reduced quadrupole oscillates at 2ωorb2\omega_{\mathrm{orb}} with magnitude Q¨∼μr2ωorb2\ddot Q\sim\mu r^2\omega_{\mathrm{orb}}^2. Why: the factor of 22 in the oscillation frequency is precisely why the wave frequency is twice the orbital frequency, and the μr2\mu r^2 scaling carries the mass and size of the orbit.

Step 2. Eliminate rr with Kepler's law r=(GM/ωorb2)1/3r=(GM/\omega_{\mathrm{orb}}^2)^{1/3} and set ωorb=πf\omega_{\mathrm{orb}}=\pi f; the mass factors combine as μM2/3=M5/3\mu M^{2/3}=\mathcal{M}^{5/3}, giving h∝(GM/c2)5/3(πf/c)2/3/Dh\propto(G\mathcal{M}/c^2)^{5/3}(\pi f/c)^{2/3}/D. Why: the same collapse to M\mathcal{M} that governed the frequency sweep governs the amplitude, so one parameter sets both tempo and volume.

Step 3. Read off the physical content: with M\mathcal{M} fixed by the frequency evolution and ff measured directly, the amplitude in equation (5) is set by the distance DD and the binary's inclination. Inverting it gives the luminosity distance with no external calibrator, up to the distance-inclination degeneracy that a single detector cannot break — the dominant systematic in the GW170817 standard-siren measurement of H0H_0. Why: this is the standard-siren principle — an inspiral advertises its own intrinsic loudness, so its apparent loudness measures distance, bypassing the cosmic distance ladder entirely. *The amplitude scaling h∝M5/3f2/3/Dh\propto\mathcal{M}^{5/3}f^{2/3}/D derived here is the basis of the standard-siren inversion used in the GW170817 Hubble-constant measurement.*

Limiting cases

The chirp mass and its consequences take instructive forms in the corners of parameter space. In the table below, MM denotes the total mass throughout.

CaseResultMeaning
Equal massesM=M/26/5≈0.435 M\mathcal{M}=M/2^{6/5}\approx0.435\,MWith m1=m2m_1=m_2 so η=14\eta=\tfrac14, the chirp mass is 0.4350.435 of the total — the maximum of M/M\mathcal{M}/M at fixed total mass.
Extreme mass ratioM≈m2(m2/m1)2/5≪m2\mathcal{M}\approx m_2(m_2/m_1)^{2/5}\ll m_2For m1≫m2m_1\gg m_2 the chirp mass is far below even the small mass; the inspiral is slow and the chirp faint.
Near coalescencef∝(tc−t)−3/8f\propto(t_c-t)^{-3/8}The frequency formally diverges as the leading post-Newtonian treatment breaks down; the real system merges and rings down first.
Cosmological sourceMz=(1+z)M\mathcal{M}_z=(1+z)\mathcal{M}The detector-frame chirp mass is redshifted; only Mz\mathcal{M}_z is directly observable from the waveform.
Amplitude scalingh∝M5/3f2/3/Dh\propto\mathcal{M}^{5/3}f^{2/3}/DStrain and sweep rate both grow as M5/3\mathcal{M}^{5/3}, tying tempo and loudness to one parameter.

Worked examples

1 · Chirp mass of GW150914 from its components (confidence builder)

The first detected black-hole merger had component masses m1≈36 M⊙m_1\approx36\,M_\odot and m2≈29 M⊙m_2\approx29\,M_\odot and lay at redshift z≈0.09z\approx0.09. Find its source-frame chirp mass, its symmetric mass ratio, and the detector-frame chirp mass an observatory would report.

Step 1. Form the products: m1m2=36×29=1044 M⊙2m_1m_2=36\times29=1044\,M_\odot^2 and M=65 M⊙M=65\,M_\odot. Then M=(1044)3/5/(65)1/5\mathcal{M}=(1044)^{3/5}/(65)^{1/5}; evaluating, (1044)3/5=64.75(1044)^{3/5}=64.75 and (65)1/5=2.305(65)^{1/5}=2.305, so M=64.75/2.305=28.1 M⊙\mathcal{M}=64.75/2.305=\mathbf{28.1\,M_\odot}.

Step 2. The symmetric mass ratio is η=1044/652=0.247\eta=1044/65^2=\mathbf{0.247}, just below the equal-mass ceiling 14\tfrac14. Consistently, the equal-mass estimate 0.435 M=28.3 M⊙0.435\,M=28.3\,M_\odot sits barely above the exact 28.1 M⊙28.1\,M_\odot, confirming a nearly equal-mass binary.

Step 3. Redshift to the detector frame: Mz=(1+z)M=1.09×28.1=30.6 M⊙\mathcal{M}_z=(1+z)\mathcal{M}=1.09\times28.1=\mathbf{30.6\,M_\odot}. This is the quantity the waveform directly delivers, matching the reported value of about 30 M⊙30\,M_\odot. Note how sharply M\mathcal{M} is fixed while the individual masses each carry several-solar-mass uncertainty — the phase constrains the combination, not the parts. M=28.1 M⊙\mathcal{M}=28.1\,M_\odot source-frame, η=0.247\eta=0.247, Mz=30.6 M⊙\mathcal{M}_z=30.6\,M_\odot detector-frame.

2 · The countdown clock of a neutron-star binary

A binary neutron star like GW170817 has source-frame chirp mass M≈1.188 M⊙\mathcal{M}\approx1.188\,M_\odot. When its gravitational-wave frequency reaches 3030 Hz, how long remains until merger, and roughly how many wave cycles does it accumulate from there to coalescence?

Step 1. Convert the chirp mass to a time using GM⊙/c3=4.926×10−6 sG M_\odot/c^3=4.926\times10^{-6}\,\mathrm{s}, so GM/c3=1.188×4.926×10−6=5.85×10−6 sG\mathcal{M}/c^3=1.188\times4.926\times10^{-6}=5.85\times10^{-6}\,\mathrm{s}. This one number, together with the frequency, sets every timescale.

Step 2. Apply the time-to-coalescence formula (4): τ=5256(GM/c3)−5/3(πf)−8/3\tau=\tfrac{5}{256}(G\mathcal{M}/c^3)^{-5/3}(\pi f)^{-8/3}. With πf=π×30=94.2 s−1\pi f=\pi\times30=94.2\,\mathrm{s^{-1}}, evaluate (5.85×10−6)−5/3=5.26×108(5.85\times10^{-6})^{-5/3}=5.26\times10^{8} and (94.2)−8/3=5.43×10−6(94.2)^{-8/3}=5.43\times10^{-6}, giving τ=0.0195×5.26×108×5.43×10−6=56 s\tau=0.0195\times5.26\times10^8\times5.43\times10^{-6}=\mathbf{56\ s}. A neutron-star binary lingers for the better part of a minute in the audio band — ample time for a low-latency alert.

Step 3. Estimate the cycle count from the phase (6): Φ=2(c3τ/5GM)5/8\Phi=2\big(c^3\tau/5G\mathcal{M}\big)^{5/8} with τ=56 s\tau=56\,\mathrm{s} and 5GM/c3=2.93×10−5 s5G\mathcal{M}/c^3=2.93\times10^{-5}\,\mathrm{s}, so c3τ/5GM=1.9×106c^3\tau/5G\mathcal{M}=1.9\times10^{6} and Φ=2×(1.9×106)5/8=1.7×104\Phi=2\times(1.9\times10^6)^{5/8}=1.7\times10^{4} rad. Dividing by 2π2\pi gives N≈2700 cyclesN\approx\mathbf{2700\ cycles}. That enormous cycle count, each cycle sensitive to M\mathcal{M}, is exactly why a neutron-star binary's chirp mass can be pinned to four significant figures. τ≈56\tau\approx56 s and ≈2700\approx2700 wave cycles from 3030 Hz to merger.

Physics of the portrait

The chirp mass is where a century of theory meets a decade of observation: a quantity that emerged from Einstein's quadrupole formula, was made concrete by Peters and Mathews in 1963, was silently confirmed by a dying pulsar, and finally became the loudest number in a new astronomy.

Origins · 1 minchirp mass and how this was discovered
CM
chirp mass
the combination ℳ = (m₁m₂)^{3/5}/M^{1/5} that governs an inspiral · post-Newtonian gravitational-wave astronomy

The chirp mass is one of the youngest quantitative concepts in gravitation and one of the most consequential. Its lineage begins with Albert Einstein, who in 1918 corrected his own earlier work and derived the quadrupole formula, showing that an accelerating mass distribution radiates gravitational energy at a rate set by the third time derivative of its mass quadrupole moment — the engine that drains an orbit. For decades the effect was thought unmeasurably small and, to some, physically suspect. The decisive calculation came in 1963 from Philip Peters and Jon Mathews, who applied the quadrupole formula to a Keplerian two-body system and computed the radiated power and orbital decay; their work contains, implicitly, the collapse of the two masses into the single combination now named the chirp mass, along with the coefficient 96/596/5 that governs the frequency sweep. Nature confirmed it before any detector could: in 1974 Russell Hulse and Joseph Taylor discovered the binary pulsar PSR B1913+16, whose orbit shrank in exact agreement with the Peters–Mathews prediction, earning the 1993 Nobel Prize and proving that gravitational radiation carries energy just as the quadrupole formula says. The final act was direct detection: when LIGO recorded GW150914 in 2015, the very first thing extracted from the data was the chirp mass, read straight off the frequency sweep. What Einstein wrote as an abstract radiation formula had become a number you could measure from a tone rising through an audio band — a black-hole pair weighed by the sound of its collision.

History
1918
— Einstein derives the quadrupole formula for gravitational radiation
1963
— Peters and Mathews compute the radiated power and orbital decay of a binary
1974
— Hulse and Taylor discover the binary pulsar PSR B1913+16
2015
— LIGO records GW150914, the first directly observed chirp
2017
— GW170817 delivers the first gravitational-wave standard-siren measurement of H₀

Going further

Texts. Misner, Thorne and Wheeler, Gravitation (1973), section 36 · Wald, General Relativity (1984), section 4.4 · Carroll, Spacetime and Geometry (2004), sections 7.5–7.6 · Hartle, Gravity (2003), chapter 23 · Schutz, A First Course in General Relativity, 2nd ed. (2009), sections 9.3–9.4 · Maggiore, Gravitational Waves, Vol. 1 (2008), section 4.1. The energy-balance derivation (orbital energy + Kepler + quadrupole luminosity) yields f˙∝M5/3f11/3\dot f\propto\mathcal{M}^{5/3}f^{11/3}; integrating gives τ(f)\tau(f) and Φ(t)\Phi(t). Key numerical checks: the equal-mass value M=0.435 M\mathcal{M}=0.435\,M and the coefficient identity 85(5/256)3/8=(1/5)5/8\tfrac85(5/256)^{3/8}=(1/5)^{5/8}. The GW150914 figures use component masses from Phys. Rev. Lett. 116, 061102 (2016) and the GW170817 standard-siren result from Nature 551, 85 (2017). The figures adopt GM⊙/c3=4.926×10−6G M_\odot/c^3=4.926\times10^{-6} s and a detector-frame M≈30 M⊙\mathcal{M}\approx30\,M_\odot.

Prerequisites. Quadrupole Formula · Plus and Cross Polarizations · Multipole Expansion of Gravitational Radiation · Gravitational Wave Equation. Leads to. Isaacson Stress-Energy of Gravitational Waves · Transverse-Traceless Gauge · Geodesic Deviation and Tidal Action of Waves.

Ask the tutor: Ask me to weigh a real detection — give me two component masses or a measured frequency and its sweep rate, and we will read off the chirp mass, the time to merger, and the standard-siren distance together.

Why is the gravitational-wave frequency twice the orbital frequency?

Because the source of radiation is the mass quadrupole moment, not the dipole. A binary's mass distribution looks identical after the two bodies swap places, which happens after only half an orbit, so the quadrupole moment returns to its original configuration twice per orbital period. Radiation is driven by the second time derivative of that moment, so the dominant wave oscillates at twice the orbital frequency, f=2forbf=2f_{\mathrm{orb}}. (There is no mass-dipole radiation at all, since momentum conservation makes the mass dipole's second derivative vanish.) This factor of two threads through every formula here, from the energy balance to the strain.

Why does the chirp mass and not the total mass govern the inspiral?

It falls out of energy balance, but there is a clean way to see it. The radiated power depends on the reduced mass μ\mu (how much mass is actually orbiting and shaking) while the orbital dynamics depend on the total mass MM (what holds the orbit together). Balancing the energy loss against the frequency evolution combines these two in exactly the proportion μM2/3\mu M^{2/3}, which is M5/3\mathcal{M}^{5/3}. So the chirp mass is the specific blend of "how much is radiating" and "how strong is the binding" that sets the sweep. The total mass and mass ratio do matter — but only at higher post-Newtonian order, as corrections to this leading behaviour.

If the chirp mass is so well measured, why are the individual masses uncertain?

Because the chirp mass is a single symmetric combination of two numbers, and one equation cannot solve for two unknowns. At leading order the waveform depends only on M\mathcal{M}, so the data constrain that combination beautifully but say nothing about how the mass is split between the bodies. The split enters only through higher post-Newtonian phase terms, which are smaller, accumulate over fewer effective cycles, and are partly degenerate with the bodies' spins. The upshot is a long, thin "banana" in the mass plane: M\mathcal{M} is nailed, while the individual masses drift along a curve of nearly constant M\mathcal{M}.

What does a standard siren measure, and why is it "standard"?

A standard siren is the gravitational-wave analogue of a standard candle. A Type Ia supernova has a known intrinsic brightness, so its apparent brightness gives its distance. An inspiral has a known intrinsic strain: once M\mathcal{M} is read from the frequency sweep and ff is measured, general relativity predicts exactly how loud the wave should be at the source. Comparing that to the observed strain yields the luminosity distance DD directly — no ladder of calibrators, no assumptions about stellar physics. It is "standard" because the calibration comes from the theory itself. Pairing DD with a redshift from an electromagnetic counterpart gives the Hubble constant in one clean step.

Does the frequency really diverge at coalescence?

Only in the leading-order formula, and that divergence is a signpost, not a physical infinity. The relation f∝(tc−t)−3/8f\propto(t_c-t)^{-3/8} comes from a Newtonian orbit corrected by the lowest-order radiation reaction; it assumes point masses on a shrinking circular orbit and knows nothing of the strong-field dynamics of merger. Long before ff could blow up, the bodies reach their innermost stable circular orbit, plunge, and merge — at which point the inspiral description hands off to numerical relativity and the signal becomes the ringdown of the final object. The formal divergence simply marks where the leading-order expansion runs out of validity.

Can you get the chirp mass from just the frequency and its rate of change?

Yes, and this is the cleanest way to see why M\mathcal{M} is the primary observable. The chirp equation f˙=965π8/3(GM/c3)5/3f11/3\dot f=\tfrac{96}{5}\pi^{8/3}(G\mathcal{M}/c^3)^{5/3}f^{11/3} can be solved algebraically for the chirp mass: measure the instantaneous frequency ff and its rate of increase f˙\dot f at any moment during the inspiral, and M\mathcal{M} follows immediately, with no reference to distance, mass ratio, or spins. In practice one does better by matched-filtering the whole phase evolution than by using a single (f,f˙)(f,\dot f) pair, but the principle is the same: the chirp mass is written directly into the tempo of the rising tone.

Before this

Where this leads