Pith. sign in

REVIEW 3 major objections 4 minor 114 references

General relativity has a conserved curvature measure—magnetic Weyl enstrophy—that forces nonlinear wave energy toward lower frequencies.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

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-extremal Kerr and AdS.

T0 review reviewed 2026-08-05 challenge →

load-bearing objection Candid, clean derivation of a magnetic Weyl balance law, but the advertised near-extremal inverse cascade is not established: Eq. (27) admits that the fastest channel breaks the conserved quantity, leaving only a one-sided constraint that does not force lower-frequency transfer. the 3 major comments →

arxiv 2608.03697 v1 pith:DMQ6FUJ5 submitted 2026-08-04 gr-qc astro-ph.HEhep-th

Gravitational Enstrophy: Local Geometric Origin and Inverse-Cascade Constraints

Luis Lehner This is my paper
classification gr-qc astro-ph.HEhep-th MSC 83C0583C3583C57
keywords gravitational enstrophymagnetic Weyl tensorFjørtoft constraintKerr black holesinverse cascadequasinormal modesAdS/CFT holographygravitational turbulence
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The reading

This paper sets out to establish that general relativity has a conserved quantity playing the role enstrophy plays in two-dimensional fluids: the integrated square of the magnetic Weyl curvature, Z = ∫ B_ab B^ab √γ d³x. For linear radiative perturbations of Kerr and Kerr–AdS, this quantity is approximately conserved in the zero-angular-momentum frame, and mode by mode it carries exactly the right spectral weight, δZ_k/W_k = ω_k². Because energy and this enstrophy are both conserved to quadratic order while nonlinear interactions redistribute mode energy, the two moments force a gravitational Fjørtoft constraint: energy cannot cascade purely toward high frequencies; the upward part must be balanced by a larger transfer downward. The paper shows the constraint is dynamically relevant where damping is slow—near-extremal Kerr and reflecting AdS—and suppressed in generic ringdown, and it maps the bulk quantity to boundary fluid enstrophy in the fluid/gravity setting.

Core claim

Central claim: the magnetic part of the Weyl tensor—the frame-dragging sector of spacetime curvature—is the gravitational analogue of fluid vorticity, and its integrated square is an approximately conserved second moment for radiative perturbations. In the radiation zone of type D vacuums, δZ_k/W_k = ω_k² mode by mode, with near-horizon corrections O(√(1−χ)) for zero-damping modes. In the ZAMO frame the Bianchi evolution has vanishing vorticity coupling and cancelling curl exchange; remaining terms are controlled, so dδZ/dt = O(ε³). With energy conservation this yields a gravitational Fjørtoft theorem: more energy flows to lower frequencies than to higher. Active regimes are near-extremal Ke

What carries the argument

The supporting object is B_ab = ⋆C_acbd u^c u^d, the magnetic (frame-dragging) part of the Weyl tensor as seen by an observer u^a; its integrated square Z is the enstrophy. Three mechanisms carry the argument: the two-derivative relation between metric amplitude and radiative curvature, which fixes η_k = δZ_k/(ω_k²W_k) = 1; the ZAMO-frame cancellations—symmetric-times-antisymmetric vorticity contraction vanishes pointwise, and δE²=δB² makes curl exchange a harmless transfer between sectors; and the Fjørtoft algebra of two conserved positive moments with weights 1 and ω², which in Theorem 1 forces lower-frequency-dominant redistribution. Conditions (C1)–(C4) delimit when the balance law is co

Load-bearing premise

The load-bearing premise is that near the horizon, in the zero-angular-momentum frame, the radiative curvature falls off faster than the redshift factor (|δB|² ∼ α^(2p) with p > 0); the paper states this exponent is not fixed by background geometry and must be verified numerically.

What would settle it

Measure p from a second-order Teukolsky or fully nonlinear evolution at χ ≈ 0.99: if p ≤ 0, the integrated acceleration source diverges and Z is not conserved. Independently, in a global AdS pure-gravity run with multi-mode data, track ⟨ω²⟩ = Z/W: a drift on the nonlinear timescale shows the resonant interactions do not conserve Z, so the constraint does not bind there.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • In near-extremal Kerr, the constraint predicts nonlinear transfer piles energy near ω ≈ mΩ_H, an observable spectral excess over linear quasinormal-mode predictions.
  • In reflecting AdS, energy and enstrophy are both conserved to O(ε³) for transverse radiative modes, giving a closed vacuum laboratory for the inverse cascade.
  • Generic Kerr ringdown remains cascade-suppressed because τ_damp ≲ τ_nl; there, Z is a diagnostic of the spectral direction, not a dynamical driver.
  • In the fluid/gravity correspondence, the bulk magnetic Weyl enstrophy maps to boundary fluid enstrophy with coefficient (2624π⁵/189)T⁵, so the ratio R_holo = 1 diagnoses the hydrodynamic regime in holographic turbulence.
  • The (W, Z) pair is one of several conserved two-moment pairs; the scalar AdS instability conserves (E, N) instead, so its direct cascade does not contradict the Fjørtoft logic.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Inference: if the near-horizon fall-off exponent is confirmed positive, the enstrophy balance becomes a hard conservation law on the cascade timescale, and high-signal near-extremal ringdowns should show excess power near ω ≈ mΩ_H—a clean separation from linear quasinormal-mode templates.
  • Inference: the same two-moment algebra may work for other positive quadratic Weyl functionals, such as the electric enstrophy or the Bel–Robinson super-energy; the paper chooses B² as the cleanest representative, but the proof structure suggests a family of conserved moments.
  • Inference: whether a third conserved quantity exists in the resonant mode-coupling of near-extremal Kerr would decide between recurrent cascades and monotone attraction to ω = mΩ_H; this can be settled by computing the sign and magnitude of the second-order Teukolsky coupling coefficients.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper proposes the magnetic Weyl functional Z = ∫ B_ab B^ab √γ d^3x as a gravitational analogue of fluid enstrophy. For linear radiative perturbations of Petrov type D backgrounds (Kerr, Kerr–AdS), it argues that the modal spectral ratio δZ_k/W_k equals ω_k^2 (Proposition 1); that the local Bianchi balance law for δZ reduces, in the ZAMO frame, to an approximate conservation law dδZ/dt = O(ε^3) under conditions (C1)–(C4); and that the pair (W, Z) then yields a gravitational Fjørtoft theorem forcing preferential transfer toward lower frequencies. The paper applies this to near-extremal Kerr, Kerr–AdS and horizonless AdS, and gives a fluid/gravity map between the bulk magnetic Weyl norm and boundary fluid enstrophy.

