Pith. sign in

REVIEW 2 major objections 5 minor 49 references

Neutron stars with a dark-energy core from the Chaplygin gas

T0 review · 2 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Two-phase stars with a Chaplygin dark-energy core and an ordinary-matter crust can be dynamically stable and consistent with measured mass-radius and tidal-deformability data.

desk verdict A modest, honest toy-model study of neutron stars with a Chaplygin dark-energy core; the new tidal deformability numbers are fine, but the stability claim rests on an unargued choice of rapid phase transition. read the letter →

arxiv 2412.13568 v1 pith:2BT2QW6J submitted 2024-12-18 gr-qc astro-ph.HE

Juan M. Z. Pretel , Mariana Dutra , Sergio B. Duarte This is my paper · ORCID
classification gr-qcastro-ph.HE MSC 83C5585A15 PACS 04.40.Dg95.36.+x97.60.Jd
keywords neutronstarsdarkenergycoreChaplygingashybridradialstabilitytidaldeformabilityequationofstategeneralrelativity
topics Dark Energy
open problems Dark Energy
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper asks whether a neutron star can hide a core of dark energy, modeled by the Chaplygin gas equation of state, beneath a crust of normal matter and still behave like the neutron stars we observe. The authors build two-phase relativistic stars with a Chaplygin dark-fluid core and a polytropic ordinary-matter crust, then compute mass-radius curves, radial oscillation frequencies, and tidal deformability. Comparing with pulsar mass-radius measurements and the GW170817 tidal constraint, they find a family of models that is dynamically stable and observationally allowed. The central claim is that neutron stars with a dark-energy core are possible, in the sense that they can be radially stable and consistent with recent astrophysical data, not that observed neutron stars necessarily have such cores.

What carries the argument

The central object is the discontinuous two-phase equation of state $p(\rho) = A\rho - B/\rho$ for the Chaplygin dark-fluid core and $p = \kappa\rho^{1+1/\eta}$ for the polytropic crust, with pressure continuity at the splitting surface fixing $B$ in terms of the density jump parameter $\alpha = \rho^-_{\rm dis}/\rho^+_{\rm dis}$. The argument is carried by solving the TOV equations for equilibrium, the Gondek et al. radial-pulsation equations with junction conditions that depend on whether the phase transition is slow or rapid, and the tidal Love-number equation with its own interface junction condition. The stability and tidal results are then compared with observational mass-radius and tidal-deformability constraints.

What would settle it

A precise measurement of the radius of a neutron star at known mass that falls outside the allowed mass-radius band in Fig. 6 would settle the matter: for example, a 1.4 solar-mass star with a radius below about 10 km would lie outside the curves for the parameter values studied, ruling out these dark-energy-core configurations.

Watch

Extended reading notes

Core claim

The paper argues that neutron stars with a dark-energy core are viable: configurations with a Chaplygin-fluid core and a polytropic crust satisfy the Tolman-Oppenheimer-Volkoff equations, remain stable against radial pulsations up to the maximum-mass point, and fall inside current observational constraints. The fundamental-mode squared frequency vanishes exactly at the maximum-mass turning point, matching the standard criterion $dM/d\rho_c > 0$, while the choice between slow and rapid phase transitions at the core-crust interface changes the low-mass stability behavior. For the parameter values studied, the mass-radius curves can reach the high masses suggested by the GW190814 secondary companion, and the tidal deformabilities for $\alpha \in [0.6, 1.0]$ lie within the GW170817 bound.

Load-bearing premise

The load-bearing premise is that the core-crust boundary is a sharp discontinuity whose phase-transition speed (slow or rapid) is known, but the paper does not justify which regime a Chaplygin core would physically follow.

Editorial extensions

