The baryon content of magnetically arrested black hole disks and jets

Using general relativistic magnetohydrodynamic simulations with a passive tracer to isolate physical baryon density, this study reveals that baryon loading in magnetically arrested black hole jets is an episodic process driven by magnetic flux eruptions and shear-driven waves, leading to charge-starved, baryon-poor plasma in spinning systems with implications for particle acceleration and non-thermal emission.

The Gist

Imagine a supermassive black hole not as a cosmic vacuum cleaner, but as a chaotic, high-speed kitchen where a chef is trying to make a very specific, high-energy smoothie (a relativistic jet) to shoot out into space.

The problem? The chef needs to know exactly what ingredients are going into the smoothie. Is it mostly pure energy (photons and electrons/positrons), or is it "heavy stuff" like protons and ions (baryons)?

This paper is about a new way to track those ingredients in a computer simulation, revealing that the black hole's "smoothie" is actually much emptier of heavy stuff than we thought, and that the process is messy, rhythmic, and full of explosions.

Here is the breakdown using everyday analogies:

1. The Problem: The "Fake Ingredient" Glitch

In previous computer simulations of black holes, scientists had a problem. To keep the computer from crashing when the area around the black hole became too empty (too low density), they had to secretly inject "dummy" mass—like adding water to a soup just so the pot doesn't boil dry.

The issue? This "dummy water" mixed with the real ingredients (the heavy protons from the accretion disk). By the time the soup was ready to shoot out as a jet, scientists couldn't tell how much was real heavy matter and how much was just the fake water added to keep the simulation stable. It was like trying to taste a cake but not knowing if the sugar was real or just a sweetener added to make the batter look right.

The Solution: The authors invented a "Passive Tracer."
Think of this as a glow-in-the-dark dye.

• They only put the dye in the real ingredients (the disk) at the start.
• When the computer adds "dummy water" (fake mass) to keep things stable, it doesn't get the dye.
• Now, the scientists can look at the jet and see exactly where the real, heavy ingredients are and where the fake stuff is. They can finally taste the real soup.

2. The Process: The "Magnetic Pressure Cooker"

The black hole is surrounded by a disk of swirling gas (the accretion disk). It acts like a pressure cooker for magnetic fields.

• The Build-up: As gas swirls in, it drags magnetic field lines with it, twisting them tighter and tighter, like winding a rubber band.
• The Eruption: Eventually, the magnetic pressure gets too high. The rubber band snaps! This is called a magnetic flux eruption.
• The Result: When the rubber band snaps, it violently pushes everything in the center outward. It creates a vacuum in the middle, sucking the heavy ingredients out of the center and shooting them up the sides of the "jet funnel."

3. The Discovery: It's a Rhythmic, Episodic Cycle

The paper found that the jet isn't a steady stream of heavy matter. It's episodic (happening in bursts).

• The Cycle:
  1. Accretion: Heavy matter flows in.
  2. Eruption: The magnetic field snaps, clearing out the center and pushing a burst of heavy matter up the sides of the jet.
  3. Starvation: The center becomes empty of heavy matter (charge-starved). The jet is now mostly pure energy.
  4. Refill: The magnetic pressure drops, and heavy matter flows back in to start the cycle again.

It's like a garden hose with a kink. You turn it on, the water builds up pressure, the kink releases with a splash (erupting water), the hose goes dry for a second, and then the water flows again.

4. The Spin Factor: The "Whirlpool" Effect

The authors tested three types of black holes: spinning forward, spinning backward, and not spinning at all.

• Spinning Black Holes: The spin creates a shear (like two layers of water sliding past each other at different speeds). This creates vortices (swirls) along the edge of the jet. These swirls act like a mixer, grabbing extra heavy matter from the sides and pulling it into the jet.
• Non-Spinning Black Holes: Without the spin, there are no swirls. The jet stays very clean and empty of heavy matter. The "mixer" is off.

5. The Big Reveal: The "Charge-Starved" Jet

The most important finding is about the Goldreich-Julian (GJ) limit.

