{"id":"ac2477ae-e446-4e69-b884-4b16f42ac4a0","arxiv_id":"2607.19492","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A bilinear-form framework computes black-hole quasinormal-mode frequency shifts to any order, but the mode-sum expansion of the first-order mode shift diverges and needs a continuum piece.","lead":"This paper develops a perturbation-theory framework for computing how black-hole quasinormal-mode frequencies shift when gravity is modified or the black hole is embedded in an environment. It gives formulas for shifts to arbitrary order and shows that the natural mode-sum expansion of the first-order mode shift diverges, so a continuous-spectrum piece is needed.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (40)'s 'arbitrary order' claim rests on the unproven structural ansatz Eq. (29): that every higher-order modified Teukolsky source is a ζ-independent linear operator acting on lower-order field shifts. No concrete theory is shown to satisfy this beyond first order, so the central claim is not ye","rationale":"The reader's weakest-assumption analysis identifies the same load-bearing point: Eq. (29) is introduced as an expectation, not a proven consequence, and Eq. (40) inherits that assumption. My read agrees with the conditional verdict. I would not strengthen to rejection because the first- and second-order machinery is explicit, the bilinear-form orthogonality is independently grounded in Ref. [30], and the numerical checks are internally consistent. However, the paper's central novelty is the 'to any order' claim, and that claim is only as secure as the structural ansatz. The manuscript itself flags the relevant limitations: Sec. III A says the higher-order structure is 'expected', sets aside mirror-mode/degenerate complications, and the second-order validations are consistency checks using known mode shifts rather than closed predictions from Eq. (29). A concrete third-order derivation in a real modified-gravity model would settle whether the ansatz actually holds. Until then, the proper verdict is conditional: the framework is promising and well-motivated, but the arbitrary-order generalization is not yet established.","tokens_in":28860,"tokens_out":9831,"duration_ms":105515,"concrete_test":"Derive the third-order sourced Teukolsky equation for a concrete beyond-GR theory continuously connected to GR, e.g., Einstein-scalar-Gauss-Bonnet or dynamical Chern-Simons, using full metric reconstruction (CCK at first order, GHZ at second order). Write the k=3 source explicitly and check whether it is exactly -3S^(3)ψ^(0) - 6S^(2)ψ^(1) - 3S^(1)ψ^(2), with each S^(i) a linear ζ-independent operator and no additional terms depending on ω^(1), ω^(2), or mixed mirror modes. If extra terms appear, recompute ω^(3)_n directly from the expanded field equation and compare with Eq. (40); disagreement would falsify the arbitrary-order claim. If the form holds through third order, the structural ansatz gains concrete support.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central assertion is that QNM frequency shifts are computable to any order via Eq. (40). That formula follows only if the k-th order modified Teukolsky equation has the exact binomial form of Eq. (29), with each S^(i) a linear, ζ-independent operator acting on ψ^(j). The text introduces Eq. (29) explicitly as 'expected' rather than derived: Sec. III A states that for a generic k-th order modified Teukolsky equation the form 'is then expected', and it separately sets aside the known complications of complex conjugation, mirror-mode coupling, and degenerate perturbation theory. No construction of S^(2) or S^(3) is given for any concrete beyond-GR or environmental model. If realistic higher-order sources contain terms not of this form—e.g., frequency-dependent reconstruction kernels, mixed mirror-mode operators, or effective nonlinearities entering through the sourced metric reconstruction—Eq. (40) omits contributions and the 'any order' claim is unsupported. The first- and second-order checks in Secs. V and VI do not test this assumption: they insert known analytic or Leaver-based mode shifts, so they never exercise Eq. (29) at order k ≥ 3. The paper's own caveats therefore leave the arbitrary-order step as the weakest load-bearing link.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The article develops a perturbation-theory framework for quasinormal-mode (QNM) frequency shifts of the Teukolsky equation, based on a conserved bilinear form under which Kerr QNMs are orthogonal. It derives the first-order frequency shift (36), the second-order shift (38), and a general k-th-order formula (40) in terms of lower-order mode shifts and linear source operators S^(j). It also constructs a formal spectral decomposition of the first-order mode shift, Eq. (60), into QNM projections plus a continuum contribution, and presents numerical evidence—in Pöschl–Teller and slowly spinning Kerr toy models—that the QNM-only sum diverges. The first- and second-order formulas are benchmarked against independent analytic and numerical results with no fitted parameters.","tokens_in":29177,"tokens_out":12071,"duration_ms":117328,"significance":"If the general formula is established, it closes an important gap: no systematic higher-order QNM perturbation theory currently exists. The paper has concrete strengths: the bilinear-form derivation of the first-order formula is explicit and reproduces the known eigenvalue-perturbation result; the second-order formula is checked against the exact Pöschl–Teller shift and against numerical Kerr data; and the divergence of the QNM sum is a sharp, testable manifestation of QNM incompleteness. The central caveat is that the advertised \"any order\" result rests on an unproved structural ansatz for the source operators, and the paper's own caveats limit the scope of the theorem.","major_comments":[{"comment":"The \"any order\" formula (40) is derived from Eq. (29), but Eq. (29) is introduced as \"expected\" rather than derived. The text assumes each S^(j) is a linear, ζ-independent operator acting on ψ^(i), and no concrete construction of S^(2) or S^(3) is given for the modified-gravity/environment class discussed in the Introduction. The checks in Secs. V and VI use known full solutions to produce ψ^(1) and verify only k = 1, 2; they never exercise the binomial ansatz at k ≥ 3. Thus the central claim \"frequency shifts to any order\" is conditional on an unproved structural hypothesis. The authors should either prove Eq. (29) for a defined class of theories (e.g., when the modified equation is O†(ζ)ψ(ζ)=0, in which case S^(j)=O†(j) follows) or explicitly state the theorem with Eq. (29) as a hypothesis and provide an example that satisfies it at order 3.","section":"Sec. III A, Eq. (29)"},{"comment":"In Eq. (57) the constant α ∈ C is left arbitrary. Replacing the sum over n′ ≠ n in Eq. (59) by the full sum over n′ in Eq. (60) subtracts the n′ = n term ⟨⟨ψ^(0)_n, ψ^(1)_n⟩⟩/⟨⟨ψ^(0)_n, ψ^(0)_n⟩⟩ ψ^(0)_n = (−iω^(1)_n t + α)ψ^(0)_n. After adding the explicit secular term −iω^(1)_n t ψ^(0)_n, an αψ^(0)_n term remains. Equation (60) is therefore an identity only after one sets α = 0, or explicitly absorbs αψ^(0)_n into the normalization of ψ^(1)_n. As written, the displayed spectral decomposition is not correct; this should be fixed and stated.","section":"Eq. (60)"},{"comment":"The framework is developed under the assumption that S^(1) is a linear operator acting on ψ^(0)_n alone, explicitly setting aside complex conjugation, mirror-mode coupling, and degenerate perturbation theory. These are not exotic: for real metric perturbations S^(1) maps ψ^(0)_lmnp to ψ^(0)_l,−m,n with frequency −ω* and the first-order problem is degenerate. Thus Eq. (40), and even Eq. (36) as presented, do not cover all of the physical applications cited in the Introduction (e.g., bGR theories with real metric perturbations). The authors should state this scope limitation in the abstract/conclusions or extend the derivation to the antilinear/mirror-mode case following Ref. [18].","section":"Sec. III A, after Eq. (27)"}],"minor_comments":[{"comment":"The caption says the true mode shift \"only appears constant due to the plot scale\"; please quantify the true value and the rate of divergence of the partial sums, so the reader can see the contrast quantitatively.","section":"Fig. 3"},{"comment":"The notation Jt is introduced as \"simply reduces to Jt: t→−t,\" but the spin-weight action of J in the 1+1 toy model should be clarified, since in the main text J includes the GHP prime operation.","section":"Sec. V B"},{"comment":"There are several typographical artifacts: \"P¨ oschl-Teller\" with a misplaced dieresis, \"l’Hˆ opital\" in App. B, and inconsistent use of \"Schr¨ odinger\". These should be cleaned up.","section":"Throughout"},{"comment":"The equivalence to the Sturm-Liouville product is stated in Eq. (21) and then used heavily; it would help to explicitly mention that the limit ω2→ω1 reproduces Eq. (19), since this is used in deriving Eq. (36).","section":"Sec. II E"}],"recommendation":"major_revision","confidential_remarks":"This is a solid formal contribution, but the advertised \"any order\" generality exceeds what is currently proved. The key revision is to turn Eq. (29) into a theorem with explicit hypotheses, or to clearly state the result as conditional on that ansatz and provide at least one nontrivial order-3 example. Eq. (60) also has a small but real missing α-term that should be fixed by choosing a normalization. I see no grounds for rejection, but the central claim needs tightening."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful core of this paper is real. It takes the bilinear-form orthogonality result from the authors' earlier work and builds a perturbation theory for Teukolsky QNM frequency shifts beyond first order. The genuinely new pieces are the k-th order formula, Eq. (40), the spectral decomposition of the first-order mode shift, Eq. (60), and the observation that the QNM-only sum diverges. First order reduces cleanly to the known eigenvalue-perturbation result, which is the right sanity check. The Pöschl-Teller and slowly-spinning Kerr validations compare against independent analytic or Leaver-based results, so they are not circular. The paper also deserves credit for being explicit that the spectral decomposition is formal and that the practical route is to solve for the mode shift numerically or via the regularized Green's function in App. B. Self-citation to the bilinear-form theorem is legitimate here; that theorem is the foundation, and it is a published, checkable result.\n\nThe main soft spot is exactly the one the stress-test note identifies: the \"arbitrary order\" claim rests on Eq. (29), the assumed binomial form of higher-order modified Teukolsky equations with linear, frequency-independent source operators. The text says \"expected,\" not derived, and it explicitly sets aside complex conjugation, mirror-mode coupling, and degenerate perturbation theory. No concrete beyond-GR or environmental theory is shown to satisfy this structure at order three or higher, and the numerical checks never exercise that assumption. That is a load-bearing caveat, but it is a caveat the authors themselves flag. It does not undermine the first- and second-order results, which are explicit, checked, and already useful. The divergence of the QNM sum is also a genuine finding; it makes incompleteness concrete and should catalyze work on regularization.\n\nI would not call the central argument flawed. I would call the \"any order\" headline stronger than the evidence. In the published version that claim should be framed as conditional on the structural ansatz of Eq. (29), with a clear statement that verifying that ansatz for specific modified theories is open work. The mode-shift spectral decomposition likewise should stay framed as formal, as it currently is.\n\nWho gets value from this: anyone computing ringdown spectral shifts in modified gravity or environmental scenarios, and anyone working on QNM completeness. It deserves a serious referee. My recommendation is to send it to peer review, with a request that the authors either prove or visibly soften the arbitrary-order claim.","headline":"Genuinely useful derivation of higher-order QNM frequency shifts with honest caveats; the \"any order\" claim rests on an explicit but unproven structural assumption, so it deserves a serious referee rather than a desk reject.","tokens_in":29667,"tokens_out":1587,"would_cite":true,"duration_ms":21312,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","83C35"],"pacs":["04.30.-w","04.70.-s"],"model":"deepseek-v4-flash","headline":"A bilinear form lifts Schrödinger perturbation theory to black hole quasinormal modes, yielding frequency shifts to arbitrary order.","keywords":["quasinormal modes","black hole spectroscopy","Teukolsky equation","bilinear form","perturbation theory","ringdown","modified gravity","spectral decomposition"],"falsifier":"Compute the third-order frequency shift for the Pöschl-Teller width perturbation by differentiating the exact frequency ω(ζ) and compare it with the third-order version of the paper's formula using numerically obtained mode shifts; a disagreement would falsify the claim that the formula holds to arbitrary order.","tokens_in":28725,"feed_emoji":"🕳️","tokens_out":5794,"duration_ms":59604,"temperature":0.7,"pith_summary":"This paper aims to give black-hole spectroscopy the same tool that quantum mechanics has had for a century: a systematic perturbative expansion for how quasinormal-mode frequencies shift when physics departs from vacuum general relativity. Using a symplectic-current bilinear form under which quasinormal modes are orthogonal—despite being incomplete as a basis—the authors derive a formula for the frequency shift at any order in a small coupling, expressed in terms of lower-order mode shifts. They also derive a formal spectral decomposition of the first-order mode shift into projections onto unperturbed modes plus a continuous-spectrum integral, and show in two examples that the mode-only sum diverges. If correct, this removes the main obstacle to computing second- and higher-order ringdown frequency shifts for modified-gravity or environmental perturbations, which will matter as gravitational-wave observations sharpen.","feed_headline":"Black hole ringdown shifts now computable to any order","feed_subtitle":"A bilinear form tames quasinormal-mode incompleteness, enabling higher-order frequency corrections for precision gravity tests.","key_machinery":"The central object is the bilinear form built from the Teukolsky symplectic current with t–φ reflection instead of complex conjugation, regularized on a complex radial contour. It makes quasinormal modes orthogonal, and its norm coincides with the derivative of the Wronskian. Combined with a modified Teukolsky equation of assumed binomial structure, the balance law yields the arbitrary-order frequency-shift formula. The mode-shift decomposition uses the Teukolsky Green's function and contour deformation into QNM poles, a branch cut, and a high-frequency arc, making the continuous spectrum explicit.","core_discovery":"The paper's central claim is that quasinormal-mode frequency shifts under small deformations of the Teukolsky equation can be computed to arbitrary order. The key formula gives the k-th-order frequency shift as a combination of matrix elements of source operators acting on lower-order mode shifts and of bilinear-form projections, so the hierarchy closes order by order. At first order it reproduces the known eigenvalue-perturbation formula; at second order it reproduces the exact Pöschl-Teller shift and the slow-spin Kerr shift. The paper also establishes that the first-order mode shift admits only a formal spectral decomposition: a sum over unperturbed QNMs plus a continuum piece, where the","pith_inferences":["If QNM-sum divergence is generic, then regularization schemes that fold the branch-cut and high-frequency-arc continuum into the mode-shift expansion could turn the formal spectral decomposition into a practical computational tool; the paper hints at but does not demonstrate this.","The same bilinear-form perturbation theory should apply to boson-cloud quasibound states, where the nonrelativistic hydrogen analogy suggests the discrete sum may converge—a concrete next target.","The arbitrary-order claim could be stress-tested by deriving explicit source operators for a specific higher-derivative gravity and checking whether the binomial structure survives metric reconstruction; if it does not, the all-orders formula needs modification.","The time-dependent n' = n projection and secular term suggest that even the notion of 'mode shift' may need a freedom-fixing prescription at higher orders, since the second-order frequency shift was proven invariant under adding a multiple of the background mode."],"forward_implications":["Second- and higher-order frequency shifts are now in principle computable once the first-order mode shift is obtained by solving the modified Teukolsky equation numerically or via the regularized Green's function.","The first-order frequency-shift formula coincides with the established eigenvalue-perturbation result, placing the new framework on known footing.","The QNM-only mode-shift expansion diverges in both worked examples, so the continuous-spectrum contribution is not a technicality; practical mode-shift computations must include it or use a representation that avoids the expansion.","The framework covers modified theories and environmental effects that can be cast as modified Teukolsky equations with linear source operators, including axisymmetry-breaking perturbations that mix different m-modes.","The degenerate case, where metric reconstruction mixes modes through complex conjugation, is left open but expected to generalize following the established treatment."],"fun_headline_variants":["QNM shifts to any order via bilinear-form perturbation","Arbitrary-order black hole ringdown shifts from bilinear form","Bilinear form computes black hole QNM shifts to any order","From Schrödinger to black holes: arbitrary-order QNM shifts"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The arbitrary-order formula rests on the unproven assumption that the k-th-order modified Teukolsky equation always takes the binomial form with linear, frequency-independent source operators acting on lower-order fields; if a real modification produces nonlinear or frequency-dependent sources, the all-orders formula does not follow.","fun_headline_variants_meta":{"raw":{"variants":["QNM shifts to any order via bilinear-form perturbation","Arbitrary-order black hole ringdown shifts from bilinear form","Bilinear form computes black hole QNM shifts to any order","From Schrödinger to black holes: arbitrary-order QNM shifts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000814,"raw_usage":{"total_tokens":3405,"prompt_tokens":747,"completion_tokens":2658,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":491,"completion_tokens_details":{"reasoning_tokens":2587}},"tokens_in":491,"tokens_out":2658,"duration_ms":18775,"temperature":1.0,"reasoning_tokens":2587,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T12:34:30.135842+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the third-order frequency shift for the Pöschl-Teller width perturbation by differentiating the exact frequency ω(ζ) and compare it with the third-order version of the paper's formula using numerically obtained mode shifts; a disagreement would falsify the claim that the formula holds to arbitrary order.","supporting_citations":[],"review_version":1}