research24 min studyLesson 909 in Physics

Translate tidal response into binary constraints

Normalize Love numbers and dimensionless deformability, then distinguish a waveform constraint from a unique stellar radius.

Before you begin

On this page

The question for today

What you’ll understand

Define dimensional and dimensionless tidal deformability.
Compute a compactness-dependent tidal response.
Identify waveform, spin, common-equation-of-state and prior assumptions in inference.

Before two neutron stars merge, each responds to the gravitational field of its companion. That response modifies the binary’s orbital evolution and gravitational-wave phase. Its leading adiabatic quadrupolar contribution is summarized by tidal deformability. Extracting a number from a waveform still requires masses, spin assumptions, and a waveform model; mapping that number to a radius or core composition introduces further stellar-structure information.

Normalize the induced quadrupole

In the usual slowly varying electric-type quadrupolar convention, write the induced mass quadrupole as minus lambda times the external tidal tensor. The sign is a tensor convention, while positive lambda measures the strength of the equilibrium response. The dimensionless Love number k2 is determined by integrating tidal perturbations together with a background stellar model and matching to the exterior. It is not fixed solely by radius or compactness.

Restore SI constants to distinguish lambda from the dimensionless deformability capital Lambda. The fifth power of radius appears explicitly, giving strong radius sensitivity at fixed mass and Love number. However, the Love number also changes when the internal density profile changes. A numerical radius estimate from a measured deformability therefore needs a background equation of state or a calibrated relation with its applicable assumptions and error, rather than taking a fifth root with an arbitrary fixed k2.

Q_{ij}=-\lambda\mathcal E_{ij},\qquad \lambda=\frac{2k_2R^5}{3G},\qquad \Lambda=\frac{\lambda c^{10}}{G^4M^5}=\frac{2k_2}{3C_*^5},\quad C_*=\frac{GM}{Rc^2}

The Newtonian-normalized tidal tensor has units s⁻² and mass quadrupole kg m², so λ has units kg m² s². Λ and k2 are dimensionless. The relation uses the standard static electric quadrupolar normalization for a compact star.↔ Scroll wide equations horizontally.

Form the binary combination measured most directly

The leading binary phasing depends on a mass-weighted combination of the two stellar deformabilities. The weights depend on the two gravitational masses, so a measurement of that combination is not two independent measurements of individual deformabilities. Equal masses with the same response give the same combined value. Unequal masses and measurement correlations can admit many pairs of individual values, especially when the waveform has limited tidal signal-to-noise.

The adiabatic approximation assumes the external field changes slowly relative to the relevant stellar response, away from important resonances. Late inspiral, dynamical tides, rotation, and possible phase conversion can require additional response physics. Extracting a static parameter does not establish that a star remains exactly spherical or cold throughout merger. The inspiral baseline is usually a cold, nearly equilibrium stellar sequence, whereas the remnant requires finite-temperature matter and a different dynamical calculation.

\widetilde\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}

Masses use any common unit and cancel dimensionally. This is the conventional leading mass-weighted adiabatic quadrupolar combination. It does not include every higher-order, spin-tidal, resonant or postmerger effect.↔ Scroll wide equations horizontally.

Keep a dated inference attached to its assumptions

The LIGO–Virgo 2018 GW170817 analysis illustrates the distinction. With its specified common equation-of-state parameterization and a requirement to support masses above $1.97M_\odot$ , it reported radii of $11.9\pm1.4$ kilometers for both components at 90 percent credibility. Its alternative equation-of-state-insensitive analysis produced different radius intervals. These are dated conditional inferences, not direct ruler measurements or a current universal neutron-star radius.

Joint analyses can combine tides with radio masses, NICER pulse profiles, and nuclear theory. They must avoid counting a shared likelihood twice and should test common-equation-of-state assumptions, priors, waveform systematics, spin restrictions, and density coverage. A very precise combined radius may reflect informative assumptions as well as data. Robust dense-matter statements report which model alternatives remain compatible, especially when a sharp phase transition could alter the relation between mass, radius, and tidal response.

Work it through

Compute a dimensionless response and an equal-mass combination

A spherical stellar model has $C_*=0.2$ and $k_2=0.096$ . Find Lambda. Then place two identical copies in an equal-mass adiabatic binary and evaluate the combined deformability.

The compactness power is $0.2^5=0.00032$ , and $2k_2/3=0.064$ . Their ratio gives $\Lambda=200$ .
For equal masses M, each weighted term in the numerator contains $13M^5\Lambda$ . The sum is $26M^5\Lambda$ , while the denominator is $32M^5$ .
The prefactor gives $(16/13)(26/32)=1$ , so $\widetilde\Lambda=200$ . This algebra checks the binary normalization without adding a new radius measurement.

The result

The single-star and equal-mass binary deformabilities are both two hundred in this stipulated model. The Love number was supplied, not inferred from compactness alone.

The living notebook

Explore the model

Interactive lab

Interactive lab · matrix

Inspect deformation without an area change

Reset and set a = 1 and b = 0, predicting final determinant and both singular values equal to 1. Compare after a fresh Reset with a = 1 and b = 0.5: determinant remains 1, while singular values are about 1.2808 and 0.7808. Play interpolates geometry from identity; readouts always describe the selected final matrix. This planar shear is not a relativistic quadrupole, a Love-number integration, or a binary waveform.

Interpolated geometryOriginal unit squareSelected final transform

Area scale · det A: 1.4
Singular values · σmax, σmin: 1.597, 0.876
Condition number · κ₂(A): 1.823

Diagonal scale a1.4Shear b0.6

t = 0 s

The selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forLove number and compactness jointly determine tidal deformability.

Pause & reason

Check your understanding

Interactive checks need JavaScript. You can still read the explanation.

Read the checkpoint explanation

The waveform combination mixes both stars and depends on their masses. A radius inversion needs the internal-profile dependence of k2 and the assumptions of the chosen structural and waveform model.

Your turn

Try it yourself

Hold mass and k2 fixed in a controlled algebraic comparison, and increase a model radius by ten percent from a state with $\Lambda=200$ . Find the new Lambda. If k2 instead falls by ten percent in the comparison, find the revised result and explain the inference lesson.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

The explicit radius dependence is fifth power, but the Love number multiplies it.

Compare with the solution

At fixed k2 the ratio is $1.1^5=1.61051$ , giving $\Lambda=322.102$ . If k2 is also multiplied by 0.9, the new value is $289.8918$ . These changes are an algebraic sensitivity comparison, not a self-consistent sequence of two equations of state. They show why a radius conversion cannot ignore the internal-profile dependence of the Love number, even though the explicit fifth power makes deformability strongly radius sensitive.

Keep the idea

Take it with you

Love number and compactness jointly determine tidal deformability.
Binary phasing primarily constrains a mass-weighted response.
Dated radius inferences retain waveform, prior, and equation-of-state conditions.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.