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

Probleemstelling

Er is een groeiende interesse in de rol van elektron-positronparen (e-paren) in de astrofysica van zwarte gaten, met name in de context van accretievloeiingen en jets. Hoewel bekend is dat paren een cruciale rol spelen bij het afschermen van elektrische velden en het geleiden van stromen in Blandford-Znajek-jets, blijft de mechanisme voor de toevoer van deze plasma's onduidelijk.

• Bestaande modellen: De "gap acceleration"-modellen (partikel-in-cell simulaties) suggereren dat ongescreende elektrische velden paren genereren via cascades. De "drizzle"-modellen (post-processing van GRMHD) suggereren dat paren ontstaan uit botsingen van opgeblazen fotonen in de jet.
• Het gat in het onderzoek: De impact van paar-fysica op de dynamica en thermodynamica van de accretievloeiing zelf is relatief onontgonnen. Voorgaande studies waren vaak analytisch, lokaal, of gebruikten één-component vloeistofmodellen die het aantal paren niet expliciet konden volgen zonder kunstmatige vloeren (floors) die de dichtheid verstoren.
• De centrale vraag: Kunnen elektron-positronparen de temperatuur van het accretieplasma reguleren via een thermostaat-mechanisme, en hoe wordt de ruimtelijke verdeling van paren bepaald door de concurrentie tussen lokale paar-processen en globale advektie?

Methodologie

De auteurs voeren globale, driedimensionale General Relativistic Magnetohydrodynamic (GRMHD) simulaties uit met de open-source code iharm3d.

• Simulatieopzet: Ze modelleren accretievloeiingen rondom Kerr-zwarte gaten (spin $a=0.9375$) in de SANE-toestand (Stable, Normal, Accreting). De initiële conditie is een Fishbone-Moncrief torus.
• Koppelingsmechanisme: In plaats van een één-component vloeistof, introduceren ze positronen als een passieve scalair ($\rho_+$) die wordt meegevoerd door het snelheidsveld van het plasma. Ze lossen een extra continuïteitsvergelijking op voor de positronmassadichtheid.
• Paarproductie en -annihilatie: De evolutie van het aantal paren wordt bepaald door een operator-splitting methode die rekening houdt met:
  • Productie via foton-foton ($\gamma\gamma$), foton-deeltje ($\gamma e$) en deeltje-deeltje ($ee$) botsingen.
  • Annihilatie.
  • De berekening is lokaal gebaseerd op de lokale Thomson-optische diepte ($\tau_T$) en de schaalhoogte ($H$).
• Beperkingen en aannames:
  • De elektron- en positrontemperatuur ($T_e$) wordt constant gehouden (verschillende modellen lopen van $10^9$ tot $10^{10.5}$ K).
  • Er wordt een "floor" (ondergrens) toegepast op de dichtheid om numerieke instabiliteiten te voorkomen.
  • In gebieden met zeer hoge dichtheid (waar $\tau_p > \tau_{crit}$) wordt paarproductie handmatig uitgeschakeld. Dit is gebaseerd op de aanname dat in deze optisch dikke gebieden Coulomb-kolli ties de elektronen afkoelen tot $\sim 10^7$ K, waardoor paarproductie inefficiënt wordt en een "pair runaway" (onbeperkte groei) wordt voorkomen.
• Parameters: Er worden modellen uitgevoerd voor verschillende accretiepercentages (van ADAF tot hoge accretie, $\dot{M} \sim 0.05 - 0.4 \dot{M}_{Edd}$) en schaalhoogtes ($H/r = 0.1, 0.2, 0.3$).

Belangrijkste Resultaten

1. Tijdschalen en Evenwicht

• De paar-evenwichtstijdschaal ($t_{eq}$) varieert sterk. In gebieden met hoge optische diepte ($\tau_p \sim 1$) is $t_{eq}$ kort ($\sim O(1) GM/c^3$), waardoor het plasma in lokaal evenwicht kan komen.
• In gebieden met lage optische diepte (korona en jets) is $t_{eq}$ zeer lang, veel langer dan de dynamische tijdschalen (inval, rotatie). Hier wordt de parenverdeling dus niet bepaald door lokale productie, maar door advektie vanuit het schijfgebied.

2. Ruimtelijke Verdeling: De "Pair Void" en de "Pair Strip"

