Pith. sign in

REVIEW 3 major objections 6 minor 2 cited by

Extremal Black Hole Weather

T0 review · 3 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash

Pith's one-line read Near-extremal Kerr may have gravitational 'weather' with dyadically suppressed high-ell modes

desk verdict A serious, novel derivation of a QNM-amplitude dynamical system for near-extremal Kerr; the dyadic equilibrium spectrum is suggestive but rests on an untested regulator assumption. read the letter →

arxiv 2412.02821 v1 pith:ANAAYF5N submitted 2024-12-03 gr-qc

classification gr-qc PACS 04.70.-s04.30.-w
keywords near-extremalKerrzero-dampedquasinormalmodesHertzpotentialsecond-orderEinsteinequationinversecascadeweakturbulencenNHEKlimitblackholeweather
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper asks whether a sufficiently rapidly spinning (near-extremal) Kerr black hole can support a weakly nonlinear 'weather' regime, in which long-lived, almost-resonant gravitational modes exchange energy before they decay. Starting from the second-order vacuum Einstein equation and a Hertz-potential description of metric perturbations, it derives an infinite-dimensional dynamical system for the amplitudes of the zero-damped quasinormal modes. In the near-horizon scaling limit, with angular momenta split into a low-$\ell$ and a high-$\ell$ eikonal sector, the system admits a time-independent equilibrium supported on axisymmetric, zero-overtone modes with amplitudes $c_\ell \sim C^{\rm low}\,2^{-\ell/2}\ell^{-7/2}$. The authors interpret this as evidence that weakly nonlinear self-interaction alone drives the perturbation toward a smooth, large-scale endpoint: an inverse cascade, with the low-$\ell$ modes acting as a pump. If the interpretation survives, it gives a gravitational counterpart of the 'weather' patterns seen in simplified models of rotating atmospheres, with no external driving and no reliance on the linear decay of the modes.

What carries the argument

The machinery is a QNM-amplitude dynamical system built for near-extremal Kerr. In the corrector-tensor scheme a metric perturbation is encoded by a Hertz potential $\Phi$, and the second-order Einstein equation becomes a sourced Teukolsky equation for $\Phi$; a conserved bilinear product for Teukolsky solutions converts the sourced equation into $\dot c_1=\alpha\sum_{2,3}(U_{123}c_2c_3+V_{123}c_2c_3^*)$. The coefficients $U_{123},V_{123}$ are three-wave overlap integrals. In the nNHEK scaling limit (the near-near-horizon extremal Kerr limit), with $\ell$ restricted to $m^2\ll\ell\ll\varepsilon^{-1}$, they factor into angular integrals over spin-weighted spherical harmonics and regulated radial triple products; the large-$\ell$ angular integrals are evaluated by stationary phase, yielding triangle selection rules, and the radial integrals reduce to three-point objects of the near-horizon $\mathrm{SL}(2,\mathbb{R})$ symmetry. The equilibrium is found by rescaling $\ell=L\bar\ell$, passing from sums to integrals, and solving the resulting scale-invariant equation by the power-law ansatz $c_{\bar\ell}\propto\bar\ell^{-2}$; after inverting the amplitude renormalization this gives the $2^{-\ell/2}\ell^{-7/2}$ spectrum.

What would settle it

A concrete test is to integrate the truncated dynamical system numerically from generic initial data for a finite but large mode set (the paper states this is in principle possible) and check whether high-$\ell$ amplitudes approach $c_\ell\propto 2^{-\ell/2}\ell^{-7/2}$ or grow instead. A second check is to include modes with $\ell\sim\varepsilon^{-1/2}$ and $\ell\gtrsim\varepsilon^{-q}$: if the regularized double integral's value, and hence the exponent $p=-2$, changes when those sectors are added, the equilibrium is an artifact of the truncation.

Watch

Extended reading notes

Core claim

On its own terms, the paper claims that the leading self-interaction of gravitational perturbations around near-extremal Kerr has a nontrivial stationary state in the space of zero-damped QNM amplitudes. The state is supported only on axisymmetric modes ($m=0$) and zero overtones ($N=0$), and at large polar index it obeys the dyadic power law $c_\ell\propto C^{\rm low}2^{-\ell/2}\ell^{-7/2}$, where $C^{\rm low}$ is a weighted sum of the low-$\ell$ amplitudes that pumps the high-$\ell$ sector. Because the linear damping factors are essentially frozen during the relevant Boyer-Lindquist time, the suppression is produced entirely by quadratic nonlinearities. The corresponding metric, the Kerr background plus the reconstructed Hertz-potential perturbation and the corrector tensor, is time-independent over a parametrically long interval that diverges as extremality is approached. The paper does not claim to prove that this equilibrium is an attractor; the claim is that it is a solution of the derived dynamical system and evidence for a high-$\ell$ inverse cascade.

Load-bearing premise

The equilibrium result depends on a specific way of cutting off divergent mode sums, which is an unproven assumption about how the omitted middle-size and very large angular-momentum modes behave; if that assumption is wrong, the claimed spectrum is not established.

Editorial extensions

