{"id":"9b90e666-45ee-4383-9c79-31d7a81059e5","arxiv_id":"2411.09032","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Excited Proca stars are always unstable: under perturbations they collapse, disperse, or migrate to the ground state, with migration possible only for a narrow set of negative-binding-energy configurations.","lead":"This paper studies excited states of Proca stars, dense objects made of a massive vector field, and finds that these excited states are always unstable. In nonlinear simulations, perturbed excited Proca stars either collapse to a black hole, disperse, or migrate to the ground state, with migration limited to a narrow parameter region.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim of instability against 'even very small perturbations' is not demonstrated: the evolutions impose a fixed 5% Gaussian perturbation, and no linear stability analysis is reported.","rationale":"The reader's conditional verdict is appropriate. My stress-test identifies a somewhat different but closely related weak point: the stability claim is stronger than the evidence because only finite-amplitude 5% perturbations are simulated and no linear perturbation analysis is presented. This is a correctness risk for the central 'always unstable' statement, not merely a presentation issue. The paper's own caveat that only spherical symmetry is considered, and that ref. [25] shows the spherical ground state is unstable against non-spherical perturbations, already limits the scope; within that scope, the missing small-amplitude limit is the most load-bearing gap. I do not think the paper should be rejected on this basis: the observed collapse/migration behavior and previous bosonic-star results make the qualitative conclusion plausible, and the authors state they checked smaller perturbations. However, the headline claim, as worded, should be conditional on a proper small-amplitude verification. The migration-endpoint identification and the asserted but not shown 'dissipation' final state are additional reasons for caution, but the finite-perturbation gap is the single most load-bearing issue because it directly addresses the quantified 'against even very small perturbations' statement. Verdict stays CONDITIONAL.","tokens_in":16883,"tokens_out":6057,"duration_ms":64872,"concrete_test":"Perform an amplitude-convergence study for a representative first-excited model (e.g., φ0 = 0.006 and φ0 = 0.081): evolve the same spherical initial data with perturbation amplitudes ε = 5%, 1%, 0.1%, and 0.01% of the maximum scalar potential, using identical gauge conditions and resolutions refined consistently. Record whether the outcome (collapse, migration, or return to the excited configuration) persists and whether the instability timescale grows systematically (e.g., roughly as ε^{-2}) down to ε = 0.01%. If the 5% outcome persists and the timescale grows systematically, the 'very small perturbations' claim is supported; if the outcome changes or the star remains near the excited configuration, the claimed 'always unstable' result is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline result is that excited Proca stars are 'always unstable against even very small perturbations' (Abstract), but the reported evidence is a finite set of numerical evolutions, each started from a Gaussian perturbation of the scalar potential with amplitude fixed at 5% of the maximum of the stationary solution (Sec. IV.B, first paragraph). A 5% perturbation is finite, not 'very small'. It can trigger collapse or migration even for a configuration that is linearly stable within a finite basin, so the observations do not by themselves establish instability in the limit of infinitesimal perturbations. The statement that 'smaller perturbations give similar results' is not quantified or shown, and the observation that truncation error eventually triggers the instability does not provide a controlled amplitude dependence; at higher resolution the effective perturbation is smaller, so without an amplitude study one cannot distinguish exponential instability from a long-lived metastable process. Because the abstract quantifies over all excited configurations and all small perturbations, the missing small-amplitude limit is directly load-bearing for the central claim. A secondary weakness is that the migration endpoint is identified by matching (ωf, Mf) to the nearest sampled ground-state configuration (Table II); for several rows the mass mismatch is larger than the stated uncertainty, so the endpoint identification is not exact.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs stationary Proca star solutions in the first two excited states under spherical symmetry and studies their nonlinear stability by adding Gaussian perturbations to the scalar potential and evolving the system with the OllinSphere BSSN code. The main findings are that the first two excited families have no stable branch, and that perturbed configurations either migrate to a ground-state configuration, collapse to a black hole, or disperse. The abstract states that excited Proca stars are always unstable against even very small perturbations.","tokens_in":17102,"tokens_out":5019,"duration_ms":49197,"significance":"If the instability claim holds, the paper fills a gap in the bosonic star literature by showing that excited Proca star families lack a stable spherical branch, in contrast to the ground state. The numerical methodology is standard and carefully tested: stationary solutions are obtained by a shooting method with Hamiltonian and Gauss constraint violations around 1e-9 and fourth-order convergence, and dynamical outcomes are diagnosed by lapse collapse, apparent horizon formation, mass loss at the boundary, and FFT frequency extraction. The paper is transparent about its restriction to spherical symmetry, and it explicitly cites the result of ref. [25] that the spherical ground state is unstable against non-spherical perturbations.","major_comments":[{"comment":"All dynamical evolutions use a Gaussian perturbation with amplitude fixed at 5% of the maximum of the scalar potential, and the statement that 'smaller perturbations give similar results' is not quantified or shown. A 5% perturbation is finite, not 'very small,' so the abstract's claim that excited Proca stars are 'always unstable against even very small perturbations' is not supported by the presented evidence. The observation that truncation error eventually triggers the instability does not provide a controlled amplitude dependence: as resolution changes, the effective perturbation size changes, so the results cannot distinguish exponential instability from a long-lived metastable process kicked by a finite perturbation. I recommend adding a systematic amplitude study (e.g., amplitudes of 0.5%, 1%, 2%, and 5% for representative models) or a linear stability analysis of the radial perturbation equations; otherwise the abstract and conclusions should be reworded to state instability under finite perturbations.","section":"Sec. IV.B and Abstract"},{"comment":"The identification of the migration endpoint with a particular ground-state configuration is not quantitatively consistent for several rows. In the first excited state, the row with φ0=0.009 reports (ωf, Mf) = (0.95407±0.00024, 0.839±0.005) compared with the ground-state values (ω, M) = (0.953, 0.822); the mass differs by 0.017, more than three times the stated uncertainty, and the frequency differs by about 4.5σ. Similar discrepancies appear for the row with φ0=0.017, where Mf=0.955 versus M=0.935. The appendix text says the agreement is 'remarkable,' which overstates the actual agreement. The authors should either quantify the systematic errors in the FFT frequency and the asymptotic mass, interpolate the ground-state family to the measured (ωf, Mf), or present the endpoint as consistent within estimated systematic uncertainties if that can be justified.","section":"Table II and Appendix A"}],"minor_comments":[{"comment":"The paper notes (citing ref. [25]) that the spherical ground state is unstable against non-spherical perturbations; the conclusion's phrase 'stable branch of the ground state' should be qualified as stable within spherical symmetry to avoid misleading readers outside that restricted context.","section":"Secs. I and V"},{"comment":"There are several typos: the Figure 6 caption duplicates 'time,' the Figure 14 caption reads 'in for the models with with,' Section III uses 'anzats' instead of 'ansatz,' and page 12 has 'the the minimum.'","section":"Throughout"},{"comment":"The perturbation description says 'unit width' but does not specify whether this refers to the standard deviation or the full width at half maximum of the Gaussian; please define the width explicitly.","section":"Sec. IV.B"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of a general relativity journal and the numerical work appears technically sound. The main concern is the mismatch between the paper's strongest wording ('always unstable against even very small perturbations') and the finite-amplitude evidence shown; this should be addressed before publication. The migration-endpoint identification also needs to be presented with more careful uncertainty accounting."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a solid, within-subfield numerical relativity paper. The genuinely new piece is the construction of the first two excited Proca star families and their nonlinear evolutions. The stationary equations are derived carefully, the shooting method is standard, and the numerical checks are convincing: constraint violations around 1e-9, fourth-order convergence, and a code (OllinSphere) that has been tested in prior work. The dynamical diagnostics are also sound—lapse collapse coincides with apparent horizon formation, and the FFT-based identification of the final frequency is a nice touch. I believe the central result: in spherical symmetry, excited Proca stars are unstable and either collapse, disperse, or migrate toward the spherical ground state.\n\nThe soft spots are real but not fatal. First, the abstract claims instability against 'even very small perturbations,' but every evolution uses a 5% Gaussian perturbation. The paper says smaller perturbations give similar results and that truncation error alone triggers the instability, but there is no controlled amplitude study or linear stability analysis. That is enough to support finite-amplitude instability; it does not strictly demonstrate the infinitesimal limit. A referee should ask for an amplitude scan or a linear perturbation analysis. Second, the spherical symmetry restriction is load-bearing. The paper itself cites [25] showing that the spherical ground state is unstable against non-spherical perturbations, so the migration endpoint is only an endpoint within the spherical ansatz. In full 3D the endpoint could differ. The paper is transparent about this, but the abstract should say 'in spherical symmetry' and soften 'always.' Third, the migration endpoint matching in Table II has some slop: for several rows the final mass Mf differs from the matched ground-state mass by more than the quoted ±0.005 (e.g., 0.839 vs 0.822, 0.655 vs 0.607). The paper acknowledges the sampling limitation, so this is minor, but it means the endpoint identification is plausible rather than exact. No code or data are released, which is common but limits independent reproduction.\n\nOverall, this is a credible extension of the bosonic star program. It deserves peer review. A referee should push for an amplitude study, a more systematic parameter scan, and a more careful abstract. I would cite it as a reference for excited Proca star instability, and I would bring it to a reading group focused on exotic compact objects.","headline":"Solid spherical-symmetry result on excited Proca stars; the instability conclusion is credible, but the abstract oversells the 'even very small perturbations' and 'always' claims.","tokens_in":17612,"tokens_out":3266,"would_cite":true,"duration_ms":33696,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.20.Ex","04.25.Dm","95.30.Sf"],"model":"deepseek-v4-flash","headline":"The paper claims that the first two excited families of spherically symmetric Proca stars are always unstable under small perturbations, and that their evolutions end in black-hole collapse, dissipation, or migration to a lower-mass…","keywords":["Proca stars","excited states","bosonic stars","numerical relativity","spherical symmetry","dynamical stability","black hole formation","binding energy"],"falsifier":"Run the same perturbed first- and second-excited configurations in a full 3+1 evolution without imposing spherical symmetry: if any such configuration relaxes to a stable excited Proca star, or if a migrating star settles on a non-spherical ground state instead of the spherical one, the paper's 'always unstable, migration to the spherical ground state' claim fails.","tokens_in":16675,"feed_emoji":"🌌","tokens_out":6661,"duration_ms":59380,"temperature":0.7,"pith_summary":"This paper asks whether excited states of Proca stars—self-gravitating balls of massive complex vector field—can be stable. The answer, for the first two excited families in spherical symmetry, is no: no stable branch exists, and even tiny perturbations drive the star to collapse to a black hole, disperse, or settle onto the stable branch of the ground state. The importance of the claim is that excited configurations have been proposed as intermediate stages in the formation of exotic compact objects; if they are universally unstable, only their transient dynamics matter. The paper also identifies where migration is possible: a small corner of parameter space with negative binding energy and low mass, where the star ejects roughly a quarter to a third of its mass and lands on a known ground-state solution.","feed_headline":"Excited Proca stars are always unstable, simulations show","feed_subtitle":"First two excited families have no stable branch; perturbed stars collapse, dissipate, or settle onto the ground state.","key_machinery":"The central object is the self-gravitating Proca field in spherical symmetry, reduced by a harmonic ansatz $\\phi=\\varphi(r)e^{-i\\omega t}$, $a_r=i a(r)e^{-i\\omega t}$, $E^r=e(r)e^{-i\\omega t}$ to a first-order ODE system for the metric functions $A,\\alpha$ and the field profiles $F=\\alpha\\varphi$ and $a$, solved as an eigenvalue problem for $\\omega$ by shooting. Excited states are labelled by nodes in the vector potential $a(r)$. Stability is probed by adding a small Gaussian perturbation to the scalar potential, re-solving the Hamiltonian and Gauss constraints to get consistent initial data, and evolving with a BSSN code adapted to spherical symmetry; the end state is diagnosed through the central lapse (collapse to zero signals horizon formation), the total mass at the boundary, and a fast Fourier transform of the scalar potential at the origin whose dominant late-time frequency is matched to the stationary ground-state family.","core_discovery":"Within spherical symmetry, the paper constructs stationary families of Proca stars in the ground, first-excited, and second-excited states by solving the Einstein–Proca ODE system with a shooting method, then evolves perturbed configurations with a fully non-linear BSSN code. It finds that the first and second excited families have no stable branch analogous to the ground state: configurations with negative binding energy to the left of the maximum mass, which one might expect to be stable, are only metastable. Under a 5% Gaussian perturbation, the low-mass ones migrate to the stable branch of the ground state, identified by matching the final dominant frequency $\\omega_f$ and final total mass $M_f$ against the ground-state family, while higher-mass ones collapse to a black hole; configurations with positive binding energy dissipate. The paper therefore concludes that excited Proca stars are always unstable against small perturbations in spherical symmetry, with a three-way final fate determined by where the initial configuration sits in the family.","pith_inferences":["If this result extends beyond the first two families, excited Proca stars in general would be transient, and searches for stable bosonic dark-matter clumps should focus on ground-state configurations.","Because the paper restricts to spherical symmetry and the paper itself notes (via reference [25]) that the spherical ground state is unstable to non-spherical perturbations, the migration endpoint in a full 3D evolution could be a non-spherical ground state rather than the spherical one; the 'always unstable' verdict would then be strengthened but the migration destination would need revision.","A testable extension is to include vector self-interactions: the paper notes these can stabilize excited boson stars, but for vector fields they risk loss of hyperbolicity, so the net effect on excited Proca stars is genuinely open."],"forward_implications":["No stable branch exists for the first two excited Proca-star families; every spherically symmetric configuration tested is unstable, so these objects cannot be long-lived equilibrium states.","A perturbed excited Proca star has three possible fates: collapse to a black hole, dispersal of the field, or migration to a lower-mass configuration on the stable ground-state branch.","Migration happens only in a narrow low-mass, negative-binding-energy region; first-excited migrants lose about 25% of their mass, second-excited migrants about 35–39%, with the ejected field carrying away the excess.","The final migrated state is not arbitrary: its dominant frequency and total mass match a specific ground-state solution, so the endpoint is predictable from the family curves.","Higher excitation makes collapse more likely; in the second excited state only very low-mass configurations show migration, and the metastable window shrinks."],"supporting_citations":[{"why":"Introduces Proca stars as self-gravitating solutions of a massive complex vector field, the objects under study.","marker":"[3]"},{"why":"Performs numerical evolutions of ground-state Proca stars and studies their stability, providing the evolution methodology.","marker":"[17]"},{"why":"Establishes ground-state stability results and collapse/migration behavior used as the baseline comparison for excited states.","marker":"[18]"},{"why":"Studies excited boson stars with shell structure and their dynamics, the scalar analogue motivating this work.","marker":"[24]"},{"why":"Shows the spherically symmetric ground state is unstable against non-spherical perturbations, the key caveat for the spherical analysis.","marker":"[25]"},{"why":"Supplies the BSSN spherical-symmetry evolution code used for the dynamical simulations.","marker":"[29]"},{"why":"Shows excited boson stars decay to the fundamental state or collapse to a black hole, the scalar counterpart of the three-way fate.","marker":"[30]"},{"why":"Provides ground-state Proca star data used to validate the stationary solutions and to identify migration endpoints.","marker":"[43]"}],"fun_headline_variants":["Excited Proca stars always unstable, even to tiny nudges","Excited Proca stars collapse, dissipate, or migrate - never stable","Excited Proca stars: always unstable, three possible fates","Even tiny perturbations doom excited Proca stars","No stable excited Proca stars: they collapse, dissipate, or migrate"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that spherical symmetry is preserved throughout the evolution, so the three-way classification (black-hole collapse, dissipation, migration to the spherical ground state) captures the true dynamical attractors; the paper itself notes, citing reference [25], that the spherical ground state is unstable to non-spherical perturbations, so the endpoint in the full theory could differ.","fun_headline_variants_meta":{"raw":{"variants":["Excited Proca stars always unstable, even to tiny nudges","Excited Proca stars collapse, dissipate, or migrate - never stable","Excited Proca stars: always unstable, three possible fates","Even tiny perturbations doom excited Proca stars","No stable excited Proca stars: they collapse, dissipate, or migrate"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000719,"raw_usage":{"total_tokens":3171,"prompt_tokens":827,"completion_tokens":2344,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":443,"completion_tokens_details":{"reasoning_tokens":2254}},"tokens_in":443,"tokens_out":2344,"duration_ms":15174,"temperature":1.0,"reasoning_tokens":2254,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:08:22.677677+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same perturbed first- and second-excited configurations in a full 3+1 evolution without imposing spherical symmetry: if any such configuration relaxes to a stable excited Proca star, or if a migrating star settles on a non-spherical ground state instead of the spherical one, the paper's 'always unstable, migration to the spherical ground state' claim fails.","supporting_citations":[{"cited_title":"In the figure we also indicate 7 different models that we considered for our dynamical evolutions of perturbed initial data","cited_arxiv_id":null,"evidence_quote":"Introduces Proca stars as self-gravitating solutions of a massive complex vector field, the objects under study."},{"cited_title":"Sanchis-Gual, C","cited_arxiv_id":null,"evidence_quote":"Performs numerical evolutions of ground-state Proca stars and studies their stability, providing the evolution methodology."},{"cited_title":"Lazarte and M","cited_arxiv_id":null,"evidence_quote":"Shows the spherically symmetric ground state is unstable against non-spherical perturbations, the key caveat for the spherical analysis."},{"cited_title":"Herdeiro, E","cited_arxiv_id":null,"evidence_quote":"Supplies the BSSN spherical-symmetry evolution code used for the dynamical simulations."},{"cited_title":"Alcubierre, Introduction to 3 + 1 Numerical Relativity (Oxford Univ","cited_arxiv_id":null,"evidence_quote":"Shows excited boson stars decay to the fundamental state or collapse to a black hole, the scalar counterpart of the three-way fate."}],"review_version":1}