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

Technische Samenvatting: Studie van Flat Spectrum Radio Quasars en BL Lacertae-objecten als Bronnen van Diffuse Ultra-Hoge-Energie Kosmische Straling

Probleemstelling
De oorsprong van ultra-hoge-energie kosmische straling (UHECR's, $E \gtrsim 10^{17}$ eV) blijft een significant onopgelost vraagstuk in de moderne astrofysica. Hoewel extragalactische oorsprong wordt gesuggereerd door spectrale kenmerken, zijn de precieze bronnen moeilijk te identificeren vanwege afwijkingen door galactische en extragalactische magnetische velden. Recente hoogprecisiegegevens van de Pierre Auger Observatory en het Telescope Array (TA) onthullen specifieke spectrale kenmerken (hardening bij de "enkel", verzachting en GZK-onderdrukking) en samenstellings-trends (overgang van lichte naar zwaardere kernen met toenemende energie). Hoewel nabije radiogalaxieën sterke kandidaten zijn, vereist de potentiële bijdrage van Actieve Galactische Kernen (AGN), specifiek Flat Spectrum Radio Quasars (FSRQ's) en BL Lacertae (BL Lac)-objecten, een rigoureuze toetsing tegen observationele beperkingen met betrekking tot flux, samenstelling en anisotropie.

Methodologie
De auteurs onderzoeken FSRQ's en BL Lacs als waarschijnlijke bronnen van diffuse UHECR's met behulp van een uitgebreid propagatiemodel binnen de standaard $\Lambda$CDM-cosmologie.

1. Bron-evolutiemodellen:

   • De studie maakt gebruik van realistische Luminositeit-afhankelijke Dichtheids-evolutie (LDDE)-functies, afgeleid van waargenomen gammastraal-luminositeitsfuncties (LF) uit Fermi-LAT-gegevens.
   • FSRQ's: Gemodelleerd met een LDDE-functie waarbij de roodverschuivingspiek van de bronnedichtheid afhankelijk is van de luminositeit ($z_c \propto L^\alpha$).
   • BL Lacs: Gemodelleerd met een intrinsieke verdeling van fotonindices (Gaussisch) die gecorreleerd is met luminositeit, gecombineerd met dezelfde LDDE-formulering.
   • Deze modellen vervangen eenvoudige fenomenologische evolutievormen (bijv. machtswetten in $1+z$) door fysisch gemotiveerde evoluties gebaseerd op waargenomen gammastraalgegevens.

2. Propagatiefysica:

   • UHECR's worden door extragalactische ruimte gepropageerd, waarbij energieverliezen door interacties met de Kosmische Microgolfachtergrond (CMB) en Extragalactische Achtergrondlicht (EBL) worden overwogen.
   • Opgenomen processen: Bethe-Heitler-paarproductie, fotopionproductie (voor protonen) en fotodisintegratie (voor kernen). De studie maakt gebruik van de Continuous Energy Loss (CEL)-benadering voor fotodisintegratie, waarbij de overlevingskans van primaire kernen (p, He, N, Si, Fe) wordt gevolgd zonder expliciet secundaire cascades op te lossen.
   • Magnetische Diffusie: Propagatie vindt plaats in een uniform, turbulent intergalactisch magnetisch veld (IGMF) met een RMS-strength $B_{rms} = 1$ nG en een coherentielengte $l_c = 0.5$ Mpc. De energie-afhankelijke diffusiecoëfficiënt $D(E)$ wordt berekend op basis van Kolmogorov-turbulentie, waarbij resonante en niet-resonante diffusieregimes worden gescheiden.

3. Fluxberekening en Fitten:

   • De diffuse flux wordt berekend door de bronemissiviteit te integreren over roodverschuiving en luminositeit, waarbij een propagatiestelling wordt toegepast om bijdragen op te tellen van een uniforme verdeling van synthetische bronnen.
   • Het model omvat een gemengde nucleaire samenstelling (p, He, N, Si, Fe) waarbij fracties als vrije parameters worden behandeld.
   • De normalisatie van de bronemissiviteit is gekoppeld aan de gammastraal-luminositeit via een baryonbeladingsfactor $\kappa = L_{CR}/L_\gamma$.
   • Parameters (bronafstand $d_s$, versnellings-efficiëntie $\eta$ en nucleaire fracties $\epsilon_Z$) worden gefit aan observationele gegevens van Auger en TA (flux en gemiddelde stortdiepte $\langle X_{max} \rangle$).
   • Statistische modelvergelijking wordt uitgevoerd met behulp van gereduceerde $\chi^2$, Akaike Informatiecriterium (AIC) en Bayesiaans Informatiecriterium (BIC).
   • Volledige-hemelsferische kaarten van integrale flux ($10-100$ EeV) en dipoolanisotropie worden gegenereerd met behulp van het HEALPix-raamwerk.

