arxiv: 2605.05816 · v1 · submitted 2026-05-07 · ✦ hep-ph · astro-ph.HE· nucl-th

Recognition: unknown

Massive hybrid stars within the extended three-flavor quark-meson diquark model

Jens O. Andersen , Mathias P. N{\o}dtvedt

Authors on Pith no claims yet

Pith reviewed 2026-05-08 08:34 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.HEnucl-th

keywords hybrid starsquark-meson diquark modelequation of stateneutron starsperturbative QCDvector mesons

0 comments

The pith

The extended three-flavor quark-meson diquark model produces hybrid stars with masses and radii matching observations when vector and axial-vector mesons are included.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper applies the extended three-flavor quark-meson diquark model to the interior of hybrid stars by computing the equation of state in the mean-field approximation at zero temperature. The model incorporates quarks, scalar and pseudoscalar mesons, diquarks, and vector and axial-vector mesons as effective degrees of freedom, with charge neutrality imposed. The equation of state is matched to a low-density nuclear equation of state, and parameters are adjusted to also satisfy perturbative QCD constraints at high density. The vector and axial-vector mesons stiffen the equation of state at intermediate densities, enabling the construction of stars whose masses and radii agree with astrophysical data. The calculation shows that stars heavier than about two solar masses contain a quark core with central density at least 3.9 times nuclear saturation density, and the speed of sound exhibits a double-peak structure before approaching the conformal limit from above.

Core claim

Within the extended three-flavor quark-meson diquark model the addition of vector and axial-vector mesons makes the equation of state sufficiently stiff to support hybrid stars whose masses and radii are in agreement with recent astrophysical observations and perturbative QCD. Stars with a mass larger than approximately 2 solar masses have a quark core with central baryon density of at least 3.9 times the nuclear saturation density. The speed of sound has a double-peak structure caused by the decrease in the strange quark mass with increasing baryon chemical potential and relaxes to the conformal limit from above at large chemical potentials.

What carries the argument

The equation of state computed in the mean-field approximation at zero temperature from the extended three-flavor quark-meson diquark model that includes quarks, scalar and pseudoscalar mesons, diquarks, and vector and axial-vector mesons.

If this is right

Hybrid stars can be constructed whose masses and radii agree with recent astrophysical observations.
Stars with masses larger than about two solar masses contain a quark core at central densities of at least 3.9 times nuclear saturation density.
The speed of sound in the dense matter exhibits a double-peak structure and approaches the conformal limit from above at high baryon chemical potentials.
The decrease of the strange quark mass with increasing baryon chemical potential produces the second peak in the speed of sound.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The model supplies a single effective framework that connects low-density nuclear matter to high-density perturbative QCD.
Very massive neutron stars are expected to contain deconfined quark matter in their cores.
Future radius measurements of stars near two solar masses can further constrain the required stiffness provided by the vector mesons.

Load-bearing premise

The mean-field approximation at zero temperature together with the specific choice of model parameters adjusted to match the low-density nuclear equation of state, perturbative QCD at high density, and astrophysical observations.

What would settle it

A precise mass-radius measurement for a star heavier than two solar masses that requires a central density below 3.9 times nuclear saturation density, or an equation of state lacking the double-peak structure in the speed of sound.

Figures

Figures reproduced from arXiv: 2605.05816 by Jens O. Andersen, Mathias P. N{\o}dtvedt.

**Figure 1.** Figure 1: Quark-, diquark- and vector condensates, and chemical potentials as view at source ↗

**Figure 2.** Figure 2: Number density as function of baryon chemical potential for all three view at source ↗

**Figure 3.** Figure 3: Normalized pressure as a function of µB for three parameter sets. The purple band is pQCD. The black solid line shows the nuclear EoS [53]. The vertical lines have the same meaning as in view at source ↗

**Figure 5.** Figure 5: Mass-radius relationship for three different parameter sets. The black line is the nuclear equation of state GMSR [53] and the dashes lines are pure quark stars with the same parameters. See main text for details. We close this section with a few remarks on future directions. So far, we have worked exclusively at T = 0. The thermodynamic potential in the mean-field approximation can be straightforwardly … view at source ↗

