Tidal Deformability of Neutron Stars

A neutron star in a binary is not a rigid point mass. Its companion's gravity varies across the star's diameter, and that gradient — the tidal field — pulls the near side harder than the far while squeezing the star transverse. The matter rearranges until pressure balances the strain, and the star settles into a slightly prolate shape carrying an induced mass quadrupole. To linear order the whole response collapses to one number, the tidal deformability $\lambda$ , defined by $Q_{ij}=-\lambda\,\mathcal{E}_{ij}$ : quadrupole equals compliance times tidal field, exactly Hooke's law for a star. The dimensionless face of $\lambda$ is the second tidal Love number $k_2$ , invented by Love in 1911 for the elastic Earth and carried into general relativity by Hinderer, Damour and Flanagan around 2008. The relativistic twist: $k_2$ is fixed by the star's entire density profile through its compactness $C=GM/Rc^2$ , so strong gravity stiffens the star and $k_2$ falls as $C$ rises, reaching exactly zero for a black hole — a nonzero deformability is the fingerprint of a material surface. What a gravitational-wave detector actually measures is $\Lambda=\tfrac{2}{3}k_2\,C^{-5}$ , and that fifth power is the whole story: because $\Lambda$ scales as $(R/M)^5$ , a one-kilometre change in radius at fixed mass moves $\Lambda$ by tens of percent — an amplification of radius differences that no other measurement resolves. The induced quadrupole feeds back on the orbit, advancing the inspiral phase by a term that enters formally at fifth post-Newtonian order but is enhanced by the large value of $\Lambda$ . When GW170817 arrived in 2017 — the first binary neutron-star merger — this phase advance turned $\Lambda$ from a theorist's parameter into a measured probe of matter above nuclear density.

Key takeaways

A companion's tidal field $\mathcal{E}_{ij}$ induces a mass quadrupole; the linear-response coefficient is the tidal deformability, $Q_{ij}=-\lambda\,\mathcal{E}_{ij}$ , or dimensionlessly the Love number $k_2$ .
The gravitational-wave observable is $\Lambda=\tfrac{2}{3}k_2\,(Rc^2/GM)^5=\tfrac{2}{3}k_2\,C^{-5}$ , whose fifth power in $1/C$ makes it exquisitely sensitive to stellar radius and hence the dense-matter equation of state.
Relativity always reduces deformability: $k_2\to0$ as $C\to\tfrac12$ , so a black hole has $\Lambda=0$ exactly and any measured $\Lambda>0$ signals extended matter.
Computing $k_2$ means integrating one $\ell=2$ perturbation equation alongside the TOV structure equations, from $y(0)=2$ at the centre to the surface.
During inspiral the induced quadrupole advances the wave phase at 5PN order, dominated by the mass-weighted combination $\tilde\Lambda$ .
GW170817 gave $\tilde\Lambda\lesssim800$ and $\Lambda_{1.4}=190^{+390}_{-120}$ , disfavouring the stiffest equations of state and, with NICER radii, anchoring multi-messenger dense-matter constraints.
Key idea: the tidal deformability measures how much a star's mass deforms into a quadrupole under an external tidal field; through $\Lambda=\tfrac{2}{3}k_2\,(Rc^2/GM)^5$ it climbs steeply with radius, converting a gravitational-wave phase measurement into a constraint on the dense-matter equation of state.

The core insight

The tidal deformability measures how much a star's mass distribution deforms into a quadrupole under an external tidal field; through the dimensionless combination $\Lambda=\tfrac{2}{3}k_2\,(Rc^2/GM)^5$ it is a steeply increasing function of stellar radius, converting a gravitational-wave phase measurement into a constraint on the dense-matter equation of state.

Where this lives today

