Dyadoregion in Kerr Black Hole Environments

• Dyadoregion is the region outside a black-hole event horizon where the proper electric field exceeds the Schwinger critical value, triggering electron-positron pair creation.
• It is defined invariantly in Kerr spacetimes using electromagnetic invariants instead of observer-dependent fields, ensuring a coordinate-independent characterization.
• Its morphology varies from axisymmetric polar lobes to non-axisymmetric multi-lobe structures, affecting beaming, energetics, and the initial conditions of plasma outflows.

Dyadoregion denotes the subset of spacetime outside a black-hole event horizon where the locally measured proper electric field reaches the Schwinger critical value and the QED vacuum breaks down into $e^\pm$ pairs. In Kerr spacetimes immersed in asymptotically uniform magnetic fields, the concept is formulated invariantly through the electromagnetic invariants rather than through any observer-dependent electric field, so the defining condition is $\tilde E \ge E_c$ in the local frame where the electric and magnetic fields are parallel. Recent analyses have treated both the aligned Wald configuration and the fully misaligned Bičák–Dvořák/Janiš solution, showing that frame dragging can induce overcritical electric fields without requiring the black hole to be charged and that the resulting dyadoregion can range from axisymmetric polar lobes to a non-axisymmetric multi-lobe structure with inclination-dependent energetics, beaming, and plasma initial conditions (Cherubini et al., 27 Oct 2025, Yang et al., 21 May 2026).

1. Invariant definition and field diagnostics

The physical definition is the same in both treatments: the dyadoregion is the spacetime region outside the event horizon where the locally measured electric field is large enough to trigger Schwinger pair creation. The critical field is

$$E_c=\frac{m_e^2 c^3}{e\hbar}\approx 1.32\times 10^{16}\ \mathrm{V\,cm^{-1}},$$

equivalently $B_c=E_c\approx 4.41\times 10^{13}\ \mathrm{G}$ in Gaussian units (Cherubini et al., 27 Oct 2025).

Its invariant characterization is based on

$$I_1 \equiv F_{\mu\nu}F^{\mu\nu}=2\,(B^2-E^2),\qquad I_2 \equiv F_{\mu\nu}{}^{\star}\!F^{\mu\nu}=-4\,\mathbf{E}\cdot\mathbf{B},$$

or, equivalently,

$$E_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}-I_1\big),\qquad B_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}+I_1\big).$$

In the notation of the cited works, these proper magnitudes are also written as $\tilde E$ and $\tilde B$. The dyadoregion is therefore the locus where $\tilde E(x)\ge E_c$, which is an observer-independent statement in curved spacetime (Cherubini et al., 27 Oct 2025, Yang et al., 21 May 2026).

For explicit local-field calculations, both papers use a ZAMO/LNRF tetrad. The ZAMO electric and magnetic fields are obtained from

$$E_\mu=F_{\mu\nu}Z^\nu,\qquad B_\mu=-{}^{\star}\!F_{\mu\nu}Z^\nu,$$

with physical components projected onto the orthonormal spatial basis. When the local fields are not parallel, one reaches the parallel-field frame by a local Lorentz boost along $\tilde E \ge E_c$0. In the aligned Wald configuration, the required boost is weakly relativistic in the astrophysically relevant domain and vanishes on the symmetry axis, where $\tilde E \ge E_c$1 (Cherubini et al., 27 Oct 2025).

This invariant construction is central because the electric and magnetic fields themselves are observer-dependent, whereas the dyadoregion criterion depends only on the invariants of $\tilde E \ge E_c$2. In that sense, the dyadoregion is not a coordinate artifact but a covariant characterization of overcritical electromagnetic structure near a rotating black hole.

2. Kerr background and external magnetic-field configurations

The background spacetime is Kerr in Boyer–Lindquist coordinates. In the misaligned treatment, the metric is written as

$\tilde E \ge E_c$3

with

$\tilde E \ge E_c$4

The event horizon is at

