⚛️ general relativity

Charging of rotating black holes: kinetic simulations of black hole magnetospheres

Kinetic general-relativistic particle-in-cell simulations demonstrate that rotating black holes immersed in external magnetic fields dynamically relax to an equilibrium charge significantly lower than the theoretical Wald limit, driven by plasma processes rather than initial conditions, thereby preserving the potential drops necessary for energy-extraction mechanisms like the Blandford-Znajek process.

Original authors: Martin Kološ, Farukh Abdulkhamidov, Arman Tursunov, Jorge A. Rueda, Benoît Cerutti

Published 2026-09-16

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Martin Kološ, Farukh Abdulkhamidov, Arman Tursunov, Jorge A. Rueda, Benoît Cerutti

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Black holes are often imagined as the universe's ultimate vacuum cleaners, swallowing everything that comes too close. But when a black hole spins and sits inside a magnetic field, it becomes something far more complex: a cosmic engine. In the space around these spinning giants, magnetic field lines get twisted by the rotation, creating powerful electric fields. These fields act like a giant separator, pulling positive charges one way and negative charges the other. For decades, physicists have debated what happens next. Does the black hole simply keep absorbing charges until it reaches a specific, predictable limit where the electric forces balance out? Or does the chaotic, living plasma surrounding the black hole change the rules entirely? This question matters because the amount of electric charge a black hole holds influences how it powers the massive jets of energy that shoot out from its poles, which in turn affects how galaxies evolve.

A team of researchers has now settled this debate by running the most detailed simulations of their kind. They used a supercomputer to model the behavior of a spinning black hole surrounded by a magnetic field, but with a crucial twist: they did not just assume the space around the hole was empty or filled with a simple, smooth fluid. Instead, they simulated the actual particles—electrons and their antimatter twins, positrons—that fill the space. These particles are created in pairs by high-energy radiation, a process that happens naturally near black holes. The researchers started their simulations in two different ways. In one set of runs, they began with a black hole that had no electric charge at all. In the other set, they started with a black hole that already held the maximum theoretical charge predicted by older, simpler theories. They then watched to see what happened over time.

The results were striking. Regardless of whether the black hole started out neutral or heavily charged, the system quickly forgot its beginning. Within a short period, both scenarios evolved into the exact same state. The black hole settled into a stable, positive electric charge that was much lower than the maximum limit previously thought to be the natural endpoint. The researchers found that for most spinning black holes, the final charge is only about thirty percent of the theoretical maximum. As the black hole's spin gets closer to the fastest possible speed, this charge drops even further, eventually fading to almost nothing. This finding proves that the charge of a black hole is not a fixed value determined by the initial conditions or a simple vacuum formula. Instead, it is a dynamic quantity, constantly regulated by the flow of particles and the creation of new matter in the surrounding plasma.

The study also clarified why the black hole does not reach that theoretical maximum. The electric field around the hole acts like a valve, directing positive charges toward the poles and negative charges toward the equator. The system reaches a balance when the magnetic field lines carrying these charges are split evenly between the two regions. The researchers derived a simple rule based on this balance, which matched their complex computer simulations almost perfectly. This confirms that the charge is determined by how the magnetic field threads through the black hole's surface, rather than by the strength of the electric field alone.

Perhaps most importantly, the fact that the black hole settles at a lower charge than the maximum means that the powerful electric fields driving energy extraction are never completely switched off. Even in a state of equilibrium, there remains a voltage difference between the black hole and the distant universe. This ensures that the mechanisms that power the high-speed jets of energy, such as the Blandford-Znajek process, continue to function. The black hole remains a potent engine, capable of accelerating particles and launching jets, because the kinetic dance of the plasma prevents it from ever becoming fully neutralized or fully saturated. The work provides a new, concrete understanding of how these cosmic engines are charged, replacing old theoretical limits with a picture driven by the real, messy physics of particle creation and magnetic fields.

Technical Summary: Charging of Rotating Black Holes via Kinetic Simulations

Problem Statement
The paper addresses a long-standing open question in black hole (BH) electrodynamics: whether a rotating black hole immersed in an external magnetic field charges up to the Wald equilibrium value, $Q_W = 2aMB$ (where $M$ is mass, $a$ is the spin parameter, and $B$ is the magnetic field strength). While the Wald solution suggests $Q_W$ is the equilibrium state where the electric potential difference between the horizon and infinity vanishes, previous theoretical approaches have yielded conflicting expectations. Some models suggest a sub-Wald equilibrium based on energy minimization or selective accretion, while others argue for complete screening ( $Q \approx Q_W$ ) or oscillatory behavior. Crucially, existing force-free or ideal magnetohydrodynamic (MHD) models cannot address this problem because they enforce $\vec{E} \cdot \vec{B} = 0$ , thereby removing the parallel electric fields necessary for charge separation. The central question is whether an astrophysical BH, embedded in a plasma-filled magnetosphere with self-consistent pair creation, evolves toward $Q_W$ or a different value, and whether this final state depends on initial conditions.

