Extending the theory-agnostic Blandford-Znajek luminosity to near-extremal rotation and validating it against nonlinear simulations in Konoplya-Rezzolla-Zhidenko black-hole spacetimes

The Blandford-Znajek (BZ) mechanism, in which magnetic fields threading a spinning black hole extract its rotational energy and power relativistic jets, has recently been given a theory-agnostic formulation: Camilloni and Rezzolla (Phys. Rev. D 114, 023014, 2026) showed, within the horizon-penetrating Konoplya-Rezzolla-Zhidenko (KRZ) parametrization of generic stationary and axisymmetric black-hole spacetimes, that the leading-order jet luminosity scales universally as $\Omega_h^2$, independent of the spacetime-deviation parameters, while the next-leading-order corrections depend on the KRZ coefficients, breaking the degeneracy for rapidly rotating holes. The analysis, however, is restricted to a split-monopole magnetic topology, is polynomial through $O(\Omega_h^6)$ and so may fail near extremal spin $a_*=1$, and has not been validated against nonlinear simulations. This record poses the completion of that theory-agnostic BZ framework: extend the luminosity formulae to the near-extremal regime with non-polynomial matched-asymptotic terms, validate them against fully nonlinear force-free or general-relativistic magnetohydrodynamic (GRMHD) simulations in generic KRZ spacetimes, and establish the horizon-scale magnetic-field topology beyond the split monopole. These are the joint theoretical prerequisites for using BZ jet power as a probe of strong-field gravity and a test of the Kerr paradigm.

Background

The Blandford-Znajek (BZ) mechanism is the standard model for how spinning black holes launch relativistic jets: magnetic field lines threading the event horizon are twisted by the hole's rotation and extract its rotational energy as an outgoing Poynting flux. The original calculation (Blandford & Znajek 1977) treated the black-hole magnetosphere as force-free, meaning the electromagnetic stresses dominate the plasma inertia so that the electric and magnetic fields are perpendicular, and adopted a split-monopole geometry, a purely radial magnetic field with an equatorial current sheet. In this setting the governing equation for the magnetic flux function is the Grad-Shafranov equation, a nonlinear second-order partial differential equation, which Blandford and Znajek solved perturbatively in the horizon angular velocity $\Omega_h$ (the angular velocity at which the horizon is dragged around by frame dragging), finding a jet luminosity of schematic form $P_{\rm BZ}\sim\Omega_h^2 f(\Omega_h)$ with $f\simeq 1$ for slowly rotating holes.

Whether this picture survives beyond the Kerr metric of general relativity matters for the Kerr paradigm, the assumption that astrophysical black holes are described by the Kerr solution. The Konoplya-Rezzolla-Zhidenko (KRZ) parametrization (Konoplya, Rezzolla & Zhidenko 2016) provides a general expansion of any stationary, axisymmetric black-hole metric as a continued-fraction series around the horizon whose dimensionless coefficients — at leading order $\varrho$, $a_1$, $b_1$ — measure deviations from Kerr, with Kerr recovered when all coefficients vanish. A horizon-penetrating (HP) form of these metrics, in which the coordinates cross the horizon smoothly so that numerical simulations can evolve fields through it (Ma & Rezzolla 2024), is the version adapted to magnetospheric work.

Camilloni & Rezzolla (2026) gave the first analytic extension of the BZ calculation to generic stationary and axisymmetric spacetimes using the HP-KRZ framework (improving on an earlier second-order attempt by Konoplya, Kunz & Zhidenko 2021). They found a universal leading order, $P_{\rm BZ}\sim\Omega_h^2$ with a spin-dependence function $f_{\rm KRZ}(\Omega_h)\simeq 1$ for $M\Omega_h\ll 1$ ($M$ being the black-hole mass) independent of the deviation parameters, so that slowly rotating black holes cannot be distinguished by their jet power; at next-leading order, $f_{\rm KRZ}$ — computed through $O(\Omega_h^4)$, giving the luminosity through $O(\Omega_h^6)$ — depends on $\varrho$, $a_1$ and $b_1$, breaking the degeneracy for rapidly rotating holes ($\Omega_h\gtrsim 0.3/M$). They also derived spin-dependent physical bounds on the KRZ parameters: not every HP-KRZ black hole can reach the extremality condition $a_*=1$ (with $a_*$ the dimensionless spin parameter $J/M^2$), and super-Kerr solutions can exceed the Kerr limit $2M\Omega_h=1$. For Kerr specifically, the analogous near-extremal loss of accuracy of the polynomial expansion was cured by discovering non-analytic contributions via matched asymptotic expansions (Camilloni et al. 2022) — the technique of separately expanding the solution in an inner near-horizon region and an outer region and matching the overlap — and by near-horizon force-free attractor solutions (Camilloni et al. 2021). These Kerr results do not transfer automatically, because in the KRZ family the extremal limit is a parameter-dependent surface rather than a point.