Tidal deformability is now a standard axis of multi-messenger dense-matter physics. The LIGO–Virgo analysis of GW170817 remains the anchor: $\tilde\Lambda\lesssim800$ (90% credible, low-spin prior) and the equation-of-state-informed $\Lambda_{1.4}=190^{+390}_{-120}$ translated, through the steep $\Lambda(R)$ relation, into $R_{1.4}\lesssim13.5\,\mathrm{km}$ , ruling out the stiffest equations of state. It does not stand alone: NICER measures neutron-star radii independently from X-ray pulse profiles — PSR J0030+0451 and the $2\,M_\odot$ PSR J0740+6620 — and the two kinds of constraint are combined into joint inferences of the pressure near nuclear saturation density. A third leg is the quasi-universal "I-Love-Q" relations, which let a measured $\Lambda$ predict a moment of inertia and help break spin–tidal degeneracies. Next-generation detectors — Cosmic Explorer and the Einstein Telescope — aim to measure $\tilde\Lambda$ to a few percent per event and to catch the remnant's post-merger oscillations, where the equation of state is imprinted even more directly.

The physical picture

Hold two ideas side by side. The first is mechanical: a neutron star is a spring, the tidal field is the applied stress, and $\lambda$ is the compliance — how far the star yields per unit pull. A soft equation of state supports little pressure and gives a large, easily distorted star and a big $\Lambda$ ; a stiff one holds the star compact and rigid, giving a small $\Lambda$ . The second is geometric: the response is a quadrupole — elongation toward the companion, compression across it — because the tidal field is itself quadrupolar. The first figure fixes both, with numbers for a canonical star.

Figure

Q = −λ E: a companion's tidal field E stretches a neutron star along the line of centres and squeezes it transverse, inducing a mass quadrupole Q. The compliance is λ = (2/3)k₂R⁵. For a canonical 1.4 M⊙ star of radius R ≈ 12 km the compactness is C = GM/Rc² ≈ 0.17, the Love number is k₂ ≈ 0.08, and the dimensionless deformability is Λ = (2/3)k₂C⁻⁵ ≈ 350

Now watch what the compliance does to an observation. The induced quadrupole makes each star's gravity slightly stronger at its companion, accelerating the inspiral and advancing the gravitational-wave phase — an effect that grows steeply as the stars close in, so it is read from the last orbits before merger. The magnitude a detector infers is $\Lambda$ , which depends on radius through the fifth power of $1/C$ : that single exponent is why tidal measurements are such a sharp radius probe. The second figure is the payoff plot — the $\Lambda$ –radius relation for a $1.4\,M_\odot$ star, with GW170817 laid on top.

Figure

Representative Λ for a 1.4 M⊙ neutron star versus radius R: Λ climbs as R⁵ through the C⁻⁵ factor, so soft small-radius stars give Λ ∼ 100 while stiff large-radius stars give Λ ∼ 900. GW170817 measured Λ = 190 with a 90% range of about 70 to 580 (shaded band), bounding R ≲ 13.5 km and disfavouring the stiffest equations of state

The core equations

The defining statement is a linear-response law: the star's quadrupole is proportional to the tidal field that produces it, with a minus sign because the star deforms so as to lower its energy in the field.

Q_{ij}=-\lambda\,\mathcal{E}_{ij} \qquad (1)

where: $Q_{ij}$ — induced traceless mass-quadrupole moment; $\mathcal{E}_{ij}$ — external (electric-type) tidal tensor from the companion; $\lambda$ — tidal deformability, the linear-response coefficient. Reading: a larger $\lambda$ means a more compliant star that develops a bigger quadrupole for the same tidal pull. The dimensionful $\lambda$ is repackaged into a dimensionless Love number by pulling out five powers of the radius — the natural length in the problem. This defines $k_2$ , the second (quadrupole, $\ell=2$ ) tidal Love number.

\lambda=\frac{2}{3}k_2 R^5 \qquad (2)

where: $k_2$ — dimensionless quadrupole tidal Love number; $R$ — stellar areal radius. Reading: all the internal physics is compressed into $k_2$ ; the geometric size sits in $R^5$ . For the waveform even $R^5$ is inconvenient, because it knows mass and compactness, not kilometres. Dividing by $M^5$ (geometrized units, $G=c=1$ ) gives the observable and exposes the fifth-power compactness sensitivity that dominates everything.

\Lambda\equiv\frac{\lambda}{M^5}=\frac{2}{3}k_2\left(\frac{Rc^2}{GM}\right)^5=\frac{2}{3}k_2\,C^{-5} \qquad (3)