If this is right

  • If the central claim is correct, dark energy is not necessarily confined to cosmological scales; it could reside inside compact stars without contradicting current observations.
  • The maximum mass increases with both $\alpha$ and $A$, so a Chaplygin core with $A = 0.48$ and $\rho^+_{\rm dis} = 0.5 \times 10^{15}\,{\rm g/cm^3}$ can reach masses compatible with the GW190814 secondary component.
  • Tidal deformabilities for $\alpha \in [0.6, 1.0]$ satisfy the GW170817 bound, so future gravitational-wave events with tighter constraints can discriminate among values of $\alpha$.
  • At central densities above $10^{15}\,{\rm g/cm^3}$, the dark-energy core occupies more than 60 percent of the stellar radius, making observable properties such as radius and tidal deformability sensitive to the dark-energy sector.
  • The fundamental-mode frequency vanishes at the maximum-mass point, so the standard $dM/d\rho_c > 0$ stability criterion and the pulsation stability analysis agree for these hybrid stars.

Reading between the lines

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

  • The existence claim is conditional on which phase-transition regime really occurs at the core-crust interface; a microphysical model of the Chaplygin-to-hadronic transition would select between the slow and rapid junction conditions and could change the allowed configurations.
  • The same two-phase matching procedure could be applied to other exotic cores or to the inverted configuration of a normal-matter core with a dark-energy crust, which the paper lists as future work.
  • Tighter radius measurements from pulsar timing would directly constrain $\alpha$, because the low-mass radius depends strongly on the density jump in these models.
  • The Chaplygin core may also alter universal relations among compactness, moment of inertia, and tidal deformability, a connection the paper leaves unexplored.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The paper studies hybrid compact stars consisting of a Chaplygin dark-fluid (CDF) core and an ordinary-matter polytropic crust. The authors solve the TOV equations, the radial pulsation equations with junction conditions at the phase-splitting surface, and the tidal deformability equations for this two-phase stellar model. They examine how the density-jump ratio α = ρ_dis^-/ρ_dis^+, the CDF parameter A, and the core-edge density ρ_dis^+ affect mass-radius relations, oscillation spectra, and tidal deformability. They compare selected model curves with observations (PSR J0952-0607, GW190814, and the GW170817 tidal deformability constraint) and conclude that neutron stars with a dark-energy core are dynamically stable and consistent with these measurements, thus claiming that the existence of such stars is possible.

Significance. If the central claim holds, the paper offers a simple proof-of-principle that a two-phase star with a Chaplygin dark-energy core is a viable, if not unique, explanation for the masses, radii, and tidal deformabilities of observed compact stars. The manuscript correctly applies standard TOV, radial pulsation, and tidal deformability machinery, and it is transparent about the toy-model nature of the crust EoS. The main value is as a demonstration of possibility within a phenomenological model, rather than a realistic EoS construction. The credibility of the claim, however, rests on two load-bearing assumptions: the choice of phase-transition regime (rapid vs slow) at the interface, and the interpretation of parameter-scan agreement as observational consistency. Both issues are addressable within the manuscript's scope.

major comments (2)
  1. [§2.2, Eqs. (10)-(11), Fig. 4] The junction conditions for radial pulsations at r = R_dis are taken from Pereira et al. [42] and depend on whether the phase transition is slow or rapid. The manuscript never establishes which regime is physically realized for a Chaplygin-gas core; the interface is introduced as a phenomenological density discontinuity without a microphysical conversion model. Figure 4 shows that the two conditions disagree qualitatively at low masses: with slow transitions the low-density dM/dρ_c < 0 branch remains stable, while with rapid transitions it is unstable. Since the observational comparisons in Fig. 6 use the rapid condition, the claim that the considered stars are dynamically stable and hence can exist is not robust against an alternative but equally plausible choice of junction conditions. The authors should justify the rapid-transition assumption for a Chaplygin core or repeat the stability and observational analysis under both conditions and show that the existence claim is independent of this choice.
  2. [§3, Figs. 6-7, Conclusions] The conclusion that the model is consistent with recent astrophysical measurements rests on curves selected from a grid of free parameters {A, α, ρ_dis^+}, specifically A = 0.48, α ∈ [0.6, 1.0], and ρ_dis^+ = 0.5 and 0.8 × 10^15 g/cm^3. No statistical measure of agreement is provided: the GW170817 constraint is a broad 90% credible interval, and the pulsar and GW190814 points are shown without a quantitative distance metric in Fig. 6. Because the parameters are scanned rather than independently predicted, the stated agreement is partly built into the parameter choice. Please quantify the quality of the match (e.g., a likelihood or reduced chi-square over the relevant observational constraints) or restate the conclusion as 'can be made consistent' rather than 'is consistent'.
