IceCube Real-time Searches for High-energy Neutrinos Coincident with LIGO/Virgo/KAGRA Gravitational-Wave Alerts in O4a

Technische Zusammenfassung: Echtzeit-Suche von IceCube nach hochenergetischen Neutrinos koinzident mit LIGO/Virgo/KAGRA-Gravitationswellen-Alarmen in O4a

Problem und Motivation
Kompakte Binärverschmelzungen (Compact Binary Coalescences, CBC), einschließlich binärer Neutronensterne (BNS), Neutronenstern–Schwarzes-Loch-Merger (NSBH) und binärer Schwarzer-Loch-Merger (BBH), sind vorhergesagte Quellen hochenergetischer (TeV–PeV) Neutrinos. Es wird erwartet, dass diese Ereignisse relativistische Jets erzeugen, die in der Lage sind, kosmische Strahlung zu beschleunigen, welche wiederum durch Wechselwirkungen Pionen erzeugt und in Neutrinos zerfällt. Während die Detektion von GW170817 und seinen elektromagnetischen Gegenstücken die Lebensfähigkeit der Multi-Messenger-Astronomie etabliert hat, bleiben die Quellen der Mehrheit der astrophysikalischen Neutrinos unidentifiziert. IceCube hat in den Beobachtungsperioden O1, O2 und O3 bereits nach Neutrino-Gegenstücken zu Gravitationswellen-Ereignissen (GW) gesucht, konnte jedoch keine statistisch signifikante Emission finden. Diese Arbeit berichtet über Echtzeit-Suchen, die während des ersten Teils der vierten Beobachtungsperiode (O4a) des LIGO-Virgo-KAGRA-Netzwerks (LVK) durchgeführt wurden, mit dem Ziel, koinzidente hochenergetische Neutrinos zu identifizieren und die Latenz der Reaktions-Pipeline zu verbessern.

Methodik
Die Analyse nutzt Daten des IceCube Neutrino Observatory, spezifisch den „Gamma-ray Follow-up“ (GFU)-Datensatz, der aus hochenergetischen, track-ähnlichen Ereignissen besteht, die mit geringer Latenz rekonstruiert wurden. Die Studie umfasst 85 signifikante und 945 Ereignisse mit geringer Signifikanz von CBC-Kandidaten, die während O4a und der vorangegangenen Engineering-Run (ER15) in Echtzeit von LVK veröffentlicht wurden. Es werden zwei unterschiedliche Analysemethoden angewandt:

1. Unbinned Maximum Likelihood (UML) Analyse: Eine generische Transienten-Suche, die für signifikante GW-Alarme verwendet wird. Diese Methode konstruiert eine Likelihood-Funktion, die die Signal- und Hintergrund-Wahrscheinlichkeitsdichtefunktionen (PDFs) über die beobachteten Neutrino-Ereignisse summiert. Die Signal-PDF nimmt eine Punktquelle mit einem Potenzgesetz-Energiespektrum ($E^{-\gamma}$) an und bezieht die GW-Skymap-Wahrscheinlichkeit ($w(\Omega)$) ein. Die Teststatistik (TS) wird über den Ort der Quelle, die Signalanzahl und den Spektralindex maximiert. P-Werte werden abgeleitet, indem die beobachtete TS gegen eine hintergrund-dominierte Verteilung verglichen wird, die durch Pseudo-Experimente generiert wurde, bei denen die Ereignisse in der Rektaszension und Zeit gestreut wurden.
2. Low-Latency Algorithm for Multimessenger Astrophysics (LLAMA): Ein modellabhängiger Bayesscher Ansatz, der sowohl auf signifikante als auch auf Ereignisse mit geringer Signifikanz angewendet wird. Diese Pipeline bezieht astrophysikalische Priors bezüglich der Verteilung der GW-Emissionsenergien, der Wahrscheinlichkeit, ob Ereignisse astrophysikalisch oder Rauschen sind, sowie der isotropen äquivalenten Energie der Neutrinoemission ein. Die Teststatistik gewichtet die Wahrscheinlichkeit von vier Hypothesen: echte Koinzidenz, Zufallskoinzidenz mit einem astrophysikalischen GW, Zufallskoinzidenz mit einem astrophysikalischen Neutrino und reiner Hintergrund.