Significance. If the approximate conservation law held as stated, this would be an important addition to gravitational perturbation theory: it identifies a curvature-level second moment that organizes spectral transfer in a way analogous to 2D turbulence, and it gives explicit falsifiable targets (spectral excess near ω = mΩ_H; R_holo = 1 in holographic turbulence). The paper's algebraic machinery is largely transparent and is a genuine strength: the pointwise cancellation of the ZAMO vorticity coupling (Sec. IV A), the self-adjoint curl integration by parts (Appendix A 3), the shear decoupling statement (Prop. 2), the explicit numerical checks (Tables III and V), and the closed-form fluid/gravity coefficient Eq. (52) are useful concrete results. The kinematic Fjørtoft theorem itself (Theorem 1) is standard and clean. However, the dynamical premise needed to apply it in the advertised regimes is not established: the fastest near-extremal nonlinearity violates Z conservation, and the near-horizon acceleration convergence rests on an unverified falloff exponent. The headline physical claim is therefore stronger than the paper's own equations support.

major comments (3)
  1. [Sec. V C and Sec. VI C] The equality-form Fjørtoft constraint is not available in the advertised near-extremal regime. Eq. (27) shows that each resonant three-wave interaction changes Z by 3ω1ω2ω3 δn3, i.e. an O(1) change relative to the modal enstrophy, while the paper's own nonlinearity analysis (Eq. (28), Sec. VI C, Fig. 2) identifies the three-wave channel as the fastest near-extremal transfer, τ_nl^{-1} ∼ v ε ω_R, and Table V confirms the triads are kinematically open. Thus ΔZ/Z ∼ O(1) over one nonlinear time, not the O(ε^3) remainder asserted in Eq. (14); premise (20) fails precisely where the cascade is claimed active. What survives is the one-sided Corollary 2, which rules out a purely direct cascade but does not force preferential lower-frequency transfer; the direction is set by coupling coefficients that the paper defers. This limitation must be moved into the abstract and main theorem, or the three-
  2. [Appendix E 7] Condition (C3) depends on convergence of the ZAMO acceleration source. Equation (E38) gives a_hatr |δB|^2 sqrt(g_rr) ∼ α^{2p−2}, so the horizon integral converges only if p > 0, where p is the near-horizon falloff exponent |δB|^2 ∼ α^{2p}. The paper explicitly states that p is not fixed by the background geometry and defers its verification to future numerics (Sec. VIII B). Since the same appendix notes that a perturbation regular in a horizon-penetrating frame can have p ≤ 0 in the singular ZAMO frame, the conservation law is conditional on an unverified input. The abstract and Sec. VI C present the near-extremal constraint as operative; they should either carry the explicit p > 0 condition as a stated assumption or provide numerical evidence for physical ZDMs.
  3. [Sec. III, Prop. 1, Remark 3(i)] The central ratio η_k = 1 is definitional. With W_k defined in Eq. (6) as (1/(2ω_k^2)) ∫ |Ψ4|^2 and δZ_k defined as (1/2)∫ |Ψ4|^2, the cancellation δZ_k/(ω_k^2 W_k) = 1 is built into the definitions; no dynamics enters. Remark 3(i) concedes this, and the genuine content is the gauge/frame reduction δB^2 = 1/2 |Ψ4|^2 and the mode-independence of the ratio. The abstract and introduction, however, present η_k = 1 as a demonstrated physical result ('we show that η_k = 1'), which overstates the status. The paper should be rephrased so that the spectral weighting appears as a normalization/definitional identity plus a frame-reduction statement, not as an independent dynamical prediction.
minor comments (4)
  1. [Abstract] Typo: 'constraintt implyging' should read 'constraint implying'.
  2. [Fig. 1] The caption notation 'n = 3 R220, −2S22' is unexplained; please define the spheroidal harmonic and overtone labels.
  3. [Sec. VIII C] The observational signatures are speculative and depend on Conjecture 1, which the paper itself identifies as unproven. Consider labeling this subsection explicitly as conditional on the conjecture.
  4. [Sec. V C] The sentence 'This change is the secular content of the O(ε^3) remainder' is potentially confusing: Eq. (27) shows the change is not O(ε^3) relative to Z but O(1) per transfer event. Clarify the distinction between field-order counting and relative change of the conserved moment.

Circularity Check

1 steps flagged

The exact spectral ratio η_k = 1 is a normalization identity: Eq. (6) defines W_k and δZ_k from the same |Ψ_4|^2 integral, so the ω^2 weighting that drives the Fjørtoft theorem is put in by construction. The independent content is the approximate conservation of δZ, which is conditional and not circular.

specific steps
  1. self definitional [Sec. III B, Eq. (6)–(7), Proposition 1; Remark 3(i)]
    "W_k[D] = 1/(2ω_k^2) ∫_D |Ψ_4^(1)|^2 √γ d^3x, δZ_k[D] = ∫_D δB_ab^(k)δB^ab_(k) √γ d^3x ... η_k ≡ δZ_k[D]/(ω_k^2 W_k[D]) = 1. (Remark 3(i)): 'With the Isaacson normalization W_k ≡ ω_k^{-2} ∫ δB_ab δB^ab √γ d^3x, the radiation-zone equality η_k = 1 is an identity between two functionals built from the same curvature density.'"

    Proposition 1's exact equality is not derived from independent definitions of energy and enstrophy: W_k and δZ_k are both defined from the same ∫|Ψ_4|^2 integral, with W_k carrying a 1/ω_k^2 factor chosen so that the ratio equals 1 identically. The paper's own Remark 3(i) concedes this. The subsequent spectral weighting Z ≈ Σ_k ω_k^2 W_k in Eq. (17), and hence the two-moment algebra of Theorem 1, inherits this constructed ratio. What is not circular is the separate dynamical claim that δZ is approximately conserved on τ_nl (Eq. (14) under C1–C4); that claim rests on the Bianchi-identity balance and on unverified conditions such as the near-horizon falloff exponent p>0, not on the normalization of W_k. So the circularity is partial: the 'prediction' of ω_k^2 spectral weighting is a normaliz

full rationale

