Search for neutrino emission from blazar $γ$-ray flares accounting for possible neutrino time delays

Technische Samenvatting: Zoektocht naar neutrino-emissie van blazar $\gamma$-straalflares rekening houdend met mogelijke neutrino-tijdvertragingen

Probleemstelling
Hoewel de IceCube-neutrino-observator het bestaan van astrofysische hoogenergetische neutrino's ($E_\nu \gtrsim 1$ TeV) heeft vastgesteld, blijven de specifieke bronnen die bijdragen aan de diffuse flux grotendeels onbekend, ondanks bewijs dat specifieke gebeurtenissen koppelt aan blazars (bijv. TXS 0506+056) en actieve sterrenstelsels (bijv. NGC 1068). Eerdere zoektochten naar correlaties tussen blazars en IceCube-neutrino's vertrouwden doorgaans op temporele samenval tussen de aankomst van het neutrino en belangrijke elektromagnetische (specifiek $\gamma$-straal) flares. Theoretische overwegingen suggereren echter dat deze samenval niet standhoudt. Als neutrino's worden geproduceerd via $p\gamma$-interacties, kunnen de tijdschalen voor protonversnelling ($t'_{acc}$) en protonenergieverliezen ($t'_{p\gamma}$) in het jet-frame orders van grootte langer zijn dan die van de elektronen die verantwoordelijk zijn voor de elektromagnetische flares. Gevolgelijk kunnen hoogenergetische neutrino's ($\gtrsim 100$ TeV) maanden tot jaren na de piek van de geassocieerde $\gamma$-straalflare arriveren. Dit artikel adresseert het gebrek aan zoektochten die rekening houden met dergelijke potentiële tijdvertragingen.

Methodologie
De auteurs voerden een cross-matching analyse uit tussen hoogenergetische neutrino-alerts uit de IceCat-1 publieke catalogus en twee verschillende blazar-bronlijsten, waarbij zij zochten naar correlaties over een breed bereik van jet-frame tijdvertragingen ($t'_{delay}$).

1. Datasets:

   • Neutrino's: 348 alerts uit de IceCat-1 catalogus ($E_\nu \gtrsim 100$ TeV), waarbij alerts met gelijktijdige IceTop-signalen zijn uitgesloten.
   • Bronlijsten:
     • H21+: Een steekproef van 294 AGN's uit de MOJAVE-database (Homan et al. 2021) met bekende Dopplerfactoren ($D$), roodverschuivingen ($z$) en Fermi-LAT lichtcurves.
     • R24+: Een steekproef van 295 blazars uit de CGRaBS-catalogus (Rodrigues et al. 2024b) met Dopplerfactoren afgeleid van spectrale energieverdeling-modellering, eveneens gekoppeld aan Fermi-LAT lichtcurves.
   • Lichtcurves: Maandelijks gebindeerde fotonflux-lichtcurves uit de Fermi-LAT Light Curve Repository (LCR).

2. Statistisch Kader:

   • Ruimtelijke Associatie: Neutrino's werden geassocieerd met blazars als de bron binnen een vergrote hoekonzekerheidsregio viel (inclusief een systematische onzekerheid $\Delta = 0,78^\circ$).
   • Wegingsschema: Voor elke alert werd een gewicht $w_\nu(t'_{delay})$ berekend op basis van de signaalsterkte van de alert, een straffactor voor grote hoekonzekerheden en een temporele factor $G_{\nu a}$.
   • Temporele Factor: De temporele factor koppelde de aankomsttijd van het neutrino ($t_\nu$) aan het globale maximum van de $\gamma$-straallichtcurve van de blazar voorafgaand aan de aankomst ($t^{\gamma \max}_a$). De factor kende een maximaal gewicht toe als het tijdsverschil overeenkwam met de verwachte waarnemer-frame vertraging ($t_{delay} = t'_{delay}(1+z)/D$) binnen een bin-grootte, waarbij een Gaussische straf werd toegepast bij afwijkingen.
   • Teststatistiek (TS): De globale TS werd gedefinieerd als de som van individuele gewichten over alle alerts.
   • Significantietoetsing: De TS voor de echte data werd vergeleken met een distributie van TS-waarden gegenereerd uit 49.999 mock-datasets waarbij de rechte ascensie en aankomsttijden van de neutrino's werden gerandomiseerd. Pre-trial $p$-waarden werden berekend en vervolgens gecorrigeerd voor trials (het scannen over $t'_{delay}$).

Belangrijkste Bijdragen

• Nieuwe Zoekstrategie: Dit werk vormt de eerste poging om te zoeken naar associaties tussen IceCube-neutrino's en blazar $\gamma$-straalflares door specifiek rekening te houden met een universele jet-frame tijdvertraging ($t'_{delay}$) veroorzaakt door trage protonversnelling of energieverliezen.
• Duo Bronselecties: De studie maakt gebruik van twee onafhankelijke blazar-catalogi met verschillende methoden voor het bepalen van Dopplerfactoren (helderheidstemperatuur-argumenten versus SED-modellering) om de robuustheid van potentiële correlaties te testen tegen onzekerheden in jet-kinematica.
• Systematisch Vertragingsscannen: De analyse scant $t'_{delay}$ van 1 maand tot $\sim 3,5 \times 10^4$ dagen, wat het theoretische bereik dekt waar neutrino-productie mogelijk achterloopt op de elektromagnetische emissie.

Resultaten

• Statistische Significantie: De zoektocht leverde geen statistisch significante bewijslast voor een correlatie tussen IceCat-1 neutrino's en de onderzochte AGN's in de H21+ of R24+ bronlijsten.
  • Voor de H21+ lijst was de minimale pre-trial $p$-waarde 0,037 (bij $t'_{delay} \approx 7,7 \times 10^3$ d).
  • Voor de R24+ lijst was de minimale pre-trial $p$-waarde 0,026 (bij $t'_{delay} \approx 1,9 \times 10^3$ d).
  • Na correctie voor trials kwamen deze minima overeen met post-trial $p$-waarden van respectievelijk 0,14 (1,06$\sigma$) en 0,11 (1,21$\sigma$), wat consistent is met de nulhypothese.
• Top Kandidaten: De meest significante associaties werden gedreven door de alert IC120523B (gelinkt aan de 2010 flare van 3C 454.3) en IC150926A (gelinkt aan een flare van 3C 279). Hoewel deze specifieke gebeurtenissen een temporele uitlijning vertoonden met de verwachte vertragingen, waren ze onvoldoende om een globale correlatie over de gehele dataset vast te stellen.
• Discrepantie in Dopplerfactoren: De optimale vertagingstijden verschilden aanzienlijk tussen de twee bronlijsten (een factor 4), grotendeels gedreven door discrepanties in de geschatte Dopplerfactoren voor 3C 454.3 en 3C 279 tussen de H21 en R24 methodologieën.

Betekenis en Claims
De auteurs concluderen bescheiden dat hun resultaten geen significante correlatie ondersteunen tussen de specifieke onderzochte blazarpopulaties en IceCat-1 neutrino's, zelfs niet wanneer rekening wordt gehouden met tijdvertragingen. Zij suggereren verschillende, niet-onderling uitsluitende redenen voor deze nulresultaten:

1. Inefficiëntie: Blazars kunnen inefficiënte bronnen van $\sim 100$ TeV neutrino's zijn, die slechts een klein deel ($\sim 0,5\%$) van de IceCube-flux bijdragen.
2. Energieschaal: Als blazars efficiënte neutrino-bronnen zijn, kan het piek van hun emissie bij veel hogere energieën liggen ($\gtrsim 10$ PeV), buiten de primaire gevoeligheid van de huidige IceCube-alert sample.
3. Model-aannames: De aanname van een universele jet-frame tijdvertraging ($t'_{delay}$) kan ongeldig zijn vanwege variaties tussen bronnen in versnellingsmechanismen en magnetische velden.
4. Ontkoppeling: Protonversnelling en energiedissipatie kunnen ruimtelijk of temporeel ontkoppeld zijn van elektronprocessen, wat betekent dat neutrino-productie niet causaal verbonden is met de waargenomen $\gamma$-straalflares.

Het artikel stelt dat toekomstige verbeteringen in gevoeligheid waarschijnlijk grotere datasets zullen vereisen (bijv. de komende IceCat-2 en CAZ catalogi) en de inclusie van radio-lichtcurves, die bekend staan om het achterlopen op $\gamma$-straalflares en mogelijk beter correleren met de tijdschalen van neutrino-productie.