Pith. sign in

REVIEW 4 major objections 5 minor 35 references

At a critical magnetic charge q_c, neutron stars with nonlinear magnetic monopoles enter frozen states bounded by a critical horizon, beyond which no static solutions exist.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review

2026-08-04 19:15 UTC pith:MTXZDVTV

load-bearing objection A plausible extension of frozen-state boson stars to fluid neutron stars, whose central claim rests on an unexamined singular limit. the 4 major comments →

arxiv 2509.09338 v1 pith:MTXZDVTV submitted 2025-09-11 gr-qc astro-ph.HEhep-th

Frozen Neutron Stars

classification gr-qc astro-ph.HEhep-th
keywords neutron starsnonlinear electrodynamicsmagnetic monopolesfrozen statesBardeen modelHayward modelcritical horizonmodified TOV equations
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper tries to establish that a neutron star containing a nonlinear magnetic monopole cannot hold arbitrarily large magnetic charge. In the Bardeen and Hayward nonlinear electrodynamics models, the modified hydrostatic-equilibrium equations have static perfect-fluid solutions only up to a critical charge q_c; at q_c the metric component 1/g_rr drops to nearly zero at the stellar surface and -g_tt also approaches zero inside, so the surface acts like a critical horizon and all matter is confined within it. The authors call this endpoint a frozen neutron star, because a distant observer would see the surface redshift and time dilation diverge, much like an extremal black hole, but no event horizon forms. If correct, the result extends frozen states from purely bosonic stars to ordinary baryonic matter, and it identifies a new possible endpoint for neutron stars that accumulate magnetic monopoles.

Core claim

The central claim is that the family of static spherical neutron-star solutions in the Einstein-Bardeen and Einstein-Hayward models terminates at a critical magnetic charge q_c. As q approaches q_c from below, the star contracts, a dense 'hard-candy' surface layer forms, and the minimum of 1/g_rr moves to the stellar surface while its value falls to 10^-9 to 10^-12; concurrently -g_tt approaches zero throughout the interior. The authors identify this degenerate surface as a critical horizon and interpret the configuration as a frozen state. They verify the effect for three equations of state (BSk19, SLy4, AP4), finding that softer EOSs and higher central densities lower q_c, and that the cau

What carries the argument

The key object is the modified TOV system: the Einstein equations sourced by a perfect-fluid stress tensor plus the stress tensor of a nonlinear magnetic monopole in the Bardeen or Hayward Lagrangian, together with the vector ansatz A = q cos(theta) dphi. In these models the magnetic charge provides a negative pressure that compresses the star and an extra gravitational potential that deepens the minimum of the radial metric function. The critical horizon emerges when that minimum coincides with the stellar surface and -g_tt simultaneously vanishes there, an endpoint controlled by the equation of state via the causality limit on pressure.

Load-bearing premise

The argument rests on treating the numerical breakdown of the static TOV equations at q_c, where the metric becomes nearly degenerate, as a genuine equilibrium state of matter plus nonlinear field rather than a coordinate artifact or an artifact of the static perfect-fluid idealization.

What would settle it

A radial linear stability calculation for the q_c Bardeen and Hayward solutions: an unstable fundamental mode, or any regular static solution with q slightly above q_c, would show that the frozen endpoint is not physical. Alternatively, a fully relativistic simulation approaching q_c could reveal whether the near-degenerate metric persists or collapses.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

Share X Bluesky LinkedIn Reddit HN

If this is right

  • At q_c a neutron star becomes a frozen star with a critical horizon; above q_c no static perfect-fluid equilibrium exists, so q_c is the endpoint of the static branch.
  • To a distant observer the frozen star mimics an extremal black hole because surface time dilation diverges, but it has no event horizon, offering a horizonless alternative for ultracompact objects.
  • The critical charge is not universal: softer equations of state and higher central densities make q_c smaller, while stiff EOSs can fail to reach it altogether.
  • In the frozen limit the total ADM mass is dominated by the nonlinear magnetic field's charge contribution, not the baryonic mass of the star.
  • Magnetic charge deforms the mass-radius and ADM-mass-radius relations, so a frozen neutron star could masquerade as a more massive or more compact object in observations.