Belangrijkste Resultaten

1. Flux en Samenstelling:

   • Zowel FSRQ- als BL Lac-scenario's kunnen het waargenomen UHECR-energiespectrum en $\langle X_{max} \rangle$-gegevens reproduceren met gemengde samenstellingen.
   • De beste-fit versnellings-efficiëntie is $\eta \approx 0.1$, consistent met numerieke simulaties van relativistische jet-versnelling.
   • De vereiste baryonbeladingsfactor is $\kappa \approx 5.65$, wat aangeeft dat het CR-vermogen energetisch haalbaar is (slechts een paar keer de gammastraal-luminositeit).
   • FSRQ-scenario: Vereist een bronafstand $d_s \approx 8 \pm 1$ Mpc om de flux te passen. Dit impliceert een bronnedichtheid die aanzienlijk hoger is dan observationele schattingen uit Fermi-LAT-luminositeitsfuncties.
   • BL Lac-scenario: Vereist $d_s \approx 160 \pm 21$ Mpc, wat consistent is met observationele schattingen van BL Lac-dichtheid.

2. Statistische Vergelijking:

   • De gereduceerde $\chi^2$-waarden voor de gecombineerde Auger- en TA-gegevens zijn 2.37 (FSRQ) en 2.55 (BL Lac).
   • Modelselectiecriteria (AIC en BIC) geven de voorkeur aan het BL Lac-model ($\Delta$AIC = 4.17, $\Delta$BIC = 4.83), wat matig bewijs biedt dat BL Lacs een betere beschrijving van de gegevens bieden dan FSRQ's.

3. Anisotropiebeperkingen:

   • BL Lacs: De voorspelde dipoolanisotropie met de op flux gefitte parameters ($z_{max}=0.5$) stemt goed overeen met Auger-waarnemingen.
   • FSRQ's: De voorspelde dipoolamplitude is significant hoger dan waargenomen bij gebruik van dezelfde parameters. Om de waargenomen anisotropie te matchen, moet het roodverschuivingsbereik kunstmatig worden beperkt tot $z_{max} \approx 0.01$, wat in strijd is met de flux-fitvereisten en de bekende sterke kosmologische evolutie van FSRQ's.

4. Fotodisintegratie:

   • Berekeningen van optische diepte en geabsorbeerde fracties tonen aan dat zware kernen (bijv. Fe) aanzienlijke fotodisintegratie ondergaan bij lagere energieën in vergelijking met lichtere kernen (bijv. He), wat leidt tot een progressieve onderdrukking van zware componenten bij de hoogste energieën.

Betekenis en Conclusies
Het artikel concludeert dat hoewel zowel FSRQ's als BL Lacs in staat zijn de waargenomen UHECR-flux te reproduceren, het BL Lac-scenario sterk de voorkeur heeft wanneer de volledige set observationele beperkingen wordt overwogen.

• FSRQ's worden ongunstig beoordeeld als dominante bronnen van diffuse UHECR's, omdat het scenario een onrealistisch hoge bronnedichtheid (kleine $d_s$) vereist om de flux te passen en een dipoolanisotropie voorspelt die observationele limieten overschrijdt, tenzij de bronpopulatie kunstmatig wordt afgekapt bij lage roodverschuivingen.
• BL Lacs blijven consistent met waarnemingen, en bieden een zelfconsistent beschrijving van flux, samenstelling en dipoolanisotropie met een bronnedichtheid en roodverschuivings-evolutie die overeenkomen met Fermi-LAT-gegevens.

De auteurs suggereren dat hoewel FSRQ's mogelijk nog bijdragen aan multi-messenger kanalen (neutrino's en gammastralen), ze waarschijnlijk ondergeschikte bijdragers zijn aan de diffuse UHECR-flux in vergelijking met BL Lac-objecten. De studie benadrukt het belang van het gebruik van realistische, luminositeit-afhankelijke bron-evolutiemodellen in plaats van vereenvoudigde fenomenologische vormen bij het beperken van UHECR-oorsprongen.