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
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
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
- 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.
The abstract, sentence by sentence
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
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| 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 was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
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.