If this is right

  • If the equilibrium is reached, the endpoint is a smooth near-Kerr geometry with no large-$\ell$ content, so gravitational radiation in this regime would be dominated by low-$\ell$ axisymmetric modes rather than a turbulent spread to high harmonics.
  • The equilibration time can be parametrically long in $\varepsilon$: the dynamical system is written in slow time $\bar t=\varepsilon t/(2M)$, and the equilibrium persists for Boyer-Lindquist times of order $M\varepsilon^{q-1}$ with $1/2<q<1$.
  • Linear decay of the zero-damped QNMs is not responsible for the suppression; the result is due to the quadratic self-interaction of the modes alone.
  • The explicit overlap coefficients for the (high,high)$\to$(high), (high,low)$\to$(high), and related channels give a quantitative replacement for earlier truncated two-mode models of resonant mode coupling, treating all QNMs on the same footing.
  • Higher-order cubic interactions are not needed for the leading cascade, but the authors note that resonant quartets are generic and would be the next correction to this picture.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The existence of the equilibrium does not by itself establish a cascade; the paper's interpretation is conditional on the equilibrium being attracting, and a decisive numerical experiment integrating the finite truncation from generic data would turn the 'weather' conjecture into a dynamical claim.
  • The equilibrium is restricted to $m=0$, and the authors explicitly note that reducing dimensionality can change cascade direction; whether non-axisymmetric modes are depleted or fed by the axisymmetric cascade is the first question to resolve.
  • The dimensional-renormalization cutoff assumption about intermediate and ultra-large-$\ell$ modes is a physical claim about how those sectors are slaved to the low/high sectors; a cleaner check would compute next-order corrections to the overlap coefficients and see whether the $p=-2$ cancellation survives.
  • The same three-wave machinery could be adapted to other long-lived resonant mode families, such as near-extremal charged black holes, which would show whether dyadic suppression is a generic feature of near-extremal horizons or specific to Kerr.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. This paper develops a formalism for weakly nonlinear gravitational perturbations of a near-extremal Kerr black hole at second order in the vacuum Einstein equation. Using the GHZ metric-reconstruction framework, the perturbation is represented through a Hertz potential and expanded in zero-damped quasinormal modes with time-dependent amplitudes; projecting the sourced Teukolsky equation with the conserved bilinear form of Green-Hollands-Sberna-Toomani-Zimmerman yields an infinite-dimensional quadratic dynamical system (Eq. 92). In the near-extremal nNHEK scaling limit, with angular momenta split into low (O(1)) and high (m^2 << ell << epsilon^{-1}) sectors, the overlap coefficients are evaluated in matched asymptotic expansions into explicit factorized angular and radial integrals (Sec. VI and Apps. H-M). The paper then seeks an equilibrium of the truncated system in the sector m = 0, N = 0, odd ell. Using a scaling ansatz c_ellbar ~ ellbar^p and replacing the finite integration bounds of Eq. (167) by (0, infinity) with a 'dimensional renormalization prescription,' the two nonlinear terms are claimed to cancel at p = -2, giving an equilibrium amplitude spectrum c_ell ~ C_low 2^{-ell/2} ell^{-7/2} (Eqs. 25 and 169-171). The authors interpret this suppression of large-ell axisymmetric amplitudes as an inverse cascade, which they call 'black hole weather.'

Significance. If the equilibrium claim survives scrutiny, this is a potentially important step beyond the heuristic single-triad model of Yang, Zimmerman and Lehner: it would be the first derivation from the vacuum Einstein equation of a concrete large-ell spectrum for zero-damped QNM self-interactions in near-extremal Kerr. The paper's strengths are real: the dynamical system (92) is derived from the second-order Einstein equation rather than postulated; the overlap coefficients are computed explicitly in a matched asymptotic expansion with substantial appendices; the ladder-operator structure of Þ and Þ' and the SL(2,R)-module interpretation (Apps. F, L) are elegant; the asymptotic formulas for angular integrals are checked numerically (App. H); and the authors are candid about every truncation, including the non-self-consistency of the N = 0 restriction and the regulator dependence of the equilibrium. The truncated system is explicit enough that the equilibrium claim is numerically falsifiable, as the authors themselves note (Sec. VIII).

major comments (3)
  1. [Sec. VII B, Eqs. (167)-(169).] The equilibrium condition must be satisfied by the actual truncated dynamical system, i.e., by Eq. (167) with the finite integration domain c <= ellbar_i <= epsilon^{1/2-q} (L ~ epsilon^{-1/2}, 1/2 < q < 1). The text replaces these finite boundaries by (0, infinity), states that the two nonlinear terms cancel for the ansatz c_ellbar ~ ellbar^p at p = -2, and regulates the resulting divergences by a 'dimensional renormalization prescription' that the paper itself describes as 'an implicit assumption on the nature of the intermediate QNMs having ell <~ epsilon^{-1/2} and QNMs with very large ell >~ epsilon^{-q}.' This is load-bearing: for p = -2 the double integral is divergent at both endpoints, so the infinite-boundary cancellation is not a controlled limit of the finite-domain integral, and the change of domain is precisely an assumption about modes for which the approximations leading to Eq. (167) are not valid. The equation that must vanish for an equilibrium is the finite-boundary version; unless its residual is shown to vanish, or to be subdominant relative to the ellbar_1 c_high,1 Re(C_low) term in Eq. (167), the cancellation fixing p = -2, and with it the dyadic spectrum (25) and (169), is imposed by the regulator rather than established for the retained modes. Please provide one of the following: (i) an estimate or numerical evaluation of the finite-boundary residual at p = -2; (ii) a demonstration that the cancellation is independent of the renormalization scheme, with the continuation parameter and the subtraction prescription specified; or (iii) a numerical solution of a finite but large truncation of Eq. (167) showing a stable fixed point with the claimed scaling.
  2. [Secs. IV B and VI A.] The dynamical system (131) rests on two truncations that are load-bearing for the central claim: the retarded Green's function is replaced by its QNM part (G_ret ~ G_qnm, with G_cut and G_arc neglected, Eqs. (68)-(70)), and the far-zone contributions to the overlap coefficients are discarded on the basis of the 'central hypothesis' in Sec. VI A that they are parametrically small in epsilon. The supporting evidence, namely the decay of the near-zone radial integrands as xbar -> infinity (App. I) and the analogy with the scalar products computed in App. D, is suggestive but does not quantify the neglected terms. Because the equilibrium solution lives entirely in the near-zone QNM sector, any far-zone or cut/arc contribution at the same order as the retained nonlinear terms would modify the equilibrium. I request a quantitative estimate, or a targeted numerical check in the spirit of the numerical tests of the angular asymptotics reported in App. H, for the far-zone portions of at least the dominant (high,high) -> (high) overlap coefficients U^near_123 and V^near_123.
  3. [Sec. VII A.] The equilibrium is solved within the sector m = 0, N = 0, odd ell, and the paper states that the N = 0 restriction is 'not self-consistent (a pair of N = 0 QNMs can excite an N > 0 QNM)' and that the odd-ell restriction is made for convenience. App. J shows that the N_S = 2 and N_S = 4 overlap coefficients share the homogeneous scaling structure of the N_S = 0 coefficients, which supports the authors' conjecture that the N = 0 restriction is inessential, but the equilibrium condition (164) is not checked against the N > 0 or even-ell channels. Since the physical interpretation, that weakly nonlinear dynamics drives the system toward an axisymmetric, zero-overtone, dyadically suppressed spectrum, is a statement about the full QNM system, the present result is an equilibrium of a heavily truncated subsystem. Please either extend the equilibrium analysis to include the leading N > 0 and even-ell channels, or argue their decoupling, or state explicitly in the abstract and conclusions that the equilibrium and the cascade interpretation are established only for the N = 0, m = 0, odd-ell subsystem.