where: $\Lambda$ — dimensionless deformability, the primary per-star observable; $C=GM/(Rc^2)$ — compactness; $M$ — gravitational mass. Reading: two stars with almost equal $k_2$ can have very different $\Lambda$ because $C^{-5}$ ranges over more than an order of magnitude across realistic radii. The Love number is not free — it is fixed by structure. Matching interior to exterior (Derivation A) yields the closed form of Hinderer (2008), reducing everything to the compactness $C$ and one boundary number $y=R\,H'(R)/H(R)$ .

k_2=\frac{8C^5}{5}(1-2C)^2\big[2+2C(y-1)-y\big]\,D^{-1} \qquad (4)

where: $y$ — logarithmic derivative of the perturbation at the surface; $C$ — compactness; $D$ — the denominator below. Reading: the prefactor $(1-2C)^2$ and the structure of $D$ conspire so that $k_2\to0$ as $C\to\tfrac12$ — the black-hole limit is built in. The denominator carries the strong-field content. Its logarithm is the relativistic correction; expanded for small $C$ it recovers the Newtonian Love number, and at $C\to\tfrac12$ it drives $k_2$ to zero.

D=2C\big[6-3y+3C(5y-8)\big]+4C^3\big[13-11y+C(3y-2)+2C^2(1+y)\big]+3(1-2C)^2\big[2-y+2C(y-1)\big]\ln(1-2C) \qquad (5)

where: each bracket is a polynomial in $C$ and $y$ ; $\ln(1-2C)$ — the strong-field relativistic term. Reading: you do not memorize $D$ — you note that $k_2$ needs only $C$ and $y$ , and $y$ comes from integrating one extra equation through the star. That extra equation is a static $\ell=2$ , even-parity perturbation. For $H(r)$ it is a second-order ODE; the clean version is its Riccati reduction for $y=rH'/H$ , integrated outward with the TOV background (Derivation B).

r\,y'+y^2+y\,e^{\lambda(r)}\!\left[1+4\pi G r^2(p-\rho)c^{-2}\right]+r^2 Q(r)=0 \qquad (6)

where: $y=rH'/H$ ; $e^{\lambda(r)}$ — the $g_{rr}$ TOV metric function; $p,\rho$ — background pressure and density; $Q(r)$ — a coefficient from the metric, pressure, density and sound speed $\mathrm{d}p/\mathrm{d}\rho$ . Reading: integrate from $y(0)=2$ to the surface, read off $y(R)$ , and feed it into equations (4)–(5). Finally the observable enters the waveform. The two deformabilities combine into a single mass-weighted $\tilde\Lambda$ controlling the leading tidal phase, which appears at fifth post-Newtonian order (Derivation C).

\tilde\Lambda=\frac{16}{13}\frac{(m_1+12m_2)\,m_1^4\Lambda_1+(m_2+12m_1)\,m_2^4\Lambda_2}{(m_1+m_2)^5},\qquad \Psi_{\rm tidal}(f)\propto-\tilde\Lambda\,x^{5/2} \qquad (7)

where: $\tilde\Lambda$ — mass-weighted combined deformability; $m_1,m_2$ — component masses; $\Lambda_1,\Lambda_2$ — individual deformabilities; $x=(\pi GMf/c^3)^{2/3}$ — velocity parameter. Reading: for equal masses and equal $\Lambda$ the weighting collapses to $\tilde\Lambda=\Lambda$ ; the phase rides at relative order $x^5$ but is enhanced by the large value of $\tilde\Lambda$ .

Derivation ATidal deformability from matching the exterior metric

The goal is to earn the relation $\lambda=\tfrac{2}{3}k_2R^5$ rather than assert it, by placing the star in a weak external tidal field and reading the two independent pieces of its exterior potential. In the buffer region $R\ll r\ll b$ (with $b$ the separation) the metric expands in powers of $r$ : the imposed tidal field grows as $r^2$ while the induced response decays as $1/r^3$ . Isolating and comparing these two coefficients defines the linear-response constant, and nondimensionalizing by the stellar radius converts it into the Love number.

Step 1. In the star's asymptotic rest frame write $-g_{tt}/c^2-1=2\Phi/c^2$ and expand $\Phi$ in the buffer zone: a monopole $-GM/r$ , an imposed tidal piece growing as $r^2$ , and an induced quadrupole falling as $1/r^3$ . Why: three distinct pieces — the star's mass, the companion pushing in, the star pushing back — separated purely by radial power.

