Vacuum breakdown in a misaligned magnetized Kerr spacetime

Abstract: Electron-positron ($e^{+}e^{-}$) pair creation by vacuum breakdown around compact objects is believed to power high-energy astrophysical transients like gamma-ray bursts (GRBs). In this work, we focus on vacuum breakdown around a Kerr black hole (BH) immersed in an asymptotically uniform magnetic field that is inclined with respect to the BH spin axis. The dyadoregion, the region where the induced electric field exceeds the critical value $E_{\text{c}}=m_{e}^{2}c^{3}/(e\hbar)$, is identified via the electromagnetic invariants. It is found that the dyadoregion consists of several lobes whose number, size, and orientation vary with the inclination. We also estimate the electromagnetic energy available for pair creation and derive a beaming factor that allows a conversion between the intrinsic dyadoregion energy and the observed isotropic energy. The thermodynamic properties of the resulting electron-positron-photon ($e^{+}e^{-}γ$) plasma are included, revealing an initial magnetic dominance. The evaluation of the minimum magnetic field required shows that misaligned magnetic fields generally favor pair creation more than aligned ones.

• The paper presents an analytical and numerical investigation of electron-positron pair creation via vacuum breakdown near Kerr black holes in misaligned magnetic fields, establishing critical conditions for the process.
• It details the dyadoregion morphology, showing how misalignment induces complex multi-lobed structures compared to the axisymmetric case, with quantifiable effects on energy distribution and beaming.
• Results indicate that misalignment lowers the threshold magnetic field for breakdown, providing insights into gamma-ray burst energetics and the dynamics of relativistic plasma outflows.

Vacuum Breakdown and Dyadoregion Structure in Misaligned Magnetized Kerr Spacetimes

Introduction

The study investigates electron-positron ($e^+e^-$) pair creation via vacuum breakdown in the vicinity of Kerr black holes immersed in asymptotically uniform magnetic fields misaligned with the black hole (BH) spin axis (2605.21910). This mechanism, relevant for astrophysical transients such as gamma-ray bursts (GRBs), contrasts with neutrino annihilation models, offering rapid and efficient conversion of electromagnetic energy into a hot, dense $e^+e^-\gamma$ plasma. The focus is a systematic analysis of the "dyadoregion"—the spatial domain where the electric field surpasses the critical Schwinger threshold $E_c = m_e^2 c^3 / (e\hbar)$—for misaligned field configurations, along with associated energetics and plasma thermodynamics.

Electromagnetic Field Configuration in Misaligned Kerr Spacetimes

The electromagnetic field solution, extending Bičák and Dvořák's construction, describes a Kerr BH within a uniform but arbitrarily inclined magnetic field. Unlike the axisymmetric Wald solution, the misalignment (inclination $i$ between field and spin axis) breaks symmetry, yielding non-axisymmetric Maxwell fields. The asymptotic field components $B_0$, $B_1$ encode the parallel and perpendicular projections relative to the spin axis, with $B = \sqrt{B_0^2 + B_1^2}$ and $B_0 = B\cos i$, $B_1 = B\sin i$. For observers, the electromagnetic field is characterized using the ZAMO (zero angular momentum observer) tetrad, determining the electric and magnetic field projections relevant for local pair production analysis.

The spatial structure of electric fields for various inclination angles, as measured by ZAMOs, is illustrated in (Figure 1).

Figure 1: Schematic representation of the electric fields measured by ZAMOs for selected inclinations $i$, with the field strength increasing toward the horizon. Axisymmetry is recovered only at $e^+e^-\gamma$0.

Dyadoregion Identification and Morphology

The dyadoregion's spatial extent, where vacuum breakdown proceeds, is set by the electric field exceeding $e^+e^-\gamma$1 in the frame where the local electric and magnetic fields are parallel. This is computed from electromagnetic invariants, yielding the effective field

$e^+e^-\gamma$2

with $e^+e^-\gamma$3 and $e^+e^-\gamma$4 as the two Lorentz invariants.

