2025-11-13T08:52:11.146722

An observational signature for extremal black holes

Aretakis, Khanna, Sabharwal
We consider scalar perturbations of the Reissner--Nordström family and the Kerr family. We derive a characteristic expression of the radiation field, at any given unit solid angle of future null infinity, and numerically show that its amplitude gets excited only in the extremal case. Our work, therefore, identifies an observational signature for extremal black holes. Moreover, we show that the source of the excitation is the extremal horizon instability and its magnitude is exactly equal to the conserved horizon charge.
academic

An observational signature for extremal black holes

Basic Information

  • Paper ID: 2307.03963
  • Title: An observational signature for extremal black holes
  • Authors: Stefanos Aretakis (University of Toronto), Gaurav Khanna (University of Rhode Island & UMass Dartmouth), Subir Sabharwal (University of Rhode Island)
  • Classification: gr-qc (General Relativity and Quantum Cosmology), hep-th (High Energy Physics - Theory), math-ph (Mathematical Physics)
  • Publication Date: July 2023 (arXiv v3: November 7, 2025)
  • Paper Link: https://arxiv.org/abs/2307.03963

Abstract

This paper investigates scalar perturbations of Reissner-Nordström and Kerr family black holes. The authors derive characteristic expressions for the radiation field at arbitrary fixed angles on future null infinity I⁺, and demonstrate through numerical methods that the amplitude is excited only in the extremal black hole case. Consequently, this work identifies an observational signature for extremal black holes. Furthermore, the authors prove that the excitation source is the extremal horizon instability, with magnitude precisely equal to the conserved horizon charge.

Research Background and Motivation

1. Core Problem to be Addressed

This paper aims to solve the problem of identifying observational signatures of extremal black holes, namely how to distinguish extremal from sub-extremal black holes through perturbation dynamics measurable by distant observers.

2. Problem Significance

Extremal black holes hold special status in theoretical physics:

  • Theoretical Significance: They saturate geometric inequalities for mass, angular momentum, and charge, playing important roles in Hawking radiation and quantum gravity research
  • Astrophysical Significance: Observational evidence suggests near-extremal black holes are ubiquitous in the universe, including stellar-mass and supermassive near-extremal black holes
  • Fundamental Physics Significance: Extremal black holes have zero temperature and possess conformally invariant near-horizon geometry, yielding new solutions to Einstein's equations

3. Limitations of Existing Methods

Previous work (Angelopoulos, Aretakis & Gajic 2018) proposed an observational signature requiring integration over the entire S² sphere, which is difficult to implement in practical observations. This paper proposes a signature requiring only a single detector at a fixed angle on null infinity, more closely approximating actual physical observation conditions.

4. Research Motivation

  • Extremal black holes exhibit horizon instability with transverse derivatives growing along the event horizon
  • Existence of a conserved horizon charge Hψ, which challenges the "no-hair theorem"
  • Need to verify whether the horizon charge Hψ can be measured by external observers, potentially falsifying the no-hair assumption for extremal black holes

Core Contributions

  1. Proposes new observational signature expression sI+(ϑ): Measuring the radiation field at a fixed angle ϑ on future null infinity I⁺ alone suffices to distinguish extremal from sub-extremal black holes
  2. Establishes theoretical conjecture: For extremal Reissner-Nordström (ERN) black holes, derives explicit expression for subleading coefficient C₂(ϑ):
    • C₂(ϑ) = -32M²Hψ + 8MI⁽¹⁾ψ
  3. Numerical verification: Through high-precision numerical simulations verifies:
    • Sub-extremal case: sI+(ϑ) → 0
    • Extremal case: sI+(ϑ) = Hψ
  4. Generalization to Kerr black holes: Successfully extends results to extremal Kerr (EK) black holes, verifying similar observational signatures
  5. Establishes physical connection: Proves the observational signature originates from extremal horizon instability, with magnitude precisely equal to the conserved horizon charge

Methodology Details

Task Definition

Input: Initial data ψ (scalar field) for black hole perturbations
Output: Observational signature sI+(ϑ) for determining whether the black hole is extremal
Constraint: Observer located at fixed angle ϑ on future null infinity I⁺, measuring only radiation field in that direction

Theoretical Framework

1. Geometric Setup

  • Future null infinity I⁺ with coordinates: I⁺ = (τ, ϑ)
    • τ: retarded time (on I⁺) or advanced time (on H⁺)
    • ϑ = (θ₁, θ₂): fixed angular direction of observer
  • Null generators parameterized as: {(τ, ϑ₀), τ ∈ ℝ}