Step 2. Keep the two quadrupolar pieces, $\Phi=-\tfrac12\mathcal{E}_{ij}x^ix^j+\tfrac{3G}{2}\,Q_{ij}n^in^j/r^3$ with $n^i=x^i/r$ . The $r^2$ coefficient is the tidal tensor $\mathcal{E}_{ij}$ (from the companion); the $1/r^3$ coefficient is the induced quadrupole $Q_{ij}$ (from the star). Why: this is the operational definition of both — the field applied and the moment returned read off the two coefficients.

Step 3. Solve the interior perturbation equation for $H(r)$ , continue into the exterior, and demand continuity of $H$ and $H'$ at the surface. This fixes the ratio of the $1/r^3$ to the $r^2$ coefficient, i.e. $Q_{ij}/\mathcal{E}_{ij}$ ; dividing out $R^5$ makes it the dimensionless $k_2$ . Why: the deformability is an integral over the whole star, and the matching is how that interior integral becomes an exterior-observable ratio.

Step 4. Comparing $Q_{ij}=-\lambda\,\mathcal{E}_{ij}$ with the matched coefficients gives $\lambda=\tfrac23 k_2R^5$ , and dividing by $M^5$ gives $\Lambda=\tfrac23 k_2 C^{-5}$ ; taking $C\to\tfrac12$ sends $k_2\to0$ . Why: the black-hole limit is the crucial check — an object with a horizon and no material surface has nothing to deform, so its Love number must vanish.

Go deeper · Derivation B · 1 minThe perturbation equation and its Riccati reduction

Derivation BThe perturbation equation and its Riccati reduction

The aim is to reduce the linearized Einstein equations for a tidally perturbed star to a single first-order ODE integrable alongside the TOV equations in one pass, without ever tracking the perturbation's amplitude. A static, even-parity, $\ell=2$ perturbation is one radial function $H(r)$ obeying a second-order linear equation whose coefficients come from the background pressure, density, metric functions and sound speed. Because $k_2$ needs only the logarithmic derivative of $H$ at the surface, we never need $H$ itself — only its ratio to its derivative — which is exactly what a Riccati substitution isolates, replacing a second-order linear problem with a first-order nonlinear one.

Step 1. Impose the ansatz: a single function $H(r)$ multiplies the $\ell=2$ spherical harmonic. Linearizing the Einstein and fluid equations about the TOV background collapses to one second-order ODE, $H''+H'\,C_1(r)+H\,C_0(r)=0$ . Why: for a static perturbation the $\ell=2$ , even-parity sector has one physical degree of freedom, so one function suffices.

Step 2. The coefficients $C_1$ , $C_0$ are built from $p(r)$ , $\rho(r)$ , the metric functions $e^{\lambda}$ , $\nu$ , and the sound speed $\mathrm{d}p/\mathrm{d}\rho$ — all supplied by the TOV solution. Why: the same $p(\rho)$ that sets the mass–radius curve sets the deformability; no new microphysics is needed.

Step 3. Substitute $y(r)=rH'(r)/H(r)$ : the linear second-order equation becomes the nonlinear first-order Riccati equation $r\,y'+y^2+y\,F(r)+r^2Q(r)=0$ . Why: $k_2$ needs only $y(R)$ , and $y$ is invariant under $H\to cH$ , so this throws away the physically irrelevant amplitude.

Step 4. Apply boundary conditions. Regularity at the centre forces $y(0)=2$ ; integrate outward and read $y(R)$ . If the density drops discontinuously to zero at the surface, subtract $y(R)\to y(R)-4\pi R^3\rho_s/M$ before using equations (4)–(5). Why: the central value follows from $H\sim r^2$ at small $r$ , and the surface term accounts for the density jump at a sharp stellar edge.

Go deeper · Derivation C · 1 minWhy the tidal signal enters at fifth post-Newtonian order

Derivation CWhy the tidal signal enters at fifth post-Newtonian order