Where Pith is reading between the lines

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

  • The paper leaves radial stability unexamined; a linear perturbation analysis is the direct way to tell whether the frozen state is a stable equilibrium or a transient that collapses.
  • The absence of static solutions above q_c could reflect a change of topology or a dynamical collapse rather than a fundamental no-go; a time-dependent simulation would determine which.
  • Because a frozen star mimics an extremal black hole externally, gravitational-wave measurements of tidal deformability or post-merger signals could be a practical way to distinguish horizonless frozen stars from black holes.
  • If the effect is generic, similar critical horizons may appear for other nonlinear electrodynamics or in modified gravity, but that is an extension the paper only gestures toward.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

Summary. The paper studies static, spherically symmetric neutron stars in Einstein gravity coupled to nonlinear electrodynamics (Bardeen and Hayward models) and to a perfect fluid, using three nuclear equations of state (BSk19, SLy4, AP4). The authors derive modified TOV equations, integrate them numerically for fixed coupling schemes, and report that beyond a critical magnetic charge q_c no physically meaningful solutions exist. At q_c, the metric functions 1/g_rr and -g_tt approach zero at the stellar surface; they interpret this as the formation of a 'critical horizon' and a transition to a 'frozen state.' The paper also presents radial pressure profiles, compactness and average density versus q, mass-radius relations, and the ADM mass split between matter and magnetic charge.

Significance. If the central claim were firmly established, this paper would meaningfully extend the frozen-star concept from boson stars to ordinary-matter compact objects and open a new direction in neutron-star structure with nonlinear magnetic monopoles. The authors' construction of the modified TOV equations from a well-defined action, the use of three realistic equations of state, and the explicit numerical exploration of the q parameter space are strengths. However, the frozen-state interpretation rests on a singular numerical limit: the endpoint at q_c is identified by the solver ceasing to return solutions, and the paper explicitly acknowledges that stability under radial perturbations is unexamined. These gaps leave the physical reality of the frozen state as an unverified inference rather than a demonstrated equilibrium configuration.

major comments (4)
  1. [Sec. III, Figs. 3 and 4] The claim that q_c corresponds to a frozen state is based on the observation that, at the critical charge, 1/g_rr and -g_tt reach values as low as 10^-9 to 10^-12 at the surface. These are near-zero values along a sequence, not an exact solution at q_c. The paper does not construct the limiting configuration directly or show that the metric functions and matter variables converge to a regular solution of the field equations as q→q_c. The numerical failure beyond q_c could equally signal that the static perfect-fluid branch terminates in a singular or unstable configuration, as in the usual TOV endpoint. A direct construction of the limiting solution, or at least a clear verification that the q→q_c limit satisfies the Einstein equations in a distributional or regularized sense, is needed before identifying q_c with a physical transition.
  2. [Eq. (19)] The modified TOV equation contains a denominator 2e^{-2β}. As e^{-2β(R)}→0 at the surface, p'(r) will diverge at r=R unless the numerator in Eq. (19) vanishes to higher order in the same limit. The paper does not check this regularity condition. Without such a check, the near-singular sequence at q_c is not demonstrated to be a solution of the assumed perfect-fluid equations. This is load-bearing because the frozen-state interpretation requires the endpoint to be a valid hydrostatic equilibrium, not merely a place where the solver stops.
  3. [Sec. IV] The authors state that 'a rigorous analysis of the stability of these frozen neutron stars under radial perturbations is essential.' For a claimed new equilibrium phase, linear stability is a necessary condition for physical relevance. The paper does not provide even a preliminary stability analysis, nor does it discuss whether the degenerate surface at g_tt=0 is compatible with a static perfect fluid that is required to have a timelike 4-velocity. The 4-velocity remains normalized by construction, but the surface is a Killing horizon, and the behavior of the fluid at that horizon is not addressed. This is a central gap, not a peripheral omission.
  4. [Tables I and II] The critical charges q_c are reported as definite numbers determined by the existence of numerical solutions. No convergence criterion, tolerance, or error estimate is stated. Since the frozen endpoint is reached only asymptotically as e^{-2β}→0, the quoted values depend on how 'no physical solution' is operationally defined. A precise definition of the numerical critical charge and a sensitivity check with respect to solver tolerances would strengthen the claim that q_c is a property of the solution space rather than an artifact of the integration scheme.