In aligned fields ($e^+e^-\gamma$5), the dyadoregion is axisymmetric, with two dominant polar lobes. For misalignment, the region's morphology becomes complex: several lobes emerge, whose number, size, and orientation evolve nontrivially with $e^+e^-\gamma$6. These features are visualized for a range of inclinations in (Figure 2).

Figure 2: The dyadoregion in the $e^+e^-\gamma$7-$e^+e^-\gamma$8 plane for different inclination angles $e^+e^-\gamma$9, showing lobe number and size variation as $E_c = m_e^2 c^3 / (e\hbar)$0 increases from $E_c = m_e^2 c^3 / (e\hbar)$1. For $E_c = m_e^2 c^3 / (e\hbar)$2, $E_c = m_e^2 c^3 / (e\hbar)$3, the pattern is axisymmetric only for $E_c = m_e^2 c^3 / (e\hbar)$4.

This non-axisymmetric dyadoregion reflects the underlying lack of symmetry in the electromagnetic field, with lobe transitions and disappearances at critical inclination values. An analytic large-distance approximation provides a means to quickly estimate dyadoregion structure as a function of parameters $E_c = m_e^2 c^3 / (e\hbar)$5.

Minimum Magnetic Field for Vacuum Breakdown

Not all field configurations drive vacuum breakdown outside the horizon. For each spin $E_c = m_e^2 c^3 / (e\hbar)$6 and inclination $E_c = m_e^2 c^3 / (e\hbar)$7, there is a minimum field strength $E_c = m_e^2 c^3 / (e\hbar)$8 such that the maximum induced $E_c = m_e^2 c^3 / (e\hbar)$9 exceeds $i$0 somewhere outside the BH. In the aligned case (analytically tractable), $i$1 increases at high spin but is otherwise determined by geometry. When the field is misaligned, $i$2 generically decreases for most values of $i$3, indicating that misalignment favors vacuum breakdown.

This dependence is shown in (Figure 3).

Figure 3: The minimum magnetic field $i$4 as a function of inclination $i$5, for various BH spins $i$6. Vacuum breakdown is more readily triggered for a wide range of misalignments, especially near $i$7.

A noteworthy exception is that, for extremal Kerr BHs, even perfectly aligned fields can permit vacuum breakdown, with a threshold field $i$8 (in units of $i$9), unaffected by the classic Meissner expulsion effect for axisymmetric stationary fields.

Dyadoregion Energetics and Beaming

The electromagnetic energy content in the dyadoregion, $B_0$0, is computed through the conserved Killing integral restricted to the region where $B_0$1. Numerically, $B_0$2 grows monotonically with both spin $B_0$3 and field strength $B_0$4. As a function of inclination, the total energy displays a $B_0$5-periodic cosine-like oscillation, with maximum for $B_0$6 (alignment) and minimum near $B_0$7 (perpendicular).

This parameter dependence is depicted in (Figure 4).

Figure 4: Left: Dyadoregion electromagnetic energy $B_0$8 vs. magnetic field strength $B_0$9 for $B_1$0; Right: $B_1$1 vs. inclination $B_1$2 for $B_1$3, with different curves for spins $B_1$4. Dashed lines are analytic approximations valid for strong fields or high spin.

A scaling law

$B_1$5

is established, with $B_1$6 encoding the inclination dependence; its variation is much less than that due to changes in $B_1$7 or $B_1$8.

Astrophysically, observed energies are often inferred assuming isotropic emission. The actual beaming can be quantified by calculating the solid angle enclosing, e.g., 90% of the total dyadoregion energy. The "beaming factor" $B_1$9, the ratio of this energy region’s solid angle to $B = \sqrt{B_0^2 + B_1^2}$0, is inclination dependent and generally increases from $B = \sqrt{B_0^2 + B_1^2}$1 at $B = \sqrt{B_0^2 + B_1^2}$2 to $B = \sqrt{B_0^2 + B_1^2}$3 at $B = \sqrt{B_0^2 + B_1^2}$4, as shown in (Figure 5).

Figure 5: Plot of the beaming factor $B = \sqrt{B_0^2 + B_1^2}$5 vs. inclination $B = \sqrt{B_0^2 + B_1^2}$6. The beaming is most concentrated (smallest $B = \sqrt{B_0^2 + B_1^2}$7) for aligned fields, and most isotropic for perpendicular configurations.