The 2026 analysis leaves the following gaps, all stated in the paper itself: the magnetic geometry is assumed to be a split monopole; the polynomial expansion through $O(\Omega_h^6)$ may be inaccurate as $a_*\to 1$; the analytic expressions have not been contrasted with fully nonlinear simulations; and, observationally, using BZ luminosity as a strong-gravity probe requires independent high-precision measurements of $\Omega_h$ together with reliable knowledge of the horizon-scale magnetic topology.

Problem Statement

Complete the theory-agnostic description of Blandford-Znajek energy extraction in generic stationary and axisymmetric black-hole spacetimes, as initiated by Camilloni and Rezzolla (2026) in the horizon-penetrating Konoplya-Rezzolla-Zhidenko (KRZ) framework, by resolving the limitations that the paper itself identifies. Specifically:

1. Near-extremal extension. Extend the perturbative BZ luminosity formulae — presently polynomial through $O(\Omega_h^6)$ and possibly inaccurate as the spin approaches $a_*=1$ — to the near-extremal regime: determine the near-extremal behaviour of the KRZ luminosity and whether the polynomial expansion must be supplemented by non-polynomial (non-analytic) contributions, candidates being the matched-asymptotic corrections known for Kerr and their generalization to the KRZ family. The extension must remain valid within the spin-dependent physical bounds of the $(\varrho, a_1, b_1)$ parameter space, where the extremal limit is itself a parameter-dependent surface and some KRZ black holes saturate the bounds before reaching $a_*=1$.

2. Nonlinear validation. Contrast the analytic expressions with fully nonlinear simulations — force-free, GRMHD (general-relativistic magnetohydrodynamics), or particle-in-cell (a kinetic plasma-simulation method) — in generic KRZ spacetimes, across a representative portion of the KRZ parameter space and spin range, thereby testing the perturbative framework and, where the discrepancies require it, improving the modeling in the regime where the present expansion becomes unreliable. The analytic expressions being tested and the acceptance tolerance must be fixed before the comparison is made, so that any such refinement is reported separately and does not enter the prediction that is being validated.

3. Magnetic topology. Establish, with high-accuracy nonlinear simulations, the actual magnetic-field topology on horizon scales in KRZ black-hole spacetimes, going beyond the assumed split monopole, and quantify its impact on the BZ efficiency relative to the split-monopole benchmark.

The common objective uniting these three deliverables is a spin-complete, simulation-validated, theory-agnostic model of BZ jet power. More specifically, the three deliverables are jointly required rather than alternatives: the programme is resolved when the near-extremal luminosity of (1) and the nonlinear simulations of (2) are shown to agree over at least the KRZ subregion $|\varrho|, |a_1|, |b_1| \le 0.1$, with a spin range that reaches the parameter-dependent extremal bound of that subregion, to a fractional tolerance fixed before the comparison and no looser than $10\%$, and when the horizon-scale topology of (3) is quantified by its impact on the BZ efficiency relative to the split-monopole benchmark. Together with independent high-precision measurements of the black-hole angular velocity $\Omega_h$ and of the horizon-scale magnetic topology — observational requirements acknowledged as still open — such a model would let the jet luminosity of rapidly rotating black holes serve as a probe of the strong-field spacetime and hence as a test of the Kerr paradigm.

References