minor comments (5)
  1. [Eq. (9)] The fluid 4-velocity should be U^μ=(e^{-α},0,0,0) for the metric signature in Eq. (7); the text appears to omit the minus sign in the exponent. Please check and correct.
  2. [Introduction] Typo: 'gracitational' should be 'gravitational'.
  3. [Sec. III C] Typo: 'ttherefore' should be 'therefore'.
  4. [References] Reference [33] is incomplete: it lacks a journal, volume, and arXiv identifier. Also, the in-text references to Refs. [19-21] for the frozen-boson-star concept would be easier to follow if the connection to the present neutron-star case were stated explicitly in Sec. II.
  5. [Fig. 1 and Fig. 2] In the fixed-s panels, the text says 'as q decreases ... the maximum central pressure increases while the stellar radius decreases,' but this behavior is opposite to the fixed-sq^2 case. Please ensure the caption and the main text clearly distinguish the two schemes and that the direction of the horizontal axis in each panel is unambiguous.

Circularity Check

1 steps flagged

Central numerical finding (critical q_c with degenerate metric limit) is computed from the paper's own modified TOV equations and is not a fitted quantity; the 'frozen state / critical horizon' interpretation imports terminology from the authors' own boson-star papers, but the derivation itself is not reduced to its inputs.

specific steps
  1. self citation load bearing [Section III.C (Radial Metric Profiles), paragraph following Fig. 4]
    "In summary, the position of this 1/grr minimum exhibits properties similar to those of the critical horizon (as discussed in Ref. [19–21]), and can be identified as the critical horizon of frozen neutron stars."

    The paper's interpretive conclusion that the degenerate-metric limit at q_c constitutes a 'frozen state' with a 'critical horizon' is imported from the authors' own prior results on Bardeen/Hayward boson stars (Refs. [19,21], which include present author Y.Q. Wang). The underlying numerical content — 1/g_rr reaching 10^-9 to 10^-12 and -g_tt approaching zero at the neutron star surface — is computed in this paper from Eqs. (14)-(19), so the self-citation is not the evidence for the metric behavior. The circularity is confined to framing/classification: the solution-branch endpoint is labeled with the authors' earlier terminology rather than independently derived. This is minor and does not reduce the core computation to its inputs.

full rationale

The derivation chain is: choose Bardeen/Hayward NED Lagrangians (Eqs. 2-3), static metric (Eq. 7), monopole ansatz (Eq. 12); integrate the modified TOV system (Eqs. 18-19) with external EOS (BSk19, SLy4, AP4) for fixed s or sq^2. The parameters s, q, and the EOS are inputs; q_c emerges as the boundary of the static solution space, and the frozen-state signature (metric degeneracy at the surface) is read off the computed metric. No parameter is fitted to data and then renamed a prediction; the central result — termination of the solution branch with 1/g_rr and -g_tt approaching zero at the surface — is an output of the stated equations, not a tautology. The acknowledged non-circular gaps are physical: the paper itself states in Sec. IV that 'a rigorous analysis of the stability of these frozen neutron stars under radial perturbations is essential,' and the endpoint is approached only to 10^-9 to 10^-12, with U^mu=(e^{-alpha},0,0,0) becoming null and the 1/(2e^{-2beta}) factor in Eq. (19) singular in the exact limit. These are correctness/interpretation concerns about whether the limit is a genuine equilibrium, not circularity. The only circularity-adjacent element is the import of the 'frozen state / critical horizon' concept from the authors' own boson-star papers [19-23]; since that citation supplies the name and interpretation rather than the numerical proof, it is a minor framing self-citation, not a load-bearing reduction. Verdict: no significant circularity; score 2.

Axiom & Free-Parameter Ledger

2 free parameters · 5 axioms · 0 invented entities

No new entities are introduced; the nonlinear magnetic monopole is a pre-existing hypothetical particle whose gravitational field is modeled via known NED Lagrangians. The ledger instead captures the hand-picked coupling s, the chosen central densities, the assumed validity of the NED Lagrangians and EOS extrapolations, and the interpretative step from numerical boundary to physical frozen state.

