{"id":"5b2bb33b-65c9-46c7-8f60-cca11fa14611","arxiv_id":"1908.09389","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"High-order WKB series summed with Padé approximants reproduce black hole quasinormal frequencies, including overtones, to up to 24 decimal places.","lead":"This paper pushes a WKB approximation for black hole ringdown frequencies to hundreds of terms and sums the series with Padé approximants, matching known numerical answers to many decimal places. It offers a fast, simple recipe for computing ringdown frequencies of non-rotating black holes without heavy numerical codes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Padé stabilization is shown only for diagonal approximants; no check that subdiagonal approximants converge to the same limit, leaving the wrong-analytic-continuation risk open.","rationale":"Agree with the reader that the weakest assumption is the unproven convergence of the Padé-summed WKB series. The paper's empirical evidence is strong: low-order Λ_k are verified against the Iyer-Will construction up to k=16, results match continued-fraction benchmarks for many modes, and the RN l=2,n=0 mode reproduces Hatsuda's value to 24 digits. Independent support includes the Borel-Padé cross-checks and the Wynn epsilon algorithm's connection to Padé approximants. However, the central claim is a method claim: it should work 'practically without any changes' for a wide class of metrics. For that claim, stabilization of one sequence is insufficient. The paper itself concedes that stabilization can be slow and that equality of consecutive results does not guarantee approach to the final value; the (0,0) mode episode shows that a stabilized diagonal sequence can disagree with a published benchmark, and the resolution required an external CF calculation. The proposed off-diagonal test directly probes whether the Padé table is convergent in the relevant sense; if diagonal and subdiagonal approximants disagree, the method's black-box use on new potentials is unsafe. Therefore the reader's CONDITIONAL verdict is appropriate; no adjustment.","tokens_in":14976,"tokens_out":13626,"duration_ms":133956,"concrete_test":"For a Schwarzschild mode without an independent high-precision benchmark, e.g. scalar l=0, n=2, generate the WKB-anharmonic series coefficients Λ_k for k up to 700 using the method described in §II, then compute three families of Padé approximants at ε=1: diagonal P_k^k, and subdiagonal P_{k-1}^k and P_k^{k-1}, for k = 150, 250, 350. If all three converge to the same complex value, the wrong-analytic-continuation concern is mitigated for that mode. If the subdiagonal families converge to a different value (or fail to converge), the diagonal stabilization is spurious and the method cannot be trusted for unbenchmarked modes.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim — that high-order Padé transforms of the WKB series reproduce known quasinormal frequencies to many digits — rests on the assumption that the diagonal Padé sequence P_k^k(1) converges to the physically correct value of ω². The formal WKB series (Eq. 3) is asymptotic and generically divergent; Padé summation of a divergent series has no general convergence theorem unless the series belongs to a special class such as Stieltjes. The paper supplies empirical stabilization for k up to 350 (Figs. 1–2) and table comparisons, but no argument that the diagonal sequence converges to the true value rather than to a different analytic continuation. That this is a real risk, not a pedantic one, is illustrated by the paper itself: for the scalar (0,0) mode the Padé result stabilizes to 16 decimal places yet differs from Andersson's published value; the authors resolve the discrepancy by arguing Andersson's value is wrong, based on their own continued-fraction code. Stabilization alone is thus not a correctness diagnostic. The paper does not compute the neighboring Padé approximants P_{k-1}^k or P_k^{k-1} (apart from low-order cases in Paper I); for series with a well-behaved Padé table these must converge to the same limit. Without such a check, the method's reliability for modes lacking independent benchmarks is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the WKB-based method for computing quasinormal modes introduced in the authors' previous work, Paper I. The idea is to generate high-order terms in the formal expansion for ω², up to order 700, using the equivalence with an anharmonic-oscillator perturbation problem and the Bender–Wu/Sulejmanpasic–Ünsal toolkit, and then to sum the resulting formal series with diagonal Padé approximants. The authors compare their results for scalar, electromagnetic, and gravitational perturbations of the Schwarzschild and Reissner–Nordström black holes against continued-fraction results, reporting agreement to up to 24 decimal places for some modes and correcting two values from Andersson's tables. They also discuss the use of the Wynn epsilon algorithm and the Borel–Padé method.","tokens_in":15235,"tokens_out":5647,"duration_ms":56193,"significance":"The strength of the paper is its systematic numerical demonstration: Tables II–XIV show that high-order diagonal Padé transforms of WKB series reproduce known continued-fraction quasinormal frequencies, including overtones, with no free parameters fitted to the target frequencies. The method is simple and black-box in structure, and the paper's cross-checks with Borel–Padé summation and continued fractions are valuable. The main gap is that the convergence of the Padé-summed series to the physically correct frequency is not established, and the paper itself flags a case where stabilization agrees with the authors' own code but disagrees with a published value. Because the method is proposed as a tool for modes that may lack independent benchmarks, this gap is load-bearing.","major_comments":[{"comment":"The stabilization of the diagonal sequence P_k^k is the central evidence for the method, but the paper reports no neighboring Padé approximants, such as entries with numerator and denominator degrees differing by one. The appendix explicitly notes that the Wynn algorithm yields only one family of approximants and that other entries of the Padé table must be constructed independently. Without checking that neighboring approximants converge to the same limit, the possibility that the diagonal sequence converges to a different analytic continuation of the divergent series remains open. Please add such a check for at least one low-overtone and one high-overtone case.","section":"Section III.A, Figs. 1–2 and Appendix Eq. (A.4)"},{"comment":"The resolution of the discrepancy with Andersson's published value rests on the authors' own continued-fraction code. The agreement between the Padé result and their continued-fraction result to 16 decimal places is impressive, but no implementation details, convergence criteria, or error estimates for that code are provided. To make the correction convincing to readers, an independent high-precision value from a published and independently implemented method, or a reproducible notebook with machine-precision settings, should be supplied.","section":"Section III.A, scalar (0,0) mode, Table II"},{"comment":"The black-box claim in the abstract and in the Final Remarks presupposes that the Padé-summed series converges to the true ω². The paper provides no argument from the structure of the series, such as Stieltjes-type properties, and no criterion for detecting when stabilization is trustworthy. Figs. 1–2 show that apparent stabilization can be slow and that low-order results are unreliable, and the paper admits that the behavior for highly damped modes is unknown. The manuscript should either narrow the claimed domain of applicability or provide a practical consistency test, such as agreement between neighboring Padé approximants and between Padé and Borel–Padé sums.","section":"Section II, Eq. (3) and Section IV"}],"minor_comments":[{"comment":"The word 'quasiormal' appears in the abstract and Introduction; it should be 'quasinormal'.","section":"Abstract and Introduction"},{"comment":"In Eq. (3), the summation index is i in one place but the term is written with ε^j and Λ_j; please make the index notation consistent.","section":"Section II, Eq. (3)"},{"comment":"In the left column for l = 0, n = 2, the entry '45907867− 1.08026685i' appears to be missing a leading '0'.","section":"Table IV"},{"comment":"The sentence 'Sulejmanpasic and and ¨Unsal' contains a duplicated 'and'.","section":"Section II"},{"comment":"'BenderWu' should be typeset as 'Bender–Wu' for consistency with the reference list.","section":"Footnote 1"},{"comment":"The abstract's claim that the l = 2, n = 0 gravitational mode agrees with the continued-fraction method to 24 decimal places should identify the comparison source; Table XIV shows agreement between Padé and Borel–Padé values, while the corresponding continued-fraction value is not listed there.","section":"Abstract and Section III.B"}],"recommendation":"major_revision","confidential_remarks":"This is a useful technical paper with strong empirical evidence. I think it is likely acceptable after the authors address the Padé-table consistency check and provide stronger substantiation for the corrected (0,0) mode value. The main risk is overclaiming black-box reliability without a convergence diagnostic."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the Matyjasek-Telecka paper. Short version: it is a solid incremental follow-up to their Paper I, and the headline claim holds up for the modes they test. The method—constructing the WKB expansion for ω² to order up to 700 via the anharmonic-oscillator mapping, then taking diagonal Padé approximants—reproduces continued-fraction results for Schwarzschild scalar, vector, and gravitational modes to 9–16 digits across dozens of cases, and for a Reissner-Nordström example it stabilizes to 24 digits. That is real evidence. They also document slow stabilization, which is genuinely useful, and they correct a couple of Andersson’s table values with their own continued-fraction code. No free parameters are fitted to the target frequencies; the WKB coefficients are independent of the QNM results. The circularity burden is low. The citation pattern is appropriate: Leaver, Konoplya, Hatsuda, Sulejmanpasic–Ünsal, Andersson are all there.\n\nThe main soft spot is the one the stress-test note flags: nothing in the paper establishes that the diagonal Padé sequence is converging to the right analytic continuation rather than to some other limit. The stabilization plots are pretty, but stabilization alone is not a correctness proof. The (0,0) scalar mode is the clearest warning: their Padé result stabilizes to 16 decimals but disagrees with Andersson; they resolve it by trusting their own CF code. That may be right, and the Borel-Padé agreement makes it plausible, but it is not an independent check. They never look at subdiagonal approximants P_{k-1}^k or P_k^{k-1}, which would be a cheap consistency test for Padé-table regularity. So for modes without an independent benchmark, the method’s reliability is unproven. That is a genuine limitation, but it is shared by most resummation-based numerical methods.\n\nOther soft spots are minor: the 24-digit agreement claimed in the abstract is not shown side by side with a CF value anywhere in the text; Table XIV shows only P and Borel-Padé. And the Mathematica notebooks are available only on request, meaning no archived code or data for independent verification.\n\nWho is this for? Anyone working on semianalytic QNM tools or testing parametrized black-hole ringdown. It deserves a serious referee: the benchmarks are extensive, the claims are specific, and the limitations are honestly discussed. The right referee will ask for the code, a subdiagonal check, and a side-by-side CF comparison for the headline number.\n\nRecommendation: send it to peer review, conditionally, with those additions requested.","headline":"A serious, well-benchmarked extension of the WKB-Padé program to order 700; the machinery is prior work, but the high-order results and corrections are new, and the paper deserves refereeing despite an unproven convergence claim and no archived code.","tokens_in":15743,"tokens_out":3682,"would_cite":true,"duration_ms":36177,"reading_group":"maybe","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":"High-order Padé summation matches black-hole frequencies to 24 digits.","keywords":["quasinormal modes","black hole perturbation theory","WKB approximation","Padé approximants","anharmonic oscillator","Schwarzschild black hole","Reissner-Nordström black hole","overtones"],"falsifier":"Compute the Schwarzschild scalar $l=0$, $n=0$ mode with an independent high-precision method—say a pseudospectral discretization of the perturbation equation—and compare with the Padé-stabilized value; a disagreement beyond the quoted 16 digits would show the stabilization is spurious. The same check applies to any $(l,n)$ pair: a single mode whose Padé plateau differs from a reliable independent calculation in more than one decimal place would falsify the universal claim.","tokens_in":14767,"feed_emoji":"🕳️","tokens_out":8183,"duration_ms":73778,"temperature":0.7,"pith_summary":"This paper claims that the WKB approach to black-hole quasinormal modes—long treated as only approximate—can be promoted to a highly accurate semianalytic method by diagonal Padé summation of a very long perturbation series for $\\omega^2$. The authors extend the standard WKB expansion to hundreds of terms (up to $k=700$) by converting the scattering problem into an anharmonic-oscillator bound-state problem, then build diagonal Padé approximants instead of summing the raw series. They show that for Schwarzschild and Reissner-Nordström black holes the resulting frequencies agree with continued-fraction numerical results to many decimal places, 24 for the gravitational $l=2$, $n=0$ mode, and that slow stabilization means low-order WKB results can mislead. The demonstration matters because it offers a black-box, potential-in/frequencies-out route to accurate quasinormal frequencies, including overtones, without a dedicated numerical eigenvalue solver.","feed_headline":"High-order Padé summation matches black-hole frequencies to 24 digits","feed_subtitle":"Diagonal Padé transforms of a 700-term WKB series reproduce continued-fraction quasinormal-mode results.","key_machinery":"The load-bearing object is the diagonal Padé approximant $P_k^k$ formed from the formal WKB series $\\omega^2 = V(x_0) + \\sum_j \\varepsilon^j \\tilde{\\Lambda}_j$, where the $\\tilde{\\Lambda}_j$ are WKB correction terms built from derivatives of the effective potential at its maximum. To supply the hundreds of terms needed, the authors use the reduction of the resonance problem to a one-dimensional anharmonic oscillator and generate high-order Rayleigh-Schrödinger corrections numerically for the chosen multipole and overtone numbers. The diagonal Padé transform replaces direct summing of the series, which is a poor strategy, and an epsilon-acceleration algorithm yields the same diagonal approximants at lower computational cost.","core_discovery":"The central discovery is that the formal WKB series for the squared quasinormal frequency, once computed to sufficiently high order, is a good input to Padé summation: the diagonal transforms $P_k^k$ stabilize to the true complex frequency, and do so even for overtones that standard low-order WKB handles poorly. For the Regge-Wheeler potential the calculations cover scalar, electromagnetic, and gravitational perturbations of Schwarzschild, and for Reissner-Nordström the results match both Borel-Padé and continued-fraction computations. The paper also finds that the longstanding benchmark value for the scalar $l=0$, $n=0$ mode is wrong: the Padé and continued-fraction codes agree on a corrected value to 16 decimal places. The convergence behavior is demonstrated empirically, not proved.","pith_inferences":["A natural testable extension is to push the same high-order Padé machinery toward the algebraically special modes the paper says WKB cannot reach; whether any plateau appears with even longer series would reveal whether that limitation is quantitative or structural.","The paper computes Borel transforms only for selected cases; a systematic comparison of Padé and Borel-Padé stabilization across modes could show whether the cheaper Padé-only workflow is sufficient whenever a plateau appears.","The anharmonic-oscillator equivalence suggests the resummation should transfer to other one-dimensional scattering problems with locally harmonic, asymptotically constant potentials, such as effective potentials in modified gravity or wormhole spacetimes, though that transfer is not demonstrated here."],"forward_implications":["If the stabilization is genuine, low-order WKB results cannot be trusted on their own; only high-order Padé sequences that show a clear plateau are reliable.","The method extends with little change to Reissner-Nordström and likely to other spherically symmetric potentials, giving a single black-box pipeline for accurate mode frequencies.","Accurate overtones up to at least $n=5$ for low $l$ become accessible to semianalytic WKB techniques, relaxing the usual $n \\lesssim l$ rule of thumb.","The corrected scalar $(0,0)$ and $(1,0)$ frequencies show that high-order Padé plus continued-fraction agreement can serve as a cross-check on published numerical values."],"supporting_citations":[{"why":"Predecessor paper that introduced Padé summation of the WKB series for ω²; the present work extends it.","marker":"[1]"},{"why":"Continued-fraction method used as the high-precision comparison and verification standard.","marker":"[3]"},{"why":"Supplied the underlying WKB matching technique and the Λ correction functions.","marker":"[13]"},{"why":"Reduced the scattering problem to an anharmonic-oscillator bound-state problem.","marker":"[19]"},{"why":"Borel summation of analytically continued Padé transforms; results reproduced for comparison.","marker":"[25]"},{"why":"Provided the efficient computer-algebra construction of high-order perturbative corrections.","marker":"[26]"},{"why":"Provided the original perturbative method for anharmonic oscillators on which the high-order computation rests.","marker":"[27]"},{"why":"Epsilon-acceleration algorithm that produces the diagonal Padé approximants efficiently.","marker":"[29]"},{"why":"Numerical quasinormal frequencies used as the baseline for the tables; one value is corrected here.","marker":"[32]"}],"fun_headline_variants":["Padé with 700-term WKB matches black-hole quasinormal modes to 24 digits","24-digit black-hole quasinormal modes from high-order Padé-WKB","WKB-Padé reproduces black-hole quasinormal frequencies to 24 decimals","Corrected benchmark: Padé-WKB black-hole quasinormal mode frequencies","High-order Padé-WKB stabilizes quasinormal modes to 24 digits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method assumes that the diagonal Padé transforms of the formal high-order WKB series converge to the true quasinormal frequency; the paper demonstrates this empirically for the cases it studies but does not prove it.","fun_headline_variants_meta":{"raw":{"variants":["Padé with 700-term WKB matches black-hole quasinormal modes to 24 digits","24-digit black-hole quasinormal modes from high-order Padé-WKB","WKB-Padé reproduces black-hole quasinormal frequencies to 24 decimals","Corrected benchmark: Padé-WKB black-hole quasinormal mode frequencies","High-order Padé-WKB stabilizes quasinormal modes to 24 digits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001036,"raw_usage":{"total_tokens":4364,"prompt_tokens":950,"completion_tokens":3414,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":566,"completion_tokens_details":{"reasoning_tokens":3306}},"tokens_in":566,"tokens_out":3414,"duration_ms":24420,"temperature":1.0,"reasoning_tokens":3306,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:13:37.776263+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the Schwarzschild scalar $l=0$, $n=0$ mode with an independent high-precision method—say a pseudospectral discretization of the perturbation equation—and compare with the Padé-stabilized value; a disagreement beyond the quoted 16 digits would show the stabilization is spurious. The same check applies to any $(l,n)$ pair: a single mode whose Padé plateau differs from a reliable independent calculation in more than one decimal place would falsify the universal claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reduced the scattering problem to an anharmonic-oscillator bound-state problem."},{"cited_title":"Wynn, Mathematical Tables and Other Aids to Computation 10, 91 (1956)","cited_arxiv_id":null,"evidence_quote":"Epsilon-acceleration algorithm that produces the diagonal Padé approximants efficiently."},{"cited_title":"Andersson, Proc","cited_arxiv_id":null,"evidence_quote":"Numerical quasinormal frequencies used as the baseline for the tables; one value is corrected here."}],"review_version":1}