A Hadron-to-Quark Transition Would Leave Two Gravitational-Wave Signatures in a Neutron-Star Merger

Somewhere inside a neutron star, matter is compressed to densities several times greater than that of an atomic nucleus. At those extremes, the familiar protons and neutrons may give way to a sea of deconfined quarks and gluons. When two neutron stars spiral together and merge, the collision can push the remnant into exactly that regime. The question is whether the transition from hadronic matter to quark matter leaves a trace we can detect from Earth.

A March 2025 preprint by Sagnik Chatterjee, Shamim Haque, and colleagues uses numerical relativity to simulate that moment. Their answer: in the simulations, yes—but the signal depends on how smooth the transition is. A sharp first-order change and a gradual mixed-phase transition would write different patterns into the gravitational waves emitted after the merger.

Why Merging Neutron Stars Are Quark Laboratories

Neutron stars are not made of ordinary atomic matter. Their cores are so dense that quantum chromodynamics, the theory of the strong force, predicts quarks and gluons may exist in a deconfined state. No terrestrial observation has yet probed those densities. Binary neutron-star mergers offer a natural laboratory because the collision can briefly create an even more extreme, hot, and rapidly rotating object.

The first such merger ever detected, GW170817, gave scientists a way to study neutron-star matter through gravitational waves. Its inspiral phase helped constrain the equation of state—the relationship between pressure, density, and other properties that tells matter how to behave. But a particularly revealing moment may come after the stars collide, when the merged remnant rings down and possibly transforms inside.

A Dial Between Sharp and Smooth

The team built nine equations of state that all share the same hadronic and quark ingredients but differ in one control parameter, Δp. At Δp = 0, the phase transition is a sharp Maxwell construction: pressure stays constant while hadronic matter converts to quark matter. As Δp increases toward 0.08, the mixed phase becomes smoother, approaching a Gibbs construction in which pressure rises gradually through the transition.

The researchers then simulated equal-mass binaries with each neutron star having a mass of 1.34, 1.36, or 1.38 solar masses. In all cases, the merging stars begin as purely hadronic objects. The phase transition, if it happens, occurs during or after the merger as the remnant is compressed further.

Extra Peaks From a Hybrid Remnant

The simulations reveal a two-part gravitational-wave signature. A post-merger remnant that stays purely hadronic rings at a dominant frequency peak, labeled f2^h. If the remnant instead converts into a hypermassive hybrid star—a rapidly rotating remnant whose core contains quark matter—additional peaks appear at higher frequencies, labeled f2^q.

In the lightest mergers, the hybrid remnant produced two extra peak pairs. In intermediate-mass mergers, one extra peak appeared. In the heaviest mergers, the remnant spent less time in the hybrid state, so the extra peaks were generally weaker; the higher-Δp case remained potentially detectable in Advanced LIGO, while lower-Δp peaks might require a next-generation observatory such as the Einstein Telescope. The paper presents these as simulated signatures for future searches, not as an observation.

The Spectrogram Remembers How the Transition Happened

The nature of the transition also leaves a mark on the spectrogram, a plot of gravitational-wave frequency over time. When Δp is large and the transition is smooth, the remnant's maximum density rises gradually. For a non-collapsing remnant, the dominant frequency climbs steadily as the core moves through the mixed phase and then saturates.

When Δp = 0 and the transition is sharp, the behavior is more dramatic. The remnant undergoes a quick mini-collapse into a hybrid star, then in some cases reverses back toward hadronic matter. That brief transformation produces a hotspot in the spectrogram, a short-lived jump above the frequency where the remnant normally rings. A steady frequency rise would favor a smoother, Gibbs-like transition. A hotspot could be a signature of a sharp first-order transition.

What GW170817 Might Say About Smoothness

The simulations also link Δp to the threshold mass for prompt collapse. Above a certain total mass, the merged remnant collapses directly into a black hole rather than lingering as a neutron star. The team found that this threshold mass depends on Δp, following a hyperbolic-tangent-like trend.

If GW170817 was an equal-mass merger that promptly collapsed into a black hole, the simulations suggest the equation of state required Δp ≳ 0.04. That would disfavor a sharp Maxwell-type phase transition for that event. But the inference depends on the prompt-collapse assumption. If GW170817's remnant collapsed later, or not at all, a sharp transition could still be allowed.

The smoothness parameter also affects how much matter is flung outward. In light mergers that do not collapse, higher Δp produces more ejecta. In heavy mergers that do collapse, higher Δp produces less. The paper cautions that intermediate-mass cases are less trivial because some remnants collapse and some do not. Since ejecta can affect the kilonova afterglow, this could offer another clue in multimessenger observations.

The Caveats Around a Simulated Signal

This is a preprint based on computer simulations, not a direct detection. The authors use one fixed pair of hadronic and quark equations of state, simulate equal-mass mergers only, and leave out magnetic fields and neutrino cooling. Their treatment of thermal effects is approximate. The collapse markers are sensitive to grid resolution and additional physics, so they should be read qualitatively. The threshold-mass relation could be made more robust by sampling Δp and the threshold mass more finely.

The threshold-mass fit also depends on the underlying equation of state, so its coefficients would change with different hadronic/quark equations of state, onset of the transition, or transition strength. The results are a guide to what future detectors might see, not a claim that a phase transition has already been observed.

Listening for the Moment Matter Changes Identity

Still, the work turns a difficult question into a concrete listening strategy. If a future binary neutron-star merger produces extra high-frequency peaks or a brief frequency hotspot, scientists may be hearing matter cross a threshold that no Earth-based observation has yet reached. The signal would not be a simple bang. It would be a fingerprint of how quarks and hadrons rearrange themselves under conditions found only in the universe's most extreme stars.