minor comments (6)
  1. [Abstract and Sec. I D.] The term 'dyadically exponentially suppressed' is potentially misleading: a dyadic cascade conventionally denotes a discrete hierarchy of scales with ratio two (e.g., ell -> 2ell), whereas the result (25) is an exponential decay in ell with base 2^{-1/2}. Please clarify the intended meaning or justify the terminology.
  2. [Abstract and Sec. VII.] The abstract states that 'our dynamical system has a time-independent solution' and reads the result as evidence for an inverse cascade. As Secs. VII A-B make clear, the statement holds only for the truncated system (N = 0, m = 0, odd ell), inside a finite time window t <~ M epsilon^{q-1}, and after the renormalization prescription applied to Eq. (167); the abstract should carry these qualifications.
  3. [Sec. VII B, footnote 28.] Footnote 28 recommends passing to a medium time ttilde = epsilon^{-q} tbar to describe the validity window of Eq. (167); this point is essential for interpreting the equilibrium and should be promoted to the main text.
  4. [App. J and App. D.] In the paragraph below Eq. (J2), the label 'F NS=2_123' should read 'F NS=4_123' (the polynomial defined by Eq. (I16) belongs to N_S = 4), and in App. D, just above Eq. (D6), the mode label '-2Upsilon_N2ell2m3' should be '-2Upsilon_N2ell2m2'.
  5. [Secs. II and V.] The symbol x is used both for the scaled near-zone coordinate (29) and for the far-zone coordinate in Eq. (113) and in the overlap and matching conditions (99)-(100); a compact table or a change of symbol for one of the two coordinates would remove a recurring source of confusion.
  6. [Sec. IV B, Eq. (70).] The compact representation of G_qnm in Eq. (70) is stated to follow from [41, Lem. 5], and it is the linchpin of the projection argument leading to Eqs. (84)-(88); a one-sentence statement of the lemma's content and normalization conventions would help the reader verify the sign and the factor 1/A_q in Eq. (84).

Circularity Check

0 steps flagged · score 2.0 of 10

No construction-level circularity: Eq. (25) is a genuine fixed point of the Einstein-derived system (92), with p=-2 selected by a cancellation condition rather than fitted or imposed; the flagged Sec. VII B regulator step and the imposed m=0/N=0 truncations are validity risks the paper itself discloses, not by-construction inputs.

full rationale

The derivation chain is: second-order vacuum Einstein equation -> sourced Teukolsky equation (81) via the GHZ decomposition (refs 24, 35, 40) -> projection with the conserved scalar product of ref 41 -> dynamical system (92) -> nNHEK/eikonal overlap coefficients (134, 159, App. K) -> truncated axisymmetric, N=0, odd-ell system (167) -> equilibrium ansatz (168) with exponent p fixed by the cancellation condition -> scaling (169) -> spectrum (25)/(170) via normalizations (156, D25). No claimed result is inserted as an input: the equilibrium is defined as a fixed point of (164) and solved, not assumed; p is "to be determined" (Sec. VII B) and is selected, as the paper states, by the requirement that the two terms in the second equation of (167) cancel for a proportionality constant of order C_low; there is no dataset, fitted parameter, or uniqueness theorem that could force (25). The central caveat, Sec. VII B, is that the finite domain (c, epsilon^{1/2-q}) is replaced by (0,infinity) and the log-divergent double integral is evaluated by a "dimensional renormalization prescription" that the authors themselves call "an implicit assumption on the nature of the intermediate QNMs having ell ~ epsilon^{-1/2} and QNMs with very large ell >= epsilon^{-q}, for which the approximations used to arrive at (167) do not necessarily hold." If excluded modes behave differently, the cancellation fixing p=-2 fails; this is a robustness/validity gap, honestly disclosed, not an equation-level reduction of output to input, since the regulator is not constructed to produce the spectrum and the paper does not claim the prescription is derived. The imposed restrictions m=0 and N=0 feed the delta_{m,0} delta_{N,0} support of (25); the paper explicitly says the N=0 restriction "is not self-consistent (a pair of N=0 QNMs can excite an N>0 QNM), but is made to arrive at a manageable truncated dynamical system," and Sec. VIII concedes that no attractor proof is given and that axisymmetry may affect cascade direction. These are truncations whose consequences are weighed in the paper, not renamed derivations; the genuinely novel ell-dependence is computed, and App. J's extension to N>0 is offered as conjecture.

Assumptions & free parameters 4 free parameters · 7 assumptions · 1 invented entities

The central result rests on a nested set of approximations: QNM dominance of the Green's function, neglect of far-zone contributions, decoupling of intermediate and very large ell modes, a power-law ansatz for the equilibrium, and a dimensional renormalization of divergent integrals. The first three are domain assumptions typical of weak turbulence analyses; the last two are ad hoc to this paper and are the least externally anchored. No fundamentally new physical entity is introduced beyond the conceptual label rQQNMs.

