{"id":"a698604e-8a46-4a03-a31c-e85068f39646","arxiv_id":"2506.04326","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The authors formulate Kerr quasinormal modes as eigenvalues of a two-dimensional spectral problem, enabling simultaneous extraction of many modes and showing that near-horizon gradients in the extremal limit are coordinate artifacts.","lead":"This paper presents a new numerical method for computing quasinormal mode frequencies of Kerr black holes by treating the Teukolsky equation as a two-dimensional eigenvalue problem, bypassing the traditional separation into radial and angular equations. The method extracts many modes at once without initial guesses, which could speed up black hole spectroscopy and studies of mode stability.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim lacks independent validation: the convergence tests compare the solver against itself, so the filtered eigenvalues could converge to the wrong QNM spectrum if the boundary-condition encoding or the filter is flawed.","rationale":"The reader's weakest assumption is precisely that the filtered eigenvalues of the discretized hyperboloidal m-mode operator coincide with the physical QNMs. My analysis agrees: the paper's evidence for this is internal convergence to a self-defined reference, not independent validation. The lack of a quantitative benchmark against established Kerr QNM values is therefore the most load-bearing gap. The proposed concrete test directly addresses this gap by comparing against an independent solver across spins and overtones, including the near-extremal regime where the filter is most fragile. Since this concern does not contradict the reader's conditional verdict—it strengthens the reason for conditionality—the verdict should remain unchanged. I see no additional, more severe internal inconsistency that would warrant rejection; the method is well-posed, the operator construction is explicit, and the internal convergence plots are credible. The concern is about validation, not about an identified error in the derivation.","tokens_in":29080,"tokens_out":5325,"duration_ms":50452,"concrete_test":"Implement an independent reference solver, e.g., Leaver's continued fraction method or a published high-accuracy Kerr QNM code, and compute the first several QNMs for a representative set of parameters: s = -2, m = 2 and m = 1, (ℓ, n) = (2,0), (2,1), (3,0), (4,0), for a/M = 0, 0.5, 0.8, 0.9, 0.99. Then run the m-mode solver with the same parameters and the paper's filtering criterion (n_low1=50, n_high1=55, n_low2=10, n_high2=11, TOL=1e-3), and compute relative errors |ω_solver - ω_independent|/|ω_independent| for every retained eigenvalue. If any retained eigenvalue differs by more than, say, 10^-6 (or by more than the claimed convergence accuracy in Fig. 4), the boundary-condition encoding or the filter is admitting or discarding the wrong modes, and the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that the eigenvalues of the discretized hyperboloidal m-mode operator (Eq. 20), after filtering, are the physical Kerr QNMs. The load-bearing assumption is that the coordinate transformation (7), the regularity factor (15), and the convergence filter (Sec. 4, TOL ~ 1e-3) together correctly encode the ingoing/outgoing boundary conditions and select exactly the physical modes. This assumption is not independently tested: the convergence analysis in Sec. 4.2 defines the reference values ω_ref using the same solver at N_ref = 55, so the reported error ϵ measures self-convergence, not accuracy. A method can converge exponentially to the wrong eigenvalues if the boundary conditions are mis-specified or if the filter admits spurious eigenvalues that happen to stabilize with resolution. The only external anchor is the qualitative Schwarzschild benchmark in Fig. 2, which is not quantified. The paper even notes that retrograde modes converge slowly near the branch cut as a/M approaches 1, so the filter's reliability is least secure exactly in the near-extremal regime where the method is promoted. Without a direct quantitative comparison to established QNM values (e.g., Leaver's continued fraction results), the central claim that the extracted spectrum is the true Kerr QNM spectrum is not fully established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a numerical method for computing quasinormal modes (QNMs) of the Kerr spacetime by formulating the frequency-domain Teukolsky equation as a two-dimensional eigenvalue problem. The authors combine a hyperboloidal compactification with an azimuthal m-mode decomposition, obtaining an operator whose eigenvalues directly yield QNM frequencies without a separation constant or root-finding. They present the discretized operator, a filtering procedure to remove spurious eigenvalues, convergence tests for two hyperboloidal gauges, a proposed re-labelling of overtones by a single index q, and studies of the angular and radial structure of the eigenfunctions. The central advertised results are that the method extracts the Kerr QNM spectrum directly as eigenvalues, that both gauges perform comparably, and that strong near-horizon gradients seen in one gauge are coordinate artefacts.","tokens_in":29338,"tokens_out":4557,"duration_ms":46933,"significance":"If the central claim is correct, the paper provides a practical and conceptually clean route to Kerr QNM spectra: for each azimuthal number m, all overtones and angular harmonics are obtained simultaneously as eigenvalues of a single discretized operator, with no need for initial guesses. This would be a useful tool for QNM stability studies, pseudospectrum computations, and mode-excitation problems, and the paper's demonstration of equivalent performance between two hyperboloidal gauges is a valuable practical message. The m-mode projection onto both spheroidal and spherical harmonic bases is also potentially useful for gravitational-wave data analysis. However, the paper's validation is entirely self-referential: the convergence study compares the solver to its own high-resolution results, and the only independent anchor is a qualitative Schwarzschild plot. Because the correctness of the extracted spectrum is the load-bearing claim, the lack of a direct quantitative comparison with established Kerr QNM values leaves the central result not fully established.","major_comments":[{"comment":"The convergence error ϵ in Eq. (59) is computed relative to ω_Ref obtained from the same solver at N_ref = 55. This is a self-consistency check and cannot detect a systematic error in the boundary-condition encoding or in the filtering procedure. The only independent anchor is the qualitative Schwarzschild benchmark in Fig. 2, which is not quantified. Since the paper's central claim is that the filtered eigenvalues are the physical Kerr QNM spectrum, please add a quantitative comparison against independent published values—for example Leaver's continued-fraction results or the high-accuracy data of Cook and Zalutskiy—for representative modes, for spins including a/M = 0.8, 0.9, and 0.99, and for both prograde and retrograde branches.","section":"Sec. 4.2, Eq. (59)"},{"comment":"The filtering procedure retains eigenvalues satisfying |1 − ω_low/ω_high| < TOL, but the paper reports no sensitivity analysis with respect to TOL or the truncation pair (n_low, n_high), nor a demonstration that spurious eigenvalues near the branch cut are always discarded. The concern is most acute precisely where the paper itself notes slow convergence, namely retrograde modes as a/M approaches 1 and spurious eigenvalues cluster near the branch cut. Without a robustness study, the filter could in principle admit or discard the wrong eigenvalues. Please quantify the dependence of the retained spectrum on TOL and resolution, and cross-check the filtered values against independent Kerr QNM data.","section":"Sec. 4, filtering criterion with TOL ≈ 1e-3"},{"comment":"The conclusion that the near-horizon gradients in the radial-fixing gauge are coordinate artefacts is inferred from comparing eigenfunctions in the radial-fixing and Cauchy-horizon-fixing gauges. These gauges use different radial coordinate choices (ρ_o = 0 versus ρ_o = κ²), so the same physical field is represented by different functions of σ; the absence of steep gradients in one gauge does not by itself prove that the gradients are unphysical. Please compare an invariant quantity, such as the original Teukolsky master function or a suitably normalized covariant quantity, or explicitly justify why the comparison of the two gauges is conclusive.","section":"Sec. 4.3.1 and Fig. 6"}],"minor_comments":[{"comment":"The text states \"For the radial fixing gauge (σc = κ2)\", but Eq. (8) with ρ_o = 0 gives σ_c = κ^{-2}, and the factor (1 − κ²σ) in Eq. (22) is consistent with σ_c = κ^{-2}; please correct this typo.","section":"Sec. 2.1.1, before Eq. (21)"},{"comment":"There are several typographical errors that should be corrected, including \"respecvelty\", \"Teulkolsky\", \"straightfoward\", \"the the\", and \"discretized operador\".","section":"Throughout"},{"comment":"The horizontal axis is described as the \"total size\" of the discretized operator; please clarify that it is the total number of grid points n_total = n1 × n2, since the truncation is parameterized by N1 = 5N2 = N.","section":"Fig. 4 caption"},{"comment":"The proposed ordering by decay rate is natural, but the statement that even q corresponds to prograde and odd q to retrograde modes requires the caveat, already partly given, that for m = 0 the q = 0 and q = 1 modes are degenerate in frequency; please make this caveat explicit in the definition of the ordering.","section":"Sec. 4.1, Eq. (56)"}],"recommendation":"major_revision","confidential_remarks":"This is a methods paper whose main novelty is the direct 2D eigenvalue formulation of the Kerr QNM problem. The derivation and numerical scheme appear sound, and the paper is well within the scope of the journal. The primary obstacle to acceptance is the validation gap: all convergence tests are self-referential and no independent quantitative benchmark is provided. I would require the authors to add a direct comparison with established Kerr QNM data, including the near-extremal retrograde regime, and to document the sensitivity of the filtering to its parameters, before this can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing is the 2D eigenvalue formulation for Kerr QNMs: no separation constant, no root-finding, all overtones for fixed m come out of one operator spectrum. That is a real methodological step, and the paper does it cleanly. The explicit operator coefficients for both hyperboloidal gauges are given, convergence is exponential, and the radial-fixing vs Cauchy-horizon-fixing comparison is useful. The eigenfunctions are regular across the domain, the angular projection onto spheroidal and spherical harmonics works as advertised, and the proposed q notation is a sensible way to label overtones by decay rate rather than by the regular/mirror convention.\n\nThe main weakness is validation. The convergence study compares the solver against its own high-resolution result, not against independent published values. A spectral method can in principle converge to the wrong eigenvalues if the boundary-condition encoding or the filter is flawed. That risk is mitigated here: the hyperboloidal framework is standard, the eigenvalues are not fitted parameters, and the Schwarzschild spectrum shows up in Fig. 2. But the Schwarzschild check is not quantified. A table comparing against Leaver's continued-fraction values for a few spins and several (ell,n) modes would settle it. The claim that the near-horizon gradients are coordinate artifacts is plausible but rests on one mode and two gauges; that part is thinner than the rest of the paper. The paper also notes that retrograde modes converge slowly near extremality, so the filter tolerance is least certain exactly where the method is being promoted.\n\nThis is a methods paper for QNM practitioners and anyone doing black-hole spectroscopy. The formulation is clear, the math is explicit, and the gaps are fixable. I would send it to a serious referee with a request for an independent benchmark table and a more careful discussion of filter reliability in the near-extremal regime. It should not be desk-rejected.","headline":"The 2D eigenvalue formulation is a real step forward, but validation is self-referential and needs an independent benchmark before the method is fully trusted.","tokens_in":29846,"tokens_out":2575,"would_cite":true,"duration_ms":34458,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","83C35","65N25"],"pacs":["04.70.-s","04.30.-w"],"model":"deepseek-v4-flash","headline":"The Kerr quasinormal-mode spectrum for each azimuthal number $m$ can be computed as a direct eigenvalue problem on a two-dimensional hyperboloidal slice, eliminating root-finding and initial guesses.","keywords":["quasinormal modes","Kerr black holes","Teukolsky equation","hyperboloidal compactification","eigenvalue problem","m-mode decomposition","black hole spectroscopy","spheroidal harmonics"],"falsifier":"Run the solver at several resolutions for a fixed $m$ and compare every filtered eigenvalue against independent high-accuracy continued-fraction values across many spins and overtones: if any retained eigenvalue fails to match a known mode, or any known mode is missing at tolerance $10^{-3}$, the central claim is refuted. A complementary check is to substitute a computed eigenpair back into the original Teukolsky equation and verify that the residual and the boundary behaviour decay at the expected rates.","tokens_in":28876,"feed_emoji":"🕳️","tokens_out":8200,"duration_ms":71126,"temperature":0.7,"pith_summary":"The paper claims that the quasinormal-mode problem for a Kerr black hole, normally solved by treating the radial and angular equations separately with a coupling separation constant and an initial guess for every mode, can be recast as one genuine eigenvalue problem for each azimuthal number $m$. The key move is to keep the Teukolsky equation as a two-dimensional elliptic PDE on a hyperboloidal slice, so the outgoing and ingoing boundary conditions are built into the coordinates rather than imposed by hand. If that works, an entire tower of overtones and angular harmonics for a fixed $m$ is obtained from a single matrix diagonalization, with no seed values. The paper also uses the complete spectrum to propose a cleaner single-index overtone label $q$, and to show that steep near-horizon gradients found by earlier radial hyperboloidal solvers are coordinate artifacts rather than physical features of extremal Kerr. This matters because quick, assumption-free spectrum construction is the natural input to black hole spectroscopy and to studies of mode stability.","feed_headline":"Kerr quasinormal modes without initial guesses","feed_subtitle":"A 2D hyperboloidal reformulation of the Teukolsky equation returns the full per-azimuth spectrum at once.","key_machinery":"The load-bearing object is the hyperboloidal $m$-mode Teukolsky operator: after the coordinate transformation (7), the regularity factor (15), and the azimuthal Fourier decomposition (16), the master equation becomes the 2D elliptic PDE (17), written as $s\\bar D_{m;\\bar\\omega}\\bar\\Phi_{m;\\bar\\omega}=0$ with spatial operators $L_1$ and $L_2$. Introducing the auxiliary function $\\bar\\Upsilon_{m;\\bar\\omega}=s\\bar\\Phi_{m;\\bar\\omega}$ converts it into the linear eigenvalue problem $Lu=s u$ of Eq. (20). Discretizing $u$ on a Chebyshev tensor grid renders the spectrum as eigenvalues of a finite matrix, with boundary conditions at null infinity and the horizon enforced geometrically by the hyperboloidal coordinates and by the regularity prefactors in Eq. (15).","core_discovery":"The central discovery is that separating variables is unnecessary for computing Kerr quasinormal modes. Combining the hyperboloidal compactification with a decomposition into azimuthal modes $m$ turns the frequency-domain Teukolsky equation into a two-dimensional linear operator whose eigenvalues are the QNM frequencies, with the auxiliary function $s\\bar\\Upsilon_{m;\\bar\\omega}=s\\,s\\bar\\Phi_{m;\\bar\\omega}$ making the problem first order in the spectral parameter. The paper demonstrates the resulting spectra for gravitational perturbations $s=-2$ across spins from Schwarzschild to $a/M=0.99$, reproduces known values, and shows that both hyperboloidal gauges converge exponentially to the same frequencies. It then exploits the directly available eigenfunctions to show that the near-extremal gradients are gauge-dependent, and to project the angular profile onto both spheroidal and spherical harmonic bases. In the authors' reading, the QNM problem is not two coupled ODE eigenvalue problems but a single 2D spectral problem for each $m$.","pith_inferences":["Beyond the paper, the same no-seed eigenvalue setup should make pseudospectrum and QNM-instability analyses in Kerr substantially cheaper, because the full eigensystem is available in one solve rather than mode by mode.","The $q$-index ordering by decay rate could serve as a standard for ringdown comparisons, removing the ambiguity in which one overtone label $n$ covers four distinct modes.","Because near-horizon gradient formation is gauge-dependent, any physical claim tied to eigenfunction steepness near extremality, such as an instability of the horizon, should be checked in a gauge-invariant quantity rather than in a fixed slicing.","A direct extension is to Kerr-Newman or other non-separable backgrounds, where the same 2D diagonalization could replace the Newton-Raphson root searches currently used."],"forward_implications":["For a fixed $m$, one matrix eigenvalue problem returns many QNMs at once, so constructing the Kerr spectrum no longer needs a root finder with a good initial guess.","The single overtone index $q$, ordered by decay rate, replaces the four-way labels and exposes the prograde/retrograde ordering directly.","Radial-fixing and Cauchy-horizon-fixing gauges give the same frequencies with the same exponential convergence, so the gauge choice does not affect spectral accuracy.","Steep near-horizon gradients in extremal-Kerr eigenfunctions are slicing artifacts, not physical features, resolving a question raised by earlier radial computations.","QNM eigenfunctions can be projected onto either spin-weighted spheroidal or spherical harmonics, easing direct comparison with gravitational-wave templates."],"supporting_citations":[{"why":"Supplies the minimal-gauge hyperboloidal coordinates and the m-mode Teukolsky operator on which the eigenvalue problem is built.","marker":"[41]"},{"why":"Provides the continued-fraction QNM baseline that the Cauchy-horizon-fixing gauge reproduces, used for benchmarking.","marker":"[23]"},{"why":"Earlier hyperboloidal radial-eigenvalue solver whose reported near-horizon gradients the paper shows to be gauge artifacts.","marker":"[42]"},{"why":"Establishes the genuine-eigenvalue treatment and spurious-eigenvalue filtering for QNM problems that the 2D method adapts.","marker":"[30]"},{"why":"Prior m-mode hyperboloidal pseudospectrum calculation showing the strategy is viable for Kerr; motivation for this systematic study.","marker":"[48]"},{"why":"Documents the two gauge choices and their Carter-Penrose geometry, including notation and Cauchy-horizon behavior.","marker":"[22]"},{"why":"PDE treatment of Kerr-Newman QNMs that still needs Newton-Raphson root finding; the contrast motivating direct eigenvalue formulation.","marker":"[46]"},{"why":"Provides high-accuracy angular spectral methods and spheroidal/spherical harmonic data used for comparison and projection tests.","marker":"[25]"}],"fun_headline_variants":["Kerr quasinormal modes as 2D eigenvalue problem","No separation needed: Kerr QNMs from 2D operator","Hyperboloidal 2D reformulation yields entire Kerr spectrum","Kerr spectra without angular separation","2D eigensolver extracts Kerr quasinormal modes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the eigenvalues of the discretized hyperboloidal operator that survive the tolerance filter are the actual quasinormal modes of the Teukolsky equation with the correct boundary conditions, rather than numerical artifacts of the discretization.","fun_headline_variants_meta":{"raw":{"variants":["Kerr quasinormal modes as 2D eigenvalue problem","No separation needed: Kerr QNMs from 2D operator","Hyperboloidal 2D reformulation yields entire Kerr spectrum","Kerr spectra without angular separation","2D eigensolver extracts Kerr quasinormal modes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000251,"raw_usage":{"total_tokens":1599,"prompt_tokens":1030,"completion_tokens":569,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":646,"completion_tokens_details":{"reasoning_tokens":489}},"tokens_in":646,"tokens_out":569,"duration_ms":5347,"temperature":1.0,"reasoning_tokens":489,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:43:42.526819+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the solver at several resolutions for a fixed $m$ and compare every filtered eigenvalue against independent high-accuracy continued-fraction values across many spins and overtones: if any retained eigenvalue fails to match a known mode, or any known mode is missing at tolerance $10^{-3}$, the central claim is refuted. A complementary check is to substitute a computed eigenpair back into the original Teukolsky equation and verify that the residual and the boundary behaviour decay at the expected rates.","supporting_citations":[{"cited_title":"Hyperboloidal framework for the Kerr spacetime","cited_arxiv_id":null,"evidence_quote":"Supplies the minimal-gauge hyperboloidal coordinates and the m-mode Teukolsky operator on which the eigenvalue problem is built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the continued-fraction QNM baseline that the Cauchy-horizon-fixing gauge reproduces, used for benchmarking."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier hyperboloidal radial-eigenvalue solver whose reported near-horizon gradients the paper shows to be gauge artifacts."},{"cited_title":"Pseudospectrum and black hole quasinormal mode instability","cited_arxiv_id":null,"evidence_quote":"Establishes the genuine-eigenvalue treatment and spurious-eigenvalue filtering for QNM problems that the 2D method adapts."},{"cited_title":"The pseudospectrum for the Kerr black hole: spin s = 0case","cited_arxiv_id":null,"evidence_quote":"Prior m-mode hyperboloidal pseudospectrum calculation showing the strategy is viable for Kerr; motivation for this systematic study."},{"cited_title":"The confluent Heun functions in black hole perturbation theory: a spacetime interpretation","cited_arxiv_id":null,"evidence_quote":"Documents the two gauge choices and their Carter-Penrose geometry, including notation and Cauchy-horizon behavior."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"PDE treatment of Kerr-Newman QNMs that still needs Newton-Raphson root finding; the contrast motivating direct eigenvalue formulation."},{"cited_title":"Cook and Maxim Zalutskiy","cited_arxiv_id":null,"evidence_quote":"Provides high-accuracy angular spectral methods and spheroidal/spherical harmonic data used for comparison and projection tests."}],"review_version":1}