{"id":"08cd1d99-4471-483d-beac-230f29a4bfad","arxiv_id":"2609.14979","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":5,"one_line_summary":"A bound state method yields accurate low-lying quasinormal mode frequencies for Reissner-Nordström black holes, with the same procedure working for both nonextremal and extremal cases.","lead":"This paper applies a recently proposed bound state method to compute quasinormal mode frequencies of Reissner-Nordström black holes. The method simplifies calculations for extremal black holes, where standard techniques require special treatment.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek/deepseek-v4-flash","headline":"No significant objection identified: the central claim for low-lying modes (n≤2) is well-supported by validation with Leaver's method; limitations for higher overtones are honestly discussed.","rationale":"The paper's central claim is empirically verified with high-accuracy comparisons to Leaver's continued fraction method for the stated mode range (n≤2). The main limitation—degradation of accuracy for higher overtones and reliance on a tuned homotopy parameter—is openly discussed and does not contradict the headline results. The reader's weakest_assumption correctly identifies the analytic continuation as the most delicate step, but the validation within the claimed operating window is sufficient. Therefore, no revision to the ACCEPT verdict is warranted.","tokens_in":12382,"tokens_out":2858,"duration_ms":29963,"concrete_test":"Compute extremal RNBH QNM frequencies for gravitational ℓ=2, n=3,4 using homotopy λ=4 with a denser α-grid (61 points, α ∈ [0.1, 0.9]) and compare to a reference method (e.g., time-domain evolution or a matrix method) to test whether the bound state method's reliability extends beyond n=2 for the extremal case.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that the bound state method yields QNM frequencies with relative errors below 10⁻⁵ for n≤2 for Reissner-Nordström black holes, and that the extremal case requires only a change of the tortoise coordinate—is robustly supported by the validation in Tables I–IV and VII. The reader's weakest_assumption (accuracy of Padé continuation from a finite α-grid) is a genuine limitation, but it is openly acknowledged and does not undermine the validated results for the claimed mode range. The homotopy deformation is an ad hoc workaround, but the paper does not claim universal applicability beyond the demonstrated cases. Overall, the evidence matches the stated claims, and no additional load-bearing concern emerges that would alter the verdict.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek/deepseek-v4-flash","summary":"The paper applies a recently proposed bound state method to compute quasinormal mode (QNM) frequencies for gravitational and electromagnetic perturbations of Reissner-Nordström black holes (RNBHs).  By a two-step coordinate transformation and introduction of a real scaling parameter α, the QNM problem is mapped to a bound state eigenvalue problem.  The bound state energies are computed numerically on a grid of real α values and analytically continued to α = −i via Padé approximants to obtain the QNM frequencies.  The method is validated against Leaver's continued fraction method for nonextremal RNBHs and against results of Onozawa et al. for the extremal case.  For low overtones (n ≤ 2) relative errors are typically below 10⁻⁵; for n = 3 errors are at the ~10⁻² to 10⁻³ level.  A homotopy deformation of the potential is shown to extend the range of reliably computable overtones by one or two steps.  A key advertised advantage is that the extremal RNBH case requires only the change of the tortoise coordinate, whereas the rest of the procedure is identical to the nonextremal case, in contrast to the significant reformulation needed in the continued fraction method.","tokens_in":12687,"tokens_out":2995,"duration_ms":29552,"significance":"If the results hold, this work provides a practical and conceptually clean alternative to the continued fraction method for spherically symmetric black holes with complicated singularity structures.  The demonstration that extremal and nonextremal RNBHs can be treated uniformly without modifying the core algorithm is a genuine simplification that could facilitate future applications to other charged or exotic black holes.  The method's accuracy for low-lying modes (n ≤ 2) is competitive with established techniques, and the explicit error reporting (Tables III and IV) is commendable.  The homotopy deformation, though ad hoc, offers a workable route to access slightly higher overtones that are difficult to reach with the basic version.  The paper does not claim universal applicability beyond the demonstrated cases, and its limitations are honestly discussed.","major_comments":[{"comment":"The abstract and text state that the method yields QNM frequencies of 'high accuracy for low-lying modes with overtones n ≤ 3.'  However, for gravitational perturbations with n = 3, the relative errors in Table III range from 1×10⁻³ to 1×10⁻² (i.e., up to ~1%).  While these values are still acceptable for many purposes, they are orders of magnitude larger than the errors for n ≤ 2 (≲2×10⁻⁵).  The authors should qualify the n = 3 accuracy as 'good' or 'reliable' rather than 'high accuracy' to avoid overstatement, or alternatively provide a clear threshold that defines 'high accuracy' and apply it consistently.","section":"Abstract and §III.A"},{"comment":"The homotopy parameter λ is chosen by trial (λ = 4 works for the cases shown), and the manuscript notes that 'λ is not necessarily the larger the better' and that singularities approach the target point α = −i as λ increases.  No systematic prescription for selecting λ is given.  While this does not invalidate the presented results, the method's reliance on manual tuning for λ weakens its claim to be a 'method' for higher overtones.  The authors should at least discuss criteria or heuristics (e.g., a scan over λ with a convergence diagnostic) that a future user could apply, or explicitly state that the choice remains an open issue.","section":"§III.B, Eq. (10) and surrounding text"}],"minor_comments":[{"comment":"The transformation from Eq. (8) to Eq. (9) substitutes −ix → αx.  It would be helpful to state explicitly that α is a real parameter that will later be analytically continued to α = −i (the target value).  Currently this is implied but could be clarified for a reader unfamiliar with [47].","section":"§II, Eq. (9)"},{"comment":"The bound state method results for n = 0,1,2 are given to five decimal places, matching the reference values exactly in many cases.  However, the relative errors are reported in separate tables.  For quick reference, consider merging the error information into the main tables (e.g., parenthetical errors) or at least stating the range of errors for each n in the text.","section":"§III.A, Tables I–IV"},{"comment":"Typographical issues: 'high accurate calculation' should be 'highly accurate calculation'; 'complicate singularity structures' should be 'complicated singularity structures'.  Also 'particular useful' should be 'particularly useful'.","section":"§IV, Conclusion"},{"comment":"The caption of Figure 1 refers to a gravitational QNM with l = 2, Q = 0.2, n = 3.  The text states that the figure demonstrates convergence with increasing Padé order.  The figure is not reproduced in the arXiv text; presumably it will be rendered in the final version.  Ensure that the axes are clearly labeled and that the reference values (dashed lines) are identified in the legend.","section":"§III.B, Fig. 1"},{"comment":"Reference [47] (the authors' previous work on the Schwarzschild case) is listed as arXiv:2608.19597.  As of this writing, this preprint is not yet published.  If it has been accepted or published by the time of submission, please update the reference accordingly.","section":"References"}],"recommendation":"minor_revision","confidential_remarks":"The paper is sound and the validation is thorough.  The main claim for low overtones (n ≤ 2) is robustly supported, and the limitation for higher overtones is openly discussed.  The two major comments above are moderate: one concerns the phrasing of accuracy for n = 3 (easily fixed), and the other points to the ad hoc nature of the homotopy parameter, which is a known limitation that does not undermine the central results.  I recommend minor revision.  The scope is appropriate for a journal like Physical Review D."},"author_rebuttal":null,"desk_editor":{"model":"deepseek/deepseek-v4-flash","letter":"This is a competent, incremental extension of the authors' own bound state method from Schwarzschild to Reissner-Nordström black holes. The main new result is that the method handles the extremal case without special machinery—just a change in the tortoise coordinate—whereas Leaver's continued fraction method requires a fundamentally different treatment. That is a real practical advantage, and the authors demonstrate it clearly with numerical results that match Leaver's to high precision for low-lying modes (n ≤ 2).\n\nWhat I like: the paper is honest about its limitations. The method fails for higher overtones (n ≥ 3–4 without homotopy, n ≥ 5–6 with it), and the authors say so. They include convergence plots and show that the Padé continuation has a finite radius of validity. The homotopy deformation is clearly an ad hoc fix (λ chosen by trial), but they flag that too and show it works for a few extra overtones. The numerical validation against Leaver's method is thorough, with relative errors reported explicitly.\n\nThe soft spots are proportionate. (1) The method is not new—it was introduced in the authors' earlier paper (ref [47]). This is a straightforward application to a different black hole, so the novelty is modest. (2) The homotopy parameter λ is not determined systematically; it is tuned by hand. The paper acknowledges this but does not offer a prescription for choosing λ a priori. (3) For non-extremal RN black holes, the bound state method offers little advantage over the well-established continued fraction method—the real benefit is for the extremal case. (4) The accuracy degrades quickly with overtone number, and the method cannot compete with Leaver's for high overtones. That said, the authors never claim otherwise.\n\nOverall, this is a solid, honest paper that delivers what it promises: a simpler way to compute low-lying QNM frequencies for extremal RN black holes, with cross-validated numbers. It is not groundbreaking, but it is useful within the subfield of black hole perturbation theory. I would bring it to a reading group focused on computational methods, but I would not cite it unless I needed those specific RN numbers. I recommend acceptance; it deserves a serious referee and will likely pass with minor revisions.","headline":"Solid incremental application of the bound state method to RN black holes, with honest validation and a clear advantage for the extremal case, but no fundamental novelty.","tokens_in":13168,"tokens_out":1184,"would_cite":false,"duration_ms":13615,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek/deepseek-v4-flash","headline":"A single bound-state computation, analytically continued to the imaginary axis, reproduces Reissner-Nordström black hole quasinormal modes to roughly one part in 100,000 and treats the extremal charge limit without changing the algorithm.","keywords":["quasinormal modes","Reissner-Nordström black holes","bound state method","analytic continuation","Padé approximants","homotopy deformation","extremal black holes","tortoise coordinate"],"falsifier":"For an $\\ell=2$ mode with charge $Q=0.4$, compute the $n=4$ gravitational and $n=5$ electromagnetic frequencies by an independent method, such as extracting them from a time-domain evolution of the perturbation equations. If these independent values disagree with the homotopy-deformed bound-state results beyond the tabulated relative errors, the rational-function continuation from the 31-point grid is the failing step; if they agree, the failure without deformation is a limitation of the fitting scheme rather than of the bound-state mapping. Repeating the test at $Q=M$ for $n=3$ would extend the check to the extremal endpoint, where only $n\\le2$ is compared.","tokens_in":14155,"feed_emoji":"🕳️","tokens_out":13501,"duration_ms":109941,"temperature":0.7,"pith_summary":"The paper claims that a recently introduced single-parameter bound-state method, previously tested on Schwarzschild black holes, also computes charged Reissner-Nordström quasinormal-mode frequencies accurately. For low-lying overtones, the results agree with the standard continued-fraction calculation to about $10^{-5}$ or better, and a homotopy deformation of the potential well extends the reliable range by one or two further overtones. The central payoff is uniformity: the extremal case $Q=M$ is handled by the same pipeline as the nonextremal case, with only the tortoise coordinate replaced by its extremal form, whereas the continued-fraction method needs a separate and more complicated construction in that limit. The paper itself reports that the undeformed method fails for higher overtones and that the analytic continuation of the bound-state energies is the limiting ingredient. Taking the argument in good faith, the contribution is a demonstration that quasinormal spectra for this whole black hole family can be extracted from ordinary bound-state data.","feed_headline":"Bound states compute black hole ringdown to 1e-5","feed_subtitle":"Deforming the potential well unlocks higher overtones; maximally charged black holes need no separate algorithm.","key_machinery":"The load-bearing object is the auxiliary scaling parameter $\\alpha$: after the quasinormal-mode-to-bound-state map $x\\to -ix$, the substitution $-ix\\to \\alpha x$ makes the inverted potential depend on the real parameter $\\alpha$, so that for real $\\alpha$ the differential equation is a genuine bound-state problem with discrete energies $E_n(\\alpha)$. A shooting method gives $E_n(\\alpha)$ on 31 uniformly spaced points with $\\alpha\\in[0.2,0.8]$, and diagonal or subdiagonal Padé approximants continue $E_n(\\alpha)$ to $\\alpha=-i$, where the original outgoing-wave problem is recovered. The homotopy deformation multiplies the potential well by $H(\\alpha,\\lambda)=(\\lambda-1)\\alpha^2+\\lambda$; since $H(-i,\\lambda)=1$ for every $\\lambda$, the quasinormal-mode problem is unchanged at the target point even though the real-$\\alpha$ well is made deeper, which moves the singularities of $E_n(\\alpha)$ away from the continuation path. The extremal case enters only through the form of the tortoise coordinate, the single place where the horizon structure changes the calculation.","core_discovery":"The central claim is that the quasinormal frequencies of odd-parity electromagnetic and gravitational perturbations of a Reissner-Nordström black hole follow from a bound-state spectrum. In tortoise coordinate $x$, the map $x\\to -ix$ turns the potential barrier into a well; introducing a real scaling parameter $\\alpha$ through $-ix\\to \\alpha x$ produces an ordinary Schrödinger eigenvalue problem with energies $E_n(\\alpha)$. These energies are computed numerically on a grid of real $\\alpha$, represented by Padé approximants, and continued analytically to $\\alpha=-i$, where the quasinormal frequency is $\\omega_n=\\sqrt{-E_n(-i)}$. For the tabulated nonextremal cases ($\\ell=2$, $Q$ from $0.2$ to $0.99$), the agreement with continued-fraction reference values is better than about $10^{-5}$ for the lowest overtones, with gravitational modes computed through $n=3$ and electromagnetic modes through $n=4$; the homotopy deformation $H(\\alpha,\\lambda)=(\\lambda-1)\\alpha^2+\\lambda$ with $\\lambda=4$ gives reliable higher overtones, reaching $n=4$ for gravitational modes, $n=5$ for electromagnetic modes, and $n=6$ for electromagnetic modes with $Q>0.8$. In the extremal case the same procedure, using $x=r-M^2/(r-M)+2M\\ln(r-M)$ as the tortoise coordinate, reproduces the known frequencies to five decimal places for $n\\le2$ and exhibits the isospectrality between electromagnetic modes of index $\\ell$ and gravitational modes of index $\\ell+1$, a degeneracy already visible in the bound-state energies.","pith_inferences":["One testable extension is to treat the homotopy parameter as a convergence diagnostic: plotting where the Padé approximants' singularities move as $\\lambda$ varies would give a practical warning when a claimed high-overtone frequency is about to become unreliable, since the paper identifies those singularities approaching $\\alpha=-i$ as the cause of failure.","The same real-$\\alpha$ continuation strategy should carry over to other spherically symmetric spacetimes with a single effective potential barrier, but the grid span, sampling density, and optimal homotopy value would need to be re-tuned; the paper does not claim a universal prescription.","For gravitational-wave template construction, a practical consequence would be a single code covering the whole charge range $0\\le Q\\le M$, including the maximally charged endpoint, because the extremal calculation is not algorithmically special.","A direct check of whether the limiting singularities of $E_n(\\alpha)$ coincide with the branch points of the inverse relation $\\omega^2=-E$ would separate a fundamental obstruction to high overtones from an artifact of the chosen Padé grid."],"forward_implications":["Nonextremal Reissner-Nordström low-lying quasinormal modes can be obtained from bound-state energies with relative errors near or below $10^{-5}$, giving an independent numerical cross-check of continued-fraction results.","Extremal black hole quasinormal modes can be computed with the same code path as nonextremal ones: only the tortoise coordinate must be replaced by its extremal form, with no new recurrence relations or singularity analysis required.","The homotopy parameter $\\lambda=4$ enlarges the computable range by roughly one or two overtones, extending gravitational results to $n=4$ and electromagnetic results to $n=5$, and to $n=6$ when $Q>0.8$.","The extremal isospectrality between electromagnetic modes of index $\\ell$ and gravitational modes of index $\\ell+1$ is reproduced, and it is shown to be present already at the level of the underlying bound-state energies."],"supporting_citations":[{"why":"Introduces the single-parameter bound-state method and the homotopy deformation that this paper applies to Reissner-Nordström black holes.","marker":"[47]"},{"why":"Establishes the continued-fraction reference calculation for nonextremal Reissner-Nordström quasinormal modes against which the new results are measured.","marker":"[15]"},{"why":"Supplies the extremal Reissner-Nordström continued-fraction reference values and the reported isospectrality between electromagnetic and gravitational modes.","marker":"[16]"},{"why":"Introduces the original correspondence between black hole quasinormal modes and bound states of the inverted potential, which the new method extends.","marker":"[30]"},{"why":"Shows that numerically computed bound-state energies can be continued analytically in model parameters, removing the need for potentials with closed-form bound-state spectra.","marker":"[45]"},{"why":"Provides the odd-parity perturbation equations and the two potentials for electromagnetic and gravitational perturbations that the calculation starts from.","marker":"[6]"}],"fun_headline_variants":["Bound state trick yields black hole ringdown to 1e-5","Deforming potential well reveals higher QNM overtones","Unified method for extremal and nonextremal black hole QNMs","High‑accuracy QNMs from a bound state spectrum"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The load-bearing premise is that a rational-function fit to the bound-state energies at 31 real values of the scaling parameter ($0.2\\le\\alpha\\le0.8$) continues accurately to $\\alpha=-i$; the paper itself reports this fails for gravitational overtones $n\\ge4$ and electromagnetic $n\\ge5$ before deformation, and it gives no general guarantee that the grid, fit order, or the fixed homotopy value $\\lambda=4$ is optimal.","fun_headline_variants_meta":{"raw":{"variants":["Bound state trick yields black hole ringdown to 1e-5","Deforming potential well reveals higher QNM overtones","Unified method for extremal and nonextremal black hole QNMs","High‑accuracy QNMs from a bound state spectrum"]},"model":"deepseek/deepseek-v4-flash","effort":"low","cost_usd":0.000798,"raw_usage":{"cost":0.000798448,"is_byok":false,"cost_details":{"upstream_inference_cost":0.000798448,"upstream_inference_prompt_cost":0.000120568,"upstream_inference_completions_cost":0.00067788},"total_tokens":3487,"prompt_tokens":1066,"completion_tokens":2421,"prompt_tokens_details":{"audio_tokens":0,"video_tokens":0,"cached_tokens":256,"cache_write_tokens":0},"completion_tokens_details":{"audio_tokens":0,"image_tokens":0,"reasoning_tokens":2346}},"tokens_in":1066,"tokens_out":2421,"cost_source":"openrouter_usage","duration_ms":23001,"image_pages":[],"temperature":1.0,"reasoning_tokens":2346,"transport_provider":"openrouter","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-09-15T12:01:34.651906+00:00","model_set":{"reader":"deepseek/deepseek-v4-flash"},"falsifier":"For an $\\ell=2$ mode with charge $Q=0.4$, compute the $n=4$ gravitational and $n=5$ electromagnetic frequencies by an independent method, such as extracting them from a time-domain evolution of the perturbation equations. If these independent values disagree with the homotopy-deformed bound-state results beyond the tabulated relative errors, the rational-function continuation from the 31-point grid is the failing step; if they agree, the failure without deformation is a limitation of the fitting scheme rather than of the bound-state mapping. Repeating the test at $Q=M$ for $n=3$ would extend the check to the extremal endpoint, where only $n\\le2$ is compared.","supporting_citations":[{"cited_title":"A New Method for Quasinormal Modes From Bound States and Homotopy deformations","cited_arxiv_id":"2608.19597","evidence_quote":"Introduces the single-parameter bound-state method and the homotopy deformation that this paper applies to Reissner-Nordström black holes."},{"cited_title":"Quasinormal modes of Reissner-Nordstrom black holes, Phys. Rev. D41, 2986","cited_arxiv_id":null,"evidence_quote":"Establishes the continued-fraction reference calculation for nonextremal Reissner-Nordström quasinormal modes against which the new results are measured."},{"cited_title":"Quasinormal modes of maximally charged black holes","cited_arxiv_id":"gr-qc/9603021","evidence_quote":"Supplies the extremal Reissner-Nordström continued-fraction reference values and the reported isospectrality between electromagnetic and gravitational modes."},{"cited_title":"Oscillations of a Black Hole, Phys. Rev. Lett.52, 1361","cited_arxiv_id":null,"evidence_quote":"Introduces the original correspondence between black hole quasinormal modes and bound states of the inverted potential, which the new method extends."},{"cited_title":"Quasi-Normal Modes from Bound States: The Numerical Approach","cited_arxiv_id":"2210.02069","evidence_quote":"Shows that numerically computed bound-state energies can be continued analytically in model parameters, removing the need for potentials with closed-form bound-state spectra."},{"cited_title":"Effective potential for even parity Regge- Wheeler gravitational perturbation equations, Phys. Rev. Lett.24, 737","cited_arxiv_id":null,"evidence_quote":"Provides the odd-parity perturbation equations and the two potentials for electromagnetic and gravitational perturbations that the calculation starts from."}],"review_version":1}