**Figure 4.** Figure 4: Speed of sound squared c 2 s as a function of baryon chemical potential µB for three different parameter sets. The solid horizontal line is the conformal limit c 2 s = 1 3 and the dashed line is a two-flavor EQMD example without a hybrid construction. See main text for details. Using the TOV equation with our hybrid EoS as input, we calculate the M–R relation shown in view at source ↗

read the original abstract

We discuss the properties of the extended three-flavor quark-meson diquark (EQMD) model as a renormalizable low-energy effective model for QCD. The effective degrees of freedom are quarks, scalar- and pseudoscalar mesons, diquarks, vector- and axial-vector mesons. We calculate the equation of state (EoS) in the mean-field approximation at $T=0$ imposing charge neutrality for electric and color charges. We match the EoS with a low-density nuclear equation of state. We discuss how the choice of parameters in the model affects the EoS and thereby the mass-radius for hybrid stars. We show that it is possible to construct hybrid stars whose masses and radii are in agreement with recent astrophysical observations and perturbative QCD (pQCD). The addition of vector and axial vector mesons to the quark-meson diquark is essential, since it makes the EoS sufficiently stiff for intermediate densities. Our results suggest that stars with a mass larger than $M\sim2M_{\odot}$ have a quark core with a central density $n_B\geq 3.9n_{\rm sat}$, where $n_{\rm sat}\approx0.165$fm$^{-3}$ is the saturation density. The speed of sound has a double-peak structure and relaxes to the conformal limit from above for large baryon chemical potentials $\mu_B$. This structure is caused by the decrease in the mass of the $s$ quark as $\mu_B$ increases.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The extended EQMD model with vector and axial-vector mesons yields hybrid stars above 2 solar masses and a double-peaked speed of sound, but the EoS stiffness is shaped by simultaneous fitting to nuclear, astrophysical, and pQCD constraints.

read the letter

The paper shows that adding vector and axial-vector mesons to the three-flavor quark-meson diquark model stiffens the intermediate-density EoS enough to produce hybrid stars with masses above 2 solar masses and quark cores at central densities of at least 3.9 times saturation. The speed of sound develops a double-peak structure that later approaches the conformal limit from above, driven by the dropping strange-quark mass with increasing chemical potential. This is the concrete new result: an explicit demonstration within this renormalizable effective model that such stars are possible while matching recent mass-radius observations and pQCD at high density. The authors also match the low-density part to an external nuclear EoS and impose electric and color neutrality in the mean-field T=0 calculation. That combination is useful for anyone who wants a ready-made EoS table from this class of models. The calculations are standard for the field and the logic is internally consistent. The main limitation is the parameter tuning. Meson masses and couplings are chosen to satisfy three regimes at once—low-density nuclear matter, astrophysical bounds, and high-density pQCD—so the reported stiffness and central densities are not independent predictions. Small shifts in the vector couplings or diquark gap could move the quark-core threshold or flatten the sound-speed peaks while still satisfying the boundary conditions. The abstract gives no numerical parameter set or stability checks, so the full text needs to be examined for how unique the fit is. This work is for people who already use effective QCD models for neutron-star cores and want to see how vector mesons change the hybrid-star window. A reader looking for a new mechanism or falsifiable prediction outside the fitted constraints will find less. It is coherent enough and grounded enough to deserve referee time; the numerics can be checked and the fitting procedure evaluated directly.

Referee Report

2 major / 2 minor

Summary. The manuscript presents the extended three-flavor quark-meson diquark (EQMD) model incorporating quarks, scalar/pseudoscalar mesons, diquarks, vector and axial-vector mesons. In the mean-field approximation at T=0 with electric and color charge neutrality, the EoS is computed and matched to a low-density nuclear EoS. Parameter choices are explored to construct hybrid star models with masses and radii consistent with astrophysical observations and pQCD constraints. Vector and axial-vector mesons are shown to be essential for EoS stiffness at intermediate densities, leading to the result that stars with M ≳ 2 M_⊙ have quark cores at central densities n_B ≥ 3.9 n_sat (with n_sat ≈ 0.165 fm^{-3}), and the speed of sound exhibits a double-peak structure relaxing to the conformal limit from above.