minor comments (5)
  1. [§1] In the Introduction, 'In Addition, in the light-cone parameterization' should read 'In addition, in the light-cone parameterization'.
  2. [§2.1, Eq. (7)] The notation ρ_dis^+ and ρ_dis^- is defined in the text, but the physical direction of the density jump is counterintuitive: α = ρ_dis^-/ρ_dis^+ with α ≤ 1, yet the core is described as the denser phase. Please clarify this notation in the text and in the figure captions for readers.
  3. [Figs. 4 and 6] Figures 4 and 6 are reproduced from the authors' own Ref. [32]. The captions state this, but the manuscript should also explicitly identify in the text which results are new in the present work (e.g., the tidal deformability analysis in Figs. 5 and 7) and confirm that reuse is compliant with the journal's copyright policy.
  4. [§2.1] The units of B are given as m^-4, while densities are quoted in g/cm^3; please provide the conversion between the geometric and physical unit systems used throughout the paper.
  5. [Abstract] The phrase 'impact of the rate of energy densities at the discontinuous surface' should be 'impact of the ratio of energy densities at the discontinuous surface'.

Circularity Check

1 steps flagged · score 4.0 of 10

The stability leg of the existence claim is imported from the authors' own prior paper (Ref. [32]), making the central conclusion partly self-citational.

  1. self citation load bearing [Sec. 3, Fig. 4 caption; abstract and Conclusions]
    "Fig. 4. Squared frequency of the fundamental vibration mode as a function of the central density (left panel) and of the gravitational mass (right plot) by using ρ+dis = 0.8 × 10^15 g/cm3, three values of A and α = 0.4 for both slow (solid lines) and rapid (dashed lines) phase transitions. Source: Taken from Ref. [32]. ... These comparisons together with the radial stability analysis show that the existence of NSs with a dark-energy core is possible."

    The central existence claim is explicitly justified by 'the radial stability analysis', but that analysis is not performed in this paper. The key stability spectra (Fig. 4) and the mass-radius curves used for the observational comparison (Fig. 6) are both captioned 'Source: Taken from Ref. [32]', i.e. the authors' own previous paper. The present manuscript re-presents those numerical results rather than re-deriving or independently verifying them. Therefore the stability part of the argument—the part that converts 'consistent with observations' into 'existence is possible'—is load-bearing on a self-citation chain: the conclusion inherits the decisive computation from the same group's earlier work, without an independent check in this paper.

full rationale

The mathematical framework is largely standard and externally cited: the TOV equations, the Gondek et al. pulsation equations, the Pereira et al. junction conditions (Eqs. 10-11), and the Postnikov et al. tidal junction condition (Eq. 16) all come from independent references, so the formalism is not defined in terms of the conclusion. The unresolved slow-vs-rapid phase-transition choice is a physical-assumption gap that affects the stability verdict, but it is a robustness/correctness concern rather than a circularity. The main circularity-related issue is the paper's reliance on the authors' own Ref. [32] for the decisive stability spectra and M-R curves; those figures are reproduced rather than re-derived here. Because the paper does add new tidal-deformability results and an explicit comparison with GW170817, GW190814, and PSR J0952-0607, the central claim retains some independent content. The observational 'agreement' is also obtained by scanning the free parameters A, alpha, and rho+_dis rather than by a sharp, parameter-free prediction; this weakens the evidential force of the comparison but is not a formal fitted-input-called-prediction circularity, since the paper presents the model as a toy model with free parameters and does not claim uniqueness. Overall, the self-citation is load-bearing for the stability component, warranting a score of 4 rather than 0-2, but the paper is not circular by construction.