$\tilde E \ge E_c$5

and the horizon angular frequency is

$\tilde E \ge E_c$6

Both the dimensional spin $\tilde E \ge E_c$7 and the dimensionless spin $\tilde E \ge E_c$8 are used (Yang et al., 21 May 2026).

For an aligned asymptotically uniform magnetic field, the standard vacuum solution is the Wald solution with 4-potential

$\tilde E \ge E_c$9

where $$E_c=\frac{m_e^2 c^3}{e\hbar}\approx 1.32\times 10^{16}\ \mathrm{V\,cm^{-1}},$$0 and $$E_c=\frac{m_e^2 c^3}{e\hbar}\approx 1.32\times 10^{16}\ \mathrm{V\,cm^{-1}},$$1. This yields an induced electric field through frame dragging even for an uncharged Kerr black hole. The aligned analysis explicitly does not use the charged Wald variant with $$E_c=\frac{m_e^2 c^3}{e\hbar}\approx 1.32\times 10^{16}\ \mathrm{V\,cm^{-1}},$$2 (Cherubini et al., 27 Oct 2025).

For a magnetic field inclined by an angle $$E_c=\frac{m_e^2 c^3}{e\hbar}\approx 1.32\times 10^{16}\ \mathrm{V\,cm^{-1}},$$3 relative to the spin axis, the exact vacuum solution is the stationary, non-axisymmetric Bičák–Dvořák/Janiš solution. The construction used in the misaligned study sets

$$E_c=\frac{m_e^2 c^3}{e\hbar}\approx 1.32\times 10^{16}\ \mathrm{V\,cm^{-1}},$$4

and introduces the azimuthal variable

$$E_c=\frac{m_e^2 c^3}{e\hbar}\approx 1.32\times 10^{16}\ \mathrm{V\,cm^{-1}},$$5

which encodes the twist produced by frame dragging. For $$E_c=\frac{m_e^2 c^3}{e\hbar}\approx 1.32\times 10^{16}\ \mathrm{V\,cm^{-1}},$$6, axisymmetry is broken and the physical fields become functions of $$E_c=\frac{m_e^2 c^3}{e\hbar}\approx 1.32\times 10^{16}\ \mathrm{V\,cm^{-1}},$$7 rather than only $$E_c=\frac{m_e^2 c^3}{e\hbar}\approx 1.32\times 10^{16}\ \mathrm{V\,cm^{-1}},$$8 (Yang et al., 21 May 2026).

This distinction between aligned and misaligned fields is structurally decisive. In the aligned case, the problem is axisymmetric and mirror-symmetric about the equatorial plane. In the misaligned case, the non-axisymmetric vacuum solution introduces preferred azimuthal directions and fundamentally changes the topology of the overcritical region.

3. Geometry and topology of the dyadoregion

In the aligned Wald configuration, the dyadoregion boundary is the implicit surface

$$E_c=\frac{m_e^2 c^3}{e\hbar}\approx 1.32\times 10^{16}\ \mathrm{V\,cm^{-1}},$$9

The resulting morphology is not spherical. It forms two polar lobes above and below the equatorial plane, bounded below by the horizon. The electric field changes sign across “polar caps” with semi-aperture angle

$B_c=E_c\approx 4.41\times 10^{13}\ \mathrm{G}$0

beyond which field lines reverse. Along the axis, the fields are already parallel and the dyadoregion reaches its largest radius. In the weakly relativistic boost regime and at small polar angle,

$B_c=E_c\approx 4.41\times 10^{13}\ \mathrm{G}$1

so the approximate radius is

$B_c=E_c\approx 4.41\times 10^{13}\ \mathrm{G}$2

with $B_c=E_c\approx 4.41\times 10^{13}\ \mathrm{G}$3 on the axis (Cherubini et al., 27 Oct 2025).

In the misaligned configuration, the defining condition becomes

$B_c=E_c\approx 4.41\times 10^{13}\ \mathrm{G}$4

