research24 min studyLesson 910 in Physics

Test a proposed dense-matter transition against several observables

Connect equilibrium jumps and onset stability to conditional multimessenger inference, without equating a suggestive posterior with a detected phase.

Before you begin

On this page

The question for today

What you’ll understand

Match cold phases using pressure and baryon chemical potential.
Apply a sharp-transition onset criterion within its assumptions.
Explain why multimessenger evidence may remain composition-degenerate.

Neutron-star cores may contain new degrees of freedom or a transition between phases, but a small radius or a soft pressure interval is not automatically evidence for one specific microscopic state. The forward problem maps a chosen equation of state into masses, radii, and tides. The inverse problem asks which matter models remain compatible with uncertain observations and theoretical constraints. These maps are rarely one-to-one over the measured density range.

State the phase construction before interpreting its jump

At zero temperature a sharp Maxwell transition between separately neutral phases has equal pressure and baryon chemical potential on its two sides. Total energy density and baryon density can jump. The chemical-potential relation connects those jumps, so specifying an arbitrary density ratio independently from the energy and pressure data can violate equilibrium. A mixed phase with global charge neutrality, surface energy, and Coulomb effects can behave differently from this locally neutral sharp-interface idealization.

Use pressure and energy density in the same units, such as MeV per cubic femtometer, and include rest energy. A finite-temperature phase boundary needs entropy and composition information that the cold construction lacks. Nucleation barriers and slow conversion can also delay a transition relative to equilibrium. Thus a cold equilibrium discontinuity is a model of possible stellar matter, not a prediction that a merger or an accreting star follows that interface instantaneously.

P_1=P_2=P_t,\qquad \mu_B=\frac{\varepsilon_1+P_t}{n_1}=\frac{\varepsilon_2+P_t}{n_2},\qquad \frac{n_2}{n_1}=\frac{\varepsilon_2+P_t}{\varepsilon_1+P_t}

These are zero-temperature equilibrium conditions with one conserved baryon chemical potential after the stated neutrality constraints. A sharp density jump is distinct from a finite-width mixed phase or a dynamical conversion calculation.↔ Scroll wide equations horizontally.

Distinguish onset instability from a separate stable branch

For a cold nonrotating barotropic star with a sharp equilibrium transition, a sufficiently large energy-density jump destabilizes the onset of an infinitesimal new-phase core. The standard onset criterion compares the jump with the transition pressure and lower-side energy density. It assumes the original branch is stable and the usual rapid-equilibrium radial perturbation treatment. Alternative conversion time scales or additional constraints can change which perturbations are physically relevant.

An unstable onset does not guarantee a disconnected stable family or twin stars. The high-density phase must become sufficiently supportive over a later range for stability to return, and its entire equation of state must be integrated. Alford, Han and Prakash’s 2013 constant-sound-speed study shows connected, disconnected, both, and absent hybrid branches in different parameter regions. A local jump criterion classifies the immediate onset only; it does not supply the maximum mass or the observational abundance of any branch.

\frac{\Delta\varepsilon_{\mathrm{crit}}}{\varepsilon_1}=\frac12+\frac32\frac{P_t}{\varepsilon_1},\qquad \Delta\varepsilon>\Delta\varepsilon_{\mathrm{crit}}\ \Longrightarrow\ \text{unstable infinitesimal-core onset}

The criterion applies to the declared cold spherical equilibrium problem and standard phase-equilibrating radial stability treatment. Δε=ε2−ε1 is positive. Equality is the marginal case; later restabilization requires additional high-density physics.↔ Scroll wide equations horizontally.

Report what a combined observation actually excludes

High pulsar masses restrict how soft an equation of state can remain at densities actually sampled by stable stars. Tidal and pulse-profile data constrain other weighted parts of the structure. The 2024 Rutherford and collaborators analysis combines NICER, gravitational-wave, radio, and nuclear-theory inputs with specified model families. Its posterior preference depends on those families and priors. A narrow band drawn from one parameterization is not a model-independent detection of a quark core or a measured transition pressure.

A credible inversion should compare alternatives that can mimic each other: a smooth crossover, a strong first-order jump with later stiffening, or nucleonic matter with changing interactions. It should state which central densities the data probe and whether a proposed transition lies above them. Merger information adds thermal, composition, rotation, and waveform uncertainties. Dated primary analyses supply empirical constraints, while a robust phase claim needs evidence that competing causal and stable descriptions cannot account for the same observations.

Work it through

Check an equilibrium jump and its onset criterion

A sharp cold transition has $\varepsilon_1=400$ , $P_t=100$ , and $\varepsilon_2=700$ , all in $\mathrm{MeV\,fm^{-3}}$ . Find its baryon-density ratio and test the standard infinitesimal-core onset criterion. Repeat only the onset test if the upper energy density becomes eight hundred.

Chemical-potential continuity gives $n_2/n_1=(700+100)/(400+100)=800/500=1.6$ . The energy-density ratio alone would give the wrong baryon-density jump.
The actual normalized energy jump is $300/400=0.75$ . The critical value is $1/2+(3/2)(100/400)=0.875$ , so this jump does not trigger the stated immediate onset instability on an initially stable branch.
If the upper energy density is eight hundred, the normalized jump becomes one, exceeding 0.875. The standard onset is then unstable; whether a separate stable branch later exists is not determined by this comparison.

The result

The original density jump is 1.6 and its onset passes this particular jump test. The stronger example destabilizes onset but does not prove the existence of twin stars or identify a real stellar phase.

The living notebook

Explore the model

Interactive lab

Interactive lab · sampling

Distinguish statistical narrowing from model uncertainty

Reset and set p = 0.5 and n = 25, predicting sample-mean variance 0.01. Compare after a fresh Reset with p = 0.5 and n = 100, predicting variance 0.0025 and half the standard deviation. This displays exact iid binomial masses, not randomly generated experiments or a dense-matter posterior. More data within one law do not test whether an assumed stellar phase model is correct.

Exact binomial massNormal interval massHighlighted outcome / mean

Mean sample proportion: 0.5
Sample-mean variance: 0.0125
Uniform CDF error bound: ≤ 0.1062

Success probability p0.5Independent samples n20

t = 0 s

Exact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). 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 forPhase matching links pressure, energy density and baryon density.

Pause & reason

Check your understanding

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

Read the checkpoint explanation

The criterion is local to the appearance of a tiny core. A later stable branch, its maximum mass, and observational distinguishability require the rest of the equation of state and full structural analysis.

Your turn

Try it yourself

Above a hypothetical transition with $P_t=100$ and $\varepsilon_2=800$ in $\mathrm{MeV\,fm^{-3}}$ , use $P=P_t+s(\varepsilon-\varepsilon_2)$ . At energy density one thousand, compare pressures for $s=1/3$ and $s=1$ . Which local causality test do they pass, and what stellar conclusion remains unavailable?

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

Give me a hint

s is the derivative dP/dε, with the same pressure and energy-density units. Causality bounds s by one in this effective fluid model.

Compare with the solution

The pressures are $100+200/3\approx166.7$ and $300\,\mathrm{MeV\,fm^{-3}}$ . The corresponding sound speeds are $c/\sqrt3$ and c, so both pass the simple nonnegative causal-slope test, with the latter saturating it. Their density normalizations still need cold-thermodynamic matching. Neither local segment determines a stable high-mass star or a disconnected branch without its lower-density continuation, interface treatment, and a full TOV and stability calculation.

Keep the idea

Take it with you

Phase matching links pressure, energy density and baryon density.
Onset instability is weaker than a prediction of an observable twin-star family.
Multimessenger constraints retain physical and statistical model degeneracies.

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.