• Think of the GJ limit as the minimum number of passengers needed on a bus to keep the engine running smoothly. If there are too few passengers, the bus sputters, and dangerous sparks (electric fields) fly everywhere.
• The study shows that for most of the time, the jet is charge-starved. It doesn't have enough heavy matter (passengers) to screen the electric fields.
• Why does this matter?
  • Sparks: When the bus is empty, those dangerous sparks can accelerate particles to insane speeds (Ultra-High-Energy Cosmic Rays).
  • Neutrinos: These sparks might be the source of high-energy neutrinos (ghost particles) detected by observatories like IceCube.
  • The "Spine" vs. The "Skin": The center of the jet (the spine) is almost entirely empty of heavy matter (pure energy), while the heavy matter is mostly stuck in the outer skin (the sheath) or being ejected in bursts.

Summary

This paper gave us a new "dye" to see what's really inside a black hole's jet. They discovered that the jet is a stop-and-go machine.

1. Magnetic fields build up pressure.
2. They snap, clearing the center and shooting out bursts of heavy matter.
3. The center becomes empty, creating a "charge-starved" zone where particles can be accelerated to extreme energies.
4. The spin of the black hole acts like a mixer, pulling more heavy matter into the jet's edges, but the core remains mostly empty.

This helps explain why we see such powerful bursts of energy and high-speed particles from black holes: they aren't just steady streams; they are violent, rhythmic explosions that periodically clear the deck, leaving the stage empty for the most extreme physics in the universe to happen.

Technical Summary

Here is a detailed technical summary of the paper "The baryon content of magnetically arrested black hole disks and jets" by Chow et al.

1. Problem Statement

Relativistic jets from black holes are critical for feedback mechanisms in galaxies and are potential sites for accelerating ultra-high-energy cosmic rays and producing high-energy neutrinos. A central unknown in jet physics is the baryon loading (the amount of protons and heavy ions entrained in the flow).

• The Challenge: In General Relativistic Magnetohydrodynamic (GRMHD) simulations, maintaining numerical stability in highly magnetized, low-density regions (like jet funnels and reconnection layers) requires imposing "density floors" (artificially injecting mass). This numerical necessity obscures the physical baryon content, making it impossible to distinguish between genuine disk-supplied baryons and numerically injected mass.
• Consequences: Without knowing the true baryon density, it is difficult to determine if the plasma is "charge-starved" (density below the Goldreich–Julian limit, allowing for particle acceleration gaps) or if the jet is inertia-dominated. This uncertainty hinders the modeling of non-thermal emission and particle acceleration mechanisms.

2. Methodology

The authors developed a novel diagnostic tool to reconstruct the physical baryon density in GRMHD simulations:

• Passive Eulerian Tracer: Instead of using discrete Lagrangian particles (as in previous works), they implemented a passive scalar field, $\chi_{tr}$.
  • Initialization: $\chi_{tr} = 1$ inside the accretion disk and $\chi_{tr} = 0$ in the ambient medium.
  • Evolution: The tracer is advected with the fluid velocity ($u^\mu$) according to the continuity equation $\nabla_\mu(\rho_{tr} u^\mu) = 0$, where $\rho_{tr} = \rho \chi_{tr}$.
  • Floor Handling: When the code injects numerical density floors to maintain stability, the tracer fraction is adjusted such that the tracer mass density ($\rho_{tr}$) remains conserved. This effectively filters out the artificially injected mass, allowing $\rho_{tr}$ to serve as a proxy for the physical baryon density supplied by the disk.
• Simulation Setup:
  • Code: BHAC (General Relativistic MHD code).
  • Geometry: Axisymmetric (2D) ideal GRMHD simulations.
  • System: Magnetically Arrested Disks (MAD) around black holes with three spin parameters: prograde ($a=0.9375$), non-spinning ($a=0$), and retrograde ($a=-0.9375$).
  • Parameters: Modeled to represent conditions similar to the M87 black hole accretion flow.

3. Key Contributions