Assumptions & free parameters 5 free parameters · 6 assumptions · 0 invented entities

The model rests on a small set of tunable parameters and on several domain assumptions about the equation of state and phase-transition physics. No new particles or fields are introduced. The central existence claim is thus a statement about a specific toy model, not a first-principles prediction.

free parameters (5)
  • alpha (ρ−dis/ρ+dis) = scanned in [0.4,1.0]; observational comparisons use [0.6,1.0]
    Ratio of energy densities at the core-crust interface; controls the density jump and affects mass, radius, stability, and tidal deformability.
  • A (CDF parameter) = scanned in [0.20,0.48]; compatibility cases use A=0.48
    Dimensionless constant in the Chaplygin EoS p = Aρ − B/ρ; sets the dark-energy contribution.
  • ρ+dis (core-edge energy density) = 0.8×10^15 and 0.5×10^15 g/cm³
    Energy density at the end of the dark-energy core; varied to match observational mass-radius data.
  • κ (polytropic constant) = 100 km²
    Constant in the crust EoS p = κρ^{1+1/η}; chosen as 'typical' from Refs. [33,34] to describe neutron stars.
  • η (polytropic index) = 1.0
    Exponent in the crust polytrope; set to 1.0 following Refs. [33,34].
assumptions (6)
  • standard math General relativity with a spherically symmetric perfect fluid describes neutron star structure (TOV equations, Eqs. 1-3).
    Standard framework; the paper assumes Einstein gravity and isotropic matter.
  • domain assumption The Chaplygin gas EoS p = Aρ − B/ρ is a valid effective model for dark energy in the stellar core.
    Motivated by cosmology and string theory via Refs. [4,20-22]; its applicability inside compact stars is assumed.
  • domain assumption The crust is described by a polytropic EoS p = κρ^{1+1/η} with η=1, κ=100 km².
    Adopted as a simple toy model; the paper acknowledges more realistic EoSs are needed.
  • ad hoc to paper The junction conditions at the phase-splitting interface, Eqs. (10)-(11) for slow and rapid phase transitions, are valid for radial pulsations.
    Taken from Pereira et al. [42]; the physically realized transition speed is not determined, making the stability conclusion conditional.
  • domain assumption The tidal deformability junction condition, Eq. (16), applies at the density discontinuity.
    Taken from Postnikov et al. [45] and Takátsy and Kovács [46]; standard for first-order transitions.
  • domain assumption Stability is determined by the sign of the squared fundamental eigenfrequency ω0².
    Standard criterion; the paper additionally connects ω0²=0 with the maximum-mass point.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Neutron stars with a dark-energy core from the Chaplygin gas." pith.science (2026). https://pith.science/paper/2BT2QW6J

@misc{pith2026241213568,
  author       = {Pith},
  title        = {Pith review of: Neutron stars with a dark-energy core from the Chaplygin gas},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2BT2QW6J}},
  note         = {Machine review of arXiv:2412.13568}
}
abstract

We analyze the effect of a Chaplygin dark fluid (CDF) core on neutron stars (NSs). To address this study, we focus on the relativistic structure of stellar configurations composed by a dark-energy core, described by a Chaplygin-like equation of state (EoS), and an ordinary-matter crust which is described by a polytropic EoS. We examine the impact of the rate of energy densities at the discontinuous surface, defined as $\alpha= \rho_{\rm dis}^-/\rho_{\rm dis}^+$, on the radius, total gravitational mass, oscillation spectrum and tidal deformability. Furthermore, we compare our theoretical predictions with several observational mass-radius measurements and tidal deformability constraints. These comparisons together with the radial stability analysis show that the existence of NSs with a dark-energy core is possible.

Figures

Figures reproduced from arXiv: 2412.13568 by the authors.