free parameters (2)
  • NED coupling s = 4.75e-81 m^3 J^-1 Wb^-2 (Bardeen), 7e-81 m^3 J^-1 Wb^-2 (Hayward)
    Chosen by hand, carried over from prior boson star studies (Refs [19-21]); the existence and value of q_c depend on this choice and it is not constrained by any observation.
  • central density rho_c = 0.5, 1.0, 1.5 x 10^18 kg/m^3
    Boundary condition scanned in Tables I-II; q_c varies with rho_c. This is a standard input in neutron star modeling, not fitted to a specific object.
axioms (5)
  • domain assumption Bardeen and Hayward Lagrangians (Eqs. 2 and 3) are the correct effective description of a nonlinear magnetic monopole's electromagnetic self-interaction.
    Adopted from Refs [15-18]; no experimental evidence for monopoles or these specific forms. The central claim is entirely within these models.
  • domain assumption Static, spherically symmetric metric and perfect-fluid matter (Eqs. 7-8).
    Standard TOV idealization; no rotation, anisotropy, or time dependence is considered.
  • domain assumption Magnetic monopole ansatz A = q cos(theta) dphi (Eq. 12).
    Restricts to a single monopole charge with a purely radial magnetic field; assumes the NED field is sourced only by the charge.
  • domain assumption Nuclear EOSs (BSk19, SLy4, AP4) remain valid up to the extreme pressures and densities reached near q_c, including the dense surface layer.
    These EOSs are fitted to nuclear matter near saturation; extrapolation to the frozen-state regime is untested. The paper invokes a causality pressure limit for AP4.
  • ad hoc to paper The numerical failure to find solutions beyond q_c marks a physical transition to a frozen state, not a solver artifact.
    Sec. III: 'beyond this critical magnetic charge q_c, it becomes numerically infeasible to obtain physically meaningful solutions.' This is an interpretation rather than a theorem.

reviewed 2026-08-04 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Frozen Neutron Stars." pith.science (2026). https://pith.science/paper/MTXZDVTV

@misc{pith2026250909338,
  author       = {Pith},
  title        = {Pith review of: Frozen Neutron Stars},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MTXZDVTV}},
  note         = {Machine review of arXiv:2509.09338}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

We investigate neutron stars with nonlinear magnetic monopoles in the framework of the Einstein-nonlinear electrodynamics model, specifically within the Bardeen and Hayward models. Solving the modified Tolman-Oppenheimer-Volkoff equations for three different equations of state, we find that upon reaching the critical magnetic charge $q_{c}$, neutron stars enter frozen states characterized by the critical horizon. This extends the concept of frozen states to compact objects composed of ordinary matter (non-field matter), thereby offering a new perspective for related research.

Figures

Figures reproduced from arXiv: 2509.09338 by Chen Tan, Yong-Qiang Wang.

