Electron-positron Pair Production in Global GRMHD Simulations of Black Hole Accretion Flows

This paper presents global 3D GRMHD simulations incorporating pair physics to reveal that electron-positron pairs in black hole accretion flows are primarily generated near the disk midplane and transported via advection to the corona and jets, where they can significantly exceed equilibrium values and potentially regulate coronal temperatures.

The Gist

Imagine a black hole not just as a cosmic vacuum cleaner, but as a chaotic, super-hot kitchen where the laws of physics are being tested to their limits. This paper is about what happens when we add a very specific, exotic ingredient to that kitchen: electron-positron pairs.

Think of these pairs as "ghost twins." They are particles of matter and antimatter that are born together, dance around, and then instantly destroy each other (annihilate) if they get too close. The big question astronomers have been asking is: Do these ghost twins play a major role in how black holes eat and shoot out jets of energy?

Here is the story of what the researchers found, explained simply:

1. The Setup: A Cosmic Kitchen Simulation

The scientists used a super-computer to simulate a black hole eating gas. Usually, these simulations treat the gas as a simple fluid. But in this study, they added a special "tracker" for the ghost twins (the positrons). They wanted to see if the twins could be created fast enough to change the recipe of the black hole's atmosphere.

They set up different "recipes" (models) based on how fast the black hole was eating (accretion rate) and how hot the gas was.

2. The "Ghost Nursery" (The Pair Source)

The most exciting discovery is where these twins are born.

• The Sweet Spot: The twins aren't born everywhere. They are born in a very specific, thin ring just above the black hole's "equator" (the accretion disk). Think of this as a nursery located right at the edge of the danger zone.
• The Balance: In this nursery, the gas is thick enough and hot enough that twins are born and destroyed at the same rate. It's a busy, balanced factory.
• The "Pair Void": Closer to the black hole itself, the twins disappear. The gas is so dense and cold (relatively speaking) that any twins that wander in get instantly destroyed. It's like a "ghost-free zone" right next to the black hole.

3. The Conveyor Belt: Advection

This is the paper's biggest "Aha!" moment.

• The Problem: In the upper atmosphere of the black hole (the "corona") and in the powerful jets shooting out into space, the gas is very thin. In thin gas, it takes a long time for twins to be born naturally. It's like trying to make a cake in an empty kitchen; you have to wait forever for ingredients to appear.
• The Solution: The black hole acts like a conveyor belt. The twins born in the dense nursery near the disk are swept up by the swirling gas and flung up into the thin corona and the jets.
• The Result: The upper atmosphere and the jets are filled with twins not because they were born there, but because they were transported there. The conveyor belt (advection) is moving them faster than they can be destroyed.

4. The "Goldilocks" Zone

The researchers found that the twins only stay in a stable, balanced state in one specific zone (the nursery).

• Too Thin (Jets/Upper Corona): The twins are out of balance. There are way more twins there than physics says should be there naturally, because the conveyor belt dumped them there faster than they could die.
• Too Dense (Near the Black Hole): The twins are wiped out instantly.
• Just Right (The Nursery): This is where the twins are in perfect equilibrium, born and dying at the same speed.

5. Why Does This Matter?

• Lighting the Jet: Black holes shoot out massive beams of energy (jets). To do this, they need to carry an electric charge. The study suggests that these transported ghost twins might provide the necessary "charge" to keep the jets running smoothly, acting like the fuel for the engine.
• Temperature Control: The birth and death of these twins act like a thermostat. If the gas gets too hot, the twins multiply and cool it down. If it gets too cold, they stop being born. This might explain why black hole atmospheres don't get infinitely hot.
• X-Ray Binaries: The amount of twins found in the simulation (about 1% of the total particles) matches what we see in real-life observations of smaller black holes (X-ray binaries). This suggests the simulation is realistic.

The Big Takeaway

Imagine a black hole as a giant, swirling storm.

• The twins are like raindrops.
• The nursery near the disk is the cloud where the rain is made.
• The conveyor belt (advection) is the wind that sweeps the rain up into the upper atmosphere and shoots it out into space.

Before this study, we thought the rain in the upper atmosphere had to be made right there. This paper shows that the wind is actually carrying the rain from the clouds below. This "wind transport" is a crucial, previously overlooked ingredient in how black holes power their spectacular jets and regulate their own temperature.

In short: Black holes don't just make ghost twins locally; they import them from the disk and shoot them out into the universe, and this process is essential for how these cosmic monsters work.

Technical Summary

1. Problem Statement

The role of electron-positron ($e^\pm$) pairs in black hole accretion flows remains a critical but under-explored aspect of high-energy astrophysics. While pairs are believed to be essential for screening electric fields in Blandford-Znajek jets and regulating plasma temperatures in accretion coronas, their global dynamics are difficult to model.

• The Gap: Previous studies relied on analytical one-zone models, post-processing radiative transfer on GRMHD snapshots, or particle-in-cell (PIC) simulations which are computationally expensive for global scales.
• The Challenge: It is unclear whether pairs reach local thermodynamic equilibrium (LTE) or if they are dominated by global advection. Furthermore, the mechanisms supplying sufficient pair density to jets (to exceed the Goldreich-Julian limit) and the impact of pairs on the global accretion structure (e.g., temperature regulation) require a self-consistent, global treatment.

2. Methodology

The authors performed global, three-dimensional General Relativistic Magnetohydrodynamic (GRMHD) simulations using the open-source code iharm3d.

