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

Sintesi Tecnica: Accelerazione di Protoni tramite Shock Collisionless nelle Corone di Buchi Neri Supermassivi

Definizione del Problema
Recenti osservazioni dell'Osservatorio Neutrinico IceCube hanno rilevato un eccesso significativo di neutrini ad alta energia provenienti da nuclei galattici attivi (AGN) radio-silenti nelle vicinanze, in particolare la galassia Seyfert NGC 1068. Queste osservazioni presentano una discrepanza: sebbene vengano rilevati neutrini ad alta energia, il corrispondente flusso di raggi gamma ad alta energia è significativamente inferiore a quanto previsto dai modelli di interazione adronica standard. Ciò suggerisce che i neutrini siano prodotti in una regione opaca ai raggi gamma, probabilmente la calda corona che circonda il buco nero supermassivo (SMBH) centrale, dove i raggi gamma vengono attenuati tramite la produzione di coppie $\gamma\gamma$ con fotoni X.

Per spiegare lo spettro dei neutrini osservato, i protoni genitori devono essere accelerati a energie di $\sim 100$ TeV all'interno di questa corona. Tuttavia, il meccanismo di accelerazione dominante nella corona degli SMBH rimane non stabilito. Sebbene l'Accelerazione Diffusiva da Shock (DSA) sia un candidato principale, la sua viabilità è incerta nelle condizioni coronali. Specificamente, gli shock coronali dovrebbero avere numeri di Mach sonici bassi ($M_s \sim 2-10$) e velocità trans-relativistiche, regimi tradizionalmente considerati inefficienti per l'accelerazione delle particelle. Inoltre, il comportamento dei plasmi a due temperature ($T_i \gg T_e$) e l'efficienza specifica dell'accelerazione protonica rispetto a quella elettronica in questi ambienti mancano di vincoli basati sui primi principi.

Metodologia
Per affrontare queste incertezze, gli autori impiegano simulazioni particle-in-cell (PIC) unidimensionali basate sui primi principi utilizzando il codice SMILEI. Lo studio si concentra su shock quasi-paralleli ($\theta = 30^\circ$) sintonizzati sui parametri fisici delle corone degli SMBH.

Il sondaggio delle simulazioni esplora sistematicamente un ampio spazio di parametri motivato dai vincoli osservativi per le galassie Seyfert:

• Numero di Mach Sonico ($M_s$): Da trans-sonico ($M_s \approx 2$) a valori più elevati ($M_s \approx 30$).
• Numero di Mach Alvénico ($M_A$): Da $4$a$46$.
• Velocità dello Shock ($v_{sh}$): Velocità trans-relativistiche ($v_{sh}/c \approx 0.17 - 0.66$).
• Rapporto di Temperatura Ione-Elettrone ($T_i/T_e$): Da $0.01$a$50$, che riflette gli stati termici disaccoppiati attesi nelle coronae collisionless.

Le simulazioni utilizzano un rapporto di massa iono-elettrone ridotto ($m_i/m_e = 100$) per garantire la fattibilità computazionale pur preservando i parametri adimensionali chiave. Le corse sono state evolute per un tempo sufficiente ($t \gtrsim 400 \Omega_{ci}^{-1}$) per permettere la saturazione delle instabilità elettromagnetiche e lo sviluppo di code di legge di potenza non termiche.

Risultati Chiave

