Tidal Deformability and the Love Number: Reading a Neutron Star's Squishiness in the Chirp

What it is and why it matters

A neutron star packs 1–2 M☉ into a sphere only ~11–13 km across, reaching central densities several times the nuclear saturation density ρ₀ ≈ 2.7×10¹⁴ g cm⁻³. How stiff that matter is — how much pressure it musters against gravity — is encoded in the equation of state (EoS), one of the great open problems in nuclear astrophysics. Tidal deformability is the cleanest gravitational handle on it.

When a companion's gravity varies across a star's diameter, it raises a tide: the star develops an induced mass quadrupole Q_ij proportional to the external tidal field E_ij. The constant of proportionality is the tidal deformability λ, with Q_ij = −λ E_ij. A stiff, incompressible star barely responds (small λ); a soft, centrally-condensed star deforms readily (large λ). Because λ depends on the full pressure-density profile, measuring it pins down the radius and stiffness of matter at densities no collider or reactor can produce, making it a genuinely astrophysical laboratory for cold, dense QCD.

The mechanism, step by step

Start with Newtonian intuition, then add relativity. (1) The companion of mass M' at separation r produces a tidal field E_ij ~ M'/r³ that stretches the star along the line of centers. (2) The star's fluid redistributes into a prolate shape, generating an induced quadrupole moment Q_ij = −λ E_ij. (3) The dimensionless Love number normalizes this: λ = (2/3) k₂ R⁵/G, so k₂ measures how much of the maximal fluid response the star actually achieves.

(4) In full general relativity, k₂ is computed by solving the Tolman–Oppenheimer–Volkoff structure equations together with a linear perturbation equation for a metric function y(r), integrated from center to surface. The result depends only on the compactness C = GM/Rc² and the internal profile. (5) The energy cost of raising this tide comes out of the orbit, accelerating the inspiral. Because the deformation stores and radiates energy slightly faster than two point masses would, the binary sweeps through frequency a touch more quickly near merger — the physical fingerprint we detect.

Characteristic numbers, scales, and the key relation

The central relation ties everything together:

Λ = (2/3) k₂ (Rc²/GM)⁵ = (2/3) k₂ C⁻⁵

The steep C⁻⁵ dependence is why Λ is such a sensitive radius probe: a 1 km change in R shifts Λ by tens of percent. For a canonical 1.4 M☉ neutron star, k₂ ≈ 0.05–0.15 and C ≈ 0.14–0.18, giving Λ₁.₄ of order a few hundred. Black holes are the crucial contrast: the static Love number of a Schwarzschild or Kerr black hole is exactly zero, so Λ = 0 — a black hole cannot be tidally deformed, a striking and much-studied prediction of general relativity.

For a merging binary, the two stars combine into a single mass-weighted parameter Λ̃ (tilde-Lambda), which is what the data actually constrain. The name honors Augustus E. H. Love, who introduced h, k, and (via Shida) l to describe Earth's solid tides in his 1911 Adams-Prize work Some Problems of Geodynamics; the stellar apsidal-motion constant is the same k₂.

How it is detected

Tidal deformability leaves its mark on the phase of the gravitational-wave chirp, not the amplitude. In the post-Newtonian expansion of the waveform, matter effects for non-spinning bodies first appear at 5PN order — formally a small, high-order correction — but they grow rapidly as the orbit shrinks, becoming most important above ~400–500 Hz in the last few orbits before contact. Detectors (LIGO Hanford/Livingston, Virgo, KAGRA) track this accumulating phase against point-mass templates; the tiny extra phase advance, integrated over hundreds of cycles, is measurable.

The landmark case is GW170817 (17 August 2017), a binary neutron star merger at ~40 Mpc. Its inspiral constrained Λ̃ and, for a 1.4 M☉ star, Λ₁.₄ ≲ 800 in the discovery analysis, tightened to roughly 70–580 with EoS-informed priors — ruling out the stiffest equations of state and favoring neutron-star radii near 11–12 km. GW190425, a heavier binary, added a second, weaker constraint. Foundational theory: Flanagan & Hinderer (2008) and Hinderer et al., building on Damour's relativistic tidal work.

Where it operates and how it differs from related effects

Adiabatic tidal deformability governs the late inspiral of binaries containing at least one neutron star (BNS, and neutron-star–black-hole systems where the NS can be measured). It is distinct from several cousins. Dynamical tides arise when the orbital frequency approaches a stellar oscillation mode (notably the f-mode near ~1–2 kHz), resonantly enhancing the response beyond the static Love-number value. Tidal disruption — the star being torn apart, relevant for NS–BH and set by the Roche limit — is a separate, nonlinear endpoint. And spin-induced quadrupoles from rotation also deform the star but scale with spin, not the external field.

Crucially, tidal deformability is a bulk equilibrium property: it probes cold (T ≪ 1 MeV), β-equilibrated matter at ~1–3 ρ₀, complementary to the hot, out-of-equilibrium physics of the post-merger remnant. Its near-total suppression for black holes also makes Λ a potential 'exotic-object' discriminator, flagging boson stars, gravastars, or quark stars that would carry anomalous Love numbers.

Open questions and significance

The biggest prize is nailing the dense-matter EoS: does a phase transition to deconfined quark matter occur in neutron-star cores, and would it show up as a kink or softening in Λ(M)? A single loud BNS event in the next generation of detectors — Einstein Telescope and Cosmic Explorer — could measure Λ to a few percent, resolving the ~11–12 km radius to sub-kilometer precision and testing 'twin-star' scenarios where two stars of equal mass have different radii.

Open theoretical issues include: how large dynamical-tide and f-mode resonance corrections really are (they bias EoS inference if mismodeled); whether higher-multipole (ℓ=3, 4) Love numbers add usable information; and, on the fundamental side, whether real astrophysical black holes truly have vanishing Love numbers, since any measured nonzero value would signal new physics or an exotic compact object. Combined with NICER's X-ray radius measurements and multimessenger kilonova data, tidal deformability has turned every neutron-star merger into a precision nuclear-physics experiment conducted at the speed of light.

Tidal response of compact and non-compact bodies: Love number k₂, compactness C = GM/Rc², and dimensionless deformability Λ = (2/3)k₂C⁻⁵. Larger Λ means a more easily deformed (less compact, 'squishier') object.

Object	Compactness C	Love number k₂	Dimensionless Λ
Black hole (Schwarzschild/Kerr)	0.5	0 (exactly)	0
Neutron star, 1.4 M☉ (soft EoS)	≈ 0.18	≈ 0.08	≈ 200
Neutron star, 1.4 M☉ (stiff EoS)	≈ 0.14	≈ 0.10	≈ 700
Quark/strange star (self-bound)	≈ 0.15–0.20	≈ 0.1–0.15	few hundred–1000+
White dwarf (0.6 M☉)	≈ 10⁻⁴	≈ 0.1–0.15	≈ 10¹⁸–10¹⁹
Jupiter (planet)	≈ 10⁻⁸	≈ 0.5	astronomically large