Study of Flat Spectrum Radio Quasars and BL Lacertae Objects as Sources of Diffusive Ultra High-Energy Cosmic Rays

Sintesi Tecnica: Studio dei Quasar a Spettro Radio Piatto e degli Oggetti BL Lacertae come Sorgenti di Raggi Cosmici Ultra-Energetici Diffusi

Enunciato del Problema
L'origine dei raggi cosmici ultra-energetici (UHECR, $E \gtrsim 10^{17}$ eV) rimane una questione significativa irrisolta nell'astrofisica moderna. Sebbene l'origine extragalattica sia suggerita dalle caratteristiche spettrali, le sorgenti precise sono difficili da identificare a causa delle deflessioni causate dai campi magnetici galattici ed extragalattici. Dati recenti ad alta precisione dall'Osservatorio Pierre Auger e dal Telescope Array (TA) rivelano specifiche caratteristiche spettrali (indurimento al "tallone", inaridimento e soppressione GZK) e tendenze nella composizione (transizione da nuclei leggeri a nuclei più pesanti con l'aumentare dell'energia). Sebbene le galassie radio vicine siano candidati forti, il potenziale contributo dei Nuclei Galattici Attivi (AGN), in particolare dei Quasar a Spettro Radio Piatto (FSRQ) e degli oggetti BL Lacertae (BL Lac), richiede un rigoroso test contro i vincoli osservativi riguardanti il flusso, la composizione e l'anisotropia.

Metodologia
Gli autori investigano gli FSRQ e i BL Lac come sorgenti plausibili di UHECR diffusi utilizzando un modello di propagazione completo all'interno della cosmologia standard $\Lambda$CDM.

1. Modelli di Evoluzione delle Sorgenti:

   • Lo studio utilizza funzioni realistiche di Evoluzione della Densità Dipendente dalla Luminosità (LDDE) derivate dalle funzioni di luminosità nei raggi gamma osservate (LF) dai dati di Fermi-LAT.
   • FSRQ: Modellati utilizzando una funzione LDDE in cui il picco di redshift della densità delle sorgenti dipende dalla luminosità ($z_c \propto L^\alpha$).
   • BL Lac: Modellati utilizzando una distribuzione intrinseca degli indici fotonici (Gaussiana) correlata con la luminosità, combinata con la stessa formalismo LDDE.
   • Questi modelli sostituiscono forme di evoluzione fenomenologiche semplici (ad esempio, leggi di potenza in $1+z$) con evoluzioni motivate fisicamente basate sui dati osservati nei raggi gamma.

2. Fisica della Propagazione:

   • Gli UHECR sono propagati attraverso lo spazio extragalattico considerando le perdite di energia dovute alle interazioni con la Radiazione Cosmica di Fondo (CMB) e la Luce di Fondo Extragalattica (EBL).
   • Processi inclusi: produzione di coppie Bethe-Heitler, produzione di fotopioni (per i protoni) e fotodisintegrazione (per i nuclei). Lo studio impiega l'approssimazione della Perdita di Energia Continua (CEL) per la fotodisintegrazione, tracciando la probabilità di sopravvivenza dei nuclei primari (p, He, N, Si, Fe) senza risolvere esplicitamente le cascate secondarie.
   • Diffusione Magnetica: La propagazione avviene in un campo magnetico intergalattico (IGMF) uniforme e turbolento con intensità RMS $B_{rms} = 1$ nG e lunghezza di coerenza $l_c = 0.5$ Mpc. Il coefficiente di diffusione dipendente dall'energia $D(E)$ è calcolato basandosi sulla turbolenza di Kolmogorov, separando i regimi di diffusione risonante e non risonante.

3. Calcolo del Flusso e Adattamento:

   • Il flusso diffuso è calcolato integrando l'emissività della sorgente su redshift e luminosità, applicando un teorema di propagazione per sommare i contributi da una distribuzione uniforme di sorgenti sintetiche.
   • Il modello include una composizione nucleare mista (p, He, N, Si, Fe) con frazioni trattate come parametri liberi.
   • La normalizzazione dell'emissività della sorgente è collegata alla luminosità nei raggi gamma tramite un fattore di carico barionico $\kappa = L_{CR}/L_\gamma$.
   • I parametri (distanza di separazione delle sorgenti $d_s$, efficienza di accelerazione $\eta$ e frazioni nucleari $\epsilon_Z$) sono adattati ai dati osservativi di Auger e TA (flusso e profondità media della doccia $\langle X_{max} \rangle$).
   • Il confronto statistico dei modelli è eseguito utilizzando il $\chi^2$ ridotto, il Criterio di Informazione di Akaike (AIC) e il Criterio di Informazione Bayesiano (BIC).
   • Mappe dell'intero cielo del flusso integrale ($10-100$ EeV) e dell'anisotropia dipolare sono generate utilizzando il framework HEALPix.