and the boundary is no longer axisymmetric. The large-$B_c=E_c\approx 4.41\times 10^{13}\ \mathrm{G}$5 approximation is

$B_c=E_c\approx 4.41\times 10^{13}\ \mathrm{G}$6

with

$B_c=E_c\approx 4.41\times 10^{13}\ \mathrm{G}$7

so that

$B_c=E_c\approx 4.41\times 10^{13}\ \mathrm{G}$8

For $B_c=E_c\approx 4.41\times 10^{13}\ \mathrm{G}$9,

$$I_1 \equiv F_{\mu\nu}F^{\mu\nu}=2\,(B^2-E^2),\qquad I_2 \equiv F_{\mu\nu}{}^{\star}\!F^{\mu\nu}=-4\,\mathbf{E}\cdot\mathbf{B},$$0

which has only discrete reflection symmetries with respect to the coordinate planes and is not axisymmetric (Yang et al., 21 May 2026).

Numerically, the misaligned dyadoregion consists of several discrete lobes whose number, size, and orientation vary with inclination. For representative parameters $$I_1 \equiv F_{\mu\nu}F^{\mu\nu}=2\,(B^2-E^2),\qquad I_2 \equiv F_{\mu\nu}{}^{\star}\!F^{\mu\nu}=-4\,\mathbf{E}\cdot\mathbf{B},$$1 and $$I_1 \equiv F_{\mu\nu}F^{\mu\nu}=2\,(B^2-E^2),\qquad I_2 \equiv F_{\mu\nu}{}^{\star}\!F^{\mu\nu}=-4\,\mathbf{E}\cdot\mathbf{B},$$2, the $$I_1 \equiv F_{\mu\nu}F^{\mu\nu}=2\,(B^2-E^2),\qquad I_2 \equiv F_{\mu\nu}{}^{\star}\!F^{\mu\nu}=-4\,\mathbf{E}\cdot\mathbf{B},$$3–$$I_1 \equiv F_{\mu\nu}F^{\mu\nu}=2\,(B^2-E^2),\qquad I_2 \equiv F_{\mu\nu}{}^{\star}\!F^{\mu\nu}=-4\,\mathbf{E}\cdot\mathbf{B},$$4 sections show four-to-six lobes with a periodicity of $$I_1 \equiv F_{\mu\nu}F^{\mu\nu}=2\,(B^2-E^2),\qquad I_2 \equiv F_{\mu\nu}{}^{\star}\!F^{\mu\nu}=-4\,\mathbf{E}\cdot\mathbf{B},$$5 and reflection symmetry under $$I_1 \equiv F_{\mu\nu}F^{\mu\nu}=2\,(B^2-E^2),\qquad I_2 \equiv F_{\mu\nu}{}^{\star}\!F^{\mu\nu}=-4\,\mathbf{E}\cdot\mathbf{B},$$6. For $$I_1 \equiv F_{\mu\nu}F^{\mu\nu}=2\,(B^2-E^2),\qquad I_2 \equiv F_{\mu\nu}{}^{\star}\!F^{\mu\nu}=-4\,\mathbf{E}\cdot\mathbf{B},$$7, two polar lobes dominate and four smaller symmetric off-axis lobes are present. As $$I_1 \equiv F_{\mu\nu}F^{\mu\nu}=2\,(B^2-E^2),\qquad I_2 \equiv F_{\mu\nu}{}^{\star}\!F^{\mu\nu}=-4\,\mathbf{E}\cdot\mathbf{B},$$8 increases toward $$I_1 \equiv F_{\mu\nu}F^{\mu\nu}=2\,(B^2-E^2),\qquad I_2 \equiv F_{\mu\nu}{}^{\star}\!F^{\mu\nu}=-4\,\mathbf{E}\cdot\mathbf{B},$$9, the polar lobes shrink, off-axis lobes rotate, and $$E_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}-I_1\big),\qquad B_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}+I_1\big).$$0-pairs disappear and reappear at critical angles $$E_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}-I_1\big),\qquad B_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}+I_1\big).$$1 and $$E_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}-I_1\big),\qquad B_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}+I_1\big).$$2; eventually a different pair $$E_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}-I_1\big),\qquad B_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}+I_1\big).$$3 becomes dominant. For $$E_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}-I_1\big),\qquad B_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}+I_1\big).$$4, the evolution mirrors the first half-period (Yang et al., 21 May 2026).

