Pith. Sign in

REVIEW

Effect of quark-hadron phase transitions on the neutron star mass distribution

High Energy Astrophysical Phenomena · Submitted 30 Oct 2026 · Published version · 114 references

Pith's verdict: A solid, honest 1D parameter study whose new 1.7 Msun peak claim rests on counting prematurely terminated runs as successful explosions.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on November 2, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.41242 (version 1)
Pith IDDADR3JIQ, Pith's permanent record of this paper and its review
CitationChunhua Zhu, Guoliang L\"u, Helei Liu, Renyu Luo, Shuai Zha, Sufen Guo, XuRun Huang, Zhijun Wang. "Effect of quark-hadron phase transitions on the neutron star mass distribution." arXiv:2610.41242 (2026). https://pith.science/paper/DADR3JIQ
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

This paper argues that a first-order quark-hadron phase transition inside a collapsing massive star can by itself create the high-mass peak seen in the neutron-star mass distribution.

The reading

The paper argues that a first-order quark-hadron phase transition inside a collapsing massive star can do what standard supernova models have struggled to do: produce a second, high-mass peak in the neutron-star birth mass distribution. Using one-dimensional general-relativistic collapse simulations with two hybrid equations of state that include such a transition, the authors find that successful explosions leave neutron stars clustered around $1.70\,M_\odot$ for RDF1.1 and within $1.55$--$1.70\,M_\odot$ for RDF1.2, instead of the canonical $\sim 1.4\,M_\odot$. The mechanism is a phase-transition-induced second shock: when the proto-neutron-star core crosses the transition boundary, it collapses again and bounces, and the outward shock truncates accretion at a mass set by the boundary itself. Because that boundary is a fixed feature of the equation of state, the resulting masses cluster rather than spreading smoothly, and low-compactness progenitors whose first shock fails can still explode through this second shock. If this is right, a single-star channel, without binary accretion, can account for the high-mass end of the observed mass distribution.

What carries the argument

The load-bearing object is the phase-transition boundary built into the hybrid equations of state RDF1.1 and RDF1.2. Each combines a hadronic phase, based on the DD2 relativistic mean-field description with an additional excluded-volume correction in RDF1.2, with a quark phase described by a relativistic density functional using a string-flip mechanism for confinement and vector repulsions, joined by a Gibbs construction so the transition is first order. The boundary fixes the thermodynamic state at which the proto-neutron-star core begins its second collapse; when the pure quark phase stiffens the equation of state, the infalling matter bounces and launches a second shock. That second shock does the causal work in the paper: it truncates accretion, and because the boundary is fixed, it does so at a characteristic mass, which is why the remnant masses cluster. The numerical code with its neutrino-transport scheme supplies the hydrodynamics, but the phase-transition boundary carries the explanatory weight.

Extended reading notes

Core claim

The paper's central claim is that the quark-hadron phase transition is not a marginal correction to core-collapse supernova theory but a primary agent in setting which stars explode and how massive their neutron stars are. In the simulations, the phase-transition onset density produces a threefold outcome: low-compactness progenitors explode with accretion truncated early and lower remnant masses; critical-compactness progenitors either explode or become black holes depending on whether the transition arrives before the proto-neutron-star mass exceeds the quark-core support limit; and high-compactness progenitors collapse to black holes sooner when the onset density is lower. The paper also reports a formation channel that standard neutrino-driven models lack: low-compactness progenitors whose first shock fails continue to accrete until the second collapse, and the phase-transition-driven second shock revives the explosion and leaves a relatively massive neutron star. Because the phase-transition boundary is a fixed geometric structure of the equation of state rather than a stochastic outcome, the final masses concentrate in a narrow high-mass band, producing a peak near $1.70\,M_\odot$ for RDF1.1 and $1.55$--$1.70\,M_\odot$ for RDF1.2 after weighting by the stellar initial mass function.

Load-bearing premise

The load-bearing premise is that a launched phase-transition second shock guarantees explosion and neutron-star formation, even though the simulations stop at shock radii of roughly 500 to 2000 km, before late-time fallback or delayed black-hole collapse can be followed.

Editorial extensions