The central circular step is confined to the spectral-weighting result. Proposition 1 defines the Isaacson energy in Eq. (6) as (1/2ω_k^2)∫|Ψ_4|^2√γ d^3x and the magnetic Weyl enstrophy as ∫δB^2√γ d^3x = (1/2)∫|Ψ_4|^2√γ d^3x in the transverse radiative sector; substituting these definitions makes η_k = δZ_k/(ω_k^2 W_k) = 1 an algebraic identity. The paper explicitly acknowledges this in Remark 3(i), calling the equality 'an identity between two functionals built from the same curvature density.' The gravitational Fjørtoft theorem then uses only this constructed ratio, together with approximate conservation of W and Z, so the advertised lower-frequency bias is, at the level of the weighting, a definitional input rather than a derived prediction. I do not score the approximate-conservation claim as circular: it is derived from the Bianchi identities, with vorticity cancellation, curl-exchange cancellation, and explicit smallness conditions (C1)–(C4), and it is openly conditional on the near-horizon falloff exponent p>0 (Appendix E7). That condition is unverified but not circular. Likewise, the self-citations to [3], [5], and [42] for ZDM lifetimes, nonlinear transfer rates, and numerical evidence are external, code-reproduced or independently checkable results, so they do not constitute load-bearing circular self-citation under the review rules. The skeptical objection that resonant three-wave interactions change Z at order unity (Eq. (27)) is a dynamical-failure argument against the equality-form constraint in near-extremal Kerr, not a circularity argument; the paper itself concedes that only the one-sided Corollary 2 survives in that regime. That reduces the scope of the claim but does not change the circularity score. Overall, because the paper's flagship spectral weighting and the consequent Fjørtoft biasing reduce to a normalization choice, while the genuinely independent conservation mechanism remains conditional and noncircular, the appropriate score is 6: partial circularity by construction.

Axiom & Free-Parameter Ledger

2 free parameters · 5 axioms · 1 invented entities

The central claim rests on the existence of an approximately conserved second moment; the paper introduces Z as such, but its spectral weighting is definitional while its conservation depends on several domain assumptions and the unverified near-horizon falloff.

free parameters (2)
  • Three-wave coupling coefficient v = v ~ 1.5
    Used in Fig. 2 for the resonant three-wave rate tau_nl^{-1} = v epsilon omega_R; presented as an estimate, not derived from second-order Teukolsky source.
  • Near-horizon radiative curvature falloff exponent p = unspecified (assumed > 0)
    Controls integrability of the ZAMO acceleration source (Appendix E7); the paper states it is not fixed by background geometry and leaves numerical verification to future work.
axioms (5)
  • standard math Goldberg-Sachs theorem: on Petrov type D vacuum backgrounds, the transverse radiative perturbation is carried by Psi_4 (and Psi_0) in Kinnersley gauge.
    Used to justify delta B^2 = 1/2 |Psi_4|^2 (Proposition 1, Sec. III).
  • domain assumption ZAMO congruence in stationary axisymmetric spacetimes is hypersurface-orthogonal, so vorticity and expansion vanish (Frobenius theorem).
    Sec. IV and Appendix E; basis for the algebraic cancellations in the balance law.
  • domain assumption Mode-by-mode equality delta E^2 = delta B^2 in the transverse radiative sector, with cross-mode terms vanishing by time averaging over beat times.
    Sec. IV B, Remark 4; needed for curl-exchange cancellation and for delta Z = sum omega_k^2 W_k.
  • domain assumption Nonlinear transfer time tau_nl ~ 1/(epsilon^2 omega_R) and coupling coefficients of order unity.
    Eq. (28); used to compare damping and transfer. The O(1) coefficient is estimated, especially the v ~ 1.5 three-wave coupling.
  • ad hoc to paper Near-horizon ZAMO-frame radiative curvature falls as |delta B|^2 ~ alpha^{2p} with p > 0.
    Appendix E7; required for convergence of the acceleration source; explicitly flagged as unverified.
invented entities (1)
  • Gravitational enstrophy Z (magnetic Weyl functional) independent evidence
    purpose: Quadratic curvature measure that acts as the second conserved spectral moment (alongside gravitational-wave energy) and drives the Fjørtoft constraint.
    Predicts spectral accumulation near omega = m Omega_H for near-extremal ringdown and R_holo ~ 1 in holographic turbulence; these are falsifiable outside the paper, though the omega^2 weighting itself is set by definition.

reviewed 2026-08-05 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Gravitational Enstrophy: Local Geometric Origin and Inverse-Cascade Constraints." pith.science (2026). https://pith.science/paper/DMQ6FUJ5

@misc{pith2026260803697,
  author       = {Pith},
  title        = {Pith review of: Gravitational Enstrophy: Local Geometric Origin and Inverse-Cascade Constraints},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/DMQ6FUJ5}},
  note         = {Machine review of arXiv:2608.03697}
}
Share X Bluesky LinkedIn Reddit HN
abstract

Two-dimensional fluids conserve energy and enstrophy, driving inverse energy cascades via Fj\o rtoft's argument. We show General Relativity admits an analogous structure: for linear radiative perturbations of Petrov type D backgrounds (Kerr, Kerr--AdS), the gravitational-wave energy $W = \sum_k W_k$ and magnetic Weyl enstrophy $\mathcal{Z} = \int B_{ab} B^{ab} \sqrt{\gamma} \, d^3x \approx \sum_k \omega_k^2 W_k$ are approximately conserved in the zero-angular momentum frame, where vorticity coupling vanishes identically and curl exchange cancels mode-by-mode. This yields a gravitational Fj\o rtoft constraintt implyging nonlinear energy transfer proceeds preferentially toward lower frequencies. The constraint is dynamically active in near-extremal Kerr ($\tau_{\text{damp}} \gg \tau_{\text{nl}}$) and confined geometries (AdS), but suppressed in generic ringdown. In AdS, $\mathcal{Z}$ maps holographically to the boundary fluid enstrophy.

Figures