This multi-lobe morphology replaces the “dyadotorus” known for certain charged, axisymmetric Kerr–Newman configurations and generalizes the two polar lobes, or “twin cones,” of the aligned Wald case. A plausible implication is that the angular structure of any resulting outflow should inherit this inclination-dependent lobe geometry.

4. Electromagnetic energetics and beaming

The intrinsic electromagnetic energy available for pair creation is computed by integrating the Killing energy density over the dyadoregion. In the misaligned treatment,

$$E_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}-I_1\big),\qquad B_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}+I_1\big).$$5

which reduces, in ZAMO variables, to

$$E_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}-I_1\big),\qquad B_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}+I_1\big).$$6

with

$$E_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}-I_1\big),\qquad B_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}+I_1\big).$$7

In the outer dyadoregion, $$E_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}-I_1\big),\qquad B_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}+I_1\big).$$8, so the energy is magnetically dominated. The large-$$E_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}-I_1\big),\qquad B_0^2=\tfrac{1}{2}\big(\sqrt{I_1^2+I_2^2}+I_1\big).$$9 scaling is

$\tilde E$0

where $\tilde E$1 and

$\tilde E$2

For the aligned case $\tilde E$3, $\tilde E$4 (Yang et al., 21 May 2026).

The aligned Wald analysis obtains the analogous estimate

$\tilde E$5

so at fixed spin the energy scales as $\tilde E$6. For $\tilde E$7, $\tilde E$8, and $\tilde E$9, this gives $\tilde B$0 erg (Cherubini et al., 27 Oct 2025). In the misaligned example $\tilde B$1, $\tilde B$2, $\tilde B$3, $\tilde B$4, the full integral gives $\tilde B$5, while the analytic scaling yields $\tilde B$6 (Yang et al., 21 May 2026).

Because the energy is angularly concentrated, the isotropic-equivalent energy exceeds the intrinsic dyadoregion energy. In the aligned Wald model, a simple double-cone estimate gives

$\tilde B$7

so $\tilde B$8 (Cherubini et al., 27 Oct 2025).

The misaligned study instead defines an energy-containment beaming factor

$\tilde B$9

where $\tilde E(x)\ge E_c$0 is the minimum solid angle enclosing $\tilde E(x)\ge E_c$1 of the dyadoregion energy. The resulting $\tilde E(x)\ge E_c$2 is symmetric and bell-shaped, with a minimum $\tilde E(x)\ge E_c$3 at $\tilde E(x)\ge E_c$4 and a maximum $\tilde E(x)\ge E_c$5 at $\tilde E(x)\ge E_c$6. The rise has inflection points near $\tilde E(x)\ge E_c$7, marking changes in the number of dominant lobes. For $\tilde E(x)\ge E_c$8, $\tilde E(x)\ge E_c$9, $$E_\mu=F_{\mu\nu}Z^\nu,\qquad B_\mu=-{}^{\star}\!F_{\mu\nu}Z^\nu,$$0, $$E_\mu=F_{\mu\nu}Z^\nu,\qquad B_\mu=-{}^{\star}\!F_{\mu\nu}Z^\nu,$$1, the corresponding isotropic-equivalent energy is $$E_\mu=F_{\mu\nu}Z^\nu,\qquad B_\mu=-{}^{\star}\!F_{\mu\nu}Z^\nu,$$2 (Yang et al., 21 May 2026).

