pith. sign in

arxiv: 2606.21617 · v1 · pith:UAALB5QXnew · submitted 2026-06-19 · ⚛️ nucl-th · astro-ph.HE· astro-ph.SR

Hyperonic equation of state for neutron stars: A systematic Bayesian comparison of density-dependent and non-linear relativistic mean-field models

Pith reviewed 2026-06-26 12:32 UTC · model grok-4.3

classification ⚛️ nucl-th astro-ph.HEastro-ph.SR
keywords hyperonic equation of stateneutron starsrelativistic mean-field modelsBayesian inferencehyperon-nucleon couplingsNICER observationsGW170817
0
0 comments X

The pith

Hyperons reduce neutron star maximum mass by 0.05-0.10 solar masses and increase radius at 1.4 solar masses by 0.5-0.8 km across models.

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

The paper conducts a Bayesian comparison of five relativistic mean-field models extended to include the full baryon octet. It finds that hyperons produce a systematic softening of the equation of state that lowers the maximum mass by 0.05-0.10 solar masses while raising the radius at 1.4 solar masses by 0.5-0.8 km. A characteristic dip appears in the speed of sound at two to three times nuclear saturation density where hyperons first appear. Most hyperonic equations of state remain consistent with the two-solar-mass limit and do not produce mass-radius curves with negative slope at low masses. The authors identify that a 1.8 solar mass star whose radius is larger than or similar to a 1.2 solar mass star would signal baryonic degrees of freedom beyond nucleons.

Core claim

Extending five RMF models to the baryon octet with scalar hyperon-nucleon couplings varied inside hypernuclear ranges and vector isoscalar couplings fixed by SU(6) symmetry yields a consistent reduction in maximum mass of 0.05-0.10 M_⊙ and an increase in radius at 1.4 M_⊙ of 0.5-0.8 km under simultaneous NICER, GW170817, χEFT, pQCD and nuclear saturation constraints. The speed of sound softens at 2-3 ρ_sat coinciding with hyperon onset. Inclusion of hyperons suppresses the extra flexibility of the isovector channel so that average proton fractions become independent of model choice. Only a residual number of the resulting equations of state produce a negative slope in the mass-radius relatio

What carries the argument

Bayesian sampling of density-dependent and nonlinear RMF models extended to the full baryon octet, with scalar hyperon couplings drawn from hypernuclear data ranges and vector couplings fixed by SU(6) symmetry, under joint constraints from NICER, GW170817, χEFT, pQCD and nuclear saturation properties.

If this is right

  • Hyperons cause a softening of the speed of sound at densities 2-3 times nuclear saturation density that coincides with their onset.
  • The proton fraction distribution becomes independent of isovector-channel flexibility once hyperons are included.
  • Less flexible nucleonic models show comparable or slightly higher proton fractions after hyperons are added because of overall EOS stiffening.
  • All hyperonic equations of state remain consistent with the two-solar-mass maximum mass constraint.
  • Only a residual fraction of hyperonic models produce mass-radius curves with negative slope at low masses.

Where Pith is reading between the lines

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

  • Simultaneous radius measurements at 1.2 and 1.8 solar masses could distinguish hyperonic from purely nucleonic compositions without requiring knowledge of the absolute maximum mass.
  • The sound-speed softening at 2-3 ρ_sat offers a potential target for gravitational-wave signals from neutron-star mergers that could be checked against the reported models.
  • Adding delta resonances or other exotic degrees of freedom on top of hyperons would likely shift the predicted radius difference between low- and high-mass stars and could be tested with the same Bayesian setup.

Load-bearing premise

The ranges allowed for hyperon-nucleon scalar couplings by hypernuclear data, together with SU(6) vector couplings and the listed external constraints, are sufficient to control all relevant uncertainties in the hyperonic equation of state.

What would settle it

Detection of a 1.8 solar mass neutron star whose radius is smaller than the radius of a 1.2 solar mass neutron star would contradict the claim that such a configuration signals hyperonic degrees of freedom.

Figures

Figures reproduced from arXiv: 2606.21617 by Constan\c{c}a Provid\^encia, Pedro Sanson, Tuhin Malik.