The goal is to see, by dimensional reasoning alone, why the tidal correction to the gravitational-wave phase sits at 5PN order yet remains measurable — the paradox that a formally tiny relativistic correction dominates the late inspiral. The strategy is to compare the energy stored in the induced quadrupole with the point-mass orbital energy, express the ratio in the velocity parameter $x\sim(M/r)$ , and note that the small power of $x$ is multiplied by the large number $\Lambda$ . It is an order-of-magnitude argument, not a full waveform calculation, but it correctly locates the term and explains its detectability.

Step 1. Estimate the induced quadrupole: the companion's tidal field is $\mathcal{E}\sim GM_2/r^3$ , so $Q\sim\lambda\mathcal{E}\sim\lambda GM_2/r^3$ . Why: this applies the defining law $Q=-\lambda\mathcal{E}$ with the Newtonian tidal field of a point companion at separation $r$ .

Step 2. Form the interaction energy, $\Delta E\sim-\tfrac12 Q\,\mathcal{E}\sim-\tfrac12\lambda(GM_2/r^3)^2$ . Why: the energy of an induced moment in the field that induced it carries the factor of one-half, as for an induced dipole in a polarizable atom.

Step 3. Divide by the point-mass orbital energy $E_{\rm pt}\sim -GM_1M_2/r$ and use $\lambda\sim\Lambda M^5$ : the ratio is $\Delta E/E_{\rm pt}\sim\Lambda\,(M/r)^5=\Lambda\,x^5$ . Why: five powers of $(M/r)$ is 5PN order, but $\Lambda\sim100$ – $1000$ multiplies it and $x$ climbs toward $0.1$ – $0.2$ in the final orbits, so the accumulated phase reaches a measurable fraction of a radian.

Limiting cases

Case	Result	Meaning
Black hole	$k_2=0,\ \Lambda=0$	A non-spinning black hole has exactly vanishing tidal deformability in general relativity; any measured $\Lambda>0$ implies a material surface.
Newtonian incompressible star	$k_2=\tfrac34$	A uniform-density Newtonian sphere has the maximal Love number $k_2=3/4$ ; central mass concentration only lowers it.
Realistic $1.4\,M_\odot$ star	$k_2\approx0.05$ – $0.15,\ \Lambda\sim100$ – $800$	Relativistic corrections push $k_2$ far below $3/4$ ; the large $C^{-5}$ factor then lifts $\Lambda$ into the hundreds.
Newtonian limit of the observable	$\Lambda=\tfrac23 k_2 C^{-5},\ C\to0$	As compactness falls, $\Lambda$ diverges — a nearly Newtonian star is "infinitely" deformable relative to its gravitational radius.

Common trapk₂ and Λ are not interchangeable, and neither is dimensionless λ

The most common slip is to treat the Love number $k_2$ and the deformability $\Lambda$ as the same thing. They are not: $\Lambda=\tfrac23 k_2 C^{-5}$ carries an extra $C^{-5}$ that ranges over more than a decade across realistic compactness, so two stars with nearly equal $k_2$ can differ by a factor of a few in $\Lambda$ — it is the $C^{-5}$ , not $k_2$ , that makes $\Lambda$ a sharp radius probe. A second slip is units: $k_2$ and $\Lambda$ are dimensionless, but the raw $\lambda$ carries dimensions (length $^5$ geometrized), so never quote a bare " $\lambda$ ". A third is degeneracy: in the waveform the tidal term is partly degenerate with spin and mass ratio, so a single-event bound on $\tilde\Lambda$ needs careful marginalization.

Worked examples

1 · The deformability of a canonical neutron star (fully worked)

A $1.4\,M_\odot$ neutron star has radius $R=12\,\mathrm{km}$ and Love number $k_2=0.08$ . Find its compactness and dimensionless tidal deformability $\Lambda$ .

Step 1. Compactness. The gravitational length of the Sun is $GM_\odot/c^2=1.477\,\mathrm{km}$ , so for $1.4\,M_\odot$ the gravitational radius is $r_g=GM/c^2=1.4\times1.477=\mathbf{2.068\,km}$ , and $C=r_g/R=2.068/12=\mathbf{0.172}$ — strongly relativistic (a black hole is $C=0.5$ ).

Step 2. The steep factor. Raise to the fifth power: $C^5=0.172^5=1.52\times10^{-4}$ , so $C^{-5}=6.58\times10^{3}$ . This is where the sensitivity lives.