Das primäre Suchfenster ist auf $\pm 500$ Sekunden relativ zur Verschmelzungszeit festgelegt, was der erwarteten Dauer kurzer Gammablitze entspricht. Zusätzlich wurde ein erweitertes Suchfenster von $[-0,1, +14]$ Tagen auf drei spezifische Kandidaten angewendet: zwei NSBH-Merger und ein Ereignis mit geringer Signifikanz mit einem potenziellen unterschwelligen Gamma-Ray-Gegenstück.

Wesentliche Beiträge und operative Verbesserungen
Ein bedeutender Beitrag dieser Arbeit ist die Verbesserung der Echtzeit-Infrastruktur von IceCube. Die Autoren berichten über die Implementierung eines automatisierten Systems, das die Ergebnisse sowohl der UML- als auch der LLAMA-Pipelines aufnimmt, sie zusammenführt und automatisch maschinenlesbare GCN-Notices versendet. Diese Automatisierung reduzierte die mediane Latenz der IceCube-Reaktion auf GW-Ereignisse im Vergleich zur vorherigen Beobachtungsperiode (O3) um den Faktor 2,7 und erreichte eine mediane Gesamtlatenz von 21 Minuten. Das System gibt zudem automatisch koinzidente Neutrino-Informationen (Richtung, Fehler, Timing) frei, falls die Pre-Trial-P-Werte unter 0,1 fallen, und versendet GCN Circulars bei P-Werten unter 0,01.

Ergebnisse

• Kurzdauer-Suche ($\pm 500$ s): Es wurde keine statistisch signifikante Neutrinoemission von irgendeinem der 85 signifikanten oder 945 niedrig-signifikanten GW-Kandidaten detektiert, die analysiert wurden. Die Verteilung der beobachteten P-Werte für beide UML- und LLAMA-Pipelines ist konsistent mit den Erwartungen eines reinen Hintergrunds. Der niedrigste beobachtete Pre-Trial-P-Wert lag bei 0,0011 für einen Kandidaten mit geringer Signifikanz (S231025a), welcher nach Berücksichtigung der Trials über den gesamten Datensatz hinweg keinen Hinweis auf eine wahre Koinzidenz zeigte.
• Langdauer-Suche: Die erweiterte Suche für die drei spezifischen Kandidaten (zwei NSBH-Merger und ein Ereignis mit geringer Signifikanz) ergab ebenfalls keine signifikante Emission.
• Obere Grenzwerte: Für alle analysierten Ereignisse setzen die Autoren 90%-Konfidenzniveaus-Obergrenzen für den zeitintegrierten Neutrinofluss und die isotrope äquivalente Energie ($E_{iso}$) fest, die in hochenergetischen Neutrinos emittiert wurde. Diese Grenzwerte nehmen eine isotrope Emission und ein Potenzgesetz-Spektrum an. Die $E_{iso}$-Grenzwerte werden über den Energiebereich von 500 GeV bis 50 PeV berechnet.
• Spezifisches Ereignis-Follow-up: Das Ereignis S230904n, ein wahrscheinlicher BBH-Merger, wurde in Echtzeit nachverfolgt. Während ein potenzieller Gegenpart (AT2023rkw) durch das Zwicky Transient Facility identifiziert wurde, identifizierte die anschließende Spektroskopie diesen als Typ Ia Supernova, die nicht mit dem GW-Ereignis zusammenhängt.

Bedeutung und Ansprüche
Das Paper behauptet, dass die Hinzunahme der 85 O4a-Ereignisse die Anzahl der mit diesen Suchen analysierten GW-Kandidaten im Vergleich zu früheren Durchläufen nahezu verdoppelt hat. Die primäre Bedeutung liegt in der Demonstration eines robusten, niederlatenz-basierten Echtzeit-Analyse-Frameworks, das sowohl signifikante als auch Ereignisse mit geringer Signifikanz verarbeiten kann. Die Autoren betonen, dass die Fähigkeit, diese Suchen in Echtzeit durchzuführen, obwohl keine Neutrino-Gegenteile gefunden wurden, eine sofortige Nachverfolgung durch andere Observatorien ermöglicht, was potenziell die Signifikanz zukünftiger Kandidaten erhöht und die Lokalisierung präzisiert. Die Ergebnisse begrenzen zudem theoretische Modelle der Neutrinoerzeugung in kompakten Binär-Mergern, indem sie strenge obere Grenzwerte für die in hochenergetischen Neutrinos emittierte Energie festlegen. Die Autoren merken an, dass zukünftige Verbesserungen der GW-Detektor-Empfindlichkeit und die geplante Erweiterung zu IceCube-Gen2 die Aussichten auf solche Koinzidenzen erhöhen werden.