Significance. If the results hold, this work supplies a renormalizable effective QCD model that connects nuclear matter at low density, the hybrid regime, and pQCD at high density, with explicit inclusion of vector/axial-vector mesons and diquarks. The demonstration that these degrees of freedom enable M > 2 M_⊙ hybrid stars and the reported double-peak structure in c_s² constitute concrete, potentially falsifiable outputs for neutron-star phenomenology. The framework's ability to simultaneously satisfy multiple constraints is a strength, though the multi-regime tuning identified in the stress-test note reduces the independence of the predictions.

major comments (2)

[Abstract and parameter choice section] Abstract and model-parameter discussion: The central claim that vector and axial-vector mesons are 'essential' for sufficient EoS stiffness (enabling M > 2 M_⊙ stars with n_B ≥ 3.9 n_sat) rests on a specific parameter set chosen to reproduce nuclear EoS at low density, pQCD at high density, and astrophysical observations simultaneously. Because the stiffness and double-peak in c_s² are obtained after this multi-constraint adjustment, the stellar properties reduce to quantities defined by the fitted inputs rather than emerging independently; an explicit sensitivity study varying vector couplings or diquark gaps while preserving the boundary conditions is required to establish robustness.
[Stellar properties and EoS results] Results on stellar structure: The statement that stars with M ∼ 2 M_⊙ have quark cores at n_B ≥ 3.9 n_sat is presented without reported uncertainties, quantitative comparison to specific observational constraints (e.g., NICER radius posteriors or GW170817 tidal deformability), or variation under reasonable changes to the mean-field treatment. The T=0 mean-field approximation with imposed neutrality may shift the transition density and central-density threshold if fluctuations or finite-temperature effects are included.

minor comments (2)

[Abstract] The abstract refers to 'recent astrophysical observations' without citing the specific datasets (e.g., PSR J0740+6620 or NICER results); these references should appear explicitly when the mass-radius constraints are discussed.
[EoS construction] Clarify the precise matching density and any interpolation or smoothing procedure used when joining the EQMD EoS to the external nuclear EoS; this affects the location of the hybrid transition and should be shown in a dedicated figure or table.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below, providing the strongest honest defense of our results while agreeing to strengthen the presentation where appropriate.

read point-by-point responses

Referee: [Abstract and parameter choice section] Abstract and model-parameter discussion: The central claim that vector and axial-vector mesons are 'essential' for sufficient EoS stiffness (enabling M > 2 M_⊙ stars with n_B ≥ 3.9 n_sat) rests on a specific parameter set chosen to reproduce nuclear EoS at low density, pQCD at high density, and astrophysical observations simultaneously. Because the stiffness and double-peak in c_s² are obtained after this multi-constraint adjustment, the stellar properties reduce to quantities defined by the fitted inputs rather than emerging independently; an explicit sensitivity study varying vector couplings or diquark gaps while preserving the boundary conditions is required to establish robustness.

Authors: The parameter set is not freely fitted but constrained by independent physical requirements: reproduction of nuclear saturation properties at low density, asymptotic freedom and pQCD at high density, and consistency with observed hybrid-star masses. The essential role of vector and axial-vector mesons is demonstrated by direct comparison to the model variant without them, which yields an EoS too soft to support M ≳ 2 M_⊙ stars. To address robustness explicitly, we will add a sensitivity analysis in the revised manuscript, varying vector couplings and diquark gaps within the windows still compatible with the low- and high-density boundary conditions. This will confirm that the double-peak structure in c_s² and the existence of high-mass hybrid stars with quark cores persist across these variations. revision: yes
Referee: [Stellar properties and EoS results] Results on stellar structure: The statement that stars with M ∼ 2 M_⊙ have quark cores at n_B ≥ 3.9 n_sat is presented without reported uncertainties, quantitative comparison to specific observational constraints (e.g., NICER radius posteriors or GW170817 tidal deformability), or variation under reasonable changes to the mean-field treatment. The T=0 mean-field approximation with imposed neutrality may shift the transition density and central-density threshold if fluctuations or finite-temperature effects are included.

