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

Technische Zusammenfassung: Suche nach Neutrinoemissionen aus Blazar-$\gamma$-Ray-Flares unter Berücksichtigung möglicher Neutrino-Zeitverzögerungen

Problemstellung
Obwohl das IceCube-Neutrinoobservatorium die Existenz astrophysikalischer hochenergetischer Neutrinos ($E_\nu \gtrsim 1$ TeV) nachgewiesen hat, bleiben die spezifischen Quellen, die zum diffusen Fluss beitragen, weitgehend unidentifiziert, obwohl es Hinweise auf eine Verknüpfung bestimmter Ereignisse mit Blazaren (z. B. TXS 0506+056) und aktiven Galaxien (z. B. NGC 1068) gibt. Vorherige Suchen nach Korrelationen zwischen Blazaren und IceCube-Neutrinos gingen typischerweise von einer zeitlichen Koinzidenz zwischen der Ankunft des Neutrinos und bedeutenden elektromagnetischen (speziell $\gamma$-Strahlungs-) Flares aus. Theoretische Überlegungen legen jedoch nahe, dass diese Koinzidenz nicht halten muss. Wenn Neutrinos über $p\gamma$-Wechselwirkungen produziert werden, können die Zeitskalen für die Protonenbeschleunigung ($t'_{acc}$) und die Protonenenergieverluste ($t'_{p\gamma}$) im Jet-Frame um Größenordnungen länger sein als die für die Elektronen, die für elektromagnetische Flares verantwortlich sind. Infolgedessen können hochenergetische Neutrinos ($\gtrsim 100$ TeV) Monate bis Jahre nach dem Peak des zugehörigen $\gamma$-Strahlungsflares eintreffen. Diese Arbeit adressiert das Fehlen von Suchen, die eine solche potenzielle Zeitverzögerung berücksichtigen.

