{"id":"9c77498e-e196-42c7-ba56-f66c8f27a384","arxiv_id":"2601.04892","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Black-hole resonances migrate along complex-plane flow lines toward hard-wall attractors, and repellers near the unperturbed modes explain why perturbation theory fails so early.","lead":"This paper explains how tiny, localized changes to the potential around a black hole can dramatically reshuffle its 'ringing' frequencies, framing the motion of these frequencies in the complex plane as a smooth flow toward hard-wall points. A generalist might read it to understand why black-hole ringdown 'spectroscopy' is fragile to small environmental perturbations such as orbiting matter or dark-matter clumps.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Attractor classification in §II D is internally inconsistent: hard-wall attractors are double zeros of g with algebraic ~1/ε approach, not simple hyperbolic attractors.","rationale":"The reader's weakest assumption correctly identifies the fixed-point classification inconsistency in Sec. II D: attractors are double zeros of g, not simple zeros, so the exponential approach is invalid. My independent reading confirms this: the closed-form definitions of F in terms of u⁻ and u⁺ imply that a simple zero of either function makes F′ have a double pole and g a double zero. The paper's own statement that 'simple zeros of g define fixed points' and the subsequent linearization with g′(z̄) ≠ 0 is therefore internally inconsistent. The algebraic ∼1/ϵ approach does not undermine the global flow picture or the explanation of the elephant-flea effect via repellers, but it is a real gap in the central framework and should be corrected. The reader's CONDITIONAL verdict is appropriate; I see no reason to move it. The exactness of Eq. (10) and the flow ODE (18) check out, and the qualitative migration to hard-wall frequencies is robust. Thus I recommend UNCHANGED (still CONDITIONAL), with the requirement that the authors revise the attractor classification and the associated exponential rate claim.","tokens_in":14170,"tokens_out":3464,"duration_ms":38605,"concrete_test":"For the Nariai/Pöschl–Teller case, take a zero z̄ of u⁺_{λω}(x₀) using the closed form in Eq. (25). Compute F(ω) = W/(u⁻u⁺) from Eqs. (25)–(26) and expand F′(ω) about ω = z̄. Verify that F′(ω) ∼ A/(ω−z̄)² and thus g(ω) = −1/F′(ω) ∼ C(ω−z̄)², so g′(z̄)=0. Then numerically integrate Eq. (18) for the fundamental QNM with x₀ = 10/√27 and check that |ω(ϵ)−z̄| behaves as const/ϵ for large ϵ (or, if the paper's exponential claim is correct, as exp(−const·ϵ)).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central flow picture rests on the claim (Sec. II D) that the integral curves of dz/dϵ = g(z) terminate at 'attracting points' that are simple zeros of g with Re[g′(z̄)] ≠ 0, giving |δz| ∝ exp(Re[g′(z̄)]ϵ). This is not self-consistent with the paper's own identification of the attractors as zeros of u⁺(x₀) or u⁻(x₀). If u⁺(x₀) has a simple zero at z̄, then F = W₀/(u⁻u⁺) has a simple pole there, so F′ has a double pole, and hence g = −1/F′ has a double zero: g(z) ≈ C(z−z̄)², so g′(z̄) = 0. The stated criterion Re[g′(z̄)] ≠ 0 fails, and the exponential approach is replaced by an algebraic one: solving dz/dϵ = C(z−z̄)² gives z−z̄ ∼ −1/(Cϵ) as ϵ → ∞ (for the attractor branch). The endpoints are still approached, so the topological conclusion 'resonances migrate to hard-wall frequencies' survives, but the dynamical-systems classification as hyperbolic attractors is quantitatively wrong. This matters because the paper uses the exponential law to characterize the approach and to define the 'attracting points' in the phase portrait; the double-zero degeneracy changes the local normal form and the rate of convergence. The elephant-flea mechanism (repellers near unperturbed resonances) is not affected, but the central framework's statement that the flow has simple zeros as attractors requires revision.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies how the quasinormal-mode and Regge-pole spectra of black-hole-like potentials deform under a delta-function perturbation of the potential. The central exact result is Eq. (10): for a perturbation ϵδ(x−x₀), the resonance condition reduces to F(λ,ω)+ϵ=0, with F constructed from the unperturbed radial functions and Wronskian. On this basis the authors introduce a flow ODE dz/dϵ=−1/F′(z), interpret resonance migration as a dynamical system with attracting and repelling points, and use it to explain the 'elephant and flea' spectral instability. The framework is applied first to the Nariai/Pöschl–Teller case, where closed forms allow a detailed analysis of attractors, repellers, and linear/nonlinear instability thresholds, and then to Schwarzschild, where numerical integration is used for the fundamental QNM and for Regge poles.","tokens_in":14346,"tokens_out":13933,"duration_ms":145280,"significance":"The exact resonance condition (10) is a valuable and nontrivial result, and the Nariai analysis provides a controlled laboratory for spectral instability. The paper is careful to separate linear, nonlinear, and anomalous instability and connects the linear coefficient to QNM excitation factors. The Schwarzschild Regge-pole results are cross-checked against an independent continued-fraction calculation. If the dynamical-systems picture is corrected as described below, the attractor–repeller mechanism would provide a clear, parameter-free explanation of why weak localized perturbations strongly destabilize high overtones. These strengths make the paper potentially suitable for publication after revision.","major_comments":[{"comment":"The classification of fixed points is internally inconsistent. The text states that simple zeros of g(z) define fixed points, that these correspond to simple poles of F′(z), and that they are typically associated with simple zeros of u⁺(x₀) or u⁻(x₀), with an exponential approach |δz|∝exp(Re[g′(z̄)]ε). But if u⁺(x₀) has a simple zero at z̄, then F=W₀/(u⁻u⁺) has a simple pole, F′ has a double pole, and g=−1/F′ has a double zero with g′(z̄)=0. The stated criterion Re[g′(z̄)]≠0 never applies to these hard-wall attractors. The correct local behavior is algebraic: solving dz/dε=C(z−z̄)² gives z−z̄∼−1/(Cε). The topological conclusion that resonances tend to the hard-wall frequencies survives, but the phase-portrait classification and all statements relying on exponential approach need revision. Please correct Sec. II D and adjust the associated discussion in Secs. III and V.","section":"Sec. II D, Eq. (18)"},{"comment":"The Regge-pole asymptotic is claimed to explain the power-law decay of ϵ_lin with n+1/2 observed in Fig. 7. However, substituting λ_n=ω+i(n+1/2) into Eq. (34) gives |λ_n^{2iω}|=e^{-2ω atan((n+1/2)/ω)}=e^{-πω}+O(1/n), which is n-independent at leading order. Thus the leading term in Eq. (34) does not produce a power-law decay. Either Eq. (33) is missing a factor from u⁻u⁺ (or from the gamma-function ratio), or the power-law claim requires a different derivation. Please supply the missing asymptotic steps or revise the claim.","section":"Sec. III.C, Eqs. (33)-(34)"}],"minor_comments":[{"comment":"Typos: 'refered' should be 'referred' (Introduction); 'Scwharzchild' in Fig. 9 caption; 'signification' should be 'significant' in Sec. IV.","section":"Throughout"},{"comment":"Notation is inconsistent between hatted quantities in Eq. (24) and unhatted ω,x in the surrounding text. Please state once the association x↔νx̂, ω↔ω̂/ν and use it consistently.","section":"Sec. III"},{"comment":"The sentence 'which admits ω=0 as a solution' is unclear in context. Is this a special case of the w-mode condition, or a spurious root? Please clarify.","section":"Sec. III.A, Eq. (29)"},{"comment":"The 90° branch-switching argument assumes F″(z_r)≠0. If F′ has a higher-order zero, the local normal form changes. A brief comment on this genericity assumption would be useful.","section":"Sec. II.D, Eq. (19)"}],"recommendation":"major_revision","confidential_remarks":"The exact resonance condition and the Nariai analysis are strong contributions, and the numerical Schwarzschild results appear reliable. The main technical flaw is the attractor classification in Sec. II D, which is fixable without changing the overall picture. The Regge-pole asymptotic inconsistency in Sec. III.C should be resolved carefully; it may indicate a missing factor rather than a false numerical observation. I recommend major revision, with the expectation that the corrected manuscript will be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Torres and Dolan have a clean exact result: for a delta-function potential perturbation, the resonance condition is exactly F(ω,λ)+ε=0 with F = W0/(u−u+). I checked the matching, and it's exact—no O(ε²) corrections. That gives a way to compute the entire deformed QNM and RP spectrum from unperturbed data, which is new and worth having. The flow-ODE reformulation is a nice way to see why the Taylor series breaks down: repellers sit near the unperturbed resonances for large x0, and trajectories bend around them. The Nariai/Pöschl–Teller analysis is mostly closed-form, and the Schwarzschild RP numerics are cross-checked against a continued-fraction code. I think this is a genuine contribution to the spectral-instability literature.\n\nBut there is a real mathematical gap in the attractor classification (Sec. II D). The paper calls the attractors simple zeros of g with Re[g′]≠0 and predicts exponential approach. That is not consistent with the paper's own identification of attractors as zeros of u±(x0). A simple zero of u± gives F a simple pole, so F′ has a double pole and g has a double zero. So g′(z̄)=0, and the approach is algebraic, ~1/ε, not exponential. The endpoints survive, and the elephant-flea mechanism is untouched, but the local normal form and the rate of convergence are wrong. That needs a correction before publication. It is not fatal to the main qualitative conclusions, but it is a claim the paper makes explicitly and it is false as written.\n\nSecond, Eq. (34) is supposed to explain the power-law decay of ϵ_lin for RPs, but the displayed leading term has modulus e^{−2ω arg λ_n}, which saturates to a constant as n→∞; it does not obviously give the clean power law shown in Fig. 7. Maybe the exact Gamma expression does over the plotted range, but the sentence as written overclaims. That should be re-derived or softened.\n\nMinor: no code or data shipped, and the Schwarzschild QNM part is limited to the fundamental mode—the abstract sounds more general than that, though the text is honest about the limitation.\n\nBottom line: the core exact condition and the migration picture are solid and useful. The attractor classification needs a serious fix, and the RP-threshold explanation needs a check. I would send this to a good referee. The authors are clearly capable of fixing these points.","headline":"A genuinely useful exact result on delta-perturbed black-hole spectra, but the dynamical-systems attractor classification is wrong in a way that needs fixing.","tokens_in":15075,"tokens_out":4647,"would_cite":true,"duration_ms":48229,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.70.-s","04.30.-w"],"model":"deepseek-v4-flash","headline":"This paper claims that adding a delta-function spike to a black-hole potential sends every quasinormal-mode and Regge-pole resonance along a smooth trajectory, and that the entire deformed spectrum follows from an exact condition built only","keywords":["quasinormal modes","Regge poles","spectral instability","delta-function perturbation","black hole perturbation theory","dynamical systems","attractor-repeller flow","elephant and flea phenomenon"],"falsifier":"Compute the exact perturbed quasinormal-mode frequencies of the Pöschl-Teller model by root-finding Eq. (10) at large ε and measuring the rate at which each mode approaches its attractor. If the deviation from the attractor decays as 1/ε rather than exponentially—as the double zero of g at a simple zero of u⁻ or u⁺ implies—then the exponential limit law in Sec. II D needs revision. A simpler check: verify symbolically that Eq. (9) contains no O(ε²) correction.","tokens_in":13849,"feed_emoji":"🕳️","tokens_out":7868,"duration_ms":69903,"temperature":0.7,"pith_summary":"The paper aims to show that a delta-function perturbation added to the radial potential of a black hole does not destroy the resonance spectrum but re-routes it: every quasinormal mode and Regge pole moves continuously along a curve in the complex frequency or angular-momentum plane as the perturbation strength grows. The key claim is that this curve is governed by an exact, non-perturbative resonance condition involving only unperturbed wavefunctions, and that the global migration is organized by a small set of attracting and repelling points. If true, the notorious 'elephant and flea' instability—where a tiny distant perturbation reshuffles the overtones—is not a failure of the spectrum's existence but a local failure of Taylor-series perturbation theory, driven by repellers sitting near the unperturbed resonances. A sympathetic reader would care because this turns a seemingly erratic spectral instability into a predictable flow, and explains why quasinormal-mode and Regge-pole overtones are sensitive in different ways.","feed_headline":"Delta spike makes black-hole resonances flow along exact curves","feed_subtitle":"Quasinormal modes and Regge poles slide to predictable attractors; the elephant-flea effect is a repeller artifact.","key_machinery":"The object doing the work is the meromorphic function F(λ,ω)=W(u⁻,u⁺)/(u⁻(x₀)u⁺(x₀)) built from the unperturbed left- and right-valid radial solutions. The resonance condition F+ε=0 converts spectral deformation into a first-order autonomous flow dz/dε=−1/F′(z); zeros of u⁻ or u⁺ at x₀ are attracting fixed points (hard-wall limits), while simple poles of F′ act as repelling points or junction points that bend trajectories by 90°, and the positions of these features move with x₀.","core_discovery":"The central discovery is Eq. (10): for a perturbation ε δ(x−x₀), the resonance condition is exactly F(λ,ω)+ε=0, where F is the ratio of the unperturbed Wronskian to the product of the two unperturbed radial solutions evaluated at x₀. No O(ε²) terms appear, so the deformed spectrum at any strength is computable from unperturbed quantities. From this, resonances migrate along integral curves of dz/dε = −1/F′(z), with fixed points of the flow attracting modes toward 'hard-wall' frequencies where one radial function vanishes at x₀, and repelling points near the unperturbed resonances explaining the breakdown of linear and quadratic approximations.","pith_inferences":["A direct consequence left implicit in the paper: the same exact reduction should hold for any compactly supported perturbation by decomposing it into delta spikes, suggesting a rational-function-in-ε form for the deformed resonance condition in more general settings.","The paper's numerical Schwarzschild analysis is restricted to the fundamental quasinormal mode; extending the flow integration to higher overtones would test whether the repeller picture survives quantitatively beyond the Nariai model.","At an attractor that is a simple zero of u⁻ or u⁺, the stated exponential approach cannot hold because g′ vanishes there; the approach should be algebraic (δz ∝ 1/ε). Measuring this rate numerically would sharpen the dynamical-system classification.","The flow picture suggests a practical numerical recipe: instead of root-finding the perturbed problem at each ε, integrate dz/dε = −1/F′(z) once, using unperturbed data, to map the entire spectral deformation."],"forward_implications":["For a delta spike, the full deformed spectrum—not just a first-order shift—is determined by unperturbed wavefunctions; no higher-order matching is needed.","Resonances never appear or disappear abruptly; as ε grows they slide continuously toward hard-wall frequencies, so the spectrum is globally smooth even when overtone labels swap.","The elephant-and-flea effect is explained as repellers sitting close to unperturbed overtones for distant perturbations, making the Taylor series valid only for extremely small ε.","Quasinormal-mode overtones destabilize exponentially with distance and overtone order, while Regge-pole overtones destabilize as a power law, as quantified by the threshold parameters ε_lin and ε_nonlin.","The same attractor–repeller skeleton organizes both the exactly solvable comparison model and the Schwarzschild case, suggesting the mechanism is generic rather than potential-specific."],"fun_headline_variants":["Black-hole spectra flow exactly under delta spikes","Delta perturbations bend quasinormal modes predictably","Resonance drift mapped via exact delta-function formula","Black-hole modes show repeller-driven spectral instability"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The flow picture stands on the assumption that F(z)=W/(u⁻u⁺) is meromorphic with only isolated poles and exactly one repelling point per unperturbed overtone, with no branch cuts crossing the trajectories; this assumption is nontrivial, because at a simple zero of u⁻ or u⁺ the fixed point is not simple and the stated exponential approach becomes algebraic.","fun_headline_variants_meta":{"raw":{"variants":["Black-hole spectra flow exactly under delta spikes","Delta perturbations bend quasinormal modes predictably","Resonance drift mapped via exact delta-function formula","Black-hole modes show repeller-driven spectral instability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000125,"raw_usage":{"total_tokens":976,"prompt_tokens":807,"completion_tokens":169,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":551,"completion_tokens_details":{"reasoning_tokens":119}},"tokens_in":551,"tokens_out":169,"duration_ms":2844,"temperature":1.0,"reasoning_tokens":119,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T11:54:52.758375+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact perturbed quasinormal-mode frequencies of the Pöschl-Teller model by root-finding Eq. (10) at large ε and measuring the rate at which each mode approaches its attractor. If the deviation from the attractor decays as 1/ε rather than exponentially—as the double zero of g at a simple zero of u⁻ or u⁺ implies—then the exponential limit law in Sec. II D needs revision. A simpler check: verify symbolically that Eq. (9) contains no O(ε²) correction.","supporting_citations":[],"review_version":1}