Authors: We will revise the results section to include direct, quantitative comparisons with NICER radius posteriors and GW170817 tidal-deformability bounds, together with the associated 1σ and 2σ contours on the mass-radius plane. We will also report the range of central densities obtained when the matching density to the nuclear EoS is varied within a physically motivated interval. The T=0 mean-field approximation with charge neutrality is the standard framework for this class of effective models; while fluctuations or finite-temperature corrections could in principle shift the transition density, such extensions lie outside the present scope and would require a separate study. We will add a brief discussion of these limitations and their expected direction of impact. revision: partial

Circularity Check

1 steps flagged

Parameter tuning to match observations, nuclear EoS and pQCD makes reported hybrid-star agreement a fit rather than independent derivation

specific steps

fitted input called prediction [Abstract]
"We show that it is possible to construct hybrid stars whose masses and radii are in agreement with recent astrophysical observations and perturbative QCD (pQCD). The addition of vector and axial vector mesons to the quark-meson diquark is essential, since it makes the EoS sufficiently stiff for intermediate densities. Our results suggest that stars with a mass larger than M∼2M⊙ have a quark core with a central density nB≥3.9nsat."

The model parameters are explicitly adjusted to match the low-density nuclear EoS, pQCD at high density, and the astrophysical observations simultaneously; the reported agreement with observations and the necessity of the vector mesons for sufficient stiffness are therefore achieved by construction through that multi-constraint fit rather than being derived independently from the Lagrangian.

full rationale

The paper computes the EoS in the mean-field T=0 approximation of the EQMD model after choosing parameters that simultaneously reproduce the low-density nuclear EoS, high-density pQCD constraints, and astrophysical mass-radius data. The central claim that hybrid stars with M>2M⊙ can be constructed with quark cores at nB≥3.9nsat therefore reduces to the statement that a parameter set satisfying those boundary conditions exists; the stiffness supplied by vector/axial-vector mesons and the double-peak cs² structure are direct consequences of that choice rather than emergent predictions. This is a moderate circularity (fitted-input-called-prediction) but the underlying Lagrangian and mean-field equations remain independently defined, so the score is not maximal.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 2 invented entities

The central claim rests on the mean-field approximation, charge-neutrality constraints, and multiple effective parameters whose values are selected to reproduce external data rather than derived from first principles.

free parameters (1)

meson masses and coupling constants
Multiple parameters in the EQMD model are chosen or fitted to produce an EoS that simultaneously matches astrophysical observations, the nuclear EoS, and pQCD constraints.

axioms (2)

domain assumption The mean-field approximation is valid for the EoS calculation at zero temperature.
Invoked to obtain the thermodynamic potential and pressure from the effective Lagrangian.
domain assumption Electric and color charge neutrality must hold in stellar matter.
Imposed as a constraint when minimizing the thermodynamic potential.

invented entities (2)

Diquarks no independent evidence
purpose: Effective degrees of freedom representing quark pairing in dense matter.
Introduced as part of the base quark-meson diquark model.
Vector and axial-vector mesons no independent evidence
purpose: Additional repulsion to stiffen the equation of state at intermediate densities.
Added in the present extension of the model.

pith-pipeline@v0.9.0 · 5582 in / 1775 out tokens · 41078 ms · 2026-05-08T08:34:18.852104+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Non-Parametric Equation of State Reveals Non-Conformal Behavior Beyond Neutron Star Densities
astro-ph.HE 2026-05 unverdicted novelty 6.0

Non-parametric EOS construction shows non-conformal behavior with evidence for soft quark matter and a hadron-quark phase transition in massive neutron star cores.

Reference graph

Works this paper leans on

65 extracted references · 5 canonical work pages · cited by 1 Pith paper · 1 internal anchor

[1]

Ecker et al, Nature Comm.161320 (2025)

C. Ecker et al, Nature Comm.161320 (2025)

2025
[2]

Freedman and McLerran, Phys

B.A. Freedman and McLerran, Phys. Rev. D16, 1130 (1977)

1977
[3]