Consequently, misaligned configurations yield more isotropic energy release, directly impacting observable energetics in GRBs and similar phenomena.

Plasma Thermodynamics in the Dyadoregion

Pair creation results in rapid formation of an $B = \sqrt{B_0^2 + B_1^2}$8 plasma, assumed to quickly thermalize. The temperature and density are obtained by equating the local creation rate (from the Schwinger process in curved spacetime with magnetic fields) to the equilibrium Fermi-Dirac density for $B = \sqrt{B_0^2 + B_1^2}$9 at zero chemical potential, leading to

$B_0 = B\cos i$0

For typical dyadoregion parameters ($B_0 = B\cos i$1), the temperature on the boundary is $B_0 = B\cos i$2.

The spatial temperature and pressure profiles closely mirror the dyadoregion shape, with hot, dense plasma predominantly located within the lobe structures. (Figure 6) and (Figure 7) provide further details on temperature and pressure distributions on the horizon and as a function of inclination.

Figure 6: Left panel: Normalized plasma temperature $B_0 = B\cos i$3 in the $B_0 = B\cos i$4-$B_0 = B\cos i$5 plane ($B_0 = B\cos i$6, $B_0 = B\cos i$7, $B_0 = B\cos i$8). Right: Plasma parameter $B_0 = B\cos i$9, demonstrating magnetic domination in the dyadoregion.

Figure 7: Left: North-pole horizon temperature $B_1 = B\sin i$0 vs. spin $B_1 = B\sin i$1 for selected $B_1 = B\sin i$2 at $B_1 = B\sin i$3. Right: $B_1 = B\sin i$4 vs. inclination $B_1 = B\sin i$5 for $B_1 = B\sin i$6.

The ratio $B_1 = B\sin i$7 is everywhere much less than unity (Figure 8), implying dynamics are initially magnetically dominated, favorable for launching highly relativistic outflows through magnetohydrodynamic mechanisms—one of the requirements for GRB engines.

Figure 8: Ratio $B_1 = B\sin i$8 at the north pole of the horizon across $B_1 = B\sin i$9 and $i$0. Behavior tracks that of temperature, confirming low plasma beta in the dyadoregion.

Implications and Prospects

The analysis demonstrates that realistic misaligned magnetized Kerr spacetimes generically enable vacuum breakdown at lower critical field strengths than aligned cases, increasing the likelihood of the process in astrophysical settings with nontrivial magnetic topology. The morphology and energetics of the dyadoregion exhibit rich $i$1- and $i$2-dependences, with observable consequences for the angular distribution and apparent luminosity of prompt emissions in transients. The strong initial magnetic dominance supports the development of magnetically driven ultra-relativistic plasma outflows.

Numerical results indicate:

• For stellar-mass black holes ($i$3, $i$4, $i$5), the intrinsic energy in the dyadoregion can exceed $i$6, matching observed GRB energetics after incorporating beaming corrections.
• The dyadoregion is most spatially concentrated when the field and spin axis are aligned, but practical systems are likely to favor mild to moderate misalignments, producing more isotropic emission and lowering pair creation thresholds.

Future theoretical developments could address the fully coupled Einstein-Maxwell system without test field approximations, as in the Kerr–Bertotti–Robinson spacetime, or pursue full GRMHD simulations to model the post-thermalization evolution and electromagnetic breakout.

Conclusion

This work extends previous analyses of vacuum breakdown around Kerr black holes to incorporate arbitrary magnetic field orientations. The key findings are that misalignment between the external field and spin axis generally facilitates vacuum breakdown and results in complex, multi-lobed dyadoregion morphologies, whose structure, energetics, and beaming are crucial for understanding the origins and observability of $i$7 plasma-driven astrophysical transients. Quantitative results connect the theoretical underpinning to plausible GRB scenarios, and the analytic and numerical tools developed enable robust modeling of such environments, setting the stage for refined simulation and observation-based inference in future high-energy astrophysics.