If this is right

  • The observed high-mass peak in the neutron-star mass distribution could be produced by single-star evolution, without requiring binary accretion or multidimensional turbulence.
  • The peak position becomes a measurable probe of the phase-transition onset density: RDF1.1 clusters remnants near $1.70\,M_\odot$, while RDF1.2 spreads them across $1.55$--$1.70\,M_\odot$.
  • The neutron-star-to-black-hole ratio shifts with the equation of state, because the lower onset density of RDF1.2 moves the explosion-failure boundary to higher compactness.
  • Abundant low-compactness progenitors, which are favored by the stellar initial mass function, would no longer be mainly a source of canonical $1.4\,M_\odot$ neutron stars but could form the upper end of the mass distribution through the second-shock channel.
  • If binary accretion or fallback adds mass after birth, the predicted peak would move upward and could coincide with the observed high-mass peak near $\sim 1.8\,M_\odot$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A sharp prediction implicit in the paper is that any first-order phase-transition equation of state should imprint its own characteristic mass peak, so comparing the observed distribution against several hybrid equations of state could be used to select among quark-matter models.
  • The authors argue that the clustering mechanism depends on the equation-of-state boundary rather than on fluid instabilities, so the peak's location may survive in multidimensional simulations even if the success-failure boundary shifts; the peak's height is the part that still depends on the explosion-counting criterion.
  • Because the paper draws on earlier work describing a phase-transition neutrino signal, a testable extension would be to search for an electron-antineutrino burst from a nearby supernova and check whether the remnant left behind is a massive neutron star.
  • At the population level, many observed neutron stars have been recycled by accretion, so a direct test of the birth-mass peak would require separating minimally accreted systems from the rest of the sample.
Original abstract

The neutron star (NS) mass distribution serves as a crucial probe of the dense matter equation of state (EOS) and core-collapse supernova explosion mechanisms. Although observations reveal a bimodal NS mass distribution, existing theoretical simulations often struggle to naturally reproduce the high-mass peak, and the impact of the quark-hadron phase transition (PT) in high-density nuclear matter remains largely unexplored. Employing the 1D general-relativistic hydrodynamics code GR1D with hybrid EOSs incorporating a first-order PT (RDF1.1 and RDF1.2), we systematically simulate core-collapse supernovae for progenitors with zero-age main-sequence masses ranging from 13 to 120 $M_\odot$. Our results reveal a threefold impact of the PT onset density on the final supernova fate: in low-compactness progenitors, an early PT truncates accretion and produces less massive NSs; in the critical-compactness regime, it facilitates successful explosions; and in high-compactness progenitors, it accelerates black hole formation. Compared with the $\sim 1.4\,M_\odot$ peak typical of the conventional core-collapse scenario, the inclusion of the PT produces NSs with masses centered around $1.70\,M_\odot$ for RDF1.1 and within $1.55$--$1.70\,M_\odot$ for RDF1.2. Furthermore, we identify a novel formation channel whereby low-compactness progenitors with failed initial shocks can form more massive NSs via a PT-induced second shock. Due to the fixed nature of the PT boundary, the final NS masses exhibit a pronounced clustering effect. This study elucidates the critical role of the quark-hadron PT in shaping the NS mass distribution, providing important astrophysical constraints on the properties of high-density nuclear matter and the formation mechanism of massive NSs.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-11-02: 1 request, 6,756 input tokens.

What the abstract itself states, as probabilities: claims first 0.10 · states limitation 0.28 · human participants 0.01 · controlled comparison 0.70 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.90 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