free parameters (4)
  • C_low (cumulative low-ell amplitude) = undetermined (sets overall scale of equilibrium)
    The equilibrium solution (25), (169) is proportional to C_low, a weighted sum of low-ell amplitudes that is not fixed by the dynamical system; it is a free normalization set by initial data.
  • Truncation scale L (low/high ell split) = L ~ epsilon^{-1/2} (with L of order 10^2-10^3 in practice)
    The split into low and high ell modes and the neglect of intermediate ell (Sec. VI B, VII A) relies on a chosen scale L ~ epsilon^{-1/2}; the results are claimed independent of L at leading order.
  • Upper cutoff exponent q (1/2 < q < 1) = chosen in (1/2,1), e.g. q = 3/4
    The high-ell sector is summed up to ell less than about epsilon^{-q}; the equilibrium integral cancellation assumes the domain (c, epsilon^{1/2-q}) is effectively (0, infinity) after renormalization.
  • Nonlinearity strength alpha = set to 1 in Eq. (167) ('an assumption about the overall amplitude')
    The overall amplitude of the perturbation is absorbed into alpha; the equilibrium solution is stated for alpha = 1.
assumptions (7)
  • domain assumption G_ret is dominated by its QNM contribution; branch cut and direct-absorption pieces are negligible in the regime of interest.
    Introduced in Sec. IV B ('we shall generally approximate Gret by the QNM piece Gqnm'); this is the key reduction from the full Einstein equation to a QNM amplitude system.
  • domain assumption Far-zone contributions to the overlap coefficients U123, V123 are parametrically small in epsilon.
    Stated as 'our central hypothesis' in Sec. VI A; supported only by heuristic decay arguments in Apps. D and I, not proven.
  • domain assumption Intermediate and very large ell QNMs decouple and can be discarded; only ell = O(1) and m^2 much less than ell much less than 1/epsilon are kept.
    Sec. VI B assumes this in the spirit of weak wave turbulence; the paper notes 'the hypothesis is thereby that they effectively decouple'.
  • ad hoc to paper Divergent equilibrium integrals can be regulated by dimensional renormalization, with the scheme encoding the unknown behavior of intermediate QNMs.
    Sec. VII B: the cancellation that fixes p = -2 relies on treating the integration domain as (0, infinity) and regulating divergences by a prescription that is 'an implicit assumption on the nature of the intermediate QNMs'.
  • ad hoc to paper The equilibrium is sought within the m = 0, N = 0, odd-ell sector.
    Sec. VII B: these restrictions are made 'to arrive at a manageable truncated dynamical system'; the N=0 restriction is explicitly not self-consistent.
  • ad hoc to paper The ansatz c_bar_ell proportional to bar_ell^p with p = -2 for the equilibrium.
    The power-law ansatz (168) is assumed; the specific value p = -2 is then fixed by the cancellation condition, so only the functional form is an input.
  • standard math Matched asymptotic expansion gives the QNM frequencies and mode functions to the stated order.
    Sec. V uses the standard matched asymptotic expansion (initiated by Teukolsky-Press, developed by Hod, Yang et al.); parts are corroborated by rigorous work [74-77] but the O(epsilon) errors are only estimated.
invented entities (1)
  • Resonant quadratic QNMs (rQQNMs)
    purpose: A new class of second-order modes that are driven on resonance in near-extremal Kerr and organized as a self-consistent QNM amplitude dynamics.
    Introduced in Sec. I E as a complement to ordinary quadratic QNMs; no falsifiable prediction outside the paper is given. It is a conceptual label for the resonant regime studied here.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Extremal Black Hole Weather." pith.science (2026). https://pith.science/paper/ANAAYF5N

@misc{pith2026241202821,
  author       = {Pith},
  title        = {Pith review of: Extremal Black Hole Weather},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ANAAYF5N}},
  note         = {Machine review of arXiv:2412.02821}
}
abstract

We consider weakly non-linear gravitational perturbations of a near-extremal Kerr black hole governed by the second order vacuum Einstein equation. Using the GHZ formalism [Green et al., Class. Quant. Grav. 7(7):075001, 2020], these are parameterized by a Hertz potential. We make an ansatz for the Hertz potential as a series of zero-damped quasinormal modes with time-dependent amplitudes, and derive a non-linear dynamical system for them. We find that our dynamical system has a time-independent solution within the near horizon scaling limit. This equilibrium solution is supported on axisymmetric modes, with amplitudes scaling as $c_\ell \sim C^{\rm low} 2^{-\ell/2} \ell^{-\frac{7}{2}}$ for large polar angular momentum mode number $\ell$, where $C^{\rm low}$ is a cumulative amplitude of the low $\ell$ modes. We interpret our result as evidence that the dynamical evolution will approach, for a parametrically long time as extremality is approached, a distribution of mode amplitudes dyadically exponentially suppressed in $\ell$, hence as the endpoint of an inverse cascade. It is reminiscent of weather-like phenomena in certain models of atmospheric dynamics of rotating bodies. During the timescale considered, the decay of the QNMs themselves plays no role given their parametrically long half-life. Hence, our result is due entirely to weakly non-linear effects.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Schr\"odinger perturbation theory for black hole quasinormal modes

    gr-qc 2026-07 conditional novelty 7.0 of 10

    A bilinear-form framework computes black-hole quasinormal-mode frequency shifts to any order, but the mode-sum expansion of the first-order mode shift diverges and needs a continuum piece.

  2. Gravitational Enstrophy: Local Geometric Origin and Inverse-Cascade Constraints

    gr-qc 2026-08 conditional novelty 5.0 of 10

    General relativity admits a magnetic Weyl enstrophy whose approximate conservation, together with energy, enforces a gravitational Fjørtoft constraint biasing nonlinear energy transfer toward lower frequencies in near...

Reference graph

Works this paper leans on