• Pair Void: In modellen met hoge accretie snelheid ontstaat er een gebied dicht bij de zwarte gat-horizon (in het equatoriale vlak) waar de parenverhouding ($z = n_+/n_p$) extreem laag is. Dit komt doordat de dichtheid zo hoog is dat $\tau_p$ de kritieke drempel overschrijdt, waardoor productie wordt onderdrukt en bestaande paren snel geannihileerd worden.
• Pair Strip: Direct buiten het "plunging region" (binnen de ISCO) en aan de basis van de corona, vormt zich een dunne strook waar $\tau_p \approx 1$. Hier is de paarproductie maximaal en bevindt het plasma zich in lokaal evenwicht. Dit gebied fungeert als de primaire bron van paren.
• Advektie: Paren uit deze dunne strook worden door de bulkstroom (turbulentie en uitstromen) meegevoerd naar de bovenste corona en de jets. In deze gebieden overschrijdt de paarverhouding vaak de lokale evenwichtswaarde met factoren van $10^2$ of meer, puur omdat ze daarheen zijn getransporteerd.

3. Dominante Productiemechanismen

• Dicht bij het vlak (hoge $\tau_p$): De dominante bron is foton-foton botsing (geïnduceerd door Compton-verstrooiing).
• In de jets/korona (lage $\tau_p$): De dominante bron is deeltje-deeltje botsing, hoewel de absolute productiesnelheid hier veel lager is dan in het vlak.

4. Goldreich-Julian Dichtheid en Jets

• De auteurs tonen aan dat de geadvecte paren in de magnetisch gedomineerde polaire jets een dichtheid kunnen bereiken die de Goldreich-Julian-dichtheid ($n_{GJ}$) overschrijdt. Dit suggereert dat advektie een mechanisme is om voldoende ladingsdragers te leveren om het elektrische veld in de jet af te schermen, zelfs zonder lokale paar-cascades in de jet zelf.

5. Invloed van Parameters

• Accretie snelheid: Hogere accretie snelheden verschuiven de "pair-maximized strip" dichter naar het zwarte gat en vergroten de hoek, wat leidt tot efficiëntere opname van paren in uitstromen.
• Zwarte gat massa: De resultaten zijn weinig gevoelig voor de massa van het zwarte gat (getest van $10 M_\odot$ tot $10^8 M_\odot$), wat impliceert dat de bevindingen zowel voor AGN als voor X-ray binaries gelden.
• Maximale paarverhouding: De maximale gevonden paarverhouding is $\sim O(0.01)$, wat consistent is met waarnemingen van X-ray binaries. De auteurs vinden geen bewijs voor een "pair-dominated" schijf ($z \gg 1$).

Bijdrage en Betekenis

1. Advektie als primaire bron: Het artikel benadrukt dat advektie een cruciale, maar vaak verwaarloosde bron van pareninjectie is. Paren worden niet lokaal in de jets geproduceerd, maar vanuit het accretievloeiingsvlak getransporteerd. Dit verklaart hoe jets paren kunnen bevatten zonder dat er lokale cascades nodig zijn.
2. Regulatie van Temperatuur: De studie ondersteunt het idee dat paar-fysica kan fungeren als een thermostaat voor de corona. In gebieden waar de paar-evenwichtstijdschaal kort is, kan de productie van paren de temperatuur van het plasma reguleren door extra koeling (via straling) te introduceren.
3. Validatie van Eén-Zone Modellen: De resultaten bevestigen grotendeels de voorspellingen van analytische "one-zone" evenwichtsmodellen, maar tonen ook aan dat globale dynamica (advektie) leidt tot significante afwijkingen van lokaal evenwicht in de corona en jets.
4. Observationele Implicaties: De bevindingen zijn relevant voor het interpreteren van $\gamma$-stralingssignalen en polarimetrie in X-ray binaries en AGN. Het feit dat de paarverhouding in de jets de Goldreich-Julian limiet kan overschrijden, biedt een oplossing voor het probleem van de ladingsdragerverdeling in Blandford-Znajek-jets.

Beperkingen:
De auteurs erkennen dat hun model vereenvoudigingen bevat, zoals het vasthouden aan een constante elektronentemperatuur, het gebruik van een ad-hoc koelfunctie, en het negeren van synchrotronstraling als bron van zachte fotonen. Toekomstig werk moet volledige stralingstransport en twee-temperatuur fysica integreren om deze mechanismen volledig te kwantificeren.

Conclusie:
Dit werk levert een eerste stap in het integreren van paar-fysica in globale GRMHD simulaties en toont aan dat de ruimtelijke en temporele verdeling van elektron-positronparen wordt gedicteerd door de concurrentie tussen lokale productie (gecontroleerd door optische diepte) en globaal transport (advektie). Dit heeft directe gevolgen voor ons begrip van de structuur van accretieschijven en de samenstelling van relativistische jets.