Step 3. Assemble $\Lambda$ . From equation (3), $\Lambda=\tfrac23 k_2 C^{-5}=0.6667\times0.08\times6582=\mathbf{351}$ — a textbook canonical value, comfortably inside the GW170817 90% interval. Λ ≈ 350 for a 1.4 M⊙, 12 km, k₂ = 0.08 star.

2 · One kilometre of radius, amplified

Two $1.4\,M_\odot$ stars differ in radius by $1\,\mathrm{km}$ : one at $12\,\mathrm{km}$ , one at $13\,\mathrm{km}$ . How much does $\Lambda$ change?

Step 1. The compactness factor alone. Since $\Lambda\propto C^{-5}\propto R^5$ at fixed mass, the geometric ratio is $(13/12)^5=1.083^5=\mathbf{1.49}$ — a $1\,\mathrm{km}$ (about $8\%$ ) radius increase lifts $\Lambda$ by $\sim49\%$ before $k_2$ changes.

Step 2. Add the Love-number trend. A less compact star deforms more easily, so $k_2$ rises too, from $\approx0.080$ at $12\,\mathrm{km}$ to $\approx0.088$ at $13\,\mathrm{km}$ ; equation (3) then gives $\Lambda(12\,\mathrm{km})\approx351$ and $\Lambda(13\,\mathrm{km})\approx576$ .

Step 3. The lesson. The observable swings from $\approx350$ to $\approx575$ — a $\sim64\%$ change — for a radius difference no direct imaging could resolve at $\sim100\,\mathrm{Mpc}$ . A single km in R nearly doubles Λ; this fifth-power amplification is why tidal measurements pin the radius.

3 · Combining two stars into $\tilde\Lambda$

A binary has $m_1=1.6\,M_\odot$ with $\Lambda_1=200$ (heavier, more compact) and $m_2=1.2\,M_\odot$ with $\Lambda_2=700$ (lighter, puffier). Find $\tilde\Lambda$ and check the equal-mass limit.

Step 1. Numerator. With $m_1^4=6.554$ , $m_2^4=2.074$ : the first term is $(1.6+14.4)(6.554)(200)=2.097\times10^{4}$ ; the second is $(1.2+19.2)(2.074)(700)=2.961\times10^{4}$ ; their sum is $5.058\times10^{4}$ .

Step 2. Denominator and result. $(m_1+m_2)^5=2.8^5=172.1$ , so $\tilde\Lambda=\tfrac{16}{13}\times(5.058\times10^{4})/172.1=\mathbf{362}$ . It lands between $\Lambda_1$ and $\Lambda_2$ , pulled toward the lighter star because the $(m_2+12m_1)$ weighting rewards the more deformable low-mass companion.

Step 3. Sanity check. Set $m_1=m_2$ , $\Lambda_1=\Lambda_2=\Lambda$ : each bracket becomes $13m$ , the numerator is $26m^5\Lambda$ , the denominator $(2m)^5=32m^5$ , and $\tilde\Lambda=\tfrac{16}{13}\cdot\tfrac{26}{32}\Lambda=\Lambda$ — the formula reduces to a single deformability for identical stars. Λ̃ ≈ 362 here; Λ̃ = Λ exactly for equal masses — well within the GW170817 bound Λ̃ ≲ 800.

Physics of the portrait

The tidal deformability of a neutron star is a modern relativistic quantity built on a century-old idea. The $k_2$ a detector infers is, name and definition intact, the Love number introduced in 1911 to describe how the solid Earth flexes under the Moon's pull.

Origins · 1 minAugustus Edward Hough Love and how this was discovered

Augustus Edward Hough Love

1863–1940 · elasticity theory and the tidal response of the Earth