Risultati Chiave

1. Flusso e Composizione:

   • Sia lo scenario FSRQ che quello BL Lac sono in grado di riprodurre lo spettro energetico osservato degli UHECR e i dati $\langle X_{max} \rangle$ con composizioni miste.
   • L'efficienza di accelerazione migliore è $\eta \approx 0.1$, coerente con le simulazioni numeriche dell'accelerazione nei getti relativistici.
   • Il fattore di carico barionico richiesto è $\kappa \approx 5.65$, indicando che la potenza dei raggi cosmici è energeticamente fattibile (solo alcune volte la luminosità nei raggi gamma).
   • Scenario FSRQ: Richiede una distanza di separazione delle sorgenti $d_s \approx 8 \pm 1$ Mpc per adattare il flusso. Ciò implica una densità di sorgenti significativamente più alta delle stime osservative dalle funzioni di luminosità di Fermi-LAT.
   • Scenario BL Lac: Richiede $d_s \approx 160 \pm 21$ Mpc, il quale è coerente con le stime osservative della densità dei BL Lac.

2. Confronto Statistico:

   • I valori di $\chi^2$ ridotto per i dati combinati di Auger e TA sono 2.37 (FSRQ) e 2.55 (BL Lac).
   • I criteri di selezione del modello (AIC e BIC) favoriscono il modello BL Lac ($\Delta$AIC = 4.17, $\Delta$BIC = 4.83), fornendo prove moderate che i BL Lac offrono una descrizione migliore dei dati rispetto agli FSRQ.

3. Vincoli di Anisotropia:

   • BL Lac: L'anisotropia dipolare prevista utilizzando i parametri adattati al flusso ($z_{max}=0.5$) concorda bene con le osservazioni di Auger.
   • FSRQ: L'ampiezza dipolare prevista è significativamente più alta di quella osservata quando si utilizzano gli stessi parametri. Per corrispondere all'anisotropia osservata, la gamma di redshift deve essere artificialmente ristretta a $z_{max} \approx 0.01$, il che contraddice i requisiti di adattamento del flusso e la nota forte evoluzione cosmologica degli FSRQ.

4. Fotodisintegrazione:

   • I calcoli dell'opacità ottica e delle frazioni assorbite mostrano che i nuclei pesanti (ad esempio Fe) subiscono una significativa fotodisintegrazione a energie più basse rispetto ai nuclei più leggeri (ad esempio He), portando a una soppressione progressiva dei componenti pesanti alle energie più elevate.

Significato e Conclusioni
Il paper conclude che, sebbene sia gli FSRQ che i BL Lac siano capaci di riprodurre il flusso osservato di UHECR, lo scenario BL Lac è fortemente favorito quando si considera l'insieme completo dei vincoli osservativi.

• Gli FSRQ sono sfavoriti come sorgenti dominanti di UHECR diffusi perché lo scenario richiede una densità di sorgenti irrealisticamente alta (piccolo $d_s$) per adattare il flusso e prevede un'anisotropia dipolare che supera i limiti osservativi a meno che la popolazione di sorgenti non sia artificialmente troncata a bassi redshift.
• I BL Lac rimangono coerenti con le osservazioni, fornendo una descrizione autoconsistente del flusso, della composizione e dell'anisotropia dipolare con una densità di sorgenti e un'evoluzione di redshift che si allineano ai dati di Fermi-LAT.

Gli autori suggeriscono che, sebbene gli FSRQ possano ancora contribuire ai canali multi-messaggero (neutrini e raggi gamma), è probabile che siano contributori subdominanti al flusso diffuso di UHECR rispetto agli oggetti BL Lac. Lo studio sottolinea l'importanza di utilizzare modelli di evoluzione delle sorgenti realistici e dipendenti dalla luminosità piuttosto che forme fenomenologiche semplificate quando si vincolano le origini degli UHECR.