REVIEW 4 major objections 4 minor 2 cited by
Frozen states of charged boson stars
T0 review · 4 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Charged boson stars can freeze into horizonless black-hole mimics with a de Sitter core and a black-hole exterior.
desk verdict The Horndeski extension is genuinely new and the light-ring analysis is clean, but the paper's central claim that exact frozen states exist is not actually demonstrated — the numerical branches stop at a double-zero limit that looks like an extremal horizon. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The key machinery is the approach of the metric function $N(r) = 1 - 2m(r)/r$ to a double zero at a finite radius, read as a quasi-horizon: the scalar and electric fields stay smooth there, while the energy-momentum tensor inside approaches the de Sitter form $\rho = -p_r = -p_t$ and outside approaches the Reissner-Nordström form. Two ingredients push the branches to this limit: the bounded, exponential scalar self-interaction $U(|\Psi|) = \mu^2\eta^2(1-\exp(-|\Psi|^2/\eta^2))$ at intermediate gauge coupling $q$, and, for the simple mass potential, a negative Horndeski vector-tensor coupling $\gamma$. The free parameters are the central value $C = \psi(0)$, the gauge coupling $q$, and the algebraic coupling $\alpha$; the solutions are constructed numerically by shooting.
What would settle it
Recompute the branch-B family with higher grid resolution and Richardson extrapolation to see whether the minimum of $N(r)$ converges to a positive value or to zero. If the minimum saturates above zero, or if curvature invariants diverge at the shell as the minimum tends to zero, then the frozen state is a numerical artifact or an extremal-horizon limit rather than a regular horizonless solution.
Extended reading notes
Core claim
The discovery is a branch of numerically constructed charged boson stars that terminates in a frozen state. Along branch B, decreasing the scalar field at the origin (or lowering the frequency $\omega$ in the Horndeski case) drives the metric function $N(r)$ toward a double zero at a finite radius $r_c$; the interior becomes de Sitter-like with $\rho = -p_r = -p_t$, a shell of finite thickness carries the transition, and the exterior is Reissner-Nordström. These are not black holes: no horizon and no singularity form. In standard Einstein gravity with a mass potential alone the minimum of $N$ never reaches zero, so the exponential self-interaction is the enabling ingredient; with negative Horndeski coupling the self-interaction can be dropped. The frozen-state configurations have $M/Q_N < 1$ in the Horndeski case and possess both a stable inner and an unstable outer light ring.
Load-bearing premise
The central claim depends on the numerical solutions staying accurate as the metric function's minimum approaches zero, and on that limiting dip being a genuinely smooth, horizonless shell rather than a numerical breakdown or an extremal horizon.
Editorial extensions
If this is right
- The frozen states are horizonless, globally regular alternatives to Reissner-Nordström black holes with the same exterior, so from afar their lensing and light-ring structure mimic a black hole.
- Because each has an inner stable and an outer unstable light ring, the known instability channel for ultracompact objects applies, so the states are not automatically stable as static solutions.
- In the Horndeski case with $\gamma < 0$, all solutions satisfy $M/Q_N < 1$, so they are stable against decay into $Q_N$ free scalar bosons.
- The exponential self-interaction and the Horndeski term play interchangeable roles in producing the frozen-state limit: either one can drive the quasi-horizon, whereas neither standard electrodynamics with a mass potential alone nor the ungauged model does.
Reading between the lines
- Editorial inference: If the limiting double zero is genuinely regular, the frozen-state branches provide a one-parameter family connecting ordinary boson stars to gravastar-like objects, which could be used to model gravitational-wave echoes.
- Editorial inference: The stable light ring inside the shell, combined with known light-ring instability results, suggests these frozen states may be nonlinearly unstable; evolving the full system numerically would test whether they are transient or long-lived.
- Editorial inference: A sharper check would be to verify whether curvature invariants at the shell stay finite as the minimum of $N(r)$ tends to zero; if they diverge, the frozen state is better read as an extremal-horizon limit than as a horizonless object.
- Editorial inference: The same construction should also work for other asymptotically flat charged solitons, such as gauged Q-balls, whenever charge repulsion balances gravity at intermediate coupling.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies spherically symmetric, static, charged boson stars in a U(1)-gauged scalar field model minimally coupled to gravity, with an optional Horndeski vector-tensor coupling. It claims that for intermediate gauge couplings, the self-interacting scalar potential (2.4) allows branches of solutions to approach configurations where the metric function N(r) develops a double zero at a finite radius, which the authors identify as 'frozen states'—globally regular, horizonless objects with a de Sitter interior, a Reissner-Nordström exterior, and a thin shell replacing the event horizon. For the Horndeski model with negative γ, the same behavior is claimed to occur even without scalar self-interaction. The paper also computes light-ring effective potentials and finds one stable and one unstable light ring for the near-frozen configurations.
Significance. If established, the existence of globally regular, horizonless ultracompact objects with a de Sitter core and a black-hole exterior would be a concrete field-theoretic realization of the Mazur-Mottola gravitational condensate star, and the claim that ordinary (non-linear) electrodynamics suffices in a gauged scalar model is of considerable interest. The paper also makes a falsifiable prediction about light-ring pairs and notes the possible light-ring instability, connecting to active literature on black-hole mimickers. However, the central existence claim is currently supported only by sparse numerical shooting results for a few parameter values and no convergence or error analysis, so the significance is contingent on the limiting procedure being made rigorous.
major comments (4)
- [Section 3.1, Fig. 3] The central claim that branch B terminates in a frozen state is based on the observation that the minimum of N(r) approaches zero as C decreases toward C≈11.2. However, a double zero of N(r) at r=rc makes rc a degenerate horizon: the proper radial distance ∫ dr/√N diverges logarithmically, and the coordinate r=rc is not part of the manifold in the standard sense. The paper does not demonstrate that the limiting configuration is a regular horizonless shell, and Eq. (2.9) contains terms proportional to 1/N^2 that become singular in this limit. No convergence tests, error bars, or an exact limiting solution with N(rc)=N'(rc)=0 are provided. Since the abstract and introduction explicitly claim that the shell 'replaces the event horizon', this is the load-bearing step of the paper, and it needs either a rigorous construction of the limiting solution or a careful numerical convergence analysis showing that the limit is not an extremal horizon.
- [Section 3.2, Fig. 5] The same degenerate-horizon concern applies to the γ<0 case. The branch is said to end exactly when the minimum of N(r) approaches zero, but only three values of ω are shown (ω=0.04, 0.03, 0.02), with no evidence that the sequence converges to a finite shell rather than to a singular or extremal configuration. The equations for γ<0 differ from the γ=0 case, and the scalar-field equation (2.9) still contains 1/N^2 terms; the regularity of the scalar and gauge fields at the would-be zero of N must be checked explicitly. Please provide a convergence study for these solutions and clarify the nature of the limiting configuration.
- [Introduction and Abstract] The manuscript states in the abstract that frozen states are globally regular and have a thin shell that 'replaces the event horizon'. The numerical evidence presented in Fig. 3 and Fig. 5 is for finite-C solutions with N_min>0, which are quasi-horizon ultracompact objects but not yet frozen states. The text itself uses the phrases 'approaches a double zero' and 'strongly suggest', indicating that the limiting configuration was not actually constructed. The claim that the event horizon is replaced by a shell of finite thickness is therefore an extrapolation, and the manuscript should either provide the limiting solution or temper the abstract and conclusions to what is demonstrated.
- [Section 3.1, parameter range] The demonstration of the crucial role of self-interaction relies on only two values of the gauge coupling (q=0.005 and q=0.01) at a single value of α=0.0001. The text states that the frozen state appears for 'intermediate values' of q, but no scan over q is shown for the self-interacting potential. Given that the claimed phenomenon depends sensitively on q (the branches differ qualitatively between q=0, 0.005, and 0.01), a more systematic parameter study is needed to establish the robustness of the effect.
minor comments (4)
- [Fig. 2 caption] The caption contains a typo: 'The show the value C' should be 'We show the value C'.
- [Section 4, Eq. (4.19)] The light-ring analysis is performed for finite-C solutions (Fig. 7), not for the limiting frozen state. Since the limiting configuration is not explicitly constructed, the statement that 'the frozen states possess one stable and one unstable lightring' should be phrased as a property of the near-frozen solutions, or the limiting analysis should be provided.
- [Section 2, Eq. (2.5)] The notation is mostly clear, but the boundary condition N(0)=0 in Eq. (2.10) is the standard regularity condition at the origin; it would help to state explicitly that this does not imply a horizon because the radial coordinate is not the areal radius near r=0. This is a minor clarification, not a technical error.
- [Section 3.2] The text says 'vector-tensor boson stars do exist in the limit C → 0 which corresponds to Ω → 0', but in Section 3.1 the limit Ω→0 is associated with the approach to the frozen state. The correspondence between C→0 and Ω→0 for γ<0 should be stated more carefully, since it seems to describe a different regime.
Circularity Check
No load-bearing circularity; the central derivation is self-contained, with only non-essential self-citations for context.
full rationale
The central results are obtained by numerically integrating the coupled ODEs (2.6)-(2.9) under the boundary conditions (2.10)-(2.11); no quantity is fitted to the claimed frozen-state output. The scalar-field central value C fixes omega, and branches are mapped by varying C or omega; the approach of the minimum of N(r) to zero is a computed outcome, not an imposed condition. The lightring effective potential (4.19) is derived directly from the obtained metric, so the lightring claims are genuine outputs. The comparisons between the mass potential, the self-interacting potential (2.4), and the Horndeski term (2.3) are also numerical outcomes of the same equations, not consequences of a presupposed conclusion. The paper's reliance on the authors' prior work [20,21,29] is contextual (previously discussed thin-shell solutions, the domain boundary q>0.75) and does not carry the existence argument. The main caveat is numerical rather than circular: the limiting configuration with a double zero of N is not constructed exactly, and no convergence tests are reported; the text itself says the metric function 'approaches' a double zero and that results 'strongly suggest' the limiting behavior. This is a verification/rigor limitation, not a reduction of the derivation to its inputs, and therefore does not raise the circularity score beyond the minor self-referential context.
Assumptions & free parameters
free parameters (4)
- alpha = 4 pi G eta^2 =
0.0001
- q, U(1) gauge coupling =
0.005 / 0.01 for gamma=0; 0.7 for gamma<0
- gamma, Horndeski vector-tensor coupling =
-3 and -5 for frozen states; +3 and +5 for comparison
- C = psi(0), central scalar field value =
11.2, 12, 13, 15 for gamma=0; psi(0)=1 for gamma!=0 scans
assumptions (5)
- domain assumption Spherically symmetric static ansatz (2.5) is sufficient to capture frozen states.
- domain assumption Boundary conditions (2.10) and asymptotic flatness (2.11) define the regular localized solution space.
- standard math The scalar field decays exponentially so that omega - q V_infinity < 1 holds.
- ad hoc to paper Numerical shooting solutions converge to true solutions of the boundary value problem.
- standard math Lightrings are identified from the effective potential Veff = L_z^2 N sigma^2 / r^2 for equatorial null geodesics.
Cite this review
Pith. "Pith review of Frozen states of charged boson stars." pith.science (2026). https://pith.science/paper/VFL5D4WK
@misc{pith2026250708946,
author = {Pith},
title = {Pith review of: Frozen states of charged boson stars},
year = {2026},
howpublished = {\url{https://pith.science/paper/VFL5D4WK}},
note = {Machine review of arXiv:2507.08946}
}
read the original abstract
In this paper, we study frozen states of charged boson stars. These solutions are globally regular and exist in a U(1) gauged scalar field model minimally coupled to gravity for suitable choices of the coupling constants. These configurations are field theoretical realizations of the Mazur-Mottola solution with a de Sitter interior, a black hole exterior and a thin shell that interpolates between the two and replaces the event horizon. We demonstrate that standard electrodynamics is sufficient to find these frozen states, but that the self-interaction of the scalar field is crucial. Adding Horndeski vector-tensor gravity to the model allows the frozen states to exist without self-interaction though. The frozen states possess one stable and one unstable lightring, the former inside the thin shell, the latter in the black hole exterior.
Figures
Figures from the paper (4 more)
Forward citations
Cited by 2 Pith papers
-
Frozen Neutron Stars
Solving modified Tolman-Oppenheimer-Volkoff equations with Bardeen and Hayward nonlinear electrodynamics, this paper finds that neutron stars reach 'frozen states' with a critical horizon at a critical magnetic charge.
-
Scalarization of Bardeen spacetime
For scalarization of the full Bardeen spacetime, small magnetic charges give the usual smooth scalarization threshold, while large charges end in a 'frozen' horizonless scalarized state rather than a Bardeen black hole.
Reference graph
Works this paper leans on
-
[21]
Y. Brihaye and B. Hartmann: Boson stars and black holes with wavy scalar hair, Phys. Rev. D 105 (2022) no.10, 104063
work page 2022
-
[1]
Penrose: Gravitational collapse: The role of general relativity, Riv
R. Penrose: Gravitational collapse: The role of general relativity, Riv. Nuovo Cim. 1 (1969), 252. 8 0.0000001 0.0000010 0.0000100 0.0001000 0.0010000 0.0100000 0.1000000 0 100 200 300 400 500 Veff r C=15.0 C=13.0 C=12.0 C=11.2 0.000001 0.000010 0.000100 0.001000 0 5 10 15 20 25 30 35 40 Veff r ω = 0.04 ω = 0.03 ω = 0.02 Figure 7: Left: We show the effect...
work page 1969
- [2]
-
[3]
D. J. Kaup: Klein-Gordon Geon, Phys. Rev. 172 (1968), 1331
work page 1968
-
[4]
F. E. Schunck and E. W. Mielke: Boson stars: Rotation, formation, and evolution, Gen. Rel. Grav. 31 (1999), 787
work page 1999
-
[5]
R. Friedberg, T. D. Lee and Y. Pang: MINI - SOLITON STARS, Phys. Rev. D 35 (1987), 3640
work page 1987
- [6]
-
[7]
F. E. Schunck and E. W. Mielke: General relativistic boson stars, Class. Quant. Grav. 20 (2003), R301-R356
work page 2003
Show all 32 references
-
[8]
S. L. Liebling and C. Palenzuela: Dynamical boson stars, Living Rev. Rel. 26 (2023) no.1, 1
2023
-
[9]
J. Bardeen: Nonsingular general relativistic gravitational collapse, in Proceedings of the 5th International Conference on Gravitation and the Theory of Relativity, Tbilissi University Press, 1968
1968
-
[10]
Ayon-Beato and A
E. Ayon-Beato and A. Garcia: Regular black hole in general relativity coupled to nonlinear electrody- namics, Phys. Rev. Lett. 80 (1998), 5056-5059
1998
-
[11]
Ayon-Beato and A
E. Ayon-Beato and A. Garcia: The Bardeen model as a nonlinear magnetic monopole, Phys. Lett. B 493 (2000), 149-152
2000
-
[12]
P. O. Mazur and E. Mottola: Gravitational Condensate Stars: An Alternative to Black Holes, Universe 9 (2023) no.2, 88
2023
-
[13]
Brustein, A
R. Brustein, A. J. M. Medved and T. Simhon: Black holes as frozen stars, Phys. Rev. D 105 (2022) no.2, 024019
2022
-
[14]
Brustein, A
R. Brustein, A. J. M. Medved, T. Shindelman and T. Simhon: Black Holes as Frozen Stars: Regular Interior Geometry, Fortsch. Phys. 72 (2024) no.1, 2300188 9
2024
-
[15]
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. C 85 (2025), 511
2025
-
[16]
X. E. Wang: From Bardeen-boson stars to black holes without event horizon, [arXiv:2305.19057 [gr-qc]]
-
[17]
Yue and Y
Y. Yue and Y. Q. Wang: Frozen Hayward-boson stars, JCAP 05 (2025), 066
2025
-
[18]
C. A. Herdeiro and E. Radu: Spherical electro-vacuum black holes with resonant, scalar Q-hair, Eur. Phys. J. C 80 (2020) no.5, 390
2020
-
[19]
J. P. Hong, M. Suzuki and M. Yamada: Spherically Symmetric Scalar Hair for Charged Black Holes, Phys. Rev. Lett. 125 (2020) no.11, 111104
2020
-
[20]
Brihaye and B
Y. Brihaye and B. Hartmann: Strong gravity effects of charged Q-clouds and inflating black holes, Class. Quant. Grav. 38 (2021) no.6, 06LT01
2021
-
[22]
G. W. Horndeski: Second-order scalar-tensor field equations in a four-dimensional space, Int. J. Theor. Phys. 10 (1974), 363
1974
-
[23]
G. W. Horndeski: Conservation of Charge and the Einstein-Maxwell Field Equations, J. Math. Phys. 17 (1976) 1980
1976
-
[24]
Mueller-Hoissen and R
F. Mueller-Hoissen and R. Sippel: Spherically Symmetric Solutions of the Nonminimally Coupled Einstein-maxwell Equations, Class. Quant. Grav. 5 (1988), 1473
1988
-
[25]
Verbin: Magnetic and electric black holes in the vector-tensor Horndeski theory, Phys
Y. Verbin: Magnetic and electric black holes in the vector-tensor Horndeski theory, Phys. Rev. D 106 (2022) no.2, 2
2022
-
[26]
E. J. Copeland and M. I. Tsumagari: Q-balls in flat potentials, Phys. Rev. D 80 (2009), 025016
2009
-
[27]
Pugliese, H
D. Pugliese, H. Quevedo, J. A. Rueda H. and R. Ruffini: On charged boson stars, Phys. Rev. D 88 (2013), 024053
2013
-
[28]
J. D. L´ opez and M. Alcubierre: Charged boson stars revisited, Gen. Rel. Grav. 55 (2023) no.6, 72
2023
-
[29]
Brihaye and B
Y. Brihaye and B. Hartmann: Note on super-critical charged boson stars, Gen. Rel. Grav. 56 (2024) no.1, 6
2024
-
[30]
P. V. P. Cunha, E. Berti and C. A. R. Herdeiro: Light-Ring Stability for Ultracompact Objects, Phys. Rev. Lett. 119 (2017) no.25, 251102
2017
-
[31]
Cardoso, L
V. Cardoso, L. C. B. Crispino, C. F. B. Macedo, H. Okawa and P. Pani: Light rings as observational evidence for event horizons: long-lived modes, ergoregions and nonlinear instabilities of ultracompact objects, Phys. Rev. D 90 (2014) no.4, 044069
2014
-
[32]
P. V. P. Cunha, C. Herdeiro, E. Radu and N. Sanchis-Gual: Exotic Compact Objects and the Fate of the Light-Ring Instability, Phys. Rev. Lett. 130 (2023) no.6, 061401 10
2023
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.