REVIEW 1 major objections 6 minor 1 cited by
Multi-frequency Proca stars form continuous families of solutions, and some of them are linearly stable.
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 · grok-4.5
2026-07-15 14:23 UTC pith:VZL2Z73X
load-bearing objection Solid continuum maps and stability bands for multi-frequency Proca stars; the advance is real but incremental, and the “general” stability claim still rests on a limited J-scan plus deferred nonlinear runs. the 1 major comments →
The continuum spectrum of nonrelativistic multi-frequency Proca 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
At fixed particle number, spherical multi-frequency Proca stars form continuous one- and two-parameter families that connect the discrete stationary states of constant linear polarization; within the families that include a nodeless component, finite stability bands exist against general linear perturbations, so that a non-negligible nodeless component is necessary (but not sufficient) for stability and radial stability alone is not enough for full linear stability.
What carries the argument
The spherical multi-frequency ansatz for the s=1 Schrödinger-Poisson system, together with the associated linear eigenvalue problem for complex growth rates λ of angular-momentum multipoles J, which partitions each continuous family into stable and unstable regions.
Load-bearing premise
That the absence of unstable modes for angular momenta up to J=6, plus consistency with a few forthcoming dynamical runs, is enough to claim linear stability against all perturbations.
What would settle it
A nonlinear numerical evolution of a configuration inside a reported stability band that develops a growing mode, or the discovery of an unstable eigenvalue with positive real part for some J>6 that is absent from the present spectrum.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper constructs and classifies spherical multi-frequency solutions of the nonrelativistic s=1 Schrödinger-Poisson system. Using the multi-component ansatz (3) and the nonlinear eigenvalue problem (4)/(6), it shows that, at fixed particle number N, 2- and 3-component configurations form continuous one- and two-parameter families that interpolate between discrete stationary states of constant linear polarization (Secs. III B–C, Figs. 2–4, Table I). Linear stability is analyzed via the perturbation system (16) for angular momenta J ≤ 6. Only families that include a nodeless component (n_x = 0) contain stable bands; a non-negligible nodeless component is necessary but not sufficient for stability, and pure radial (J=0) stability does not guarantee full linear stability (Sec. IV, Figs. 5–7). The authors briefly discuss possible observational signatures of particle spin in ultralight dark-matter models.
Significance. If the continuum construction and the reported stability bands hold, the work substantially enlarges the known solution space of self-gravitating spin-1 solitons beyond the discrete stationary states previously studied. The explicit maps of continuous families (Figs. 2, 4) and the identification of linearly stable multi-frequency configurations that coexist with the ground state are new and of direct interest for ultralight vector dark-matter phenomenology. The numerical methodology builds on the authors’ earlier validated discretization of the linearized system, and the necessity-but-not-sufficiency of a nodeless component is a clean, falsifiable structural result. The main limitation is that full linear stability against general perturbations is inferred from a finite J-scan together with deferred nonlinear simulations; within that scope the contribution is solid and timely.
major comments (1)
- Sec. IV B (and footnote 6) and Sec. II B: the claim of linear stability “against general perturbations” rests on the absence of eigenvalues with positive real part for J ≤ 6, plus the statement that “whenever numerical simulations were performed, the evolution was found to be consistent” with those simulations deferred to a forthcoming paper. Footnote 6 itself notes that higher-J instabilities cannot yet be ruled out. For the abstract and Table I statements to be fully supported, either (i) a clearer bound on the J-range that can be considered exhaustive for these configurations, or (ii) a short quantitative summary of the already-performed dynamical checks (even if the full study is forthcoming), should be supplied. Without one of these, the wording should be softened to “stable against perturbations with J ≤ 6 (and consistent with available dynamical runs).”
minor comments (6)
- Fig. 4 caption and surrounding text: the family labels are written (n_x, n_y, n_x) = (0,1,2) etc.; the third index should be n_z.
- Fig. 7 caption: “(0, n_x, n_y) families” should read “(0, n_y, n_z) families” for consistency with the notation of Sec. III C and Table I.
- Sec. III C, first paragraph: “0 < n_y < n_z” is stated after fixing n_x = 0; a brief reminder that the ordering n_x < n_y < n_z is maintained would avoid any ambiguity for readers who skip the earlier footnote.
- Eq. (24a): the sum index is written “∑_i” while the left-hand side is E_i; a distinct dummy index would improve readability.
- Table I: the parenthetical stability entries “(0)” for stationary and 1-component states are not defined in the caption or main text; a short legend would help.
- Reference [32] points to a movie whose URL is rendered as a string of boxes in the manuscript; the link should be checked for the final version.
Circularity Check
No by-construction reductions; continuum families and stability bands are independent numerical outputs built on a self-cited framework that is not load-bearing for the new claims.
full rationale
The paper’s central results—the continuous one- and two-parameter families of 2- and 3-component multi-frequency solutions at fixed N (Secs. III B–C, Figs. 2–4) and the associated linear-stability bands (Sec. IV, Figs. 5–7)—are obtained by direct numerical solution of the nonlinear eigenvalue problem (4)/(22) for continuously varied central amplitudes and of the linearized spectral problem (16) for J ≤ 6. These outputs do not reduce algebraically or statistically to any fitted input parameter, nor are they defined in terms of the quantities they claim to predict. The spherical multi-frequency ansatz (3), the Hamiltonian operator (7), the form of the perturbation matrices (18)/(25), and the known stability of the pure ground state are taken from the authors’ prior works [23,24]; that self-citation supplies the setup rather than the target continuum maps or the necessity-but-not-sufficiency of a nodeless component. No uniqueness theorem is imported to forbid alternatives, no parameter is fitted and then re-predicted, and no known empirical pattern is merely renamed. Once the Schrödinger-Poisson system and the spherical ansatz are granted, the derivation chain is self-contained numerical exploration. The incomplete J-scan and deferred nonlinear simulations are correctness caveats, not circularity.
Axiom & Free-Parameter Ledger
free parameters (1)
- total particle number N =
43.5 (arbitrary)
axioms (4)
- domain assumption The non-relativistic s=1 Schrödinger-Poisson system correctly describes the dynamics of a massive vector field in the weak-field, low-velocity regime.
- domain assumption Spherical multi-frequency solutions are completely captured by the Cartesian ansatz (3) with real radial profiles and constant frequencies.
- domain assumption Mode stability (all eigenvalues of the linearized operator purely imaginary) is a reliable indicator of dynamical stability for the configurations studied.
- standard math Standard spectral theory of self-adjoint and non-self-adjoint operators on radial function spaces applies to the discretized perturbation matrices.
invented entities (1)
-
multi-frequency Proca stars (2- and 3-component continuum families)
no independent evidence
read the original abstract
Multi-frequency Proca stars are excited selfgravitating solutions of the $s=1$ Schr\"odinger-Poisson system that generalize the conventional stationary states of a massive vector field. Unlike stationary states, which are characterized by a single oscillation frequency, multi-frequency configurations exhibit a quasi-periodic dynamics involving two or three distinct frequencies. In this paper, we present a systematic study of the spectrum of spherical multi-frequency Proca stars and show that, at fixed particle number, they form continuous families interpolating between discrete stationary states of constant linear polarization. Furthermore, we analyze their stability and demonstrate that a subset of these multi-frequency configurations are linearly stable against general perturbations. In particular, we show that a necessary, although not sufficient, condition for stability is the presence of a non-negligible nodeless component, and that radial stability alone is not sufficient to guarantee full linear stability. Finally, we briefly discuss the potential implications of multi-frequency states for proving the particle spin in ultralight dark matter models.
Forward citations
Cited by 1 Pith paper
-
Existence of nonrelativistic $\ell$- and multi-$\ell$-boson stars and their radial stability
Existence of an infinite class of spherically-symmetric solutions to the multi-field Schrödinger-Poisson system is established via global minimization of the energy functional on rotationally invariant H1 functions wi...
Reference graph
Works this paper leans on
-
[1]
R. Brito, V. Cardoso, C. A. R. Herdeiro, and E. Radu, Proca stars: Gravitating Bose–Einstein condensates of massive spin 1 particles, Phys. Lett. B752, 291 (2016), arXiv:1508.05395 [gr-qc]
Pith/arXiv arXiv 2016
-
[2]
I. Salazar Landea and F. Garc´ ıa, Charged Proca Stars, Phys. Rev. D94, 104006 (2016), arXiv:1608.00011 [hep- th]
Pith/arXiv arXiv 2016
-
[3]
Y. Brihaye, T. Delplace, and Y. Verbin, Proca Q Balls and their Coupling to Gravity, Phys. Rev. D96, 024057 (2017), arXiv:1704.01648 [gr-qc]
Pith/arXiv arXiv 2017
-
[4]
N. Sanchis-Gual, C. Herdeiro, E. Radu, J. C. Degollado, and J. A. Font, Numerical evolutions of spherical Proca stars, Phys. Rev. D95, 104028 (2017), arXiv:1702.04532 [gr-qc]
Pith/arXiv arXiv 2017
-
[5]
Minamitsuji, Vector boson star solutions with a quar- tic order self-interaction, Phys
M. Minamitsuji, Vector boson star solutions with a quar- tic order self-interaction, Phys. Rev. D97, 104023 (2018), arXiv:1805.09867 [gr-qc]
Pith/arXiv arXiv 2018
-
[6]
N. Sanchis-Gual, C. Herdeiro, J. A. Font, E. Radu, and F. Di Giovanni, Head-on collisions and orbital merg- ers of Proca stars, Phys. Rev. D99, 024017 (2019), arXiv:1806.07779 [gr-qc]
Pith/arXiv arXiv 2019
-
[7]
C. A. R. Herdeiro, G. Panotopoulos, and E. Radu, Tidal Love numbers of Proca stars, JCAP08, 029, arXiv:2006.11083 [gr-qc]
Pith/arXiv arXiv 2006
-
[8]
C. A. R. Herdeiro and E. Radu, Asymptotically flat, spherical, self-interacting scalar, Dirac and Proca stars, Symmetry12, 2032 (2020), arXiv:2012.03595 [gr-qc]
Pith/arXiv arXiv 2032
-
[9]
J. Calder´ on Bustillo, N. Sanchis-Gual, A. Torres-Forn´ e, J. A. Font, A. Vajpeyi, R. Smith, C. Herdeiro, E. Radu, and S. H. W. Leong, GW190521 as a Merger of Proca Stars: A Potential New Vector Boson of 8.7�10 �13 eV, Phys. Rev. Lett.126, 081101 (2021), arXiv:2009.05376 [gr-qc]
Pith/arXiv arXiv 2021
-
[10]
C. A. R. Herdeiro, A. M. Pombo, E. Radu, P. V. P. Cunha, and N. Sanchis-Gual, The imitation game: Proca stars that can mimic the Schwarzschild shadow, JCAP 04, 051, arXiv:2102.01703 [gr-qc]
-
[11]
N. Sanchis-Gual, J. Calder´ on Bustillo, C. Herdeiro, E. Radu, J. A. Font, S. H. W. Leong, and A. Torres- Forn´ e, Impact of the wavelike nature of Proca stars on their gravitational-wave emission, Phys. Rev. D106, 124011 (2022), arXiv:2208.11717 [gr-qc]
Pith/arXiv arXiv 2022
-
[12]
J. a. L. Rosa and D. Rubiera-Garcia, Shadows of boson and Proca stars with thin accretion disks, Phys. Rev. D 106, 084004 (2022), arXiv:2204.12949 [gr-qc]
Pith/arXiv arXiv 2022
-
[13]
C. A. R. Herdeiro, E. Radu, N. Sanchis-Gual, N. M. San- tos, and E. dos Santos Costa Filho, The non-spherical ground state of Proca stars, Phys. Lett. B852, 138595 (2024), arXiv:2311.14800 [gr-qc]
Pith/arXiv arXiv 2024
-
[14]
Z. Wang, T. Helfer, and M. A. Amin, General relativistic polarized Proca stars, Phys. Rev. D109, 024019 (2024), arXiv:2309.04345 [gr-qc]
Pith/arXiv arXiv 2024
-
[15]
I. Sengo, P. V. P. Cunha, C. A. R. Herdeiro, and E. Radu, The imitation game reloaded: effective shadows of dy- namically robust spinning Proca stars, JCAP05, 054, arXiv:2402.14919 [gr-qc]
-
[16]
M. Jain and M. A. Amin, Polarized solitons in higher- spin wave dark matter, Phys. Rev. D105, 056019 (2022), arXiv:2109.04892 [hep-th]
Pith/arXiv arXiv 2022
-
[17]
H.-Y. Zhang, M. Jain, and M. A. Amin, Polarized vector oscillons, Phys. Rev. D105, 096037 (2022), arXiv:2111.08700 [astro-ph.CO]
Pith/arXiv arXiv 2022
-
[18]
P. Adshead and K. D. Lozanov, Self-gravitating Vec- tor Dark Matter, Phys. Rev. D103, 103501 (2021), arXiv:2101.07265 [gr-qc]
Pith/arXiv arXiv 2021
-
[19]
M. Gorghetto, E. Hardy, J. March-Russell, N. Song, and S. M. West, Dark photon stars: formation and role as dark matter substructure, JCAP08(08), 018, arXiv:2203.10100 [hep-ph]
-
[20]
Zhang,Probing ultralight dark fields in cosmological and astrophysical systems, Ph.D
H.-Y. Zhang,Probing ultralight dark fields in cosmological and astrophysical systems, Ph.D. thesis, Rice U. (2023), arXiv:2401.00043 [hep-ph]
Pith/arXiv arXiv 2023
-
[21]
J. Chen, L. H. Nguyen, and D. J. E. Marsh, Vector Dark Matter Halo: From Polarization Dynamics to Direct De- tection, (2024), arXiv:2407.17315 [astro-ph.CO]
Pith/arXiv arXiv 2024
-
[22]
Zhang, Unified view of scalar and vector dark mat- ter solitons, (2024), arXiv:2406.05031 [hep-ph]
H.-Y. Zhang, Unified view of scalar and vector dark mat- ter solitons, (2024), arXiv:2406.05031 [hep-ph]
arXiv 2024
-
[23]
E. C. Nambo, A. Diez-Tejedor, E. Preciado-Govea, A. A. Roque, and O. Sarbach, Nonrelativistic Proca stars: Spherical stationary and multi-frequency states, (2024), arXiv:2412.06901 [gr-qc]
Pith/arXiv arXiv 2024
-
[24]
E. C. Nambo, G. Diaz-Andrade, A. Diez-Tejedor, E. Preciado-Govea, A. A. Roque, and O. Sarbach, Lin- ear stability of nonrelativistic Proca stars, (2025), arXiv:2512.04376 [gr-qc]
arXiv 2025
-
[25]
D. L´ opez Nacir and F. R. Urban, Vector Fuzzy Dark Matter, Fifth Forces, and Binary Pulsars, JCAP10, 044, arXiv:1807.10491 [astro-ph.CO]
-
[26]
M. A. Amin, M. Jain, R. Karur, and P. Mocz, Small- scale structure in vector dark matter, JCAP08(08), 014, arXiv:2203.11935 [astro-ph.CO]
-
[27]
Y.-M. Wu, Z.-C. Chen, Q.-G. Huang, X. Zhu, N. D. R. Bhat, Y. Feng, G. Hobbs, R. N. Manchester, C. J. Rus- sell, and R. M. Shannon (PPTA), Constraining ultralight vector dark matter with the Parkes Pulsar Timing Array second data release, Phys. Rev. D106, L081101 (2022), arXiv:2210.03880 [astro-ph.CO]
Pith/arXiv arXiv 2022
-
[28]
A. G. Abacet al.(KAGRA, LIGO Scientific, VIRGO), Ultralight vector dark matter search using data from the KAGRA O3GK run, Phys. Rev. D110, 042001 (2024), arXiv:2403.03004 [astro-ph.CO]
Pith/arXiv arXiv 2024
-
[29]
J. A. Dror and Q. Wei, Pulsar timing detection of ul- 12 tralight vector dark matter, Phys. Rev. D112, 075024 (2025), arXiv:2505.22719 [hep-ph]
arXiv 2025
-
[30]
T. F. Chase, D. L´ opez Nacir, and N. Yunes, Gravitational waves of quasicircular, inspiraling black hole binaries in an ultralight vector-dark-matter environment, Phys. Rev. D112, 103511 (2025), arXiv:2505.21383 [astro-ph.CO]
arXiv 2025
-
[31]
A. A. Roque, E. Ch´ avez Nambo, and O. Sarbach, Radial linear stability of nonrelativisticℓ-boson stars, Phys. Rev. D107, 084001 (2023), arXiv:2302.00717 [gr-qc]
Pith/arXiv arXiv 2023
-
[32]
Movie showing the time evolution of a 2- component multi-frequency Proca star,������ ������������������������������������������� ������������������������������������������ (2024)
2024
-
[33]
E. Ch´ avez Nambo, A. A. Roque, and O. Sarbach, Are nonrelativistic ground stateℓ-boson stars only stable forℓ=0 andℓ=1?, Phys. Rev. D108, 124065 (2023), arXiv:2310.18405 [gr-qc]
Pith/arXiv arXiv 2023
-
[34]
I. M. Moroz, R. Penrose, and P. Tod, Spherically sym- metric solutions of the Schrodinger-Newton equations, Class. Quant. Grav.15, 2733 (1998)
1998
-
[35]
J. Singh, Shooting method for solving two-point bound- ary value problems in odes numerically and applications to medical science, inArtificial Intelligence Technology in Healthcare(CRC Press, 2024) p. 121–141
2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.