Self-Consistent Modelling of Neutrino Production in Turbulent Black Hole Coronae

Probleemstelling

In 2022 rapporteerde de IceCube-samenwerking een statistisch significant overschot aan hoge-energie neutrino's (1–10 TeV) dat ruimtelijk geassocieerd is met de nabije Seyfert 2-galaxie NGC 1068. Dit suggereert dat niet-gejette actieve galactische kernen (AGN) belangrijke bronnen zijn van extragalactische neutrino's. Een cruciaal raadsel is de discrepantie tussen het waargenomen neutrinosignaal en de zeer lage flux van gammastraling; het intrinsieke gamma-signaal wordt sterk onderdrukt door $\gamma\gamma$-absorptie in een dichte stralingsomgeving.

De fysieke uitdaging ligt in het modelleren van protonversnelling in de turbulente corona van een superzwaar zwart gat. Bestaande methoden hebben beperkingen:

1. PIC-simulaties (Particle-in-Cell): Kunnen microfysische processen simuleren, maar kunnen niet de enorme dynamische bereik (minimaal 5 ordes van grootte) overbruggen die nodig is om protonen te versnellen tot de door IceCube waargenomen energieën (tot 100 TeV) binnen realistische tijdschalen.
2. Standaard transportmodellen: Vaak gebruiken ze stationaire aannames, semi-analytische spectra of negeren ze de terugkoppeling (backreaction) van versnelde deeltjes op de turbulente cascade. Dit is problematisch omdat de geïnfereerde protonluminositeit hoog is, wat suggereert dat versnelde deeltjes de turbulente energie zelf kunnen dempen en zo het versnellingsproces beïnvloeden.

Er is behoefte aan een zelfconsistent, tijdsafhankelijk model dat niet-lineaire kinetische processen, turbulente versnelling, stralingsverliezen en hadronische interacties gelijktijdig behandelt.

Methodologie: De Turb-AM3 Code

De auteurs introduceren een nieuwe hybride numerieke code, Turb-AM3, die de volgende componenten combineert:

1. AM3 (Radiative Solver): Een state-of-the-art, tijdsafhankelijke lepto-hadronische radiatieve solver die alle relevante stralings- en hadronische verliezen berekent, inclusief de ontwikkeling van elektromagnetische cascades en de productie van secundaire deeltjes (fotonen, neutrino's).
2. Stochastische Versnelling Module: Een module die de versnelling van geladen deeltjes in magnetized turbulente plasma's beschrijft via een transportvergelijking in impulsruimte.
   • Het gebruikt geavanceerde operatoren (Fokker-Planck voor quasi-lineaire interacties en een "Generalized Fermi" operator voor sterk turbulente regimes) gebaseerd op recente theoretische vooruitgang en PIC-simulaties.
   • Het modelleert injectie van protonen uit het thermische reservoir.
3. Zelfconsistent Terugkoppeling (Backreaction): Een kerninnovatie is de koppeling tussen de versnelde protonen en de turbulente cascade. De code lost een vergelijking voor de turbulente energiespectrum op die rekening houdt met de demping van turbulentie door de energie-overdracht aan de versnelde protonen. Dit zorgt voor een zelfregulerend mechanisme: als te veel energie naar protonen gaat, neemt de turbulentie af, wat de versnelling vertraagt.
4. Fysische Opstelling:
   • Doel: De corona van NGC 1068 (SMBH massa $\sim 10^{7.2} M_\odot$, corona grootte $R_{cor} \sim 15 r_g$).
   • Omgeving: Een hete, magnetische plasma met een thermische elektronenverdeling (Maxwelliaans) en een niet-thermische protoncomponent.
   • Transport: Deeltjes ondergaan diffusie, advektie (inwaartse stroming naar het zwarte gat) en stochastische versnelling.

Belangrijkste Bijdragen

• Ontwikkeling van Turb-AM3: De eerste code die stochastische versnelling in turbulente coronae volledig koppelt aan een tijdsafhankelijke radiatieve solver met inbegrip van de dynamische demping van turbulentie door versnelde deeltjes.
• Zelfconsistentie: Het model behandelt gelijktijdig de niet-lineaire aspecten van versnelling, stralingsverliezen en de evolutie van het magnetische turbulentiespectrum zonder vooraf bepaalde spectra of stationaire aannames.
• Template Spectra: Het bieden van gedetailleerde template spectrale vormen voor TeV-neutrino's die afhankelijk zijn van de fysieke parameters van de corona.

Resultaten

1. Reproductie van IceCube-data: Met een standaard set fysische parameters ($R_{cor} \approx 15 r_g$, $v_A \approx 0.25c$, $\ell_c \approx 4 r_g$, injectiefractie $\xi_p = 0.1$) reproduceert het model het waargenomen neutrinosignaal van NGC 1068 zeer goed.
   • Het neutrinospectrum toont een karakteristiek profiel: een vlakke uitbreiding bij lage energieën (gedreven door $pp$-interacties) en een piek rond 1–10 TeV (gedreven door $p\gamma$-interacties met de rijke röntgenstraling van de corona).
   • Het model blijft consistent met de strenge bovenlimieten van gammastraling (Fermi-LAT, MAGIC), omdat de hoge-energie gammastraling effectief wordt geabsorbeerd en gereprocesserd in de corona.
2. Zelfregulerend Mechanisme: De simulaties tonen aan dat de versnelde protonen turbulentie dempen zodra hun energiedichtheid vergelijkbaar wordt met de magnetische energiedichtheid. Dit beperkt de maximale protonenergie en zorgt ervoor dat de totale niet-thermische protonluminositeit onder de bolometrische luminositeit blijft. Dit verklaart de normalisatie van de flux zonder fijne afstemming van parameters.
3. Parameter Ruimte Onderzoek:
   • De spectrale vorm is robuust voor injectiefracties $\xi_p \gtrsim 10^{-2}$.
   • De verhouding tussen ontsnappingstijd en versnellingstijd ($t_{esc}/t_{acc}$) is cruciaal; alleen als $t_{esc}/t_{acc} \gtrsim 1$ kunnen protonen voldoende energie opbouwen.
   • De advektietijd ($t_{adv}$) bepaalt de hoog-energetische cutoff.
4. Toepasbaarheid op Andere Bronnen: Het model kan het neutrinosignaal van andere Seyfert-galaxieën (NGC 4151 en NGC 7469) verklaren door simpelweg de schijfluminositeit (en dus de doelfotonendichtheid) aan te passen, wat de waargenomen trend van lagere flux bij hogere neutrino-energieën verklaart.
5. Template Voorspelling: De auteurs voorspellen een universeel template-spectrum voor turbulente versnelling: een brede piek rond enkele TeV met een vlakke staart bij lagere energieën en een scherpe cutoff boven de 10 TeV door fotohadronische verliezen.

Significantie

Dit werk biedt een krachtig nieuw kader voor het interpreteren van multi-messenger signalen uit compacte astrofysische omgevingen.

• Bevestiging van het Mechanisme: Het ondersteunt het idee dat de corona van een AGN een efficiënte bron is van hoge-energie neutrino's via stochastische versnelling in turbulente magnetische velden.
• Diagnostisch Hulpmiddel: De specifieke vorm van het neutrinospectrum (vooral de verhouding tussen de $pp$- en $p\gamma$-componenten) kan dienen als een directe sonde voor de condities in de corona, zoals de dichtheid van deeltjes, de mate van paarbelasting (pair-loading) en de eigenschappen van de turbulentie.
• Toekomstige Observaties: De afgeleide template spectra helpen bij het optimaliseren van zoektochten met toekomstige neutrino-telescopen en kunnen helpen om de eigenschappen van de versnellingsmechanismen in AGN-corona's te beperken zodra er meer data beschikbaar komt.

Samenvattend bewijst dit artikel dat een zelfconsistent model van turbulente versnelling, inclusief terugkoppelingseffecten, de complexe multi-messenger waarnemingen van NGC 1068 en andere Seyfert-galaxieën succesvol kan verklaren zonder in strijd te zijn met gamma-stralingsbeperkingen.