Proton Acceleration by Collisionless Shocks in Supermassive Black Hole Coronae: Implications for High-Energy Neutrinos

Technische Samenvatting: Protonversnelling door botsingsloze schokken in de corona van supermassieve zwarte gaten

Probleemstelling
Recente waarnemingen door het IceCube Neutrino Observatory hebben een significante overmaat aan hoogenergetische neutrino's gedetecteerd van nabijgelegen radio-stille actieve galactische kernen (AGN's), met name de Seyfert-sterrenstelsel NGC 1068. Deze waarnemingen presenteren een discrepantie: hoewel er hoogenergetische neutrino's worden gedetecteerd, is de corresponderende flux van hoogenergetische gammastraling aanzienlijk lager dan voorspeld door standaard hadronische interactiemodellen. Dit suggereert dat neutrino's worden geproduceerd in een regio die opaak is voor gammastraling, waarschijnlijk de hete corona die het centrale supermassieve zwarte gat (SMBH) omringt, waar gammastraling wordt gedempt via $\gamma\gamma$ paarproductie met röntgenfotonen.

Om dit waargenomen neutrino-spectrum te verklaren, moeten moederprotonen worden versneld tot energieën van $\sim 100$ TeV binnen deze corona. Echter, het dominante versnellingsmechanisme in SMBH-corona's blijft onvastgesteld. Hoewel Diffuse Schokversnelling (DSA) een leidende kandidaat is, is de levensvatbaarheid ervan onzeker onder coronaal condities. Specifiek worden coronaal schokken verwacht te hebben met lage sonische Mach-getallen ($M_s \sim 2-10$) en trans-relativistische snelheden, regimes die traditioneel als inefficiënt voor deeltjesversnelling worden beschouwd. Voorts ontbreken eerste-principes-beperkingen voor het gedrag van twee-temperatuurplasma's ($T_i \gg T_e$) en de specifieke efficiëntie van proton- versus elektronversnelling in dergelijke omgevingen.

Methodologie
Om deze onzekerheden aan te pakken, maken de auteurs gebruik van eerste-principes, eendimensionale particle-in-cell (PIC) simulaties met de SMILEI-code. De studie richt zich op quasi-parallelle botsingsloze schokken ($\theta = 30^\circ$) afgestemd op de fysieke parameters van SMBH-corona's.

De simulatie-survey verkent systematisch een breed parameterruimte gemotiveerd door observationele beperkingen voor Seyfert-stelsels:

• Sonisch Mach-getal ($M_s$): Variërend van trans-sonisch ($M_s \approx 2$) tot hogere waarden ($M_s \approx 30$).
• Alfvénisch Mach-getal ($M_A$): Variërend van $4$ tot $46$.
• Schoksnelheid ($v_{sh}$): Trans-relativistische snelheden ($v_{sh}/c \approx 0.17 - 0.66$).
• Ionen-naar-elektrontemperatuurratio ($T_i/T_e$): Variërend van $0.01$ tot $50$, wat de ontkoppelde thermische staten reflecteert die in botsingsloze corona's worden verwacht.

De simulaties maken gebruik van een gereduceerde ion-naar-elektron massa ratio ($m_i/m_e = 100$) om computationele haalbaarheid te garanderen terwijl cruciale dimensieloze parameters behouden blijven. De runs werden geëvolueerd voor voldoende tijd ($t \gtrsim 400 \Omega_{ci}^{-1}$) om de verzadiging van elektromagnetische instabiliteiten en de ontwikkeling van niet-thermische power-law staarten toe te laten.

Belangrijkste Resultaten