Figure 1
Figure 1. Figure 1: FIG. 1. Posterior domains for the EOS (log-log), mass-radius (stable branch only), and [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Particle fractions [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Proton fraction [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Hyperon fractions [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Posterior distributions of hyperon onset densities in [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Posterior distributions of the hyperon optical potentials [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Posterior distributions of the mass-radius slope [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
read the original abstract

A systematic Bayesian inference study of the equation of state (EOS) of dense matter with strangeness is presented, extending five relativistic mean-field (RMF) models with both constant and density-dependent couplings to include the full baryon octet. The hyperon-nucleon couplings in the scalar channel are varied within ranges informed by hypernuclear data, while vector isoscalar couplings are fixed by the SU(6) symmetry quark model. Observational constraints from NICER (PSR J0030, J0437, J0740) and GW170817, theoretical constraints from chiral effective field theory ($\chi$EFT) and perturbative QCD (pQCD), and experimental constraints from nuclear saturation properties are imposed simultaneously. We find that the inclusion of hyperons systematically reduces the maximum neutron star mass by $0.05$-$0.10,M_\odot$ across all models while increasing the radius at $1.4,M_\odot$ by $0.5$-$0.8$ km. The speed of sound exhibits a characteristic softening at densities $2$-$3,\rho_{\rm sat}$ coinciding with hyperon onset. All hyperonic models remain consistent with the $2,M_\odot$ constraint. Models with a more flexible isovector channel span a larger proton fraction when only nucleons are included. However, the extra flexibility is visibly suppressed by hyperons, meaning that the average proton distribution is independent of model flexibility when hyperons are included. Less flexible models show comparable or slightly increased proton fractions due to EOS stiffening when hyperons are included. Only a residual number of hyperonic equations of state give rise to a mass-radius curve with a negative slope at low masses. A $1.8,M_\odot$ neutron star with a radius larger than or similar to the radius of a $1.2,M_\odot$ star would provide strong evidence that the star contains baryonic degrees of freedom beyond nucleons.

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.

Referee Report

3 major / 2 minor

Summary. The paper performs a systematic Bayesian inference study of hyperonic extensions to five RMF models (both density-dependent and non-linear), varying scalar hyperon-nucleon couplings within hypernuclear bounds while fixing vector isoscalar couplings via SU(6). Simultaneous constraints from NICER (J0030, J0437, J0740), GW170817, χEFT, pQCD, and nuclear saturation properties are applied. Key findings include a systematic 0.05–0.10 M⊙ reduction in maximum mass, 0.5–0.8 km increase in R(1.4), softening of the speed of sound at 2–3 ρ_sat, suppression of isovector flexibility by hyperons, and the observation that only a residual fraction of hyperonic EOS exhibit negative low-mass M-R slopes. The central claim is that an observed R(1.8) ≳ R(1.2) would constitute strong evidence for baryonic degrees of freedom beyond nucleons.

Significance. If the uniqueness of the proposed M-R diagnostic to hyperons holds under the stated constraints, the work supplies a concrete, falsifiable observational test for strangeness in neutron-star cores and demonstrates the value of multi-model Bayesian pipelines for EOS inference. The consistent numerical shifts across model classes and the explicit incorporation of χEFT/pQCD bounds are strengths that enhance the robustness of the reported trends.

major comments (3)
  1. [coupling strategy paragraph] Coupling-strategy paragraph (abstract and methods): vector isoscalar couplings are fixed to SU(6) values while only scalar couplings are varied inside hypernuclear ranges. No sensitivity analysis to reasonable deviations from SU(6) (or alternative symmetry assumptions) is reported; such variations could permit nucleonic models to reproduce comparable or inverted low-mass radius slopes after the same NICER/GW170817/χEFT/pQCD constraints, rendering the R(1.8) ≳ R(1.2) diagnostic non-unique to hyperons.
  2. [results on proton fractions] Results section on proton fractions: the statement that hyperons suppress isovector flexibility so that the average proton distribution becomes independent of model flexibility requires quantitative support (e.g., reported variances, Kolmogorov-Smirnov distances, or explicit comparison of posterior widths between nucleonic and hyperonic ensembles).
  3. [abstract and concluding paragraph] Abstract and final paragraph: the claim that R(1.8) ≳ R(1.2) provides strong evidence for degrees of freedom beyond nucleons rests on the assertion that only a residual number of hyperonic EOS exhibit negative low-mass slopes. The manuscript must explicitly quantify the corresponding fraction (or posterior probability) for the nucleonic counterparts under identical constraints to establish the claimed specificity.
minor comments (2)
  1. [abstract] The abstract states that all models satisfy the 2 M⊙ limit but does not report the posterior probability mass above 2 M⊙ or the precise form of the likelihood function used for the NICER and GW170817 data.
  2. [results] Notation for the speed-of-sound softening (2–3 ρ_sat) should be cross-referenced to the specific density range shown in the relevant figure or table.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments highlight important points for strengthening the presentation and claims. We address each major comment below and indicate planned revisions.

read point-by-point responses
  1. Referee: [coupling strategy paragraph] Coupling-strategy paragraph (abstract and methods): vector isoscalar couplings are fixed to SU(6) values while only scalar couplings are varied inside hypernuclear ranges. No sensitivity analysis to reasonable deviations from SU(6) (or alternative symmetry assumptions) is reported; such variations could permit nucleonic models to reproduce comparable or inverted low-mass radius slopes after the same NICER/GW170817/χEFT/pQCD constraints, rendering the R(1.8) ≳ R(1.2) diagnostic non-unique to hyperons.

    Authors: The choice to fix vector isoscalar couplings via SU(6) follows the standard approach in the hyperonic RMF literature to constrain the parameter space while remaining consistent with hypernuclear phenomenology. We acknowledge that a dedicated sensitivity study to deviations from SU(6) is not included and would require an expanded analysis. In the revision we will explicitly qualify the uniqueness claim as holding under the SU(6) assumption commonly adopted in the field, add a brief discussion of this modeling choice in the methods section, and note that large deviations are disfavored by existing hypernuclear constraints. We maintain that the reported trends remain representative within this framework. revision: partial

  2. Referee: [results on proton fractions] Results section on proton fractions: the statement that hyperons suppress isovector flexibility so that the average proton distribution becomes independent of model flexibility requires quantitative support (e.g., reported variances, Kolmogorov-Smirnov distances, or explicit comparison of posterior widths between nucleonic and hyperonic ensembles).

    Authors: We agree that quantitative measures will strengthen this statement. In the revised manuscript we will report the standard deviations of the proton-fraction posteriors for the nucleonic and hyperonic ensembles across all five models and include a direct comparison of posterior widths to demonstrate the suppression of isovector flexibility by hyperons. revision: yes

  3. Referee: [abstract and concluding paragraph] Abstract and final paragraph: the claim that R(1.8) ≳ R(1.2) provides strong evidence for degrees of freedom beyond nucleons rests on the assertion that only a residual number of hyperonic EOS exhibit negative low-mass slopes. The manuscript must explicitly quantify the corresponding fraction (or posterior probability) for the nucleonic counterparts under identical constraints to establish the claimed specificity.

    Authors: We accept that the specificity of the diagnostic requires explicit comparison with the nucleonic case. In the revision we will add the fraction (and, where appropriate, posterior probability) of nucleonic EOS that exhibit negative low-mass M-R slopes under the identical set of NICER, GW170817, χEFT, pQCD and saturation constraints. This quantification will be included in both the results section and the concluding discussion to support the claim. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results from external-constraint sampling

full rationale

The paper extends five RMF models to hyperons by varying scalar couplings inside hypernuclear ranges and fixing vector isoscalar couplings via SU(6), then performs Bayesian inference against independent external datasets (NICER pulsars, GW170817, χEFT, pQCD, nuclear saturation properties). Reported outcomes (mass reduction of 0.05-0.10 M_⊙, radius increase of 0.5-0.8 km at 1.4 M_⊙, suppression of isovector flexibility, and rarity of negative low-mass M-R slopes) are direct consequences of this sampling process. The diagnostic claim about R(1.8) ≳ R(1.2) follows from the sampled ensemble rather than any redefinition or refitting of the target observables themselves. No self-citation is invoked as load-bearing justification, and no step reduces a reported prediction to a fitted input by construction.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central quantitative claims rest on the validity of the chosen coupling ranges and the SU(6) fixing rule inside the RMF framework; no new particles or forces are postulated.

free parameters (1)
  • hyperon-nucleon scalar couplings
    Varied inside ranges informed by hypernuclear data to propagate uncertainty
axioms (1)
  • domain assumption Vector isoscalar hyperon-nucleon couplings fixed by SU(6) symmetry quark model
    Standard assumption invoked when extending RMF models to the baryon octet

pith-pipeline@v0.9.1-grok · 5917 in / 1441 out tokens · 40319 ms · 2026-06-26T12:32:07.978408+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

44 extracted references · 2 canonical work pages

  1. [1]

    III A for RGOPT were used to generate the corresponding plots

    inside NS that describe two solar mass stars will present with a very large probability a positive slope of the mass-radius curves for masses below∼1.5M ⊙. Moreover, a positive slope at 1.8 solar masses would be a very strong indication of the presence of non-nucleonic degrees of freedom inside NS. The radii of different mass NS given in Table IV, and dis...

  2. [2]

    Malik, M

    T. Malik, M. Ferreira, B. K. Agrawal, and C. Providˆ encia, Relativistic Description of Dense Matter Equation of State and Compatibility with Neu- tron Star Observables: A Bayesian Approach, Astrophys. J.930, 17 (2022), arXiv:2201.12552 [nucl-th]

  3. [3]

    Typel and H

    S. Typel and H. H. Wolter, Relativistic mean field calcu- lations with density dependent meson nucleon coupling, Nucl. Phys. A656, 331 (1999)

  4. [4]

    P. Char, C. Mondal, F. Gulminelli, and M. Oertel, Gener- alized description of neutron star matter with a nucleonic relativistic density functional, Phys. Rev. D108, 103045 (2023), arXiv:2307.12364 [nucl-th]

  5. [5]

    N. K. Glendenning, Neutron Stars Are Giant Hypernu- clei?, Astrophys. J.293, 470 (1985)

  6. [6]

    Schaffner and I

    J. Schaffner and I. N. Mishustin, Hyperon rich matter in neutron stars, Phys. Rev. C53, 1416 (1996), arXiv:nucl- th/9506011

  7. [7]

    Weissenborn, D

    S. Weissenborn, D. Chatterjee, and J. Schaffner-Bielich, Hyperons and massive neutron stars: the role of hyperon potentials, Nucl. Phys. A881, 62 (2012), arXiv:1111.6049 [astro-ph.HE]

  8. [8]

    Bednarek, P

    I. Bednarek, P. Haensel, J. L. Zdunik, M. Bejger, and R. Manka, Hyperons in neutron-star cores and two- solar-mass pulsar, Astron. Astrophys.543, A157 (2012), arXiv:1111.6942 [astro-ph.SR]

  9. [9]

    Providencia and A

    C. Providencia and A. Rabhi, Interplay between the sym- metry energy and the strangeness content of neutron stars, Phys. Rev. C87, 055801 (2013), arXiv:1212.5911 [nucl-th]

  10. [10]

    L. L. Lopes and D. P. Menezes, Hypernuclear matter in a complete SU(3) symmetry group, Phys. Rev. C89, 025805 (2014), arXiv:1309.4173 [nucl-th]

  11. [11]

    Fortin, C

    M. Fortin, C. Providencia, A. R. Raduta, F. Gulminelli, J. L. Zdunik, P. Haensel, and M. Bejger, Neutron star radii and crusts: uncertainties and unified equations of state, Phys. Rev. C94, 035804 (2016), arXiv:1604.01944 [astro-ph.SR]

  12. [12]

    Vidana, D

    I. Vidana, D. Logoteta, C. Providencia, A. Polls, and I. Bombaci, Estimation of the effect of hyperonic three- body forces on the maximum mass of neutron stars, EPL 94, 11002 (2011), arXiv:1006.5660 [nucl-th]

  13. [13]

    H. J. Schulze, A. Polls, A. Ramos, and I. Vidana, Max- imum mass of neutron stars, Phys. Rev. C73, 058801 (2006)

  14. [14]

    Chatterjee and I

    D. Chatterjee and I. Vida˜ na, Do hyperons exist in the interior of neutron stars?, Eur. Phys. J. A52, 29 (2016), arXiv:1510.06306 [nucl-th]

  15. [15]

    Lonardoni, A

    D. Lonardoni, A. Lovato, S. Gandolfi, and F. Ped- eriva, Hyperon Puzzle: Hints from Quantum Monte Carlo Calculations, Phys. Rev. Lett.114, 092301 (2015), arXiv:1407.4448 [nucl-th]

  16. [16]

    Logoteta, I

    D. Logoteta, I. Vidana, and I. Bombaci, Impact of chiral hyperonic three-body forces on neutron stars, Eur. Phys. J. A55, 207 (2019), arXiv:1906.11722 [nucl-th]

  17. [17]

    Gerstung, N

    D. Gerstung, N. Kaiser, and W. Weise, Hyperon–nucleon three-body forces and strangeness in neutron stars, Eur. Phys. J. A56, 175 (2020), arXiv:2001.10563 [nucl-th]

  18. [18]

    Vidana, V

    I. Vidana, V. M. Sarti, J. Haidenbauer, D. L. Mihaylov, and L. Fabbietti, Neutron Star Properties and Fem- toscopic Constraints, Eur. Phys. J. A61, 59 (2025), arXiv:2412.12729 [nucl-th]

  19. [19]

    Fonsecaet al., Refined Mass and Geometric Measure- ments of the High-mass PSR J0740+6620, Astrophys

    E. Fonsecaet al., Refined Mass and Geometric Measure- ments of the High-mass PSR J0740+6620, Astrophys. J. Lett.915, L12 (2021), arXiv:2104.00880 [astro-ph.HE]

  20. [20]

    Antoniadiset al., A Massive Pulsar in a Com- pact Relativistic Binary, Science340, 6131 (2013), arXiv:1304.6875 [astro-ph.HE]

    J. Antoniadiset al., A Massive Pulsar in a Com- pact Relativistic Binary, Science340, 6131 (2013), arXiv:1304.6875 [astro-ph.HE]

  21. [21]

    Weissenborn, D

    S. Weissenborn, D. Chatterjee, and J. Schaffner-Bielich, Hyperons and massive neutron stars: vector repul- sion and SU(3) symmetry, Phys. Rev. C85, 065802 (2012), [Erratum: Phys.Rev.C 90, 019904 (2014)], arXiv:1112.0234 [astro-ph.HE]

  22. [22]

    A. Gal, E. V. Hungerford, and D. J. Millener, Strangeness in nuclear physics, Rev. Mod. Phys.88, 035004 (2016), arXiv:1605.00557 [nucl-th]

  23. [23]

    Fortin, S

    M. Fortin, S. S. Avancini, C. Providˆ encia, and I. Vida˜ na, Hypernuclei and massive neutron stars, Phys. Rev. C95, 065803 (2017), arXiv:1701.06373 [nucl-th]

  24. [24]

    Descouvemont, Nuclear reactions of astrophysi- cal interest, Front

    C. Providˆ encia, M. Fortin, H. Pais, and A. Rabhi, Hy- peronic stars and the nuclear symmetry energy, Fron- tiers in Astronomy and Space Sciences6, 10.3389/fs- pas.2019.00013 (2019)

  25. [25]

    Fortin, A

    M. Fortin, A. R. Raduta, S. Avancini, and C. Providˆ encia, Relativistic hypernuclear compact stars with calibrated equations of state, Phys. Rev. D 101, 034017 (2020), arXiv:2001.08036 [hep-ph]

  26. [26]

    X. Sun, Z. Miao, B. Sun, and A. Li, Astrophysical Impli- cations on Hyperon Couplings and Hyperon Star Prop- erties with Relativistic Equations of States, Astrophys. J.942, 55 (2023), arXiv:2205.10631 [astro-ph.HE]

  27. [27]

    Huang, L

    C. Huang, L. Tolos, C. Providˆ encia, and A. Watts, Constraining a relativistic mean field model using neu- tron star mass–radius measurements II: hyperonic mod- els, Mon. Not. Roy. Astron. Soc.536, 3262 (2025), arXiv:2410.14572 [astro-ph.HE]

  28. [28]

    Malik and C

    T. Malik and C. Providˆ encia, Bayesian inference of sig- natures of hyperons inside neutron stars, Phys. Rev. D 106, 063024 (2022), arXiv:2205.15843 [nucl-th]

  29. [29]

    Fortin, J

    M. Fortin, J. L. Zdunik, P. Haensel, and M. Bejger, Neu- tron stars with hyperon cores: stellar radii and equation of state near nuclear density, Astron. Astrophys.576, A68 (2015), arXiv:1408.3052 [astro-ph.SR]

  30. [30]

    Schaffner, C

    J. Schaffner, C. B. Dover, A. Gal, C. Greiner, D. J. Mil- lener, and H. Stoecker, Multiply strange nuclear systems, Annals Phys.235, 35 (1994)

  31. [31]

    Buchner, A

    J. Buchner, A. Georgakakis, K. Nandra, L. Hsu, C. Rangel, M. Brightman, A. Merloni, M. Salvato, J. Donley, and D. Kocevski, X-ray spectral modelling of the AGN obscuring region in the CDFS: Bayesian model selection and catalogue, Astron. Astrophys.564, A125 12 (2014), arXiv:1402.0004 [astro-ph.HE]

  32. [32]

    Huthet al., Constraining Neutron-Star Matter with Microscopic and Macroscopic Collisions, Nature606, 276 (2022), arXiv:2107.06229 [nucl-th]

    S. Huthet al., Constraining Neutron-Star Matter with Microscopic and Macroscopic Collisions, Nature606, 276 (2022), arXiv:2107.06229 [nucl-th]

  33. [33]

    T. E. Rileyet al., AN ICERView of PSR J0030+0451: Millisecond Pulsar Parameter Estimation, Astrophys. J. Lett.887, L21 (2019), arXiv:1912.05702 [astro-ph.HE]

  34. [34]

    Choudhuryet al., A NICER View of the Nearest and Brightest Millisecond Pulsar: PSR J0437–4715, Astro- phys

    D. Choudhuryet al., A NICER View of the Nearest and Brightest Millisecond Pulsar: PSR J0437–4715, Astro- phys. J. Lett.971, L20 (2024), arXiv:2407.06789 [astro- ph.HE]

  35. [35]

    T. E. Rileyet al., A NICER View of the Massive Pul- sar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy, Astrophys. J. Lett.918, L27 (2021), arXiv:2105.06980 [astro-ph.HE]

  36. [36]

    B. P. Abbottet al.(LIGO Scientific, Virgo), GW170817: Observation of Gravitational Waves from a Binary Neu- tron Star Inspiral, Phys. Rev. Lett.119, 161101 (2017), arXiv:1710.05832 [gr-qc]

  37. [37]

    B. P. Abbottet al.(LIGO Scientific, Virgo), GW170817: Measurements of neutron star radii and equation of state, Phys. Rev. Lett.121, 161101 (2018), arXiv:1805.11581 [gr-qc]

  38. [38]

    Komoltsev and A

    O. Komoltsev and A. Kurkela, How Perturbative QCD Constrains the Equation of State at Neutron- Star Densities, Phys. Rev. Lett.128, 202701 (2022), arXiv:2111.05350 [nucl-th]

  39. [39]

    Ferreira and C

    M. Ferreira and C. Providˆ encia, Constraining neutron star matter from the slope of the mass-radius curves, Phys. Rev. D110, 063018 (2024), arXiv:2406.12582 [nucl-th]

  40. [40]

    Albino, T

    M. Albino, T. Malik, M. Ferreira, and C. Providˆ encia, Bayesian inference of hybrid stars with large quark cores, Phys. Rev. D113, 083019 (2026), arXiv:2511.02653 [nucl-th]

  41. [41]

    Oertel, F

    M. Oertel, F. Gulminelli, C. Providˆ encia, and A. R. Raduta, Hyperons in neutron stars and supernova cores, Eur. Phys. J. A52, 50 (2016), arXiv:1601.00435 [nucl-th]

  42. [42]

    Friedman and A

    E. Friedman and A. Gal, In-medium nuclear interac- tions of low-energy hadrons, Phys. Rept.452, 89 (2007), arXiv:0705.3965 [nucl-th]

  43. [43]

    Hiyama, K

    E. Hiyama, K. Sasaki, T. Miyamoto, T. Doi, T. Hat- suda, Y. Yamamoto, and T. A. Rijken, Possible lightest Ξ Hypernucleus with Modern ΞNInteractions, Phys. Rev. Lett.124, 092501 (2020), arXiv:1910.02864 [nucl-th]

  44. [44]

    Ferreira and C

    M. Ferreira and C. Providˆ encia, Learning about neu- tron star composition from the slope of the mass- radius diagram, Phys. Rev. D112, 083058 (2025), arXiv:2506.00550 [nucl-th]