Figures reproduced from arXiv: 2608.03697 by Luis Lehner.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Timescale separation for the gravitational Fjørtoft [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p035_3.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

114 extracted references · 38 canonical work pages · 7 internal anchors

  1. [1]

    The symmetric tensor field, γab =g ab +u aub ,(A4) acts as the spatial projector orthogonal tou a

    Spatial Projection and Kinematical Decomposition Letu a be a future-directed, unit timelike vector field (uaua =−1). The symmetric tensor field, γab =g ab +u aub ,(A4) acts as the spatial projector orthogonal tou a. Angular brackets denote the projected, symmetric, and trace-free (STF) part of a spatial tensor: T⟨ab⟩ = γ(a cγb) d − 1 3 γabγcd Tcd .(A5) Th...

  2. [2]

    Quadratic W eyl Invariants Using the orthogonality of the STF basis and the norms in Eq. (B21), the positive-definite Bel–Robinson energy density associated with the observeru a is: EabEab +B abBab =Q abQ ab =|Ψ 0|2 +|Ψ 4|2 + 4|Ψ1|2 + 4|Ψ3|2 + 6|Ψ2|2 .(B23) This positive-definite norm receiving contributions from all five Weyl scalars represents the total...

  3. [3]

    Applying the curl definition in Eq

    Curl Integration by Parts LetS ab andT ab be spatial, symmetric, trace-free ten- sors. Applying the curl definition in Eq. (A11) and not- ing that the contraction of any spatial symmetric tensor with the antisymmetric indices of the Levi-Civita tensor allows the projection brackets to be omitted, we have: Sab(curlT) ab =S abϵcdaDcTb d .(A14) Integrating t...

  4. [4]

    (A9) withB ab yields, withτthe proper time along ua: 1 2 d dτ BabBab =−B ab(curlE) ab −θB abBab + 3Babσc⟨aBb⟩ c +B ab Bc⟨aωb⟩ c −2B abϵcd⟨aacEb⟩ d

    Quadratic Balance and the Exact V orticity Identity Contracting the magnetic Weyl evolution equation Eq. (A9) withB ab yields, withτthe proper time along ua: 1 2 d dτ BabBab =−B ab(curlE) ab −θB abBab + 3Babσc⟨aBb⟩ c +B ab Bc⟨aωb⟩ c −2B abϵcd⟨aacEb⟩ d . (A17) As shown in IV A, the fourth term in the right hand side vanishes. This local cancellation requir...

  5. [5]

    Complex STF W eyl T ensor We combine the spatial electric and magnetic Weyl tensors into a single complex spatial STF tensor: Qab =E ab +iB ab .(B10) In the orthonormal spatial basis ( ˆθa,ˆφa,ˆra), its compo- nents are: Qθθ =−Ψ 2 + 1 2 (Ψ0 + Ψ4),(B11) Qφφ =−Ψ 2 − 1 2 (Ψ0 + Ψ4),(B12) Qrr = 2Ψ2 ,(B13) Qθφ = i 2 (Ψ0 −Ψ 4),(B14) Qθr = Ψ1 −Ψ 3 ,(B15) Qφr =i(Ψ...

  6. [6]

    Pure Radiative Sectors Let parenthesized numerical subscripts denote the re- striction to a subsector where only the indicated Weyl scalar is non-zero; thusQ ab(s),E ab(s), andB ab(s) rep- resent the complex Weyl tensor and its real/imaginary parts evaluated with Ψ s alone. For a pure outgoing transverse radiative mode (Ψ 4 ̸= 0), the complex STF tensor r...

  7. [7]

    This sector is non-radiative and does not participate in the spectral-transfer dynamics of the Fjørtoft cascade

    Coulomb Sector For a purely Coulombic field, only the Ψ 2 component is non-zero: Qab(2) = 2Ψ2Pab .(B35) Using the basis normP abP ab = 3/2, the quadratic in- variants evaluate to: Eab(2)Eab (2) +B ab(2)Bab (2) = 6|Ψ2|2 ,(B36) and Eab(2)Eab (2) −B ab(2)Bab (2) = 6 Re Ψ2 2 ,(B37) Eab(2)Bab (2) = 3 Im Ψ2 2 . This sector is non-radiative and does not particip...

  8. [8]

    III and shown to be satisfied un- der controlled assumptions in Sec

    F requency W eighting The enstrophy spectral weighting of the quadratic Weyl structure—namely, that radiative modes scale as Ψ(k) rad ∼ω 2 khk, implying R δB 2√γ d3x/Wk ∼ω 2 k—is analyzed in Sec. III and shown to be satisfied un- der controlled assumptions in Sec. H 3. We note that the helicity-like pseudoscalar functional constructed with E(k) ab Bab (k)...

  9. [9]

    B, the complex STF Weyl tensorQ ab =E ab +iB ab decomposes as in Eq

    The radiative/Coulomb split In the NP tetrad of Sec. B, the complex STF Weyl tensorQ ab =E ab +iB ab decomposes as in Eq. (B22): •Ψ 0,Ψ 4: transverse radiative data •Ψ 2: Coulombic data •Ψ 1,Ψ 3: longitudinal/mixed projections Schematically, δQab =δQ rad ab +δQ C ab +δQ L ab +δQ gauge ab ,(C1) 27 with δQrad ab ↔δΨ 0, δΨ4,(C2) δQC ab ↔δΨ 2,(C3) δQL ab ↔δΨ ...

  10. [10]

    With Ψ2 = ΨR 2 +iΨ I 2, EC ab = 2ΨR 2 Pab, B C ab = 2ΨI 2Pab.(C7) This contributes to the magnetic norm, but not as a prop- agating mode

    The Coulomb sector A purely Coulombic field has only Ψ 2:Q C ab = 2Ψ2Pab. With Ψ2 = ΨR 2 +iΨ I 2, EC ab = 2ΨR 2 Pab, B C ab = 2ΨI 2Pab.(C7) This contributes to the magnetic norm, but not as a prop- agating mode. There is no associated frequencyω k, so the spectral scalingZ k ∼ω 2 kWk does not apply. Example: in Schwarzschild with static observers, Ψ 2 is ...

  11. [11]

    OnlyδZ rad is the positive-definite enstrophy that pairs with energy

    Contamination of the magnetic norm Decompose the full perturbation as δBab =δB rad ab +δB C ab +δB L ab +δB gauge ab .(C10) The total magnetic norm expands as Z Σ δBabδB ab√¯γ d3x=δZ rad +δZ C +δZ L + 2Xrad,C + 2Xrad,L + 2XC,L +· · ·,(C11) with cross-terms likeX rad,C = R Σ δB rad ab δB ab C √¯γ d3x. OnlyδZ rad is the positive-definite enstrophy that pair...

  12. [12]

    Time averages vs. instantaneous norms A real radiative perturbation has both +ωand−ω components: δB rad ab =b abe−iωt + ¯babeiωt.(C13) The quadratic integrand is δB rad ab δB ab rad = 2bab¯bab +b abbabe−2iωt + ¯bab¯babe2iωt. (C14) The physical modal enstrophy is the stationary part, 2 R Σ bab¯bab√¯γ d3x—the time average over several periods. The oscillato...

  13. [13]

    The radiative subspace does not close exactly under the nonlinear Bianchi equations

    Coulomb terms in the balance law The radiative enstrophy balance takes the form d dt δZrad =R rad +R mix +F Z ∂Σ +O(ϵ 3),(C15) 28 whereR rad is built from radiative fields alone, whileRmix couples to Coulombic, longitudinal, and constraint sec- tors. The radiative subspace does not close exactly under the nonlinear Bianchi equations. Projecting onto Ψ 0/Ψ...

  14. [14]

    In Kerr and Kerr–AdS, the background curvature is Coulombic ( ¯Ψ2), while radiative perturbations live in δΨ0, δΨ4

    Implications for the main text regimes The role of the Coulomb sector varies by setting. In Kerr and Kerr–AdS, the background curvature is Coulombic ( ¯Ψ2), while radiative perturbations live in δΨ0, δΨ4. Dynamical mass and spin changes enter throughδΨ 2, but the relevant enstrophy is the projected radiative norm, not the raw variation of the stationary c...

  15. [15]

    The Super-Poynting V ector Applying the integration-by-parts identity for the spa- tial curl operator (derived in Appendix A 3) yields the boundary term: Z D Bab(curlE) ab √γ d3x= Z D Eab(curlB) ab √γ d3x + Z ∂D scϵcdaBabEb d√q d2x . (D2) To clarify the physical significance of the boundary in- tegrand, we define the spatial Bel–Robinson (or super- Poynti...

  16. [16]

    Magnetic Enstrophy Balance Recall the magnetic Weyl evolution equation has the schematic form: ˙B⟨ab⟩ + (curlE) ab =K B ab ,(D6) whereK B ab represents the local kinematical terms involv- ing spatial expansion, shear, vorticity, and acceleration. 29 Contracting this evolution equation with 2B ab and inte- grating over the domainD t yields: d dt ZB[Dt] =−2...

  17. [17]

    For an outgoing transverse mode, the duality rotation relates E and B

    Radiative Flux Through a Large Sphere Take a large sphere at radius r with outgoing normal sa. For an outgoing transverse mode, the duality rotation relates E and B. The radial super-Poynting flux therefore satisfies: |saP a|=E rad ab Eab rad =B rad ab Bab rad ,(D15) which, in terms of the outgoing Newman–Penrose scalar, scales as: Erad ab Eab rad =B rad ...

  18. [18]

    In a horizon-regular tetrad, the curvature flux through the horizon is governed by the normal component of the super-Poynting vector

    Horizon Flux For a black-hole exterior, the spatial domain is bounded internally by the future event horizon. In a horizon-regular tetrad, the curvature flux through the horizon is governed by the normal component of the super-Poynting vector. Schematically, this is written as: F Z H+ ∼ Z H+ |Ψhor|2 dv dA ,(D23) where Ψhor represents the Weyl scalar carry...

  19. [19]

    The local flux density remains: saP a =s aϵabcEb dBcd ,(D26) but the physical behavior is determined by the chosen boundary conditions at conformal infinity

    AdS Boundary Conditions In asymptotically Anti–de Sitter spacetimes, the con- formal boundary is timelike rather than null. The local flux density remains: saP a =s aϵabcEb dBcd ,(D26) but the physical behavior is determined by the chosen boundary conditions at conformal infinity. For reflecting AdS boundary conditions, the gravi- tational symplectic flux...

  20. [20]

    Let the boundary∂Dbe a timelike wall

    Finite Reflecting Domains The same logic applies to finite spatial domains bounded by artificial reflecting boundaries. Let the boundary∂Dbe a timelike wall. A perfect curvature- reflecting boundary condition imposes: saδP a = 0 on∂D,(D31) which yields: F Z ∂D = 0,(D32) meaning the enstrophy balance is governed entirely by the remaining bulk terms. If ins...

  21. [21]

    Numerical Evaluation The enstrophy flux can be evaluated directly fromE ij andB ij on each spatial slice. Lettings i be the outward unit normal to the extraction surface, the spatial super- Poynting vector is: δP i =ϵ ijk δEj ℓδBkℓ ,(D34) and the integrated flux is: F Z ∂D = 2 Z ∂D siδP i√q d2x ,(D35) where positive flux decreases the interior enstrophy. ...

  22. [22]

    Stationary Axisymmetric ADM F orm Consider a stationary, axisymmetric metric written in the standard ADM 3 + 1 form: ds2 =−α 2dt2 +γ ij dxi +β idt dxj +β jdt ,(E1) where indicesi, j∈ {r, θ, ϕ}. For a circular, stationary, and axisymmetric spacetime, we can choose coordinate charts such that the spatial metric is diagonal and the shift vector is purely azi...

  23. [23]

    For a stationary background spacetime, the metric sat- isfies∂ tγij = 0, reducing the extrinsic curvature to: Kij =− 1 2α (Diβj +D jβi).(E9) 32 Using the circular ADM form in Eq

    Extrinsic Curvature, Expansion, and Shear We define the extrinsic curvature of the spatial slices as: Kij =γ i aγj b∇anb .(E7) In terms of the ADM variables, this is expressed as: Kij = 1 2α (∂tγij −D iβj −D jβi).(E8) This definition differs by an overall sign from conven- tions often adopted in numerical relativity, but the rel- ative signs of the result...

  24. [24]

    For regular calculations at the horizon, a horizon- penetrating coordinate system must be employed

    Kerr Example For a Kerr black hole in Boyer–Lindquist coordinates, we define the standard functions: ρ2 =r 2 +a 2 cos2 θ ,∆ =r 2 −2M r+a 2 ,(E19) and A= (r 2 +a 2)2 −a 2∆ sin2 θ .(E20) The ADM metric functions are given by: α= ∆ρ2 A 1/2 , γ rr = ρ2 ∆ , γ θθ =ρ 2, γ ϕϕ = Asin 2 θ ρ2 , (E21) and the frame-dragging profile is: Ωfd = 2aM r A .(E22) Evaluating...

  25. [25]

    Kerr–AdS The same kinematic analysis applies to Kerr–AdS spacetimes. In standard Boyer–Lindquist-type coordi- nates, the metric is: ds2 =− ∆r ρ2 dt− asin 2 θ Ξ dϕ 2 + ρ2 ∆r dr2 + ρ2 ∆θ dθ2 + ∆θ sin2 θ ρ2 a dt−r2 +a 2 Ξ dϕ 2 ,(E26) where ρ2 =r 2 +a 2 cos2 θ ,∆ θ = 1− a2 ℓ2 cos2 θ ,Ξ = 1− a2 ℓ2 , (E27) and ∆r = (r2 +a 2) 1 + r2 ℓ2 −2M r .(E28) 33 From this ...

  26. [26]

    Shear Source in theB 2 Balance The shear term appearing in the local magnetic Weyl balance is: S B σ = 6 Z Dt BijσkiBj k√γ d3x ,(E30) where the coefficient of 6 arises from contracting the sym- metric, trace-free combination 3σ c⟨iBj⟩ c with 2B ij. In the ZAMO orthonormal frame, the only non- vanishing shear components areσ ˆrˆϕ andσ ˆθ ˆϕ and the tensor ...

  27. [27]

    Acceleration Source The acceleration contribution to the magnetic Weyl balance is also non-zero for ZAMO observers. From Eq. (E6), the orthonormal components of the accelera- tion are: aˆr= 1√γrr ∂r lnα , a ˆθ = 1√γθθ ∂θ lnα , a ˆϕ = 0. (E35) The acceleration term in theB 2 balance is given by: S B a =−4 Z Dt BijϵkliakEj l√γ d3x ,(E36) which dynamically c...

  28. [28]

    (E36) is the sole bulk term that cannot be bounded purely by algebraic or kine- matic arguments, because the ZAMO acceleration di- verges at the horizon

    Convergence of the Acceleration Source The acceleration source in Eq. (E36) is the sole bulk term that cannot be bounded purely by algebraic or kine- matic arguments, because the ZAMO acceleration di- verges at the horizon. Specifically, the Kerr ZAMO kine- matics yieldα aˆr→κ + asr→r + (whereκ + is the black hole surface gravity), meaning that the local ...

  29. [29]

    A Quantitative Smallness Criterion To apply the Fjørtoft argument, the shear and accel- eration source terms do not need to vanish identically, but they must remain small over the characteristic non- linear transfer timescaleτ nl. We define the dimensionless diagnostics: χσ(t)≡ δS B σ (t) δZB(t)/τnl , χ a(t)≡ δS B a (t) δZB(t)/τnl .(E41) The local bulk dy...

  30. [30]

    Interpretation for Kerr and Near-Extremal Kerr Our ZAMO analysis clarifies the physical role of black hole rotation. The presence of background rotation does not destroy the algebraic structure of the local Weyl en- strophy; the Bianchi identities, the electric-magnetic curl exchange, and the algebraic cancellation of vorticity re- main fully intact. Rota...

  31. [31]

    The shear and expansion are computed directly asσ ij =K ij − 1 3 Kγ ij andθ=K(whereK= 0 for stationary circular back- grounds, but may be non-zero in dynamical spacetimes)

    Numerical Implementation In a standard 3 + 1 numerical simulation, the ZAMO (or Eulerian) observer corresponds to the unit normal vectorn a, and all required variables are readily available on each spatial slice: the lapseα, the shiftβ i, the spatial metricγ ij, the extrinsic curvatureK ij, and the electric and magnetic Weyl tensorsE ij andB ij. The shear...

  32. [32]

    Continuous spectra For a continuous distribution, letW(Ω)≥0 be the spectral density. Then W= Z ∞ 0 W(Ω)dΩ,Z= Z ∞ 0 ΩW(Ω)dΩ.(F4) A redistributionδW(Ω) preserving both invariants satis- fies Z ∞ 0 δW(Ω)dΩ = 0, Z ∞ 0 ΩδW(Ω)dΩ = 0.(F5) The continuous form suits scattering states; the discrete form suits confined geometries. The structure is identi- cal

  33. [33]

    In 36 a degenerate subspace, the energy and enstrophy are positive-definite quadratic forms on the mode ampli- tudes

    Degeneracies and normalization Modes can share the same Ω due to polarizations, an- gular degeneracies, or standing-wave combinations. In 36 a degenerate subspace, the energy and enstrophy are positive-definite quadratic forms on the mode ampli- tudes. Diagonalize them simultaneously. Since the en- strophy weight is the same for all modes in the subspace,...

  34. [34]

    Schematically,H= P k χkΩkWk, where χk is a polarization or phase factor that can change sign

    The helicity-like invariant The mixed functional H= Z Σ EabBab√γ d3x(F7) has the same frequency scaling asE 2 andB 2, but it’s sign-indefinite. Schematically,H= P k χkΩkWk, where χk is a polarization or phase factor that can change sign. It cannot serve as the positive-definite enstrophy for the Fjørtoft pair. If conserved, it would constrain chirality or...

  35. [35]

    Energy and enstro- phy leak through horizons, escape to infinity, or exchange with matter

    When the system is open In GR the radiative sector is open. Energy and enstro- phy leak through horizons, escape to infinity, or exchange with matter. The balance laws are ˙W=I W − DW − FW , ˙Z=I Z − DZ − FZ +R bulk. (F8) The Fjørtoft constraint is active only if these open- system terms are slow compared to the nonlinear trans- fer: ˙W W ≪τ −1 nl , ˙Z Z ...

    work page

  36. [36]

    In the super- radiant regime 0< ω < mΩH , the canonical energy flux through the horizon is negative—the black hole feeds en- ergy into the mode

    Superradiance: whenW k is not positive The Fjørtoft argument assumesW k ≥0. In the super- radiant regime 0< ω < mΩH , the canonical energy flux through the horizon is negative—the black hole feeds en- ergy into the mode. The pair (W,Z) no longer satisfies the closed algebra. The fix, described in Section VI D, is to use the co-rotating energyW (χ) associa...

    work page

  37. [37]

    (42), with the fluid-weighted tensor of Eq

    W eighted Magnetic W eyl T ensor On the ingoing EF slicesv= const of the metric Eq. (42), with the fluid-weighted tensor of Eq. (46), ˜Bab = ⋆Cacbd ucud = 1 2 ϵac ef Cef bducud ,(G1) whereu a =u µ∂µ is the bulk extension of the boundary- normalized fluid velocity, kept unnormalized in the bulk (|u|2 =r 2fon the brane). The unit-normalized alterna- tives a...

    work page

  38. [38]

    (49) together with the incompressibility iden- tity R σijσij d2x= 1 2 R ω2 fl d2x(decaying or periodic data) factorizes it into Eq

    V orticity and Shear contributions The bulk enstrophy perturbation is δ ˜Z= Z ∞ rh dr Z δ ˜Bab δ ˜Bab √σ d2x ,(G3) and substituting the two-sector decomposition Eq. (49) together with the incompressibility iden- tity R σijσij d2x= 1 2 R ω2 fl d2x(decaying or periodic data) factorizes it into Eq. (52), with ˜Cω,σ (T) = Z ∞ rh ˜Kω,σ (r;T)r 2 dr .(G4) The ra...

    work page

  39. [39]

    Spacelike slices and the observer choice With a Cauchy foliation the fluid-gravity correspon- dence has also been demonstrated [56, 68]. To link black- hole behavior to the asymptotic fluid, one can conve- niently choose a horizon-penetrating deformation of the EF slicing, τ=v−h(r),0< h ′(r)< 2 r2f(br) ,(G7) for which the induced radial metric componentγ ...

    work page

  40. [40]

    (52), and can be extracted by repeating the evaluation with the full first-order metric includingds 2 grad

    The closed forms above are thek→0 profiles; for a boundary velocity modev(x) =v 0eik·x with ωfl(k) =iϵ ijkiv0,jeik·x, finite-wavenumber correc- tions enter atO(∂ 2v), consistently with Eq. (52), and can be extracted by repeating the evaluation with the full first-order metric includingds 2 grad

    work page

  41. [41]

    Extraction at boundary timetpairs Ω CFT(t) with the bulk integral on the EF slicev=t, which reaches conformal infinity at that boundary time. Appendix H: Quantitative checks: shear, horizon suppression, and the spectral ratio The enstrophy conservation argument rests on three quantitative claims: the ZAMO shear is small, horizon absorption is suppressed n...

    work page

  42. [42]

    Proposition 2 eliminates the radial compo- nentσ ˆrˆϕ, butσ ˆθ ˆϕ remains; we now show its effect is small

    Smallness of the ZAMO shear The ZAMO shear couples to the magnetic enstrophy throughσ ˆθ ˆϕ. Proposition 2 eliminates the radial compo- nentσ ˆrˆϕ, butσ ˆθ ˆϕ remains; we now show its effect is small. In Boyer–Lindquist coordinates the frame-dragging an- gular velocity is Ωfd =− gtϕ gϕϕ = 2M ar A , A= (r 2 +a 2)2 −a 2∆ sin2 θ, (H1) with ∆ =r 2 −2M r+a 2. ...

    work page

  43. [43]

    The fractional enstrophy loss through the horizon is F Z H+ δZ/τ nl ∼ ω2 Rκ+ ω2 R ·ϵ 2ωR = κ+ ϵ2ωR → κ+ mΩH (ϵ∼O(1))

    Suppression of horizon absorption For a ZDM withω≃mΩ H , the horizon energy flux vanishes at superradiant threshold: F W H+ ∝(ω R −mΩ H )|A|2 ∼κ +|A|2,(H5) whereAis the mode amplitude andω R −mΩ H ∼κ +. The fractional enstrophy loss through the horizon is F Z H+ δZ/τ nl ∼ ω2 Rκ+ ω2 R ·ϵ 2ωR = κ+ ϵ2ωR → κ+ mΩH (ϵ∼O(1)). (H6) The exact ratio is κ+ ΩH = √ M ...

    work page

  44. [44]

    But the interpretation of the near-horizon correction needs care

    The spectral ratioδZ k/Wk =ω 2 k: where it holds and where it doesn’t The identityδZ k/Wk =ω 2 k follows from the temporal structure of the mode, not from any radiation-zone ap- proximation. But the interpretation of the near-horizon correction needs care. a. Radiation-zone proof Decompose Ψ (1) 4 =R lm(r) −2Slm aω (θ)e −iωkt+imϕ. The Isaacson energy is W...

    work page

  45. [45]

    For ZDMs,ω ℓmn ≃mΩ H − κ+[cℓm +i(n+ 1 2 )], withc ℓm of order unity [42]

    Three-wave resonance is kinematically open The three-wave resonance requiresω 3 ≈ω 1 +ω 2 within the mode linewidths. For ZDMs,ω ℓmn ≃mΩ H − κ+[cℓm +i(n+ 1 2 )], withc ℓm of order unity [42]. For an azimuthally allowed triad (m1 +m2 =m 3), the ΩH terms cancel, and both the detuning ∆ =|Re(ω 1 +ω2 −ω3)|and 40 the combined linewidth Γ = P i |Imω i|scale asκ...

    work page

  46. [46]

    Conditions for the cascade to develop The Fjørtoft constraint can act only whenτ damp ≫ τnl; the form it takes is set by the dominant interaction (Sec. V C). Forl=m= 2 ZDMs: τ −1 nl ∼ϵ 2 ωR =ϵ 2 mΩH , τ damp ∼κ −1 + ≃ M √ 2√1−χ . (H17) The separation ratio is τdamp τnl = ϵ2 mΩH κ+ ≃ ϵ2 mp 2(1−χ) .(H18) Forl=m= 2, the onset conditionτ damp/τnl ≥1 gives ϵ2 ...

    work page

  47. [47]

    Fjørtoft, Tellus5, 225 (1953)

    R. Fjørtoft, Tellus5, 225 (1953)

    work page 1953

  48. [48]

    R. H. Kraichnan and D. Montgomery, Reports on Progress in Physics43, 547 (1980)

    work page 1980

  49. [49]

    H. Yang, A. Zimmerman, and L. Lehner, Physical Review Letters114, 081101 (2015), arXiv:1402.4859 [gr-qc]

    Pith/arXiv arXiv 2015

  50. [50]

    Iuliano, S

    C. Iuliano, S. Hollands, S. R. Green, and P. Zimmerman, Extremal black hole weather (2024), arXiv:2412.02821 [gr-qc]

    Pith/arXiv arXiv 2024

  51. [51]

    S. Ma, L. Lehner, H. Yang, L. E. Kidder, H. P. Pfeiffer, and M. A. Scheel, Physical Review Letters136, 061401 (2026), arXiv:2508.13294 [gr-qc]

    arXiv 2026

  52. [52]

    Carrasco, L

    F. Carrasco, L. Lehner, R. C. Myers, O. Reula, and A. Singh, Physical Review D86, 126006 (2012), arXiv:1210.6702 [hep-th]. •χ= 0.9999:ϵ≳0.084 5 Forχ≳0.999 then= 1 modes must be tracked onto the ZDM branch directly, Imω≃ −(n+ 1 2 )κ+: continued-fraction con- tinuation in spin can hop onto an adjacent overtone-multiplet branch with Imω≃ − 5 2 κ+, which woul...

    Pith/arXiv arXiv 2012

  53. [53]

    Adams, P

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

    Pith/arXiv arXiv 2014

  54. [54]

    Eling and Y

    C. Eling and Y. Oz, JHEP11, 079, arXiv:1308.1651 [hep- th]

    Pith/arXiv arXiv

  55. [55]

    Enstrophy from symmetry

    R. Marjieh, N. Pinzani-Fokeeva, and A. Yarom, Enstro- phy from symmetry (2021), arXiv:2009.03980 [hep-th]

    work page internal anchor Pith review Pith/arXiv arXiv 2021

  56. [56]

    Donnelly, L

    W. Donnelly, L. Freidel, S. F. Moosavian, and A. J. Sper- anza, JHEP09, 008, arXiv:2012.10367 [hep-th]

    Pith/arXiv arXiv 2012

  57. [57]

    Redondo-Yuste and L

    J. Redondo-Yuste and L. Lehner, JHEP02, 240, arXiv:2212.06175 [gr-qc]

    Pith/arXiv arXiv

  58. [58]

    R. A. Isaacson, Physical Review166, 1272 (1968)

    1968

  59. [59]

    Kinnersley, J

    W. Kinnersley, J. Math. Phys.10, 1195 (1969)

    1969

  60. [60]

    Newman and R

    E. Newman and R. Penrose, J. Math. Phys.3, 566 (1962)

    1962

  61. [61]

    S. R. Green, F. Carrasco, and L. Lehner, Physical Review X4, 011001 (2014), arXiv:1309.7940 [hep-th]

    Pith/arXiv arXiv 2014

  62. [62]

    J. R. Westernacher-Schneider, L. Lehner, and Y. Oz, JHEP12, 067, arXiv:1510.00736 [hep-th]

    work page internal anchor Pith review Pith/arXiv arXiv

  63. [63]

    J. M. Bardeen, W. H. Press, and S. A. Teukolsky, The Astrophysical Journal178, 347 (1972)

    1972

  64. [64]

    Bel, Comptes Rendus de l’Acad´ emie des Sciences247, 1094 (1958)

    L. Bel, Comptes Rendus de l’Acad´ emie des Sciences247, 1094 (1958)

    work page 1958

  65. [65]

    C. W. Misner, K. S. Thorne, and J. A. Wheeler,Gravi- tation(W. H. Freeman, San Francisco, 1973)

    1973

  66. [66]

    Maartens and B

    R. Maartens and B. A. Bassett, Classical and Quantum Gravity15, 705 (1998), arXiv:gr-qc/9704059

    Pith/arXiv arXiv 1998

  67. [67]

    Friedrich, Classical and Quantum Gravity13, 1451 (1996)

    H. Friedrich, Classical and Quantum Gravity13, 1451 (1996). 41

    work page 1996

  68. [68]

    H. Yang, F. Zhang, S. R. Green, and L. Lehner, Phys. Rev. D91, 084007 (2015), arXiv:1502.08051 [gr-qc]

    Pith/arXiv arXiv 2015

  69. [69]

    J. L. Trueba and A. F. Ra˜ nada, European Journal of Physics17, 141 (1996)

    work page 1996

  70. [70]

    R. A. Isaacson, Physical Review166, 1263 (1968), part II:Nonlinear terms and the effective stress tensor, Phys. Rev.166, 1272 (1968), doi:10.1103/PhysRev.166.1272

    doi:10.1103/physrev.166.1272 1968

  71. [71]

    Maggiore,Gravitational Waves: Volume 1: Theory and Experiments(2007)

    M. Maggiore,Gravitational Waves: Volume 1: Theory and Experiments(2007)

    work page 2007

  72. [72]

    Hollands and R

    S. Hollands and R. M. Wald, Commun. Math. Phys.321, 629 (2013), arXiv:1201.0463 [gr-qc]

    Pith/arXiv arXiv 2013

  73. [73]

    Canonical Energy and Hertz Potentials for Perturbations of Schwarzschild Spacetime

    K. Prabhu and R. M. Wald, Classical and Quantum Gravity35, 235004 (2018), arXiv:1807.09883 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  74. [74]

    S. A. Teukolsky and W. H. Press, The Astrophysical Journal193, 443 (1974)

    work page 1974

  75. [75]

    J. N. Goldberg and R. K. Sachs, Acta Physica Polonica22 (Suppl.), 13 (1962), reprinted as a Golden Oldie: Gen. Relativ. Gravit.41, 433–444 (2009), doi:10.1007/s10714-008-0722-5

    work page doi:10.1007/s10714-008-0722-5 1962

  76. [76]

    Bel, General Relativity and Gravitation32, 2047 (2000), golden Oldie reprint of the 1962 original

    L. Bel, General Relativity and Gravitation32, 2047 (2000), golden Oldie reprint of the 1962 original

    work page 2047

  77. [77]

    J. M. M. Senovilla, Classical and Quantum Gravity17, 2799 (2000), arXiv:gr-qc/9906087

    Pith/arXiv arXiv 2000

  78. [78]

    Dynamical laws of superenergy in General Relativity

    A. Garc ´ ıa-Parrado G´ omez-Lobo, Classical and Quantum Gravity25, 015006 (2008), arXiv:0707.1475 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2008

  79. [79]

    D. A. Nichols, R. Owen, F. Zhang, A. Zimmerman, J. Brink, Y. Chen, J. D. Kaplan, G. Lovelace, K. D. Matthews, M. A. Scheel, and K. S. Thorne, Physical Re- view D84, 124014 (2011), arXiv:1108.5486 [gr-qc]

    Pith/arXiv arXiv 2011

  80. [80]

    Chandrasekhar,The Mathematical Theory of Black Holes(Oxford University Press, New York, 1983)

    S. Chandrasekhar,The Mathematical Theory of Black Holes(Oxford University Press, New York, 1983)

    1983

Showing first 80 references.

This paper was first reviewed by deepseek-v4-flash on August 5, 2026.