The aligned and misaligned beaming prescriptions are therefore not numerically identical. The later study explicitly notes that its aligned $$E_\mu=F_{\mu\nu}Z^\nu,\qquad B_\mu=-{}^{\star}\!F_{\mu\nu}Z^\nu,$$3 is larger than the twin-cone estimate $$E_\mu=F_{\mu\nu}Z^\nu,\qquad B_\mu=-{}^{\star}\!F_{\mu\nu}Z^\nu,$$4 because it is defined by energy containment rather than by an assumed cone geometry. This is a difference of definition rather than a contradiction in the underlying energetics.

5. Pair creation, local rates, and $$E_\mu=F_{\mu\nu}Z^\nu,\qquad B_\mu=-{}^{\star}\!F_{\mu\nu}Z^\nu,$$5 plasma thermodynamics

In the local parallel-field frame, the Schwinger pair-creation rate per unit 4-volume is written in the misaligned analysis as

$$E_\mu=F_{\mu\nu}Z^\nu,\qquad B_\mu=-{}^{\star}\!F_{\mu\nu}Z^\nu,$$6

and in the pure-electric limit $$E_\mu=F_{\mu\nu}Z^\nu,\qquad B_\mu=-{}^{\star}\!F_{\mu\nu}Z^\nu,$$7 as

$$E_\mu=F_{\mu\nu}Z^\nu,\qquad B_\mu=-{}^{\star}\!F_{\mu\nu}Z^\nu,$$8

The aligned study presents the equivalent invariant Schwinger expression in terms of $$E_\mu=F_{\mu\nu}Z^\nu,\qquad B_\mu=-{}^{\star}\!F_{\mu\nu}Z^\nu,$$9 and $\tilde E \ge E_c$00, emphasizing that in practice the $\tilde E \ge E_c$01 term is sufficient because of the exponential suppression (Yang et al., 21 May 2026, Cherubini et al., 27 Oct 2025).

Both treatments apply these flat-space QED rates locally in curved spacetime using local tetrads and the invariant volume element. The aligned work states that the locally constant field approximation is justified because the field varies on macroscopic scales relative to the Compton wavelength, while the misaligned work notes that curvature enters through the local tetrad construction and $\tilde E \ge E_c$02 (Cherubini et al., 27 Oct 2025, Yang et al., 21 May 2026).

For plasma thermodynamics, the misaligned study assumes rapid thermalization into an $\tilde E \ge E_c$03 plasma on the Compton timescale

$\tilde E \ge E_c$04

and obtains

$\tilde E \ge E_c$05

leading to

$\tilde E \ge E_c$06

On the dyadoregion boundary, where $\tilde E \ge E_c$07, the temperature is nearly constant: $\tilde E \ge E_c$08 For $\tilde E \ge E_c$09, $\tilde E \ge E_c$10 (Yang et al., 21 May 2026).

The aligned Wald analysis describes pair-photon thermalization on timescales $\tilde E \ge E_c$11 s and derives

$\tilde E \ge E_c$12

For representative parameters, $\tilde E \ge E_c$13 at the boundary, while near the pole on the horizon $\tilde E \ge E_c$14 can reach a few $\tilde E \ge E_c$15, for example $\tilde E \ge E_c$16 for $\tilde E \ge E_c$17 (Cherubini et al., 27 Oct 2025). The numerical values are therefore configuration-dependent.

In both analyses, the equilibrium equation of state is

$\tilde E \ge E_c$18

with magnetic pressure $\tilde E \ge E_c$19 or $\tilde E \ge E_c$20. A central conclusion is that $\tilde E \ge E_c$21, including near the horizon for typical GRB parameters, so the nascent plasma is initially magnetically dominated (Yang et al., 21 May 2026, Cherubini et al., 27 Oct 2025).

This magnetic dominance is physically significant because it supports Poynting-flux–driven acceleration and ultrarelativistic outflows. A plausible implication is that the dyadoregion should be interpreted not only as a pair-production volume but also as a set of initial conditions for subsequent GRMHD evolution.