116 extracted references · 46 canonical work pages · cited by 2 Pith papers

  1. [1]

    The conceptual reason for this is that the branch cut starting at ω = 0 in the complex frequency plane characterizing the contribution in the retarded Green’s function (see Sec

    Axisymmetric QNMs The treatment of QNMs in the nNHEK approximation ε ≪ 1 is qualitatively different, and in fact more subtle, for axisymmetric (m = 0) modes than for m ⁄= 0, as has 31 been noted in several works before, see e.g., [18, 93]. The conceptual reason for this is that the branch cut starting at ω = 0 in the complex frequency plane characterizing ...

  2. [2]

    We assume vanishing overtone numbers, i.e., Ni = 0

    (high),(high) → (low) In the following formulas, ℓ1 corresponds to the low ℓ, whereas ℓ2,ℓ 3 are large e.g., ¯ℓ2, ¯ℓ3 ≥ c > 0, in terms of the rescaled angular momenta ¯ℓ = ℓ/L. We assume vanishing overtone numbers, i.e., Ni = 0. 46 U near 123 =δm1,m2+m3 i−ℓ2−ℓ3(−1)m1π 5 2 Γ( −1 +h1 −im1)¯ℓ2 1 2 ¯ℓ3 1 2 csc [ π 2 (h1 +ℓ2 +ℓ3 + 2im1) ] 4 √ 2+2C0ℓ1m1M 4 3 Γ...

  3. [3]

    We assume vanishing overtone numbers, i.e., Ni = 0

    (high),(low) → (low) In the following formulas, ℓ1 and ℓ3 correspond to the low ℓ’s, whereas ℓ2 is large e.g., ¯ℓ2 ≥ c > 0, in terms of the rescaled angular momenta ¯ℓ = ℓ/L. We assume vanishing overtone numbers, i.e., Ni = 0. 47 U near 123 = −δm1,m2+m3 i−ℓ2(−1)m2π 7 4 Γ( −1 +h1 −im1) sec [ 1 2π(h1 +h3 +ℓ2 + 2im1) ] 2 3 4M 5 3 ¯ℓ2 2+2C0ℓ1m1Γ( −2 −h1 −im1)...

  4. [4]

    Lagos and L

    M. Lagos and L. Hui, Physical Review D 107, 044040 (2023)

  5. [5]

    Sberna, P

    L. Sberna, P. Bosch, W. E. East, S. R. Green, and L. Lehner, Physical Review D 105, 064046 (2022)

  6. [6]

    Bourg, R

    P. Bourg, R. P. Macedo, A. Spiers, B. Leather, B. Bonga, and A. Pound, arXiv preprint arXiv:2405.10270 (2024)

  7. [7]

    minimal subtraction scheme

    in nNHEK, we have e.g., λI =i. (C6) The Teukosky operators sO,s ≥ 0 and their adjoints have useful intertwining relations with I, J . For the case of J , this intertwining relation is [41] sOΨ 2s 3 2 J = Ψ 2s 3 2 J sO† (C7) where Ψ 2 is the non-vanishing background Weyl scalar of Kerr. For the case of I, this intertwining relation is sO†Ψ − 2s 3 2 I = Ψ −...

  8. [8]

    Andersson and K

    N. Andersson and K. Glampedakis, Physical Review Letters 84, 4537–4540 (2000)

Show all 116 references
  1. [9]

    Hod, Phys

    S. Hod, Phys. Lett. B 715, 348 (2012) , arXiv:1207.5282 [gr-qc]

  2. [10]

    H. Yang, F. Zhang, A. Zimmerman, D. A. Nichols, E. Berti, and Y. Chen, Physical Review D 87 (2013)

  3. [11]

    H. Yang, A. Zimmerman, A. Zengino˘ glu, F. Zhang, E. Berti, and Y. Chen, Phys. Rev. D 88, 044047 (2013)

  4. [12]

    G. F. Torres del Castillo, 3-D Spinors, Spin-Weighted Functions and their Applications (Birkh¨ auser Boston, Boston, MA, 2003)

  5. [13]

    R. A. Breuer, M. P. Ryan, and S. Waller, Proc R Soc London Ser A 358, 71 (1977)

  6. [14]

    Y3 is spin weighted spherical harmonic with parameters (s3,ℓ 3,m 3)

    := π∫ 0 d¯θ(sin ¯θ)S1S2Y3(¯θ), (H20) whereS1,S 2 are well behaved functions of ¯θ e.g., smooth on [0,π ] up to and including the boundary of the interval. Y3 is spin weighted spherical harmonic with parameters (s3,ℓ 3,m 3). In our computation of the overlap coeffi- cients in the...

  7. [15]

    In the following formulas, ℓ3 corresponds to the low ℓ QNM

    (high),(low) → (high) Here we give the explicit form of the overlap coeffi- cients for one low ℓ and two high ℓ for vanishing overtone numbers, i.e., Ni = 0. In the following formulas, ℓ3 corresponds to the low ℓ QNM. Since U near 123 is symmetric in (23), the case when ℓ2 is th...

  8. [16]

    Klainerman and J

    S. Klainerman and J. Szeftel, Pure Appl. Math. Quart. 19, 791 (2023) , arXiv:2104.11857 [math.AP]

  9. [17]

    Dafermos, G

    M. Dafermos, G. Holzegel, I. Rodnianski, and M. Taylor, arXiv preprint arXiv:2104.08222 (2021)

  10. [18]

    Casals, S

    M. Casals, S. E. Gralla, and P. Zimmerman, Physical Review D 94, 064003 (2016)

  11. [19]

    (high),(high) → (high)

    for the definition of h+ ≡ sh+ℓm], which in view of the struc- ture ( L1) of H a 0,H a ± at once implies the following. For any given (s,N,ℓ,m ), sΥ near N ℓm must be an eigenfunction of LH0 , in fact with eigenvalue N +h. Due to the factors of e±¯t in H a ±, LH± map sΥ near N ...

  12. [20]

    Mitman, M

    K. Mitman, M. Lagos, L. C. Stein, S. Ma, L. Hui, Y. Chen, N. Deppe, F. H´ ebert, L. E. Kidder, J. Moxon, M. A. Scheel, S. A. Teukolsky, W. Throwe, and N. L. Vu, Physical Review Letters 130 (2023)

  13. [21]

    M. H.-Y. Cheung, V. Baibhav, E. Berti, V. Car- doso, G. Carullo, R. Cotesta, W. Del Pozzo, F. Duque, T. Helfer, E. Shukla, et al. , Physical Review Letters 130, 081401 (2023)

  14. [22]

    Ma and H

    S. Ma and H. Yang, Physical Review D 109 (2024)

  15. [23]

    H. Yang, A. Zimmerman, and L. Lehner, Physical review letters 114, 081101 (2015)

  16. [24]

    S. R. Green, S. Hollands, and P. Zimmerman, Classical and Quantum Gravity 37, 075001 (2020)

  17. [25]

    J. M. Bardeen, W. H. Press, and S. A. Teukolsky, Astrophys. J. 178, 347 (1972)

  18. [26]

    Khera, S

    N. Khera, S. Ma, and H. Yang, arXiv preprint arXiv:2410.14529 (2024)

  19. [27]

    A. J. Amsel, G. T. Horowitz, D. Marolf, and M. M. Roberts, Journal of High Energy Physics 2009, 044 (2009)

  20. [28]

    Bredberg, C

    I. Bredberg, C. Keeler, V. Lysov, and A. Strominger, Nuclear Physics B-Proceedings Supplements 216, 194 (2011)

  21. [29]

    Bizo´ n, M

    P. Bizo´ n, M. Maliborski, and A. Rost- worowski, Phys. Rev. Lett. 115, 081103 (2015) , arXiv:1506.03519 [gr-qc]

  22. [30]

    Moschidis, Invent

    G. Moschidis, Invent. Math. 231, 467 (2023) , arXiv:1812.04268 [math.AP]

  23. [31]

    Balasubramanian, A

    V. Balasubramanian, A. Buchel, S. R. Green, L. Lehner, and S. L. Liebling, Phys. Rev. Lett. 113, 071601 (2014) , arXiv:1403.6471 [hep-th]

  24. [32]

    [41], λJ = √ 2(r −ia cosθ)√ ∆

    and BL co- ordinates, we have e.g. [41], λJ = √ 2(r −ia cosθ)√ ∆ . (C3) In nNHEK in the frame ( 34) and the coordinates ( 28), we have λJ =M √ 2 (1 −i cos ¯θ√f ) (C4) instead, where f = ¯x(¯x + 2). For the case of I one may give an analogous construc- tion. In this case, we ha...

  25. [33]

    Flammer, Spheroidal wave functions (Courier Corpo- ration, 2014)

    C. Flammer, Spheroidal wave functions (Courier Corpo- ration, 2014)

  26. [34]

    and the coordinates ¯xµ = ( ¯t, ¯x, ¯θ, ¯φ) of Eq. ( 28). These spin coefficients and the non-vanishing optical scalars are, in GHP notations [55] τ = −i sin ¯θ M √ 2 ( 1 + cos2 ¯θ ), (A1a) τ ′ = − i sin ¯θ M √ 2(1 −i cos ¯θ)2, (A1b) β = cot ¯θ 2 √ 2M (1 +i cos ¯θ), (A1c) β′ = −...

  27. [35]

    H¨ afner, P

    D. H¨ afner, P. Hintz, and A. Vasy, Invent. Math. 223, 1227 (2020) , [Erra- tum: Invent.Math. 236, 477–481 (2024)], arXiv:1906.00860 [math.AP]

  28. [36]

    L. S. Kegeles and J. M. Cohen, Physical Review D 19, 1641 (1979)

  29. [37]

    P. L. Chrzanowski, Physical Review D 11, 2042 (1975)

  30. [38]

    S. E. Gralla, A. Zimmerman, and P. Zim- merman, Phys. Rev. D 94, 084017 (2016) , arXiv:1608.04739 [gr-qc]

  31. [39]

    Aretakis, Commun

    S. Aretakis, Commun. Math. Phys. 307, 17 (2011) , arXiv:1110.2007 [gr-qc]

  32. [40]

    Bantilan, P

    H. Bantilan, P. Figueras, M. Kunesch, and P. Ro- matschke, Phys. Rev. Lett. 119, 191103 (2017) , arXiv:1706.04199 [hep-th]

  33. [41]

    Angelopoulos, S

    Y. Angelopoulos, S. Aretakis, and D. Gajic, Advances in Mathematics 417, 108939 (2023)

  34. [42]

    J. M. Baarden, in Les Houches Summer School of Theo- retical Physics: Black Holes (1973) pp. 241–290

  35. [43]

    V. E. Zakharov, V. S. L’vov, and G. Falkovich, Kolmogorov spectra of turbulence I: Wave turbulence (Springer Science & Business Media, 2012)

  36. [44]

    S´ erazin, T

    G. S´ erazin, T. Penduff, B. Barnier, J.-M. Mo- lines, B. K. Arbic, M. M¨ uller, and L. Terray, Journal of Physical Oceanography 48, 1385 (2018)

  37. [45]

    Bardeen and G

    J. Bardeen and G. T. Horowitz, Physical Review D 60, 104030 (1999)

  38. [46]

    Baier, P

    R. Baier, P. Romatschke, D. T. Son, A. O. Starinets, and M. A. Stephanov, JHEP 04, 100 , arXiv:0712.2451 [hep-th]

  39. [47]

    Adams, P

    A. Adams, P. M. Chesler, and H. Liu, Phys. Rev. Lett. 112, 151602 (2014) , arXiv:1307.7267 [hep-th]

  40. [48]

    S. R. Green, F. Carrasco, and L. Lehner, Phys. Rev. X 4, 011001 (2014) , arXiv:1309.7940 [hep-th]

  41. [49]

    S. A. Teukolsky,Physical Review Letters 29, 1114 (1972)

  42. [50]

    Geroch, A

    R. Geroch, A. Held, and R. Penrose, J. Math. Phys. 14, 874 (1973)

  43. [51]

    Craps, O

    B. Craps, O. Evnin, and J. Vanhoof, JHEP 10 (10), 048, arXiv:1407.6273 [gr-qc]

  44. [52]

    These operators satisfy the same commutation relations and furthermore commute with the spins Teukolsky operators, [sO, LX ] = 0

    LX,X a ∈ {H a 0,H a ±}. These operators satisfy the same commutation relations and furthermore commute with the spins Teukolsky operators, [sO, LX ] = 0. Hence the complex linear span of the near zone QNM solutions [see Eq. ( 103a)] to the spin s Teukolsky equations is a modul...

  45. [53]

    envelope

    that m = 0 and that ωN ℓ0 = −iε(h+ +N )/(2M ) = −iε(ℓ + 1 +N )/(2M ) up to terms of order O(ε2), cf. Eqs. ( 118), ( 95). Neglecting O(ε3)-terms, the approximation to the angular Teukolsky equation replacing Eq. ( 101) becomes [ 1 sin ¯θ d d¯θ ( sin ¯θ d d¯θ ) + ( sEℓ0 − s2 sin...

  46. [54]

    mathematical aspects of general rela- tivity

    C. Kehle, Talk at “mathematical aspects of general rela- tivity” mfo meeting id 2432 (August 2024)

  47. [55]

    Srednicki, Quantum field theory (Cambridge Univer- sity Press, 2007)

    M. Srednicki, Quantum field theory (Cambridge Univer- sity Press, 2007)

  48. [56]

    Hollands and V

    S. Hollands and V. Toomani, arXiv preprint arXiv:2405.18604 (2024)

  49. [57]

    S. A. Teukolsky, Astrophys. J. 185, 635 (1973)

  50. [58]

    W. H. Press and S. A. Teukolsky, Apj 185, 649 (1973)

  51. [59]

    Casals, S

    M. Casals, S. Hollands, A. Pound, and V. Toomani, arXiv preprint arXiv:2402.15468 (2024)

  52. [60]

    S. R. Green, S. Hollands, L. Sberna, V. Toomani, and P. Zimmerman, Physical Re- view D 10.1103/PhysRevD.107.064030 (2022), arXiv:2210.15935

  53. [61]

    Ori, Phys

    A. Ori, Phys. Rev. D 67, 124010 (2003) , arXiv:0207045 [gr-qc]

  54. [62]

    E. W. Leaver, J. Math. Phys. 27, 1238 (1986)

  55. [63]

    Gautschi, SIAM Review 9, 24 (1967) , https://doi.org/10.1137/1009002

    W. Gautschi, SIAM Review 9, 24 (1967) , https://doi.org/10.1137/1009002

  56. [64]

    Bhattacharyya, V

    S. Bhattacharyya, V. E. Hubeny, S. Minwalla, and M. Rangamani, JHEP 02, 045 , arXiv:0712.2456 [hep-th]

  57. [65]

    L. R. Price, Developments in the perturbation theory of algebraically special spacetimes , Ph.D. thesis, U. Florida (2007)

  58. [66]

    Starobinskii and S

    A. Starobinskii and S. Churilov, Zh. eksp. teor. Fiz 65, 3 (1973)

  59. [67]

    S. A. Teukolsky and W. H. Press, Astrophysical Journal, vol. 193, Oct. 15, 1974, pt. 1, p. 443-461. 193, 443 (1974)

  60. [68]

    E. W. Leaver, Physical Review D 34, 384 (1986)

  61. [69]

    Casals, C

    M. Casals, C. Kavanagh, and A. C. Ottewill, Physical Review D 94, 124053 (2016)

  62. [70]

    R. M. Wald, General relativity (University of Chicago press, 2010)

  63. [71]

    S. B. Edgar and G. Ludwig, General Relativity and Gravitation 32, 637 (2000)

  64. [72]

    Geroch, Commun.Math

    R. Geroch, Commun.Math. Phys. 13, 180 (1969)

  65. [73]

    S. E. Gralla and P. Zimmerman, 51 JHEP 06 (6), 061, arXiv:1804.04753 [gr-qc]

  66. [74]

    Comp` ere, K

    G. Comp` ere, K. Fransen, T. Hertog, and J. Long, Class. Quantum Gravity 35, 104002 (2018)

  67. [75]

    Chandrasekhar, The mathematical theory of black holes, Vol

    S. Chandrasekhar, The mathematical theory of black holes, Vol. 69 (Oxford university press, 1998)

  68. [76]

    Araneda, Class

    B. Araneda, Class. and Quantum Gravity 34, 035002 (2016)

  69. [77]

    Prabhu and R

    K. Prabhu and R. M. Wald, Class. Quantum Gravity 35, 1 (2018) , arXiv:1807.09883

  70. [78]

    non-modal stability

    corre- sponding to the equilibrium Hertz potential φeq, and where xeq ab the corresponding corrector, see Eq. ( 94). The upshot is that, by Eqs. ( 170), or ( 171), the equi- librium metric has dyadically small amplitudes for large ℓ QNMs, and thus the equilibrium metric, which...

  71. [79]

    Iyer and R

    V. Iyer and R. M. Wald, Phys. Rev. D 50, 846 (1994)

  72. [80]

    Cannizzaro, L

    E. Cannizzaro, L. Sberna, S. R. Green, and S. Hollands, Physical Review Letters 132, 051401 (2024)

  73. [81]

    Kato, Perturbation theory for linear operators , Vol

    T. Kato, Perturbation theory for linear operators , Vol. 132 (Springer Science & Business Media, 2013)

  74. [82]

    Spiers, A

    A. Spiers, A. Pound, and J. Moxon, arXiv preprint arXiv:2305.19332 (2023)

  75. [83]

    B. F. Whiting, J. Math. Phys. 30, 1301 (1989)

  76. [84]

    Kavanagh, A

    C. Kavanagh, A. C. Ottewill, and B. Wardell, Phys. Rev. D 93, 124038 (2016)

  77. [85]

    A. G. Shah and B. F. Whiting, General Relativity and Gravitation 48, 10.1007/s10714-016-2064-z (2016)

  78. [86]

    Ferrari and B

    V. Ferrari and B. Mashhoon, Physical Review D 30, 295 (1984)

  79. [87]

    H. Yang, D. A. Nichols, F. Zhang, A. Zimmerman, Z. Zhang, and Y. Chen, Phys. Rev. D 86, 104006 (2012) , arXiv:1207.4253 [gr-qc]

  80. [88]

    Berti, V

    E. Berti, V. Cardoso, and M. Casals, Phys. Rev. D 73, 024013 (2006)

  81. [89]

    R. H. Price, Physical Review D 5, 2439 (1972)

  82. [90]

    Campanelli and C

    M. Campanelli and C. O. Lousto, Phys. Rev. D 59, 124022 (1999)

  83. [91]

    Detweiler, Astrophys

    S. Detweiler, Astrophys. J. 239, 292 (1980)

  84. [92]

    Sasaki and T

    M. Sasaki and T. Nakamura, General relativity and gravitation 22, 1351 (1990)

  85. [93]

    Joykutty, in Annales Henri Poincar´ e, Vol

    J. Joykutty, in Annales Henri Poincar´ e, Vol. 23 (Springer, 2022) pp. 4343–4390

  86. [94]

    Hintz and Y

    P. Hintz and Y. Xie, Physical Review D 104, 064037 (2021)

  87. [95]

    Hintz and Y

    P. Hintz and Y. Xie, Journal of Mathematical Physics 63 (2022)

  88. [96]

    Gajic and C

    D. Gajic and C. Warnick, Commun.Math. Phys. 385, 1395 (2021)

  89. [97]

    DLMF, NIST Digital Library of Mathematical Functions , https://dlmf.nist.gov/, Release 1.2.1 of 2024-06-15, f. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller, B. V. Saunders, H. S. Cohl, and M. A. McClain, eds

  90. [98]

    Berti, V

    E. Berti, V. Cardoso, and A. O. Starinets, Class. Quant. Grav. 26, 163001 (2009) , arXiv:0905.2975 [gr-qc]

  91. [99]

    K. D. Kokkotas and B. G. Schmidt, Living Rev. Rel. 2, 2 (1999) , arXiv:gr-qc/9909058

  92. [101]

    For this, we self-consistently assume in the full angular Teukolsky equation (

    of the solutions to the angular Teukolsy equation by including terms of O(ε). For this, we self-consistently assume in the full angular Teukolsky equation (

  93. [102]

    A. G.-P. G´ omez-Lobo,The N.P. and G.H.P. formalisms. (2009)

  94. [108]

    A. O. Barut and C. Fronsdal, Proc. R. Soc. Lond. A 287, 532 (1965)

  95. [109]

    In view of Eqs

    of sh+, we find that sh+ = ℓ + 1. In view of Eqs. ( D9), ( D10) we would thereby conclude for instance that the scalar product of an axisymmetric QNM would diverge for s = ±2,m = 0,N = 0. A similar conclusion can be drawn for all axisymmetric QNMs. In order to be self-consisten...

  96. [110]

    Stigner, On tensor product of non-unitary repre- sentations of sl(2,R) , Ph.D

    C. Stigner, On tensor product of non-unitary repre- sentations of sl(2,R) , Ph.D. thesis , Karlstadt University (2007)

  97. [111]

    P. J. Schmid, Annual Review of Fluid Mechanics 39, 129 (2007)

  98. [112]

    R. M. Wald, Physical Review Letters 41, 203 (1978)

  99. [113]

    Bredberg, T

    I. Bredberg, T. Hartman, W. Song, and A. Strominger, Journal of High Energy Physics 2010, 1 (2010)

  100. [114]

    J. N. Goldberg, A. J. Macfarlane, E. T. New- man, F. Rohrlich, and E. C. G. Sudarshan, J. Math. Phys. 8, 2155 (1967)

  101. [115]

    Bakaleinikov and A

    L. Bakaleinikov and A. Silbergleit, J. Math. Phys. 61, 10.1063/5.0016206 (2020)

  102. [116]

    Gervois and H

    A. Gervois and H. Navelet, J. Math. Phys. 25, 3350 (1984)

  103. [117]

    Rosenheinrich, Tables of Some Indefinite Integrals of Bessel Functions (2016)

    W. Rosenheinrich, Tables of Some Indefinite Integrals of Bessel Functions (2016)

  104. [123]

    Y2,Y 3 are spin weighted spherical harmonics, as before

    := π∫ 0 d¯θ(sin ¯θ)S1Y2Y3(¯θ), (H15) whereS1 is some well behaved function of ¯θ e.g., smooth on [0,π ] up to and including the boundary of the interval, which is independent of ℓ2,ℓ 3. Y2,Y 3 are spin weighted spherical harmonics, as before. In our computation of the overlap ...

  105. [125]

    (D9) Fors = +2, we can use ( 64) and get +2Anear 0ℓm =M 2 3 2− 1 2 −(h++im)× Γ(2 h+)Γ( −h+ −im + 3)Γ( h+ +im + 2) |Γ( h+ −im − 2)|2

    may be further simplified for overtone num- ber N = 0: −2Anear 0ℓm =M 10 3 2 7 2 −(h++im) Γ(2 h+)Γ( −h+ −im + 3) Γ( h+ −im + 2) . (D9) Fors = +2, we can use ( 64) and get +2Anear 0ℓm =M 2 3 2− 1 2 −(h++im)× Γ(2 h+)Γ( −h+ −im + 3)Γ( h+ +im + 2) |Γ( h+ −im − 2)|2 . (D10) These no...

  106. [156]

    and the definition ( D25), this trans- lates into 29: φeq(¯xµ) ∝ (low) +Clow,eq ∑ ε−q/greaterorsimilarℓ/greaterorsimilarε−1/2 2− ℓ 2ℓ− 7 2 · −2Υ near 0ℓ0 (¯t, ¯x, ¯θ), (170) in the near horizon zone where −2Υ near 0ℓ0 (¯xµ) are the near zone, zero-damped QNMs as in Eq. ( 150). ...

  107. [165]

    is still quite intractable due to the complicated forms of the overlap coefficients as indicated by ... in Eq. ( 165). With the goal of identifying an equilibrium distribution, we now make further simplifications to Eqs. ( 165). We restrict to: • axisymmetric QNMs (m = 0), • zero...

  108. [167]

    were (0, ∞) instead of (c,ε 1/2−q), the two terms on the right side of the sec- ond equation in ( 167) would cancel for a suitable propor- tionality factor of the order of Clow,eq, the total QNM amplitude of the low ℓ QNMs, and for p = −2. For these choices, the double integra...

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.