{"id":"906fb875-d6a0-46b0-9993-4cce3e2f8b1b","arxiv_id":"2412.06901","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Nonrelativistic Proca stars have a ground state with constant polarization (linear or circular depending on the sign of the spin-spin coupling), and a symmetry-enhanced sector with λs=0 contains a continuum of multi-frequency solutions connecting stationary states.","lead":"This paper derives the nonrelativistic limit of a self-gravitating, self-interacting massive vector field and classifies its equilibrium configurations, which are called Proca stars. It proves that the lowest-energy state for fixed particle number is a spherically symmetric star with constant polarization, and finds a new family of multi-frequency stars when the spin-spin interaction is tuned to zero.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ground-state theorem is sound, but the headline multi-frequency continuum is confined to the exact λs=0 surface, where the accidental U(3) symmetry holds; any nonzero spin-spin coupling removes it exactly, and the numerical shooting evidence is not independently verified.","rationale":"The strongest analytical claim, that for λ0≥0 a ground state exists and is a stationary, spherically symmetric constant-polarization Proca star, is well argued and does not hinge on the fine-tuning of λs. The energy inequality, the polarization cases in Appendix C, and the appeal to symmetric-decreasing rearrangement and Choquard-type minimizers are credible and internally consistent. The soft spot is exactly the multi-frequency continuum, as the reader's weakest assumption states. The continuum requires the accidental U(3) symmetry, and Eq. (12) shows Q is conserved only at λs=0. For any nonzero λs, even infinitesimal, the generic-sector analysis in Section III B forces every equilibrium at fixed N to be a single-frequency stationary state; the multi-frequency ansatz cannot satisfy Eq. (2a) because the spin density becomes time-dependent when the component frequencies differ. Hence the headline discovery is tied to a codimension-one, fine-tuned surface in parameter space. In addition, the only evidence for the continuum is the numerical shooting in Section V, with no code or convergence analysis, so the one-parameter family in Fig. 10 is not independently verified. These considerations do not overturn the paper: they qualify it. The reader's CONDITIONAL verdict is appropriate, and I would not change it. The paper should clearly present the multi-frequency continuum as a symmetry-enhanced limiting case and provide reproducible numerics or a convergence study before the continuum claim is treated as established.","tokens_in":37391,"tokens_out":11643,"duration_ms":146770,"concrete_test":"Numerical robustness test: initialize the (nx,ny,nz)=(0,1,0) multi-frequency state from Fig. 10 in the full s=1 Gross-Pitaevskii-Poisson system with λs=10^-3 λn and fixed λn, using the same initial σi(r). Evolve for several periods 2π/(Ex-Ey). If no nearby equilibrium exists, the state will precess, shed energy, or fail to maintain the time-periodic form e^{-iE_i t}σi(r), whereas the same initial data at λs=0 is exactly stationary. Comparing these evolutions directly tests whether the multi-frequency continuum is a robust physical family or only a measure-zero artifact of the accidental U(3) symmetry.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III B's global-minimum theorem is well supported: the reduction to the scalar functional Escalar via λ0 in Eq. (32), the constant-polarization equality conditions, and the symmetric-decreasing rearrangement argument are standard and internally consistent. The load-bearing weak point is the continuum of multi-frequency states claimed in Sections III C, IV B, and V (Figs. 6, 10, 15). These states are equilibria only as critical points of E_Q with a matrix Lagrange multiplier Ê, which is justified only when λs=0 and the tensor Q is conserved (Eq. 12). For any λs≠0, the U(3) symmetry is broken, Q reduces to N and S, and the matrix-multiplier construction is no longer available. The paper's own generic-sector argument (Eqs. 21-22) shows every critical point of E at fixed N obeys Hψ=Eψ and hence evolves as a single-frequency stationary state. Thus an arbitrarily small λs destroys the continuum exactly, not merely perturbatively: a two-component multi-frequency state with real σx,σy and Ex≠Ey has vanishing spin at t=0, but its spin density becomes s_z(t)=2σxσy sin((Ex-Ey)t) under time evolution, so the spin-spin term cannot remain time-independent. The multi-frequency discovery is therefore tied to a measure-zero parameter surface. Separately, the one-parameter family in Fig. 10 is established only by a shooting method; no solver, residuals, error analysis, or convergence study is provided, so the existence of the continuum is not independently confirmed. The central ground-state theorem does not depend on this fine-tuning and deserves credit.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a nonrelativistic effective field theory for a selfgravitating, selfinteracting massive vector field, the s=1 Gross-Pitaevskii-Poisson system, and studies its spherically symmetric equilibrium configurations. In the generic sector (λs≠0) the admissible equilibria are stationary, single-frequency states; in the symmetry-enhanced sector (λs=0) the theory acquires an accidental U(3) symmetry and admits, in addition, multi-frequency states whose components oscillate with two or three distinct frequencies. The authors prove that for λ0≥0 (λ0=λn for λs≥0 and λn−|λs| for λs<0) a global minimum of the energy at fixed particle number exists and is a stationary, spherically symmetric, constant-polarization state, linearly polarized for λs>0 and circularly polarized for λs<0. They then numerically construct stationary and multi-frequency solutions, reporting a one-parameter family of multi-frequency states with fixed particle number that fills a region in the mass-radius and energy diagrams, connecting the nodeless and first-excited linearly polarized stationary states.","tokens_in":37694,"tokens_out":6383,"duration_ms":65451,"significance":"The analytic ground-state theorem in Sec. III B is a rigorous, well-grounded result: it reduces the vector problem to the scalar Choquard functional, uses symmetric decreasing rearrangements to enforce spherical symmetry, and identifies the polarization through the equality conditions of the spin inequality. This is a genuine variational characterization, not a fit, and it yields a falsifiable prediction (linear versus circular polarization depending on the sign of λs). The multi-frequency states are novel and connect the Proca-star literature to multi-state boson stars; the numerical results suggest a rich solution space in the symmetry-enhanced sector. However, the existence of the multi-frequency continuum is established only numerically, and the continuum itself is confined to the exactly fine-tuned surface λs=0, which substantially limits its physical robustness. If the numerical evidence is validated and the fine-tuning caveat is made prominent, the paper would be a solid contribution to the nonrelativistic vector dark matter and soliton literature.","major_comments":[{"comment":"The existence of the one-parameter continuum of multi-frequency states is a central claim of the paper, but it is supported only by a shooting method. No residuals, convergence tests, mesh-dependence checks, or independent verification of the eigenvalues are reported, and the code is not made available. Since a shooting method can in principle produce spurious families near turning points or lose accuracy for large radii, please provide a quantitative accuracy assessment (for example, ODE residuals after the numerical solution, comparison of the energy eigenvalues computed from Eq. (E2a) with those obtained from the asymptotic fit, or a convergence study under increasing integration interval and decreasing step size). Without such validation, the existence of a genuine continuum rather than a discrete set of numerically connected points remains unsubstantiated.","section":"Sec. V, Figs. 10 and 15"},{"comment":"The multi-frequency equilibria and the continuum connecting stationary states require λs=0 exactly, because the tensor Q is conserved only in that case (Eq. (12)). The paper correctly restricts the symmetry-enhanced sector to λs=0, but the abstract and conclusions present the continuum as a headline result without emphasizing that any nonzero λs, no matter how small, destroys these states exactly, not perturbatively: the generic-sector argument in Sec. III implies that all fixed-N equilibria are single-frequency stationary states for λs≠0. Please add an explicit discussion of this fine-tuning issue in the abstract and conclusions, and either temper the claim or justify why the λs=0 surface is physically relevant (for example, by a symmetry-protection argument or by showing that the continuum approximately survives for small λs in a perturbative sense).","section":"Sec. II B (Eq. (12)) and Sec. III C"}],"minor_comments":[{"comment":"The text contains two consecutive 'open question' passages at the end of the section that are redundant and could be merged into a single statement.","section":"Sec. III C"},{"comment":"The time labels t1, t2, t3 are defined with E for the stationary rows but with Ex for the multi-frequency row; please use a consistent notation and clarify that E is the single frequency for stationary states.","section":"Fig. 2 caption"},{"comment":"The symbol 'i' is used both for the imaginary unit and as an index (e.g., ψi), which is occasionally confusing; consider using a roman upright i for the imaginary unit or a different index letter.","section":"Sec. I, Eq. (1)"},{"comment":"The abbreviation 'c.f.' should be 'cf.' in several places (for example, 'c.f. Eqs. (32) in Ref. [45]').","section":"Sec. IV A, after Eq. (45)"},{"comment":"The '±' sign is introduced as referring to repulsive and attractive selfinteractions, but the relation to the earlier parameter λ0 in Eq. (32) and to λ* is not spelled out; a brief note connecting both would improve readability.","section":"Sec. V, Eq. (52a)"}],"recommendation":"major_revision","confidential_remarks":"The analytic part is the main strength and is, in my reading, correct and carefully argued. The numerical part is the weakest link: the existence of the continuum is a headline claim that rests entirely on a shooting code with no accuracy assessment, and the fine-tuning of λs=0 is underplayed in the presentation. If the authors can provide a credible numerical validation and rebalance the claims, the paper could be acceptable for publication. I would not reject on the fine-tuning issue alone, since the symmetry-enhanced sector is a clearly defined limit and the paper is explicit about its assumptions, but the current presentation oversells the multi-frequency continuum as a robust discovery rather than a property of a measure-zero parameter surface."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Good paper for the subfield, and for once the central theorem is worth the attention. The ground-state result is real: when λ0≥0, the minimizer of the energy for fixed particle number is a spherical, stationary, constant-polarization Proca star — linear for λs>0, circular for λs<0. The argument is standard but careful, using the Choquard energy, concentration-compactness, and the polarization inequality, and it settles a question the literature had left open. The authors also deserve credit for being explicit about what is not new: the λn=0 multi-frequency states are just a relabeling of known multi-state boson stars, and they map it out in Table II. The self-interacting radial-polarized stationary stars and the λn≠0 multi-frequency family are genuinely new.\n\nThe soft spot is the paper's most advertised claim. The multi-frequency continuum exists only on the λs=0 surface, where the accidental U(3) symmetry makes the full tensor Q conserved. Any nonzero λs, however small, removes that symmetry, kills the matrix Lagrange multiplier construction, and in the simple two-component example the spin density s_z(t) = 2σxσy sin((Ex−Ey)t) becomes time-dependent, so the configuration is not an equilibrium at all. The paper states the λs=0 condition but never discusses the fate of these states under small explicit breaking. That makes the continuum a genuine but measure-zero feature of the effective theory. The numerical evidence is also thinner than the rest of the paper: a shooting method with asymptotic matching, no residuals, no convergence study. Fine for illustrative plots, not quite enough to call the continuum verified. A referee should ask for a short section on the λs→0 limit and either the code or a convergence analysis.\n\nThe ground-state theorem does not depend on that fine-tuning, and the stationary-state numerics line up with it. I also checked the polarization appendix; the equality conditions are correct. The paper is honest about its open questions, including ground-state uniqueness and linear stability.\n\nRecommendation: send it to peer review. It deserves referee time. The theorem and the effective-theory mapping are solid, and the multi-frequency claim, once the measure-zero caveat is made explicit, is a legitimate mathematical observation. I'd cite the ground-state theorem; the continuum I'd cite only with the caveat attached.","headline":"Solid variational ground-state theorem, honest mapping to known boson stars, but the multi-frequency continuum is confined to the λs=0 fine-tuned surface and not yet numerically verified.","tokens_in":38269,"tokens_out":3727,"would_cite":true,"duration_ms":38368,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Proca star ground states are spherical and polarized.","keywords":["Proca stars","nonrelativistic boson stars","Gross-Pitaevskii-Poisson system","vector dark matter","multi-frequency states","ground state","polarization","U(3) symmetry"],"falsifier":"Adapt the shooting code to a small nonzero $\\lambda_s$ (for instance $\\lambda_s=10^{-6}\\lambda_n$) and search for a branch of multi-frequency solutions of Eq. (48) with distinct frequencies $E_x\\neq E_y$; a persistent branch would falsify the claim that multi-frequency states require $\\lambda_s=0$. Alternatively, find a global minimizer for $\\lambda_0\\ge 0$ that is not a stationary, spherically symmetric constant-polarization state, for example a multi-frequency or radially polarized configuration with energy below the constantly polarized nodeless state at the same $N$.","tokens_in":37201,"feed_emoji":"🌟","tokens_out":11720,"duration_ms":104981,"temperature":0.7,"pith_summary":"This paper studies the nonrelativistic limit of a self-gravitating, self-interacting spin-1 (vector) field, the setting in which ultralight vector dark matter would form galactic halos. Its central result is an existence and characterization theorem: when the effective coupling $\\lambda_0 = \\lambda_n$ for $\\lambda_s \\ge 0$ and $\\lambda_0 = \\lambda_n - |\\lambda_s|$ for $\\lambda_s < 0$ is nonnegative, the energy at fixed particle number is bounded below and has a global minimum, and every such minimum is a stationary, spherically symmetric Proca star of constant polarization (linear if $\\lambda_s>0$, circular if $\\lambda_s<0$). The paper also identifies a symmetry-enhanced sector $\\lambda_s=0$ in which the theory has an accidental U(3) symmetry, allowing equilibrium configurations whose wave function oscillates at two or three distinct frequencies. These multi-frequency states form a continuum in solution space that connects stationary states of constant polarization, whereas stationary states alone form a discrete set for fixed particle number. A sympathetic reader should care because this fixes the expected ground-state shape of vector dark matter halos and predicts a new family of equilibria if spin-spin couplings are absent.","feed_headline":"Proca star ground states are spherical and polarized","feed_subtitle":"New proof fixes the lowest-energy shape of vector dark matter halos and uncovers a continuum of multi-frequency states.","key_machinery":"The load-bearing object is the $s=1$ Gross-Pitaevskii-Poisson system, a nonlinear Schrodinger-type equation for a three-component complex vector field $\\vec{\\psi}(t,\\vec{x})$ coupled to a Newtonian potential $U$ satisfying Poisson's equation, with two self-interaction parameters $\\lambda_n$ (density-density) and $\\lambda_s$ (spin-spin). The argument runs through the energy functional $E=T+\\lambda_n F_n+\\lambda_s F_s-D$, its behaviour under the scaling $\\vec{\\psi}\\mapsto \\nu^{3/2}\\vec{\\psi}(\\nu\\vec{x})$, and the effective coupling $\\lambda_0$ defined in Eq. (32), which decides whether $E$ is bounded below at fixed $N$. For the minimizer claim, the paper uses a polarization decomposition with the inequality $\\lambda_n F_n+\\lambda_s F_s\\ge \\lambda_0 F_n$ and the symmetric decreasing rearrangement to reduce any minimizer to a spherical, constantly polarized profile. For the multi-frequency sector, the machinery is the accidental U(3) symmetry present when $\\lambda_s=0$, with the conserved tensor $\\hat{Q}=\\int \\vec{\\psi}^*\\otimes\\vec{\\psi}\\,dV$, whose diagonalization turns the equilibrium equations into a nonlinear multi-eigenvalue problem with frequencies $E_\\lambda$; the Sturm oscillation theorem then orders the radial components by node number and shows that components with equal node numbers are proportional, hence stationary.","core_discovery":"On the paper's own terms: for the $s=1$ Gross-Pitaevskii-Poisson system describing a nonrelativistic self-gravitating vector field, all equilibrium configurations in the generic sector ($\\lambda_s\\neq 0$) are stationary states $\\vec{\\psi}(t,\\vec{x})=e^{-iEt}\\sigma(r)\\hat{\\epsilon}$, while in the symmetry-enhanced sector ($\\lambda_s=0$) there are also multi-frequency states $\\vec{\\psi}(t,\\vec{x})=\\sum_{\\lambda=1}^3 e^{-iE_\\lambda t}\\sigma_\\lambda(r)\\hat{e}_\\lambda$. The proof shows that when $\\lambda_0\\ge 0$ the energy functional at fixed $N$ is bounded below and attains its minimum at a stationary, spherically symmetric state of constant polarization with monotonically decreasing positive radial profile; for $\\lambda_s>0$ the polarization is linear, for $\\lambda_s<0$ it is circular, and in the free theory ($\\lambda_n=\\lambda_s=0$) this ground state is unique up to translations and rigid unitary transformations. In the symmetry-enhanced sector, multi-frequency solutions are shown to fill regions of the mass-radius and energy-particle-number diagrams, bounded by the $n=0$ and $n=1$ constant-polarization stationary states; the paper constructs examples numerically, computes their charges and energies, and verifies that the constantly polarized nodeless state has the lowest energy.","pith_inferences":["Beyond the paper: if ultralight vector dark matter has any nonzero spin-spin coupling $\\lambda_s$, multi-frequency halos are forbidden, so detecting a stably oscillating multi-frequency core would place an upper bound on $|\\lambda_s|$ and effectively confirm an accidental U(3) symmetry in the low-energy theory.","Beyond the paper: the continuous family of multi-frequency states that links ground and excited stationary states suggests that slow dynamical processes could drive a Proca star from one stationary branch to another; the announced linear-stability analysis should reveal whether these states are orbitally stable or decay toward the ground state.","Beyond the paper: because constantly polarized Proca stars coincide with $\\ell=0$ boson stars (with coupling $\\lambda_0$) and radial Proca stars coincide with $\\ell=1$ boson stars, existing results on scalar boson star stability, collisions, and gravitational-wave signals can be imported to the vector case at least in the free and $\\lambda_s=0$ sectors.","Beyond the paper: the proof leaves open whether the ground state is unique for $\\lambda_0>0$; a systematic numerical search for minimizers with different radial profiles but the same particle number and energy would settle the uniqueness question left open in the paper."],"forward_implications":["For $\\lambda_0\\ge 0$, a Proca star with fixed particle number always possesses a ground state, and that ground state is spherical, stationary, and constantly polarized with negative energy, providing a candidate stable endpoint for vector dark matter condensation.","If spin-spin self-interactions are absent ($\\lambda_s=0$), equilibrium solutions are not isolated: multi-frequency states form continuous families that interpolate between the $n=0$ and $n=1$ constant-polarization stationary states in the mass-radius and energy diagrams.","Radially polarized Proca stars are excited states relative to constantly polarized ones when $\\lambda_0\\ge 0$, so the nonrelativistic limit of the original spherically symmetric relativistic Proca stars is an excited configuration rather than a ground state.","In the generic sector $\\lambda_s\\neq 0$, no multi-frequency states exist; only stationary linear, circular, or radial polarizations are possible, and the degeneracy among constant-polarization states is broken.","Circularly polarized Proca stars carry macroscopic spin angular momentum $\\vec{S}=\\alpha N\\hat{e}_z$ despite a spherically symmetric density, implying the corresponding relativistic rotating stars should be non-spherical."],"supporting_citations":[{"why":"Introduces Proca stars as gravitating Bose-Einstein condensates of massive spin-1 particles and supplies the relativistic spherical solutions whose nonrelativistic limit corresponds to radial-polarization states.","marker":"[10]"},{"why":"Establishes polarized solitons in higher-spin wave dark matter, the earlier nonrelativistic framework whose stationary polarized states this paper extends.","marker":"[31]"},{"why":"Derives the nonrelativistic limit of a massive vector field and constructs polarized vector oscillons, providing the baseline parameterization with $\\lambda_n$ and $\\lambda_s$.","marker":"[32]"},{"why":"Supplies the multi-scalar $\\ell$-boson star framework and the shooting method used to construct spherical stationary solutions and label them by node number.","marker":"[42]"},{"why":"Provides the multi-state boson star equations that, in the $\\lambda_s=0$ sector, map exactly onto the multi-frequency Proca star equations.","marker":"[44]"},{"why":"Gives the self-interacting $s=0$ Gross-Pitaevskii-Poisson ground-state analysis and numerical curves that linearly and circularly polarized Proca stars reproduce.","marker":"[45]"},{"why":"Proves existence and uniqueness of the minimizer for Choquard's nonlinear equation, the scalar minimization problem that underwrites the Proca ground-state theorem.","marker":"[46]"},{"why":"Provides the concentration-compactness principle used to establish existence of minimizers of the scalar energy functional.","marker":"[48]"},{"why":"Complements Ref. [48] in the existence argument for the same minimization problem.","marker":"[49]"},{"why":"Supplies the Sturm oscillation and comparison theorems used to order radial eigenfunctions by node number and to show equal-node components must be proportional.","marker":"[50]"}],"fun_headline_variants":["Proca star ground state proven spherical and polarized","Nonrelativistic Proca stars: multi-frequency states and proof","Vector field halos: spherical ground states and continuum of states","Proca stars: stationary and multi-frequency equilibrium states"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The continuum of multi-frequency Proca stars rests on the spin-spin coupling $\\lambda_s$ being exactly zero; any nonzero $\\lambda_s$, no matter how small, breaks the U(3) symmetry and eliminates those states, although the ground-state theorem itself does not depend on this fine-tuning.","fun_headline_variants_meta":{"raw":{"variants":["Proca star ground state proven spherical and polarized","Nonrelativistic Proca stars: multi-frequency states and proof","Vector field halos: spherical ground states and continuum of states","Proca stars: stationary and multi-frequency equilibrium states"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000848,"raw_usage":{"total_tokens":3781,"prompt_tokens":1126,"completion_tokens":2655,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":742,"completion_tokens_details":{"reasoning_tokens":2588}},"tokens_in":742,"tokens_out":2655,"duration_ms":20709,"temperature":1.0,"reasoning_tokens":2588,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T19:20:14.840095+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Adapt the shooting code to a small nonzero $\\lambda_s$ (for instance $\\lambda_s=10^{-6}\\lambda_n$) and search for a branch of multi-frequency solutions of Eq. (48) with distinct frequencies $E_x\\neq E_y$; a persistent branch would falsify the claim that multi-frequency states require $\\lambda_s=0$. Alternatively, find a global minimizer for $\\lambda_0\\ge 0$ that is not a stationary, spherically symmetric constant-polarization state, for example a multi-frequency or radially polarized configuration with energy below the constantly polarized nodeless state at the same $N$.","supporting_citations":[{"cited_title":"Probing ultralight dark fields in cosmological and astrophysical systems","cited_arxiv_id":"2401.00043","evidence_quote":"Provides the multi-state boson star equations that, in the $\\lambda_s=0$ sector, map exactly onto the multi-frequency Proca star equations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Proves existence and uniqueness of the minimizer for Choquard's nonlinear equation, the scalar minimization problem that underwrites the Proca ground-state theorem."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the concentration-compactness principle used to establish existence of minimizers of the scalar energy functional."},{"cited_title":"Matos and L","cited_arxiv_id":null,"evidence_quote":"Complements Ref. [48] in the existence argument for the same minimization problem."}],"review_version":1}