6. Threshold field strength, assumptions, and astrophysical interpretation

For fixed spin and inclination, vacuum breakdown requires that the peak proper electric field reach $\tilde E \ge E_c$22 outside the horizon. In the aligned Wald treatment, the minimum field is

$\tilde E \ge E_c$23

For $\tilde E \ge E_c$24, this gives $\tilde E \ge E_c$25 (Cherubini et al., 27 Oct 2025).

The broader misaligned analysis gives, for the aligned case $\tilde E \ge E_c$26,

$\tilde E \ge E_c$27

An extremal Kerr black hole with $\tilde E \ge E_c$28 therefore has $\tilde E \ge E_c$29, corresponding to

$\tilde E \ge E_c$30

For misalignment, $\tilde E \ge E_c$31 is obtained numerically and is lower than the aligned value for most inclinations, especially near $\tilde E \ge E_c$32. For $\tilde E \ge E_c$33, however, there exist four angles per $\tilde E \ge E_c$34 where breakdown is slightly harder; at $\tilde E \ge E_c$35 and $\tilde E \ge E_c$36, one finds $\tilde E \ge E_c$37, or $\tilde E \ge E_c$38 (Yang et al., 21 May 2026).

This result clarifies an important point in the literature. The aligned Wald study notes that detailed misaligned configurations are beyond its analytic model and states that misalignment would reduce the peak $\tilde E \ge E_c$39 and shrink or suppress the dyadoregion. The later explicit calculation with the full non-axisymmetric vacuum solution instead finds that misaligned magnetic fields generally favor pair creation more than aligned ones. The contrast reflects the difference between an extrapolation from the aligned model and a direct computation in the misaligned geometry (Cherubini et al., 27 Oct 2025, Yang et al., 21 May 2026).

Both analyses adopt the test-field approximation,

$\tilde E \ge E_c$40

so electromagnetic backreaction on the Kerr metric is neglected. They also assume vacuum Maxwell fields before breakdown, neglect plasma currents, resistive effects, and screening during the initial phase, and apply local flat-space Schwinger formulas in tetrads under the assumption that curvature is negligible on microphysical scales. The misaligned treatment further notes that radiative transfer and photon escape are not modeled (Yang et al., 21 May 2026, Cherubini et al., 27 Oct 2025).

Astrophysically, the parameter range emphasized is that of stellar-mass black holes with $\tilde E \ge E_c$41–$\tilde E \ge E_c$42, high spins, and magnetic fields $\tilde E \ge E_c$43–$\tilde E \ge E_c$44. In this regime, the misaligned model predicts dyadoregion energies $\tilde E \ge E_c$45–$\tilde E \ge E_c$46, with isotropic-equivalent energies $\tilde E \ge E_c$47–$\tilde E \ge E_c$48 after dividing by $\tilde E \ge E_c$49, consistent with GRB prompt energetics (Yang et al., 21 May 2026). The aligned paper similarly finds that stellar-mass systems can comfortably exceed threshold, while supermassive black holes with AGN-like fields, such as $\tilde E \ge E_c$50 and $\tilde E \ge E_c$51, have $\tilde E \ge E_c$52 and do not form a dyadoregion (Cherubini et al., 27 Oct 2025).

Relative to earlier dyadosphere and dyadotorus concepts, the dyadoregion around a Kerr black hole in an external magnetic field establishes that vacuum breakdown need not rely on a black-hole monopolar charge. In the aligned case the induced field is globally quadrupolar and produces polar lobes; in the misaligned case axial symmetry is broken and the overcritical region fragments into several discrete lobes. This geometry, together with the initial magnetic dominance of the $\tilde E \ge E_c$53 plasma, supports ultrarelativistic outflows whose angular structure can reflect the lobe pattern and potentially imprint inclination-dependent polarization and jet morphology signatures (Cherubini et al., 27 Oct 2025, Yang et al., 21 May 2026).