1. Novel Diagnostic: Introduction of a passive Eulerian tracer method that successfully decouples physical baryon transport from numerical density floor injections in GRMHD simulations.
2. Mechanism Identification: Identification of magnetic flux eruption cycles as the primary regulator of baryon loading in MAD systems.
3. Spin Dependence: Demonstration that black hole spin significantly alters the baryon entrainment mechanism, specifically through shear-driven vortices at the jet boundary.
4. Charge Starvation Mapping: Quantitative mapping of the Goldreich–Julian (GJ) screening boundary, revealing that jet power is predominantly carried by charge-starved plasma.

4. Key Results

A. Episodic Baryon Loading via Flux Eruptions

• Mechanism: Baryon loading is not continuous but episodic, regulated by magnetic flux eruption cycles.
• Process:
  1. Magnetic flux accumulates near the horizon until magnetic pressure exceeds ram pressure.
  2. A reconnection event occurs in the inner equatorial current sheet, violently ejecting plasma.
  3. This eruption evacuates baryons from the innermost equatorial region, creating a baryon-poor, magnetically dominated cavity (charge-starved region).
  4. Moderately magnetized, baryon-rich disk material is expelled along the jet funnel wall, establishing a transient mass-loading channel.
• Observation: During eruptions, the total mass accretion rate ($\dot{M}$) diverges significantly from the tracer accretion rate ($\dot{M}_{tr}$), confirming that most mass accreted during these phases is numerical floor mass, not physical baryons.

B. Role of Black Hole Spin

• Rotating Black Holes (Prograde/Retrograde):
  • Shear-driven waves develop along the jet boundary (funnel wall).
  • These waves generate vortices that propagate outward, entraining additional ambient baryons and thickening the mixing layer.
  • This mechanism is absent in the non-spinning case.
• Non-Spinning Black Hole:
  • Lacks frame-dragging effects; flux eruptions are delayed and less frequent.
  • The jet boundary remains largely laminar with no shear-driven vortices.
  • The charge-starved region is narrower and more steady compared to rotating cases.

C. Evolution of the Current Sheet

• The plasma feeding the equatorial current sheet evolves from baryon-rich (early eruption phase) to baryon-poor (late eruption phase).
• This transition implies that hadronic emission models (e.g., proton synchrotron) must account for time-dependent composition, as the reconnection layer may be dominated by pairs during the peak of an eruption.

D. Goldreich–Julian (GJ) Screening Boundary

• The authors mapped the surface where plasma density drops below the GJ limit ($\sigma_{tr} > \sigma_{GJ}$).
• Findings:
  • Large flux eruptions temporarily evacuate the inner magnetosphere, extending the charge-starved region from the jet spine down to the equatorial plane.
  • Jet Power Composition: For both prograde and retrograde spins, the majority of the jet's electromagnetic power is carried by baryon-poor plasma with densities below the GJ threshold.
  • In the prograde case, 50% of the EM power is carried by plasma with $\sigma_{tr} < \sigma_{GJ}$ only 16% of the time, but the bulk of the power is consistently carried by low-density plasma.

5. Significance and Implications

• Jet Composition: The study provides strong evidence that relativistic jets in MAD systems are predominantly composed of low-density, charge-starved plasma, even when baryons are present at the boundaries. This supports models where particle acceleration occurs in vacuum gaps (spark gaps) near the jet base.
• High-Energy Emission: The time-dependent baryon content of the current sheet suggests that proton synchrotron emission (a candidate for GeV emission in M87) may be highly variable, peaking during the baryon-rich phases of flux eruptions.
• Neutrino Production: The identification of charge-starved regions and reconnection layers as sites for proton acceleration reinforces the potential of these systems as sources of high-energy neutrinos.
• Future Directions: The authors acknowledge the limitations of 2D axisymmetry (which suppresses 3D turbulence and interchange instabilities) and call for 3D simulations incorporating radiative cooling and kinetic effects (resistivity) to refine these findings.

In summary, this paper establishes a robust framework for diagnosing jet composition in GRMHD simulations, revealing that magnetic flux eruptions drive an episodic, spin-dependent cycle of baryon loading that leaves the jet spine largely charge-starved and capable of accelerating particles to ultra-high energies.