Augustus Edward Hough Love was an English mathematician who spent his career on the theory of elasticity, and it is from his 1911 monograph Some Problems of Geodynamics that the whole vocabulary of tidal response descends. Studying how the solid Earth deforms under the pull of the Moon and Sun, Love introduced three dimensionless numbers — later the Love numbers $h$ , $k$ and $l$ — quantifying the radial rise of the surface, the extra gravitational potential of the redistributed mass, and the horizontal displacement. The second, $k$ (the potential Love number), is the direct ancestor of the neutron-star $k_2$ : both measure how much a self-gravitating body's external field is altered by its own tidal deformation. For most of the twentieth century this stayed a topic in geophysics and celestial mechanics. The leap to compact objects came around 2008, when Tanja Hinderer computed the fully relativistic $k_2$ of a neutron star, Éanna Flanagan and Hinderer showed the deformability imprints on a binary's gravitational-wave phase, and Thibault Damour and Alessandro Nagar completed the strong-field formalism. A quantity Love defined for a slightly flexing Earth had become the sharpest probe of matter at several times nuclear density — realized nine years later when GW170817 measured it directly.

History

1911

— Love defines the tidal Love numbers h, k, l for the elastic Earth

1939

— Oppenheimer and Volkoff derive the relativistic equations of neutron-star structure

2008

— Flanagan and Hinderer show tidal deformability imprints on the inspiral; Hinderer computes the relativistic k₂

2009

— Damour and Nagar complete the relativistic tidal Love-number formalism

2017

— GW170817, the first binary neutron-star merger, measures the tidal signature directly

Going further

Texts. Carroll and Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 16 · Shapiro and Teukolsky, Black Holes, White Dwarfs and Neutron Stars, Ch. 9 · Glendenning, Compact Stars, Ch. 4 · Maggiore, Gravitational Waves, Vol. 1, Ch. 5. Key numerical checks: $\Lambda\approx350$ for $M=1.4\,M_\odot$ , $R=12\,\mathrm{km}$ , $k_2=0.08$ ; $(13/12)^5=1.49$ ; the equal-mass formula reduces to $\tilde\Lambda=\Lambda$ ; and the limits $k_2\to0$ as $C\to\tfrac12$ and Newtonian $k_2=3/4$ both follow from the closed-form expressions. The measured values $\tilde\Lambda\lesssim800$ , $\Lambda_{1.4}=190^{+390}_{-120}$ , and $R_{1.4}\lesssim13.5\,\mathrm{km}$ are from Abbott et al. (2018), Phys. Rev. Lett. 121, 161101.

Prerequisites. Neutron Stars · Tolman–Oppenheimer–Volkoff Equation · Equation of State of Dense Matter · Neutron Star Structure and Interior. Leads to. Equation of State of Dense Matter · Neutron Star Structure and Interior.

Ask the tutor: ask me to integrate the perturbation equation for a specific equation of state, walk through the closed-form $k_2$ term by term, or show exactly how GW170817's phase measurement became a radius bound.

Is tidal deformability the same as the star being tidally disrupted?

No — they are different regimes. Tidal deformability is a linear elastic response: the star develops a small quadrupole while remaining intact and in equilibrium, captured by the single number $\lambda$ (or $k_2$ , or $\Lambda$ ). Tidal disruption is the nonlinear endpoint, when the tidal field overwhelms self-gravity and the star is shredded. For two neutron stars the deformability governs the whole inspiral and imprints the measurable phase advance; because neutron stars are so compact, they mostly merge before disrupting. That is exactly why the clean linear deformability, not messy disruption physics, is what detectors read.

Why does relativity make k₂ smaller than the Newtonian value?

Strong-field gravity effectively stiffens the star against tidal distortion. In the Newtonian problem a uniform-density sphere reaches the maximal $k_2=3/4$ , and central mass concentration lowers it. Relativity adds further suppression: the deeper potential well binds the matter more tightly, and the closed-form $k_2$ carries a factor $(1-2C)^2$ with a $\ln(1-2C)$ term that both drive $k_2\to0$ as the compactness climbs toward the black-hole value $1/2$ . So realistic $1.4\,M_\odot$ stars have $k_2\approx0.05$ – $0.15$ , an order of magnitude below the Newtonian ceiling — yet their observable $\Lambda$ is large because of the $C^{-5}$ enhancement.

If a black hole has Λ = 0, can we tell neutron stars from black holes this way?