1. F. Camilloni and L. Rezzolla (2026). Universality of the Blandford-Znajek emission in stationary and axisymmetric spacetimes. Physical Review D 114, 023014. doi:10.1103/8b1s-5mmm, arXiv:2602.23417. https://arxiv.org/abs/2602.23417
2. R. D. Blandford and R. L. Znajek (1977). Electromagnetic extraction of energy from Kerr black holes. Monthly Notices of the Royal Astronomical Society 179, 433-456. doi:10.1093/mnras/179.3.433. https://doi.org/10.1093/mnras/179.3.433
3. R. Konoplya, L. Rezzolla, and A. Zhidenko (2016). General parametrization of axisymmetric black holes in metric theories of gravity. Physical Review D 93, 064015. doi:10.1103/PhysRevD.93.064015, arXiv:1602.02378. https://arxiv.org/abs/1602.02378
4. F. Camilloni, O. J. C. Dias, G. Grignani, T. Harmark, R. Oliveri, M. Orselli, A. Placidi, and J. E. Santos (2022). Blandford-Znajek monopole expansion revisited: novel non-analytic contributions to the power emission. Journal of Cosmology and Astroparticle Physics 2022(07), 032. doi:10.1088/1475-7516/2022/07/032, arXiv:2201.11068. https://arxiv.org/abs/2201.11068
5. F. Camilloni, G. Grignani, T. Harmark, R. Oliveri, and M. Orselli (2021). Force-free magnetosphere attractors for near-horizon extreme and near-extreme limits of Kerr black hole. Classical and Quantum Gravity 38, 075022. doi:10.1088/1361-6382/abdf70, arXiv:2007.15662. https://arxiv.org/abs/2007.15662
6. R. A. Konoplya, J. Kunz, and A. Zhidenko (2021). Blandford-Znajek mechanism in the general stationary axially-symmetric black-hole spacetime. Journal of Cosmology and Astroparticle Physics 2021(12), 002. doi:10.1088/1475-7516/2021/12/002, arXiv:2102.10649. https://arxiv.org/abs/2102.10649
7. Y. Ma and L. Rezzolla (2024). Horizon-penetrating form of parametrized metrics for static and stationary black holes. Physical Review D 110, 024032. doi:10.1103/PhysRevD.110.024032, arXiv:2404.06509. https://arxiv.org/abs/2404.06509
8. C. Meringolo, F. Camilloni, and L. Rezzolla (2025). Electromagnetic Energy Extraction from Kerr Black Holes: Ab Initio Calculations. The Astrophysical Journal Letters 992, L8. doi:10.3847/2041-8213/ae06a6, arXiv:2507.08942. https://arxiv.org/abs/2507.08942
9. K. Chatterjee, Z. Younsi, P. Kocherlakota, and R. Narayan (2025). On the Universality of Energy Extraction from Black Hole Spacetimes. The Astrophysical Journal Letters 991, L58. doi:10.3847/2041-8213/ae0740, arXiv:2310.20043. https://arxiv.org/abs/2310.20043
10. J. Wu, L. Luo, Z. Lei, and X. Heng (2026). Model dependent analytic spin torsion corrections to Blandford Znajek energy extraction in Einstein Cartan gravity. arXiv:2608.23174. https://arxiv.org/abs/2608.23174
11. S. Murodov, O. Kholturayev, B. Rahmatov, J. Rayimbaev, I. Egamberdiev, and S. Karshiboev (2026). Blandford-Znajek Scaling in a Power-Law Rotating Kalb-Ramond Geometry: Magnetic-Flux Systematics and Bayesian Identifiability. arXiv:2608.11962. https://arxiv.org/abs/2608.11962
12. H. Feng, Z. Cai, R.-J. Yang, and J. Zhang (2026). On the perturbation solutions to the Blandford-Znajek mechanism in the Kerr-Sen black hole. Chinese Physics C 50, 065103. doi:10.1088/1674-1137/ae4327, arXiv:2603.03391. https://arxiv.org/abs/2603.03391
13. R. Yang, S. Chen, and J. Jing (2026). Effects of tidal charge on Blandford-Znajek process around braneworld black holes. The European Physical Journal C 86, 28. doi:10.1140/epjc/s10052-026-15288-w, arXiv:2510.24143. https://arxiv.org/abs/2510.24143