B. A. Freedman and L. D. McLerran, Phys. Rev. D16, 1147 (1977)

1977
[4]

B. A. Freedman and L. D. McLerran, Phys. Rev. D16, 1169 (1977)

1977
[5]

Baluni, Phys

V . Baluni, Phys. Rev. D17, 2092 (1978)

2092
[6]

Kurkela, P

A. Kurkela, P. Romatschke and A. Vuorinen, Phys. Rev. D81105021 (2010)

2010
[7]

Kneur and L

J.-L. Kneur and L. Fernandez, Phys. Rev. D111, 034020 (2025)

2025
[8]

Gorda, R

T. Gorda, R. Paatelainen, S. Sappi and K. Seppanen, Phys. Rev. Lett.131181902 (2023)

2023
[9]

Kärkkäinen et al, Phys

A. Kärkkäinen et al, Phys. Rev. Lett.135, 021901 (2025)

2025
[10]

Bailin and A

D. Bailin and A. Love, Phys. Rept.107, 325 (1984)

1984
[11]

M. G. Alford, A. Schmitt, K. Rajagopal and T. Schäfer, Rev. Mod. Phys.80, 1455 (2008)

2008
[12]

Fukushima and T

K. Fukushima and T. Hatsuda, Rept. Prog. Phys.74, 014001 (2011)

2011
[13]

Kurkela, K

A. Kurkela, K. Rajagopal and R. Steinhorst, Phys. Rev. Lett.132, 262701 (2024)

2024
[14]

Geißel, T

A. Geißel, T. Gorda and J. Braun, Phys. Rev. Lett.135, 211901 (2025)

2025
[15]

Weinberg, Phys

S. Weinberg, Phys. Lett. B251, 288 (1990)

1990
[16]

Gasser and H

J. Gasser and H. Leutwyler, Ann. Phys.158, 142 (1984)

1984
[17]

Hebeler, J

K. Hebeler, J. M. Lattimer, C. J. Pethick and A. Schwenk, Astrophys. J.773, 11 (2013)

2013
[18]

I. Tews, J. Carlson, S. Gandolfi and S. Reddy, Astrophys. J.860, 149 (2018)

2018
[19]

Drischler, K

C. Drischler, K. Hebeler and A. Schwenk, Phys. Rev. Lett. 122, 042501 (2019)

2019
[20]

Drischler, R

C. Drischler, R. J. Furnstahl, J. A. Melendez and D. R. Phillips, Phys. Rev. Lett.125, 202702 (2020)

2020
[21]

Keller, K

J. Keller, K. Hebeler and A. Schwenk, Phys. Rev. Lett. 130, 072701 (2023)

2023
[22]

Oertel, M

M. Oertel, M. Hempel, T. Klähn, and S. Typel, Rev. Mod. Phys.89, 015007 (2017)

2017
[23]

R. L. S. Farias, G. Dallabona, G. Krein and O. A. Battistel, Phys. Rev. C73, 018201 (2006)

2006
[24]

Gholami, M

H. Gholami, M. Hofmann and M. Buballa, Phys. Rev. D 111, 014006 (2025)

2025
[25]

https://science.nasa.gov/mission/nicer/
[26]

Strodthoff, B.-J

N. Strodthoff, B.-J. Schaefer and L. von Smekal, Phys. Rev. D85, 074007 (2012)

2012
[27]

Braun, M

J. Braun, M. Leonhardt and J. M. Pawlowski SciPost Phys.6, 056 (2019)

2019
[28]

Braun and B

J. Braun and B. Schallmo, Phys. Rev. D105, 036003 (2022)

2022
[29]

Stoll, N

J. Stoll, N. Zorbach and J. Braun, arXiv: 2510.01066 [hep-ph]

work page arXiv
[30]

J. O. Andersen and M. P. Nødtvedt, Phys. Rev. D111, 034031 (2025)

2025
[31]

J. O. Andersen and M. P. Nødtvedt, Phys. Rev. D113, 014026 (2026)

2026
[32]

J. O. Andersen and M. P. Nødtvedt, arXiv: 2602.18256 [hep-ph]

work page arXiv
[33]

Gholami, L