In principle yes — it is one of the sharpest discriminators available. A classical non-spinning black hole has $k_2=\Lambda=0$ exactly, so a binary with vanishing tidal signature is consistent with two black holes, while a clearly nonzero $\tilde\Lambda$ demands at least one material surface. In practice the measurement is hard: $\Lambda$ is small for the most compact neutron stars, the tidal phase accumulates only in the last orbits, and it is partly degenerate with spin. GW170817 favoured nonzero $\tilde\Lambda$ but was also consistent with fairly small values; future detectors with better high-frequency sensitivity will sharpen the test considerably.

Why does the tidal term enter at 5PN yet still get measured?

"5PN" describes how the tidal correction scales with orbital velocity: it is suppressed by five extra powers of $x\sim(v/c)^2$ relative to the leading point-mass dynamics, which sounds hopelessly small. Two facts rescue it. First, the coefficient is not order unity but order $\Lambda\sim100$ – $1000$ , so the term is enhanced by two-to-three orders of magnitude. Second, $x$ is not small near merger — it climbs toward $0.1$ – $0.2$ in the last several orbits, right in the detector's most sensitive band. The accumulated tidal phase then reaches a fraction of a radian, comparable to the phase resolution of the instrument. So the formal PN order undersells the physics: a high-order term with a huge coefficient, evaluated where $x$ is largest, is exactly what becomes observable.

What is the difference between λ, k₂, Λ and Λ̃?

They are four repackagings of one response. The raw $\lambda$ is the dimensionful compliance in $Q_{ij}=-\lambda\mathcal{E}_{ij}$ . The Love number $k_2=3\lambda/(2R^5)$ strips out the geometric size to give a dimensionless structure factor between $0$ (black hole) and $3/4$ (Newtonian uniform sphere). The deformability $\Lambda=\lambda/M^5=\tfrac23 k_2 C^{-5}$ is what a single star contributes to the waveform, carrying the crucial $C^{-5}$ radius amplification. And $\tilde\Lambda$ is the specific mass-weighted average of the two stars' $\Lambda$ values that appears in the leading tidal phase. Gravitational-wave papers usually quote $\tilde\Lambda$ for the system and $\Lambda_{1.4}$ for a canonical star.

Does the equation of state enter only through k₂, or also through the radius?

Through both, and that is the point. The equation of state — the pressure-density relation $p(\rho)$ above nuclear saturation — sets the entire stellar structure through the TOV equations, fixing the radius $R$ at a given mass: a stiffer (higher-pressure) equation of state gives a larger $R$ , a softer one a smaller $R$ . It also sets the internal density profile that fixes $k_2$ . But because $\Lambda=\tfrac23 k_2 C^{-5}$ with $C^{-5}\propto R^5$ , the radius dependence dominates — $\Lambda$ is far more sensitive to $R$ than to $k_2$ . That is why a measured $\Lambda$ is reported primarily as a radius constraint, and why $\Lambda$ and NICER's direct radii probe the same quantity from two independent directions.

Why is the combined Λ̃ weighted toward the lighter star?

Because the lighter star is usually the more deformable one, and the waveform is dominated by whichever star responds most. In the $\tilde\Lambda$ formula each contribution is weighted like $(m_i+12m_j)m_i^4$ : the $m_i^4$ favours the heavier star, but the $(m_i+12m_j)$ factor, with its large coefficient on the companion's mass, pulls the balance back toward the lighter one. Physically a lighter neutron star is less compact, so its $\Lambda$ is much larger through the fifth-power scaling, and its enormous deformability makes it the louder tidal voice. Worked example 3 shows this: $\tilde\Lambda=362$ sits closer to the lighter star's $\Lambda_2=700$ than a naive mass average would suggest.

Tidal Deformability of Neutron Stars

Where this lives today

The physical picture

The core equations

Limiting cases

Worked examples

1 · The deformability of a canonical neutron star (fully worked)

2 · One kilometre of radius, amplified

3 · Combining two stars into $\tilde\Lambda$

Physics of the portrait

Going further

Before this

Where this leads

Tidal Deformability of Neutron Stars

Where this lives today

The physical picture

The core equations

Limiting cases

Worked examples

1 · The deformability of a canonical neutron star (fully worked)

2 · One kilometre of radius, amplified

3 · Combining two stars into Λ~\tilde\LambdaΛ~

Physics of the portrait

Going further

Before this

Where this leads

3 · Combining two stars into $\tilde\Lambda$