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 →
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 →
Frozen Neutron Stars
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
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.
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
- 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.
Referee Report
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)
- [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.
- [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.
- [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.
- [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)
- [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.
- [Introduction] Typo: 'gracitational' should be 'gravitational'.
- [Sec. III C] Typo: 'ttherefore' should be 'therefore'.
- [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.
- [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
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
-
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
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)
- central density rho_c =
0.5, 1.0, 1.5 x 10^18 kg/m^3
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.
- domain assumption Static, spherically symmetric metric and perfect-fluid matter (Eqs. 7-8).
- domain assumption Magnetic monopole ansatz A = q cos(theta) dphi (Eq. 12).
- 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.
- 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.
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}
}
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
Reference graph
Works this paper leans on
-
[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...
2010
-
[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]]
Pith/arXiv arXiv 2017
-
[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]]
Pith/arXiv arXiv 2017
-
[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]]
Pith/arXiv arXiv 1946
-
[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)
1992
-
[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]]
Pith/arXiv arXiv 2015
-
[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]]
Pith/arXiv arXiv 2013
-
[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)
1996
-
[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)
1995
-
[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]]
Pith/arXiv arXiv 2009
-
[11]
Magnetic Monopoles in Unified Gauge The- ories,
G. ’t Hooft, “Magnetic Monopoles in Unified Gauge The- ories,” Nucl. Phys. B79, 276-284 (1974)
1974
-
[12]
Creation of magnetic monopoles in pulsars,
M. Bonnardeau, A. K. Drukier, “Creation of magnetic monopoles in pulsars,” Nature277, 543-544 (1979)
1979
-
[13]
MONOPOLES IN NEUTRON STARS,
J. A. Harvey, “MONOPOLES IN NEUTRON STARS,” Nucl. Phys. B236, 255-268 (1984)
1984
-
[14]
Superheavy Magnetic Monopoles and Proton Decay,
V. A. Rubakov, “Superheavy Magnetic Monopoles and Proton Decay,”
-
[15]
Monopole Catalysis of Baryon Decay,
C. G. Callan, Jr., “Monopole Catalysis of Baryon Decay,” Nucl. Phys. B212, 391-400 (1983)
1983
-
[16]
Bardeen, in Proceedings of GR5, Tiflis, U.S.S.R
J. Bardeen, in Proceedings of GR5, Tiflis, U.S.S.R. (1968)
1968
-
[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]]
Pith/arXiv arXiv 2006
-
[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]]
arXiv 1998
-
[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]]
Pith/arXiv arXiv 2000
-
[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]]
Pith/arXiv arXiv 2025
-
[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]
Y. Yue and Y. Q. Wang, “Frozen Hayward-boson stars,” JCAP05, 066 (2025) [arXiv:2312.07224 [gr-qc]]
Pith/arXiv arXiv 2025
-
[23]
S. S. Chicaiza-Medina and J. C. Degollado, “Hayward Boson Stars,” [arXiv:2508.11906 [gr-qc]]
-
[24]
Frozen states of charged boson stars,
Y. Brihaye and B. Hartmann, “Frozen states of charged boson stars,” [arXiv:2507.08946 [gr-qc]]
-
[25]
On Continued grav- itational contraction,
J. R. Oppenheimer and H. Snyder, “On Continued grav- itational contraction,” Phys. Rev.56, 455-459 (1939)
1939
-
[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)
1971
-
[27]
Introducing the black hole,
R. Ruffini and J. A. Wheeler, “Introducing the black hole,” Phys. Today24, no.1, 30 (1971)
1971
-
[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)
1939
-
[29]
On massive neu- tron cores,
J. R. Oppenheimer and G. M. Volkoff, “On massive neu- tron cores,” Phys. Rev.55, 374-381 (1939)
1939
-
[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]]
Pith/arXiv arXiv 2013
-
[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]]
Pith/arXiv arXiv 2001
-
[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]]
arXiv 1997
-
[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]
Frozen gravitational stars in Einsteinian cubic gravity,
Y. Q. Wang, “Frozen gravitational stars in Einsteinian cubic gravity,”
-
[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]]
Pith/arXiv arXiv 2025
This paper was first reviewed by deepseek-v4-flash on August 4, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.