2. Radiation Field Definition

Scalar field ψ satisfies the wave equation: □_g ψ = 0

Radiation field Ψ defined as:

Ψ(τ, ϑ) = lim_{r→∞} (rψ)(τ, r, ϑ)

3. Tail Asymptotic Behavior

Sub-extremal RN (known result):

Ψ(τ, ϑ) = C₁(ϑ)/τ² + C₂(ϑ)log(τ/M)/τ³ + O(1/τ³)

where:

  • C₁(ϑ) = -2I⁽¹⁾ψ
  • C₂(ϑ) = 8MI⁽¹⁾ψ
  • C₂(ϑ)/C₁(ϑ) = -4M (key ratio)

Extremal RN (conjecture in this paper):

  • C₁(ϑ) = 4MHψ - 2I⁽¹⁾ψ (known)
  • C₂(ϑ) = -32M²Hψ + 8MI⁽¹⁾ψ (conjecture in this paper)

4. Conserved Quantity Definitions

Horizon charge:

H[ψ] := -M²/(4π) ∫_{EH} ∂_r(rψ) dΩ

Newman-Penrose constant:

I⁽¹⁾[ψ] = M/(4π) ∫_{I⁺∩{τ≥0}} Ψ(τ, ϑ)dϑdτ

Observational Signature Construction

For RN Black Holes

Extract coefficients from radiation field:

C₁(ϑ) = lim_{τ→∞} τ² · Ψ(τ, ϑ)

C₂(ϑ) = lim_{τ→∞} [τ/log(τ/M)] · [τ² · Ψ(τ, ϑ) - C₁(ϑ)]

Core observational signature:

sI+(ϑ) = -1/(16M²) (C₂(ϑ) + 4MC₁(ϑ))

Equivalent form:

sI+(ϑ) = -1/(16M²) · lim_{τ→∞} [4Mτ²Ψ(τ,ϑ) - τ/log(τ/M) · lim_{τ→∞}(τ²Ψ(τ,ϑ)) + τ³Ψ(τ,ϑ)/log(τ/M)]

Key properties:

sI+(ϑ) = { 0,      sub-extremal RN
          { H[ψ],  extremal RN

For Kerr Black Holes

Similar definition with additional factor 1/4:

s^Kerr_I+(ϑ) = -1/(64M²) (C₂(ϑ) + 4MC₁(ϑ))

Maintains properties:

s^Kerr_I+(ϑ) = { 0,      sub-extremal Kerr
                { H[ψ],  extremal Kerr

Technical Innovations

  1. Anisotropic measurement: No sphere integration needed; single-point measurement suffices
  2. Subleading term utilization: Cleverly exploits information about extremality carried by logarithmic coefficient C₂(ϑ)
  3. Universality: Applicable to both RN and Kerr black hole families
  4. Physical measurability: Directly connects horizon "hair" Hψ with far-field observations

Experimental Setup

Numerical Methods

1. Coordinate System

Employs compactified hyperbolic coordinates (τ, ρ, θ₁, θ₂):

  • Penetrates the horizon
  • Brings I⁺ to finite radial coordinate ρ(I⁺) = S < ∞
  • Coordinate transformation:
Ω = 1 - ρ/S
r = ρ/Ω(ρ)
v ≡ t + r* - r = τ + ρ/Ω(ρ) - ρ - 4M log Ω(ρ)

Parameter choices:

  • S = 19.0 (I⁺ location)
  • ρ(H⁺) = 0.95 (horizon location)

2. Numerical Solution Scheme

  • PDE rewriting: Reformulates second-order PDE as two coupled first-order differential equations
  • Spatial discretization: High-order weighted essentially non-oscillatory (WENO) finite difference scheme
  • Time integration: Explicit Shu-Osher time stepping
  • Dimensionality:
    • RN case: (1+1)D system (spherical harmonic decomposition)
    • Kerr case: (2+1)D system

3. Computational Resources

  • GPU parallel computing (GPGPU architecture)
  • RN simulations: ~1 week per simulation
  • Kerr simulations: ~1 month per simulation
  • Platforms: UMass-URI UNITY supercomputer and MIT SuperCloud

Initial Data

Type A initial data (astrophysically relevant):

  • Penetrates horizon (H⁺)
  • Compactly supported
  • No support on I⁺

Specific form: Truncated Gaussian pulse

  • Center positions: ρ = 1.0, 1.1, 1.2
  • Widths: 0.16, 0.22, 0.32
  • Non-zero region: ρ ∈ 0.95, 8
  • Angular distribution: Spherically symmetric (ℓ = 0 monopole)

Physical reasonableness:

  • Higher multipoles decay rapidly
  • Late times dominated by ℓ = 0 (consistent with Price's law)
  • Can be viewed as "generic" initial data

Data Analysis Methods

Sub-extremal Case

Uses local power index (LPI):

LPI ≡ τψ̇/ψ

Relation:

τ[τψ̇/ψ + 2] = -C₂(ϑ)/C₁(ϑ) · log(τ/M) + C₂(ϑ)/C₁(ϑ) + O(M²log(τ/M)/τ)

Obtain ratio C₂(ϑ)/C₁(ϑ) through regression of τ × (LPI + 2) versus log(τ/M)

Smoothing algorithm:

  • Simple moving average
  • Dynamic window size
  • Minimizes variance while preserving kurtosis

Extremal Case

  • Evolve to very late times: t/M = 8000
  • Late-time data: τ/M ≳ 5000
  • Directly fit C₁(ϑ) and C₂(ϑ) using Levenberg-Marquardt algorithm
  • Avoids numerical noise introduced by time derivatives

Experimental Results

Main Results

1. Sub-extremal RN Black Holes (Figure 1)

Test three charge ratios: Q/M = 0.3, 0.4, 0.5

Observations:

  • At late times (τ/M ≳ 3300), sI+(ϑ) asymptotically approaches 0
  • Verifies known theoretical result: C₂(ϑ) = -4MC₁(ϑ)
  • Numerical precision: sI+(ϑ) ~ O(10⁻²)

Significance:

  • Validates numerical method correctness
  • Establishes code credibility

2. Extremal RN Black Holes (Figure 2)

Uses 5 different initial data sets:

  • Gaussian pulses with different positions and widths
  • One set with zero initial field and Gaussian momentum

Key findings:

  • sI+(ϑ) exhibits linear relationship with Hψ
  • Best-fit line slope: 1.0 ± 0.2
  • Includes origin (0,0) as constraint (sI+(ϑ) must vanish when Hψ=0)

Numerical performance:

  • Levenberg-Marquardt fit errors negligible (imperceptible in figure)
  • Consistent results across different initial data
  • Verifies conjecture: sI+(ϑ) = Hψ

3. Sub-extremal Kerr Black Holes (Figure 3)

Test three spin ratios: a/M = 0.3, 0.4, 0.5

Observations:

  • At late times (τ/M ≳ 3300), s^Kerr_I+(ϑ) asymptotically approaches 0
  • Behavior similar to RN case
  • Validates method generalization effectiveness

4. Extremal Kerr Black Holes (Figure 4)

Uses 3 different initial data sets

Key findings:

  • s^Kerr_I+(ϑ) exhibits linear relationship with Hψ
  • Best-fit line slope: 1.0 ± 0.1
  • Error bars visible but small

Significance:

  • Successfully generalizes to rotating black holes
  • Observational signature demonstrates universality

Ablation Experiments

Initial Data Without Horizon Support (Figure 5, Appendix)

  • Gaussian pulse center: ρ = 4.0
  • Width: 0.22
  • Truncation window: 1.0, 8.0
  • Horizon location: ρ = 0.95 (no overlap)

Results:

  • sI+(ϑ) approaches 0 for both ERN and EK cases
  • Verifies observational signature vanishes when Hψ = 0
  • Completely consistent with theoretical expectations

Numerical Precision Analysis

Error sources:

  1. Accumulated truncation errors from derivatives
  2. Noise introduced by data analysis methods
  3. Finite evolution time

Precision limitations:

  • Sub-extremal case: sI+(ϑ) ~ O(10⁻²) (should be 0)
  • Reliable detection threshold: Hψ ≫ 10⁻³
  • Sufficient for "principle-of-concept" verification, but not constrained by current numerical limitations

Statistical significance:

  • Extremal case linear fit: R² close to 1
  • Slope confidence interval: ±0.1 to ±0.2
  • Highly consistent with theoretical prediction (slope=1)

Related Work

1. Extremal Black Hole Stability Theory

  • Aretakis (2011-2015) 5-7: First discovered extremal black hole horizon instability with transverse derivatives growing along horizon
  • Lucietti & Reall (2012) 19: Gravitational instability of extremal Kerr black holes
  • Casals et al. (2016) 10: Non-axisymmetric modes and enhanced growth rates

2. Asymptotic Behavior Analysis

  • Angelopoulos, Aretakis & Gajic (2019) 2: Logarithmic corrections to radiation field for sub-extremal RN black holes
  • Angelopoulos, Aretakis & Gajic (2020) 3: Late-time asymptotics of wave equation on extremal RN backgrounds
  • Price (1986) 4: Price's law, tail behavior for spherically symmetric black holes

3. Observational Signature Research

  • Angelopoulos, Aretakis & Gajic (2018) 8: Horizon "hair" and measurements at null infinity (requires sphere integration)
  • Gralla, Hughes & Warburton (2016) 16: Gravitational signatures from near-horizon geometry
  • Gralla, Lupsasca & Strominger (2018) 17: High-spin observational signatures for Event Horizon Telescope

4. Astrophysical Observations

  • Brenneman et al. (2011) 9: Spin of supermassive black hole in NGC 3783
  • Volonteri et al. (2005) 23: Spin distribution of massive black holes and cosmic evolution
  • Sadowski et al. (2011) 20: Possibility of supercritical accretion flows spinning black holes beyond Thorne limit

5. Quasinormal Modes and Superradiance

  • Yang et al. (2013) 24: Quasinormal mode spectrum bifurcation near-extremal Kerr spacetime
  • Glampedakis & Andersson (2001) 15: Late-time dynamics for rapidly rotating black holes
  • Gajic (2023) 14: Azimuthal instability on extremal Kerr

Advantages of This Work

  1. Single-point measurement: Compared to 8 requiring sphere integration, this work requires only fixed angle
  2. Numerical verification: Provides detailed numerical evidence rather than pure theoretical derivation
  3. Universality: Applicable to both RN and Kerr families
  4. Practicality: Closer to actual observational conditions

Conclusions and Discussion

Main Conclusions

  1. Theoretical contributions:
    • Proposes anisotropic observational signature sI+(ϑ) for extremal black holes
    • Conjectures and numerically verifies explicit expression for subleading coefficient C₂(ϑ)
    • Establishes direct connection between far-field observations and horizon charge Hψ
  2. Numerical verification:
    • Sub-extremal case: sI+(ϑ) = 0 (RN and Kerr)
    • Extremal case: sI+(ϑ) = Hψ (RN and Kerr)
    • Linear relationship slope 1.0 with error ±0.1 to ±0.2
  3. Physical significance:
    • Horizon "hair" Hψ can be measured by external observers
    • Potentially falsifies no-hair assumption for extremal black holes (within linear perturbation theory framework)
    • Excitation source is extremal horizon instability

Limitations

  1. Theoretical framework constraints:
    • Applicable only to black hole perturbation theory (fixed background spacetime)
    • Not applicable to full nonlinear Einstein field equations
    • Idealized isolated black hole model
  2. Numerical precision:
    • Sub-extremal case residual error ~ O(10⁻²)
    • Reliable detection threshold: Hψ ≫ 10⁻³
    • Long-time evolution accumulates errors
  3. Computational cost:
    • RN: 1 week/simulation
    • Kerr: 1 month/simulation
    • Near-extremal cases not attempted
  4. Initial data limitations:
    • Only considers axisymmetric modes (m = 0)
    • Type A initial data (astrophysically relevant)
    • Spherically symmetric angular distribution (ℓ = 0 dominated)
  5. Observational practical challenges:
    • Requires extremely long observation times (τ/M ≳ 5000)
    • Requires extremely high precision measurement of subleading terms
    • Complexity of actual astrophysical environments

Future Directions

  1. Theoretical extensions:
    • Transient behavior in near-extremal cases (charge/spin within 5×10⁻⁵M)
    • Non-axisymmetric perturbations (m ≠ 0)
    • Gravitational wave perturbations (non-scalar fields)
  2. Numerical improvements:
    • Enhance numerical precision, reduce errors
    • More efficient algorithms to lower computational cost
    • Explore near-extremal parameter space
  3. Physical applications:
    • Combine with Event Horizon Telescope observations
    • Applications in gravitational wave detection
    • Search for black holes exceeding Thorne limit
  4. Nonlinear effects:
    • Full Einstein equation numerical simulations
    • Fate of horizon "hair" in nonlinear dynamics

In-Depth Evaluation

Strengths

  1. Strong innovation:
    • First proposes single-point measurement extremality criterion
    • Cleverly exploits information in subleading logarithmic terms
    • Connects abstract horizon charge to observable quantities
  2. Rigorous methodology:
    • Theoretical derivations based on established rigorous mathematical results
    • Advanced numerical methods (WENO scheme, GPU parallelization)
    • Cross-validation with multiple initial data sets
  3. Convincing results:
    • Sub-extremal case perfectly matches known theory
    • Extremal case linear fit slope precisely 1
    • Consistent results for both RN and Kerr families
  4. Deep physical insights:
    • Reveals far-field manifestation of extremal horizon instability
    • New perspective challenging no-hair assumption
    • Connects near-horizon physics with far-field observations
  5. Clear presentation:
    • Logical structure
    • Precise mathematical expressions
    • Clear physical picture

Weaknesses

  1. Theoretical completeness:
    • C₂(ϑ) expression remains conjecture lacking rigorous mathematical proof
    • Relies on applicability assumptions of Price's law
    • Validity range of linear approximation insufficiently discussed
  2. Numerical limitations:
    • Relatively large sub-extremal residual error (10⁻² magnitude)
    • Near-extremal transient behavior not explored
    • High computational cost limits parameter space exploration
  3. Practical applicability challenges:
    • Observation timescale extremely long (thousands of M)
    • Requires measuring subleading terms (high technical difficulty)
    • Applicability in actual astrophysical environments questionable
  4. Insufficient comparative experiments:
    • Lacks detailed comparison with other observational signature schemes
    • Missing sensitivity analysis to observational errors
    • Insufficient discussion of complementarity with gravitational wave observations
  5. Physical interpretation:
    • Insufficient explanation for why Kerr requires additional 1/4 factor
    • Limited discussion of physical meaning of different initial data
    • Insufficient connection with string theory and quantum gravity

Impact

  1. Academic value:
    • Provides new tools for extremal black hole research
    • Advances study of horizon "hair" observability
    • Inspires reconsideration of no-hair assumption
  2. Practical value:
    • Provides theoretical guidance for future black hole observations
    • Potential application to searching for super-Thorne-limit black holes
    • Possibility of combining with multi-messenger astronomy
  3. Reproducibility:
    • Detailed method description
    • Numerical techniques supported by prior work
    • High computational cost limits independent verification
  4. Field contribution:
    • Successful combination of mathematical relativity and numerical relativity
    • Demonstrates power of GPU computing in gravitational physics
    • Provides new perspective for extremal black hole astrophysics

Applicable Scenarios

  1. Theoretical research:
    • Extremal black hole perturbation theory
    • Horizon dynamics research
    • No-hair theorem testing
  2. Numerical simulation:
    • Long-time evolution of black hole perturbations
    • Asymptotic behavior analysis of radiation fields
    • High-precision numerical relativity
  3. Observational astronomy (long-term potential):
    • Near-extremal black hole identification
    • Precise black hole parameter measurement
    • Extreme physics environment detection
  4. Inapplicable scenarios:
    • Short timescale dynamical processes
    • Strong nonlinear gravitational effects
    • Complex astrophysical environments (accretion disks, companions, etc.)

Selected References

2 Angelopoulos, Y., Aretakis, S., and Gajic, D. (2019). Logarithmic corrections in the asymptotic expansion for the radiation field along null infinity. Journal of Hyperbolic Differential Equations.

3 Angelopoulos, Y., Aretakis, S., and Gajic, D. (2020). Late-time asymptotics for the wave equation on extremal Reissner–Nordström backgrounds. Advances in Mathematics.

5-7 Aretakis, S. (2011-2015). Stability and instability of extreme Reissner–Nordström black hole spacetimes series. Pioneering work on extremal black hole horizon instability.

8 Angelopoulos, Y., Aretakis, S., and Gajic, D. (2018). Horizon Hair of Extremal Black Holes and Measurements at Null Infinity. Phys. Rev. Lett. Prior observational scheme requiring sphere integration.

13 Burko, L.M., Khanna, G., and Sabharwal, S. (2019). (Transient) Scalar Hair for (Nearly) Extreme Black Holes. Phys. Rev. Research. Related work on near-extremal cases.

Overall Assessment: This is an important contribution to the field of observational signatures of extremal black holes. Through combination of theoretical derivation and high-precision numerical simulation, the authors successfully propose and verify a practical extremality criterion. Despite challenges in theoretical completeness and practical application, this work provides important theoretical foundation and technical pathway for understanding extremal black hole physics and its potential observations, with significant implications for both mathematical relativity and astrophysics.