H. Gholami et al, arXiv: 2505.22542 [hep-ph]

work page arXiv
[34]

Elitzur, Phys

S. Elitzur, Phys. Rev. D12, 3978 (1975)

1975
[35]

Watanabe and H

H. Watanabe and H. Murayama, Phys. Rev. Lett.108, 251602 (2012)

2012
[36]

Hidaka, Phys

Y . Hidaka, Phys. Rev. Lett.110, 091601 (2013)

2013
[37]

Nicolis, R

A. Nicolis, R. Penco, F. Piazza and R. A. Rosen, JHEP11, 055 (2013)

2013
[38]

Watanabe, T

H. Watanabe, T. Brauner and H. Murayama, Phys. Rev. Lett.111, 021601 (2013)

2013
[39]

Parganlija et al, Phys

D. Parganlija et al, Phys. Rev. D87, 014011 (2013)

2013
[40]

Renormalization-Group Invariant Parity-Doublet Model for Nuclear and Neutron-Star Matter

M. Recchi, L. v. Smekal and J. Wambach, arXiv: 2511.07226 [nucl-th]

work page internal anchor Pith review Pith/arXiv arXiv
[41]

Typel et al, Eur

S. Typel et al, Eur. Phys. J. A58, 221 (2022)

2022
[42]

N. K. Glendenning, Phys. Rev. D46, 1274 (1992)

1992
[43]

Baym et al

G. Baym et al. Rept. Prog. Phys.81, 056902 (2018)

2018
[44]

M. P. Nødtvedt, Code available at: https://github.com/mathiaspno/EQMD_HybridStars
[45]

V . A. Novikov et al, Nucl. Phys. B191, 301 (1981)

1981
[46]

E. S. Fraga, R. da Mata, S. Pitsingkos and A. Schmitt, Phys. Rev. D106, 074018 (2022)

2022
[47]

P. B. Demorest et al, Nature467, 1081 (2010)

2010
[48]

Antoniadis et al, Science vol

J. Antoniadis et al, Science vol. 340 NO.6131(2013)

2013
[49]

Linares et al, ApJ85954 (2018)

M. Linares et al, ApJ85954 (2018)

2018
[50]

R. W. Romani et al, ApJL934L17 (2022)

2022
[51]

T. E. Riley et al ApJL887L21 (2019)

2019
[52]

Salmi et al, ApJ974294 (2024)

T. Salmi et al, ApJ974294 (2024)

2024
[53]

Grams, J

G. Grams, J. Margueron, R. Somasundaram and S. Reddy, Eur. Phys. J. A58, 56 (2022)

2022
[54]

Annala et al, Nature Comm.14, 8451 (2023)

E. Annala et al, Nature Comm.14, 8451 (2023). 7

2023
[55]

Komoltsev and A

O. Komoltsev and A. Kurkela, Phys. Rev. Lett.128 202701 (2022)

2022
[56]

Geißel, T

A. Geißel, T. Gorda and J. Braun, Phys. Rev. D110, 014034 (2024)

2024
[57]

Fukushima and S

K. Fukushima and S. Minato, Phys. Rev. D111, 094006 (2025)

2025
[58]

Chiba and T

R. Chiba and T. Kojo, Phys. Rev. D109, 076006 (2024)

2024
[59]

B. B. Brandt et al, Phys. Rev. D112, 054038 (2025)

2025
[60]

Ayala, B

A. Ayala, B. S. Lopes, R. L. S. Farias and L. C. Parra, Phys. Lett. B864, 139396 (2025)

2025
[61]

Abbott et al, Phys

R. Abbott et al, Phys. Rev. Lett.134, 011903 (2025)

2025
[62]

Gholami et al, Phys

H. Gholami et al, Phys. Rev. D111, 103034 (2025)

2025
[63]

M. G. Alford at al, arXiv 2509.04240 [nucl-th]

work page arXiv
[64]

M. G. Alford, A. Harutyunyan, A. Sedrakian and S. Tsiopelas, Phys. Rev. D110, L061303 (2024)

2024
[65]

B. P. Abbott et al, Phys. Rev. X9, 011001 (2019). 8

2019