Methodik
Die Autoren führten eine Cross-Matching-Analyse zwischen hochenergetischen Neutrino-Alerts aus dem IceCat-1 öffentlichen Katalog und zwei unterschiedlichen Blazar-Quellensammlungen durch, wobei sie nach Korrelationen über einen weiten Bereich von Jet-Frame-Zeitverzögerungen ($t'_{delay}$) scannten.

1. Datensätze:

   • Neutrinos: 348 Alerts aus dem IceCat-1 Katalog ($E_\nu \gtrsim 100$ TeV), unter Ausschluss von Alerts mit koinzidenten IceTop-Signalen.
   • Quelllisten:
     • H21+: Eine Stichprobe von 294 AGNs aus der MOJAVE-Datenbank (Homan et al. 2021) mit bekannten Doppler-Faktoren ($D$), Rotverschiebungen ($z$) und Fermi-LAT-Lichtkurven.
     • R24+: Eine Stichprobe von 295 Blazaren aus dem CGRaBS-Katalog (Rodrigues et al. 2024b) mit durch SED-Modellierung abgeleiteten Doppler-Faktoren, ebenfalls kreuzweise mit Fermi-LAT-Lichtkurven abgeglichen.
   • Lichtkurven: Monatlich gebinnte Photonenfluss-Lichtkurven aus dem Fermi-LAT Light Curve Repository (LCR).

2. Statistischer Rahmen:

   • Räumliche Assoziation: Neutrinos wurden mit Blazaren assoziiert, wenn die Quelle innerhalb einer erweiterten Winkel-Unsicherheitsregion (unter Berücksichtigung einer systematischen Unsicherheit von $\Delta = 0,78^\circ$) lag.
   • Gewichtungsschema: Für jeden Alert wurde ein Gewicht $w_\nu(t'_{delay})$ berechnet, basierend auf der Signalstärke des Alerts, einem Straffaktor für große Winkelunsicherheiten und einem temporalen Faktor $G_{\nu a}$.
   • Temporaler Faktor: Der temporale Faktor verknüpfte die Ankunftszeit des Neutrinos ($t_\nu$) mit dem globalen Maximum der $\gamma$-Strahlungs-Lichtkurve des Blazars vor der Ankunft ($t^{\gamma \max}_a$). Der Faktor wies das maximale Gewicht zu, wenn die Zeitdifferenz innerhalb einer Bin-Größe der erwarteten Beobachter-Frame-Verzögerung ($t_{delay} = t'_{delay}(1+z)/D$) entsprach, wobei eine Gauß-Strafe für Abweichungen angewendet wurde.
   • Teststatistik (TS): Die globale TS wurde als Summe der individuellen Gewichte über alle Alerts definiert.
   • Signifikanzprüfung: Die TS der realen Daten wurde mit einer Verteilung von TS-Werten verglichen, die aus 49.999 Mock-Datensätzen generiert wurden, bei denen die Rektaszentren und Ankunftszeiten der Neutrinos randomisiert wurden. Pre-Trial-$p$-Werte wurden berechnet und anschließend für die Anzahl der Versuche (Scanning über $t'_{delay}$) korrigiert.

Wesentliche Beiträge

• Neuartige Suchstrategie: Diese Arbeit stellt den ersten Versuch dar, nach Assoziationen zwischen IceCube-Neutrinos und Blazar-$\gamma$-Strahlungs-Flares zu suchen, die spezifisch eine universelle Jet-Frame-Zeitverzögerung ($t'_{delay}$) aufgrund langsamer Protonenbeschleunigung oder Energieverluste berücksichtigen.
• Duale Quellensammlungen: Die Studie nutzt zwei unabhängige Blazar-Kataloge mit unterschiedlichen Methoden zur Bestimmung der Doppler-Faktoren (Temperatur-Helligkeits-Argumente vs. SED-Modellierung), um die Robustheit potenzieller Korrelationen gegenüber Unsicherheiten in der Jet-Kinematik zu testen.
• Systematisches Delay-Scanning: Die Analyse scannt $t'_{delay}$ von 1 Monat bis zu $\sim 3,5 \times 10^4$ Tagen und deckt damit den theoretischen Bereich ab, in dem die Neutrino-Produktion dem elektromagnetischen Ausstoß nachhinken könnte.

Ergebnisse

• Statistische Signifikanz: Die Suche ergab keine statistisch signifikante Evidenz für eine Korrelation zwischen IceCat-1 Neutrinos und den untersuchten AGN-Listen (H21+ oder R24+).
  • Für die H21+-Liste betrug der minimale Pre-Trial-$p$-Wert 0,037 (bei $t'_{delay} \approx 7,7 \times 10^3$ d).
  • Für die R24+-Liste betrug der minimale Pre-Trial-$p$-Wert 0,026 (bei $t'_{delay} \approx 1,9 \times 10^3$ d).
  • Nach der Trial-Korrektur entsprechen diese Minima Post-Trial-$p$-Werten von 0,14 (1,06$\sigma$) bzw. 0,11 (1,21$\sigma$), was mit der Nullhypothese konsistent ist.
• Top-Kandidaten: Die signifikantesten Assoziationen wurden durch die Alerts IC120523B (verknüpft mit dem 2010er Flare von 3C 454.3) und IC150926A (verknüpft mit einem Flare von 3C 279) getrieben. Obwohl diese spezifischen Ereignisse eine zeitliche Übereinstimmung mit den erwarteten Verzögerungen zeigten, reichten sie nicht aus, um eine globale Korrelation über den gesamten Datensatz zu etablieren.
• Diskrepanz in den Doppler-Faktoren: Die optimalen Verzögerungszeiten unterschieden sich signifikant zwischen den beiden Quelllisten (um einen Faktor 4), was primär auf die Diskrepanzen in den geschätzten Doppler-Faktoren für 3C 454.3 und 3C 279 zwischen den H21- und R24-Methodologien zurückzuführen ist.

Bedeutung und Behauptungen
Die Autoren kommen vorsichtig zu dem Schluss, dass ihre Ergebnisse keine signifikante Korrelation zwischen den untersuchten IceCat-1 Neutrinos und den spezifischen Blazar-Populationen stützen, selbst wenn man Zeitverzögerungen berücksichtigt. Sie führen mehrere, nicht notwendigerweise gegenseitig ausschließende Gründe für diese Nullergebnisse an:

1. Ineffizienz: Blazare könnten ineffiziente Quellen für $\sim 100$ TeV Neutrinos sein und nur einen kleinen Bruchteil ($\sim 0,5\%$) des IceCube-Flusses beitragen.
2. Energieskala: Falls Blazare effiziente Neutrinoquellen sind, könnte das Maximum ihrer Emission bei viel höheren Energien ($\gtrsim 10$ PeV) liegen, jenseits der primären Sensitivität der aktuellen IceCube-Alert-Stichprobe.
3. Modellannahmen: Die Annahme einer universellen Jet-Frame-Zeitverzögerung ($t'_{delay}$) könnte aufgrund von Variationen in den Beschleunigungsmechanismen und Magnetfeldern zwischen den einzelnen Quellen ungültig sein.
4. Entkopplung: Protonenbeschleunigung und Energiedissipation könnten räumlich oder zeitlich von Elektronenprozessen entkoppelt sein, was bedeutet, dass die Neutrino-Produktion nicht kausal mit den beobachteten $\gamma$-Strahlungs-Flares verknüpft ist.

Das Paper postuliert, dass zukünftige Verbesserungen der Sensitivität wahrscheinlich größere Datensätze (z. B. die kommenden IceCat-2 und CAZ Kataloge) sowie die Einbeziehung von Radio-Lichtkurven erfordern werden, welche bekanntlich $\gamma$-Strahlungs-Flares nachhinken und möglicherweise besser mit den Neutrino-Produktions-Zeitskalen korrelieren.