1. Accelerazione Protonica Robusta: Lo studio rileva che la DSA è un meccanismo robusto ed efficiente per l'accelerazione dei protoni nelle corone degli SMBH. Nonostante i bassi numeri di Mach sonici, gli shock canalizzano costantemente circa il 10% dell'energia cinetica dello shock in ioni non termici. Questa efficienza rimane stabile in tutto il sondaggio dei parametri, inclusi i casi con $M_s \approx 2$, sfidando l'idea che gli shock a basso $M_s$ siano acceleratori inefficaci.
2. Dinamica di Accelerazione: L'energia massima dei protoni accelerati aumenta linearmente con il tempo dopo che la struttura dello shock raggiunge uno stato quasi stazionario. Il tasso di accelerazione segue la scala prevista dalla teoria DSA, $t_{acc} \propto \gamma \beta_{sh}^{-2} \Omega_{ci}^{-1}$. Gli autori derivano un coefficiente di accelerazione normalizzato, $C^{-1} \equiv d(\gamma \beta_{sh}^{-2})/d(t\Omega_{ci})$, che presenta un valore universale di $C^{-1} \approx 0.1$, variando di meno di un fattore due in tutti i parametri testati.
3. Inefficienza dell'Accelerazione Elettronica: Al contrario dei protoni, l'accelerazione elettronica è altamente inefficiente, con frazioni di energia elettronica non termica che rimangono al di sotto dell'1% (tipicamente $\sim 10^{-3}$ per i rapporti $T_i/T_e$ coronali). L'efficienza dell'accelerazione elettronica diminuisce all'aumentare del rapporto $T_i/T_e$. In specifici casi trans-relativistici con amplificazione del campo magnetico, la coda di legge di potenza elettronica non riesce a estendersi a energie più elevate nel tempo, suggerendo che la DSA sostenuta sia inibita quando i campi amplificati rendono lo shock effettivamente superluminale per gli elettroni.
4. Meccanismi di Instabilità: La turbolenza magnetica a monte necessaria per lo scattering è guidata principalmente dall'instabilità non risonante (NRI, o instabilità di Bell) nella maggior parte dello spazio dei parametri. L'instabilità risonante (RI) diventa dominante solo a numeri di Mach alvénici molto bassi ($M_A \lesssim 8$).

Implicazioni per NGC 1068 e la Produzione di Neutrini
Applicando questi risultati a NGC 1068, gli autori dimostrano che:

• Scale Temporali: I protoni possono essere accelerati a $\sim 100$ TeV su scale temporali più brevi delle relative scale di raffreddamento (Bethe-Heitler, fotomesone, $pp$) e dinamiche (caduta libera).
• Budget Energetico: La potenza protonica non termica necessaria per spiegare l'osservata luminosità dei neutrini è stimata in $\sim 3\%$ della potenza totale dello shock. Questo è ampiamente entro il $\sim 10\%$ dell'efficienza di accelerazione misurata nelle simulazioni.
• Deficit di Raggi Gamma: La bassa efficienza dell'accelerazione elettronica e l'alta opacità della corona ai raggi gamma spiegano naturalmente il deficit osservato nei raggi gamma ad alta energia, consentendo al contempo la fuga dei neutrini.

Significato e Rivendicazioni
L'articolo rivendica di fornire i primi vincoli, basati sui primi principi e auto-consistenti, sulla DSA nel contesto specifico delle corone degli SMBH. Il significato primario risiede nell'aver stabilito che gli shock collisionless sono un meccanismo vitale ed efficiente per produrre i protoni ad alta energia necessari per spiegare il segnale dei neutrini di IceCube da NGC 1068, anche nelle condizioni difficili di bassi numeri di Mach sonici e alti rapporti di temperatura iono-elettrone.

Gli autori dichiarano esplicitamente che, sebbene i loro risultati supportino fortemente il modello adronico per la produzione di neutrini nelle corone degli AGN, l'applicazione dei risultati delle simulazioni cinetiche locali a un oggetto astrofisico complesso come una corona di AGN richiede una valutazione attenta. L'analisi del budget energetico è presentata come una valutazione di ordine di grandezza basata sulla microfisica stabilita nel documento, colmando il divario tra la fisica del plasma microscopica e i fenomeni macroscopici ad alta energia. Lo studio non pretende di risolvere tutte le incertezze riguardanti l'energia massima dei protoni o la transizione dettagliata delle instabilità in 3D, notando che queste richiedono ulteriori indagini.