The neutron star (NS) mass distribution serves as a crucial probe of the dense matter equation of state (EOS) and core-collapse supernova explosion mechanisms. Although observations reveal a bimodal NS mass distribution, existing theoretical simulations often struggle to naturally reproduce the high-mass peak, and the impact of the quark-hadron phase transition (PT) in high-density nuclear matter remains largely unexplored. Employing the 1D general-relativistic hydrodynamics code GR1D with hybrid EOSs incorporating a first-order PT (RDF1.1 and RDF1.2), we systematically simulate core-collapse supernovae for progenitors with zero-age main-sequence masses ranging from 13 to 120 $M_\odot$. Our results reveal a threefold impact of the PT onset density on the final supernova fate: in low-compactness progenitors, an early PT truncates accretion and produces less massive NSs; in the critical-compactness regime, it facilitates successful explosions; and in high-compactness progenitors, it accelerates black hole formation. Compared with the $\sim 1.4\,M_\odot$ peak typical of the conventional core-collapse scenario, the inclusion of the PT produces NSs with masses centered around $1.70\,M_\odot$ for RDF1.1 and within $1.55$--$1.70\,M_\odot$ for RDF1.2. Furthermore, we identify a novel formation channel whereby low-compactness progenitors with failed initial shocks can form more massive NSs via a PT-induced second shock. Due to the fixed nature of the PT boundary, the final NS masses exhibit a pronounced clustering effect. This study elucidates the critical role of the quark-hadron PT in shaping the NS mass distribution, providing important astrophysical constraints on the properties of high-density nuclear matter and the formation mechanism of massive NSs.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 The neutron star (NS) mass distribution serves as a crucial probe of the dense matter equation of state (EOS) and core-collapse supernova explosion mechanisms. background 1.00 0.04 0.03 asserted 1.00
s1 Although observations reveal a bimodal NS mass distribution, existing theoretical simulations often struggle to naturally reproduce the high-mass peak, and the impact of the quark-hadron phase transition (PT) in high-density nuclear matter remains largely unexplored. gap 0.97 0.06 0.03 asserted 0.98
s2 Employing the 1D general-relativistic hydrodynamics code GR1D with hybrid EOSs incorporating a first-order PT (RDF1.1 and RDF1.2), we systematically simulate core-collapse supernovae for progenitors with zero-age main-sequence masses ranging from 13 to 120 $M_\odot$. method 0.99 0.96 0.03 asserted 0.93
s3 Our results reveal a threefold impact of the PT onset density on the final supernova fate: in low-compactness progenitors, an early PT truncates accretion and produces less massive NSs; in the critical-compactness regime, it facilitates successful explosions; and in high-compactness progenitors, it accelerates black hole formation. result 1.00 0.96 0.36 asserted 1.00
s4 Compared with the $\sim 1.4\,M_\odot$ peak typical of the conventional core-collapse scenario, the inclusion of the PT produces NSs with masses centered around $1.70\,M_\odot$ for RDF1.1 and within $1.55$--$1.70\,M_\odot$ for RDF1.2. result 1.00 0.95 0.23 asserted 1.00
s5 Furthermore, we identify a novel formation channel whereby low-compactness progenitors with failed initial shocks can form more massive NSs via a PT-induced second shock. result 1.00 0.95 0.10 asserted 1.00
s6 Due to the fixed nature of the PT boundary, the final NS masses exhibit a pronounced clustering effect. result 0.86 0.86 0.10 asserted 1.00
s7 This study elucidates the critical role of the quark-hadron PT in shaping the NS mass distribution, providing important astrophysical constraints on the properties of high-density nuclear matter and the formation mechanism of massive NSs. interpretation 0.99 0.84 0.89 asserted 1.00

What kind of paper the classifier reports this is: simulation study 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. FIG. 1: Properties of compact stars described by the RDF1.1 (blue) and RDF1.2 (red) EOSs. (a) Cold mass-radius
Figure 2. FIG. 2: Time evolution of shock radius (panels a, d, g), central density (panels b, e, h), and PNS mass (defined as
Figure 3. shows the explosion simulation results of pre￾supernova models with MZAMS ranging from 13 to 40 M⊙ (the complete data for 13–120 M⊙ are listed in Ap￾pendix B). The horizontal axis represents the initial mass MZAMS, and the vertical axis denotes the compactness parameter ξ1.75. In
Figure 4. FIG. 4: Probability density function (PDF) of the neutron star gravitational mass distribution for the RDF1.1 (left
Figure 5. FIG. 5: The left panel shows the gravitational mass
Figure 6. FIG. 6: Time evolution of shock radius (panels a), central density (panels b), and PNS mass (defined as the

Figures reproduced from arXiv:2610.41242 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.