Figure 1
Figure 1. Mass-radius profile (left panel) and mass-central density relation (right panel) for hybrid [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Mass-radius diagram (left) and mass-central density relation (right) as in Fig. [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Percentage ratio of the radius of the discontinuous surface to the radius of the star (top [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Squared frequency of the fundamental vibration mode as a function of the central density [PITH_FULL_IMAGE:figures/full_fig_p008_4.png]
Figure 5
Figure 5. Figure 5: Dimensionless tidal deformability Λ vs gravitational mass [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]
Figure 6
Figure 6. Figure 6: Mass-Radius diagram of neutron stars with a dark-energy core (left) and oscillation spec [PITH_FULL_IMAGE:figures/full_fig_p009_6.png]
Figure 7
Figure 7. Figure 7: Variation of tidal Love number (left) and dimensionless tidal deformability (right) with [PITH_FULL_IMAGE:figures/full_fig_p010_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages

  1. [42]

    J. P. Pereira, C. V. Flores and G. Lugones, Astrophys. J. 860 (2018) 12

    work page 2018

  2. [1]

    A. G. Riess et al., Astrophys. J. 116 (1998) 1009

    work page 1998

  3. [2]

    Perlmutter et al., Astrophys

    S. Perlmutter et al., Astrophys. J. 517 (1999) 565

    work page 1999

  4. [3]

    Ratra and P

    B. Ratra and P. J. E. Peebles, Phys. Rev. D37 (1988) 3406

    work page 1988

  5. [4]

    Kamenshchik, U

    A. Kamenshchik, U. Moschella and V. Pasquier, Phys. Lett. B511 (2001) 265

    work page 2001

  6. [5]

    E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D15 (2006) 1753

    work page 2006

  7. [6]

    L. Miao, L. Xiao-Dong, W. Shuang and W. Yi, Commun. Theor. Phys.56 (2011) 525

    work page 2011

  8. [7]

    Joyce, L

    A. Joyce, L. Lombriser and F. Schmidt, Annu. Rev. Nucl. Part. Sci.66 (2016) 95

    work page 2016

Show all 49 references
  1. [8]

    Koyama, Rep

    K. Koyama, Rep. Prog. Phys.79 (2016) 046902

    2016

  2. [9]

    Montefalcone, P

    G. Montefalcone, P. J. Steinhardt and D. H. Wesley, JHEP 06 (2020) 91

    2020

  3. [10]

    Shankaranarayanan and J

    S. Shankaranarayanan and J. P. Johnson, Gen. Relativ. Gravit.54 (2022) 44

    2022

  4. [11]

    L. Wang, R. R. Caldwell, J. P. Ostriker and P. J. Steinhardt,Astrophys. J.530 (2000) 17

    2000

  5. [12]

    Douspis, A

    M. Douspis, A. Riazuelo, Y. Zolnierowski and A. Blanchard, A&A 405 (2003) 409

    2003

  6. [13]

    R. R. Caldwell, Phys. Lett. B545 (2002) 23

    2002

  7. [14]

    Rahaman et al., Gen

    F. Rahaman et al., Gen. Relativ. Gravit.44 (2012) 107

    2012

  8. [15]

    Bhar and F

    P. Bhar and F. Rahaman, Eur. Phys. J. C75 (2015) 41

    2015

  9. [16]

    Banerjee, M

    A. Banerjee, M. K. Jasim and A. Pradhan, Mod. Phys. Lett. A35 (2020) 2050071

    2020

  10. [17]

    M. C. Bento, O. Bertolami and A. A. Sen, Phys. Rev. D66 (2002) 043507

    2002

  11. [18]

    R. R. R. Reis, I. Waga, M. O. Calv˜ ao and S. E. Jor´ as,Phys. Rev. D68 (2003) 061302

    2003

  12. [19]

    L. Xu, J. Lu and Y. Wang, Eur. Phys. J. C72 (2012) 1883

    2012

  13. [20]

    Bili´ c, G

    N. Bili´ c, G. Tupper and R. Viollier, Phys. Lett. B535 (2002) 17

    2002

  14. [21]

    Bordemann and J

    M. Bordemann and J. Hoppe, Phys. Lett. B317 (1993) 315. December 19, 2024 1:27 ws-procs961x669 ws-procs961x669 page 12 12

    1993

  15. [22]

    Ogawa, Phys

    N. Ogawa, Phys. Rev. D62 (2000) 085023

    2000

  16. [23]

    Makler, S

    M. Makler, S. Q. de Oliveira and I. Waga, Phys. Lett. B555 (2003) 1

    2003

  17. [24]

    Rahaman, S

    F. Rahaman, S. Ray, A. K. Jafry and K. Chakraborty, Phys. Rev. D82 (2010) 104055

    2010

  18. [25]

    P. Bhar, M. Govender and R. Sharma, Pramana 90 (2018) 5

    2018

  19. [26]

    Panotopoulos, A

    G. Panotopoulos, A. Rinc´ on and I. Lopes, Eur. Phys. J. Plus135 (2020) 856

    2020

  20. [27]

    Tello-Ortiz et al., Eur

    F. Tello-Ortiz et al., Eur. Phys. J. C80 (2020) 371

    2020

  21. [28]

    Panotopoulos, A

    G. Panotopoulos, A. Rinc´ on and I. Lopes, Phys. Dark Univ.34 (2021) 100885

    2021

  22. [29]

    A. K. Prasad, J. Kumar and A. Sarkar, Gen. Relativ. Gravit.53 (2021) 108

    2021

  23. [30]

    J. M. Z. Pretel, Eur. Phys. J. C83 (2023) 26

    2023

  24. [31]

    O. P. Jyothilakshmi, L. J. Naik and V. Sreekanth, Eur. Phys. J. C84 (2024) 427

    2024

  25. [32]

    J. M. Z. Pretel, M. Dutra and S. B. Duarte, Phys. Rev. D109 (2024) 023524

    2024

  26. [33]

    K. D. Kokkotas and J. Ruoff, A&A 366 (2001) 565

    2001

  27. [34]

    J. M. Z. Pretel, S. E. Jor´ as, R. R. R. Reis and J. D. V. Arba˜ nil, JCAP 08 (2021) 055

    2021

  28. [35]

    Chandrasekhar, Astrophys

    S. Chandrasekhar, Astrophys. J. 140 (1964) 417

    1964

  29. [36]

    Chandrasekhar, Phys

    S. Chandrasekhar, Phys. Rev. Lett.12 (1964) 114

    1964

  30. [37]

    H. M. Vath and G. Chanmugam, A&A 260 (1992) 250

    1992

  31. [38]

    Gondek and P

    D. Gondek and P. Haensel and J. L. Zdunik, A&A 325 (1997) 217

    1997

  32. [39]

    J. M. Z. Pretel and M. F. A. da Silva, MNRAS 495 (2020) 5027

    2020

  33. [40]

    Bora and U

    J. Bora and U. Dev Goswami, MNRAS 502 (2021) 1557

    2021

  34. [41]

    B. Hong, Z. Ren, and X.-L. Mu, Chinese Phys. C46 (2022) 065104

    2022

  35. [43]

    Chatziioannou, Gen

    K. Chatziioannou, Gen. Relativ. Gravit.52 (2020) 109

    2020

  36. [44]

    Dietrich, T

    T. Dietrich, T. Hinderer and A. Samajdar, Gen. Relativ. Gravit.53 (2021) 27

    2021

  37. [45]

    Postnikov, M

    S. Postnikov, M. Prakash and J. M. Lattimer, Phys. Rev. D82 (2010) 024016

    2010

  38. [46]

    Tak´ atsy and P

    J. Tak´ atsy and P. Kov´ acs,Phys. Rev. D102 (2020) 028501

    2020

  39. [47]

    LIGO and Virgo Collab. (B. P. Abbott et al.), Phys. Rev. Lett.121 (2018) 161101

    2018

  40. [48]

    Abbott et al., Astrophys

    R. Abbott et al., Astrophys. J. Lett.896 (2020) L44

    2020

  41. [49]

    LIGO and Virgo Collab. (B. P. Abbott et al.), Phys. Rev. Lett.119 (2017) 161101

    2017

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.