{"id":"7198e710-cacf-46de-acee-cbfda7301314","arxiv_id":"2608.06067","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"For excited boson stars, the first zero of the fundamental radial mode matches the first critical point of mass, charge, and binding energy for all tested node numbers and self-interactions.","lead":"This paper computes the lowest radial vibration mode of excited boson stars, using a new technique that stays regular even when the star's field profile has zeros. The main finding is that the mode's stability change coincides with an equilibrium turning point, offering a fast spectral way to diagnose when these compact objects become unstable.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central claim's breadth is unsupported: zero-crossing data for n=4..7 are not shown, and 'within numerical resolution' is never quantified.","rationale":"The reader's weakest assumption concerns the completeness and faithfulness of the reformulated perturbation system, which is indeed fundamental: a flawed eigenvalue problem would invalidate every numerical result. However, the paper's nodeless benchmark (n=0, critical frequency omega=0.853) matches known results, giving some confidence in the framework. The most concretely load-bearing concern for the central claim is the breadth-versus-evidence gap: the headline coincidence is asserted for n up to 7, but only n=0..3 zero-crossing data are shown, and the numerical resolution is never defined. This is a missing-support passage that the manuscript itself implicitly acknowledges. The review rules require flagging such passages explicitly, and they are directly tied to the central claim. A conditional accept remains appropriate because the concern is about verification and documentation rather than a demonstrated internal inconsistency, but the higher-node claim and the 'within numerical resolution' wording should be substantiated before the abstract can be taken at face value.","tokens_in":18271,"tokens_out":27482,"duration_ms":261134,"concrete_test":"Recompute, with the same shooting method, the zero crossing of chi^2_0,n and the first critical point for n=4,5,6,7 at eta=0 and eta=58, and report for each case: omega_zero, omega_crit from M/Q/E_b, the difference Delta_omega, and the numerical uncertainty estimated from varying r_out over the plateau and from grid resolution. Also report epsilon(0; r_out) at the claimed zero and at a nearby non-zero configuration. If any Delta_omega exceeds the estimated uncertainty, or if epsilon is not selectively small only at the claimed zero, the central claim must be restricted to n<=3 or qualified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract and conclusions claim that for all branches examined the first zero of the constrained fundamental radial eigenvalue coincides with the first critical point of M, Q, and E_b. The only zero-crossing data actually displayed are in Tables II–V for n=0..3. For n=4..7, Sec. IV.A states only that the calculation was extended and 'these results are therefore not displayed separately'—an explicit missing-support assertion. Table I contains positive eigenvalues at nonlinear dynamical thresholds, not the zero-crossing versus critical-frequency comparison needed for the central claim. In addition, the paper never quantifies 'within numerical resolution': no error bars, grid-convergence data, r_out/plateau convergence studies, or thresholds for the singular-value ratio epsilon(0; r_out) used in the rank-condition diagnostic. If the true zero-critical frequency difference for any n>=4 were larger than the displayed four-significant-digit agreement, the abstract's 'all branches' statement would be false; as written, that possibility cannot be checked from the manuscript.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a regular radial perturbation framework for spherically symmetric equilibrium boson stars with radial nodes. Introducing additive variables f and g instead of the conventional relative variables, the authors write a closed first-order system for the metric and scalar perturbations that is regular at the zeros of the background scalar field, integrate it through the nodes of excited configurations, and compute the lowest radial eigenvalue chi^2_{0,n} along fixed-node branches. For n=0,...,3 and several values of the quartic self-interaction coupling, they report a zero crossing of chi^2_{0,n} and claim that this zero coincides, within numerical resolution, with the first critical point of the ADM mass, Noether charge, and binding energy. They further evaluate chi^2_{0,n} at the nonlinear stability thresholds of Brito et al. and fit an empirical scaling relation in the node number and coupling, testing it on n=5,...,7.","tokens_in":18480,"tokens_out":6893,"duration_ms":73948,"significance":"If the central claim is correct, the paper provides a useful technical advance: the regular variable formulation permits a direct integration of the radial perturbation equations through the nodes of excited backgrounds, and it ties the constrained radial zero mode to equilibrium turning points for a class of configurations. The perturbation system is written out in full, and the tables for n=0,...,3 give a substantial amount of reproducible numerical data. The empirical correlation with nonlinear threshold models is interesting and, unusually, is tested on configurations not used in the fit. The result is a plausible and potentially useful diagnostic, but the breadth of the central claim currently exceeds the evidence displayed in the manuscript.","major_comments":[{"comment":"The abstract and conclusions state that 'for all branches examined' the first zero of the fundamental radial eigenvalue coincides with the first critical point, but the only displayed zero-crossing data are for n=0,...,3 in Tables II-V and Fig. 3. The sentence 'We have extended the calculation to higher excited states up to n=7 and found the same behavior in all cases examined; these results are therefore not displayed separately' is not sufficient support for the central claim as worded. At minimum, the paper should provide a table or figure listing the zero-crossing frequencies and the corresponding critical frequencies for n=4,...,7, or the abstract and conclusions should be restricted to the branches for which data are shown.","section":"Sec. IV.A, paragraph after Fig. 3; Tables II-V"},{"comment":"The phrase 'within numerical resolution' is used in the abstract and conclusions but is never quantified. No grid-convergence study, no statement of the shooting accuracy, no values of r_out/R99 used for the plateau, and no thresholds for the singular-value ratio epsilon(0;r_out) in Eq. (59) are given. The tabulated zero frequencies and critical frequencies agree to four or five significant digits, but without an error estimate or a bracket size for the zero crossing, the coincidence claim cannot be assessed quantitatively. Please provide a quantitative measure of the numerical uncertainty in the zero-crossing location for representative cases.","section":"Sec. III.D and Sec. IV.A"},{"comment":"The identification of chi^2_{0,n} as the fundamental (lowest) radial eigenvalue along the entire branch is asserted rather than demonstrated. The text states that the 'lowest real root' is located on a reference background and followed continuously, with selected anchor points rescanned to verify that no lower root was missed, but no details are given on the scanning density, the mode spacing, or how mode crossings or avoided crossings were excluded. Since the central claim concerns the 'first zero of the constrained fundamental radial eigenvalue,' the completeness of this search is load-bearing. Please provide explicit evidence, such as the scanning interval and the anchor-point rescan results, that no lower real root exists in the omega range traversed.","section":"Sec. III.C, paragraph on eigenvalue search"}],"minor_comments":[{"comment":"Because Eq. (28) follows from the first law, the extrema of M, Q, and E_b coincide by construction for any regular parameter along the branch. The wording 'independently identified from the ADM mass, the Noether charge, and the binding energy' therefore overstates the independence of these three diagnostics; consider rephrasing to make clear that the non-trivial comparison is between the chi^2 zero and this common equilibrium critical point.","section":"Eq. (28) and Sec. IV.A"},{"comment":"The notation '751.4556e-4' in Table I is ambiguous: it should be clear that the first column entry is eta=75 and the second is chi^2=1.4556e-4. Please separate the columns or use explicit notation.","section":"Table I"},{"comment":"The empirical fits are reported with R^2 values but without uncertainties on the fitted prefactors and exponents; given that only four data points are used for a two-parameter fit, these uncertainties should be quoted if the scaling is presented as predictive.","section":"Eqs. (60)-(61)"},{"comment":"The caption states that dashed segments are guides to the eye, but the zero crossing is read from interpolation between solid data points. Please state explicitly how the zero-crossing frequency was extracted from the plotted data.","section":"Fig. 3 caption"},{"comment":"The DOI '10.1103/2dhs-phl4' appears malformed; please verify the bibliographic data.","section":"Reference [42]"}],"recommendation":"major_revision","confidential_remarks":"The main reservation is that the central 'all branches examined' claim is broader than the displayed evidence, and the 'within numerical resolution' qualifier is not quantified. Both are fixable within the manuscript's scope by adding the n=4,...,7 data and a convergence/error analysis. The regular-variable formulation itself appears sound and is a useful contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the new thing: the additive variables f and g that stay regular at nodes solve a real pain point for excited boson star perturbation theory. Gleiser\\u2013Kain relative variables blow up at zeros of phi0; this reformulation lets you shoot straight through the nodes, and the equations are written out in full. That is a concrete advance, and the nodeless benchmark works. The systematic fundamental-mode curves for n=1..3 across eta values are also new, and the tables give enough entries to be checked.\\n\\nThe central claim is credible for the data actually shown. The zero crossing of chi^2_{0,n} is computed independently of the equilibrium critical point (which comes from M, Q, E_b), so there is no circularity. For n=0..3 the tables match to four significant digits. Good evidence that the turning-point/zero-mode correspondence extends to nodeful branches, at least for these cases.\\n\\nNow the soft spots. First, \"within numerical resolution\" is never quantified. No error bars, no grid or r_out convergence tables, no threshold for the singular-value ratio epsilon. The r_out/R99 plateau is described but not demonstrated. That is a fixable but real gap. Second, the n=4..7 extension is asserted but not shown. The sentence \"these results are therefore not displayed separately\" is doing load-bearing work for the abstract's \"all branches examined.\" Table I gives eigenvalues at nonlinear thresholds, but those are not zero-crossing comparisons against critical frequencies. If the n=4..7 behavior differs from n=0..3, the abstract would overclaim; as written, we cannot verify. Third, the empirical scaling in Eqs. (60)-(61) is honestly labeled and out-of-sample tested, which is good, but it remains an empirical fit, not a derivation. They do not oversell it.\\n\\nMinor: no code release, and the rank-condition diagnostic is reasonable but the scale sensitivity of epsilon is not discussed. Those are minor compared with the convergence quantification.\\n\\nOverall, this is a solid paper worth serious refereeing. The method is reusable, the data for n=0..3 is genuine, and the central correspondence seems plausible where actually displayed. The package is not ready as is\\u2014the convergence documentation and the n=4..7 support need addressing\\u2014but the core is sound. Who is it for? Anyone working on boson star stability or exotic compact objects; a GR group reading group could use it as an example of a regular perturbation formulation. Send it to peer review, conditional on the authors providing quantified convergence and either showing or explicitly demoting the n=4..7 claim.","headline":"Regular variable method is a genuine advance and the n=0..3 zero-crossing data look credible, but the \"all branches\" claim outruns the displayed data and \"within numerical resolution\" needs quantification.","tokens_in":749,"tokens_out":954,"would_cite":true,"duration_ms":32735,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.40.-b"],"model":"deepseek-v4-flash","headline":"The first zero of the fundamental radial mode coincides with the first critical point of the equilibrium branch for every boson star configuration examined.","keywords":["boson stars","excited states","radial perturbations","radial stability","turning points","Einstein–Klein–Gordon system","self-interacting scalar field","eigenvalue problem"],"falsifier":"Take the $n=1$, $\\eta=28$ branch near the predicted crossing $\\omega\\approx 0.858$; perturb a configuration with $\\omega$ just above and just below that value and evolve it spherically. If the onset of instability does not occur at the crossing, or if an independent integration of the original singular perturbation equations puts the zero at a different $\\omega$, the central claim is falsified.","tokens_in":18056,"feed_emoji":"🌟","tokens_out":9011,"duration_ms":83736,"temperature":0.7,"pith_summary":"This paper establishes that for spherically symmetric boson stars, including excited configurations whose scalar profile has radial nodes, the lowest radial perturbation mode crosses zero at exactly the same point along the equilibrium branch where the ADM mass, Noether charge, and binding energy reach their first simultaneous extremum. To get there, the authors reformulate the linearized Einstein–Klein–Gordon pulsation equations in additive variables that remain finite at the nodes of the background field, so the eigenvalue problem can be integrated directly through nodeful backgrounds. They verify the coincidence for ground and excited branches up to seven nodes across several values of the quartic self-interaction strength, and they further show that the fundamental radial eigenvalue evaluated at the thresholds identified in nonlinear spherical evolutions obeys a simple power law in node number and self-interaction coupling. A sympathetic reader would care because the result connects a spectral stability diagnostic with a purely equilibrium turning-point criterion, making stability boundaries of boson stars easier to locate and opening a regular perturbative window onto excited configurations.","feed_headline":"For boson stars, radial instability starts exactly at turning point","feed_subtitle":"A node-regular method finds the first radial zero at the first equilibrium critical point.","key_machinery":"The load-bearing object is the pair of regular additive perturbation variables $f \\equiv \\varphi_0\\,\\delta\\varphi_1$ and $g \\equiv -\\varphi_0\\,\\xi/\\omega$, which replace the conventional relative variables $\\delta\\varphi_1$, $\\delta\\varphi_2$ that contain inverse powers of $\\varphi_0$. In these variables the linearized Einstein–Klein–Gordon system closes as a first-order system in $(f, f', g, g', \\mathcal{N})$ with no singular terms at the background nodes, so ordinary shooting can integrate through the nodes of any excited configuration. The eigenvalue $\\chi^2$ is then extracted with a two-parameter shooting ($\\chi^2$ and $\\zeta_1$) and outer boundary conditions on $f$, $\\delta q$, and $\\Lambda$, with a zero-mode diagnostic based on the rank of the $3\\times 2$ boundary map $\\mathcal{C}(\\chi^2; r_{\\mathrm{out}})$ and its singular-value ratio $\\epsilon$.","core_discovery":"On the paper's own terms, the central discovery is that the constrained fundamental radial eigenvalue $\\chi^2_{0,n}$ of an excited boson star branch crosses zero at the same background frequency as the first simultaneous extremum of the ADM mass, Noether charge, and binding energy. Previously the correspondence between radial mode zero and turning point was established only for the nodeless branch; the obstacle for excited states was that standard perturbation variables become singular at the nodes of the background scalar field. By working with regular additive variables $f$ and $g$, the authors obtain a closed first-order system with no inverse powers of the background field, allowing a shooting method to follow the lowest radial mode continuously along fixed-node branches. For every branch examined ($n=0$ through $n=7$, with quartic couplings $\\eta=0$, $7.9577$, $28$, $58$), the zero crossing coincides with the critical point within numerical resolution. Additionally, evaluating the radial eigenvalue at the self-interaction thresholds reported in nonlinear spherical evolutions yields a power-law scaling $\\chi^2_{0,n} \\propto \\eta\\, n^{p}$, with the same scaling holding for higher-node configurations not used in the fit.","pith_inferences":["If the zero-crossing/turning-point correspondence survives outside the tested parameter range, it supports a general stability-exchange principle for boson stars that might also apply to rotating or charged variants, where critical points are still defined by mass/charge extrema.","The empirical power law might be derivable from a WKB estimate of the effective potential between nodes; an analytic derivation would turn it into a predictive formula.","A direct next test is to compute the lowest nonradial mode on the same branches; if the radial zero remains the first mode to destabilize, the correspondence would extend to higher multipoles, with possible observational signatures in gravitational-wave echoes."],"forward_implications":["If the coincidence holds generally, the radial stability boundary for a fixed-node branch can be located from the equilibrium curves alone, without a separate eigenvalue solve.","The regular treatment of nodes extends standard radial-mode technology to nodeful solitonic backgrounds, where the old relative variables fail.","The empirical relation $\\chi^2_{0,n}/\\mu^2 \\simeq C\\, \\eta\\, n^{p}$ at the nonlinear thresholds can serve as a quick estimator for the coupling at which excited configurations become dynamically long lived, before running full evolutions.","The observed softening of the fundamental mode with increasing node number implies that higher excited branches are more easily destabilized, consistent with their transient role in dynamics."],"supporting_citations":[{"why":"Review of boson-star equilibrium branches and the first law $dM = \\omega\\, dQ$ used to define critical points.","marker":"[7]"},{"why":"Provides the standard radial pulsation equations in relative variables that the regular reformulation replaces.","marker":"[23]"},{"why":"Introduces the turning-point/mass-extremum stability diagnostic used as the comparison.","marker":"[43]"},{"why":"Establishes the nodeless correspondence between radial-mode zero and the mass maximum, the result being generalized.","marker":"[45]"},{"why":"Documents a case where the nodeless correspondence fails, motivating the explicit excited-branch check.","marker":"[46]"},{"why":"Reports the nonlinear spherical-evolution thresholds whose radial eigenvalues the paper fits.","marker":"[48]"},{"why":"Earlier radial stability criteria for excited boson stars that the current computation refines.","marker":"[49]"}],"fun_headline_variants":[],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result stands only if the reformulated $f,g$ perturbation equations together with the outer boundary conditions reproduce the true physical radial spectrum without introducing spurious or gauge zero modes, and if the rank condition truly distinguishes a zero mode from a near-zero mode.","fun_headline_variants_meta":{"error":"Client error '402 Payment Required' for url 'https://api.deepseek.com/chat/completions'\nFor more information check: https://developer.mozilla.org/en-US/docs/Web/HTTP/Status/402"},"cache_creation_input_tokens":0},"created_at":"2026-08-07T15:34:10.740311+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the $n=1$, $\\eta=28$ branch near the predicted crossing $\\omega\\approx 0.858$; perturb a configuration with $\\omega$ just above and just below that value and evolve it spherically. If the onset of instability does not occur at the crossing, or if an independent integration of the original singular perturbation equations puts the zero at a different $\\omega$, the central claim is falsified.","supporting_citations":[{"cited_title":"Gravitational Stability of Boson Stars","cited_arxiv_id":"0810.0696","evidence_quote":"Introduces the turning-point/mass-extremum stability diagnostic used as the comparison."}],"review_version":1}