{"id":"5a15b57a-4399-44fe-a534-170c65c0d78f","arxiv_id":"2507.01919","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Kerr quadratic quasi-normal mode amplitudes and phases are computed analytically in the eikonal limit via the Penrose limit, giving an explicit spin-dependent nonlinearity ratio.","lead":"This paper derives an analytic formula for the size and phase of the nonlinear quadratic ringdown mode of a spinning Kerr black hole in the large-angular-momentum limit. The formula could help future gravitational-wave detectors interpret nonlinear mode mixing in black hole merger signals.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (59) ignores an O(1) frequency detuning between the sum of the linear modes and the nonlinear QNM; the retained spin term in Eq. (63) is below the controlled order, so the spin dependence is not yet quantitatively established.","rationale":"The reader's weakest_assumption (light-ring localization plus the inconsistent retention of the subleading α term) is real and well identified. My pass sharpens it: the largest unquantified error is not only the localization itself but the O(1) detuning between the sum of linear QNM frequencies and the nonlinear mode frequency, which Eq. (59) sets equal at leading order. This detuning enters the second-order amplitude through the Green's-function denominator and is therefore as important as the terms that set the prefactor of Eq. (63), while the retained α^4/ω^4 term is far smaller in the eikonal hierarchy. The paper is transparent about several limitations and the framework is plausible, so this does not call for rejection; it strengthens the case for a conditional verdict pending a calculation that keeps the finite source width and the actual QNM frequency. I therefore leave the reader's CONDITIONAL verdict unchanged.","tokens_in":18424,"tokens_out":34848,"duration_ms":399641,"concrete_test":"For Schwarzschild (a=0), where α^4/ω^4 is negligible and C_ℓ reduces to a Gaunt integral, build the quadratic Teukolsky source from the eikonal linear solutions of App. C with ℓ1=2, ℓ2=ℓ for ℓ=5,10,20,40, and solve the sourced Teukolsky equation with the exact complex QNM frequency ω_{ℓ+2} instead of the Eq. (59) identification. If the large-ℓ slope of R_{2×ℓ} shifts by more than the O(1/ℓ) size of the detuning, or if the finite-width source integral departs from the light-ring-localized value by more than a few percent, Eq. (63) is not a controlled leading-order result. Repeat at a/M=0.5 and 0.9 for ℓ1=ℓ2=10 to check whether the plotted spin dependence survives the corrected calculation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central formula (63) is controlled only at leading order in ℓ, yet Eq. (59) identifies the source frequency with the free nonlinear QNM frequency by writing ω_{ℓ1,ℓ1,0}+ω_{ℓ2,ℓ2,0} ≃ (ℓ1+ℓ2)/b0 − iα ≃ ω_{ℓ1+ℓ2,ℓ1+ℓ2,0}. The paper's own Penrose spectrum (B13) gives the fundamental damping as ω_prec/2, with ω_prec = α(r0−M)/(r0+3M); so the sum of two fundamentals has damping ω_prec, whereas the total mode has damping ω_prec/2. The mismatch is O(1) in units of M, i.e. O(1/ℓ) relative to the real frequency. Since the second-order amplitude is controlled by the Green's-function denominator ω1+ω2−ω_{ℓ1+ℓ2}, this O(1) detuning changes the prefactor at the same relative order as the neglected finite-source-width/localization corrections, while the retained iα^4/ω^4 term in (63) is O(ℓ^{−4}). Thus the spin dependence attributed to the light-ring term is not derived at a controlled order. The authors disclose the localization assumption (Sec. V B) and the α-truncation inconsistency (Sec. VI), but they do not quantify the frequency-detuning error; the comparison with the 0.11ℓ numerical result (Sec. VI) sits exactly in this uncontrolled regime.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends to Kerr black holes a Penrose-limit computation of quadratic quasi-normal-mode (QNM) amplitudes that was previously developed for Schwarzschild. The authors take the Penrose limit around equatorial circular null geodesics, solve the linear and second-order Einstein equations in the resulting pp-wave background, match the Weyl scalar Psi_4 to eikonal Teukolsky solutions, and obtain an analytic formula, Eq. (63), for the ratio R_{ell1 x ell2} of the quadratic amplitude to the product of linear amplitudes. The formula contains a subleading i alpha^4/omega^4 term that is meant to encode black-hole spin, together with spin dependence through b0 and the overlap integrals C_ell. The paper plots the ratio as a function of spin and compares the 2 x ell slope with the numerical result ~0.11 ell of Ref. [27].","tokens_in":18694,"tokens_out":22705,"duration_ms":257013,"significance":"If correct, Eq. (63) is the first fully analytic, parameter-free prediction of the spin dependence of quadratic QNM amplitudes in the eikonal limit, a result of direct relevance for ringdown tests with third-generation detectors. The paper's strengths are its self-contained algebraic derivation, the absence of fitted parameters, and the explicit statement of the main assumptions (source localization at the light ring, equal ++/-- amplitudes, retention of a subleading alpha term). The central caveat is that the spin-dependent term is kept at an order below the controlled eikonal approximation, so the quantitative spin dependence is not yet established to the accuracy claimed.","major_comments":[{"comment":"Equation (59) identifies the sum of two linear QNM frequencies with the free nonlinear QNM frequency at leading order, but the paper's own Penrose spectrum, Eq. (B13) with (B14), gives the fundamental QNM as omega_{ell,ell,0} = ell/b0 + i omega_prec/2 + omega_prec/2 (up to the sign convention in the exponential). The sum of two fundamentals therefore has damping omega_prec, whereas the total mode has damping omega_prec/2. This detuning is O(1) in units of M, i.e., O(1/ell) relative to the real part. Because the second-order amplitude is controlled by the Green's-function denominator omega_1 + omega_2 - omega_{ell1+ell2} (equivalently, in the pp-wave variables, by p_u^{src} - p_u^{free}), this O(1) detuning changes the prefactor at the same relative order as the neglected finite-source-width corrections. The retained i alpha^4/omega^4 term in Eq. (63) is O(ell^{-4}) relative to the leading terms, which is far below the controlled order. The authors should either compute the detuning factor explicitly or demonstrate that it cancels in the ratio; as it stands, the spin dependence attributed to the light-ring term is not quantitatively established.","section":"Sec. VI, Eq. (59), and Appendix B"},{"comment":"The assumption that the nonlinear source is localized exactly at the light ring is load-bearing for the final ratio. The WKB turning points in x_1 and x_2 scale as 1/sqrt(ell), so finite-source-width corrections are expected at O(1/ell) in the eikonal expansion. The paper acknowledges the localization assumption qualitatively but does not quantify the error. The comparison with the 0.11 ell numerical result of Ref. [27] in Sec. VI is made at ell = 2, which is outside the eikonal regime, so it cannot validate the controlled-order claim. The authors should provide an estimate of the finite-width correction, or explicitly restrict the prediction to leading order in ell and withdraw the quantitative comparison at ell = 2.","section":"Sec. V B and Sec. VI"},{"comment":"The identity S ~ sum_{ell1,ell2} h_{ell1} h_{ell2} = sum_{ell1>=ell2} sigma_{ell1,ell2} h_{ell1} h_{ell2}, with sigma = 1/2 for ell1=ell2 and sigma = 1 otherwise, appears algebraically incorrect. For unordered pairs, the correct coefficients relative to the ordered double sum are 1 for equal indices and 2 for unequal indices (or one should keep the ordered sum). Since sigma enters linearly in Eq. (58) and hence in Eq. (63), this changes the normalization of the predicted ratio by a factor of two. Please check whether this is a typo in the definition of sigma or whether a different summation convention is intended, and correct it explicitly.","section":"Eq. (40)"}],"minor_comments":[{"comment":"In the expression for S_u2, the first line contains two terms both written with e^{f_-}; from the pattern of S_u1, S_11, and S_22, one of them should presumably be e^{f_+}. Please check and correct this typo.","section":"Appendix A, Eq. (A2)"},{"comment":"The notation omega_{ell1,ell1,0} is confusing; in the standard notation of Eq. (C1) the second index is the azimuthal number m, so for the modes considered here one would write omega_{ell1,m=ell1,n=0} or define the notation explicitly.","section":"Eq. (59)"},{"comment":"The word \"adimensional\" should be \"dimensionless\".","section":"Fig. 2 caption"},{"comment":"The notation for the ratio R_{ell1 x ell2} and for the amplitudes A_{ell1 x ell2}^{ell1+ell2} and A_{ell_i ell_i} is not fully defined; in particular, the connection between the A_ell appearing in Eq. (60) and the A_{ell m} in the TT expansion of Eq. (61) should be stated more explicitly.","section":"Eqs. (60)-(61)"},{"comment":"The statement that the phase stabilizes around exp(-2i pi/5) would be more informative if accompanied by a numerical fit or an analytic estimate; otherwise it is difficult for the reader to assess the convergence in ell.","section":"Sec. VI, Fig. 3 discussion"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a novel and potentially useful method, but the uncontrolled frequency detuning in Eq. (59) and the symmetry-factor issue in Eq. (40) affect the central quantitative claims. If the authors can remove or properly control the subleading alpha term, or provide a detailed error estimate showing that the detuning does not alter the leading-order ratio, the paper could become suitable for publication. As it stands, the spin-dependence claim is stronger than the controlled approximation supports."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things before you open the PDF. This is the first fully analytic Kerr extension of the same group's Schwarzschild Penrose-limit calculation: they get a closed-form ratio for the quadratic-to-linear QNM amplitude in the eikonal limit, with spin entering through the light-ring quantities and the C_ell overlap integrals. And the paper is careful about disclosing its main assumptions. Both are to their credit.\n\nThe problem is that the control of the approximation is not where it appears to be. The frequency-matching step in Eq. (59) is the load-bearing soft spot. The paper's own spectrum (B13) gives the fundamental damping as ω_prec/2 per mode. Two fundamentals therefore have damping ω_prec, while the total mode has damping ω_prec/2. The detuning is O(ω_prec) ~ O(1/M), which is O(1) in the denominator that controls the second-order amplitude. That is not a small correction; it is at the same level as the finite-source-width/localization effects the paper says it is ignoring. Meanwhile the spin-dependent term they retain in (63), iα^4/ω^4, is O(ℓ^-4) relative to the leading bracket, far below the O(1/ℓ) relative uncertainty introduced by the detuning. So the quantitative spin dependence in Figure 2 is not controlled at the claimed order.\n\nThe authors are forthright about the light-ring localization and the inconsistent truncation of α, and they flag the comparison with the numerical 0.11ℓ result as outside the eikonal regime. That honesty is real. But the frequency-detuning error is not listed among the caveats, and it directly affects the magnitude of the ratio, not just a subleading phase.\n\nWhat is genuinely new and useful: the C_ell projection onto spin-weighted spheroidal harmonics, the analytic handling of the ℓ1×ℓ2 source phases, and the demonstration that the Penrose-limit method extends to Kerr. The paper is self-contained, has no fitted parameters, and is written so the approximations are easy to audit. For a reader working on nonlinear ringdown, especially in preparation for third-generation detectors, it deserves a serious read.\n\nMy recommendation: send it to peer review, but the referee should push hard on the frequency matching. The framework is likely salvageable at leading order, but the spin dependence needs to be either derived at a controlled order or explicitly labeled as a heuristic extrapolation.","headline":"A clean but uncontrolled spin-dependent extension of the Schwarzschild Penrose-limit calculation; the frequency-detuning issue in Eq. (59) is real and needs fixing before the spin dependence can be trusted.","tokens_in":19250,"tokens_out":8543,"would_cite":false,"duration_ms":87157,"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":"The paper derives a fully analytic, spin-dependent formula for the amplitude and phase of quadratic quasi-normal modes in Kerr ringdown, in the eikonal limit, using the Penrose limit around the photon ring.","keywords":["Kerr black holes","ringdown","quasi-normal modes","quadratic quasi-normal modes","Penrose limit","eikonal limit","gravitational waves","black hole spin"],"falsifier":"Numerically solve the second-order Teukolsky equation for Kerr with a source whose radial profile is a narrow but finite shell centered on the photon ring, and compare the quadratic quasi-normal-mode amplitude ratio for, say, $\\ell_1=\\ell_2=10$ and $a/M=0.9$ with Eq. (63); if the ratio moves by more than the retained $\\alpha^4/\\omega^4$ correction as the shell width is varied, the localization assumption is disproved.","tokens_in":18153,"feed_emoji":"🕳️","tokens_out":10418,"duration_ms":108898,"temperature":0.7,"pith_summary":"Ringdown gravitational waves from a Kerr black hole contain not just the linear quasi-normal modes but also 'quadratic' modes at sums and differences of the linear frequencies, and these nonlinearities are a target for future gravitational-wave detectors. This paper tries to compute them analytically in the eikonal, large-angular-momentum limit. It works in the Penrose limit around the equatorial circular null geodesic, where the Kerr metric becomes a simple plane-wave metric, solves the linear and second-order perturbation equations there, and matches the result to the asymptotic Teukolsky solution at null infinity. The payoff is a closed-form, spin-dependent expression for the ratio of the quadratic amplitude to the linear amplitude squared, including its phase. If correct, it gives an analytic prediction that numerical relativity and future detectors can check directly.","feed_headline":"Formula predicts Kerr ringdown's nonlinear amplitude and phase","feed_subtitle":"The eikonal ratio matches numerical growth and ties spin dependence to photon-ring geometry.","key_machinery":"The engine is the Penrose limit taken along the equatorial circular null geodesic of Kerr. Around that geodesic the metric becomes the plane-wave metric $ds^2 = 2\\,du\\,dv + \\alpha^2(x_1^2-x_2^2)\\,du^2 + dx_1^2 + dx_2^2$, with $\\alpha$ encoding the black hole spin through $\\alpha^2 = 12M\\Delta/[r_0^3(r_0-M)^2]$, equivalently $\\alpha^2=(b_0^2-a^2)/r_0^4$. On this background the linear perturbation separates into harmonic-oscillator equations in the transverse coordinates, giving the eikonal quasi-normal-mode frequency $p_v=\\ell/b_0$ and damping set by $\\alpha$; the second-order perturbation is then solved as a sourced version of the same equations. The pp-wave solution is matched to the eikonal Teukolsky solution through WKB principal functions, and the resulting ratio is projected onto spin-weighted spheroidal harmonics through the overlap integral $C_\\ell$. The nonlinear input is the assumption, taken from the peaked structure of the second-order Teukolsky source, that the source is localized exactly at the light ring.","core_discovery":"The central result is the eikonal-limit amplitude ratio for the quadratic quasi-normal mode produced by two linear modes with angular momenta $\\ell_1$ and $\\ell_2$ on a Kerr background. The paper obtains $$R_{\\ell_1\\times\\ell_2} = \\left[\\frac{\\$ell_1^{3}$\\ell_2+\\ell_1\\$ell_2^{3}$ - i\\sigma_{\\ell_1,\\ell_2}(\\ell_1+\\ell_2)^4 + i(\\$alpha^{4}$/\\$omega^{4}$)}{8\\ell_1\\ell_2(\\ell_1+\\ell_2)^2}\\right]\\frac{C_{\\ell_1+\\ell_2}}{C_{\\ell_1}C_{\\ell_2}},$$ where $\\sigma_{\\ell_1,\\ell_2}=1$ for $\\ell_1\\neq\\ell_2$ and $1/2$ for equal modes, $\\alpha$ is the Penrose-limit frequency scale, $\\omega\\simeq\\ell/b_0$, and $C_\\ell$ projects the WKB angular phase onto spin-weighted spheroidal harmonics. Spin enters through $\\alpha$, the impact parameter $b_0$, the photon-ring radius $r_0$, and the overlap integrals $C_\\ell$. The paper reports that for $\\ell_1=\\ell_2=\\ell$ the amplitude grows as $\\ell^{1/4}$, that for a $2\\times\\ell$ mode it grows roughly as $0.08\\ell$ (compared with a numerical result of roughly $0.11\\ell$), that the phase approaches $\\exp(-2i\\pi/5)$ at large $\\ell$, and that the prograde orbit responds much more strongly to spin than the retrograde orbit.","pith_inferences":["A natural next step the paper does not take is to use Eq. (63) as the large-$\\ell$ anchor for a resummed or next-to-leading-order prediction at $\\ell=2$, the multipole most relevant for current detectors.","Because $\\alpha$, $b_0$, and $r_0$ are all photon-ring quantities, the result implies that the spin dependence of quadratic ringdown is fixed by the photon-ring geometry; a numerical experiment that altered the photon ring while keeping the quasi-normal-mode spectrum fixed would sharply test that localization.","The same machinery could be pushed to subleading overtones or to the difference-frequency $g_\\pm$ modes, yielding analytic predictions for the full nonlinear peak structure rather than just the sum-frequency peak.","The mild spin dependence found for the retrograde orbit and the stronger dependence for the prograde orbit suggest future detectors might use the ratio's spin trend, rather than its absolute value, as a cleaner probe of black-hole spin during ringdown."],"forward_implications":["Equation (63) gives a closed-form amplitude and phase for the quadratic quasi-normal mode at any $\\ell_1,\\ell_2$ in the eikonal limit, with the spin dependence carried by $\\alpha$, $b_0$, and the overlap integrals $C_\\ell$.","For equal modes $\\ell_1=\\ell_2=\\ell$, the predicted amplitude scales as $\\ell^{1/4}$ and is only mildly spin-dependent, while the phase settles near $\\exp(-2i\\pi/5)$ for large $\\ell$.","For a $2\\times\\ell$ combination, the predicted ratio grows roughly as $0.08\\ell$, in reasonable agreement with the numerical value of about $0.11\\ell$ even though $\\ell=2$ lies outside the eikonal regime.","The prograde orbit is much more sensitive to spin than the retrograde orbit because its impact parameter shrinks sharply as $a/M\\to 1$, so the spin dependence of the ratio is mostly a photon-ring effect.","If the symmetry assumption $A_{++}=A_{--}$ is relaxed, the nonlinear ratio becomes dependent on the initial excitation amplitudes of the two polarizations, as the paper shows by keeping distinct amplitudes for the two mirror modes."],"supporting_citations":[{"why":"Gives the Schwarzschild Penrose-limit nonlinear-ratio calculation that this paper extends to Kerr and supplies the matching procedure.","marker":"[44]"},{"why":"Shows the Penrose limit around a closed null geodesic reproduces the Kerr quasi-normal-mode spectrum and fixes the frame and $\\alpha^2$ used here.","marker":"[43]"},{"why":"Provides the eikonal-limit Teukolsky and Hamilton-Jacobi solution and the principal functions used to match the pp-wave solution to asymptotic Kerr waves.","marker":"[40]"},{"why":"Argues the second-order source is peaked near the light ring, the assumption that lets the authors localize the source at $r_0$.","marker":"[56]"},{"why":"Earlier evidence that the nonlinear source in the wave equation is concentrated near the light ring, supporting the same localization.","marker":"[57]"},{"why":"Numerical result for the $2\\times\\ell$ amplitude ratio (about $0.11\\ell$) that the paper compares with its eikonal prediction (about $0.08\\ell$).","marker":"[27]"},{"why":"Shows mirror-mode excitations depend on initial conditions, motivating the separate $A_{++}$ and $A_{--}$ amplitudes kept in the calculation.","marker":"[22]"},{"why":"Suppresses the $f_+$ mirror-mode contribution at large $\\ell$, so the paper keeps only the $f_-$ term in the final ratio.","marker":"[58]"}],"fun_headline_variants":["Penrose limit yields analytic Kerr ringdown nonlinearity","Formula predicts spin-dependent Kerr ringdown nonlinearity","Eikonal method computes quadratic quasinormal modes exactly","Kerr ringdown's quadratic modes calculated via Penrose limit","Amplitude ratio for Kerr ringdown from photon-ring geometry"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the nonlinear source is localized exactly at the light ring and that only the spin-dependent $\\alpha$ term needs to be kept beyond leading order; if the source has finite radial width, the predicted spin dependence changes.","fun_headline_variants_meta":{"raw":{"variants":["Penrose limit yields analytic Kerr ringdown nonlinearity","Formula predicts spin-dependent Kerr ringdown nonlinearity","Eikonal method computes quadratic quasinormal modes exactly","Kerr ringdown's quadratic modes calculated via Penrose limit","Amplitude ratio for Kerr ringdown from photon-ring geometry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000421,"raw_usage":{"total_tokens":2173,"prompt_tokens":963,"completion_tokens":1210,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":579,"completion_tokens_details":{"reasoning_tokens":1144}},"tokens_in":579,"tokens_out":1210,"duration_ms":13112,"temperature":1.0,"reasoning_tokens":1144,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:41:24.550739+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically solve the second-order Teukolsky equation for Kerr with a source whose radial profile is a narrow but finite shell centered on the photon ring, and compare the quadratic quasi-normal-mode amplitude ratio for, say, $\\ell_1=\\ell_2=10$ and $a/M=0.9$ with Eq. (63); if the ratio moves by more than the retained $\\alpha^4/\\omega^4$ correction as the shell width is varied, the localization assumption is disproved.","supporting_citations":[{"cited_title":"overtone","cited_arxiv_id":null,"evidence_quote":"Gives the Schwarzschild Penrose-limit nonlinear-ratio calculation that this paper extends to Kerr and supplies the matching procedure."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the eikonal-limit Teukolsky and Hamilton-Jacobi solution and the principal functions used to match the pp-wave solution to asymptotic Kerr waves."},{"cited_title":"Igata, H","cited_arxiv_id":null,"evidence_quote":"Argues the second-order source is peaked near the light ring, the assumption that lets the authors localize the source at $r_0$."}],"review_version":1}