Figure 1
Figure 1. Figure 1: FIG. 1: The top subplots show the radial pressure profiles in the Einstein-Bardeen framework with [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: The top upper subplots display the variation of the compactness [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: The top subplots show [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: The top subplots show [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: The left subplots shows the dependence of the [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: The [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗

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

35 extracted references · 16 linked inside Pith

  1. [1]

    filled hard candy like

    The total gravitational mass of the compact star mat- terM M≡m(R) = 4π Z R 0 ρ(x)x2dx.(21) The total ADM mass, incorporating contributions from the nonlinear magnetic monopoles, can be derived from the 1 r term in the asymptotic expansion of the metric at spatial infinityr→ ∞. Coincidentally, at spatial infinity r→ ∞, Einstein-Bardeen and Einstein-Hayward...

  2. [2]

    Equa- tions of state for supernovae and compact stars,

    M. Oertel, M. Hempel, T. Kl¨ ahn and S. Typel, “Equa- tions of state for supernovae and compact stars,” Rev. Mod. Phys.89, no.1, 015007 (2017) [arXiv:1610.03361 [astro-ph.HE]]

  3. [3]

    Gravitational Waves and Gamma- rays from a Binary Neutron Star Merger: GW170817 and GRB 170817A,

    B. P. Abbottet al.[LIGO Scientific, Virgo, Fermi-GBM and INTEGRAL], “Gravitational Waves and Gamma- rays from a Binary Neutron Star Merger: GW170817 and GRB 170817A,” Astrophys. J. Lett.848, no.2, L13 (2017) [arXiv:1710.05834 [astro-ph.HE]]

  4. [4]

    Neutron stars as extreme labora- tories for gravity tests,

    L. Shao and K. Yagi, “Neutron stars as extreme labora- tories for gravity tests,” Sci. Bull.67, 1946-1949 (2022) [arXiv:2209.03351 [gr-qc]]

  5. [5]

    Formation of very strongly magnetized neutron stars - implications for gamma-ray bursts,

    R. C. Duncan and C. Thompson, “Formation of very strongly magnetized neutron stars - implications for gamma-ray bursts,” Astrophys. J. Lett.392, L9 (1992)

  6. [6]

    Magnetars: the physics behind observations. A review,

    R. Turolla, S. Zane and A. Watts, “Magnetars: the physics behind observations. A review,” Rept. Prog. Phys.78, no.11, 116901 (2015) [arXiv:1507.02924 [astro- ph.HE]]

  7. [7]

    Unifying the observational diversity of isolated neutron stars via magneto-thermal evolution models,

    D. Vigan` o, N. Rea, J. A. Pons, R. Perna, D. N. Aguilera and J. A. Miralles, “Unifying the observational diversity of isolated neutron stars via magneto-thermal evolution models,” Mon. Not. Roy. Astron. Soc.434, 123 (2013) [arXiv:1306.2156 [astro-ph.SR]]

  8. [8]

    The Soft gamma re- peaters as very strongly magnetized neutron stars. 2. Quiescent neutrino, x-ray, and Alfven wave emission,

    C. Thompson and R. C. Duncan, “The Soft gamma re- peaters as very strongly magnetized neutron stars. 2. Quiescent neutrino, x-ray, and Alfven wave emission,” Astrophys. J.473, 322 (1996)

  9. [9]

    The Soft gamma re- peaters as very strongly magnetized neutron stars - 1. Radiative mechanism for outbursts,

    C. Thompson and R. C. Duncan, “The Soft gamma re- peaters as very strongly magnetized neutron stars - 1. Radiative mechanism for outbursts,” Mon. Not. Roy. As- tron. Soc.275, 255-300 (1995)

  10. [10]

    Early evolu- tion of newly born magnetars with a strong toroidal field,

    S. Dall’Osso, S. N. Shore and L. Stella, “Early evolu- tion of newly born magnetars with a strong toroidal field,” Mon. Not. Roy. Astron. Soc.398, 1869 (2009) [arXiv:0811.4311 [astro-ph]]

  11. [11]

    Magnetic Monopoles in Unified Gauge The- ories,

    G. ’t Hooft, “Magnetic Monopoles in Unified Gauge The- ories,” Nucl. Phys. B79, 276-284 (1974)

  12. [12]

    Creation of magnetic monopoles in pulsars,

    M. Bonnardeau, A. K. Drukier, “Creation of magnetic monopoles in pulsars,” Nature277, 543-544 (1979)

  13. [13]

    MONOPOLES IN NEUTRON STARS,

    J. A. Harvey, “MONOPOLES IN NEUTRON STARS,” Nucl. Phys. B236, 255-268 (1984)

  14. [14]

    Superheavy Magnetic Monopoles and Proton Decay,

    V. A. Rubakov, “Superheavy Magnetic Monopoles and Proton Decay,”

  15. [15]

    Monopole Catalysis of Baryon Decay,

    C. G. Callan, Jr., “Monopole Catalysis of Baryon Decay,” Nucl. Phys. B212, 391-400 (1983)

  16. [16]

    Bardeen, in Proceedings of GR5, Tiflis, U.S.S.R

    J. Bardeen, in Proceedings of GR5, Tiflis, U.S.S.R. (1968)

  17. [17]

    Formation and evaporation of regu- lar black holes,

    S. A. Hayward, “Formation and evaporation of regu- lar black holes,” Phys. Rev. Lett.96, 031103 (2006) [arXiv:gr-qc/0506126 [gr-qc]]

  18. [18]

    Regular black hole in general relativity coupled to nonlinear electrodynam- ics,

    E. Ayon-Beato and A. Garcia, “Regular black hole in general relativity coupled to nonlinear electrodynam- ics,” Phys. Rev. Lett.80, 5056-5059 (1998) [arXiv:gr- qc/9911046 [gr-qc]]

  19. [19]

    The Bardeen model as a nonlinear magnetic monopole,

    E. Ayon-Beato and A. Garcia, “The Bardeen model as a nonlinear magnetic monopole,” Phys. Lett. B493, 149- 152 (2000) [arXiv:gr-qc/0009077 [gr-qc]]

  20. [20]

    Non-Topological Soliton Bardeen Boson Stars and Their Frozen States,

    Z. H. Zhao, Y. N. Gu, S. C. Liu, L. X. Huang and Y. Q. Wang, “Non-Topological Soliton Bardeen Boson Stars and Their Frozen States,” Eur. Phys. J. C85, 511 (2025) [arXiv:2502.14153 [gr-qc]]

  21. [21]

    From Bardeen-boson stars to black holes without event horizon,

    X. E. Wang, “From Bardeen-boson stars to black holes without event horizon,” [arXiv:2305.19057 [gr-qc]]

  22. [22]

    Frozen Hayward-boson stars,

    Y. Yue and Y. Q. Wang, “Frozen Hayward-boson stars,” JCAP05, 066 (2025) [arXiv:2312.07224 [gr-qc]]

  23. [23]

    Hayward Boson Stars,

    S. S. Chicaiza-Medina and J. C. Degollado, “Hayward Boson Stars,” [arXiv:2508.11906 [gr-qc]]

  24. [24]

    Frozen states of charged boson stars,

    Y. Brihaye and B. Hartmann, “Frozen states of charged boson stars,” [arXiv:2507.08946 [gr-qc]]

  25. [25]

    On Continued grav- itational contraction,

    J. R. Oppenheimer and H. Snyder, “On Continued grav- itational contraction,” Phys. Rev.56, 455-459 (1939)

  26. [26]

    Y. B. Zel’dovich and I. D. Novikov,Relativistic Astro- physics 1: Stars and Relativity, (University of Chicago Press, Chicago 1971), p. 369 (translation from the 1967 Russian edition)

  27. [27]

    Introducing the black hole,

    R. Ruffini and J. A. Wheeler, “Introducing the black hole,” Phys. Today24, no.1, 30 (1971)

  28. [28]

    Static solutions of Einstein’s field equa- tions for spheres of fluid,

    R. C. Tolman, “Static solutions of Einstein’s field equa- tions for spheres of fluid,” Phys. Rev.55, 364-373 (1939)

  29. [29]

    On massive neu- tron cores,

    J. R. Oppenheimer and G. M. Volkoff, “On massive neu- tron cores,” Phys. Rev.55, 374-381 (1939)

  30. [30]

    Analytical representations of uni- fied equations of state for neutron-star matter,

    A. Y. Potekhin, A. F. Fantina, N. Chamel, J. M. Pear- son and S. Goriely, “Analytical representations of uni- fied equations of state for neutron-star matter,” Astron. Astrophys.560, A48 (2013) [arXiv:1310.0049 [astro- ph.SR]]

  31. [31]

    A unified equation of state of dense matter and neutron star structure,

    F. Douchin and P. Haensel, “A unified equation of state of dense matter and neutron star structure,” Astron. As- trophys.380, 151 (2001) [arXiv:astro-ph/0111092 [astro- ph]]

  32. [32]

    Spin - isospin structure and pion condensation in nucleon mat- ter,

    A. Akmal and V. R. Pandharipande, “Spin - isospin structure and pion condensation in nucleon mat- ter,” Phys. Rev. C56, 2261-2279 (1997) [arXiv:nucl- th/9705013 [nucl-th]]

  33. [33]

    Analytical Representation for Equations of State of Dense Matter,

    C. Gungor and K. Y. Eksi, “Analytical Representation for Equations of State of Dense Matter,”

  34. [34]

    Frozen gravitational stars in Einsteinian cubic gravity,

    Y. Q. Wang, “Frozen gravitational stars in Einsteinian cubic gravity,”

  35. [35]

    Boson stars and their frozen states in an infinite tower of higher- derivative gravity,

    T. X. Ma, T. F. Fang and Y. Q. Wang, “Boson stars and their frozen states in an infinite tower of higher- derivative gravity,” Eur. Phys. J. C85, no.5, 542 (2025) [arXiv:2406.08813 [gr-qc]]

This paper was first reviewed by deepseek-v4-flash on August 4, 2026.