• Simulation Setup:

  • Geometry: Kerr black holes with spin $a = +0.9375$.
  • Accretion State: Simulations were run in the SANE (Stable and Normal Evolving) state using a Fishbone-Moncrief torus initial condition.
  • Models: A suite of models varying the accretion rate (represented by the aspect ratio $H/r = 0.1, 0.2, 0.3$) and black hole mass ($10^8 M_\odot$ and $10 M_\odot$).
  • Cooling: An ad-hoc cooling function was applied to maintain specific disk aspect ratios ($H/r$), simulating radiative losses without full radiative transfer.

• Pair Physics Implementation:

  • Passive Scalar Approach: Instead of a single-fluid approximation, the authors explicitly evolved the positron mass density ($\rho_+$) as a passive scalar advected by the fluid velocity field.
  • Charge Neutrality: Enforced via $n_e = n_+ + n_p$, allowing unique determination of electron, positron, and proton densities.
  • Reaction Rates: Pair production ($\dot{n}_C$) and annihilation ($\dot{n}_A$) rates were calculated locally based on the Thomson scattering optical depth ($\tau_p$) and electron temperature ($T_e$).
  • Temperature Assumption: Electrons and positrons were assumed to maintain a constant temperature ($T_e = 10^9 - 10^{10.5}$ K).
  • Runaway Suppression: To prevent unphysical "pair runaway" (infinite pair production) in optically thick regions where no equilibrium solution exists, pair production was manually disabled when $\tau_p > \tau_{crit}$. Annihilation was still allowed, justified by the expectation that Coulomb collisions would cool the plasma to $\sim 10^7$ K in these dense regions.

3. Key Contributions

• Global Advection of Pairs: This is one of the first studies to explicitly track the spatiotemporal distribution of pairs in a global 3D GRMHD simulation, moving beyond local equilibrium assumptions.
• Advection as a Dominant Source: The study identifies advection from the disk as a primary mechanism for injecting pairs into the upper corona and jets, challenging the view that local production is the sole source.
• Pair Void Identification: The discovery of a "pair void" extending a few gravitational radii from the black hole in high-accretion-rate models, caused by the suppression of production in the optically thick midplane.
• Goldreich-Julian Screening: Demonstration that advected pairs can populate the magnetically dominated polar regions with densities exceeding the Goldreich-Julian limit, potentially enabling electric field screening in jets.

4. Key Results

A. Timescales and Equilibrium

• Equilibrium vs. Advection: The distribution of pairs is governed by the competition between the pair equilibrium timescale ($t_{eq}$) and the dynamical (inflow) timescale ($t_{inflow}$).
  • High $\tau_p$ (Disk Midplane/Coronal Base): $t_{eq} \ll t_{inflow}$. Pairs are in near-equilibrium ($z \approx z_{eq}$).
  • Low $\tau_p$ (Upper Corona/Jets): $t_{eq} \gg t_{inflow}$. Pairs are far from equilibrium ($z \gg z_{eq}$) because they are transported faster than they can annihilate or equilibrate locally.
• Pair Fraction: In high accretion rate models ($H/r=0.1$), the maximum pair fraction reaches $z \sim 0.01$, consistent with X-ray binary observations. The ratio $z/z_{eq}$ can exceed $10^2$ in the jets, indicating significant disequilibrium.

B. Spatial Distribution

• Pair Void: In high accretion rate models, a region near the midplane close to the horizon is devoid of pairs. This is due to the high optical depth suppressing production (via the $\tau_{crit}$ cutoff) and rapid annihilation.
• Source Strip: A thin strip at the base of the corona (where $\tau_p \approx 1$) acts as the primary "pair source," maintaining near-equilibrium and generating a high flux of pairs.
• Plunging Region: Inside the ISCO, advection dominates. Pairs are rapidly swept inward, preventing local equilibrium even where $\tau_p \sim 1$.

C. Production Mechanisms

• Photon-Photon ($\gamma\gamma$): Dominates pair production in the optically thick midplane and coronal base.
• Particle-Particle ($ee$): Dominates in the optically thin upper corona and jets, though the absolute production rate is low.
• Advection: The primary transport mechanism moving pairs from the "source strip" to the jets and upper corona.

D. Parameter Sensitivity

• Accretion Rate: Higher accretion rates shift the "pair-maximized strip" closer to the horizon and increase the opening angle, enhancing pair capture by outflows.
• Black Hole Mass: Results are largely scale-free regarding the pair fraction ($z$), though the absolute number density scales with mass. The physics holds for both AGN ($10^8 M_\odot$) and X-ray binaries ($10 M_\odot$).

5. Significance and Implications

• Jet Composition: The results suggest that advection can supply sufficient pair density to magnetically dominated jets to exceed the Goldreich-Julian density, offering a mechanism to screen electric fields without requiring local pair cascades in the jet itself.
• Temperature Regulation: The study supports the hypothesis that pair physics could act as a "thermostat" for the coronal temperature. In regions where $t_{eq}$ is comparable to the heating timescale (Coulomb collisions), pairs could regulate the plasma temperature by increasing radiative cooling efficiency.
• Observational Consistency: The maximum pair fraction of $\sim 0.01$ aligns with recent inferences from X-ray binaries, suggesting that accretion disks are generally not pair-dominated, contrary to some early theoretical predictions.
• Limitations & Future Work: The authors note that the assumption of constant electron temperature and the lack of full radiative transfer (specifically synchrotron and non-thermal processes) are simplifications. Future work requires coupling these pair dynamics with two-temperature physics and full radiative transfer to fully resolve the coronal temperature structure.

In conclusion, this paper establishes that advection is a critical, previously underappreciated source of pairs in black hole accretion flows, capable of populating jets and coronae with non-equilibrium pair densities that influence global dynamics and observational signatures.