1. Robuuste Protonversnelling: De studie stelt vast dat DSA een robuust en efficiënt mechanisme is voor protonversnelling in SMBH-corona's. Ondanks de lage sonische Mach-getallen kanaliseren schokken consistent ongeveer 10% van de kinetische energie van de schok in niet-thermische ionen. Deze efficiëntie blijft stabiel over de gehele parametersurvey, inclus\n bij gevallen met $M_s \approx 2$, wat het idee uitdaagt dat lage-$M_s$ schokken ineffectieve versnellers zijn.
2. Versnellingsdynamica: De maximale energie van versnelde protonen neemt lineair toe met de tijd nadat de schokstructuur een quasi-stationaire toestand heeft bereikt. De versnellingssnelheid volgt de schaling voorspeld door DSA-theorie, $t_{acc} \propto \gamma \beta_{sh}^{-2} \Omega_{ci}^{-1}$. De auteurs leiden een genormaliseerde versnellingscoëfficiënt af, $C^{-1} \equiv d(\gamma \beta_{sh}^{-2})/d(t\Omega_{ci})$, die een universele waarde van $C^{-1} \approx 0.1$ vertoont, variërend met minder dan een factor twee over alle geteste parameters.
3. Inefficiëntie van Elektronversnelling: In tegenstelling tot protonen is de elektronversnelling zeer inefficiënt, waarbij de niet-thermische elektronenergiefracties onder de 1% blijven (typisch $\sim 10^{-3}$ voor coronaal-achtige $T_i/T_e$ ratio's). De efficiëntie van de elektronversnelling neemt af naarmate de $T_i/T_e$ ratio toeneemt. In specifieke trans-relativistische gevallen met sterke magnetische veldversterking, faalt de niet-thermische power-law staart van elektronen om naar hogere energieën te strekken over de tijd, wat suggereert dat aanhoudende DSA wordt geremd wanneer versterkte velden de schok effectief incidenteel superluminaal maken voor elektronen.
4. Instabiliteitsmechanismen: De upstream magnetische turbulentie die nodig is voor verstrooiing wordt primair gedreven door de niet-resonante instabiliteit (NRI, of Bell-instabiliteit) over het grootste deel van de parameterruimte. Resonante instabiliteit (RI) wordt dominant bij zeer lage Alfvénische Mach-getallen ($M_A \lesssim 8$).

Implicaties voor NGC 1068 en Neutrino-productie
Door de resultaten toe te passen op NGC 1068, demonstreren de auteurs dat:

• Tijdschalen: Protonen kunnen worden versneld tot $\sim 100$ TeV op tijdschalen die korter zijn dan de relevante koelings (Bethe-Heitler, fotomeson, $pp$) en dynamische (vrije val) tijdschalen.
• Energiebudget: De niet-thermische protonenergie die nodig is om de geobserveerde neutrino-luminositeit te verklaren, wordt geschat op $\sim 3\%$ van de totale schokenergie. Dit valt ruim binnen de $\sim 10\%$ versnellings-efficiëntie die in de simulaties is gemeten.
• Gamma-straaltekort: De lage efficiëntie van elektronversnelling en de hoge opaciteit van de corona voor gammastraling verklaren op natuurlijke wijze het geobserveerde tekort aan hoogenergetische gammastraling, terwijl neutrino's wel kunnen ontsnappen.

Betekenis en Claims
Het artikel beweert de eerste eerste-principes, zelfconsistente beperkingen te bieden op DSA in de specifieke context van SMBH-corona's. De primaire betekenis ligt in het vaststellen dat botsingsloze schokken een levensvatbaar en efficiënt mechanisme zijn voor het produceren van de hoogenergetische protonen die nodig zijn om het IceCube-neutrino-signaal van NGC 1068 te verklaren, zelfs onder de uitdagende condities van lage sonische Mach-getallen en hoge ion-naar-elektron temperatuurratio's.

De auteurs geven expliciet aan dat hoewel hun resultaten de hadronische modellen voor neutrino-productie sterk ondersteunen, de toepassing van lokale kinetische simulatieresultaten op een complex astrofysisch object zoals een AGN-corona zorgvuldige overweging vereist. De energiebudgetanalyse wordt gepresenteerd als een orde van grootte beoordeling gebaseerd op de microfysica die in het artikel wordt vastgesteld, waarmee de kloof tussen microscopische plasmafysica en macroscopische hoogenergetische fenomenen wordt overbrugd. De studie claimt niet alle onzekerheden met betrekking tot de maximale protonenergie of de gedetailleerde transitie van instabiliteiten in 3D op te lossen, waarbij wordt opgemerkt dat deze verder onderzoek vereisen.