Methodology
The authors employ first-principles, axisymmetric general-relativistic particle-in-cell (GRPIC) simulations using the GRZeltron code. This approach avoids assumptions about plasma conductivity or supply mechanisms.

Setup: Simulations evolve an asymptotically uniform Wald magnetic field around a Kerr black hole for various spin parameters ( $a = 0.2$ to $0.95$).
Initial Conditions: Two distinct initial states are tested for each spin: a neutral black hole ( $Q_0 = 0$ ) and a black hole pre-charged to the Wald value ( $Q_0 = Q_W$ ).
Physics Included: The simulations self-consistently model pair creation via inverse-Compton scattering of soft photons, radiation reaction forces, and the full kinetic response of the plasma.
Diagnostics: The accumulated charge is measured via Gauss's law flux through the event horizon and cross-verified by integrating the radial current flux. The system is evolved for $\approx 160$ gravitational times ( $t_g$ ) to reach a saturated state.

Key Contributions and Analytic Framework
The paper provides the first direct GRPIC determination of the equilibrium BH charge in the Wald problem. Additionally, the authors derive an analytic expression for the equilibrium charge ( $\xi_{eq} \equiv Q_{eq}/Q_W$ ) based on a magnetic flux balance criterion.

Analytic Derivation: The authors argue that equilibrium is reached when the magnetic flux threading the polar cap (where positive charges are captured) equals the flux threading the equatorial belt (where negative charges are captured). This requires the separatrix defined by $\vec{E} \cdot \vec{B} = 0$ to bisect the magnetic flux of the horizon.
Resulting Formula: This yields a closed-form expression for $\xi_{eq}$ (Eq. 22) that depends on the spin $a$ . The model predicts $\xi_{eq} \approx 0.25$ in the slow-rotation limit and a decline toward zero as $a \to 1$ due to the expulsion of magnetic flux from the horizon (the Meissner effect).

Results

Dynamical Attractor: The simulations demonstrate that the system loses memory of its initial charge state. Both the $Q_0 = 0$ and $Q_0 = Q_W$ branches converge to the same positive, saturated charge value within a few tens of gravitational times. This confirms the equilibrium charge is a dynamical attractor of the kinetic magnetosphere.
Sub-Wald Equilibrium: The equilibrium charge is consistently well below the Wald limit.
- For $a \lesssim 0.7$ , the normalized charge is $\xi_{eq} \approx 0.30\text{--}0.32$ .
- For higher spins ( $a = 0.9$ ), $\xi_{eq} \approx 0.20$ .
- For near-extreme spins ( $a = 0.95$ ), $\xi_{eq} \approx 0.14$ .
- In no case does the charge approach $\xi = 1$ or the much smaller minimum-energy values ( $\xi \lesssim 0.05$ ) proposed in other literature.
Comparison with Analytic Model: The measured values agree qualitatively with the analytic flux-balance prediction, reproducing the flat plateau at low spins and the decline at high spins. However, the simulated charges are systematically 20–30% higher than the analytic prediction at intermediate spins. The authors attribute this discrepancy to plasma currents suppressing the Meissner expulsion of magnetic flux, keeping the polar capture channel open longer than the vacuum Wald solution predicts.

Significance and Conclusions
The paper concludes that the charge state of an astrophysical black hole is driven by kinetic plasma processes rather than vacuum electrostatics or initial data.

Upper Bound: The Wald charge $Q_W$ acts as an upper bound, not a general equilibrium value.
Implications for Energy Extraction: Since $\xi_{eq} < 1$ , the potential drop between the horizon and infinity is not fully screened. Consequently, the accelerating gaps near the inner light surface persist, pair discharges continue to load the magnetosphere, and the Blandford-Znajek mechanism for energy extraction is not quenched by charging.
Limitations: The authors note that their results are restricted to axisymmetric configurations, excluding non-axisymmetric instabilities of the current sheet and tilted magnetic field geometries. They also acknowledge that the scale separation between the gravitational radius and the Larmor radius in simulations is smaller than in real astrophysical sources.

The work establishes that the equilibrium charge is a kinetic quantity determined by the competition between horizon absorption of opposite charge signs in a pair-producing plasma, providing a new benchmark for models of black hole magnetospheres and ultra-high-energy cosmic ray acceleration.

Technical Summary: Charging of Rotating Black Holes via Kinetic Simulations

More like this