Search for Long-Transient Gravitational Waves from Supernova SN2023ixf using GFH-v2 Pipeline

Technisches Resümee: Suche nach langanhaltenden transienten Gravitationswellen von SN 2023ixf mittels GFH-v2

Problem und Motivation
Es wird vorhergesagt, dass schnell rotierende, junge Magnetare, die bei Kernkollaps-Supernovae (CCSNe) entstehen, langanhaltende transiente kontinuierliche Gravitationswellen (tCWs) emittieren. Im Gegensatz zu kurzzeitigen Bursts, die mit dem Kollaps selbst einhergehen, entstehen diese Signale durch das Spin-down eines neu gebildeten, nicht-achsymmetrischen Neutronensterns oder Magnetars. Die Emission wird durch den Verlust von Rotationsenergie angetrieben, was ein Signal erzeugt, dessen Frequenz gemäß einem Potenzgesetz abnimmt und dessen Schwingungsamplitude über die Zeit abfällt. Der Nachweis dieser Signale würde direkte Einblicke in die Geburt, die interne Struktur und die Rotationsentwicklung von Neutronensternen ermöglichen.

SN 2023ixf, eine Typ-II-Kernkollaps-Supernova, die in der nahen Galaxie M101 (Entfernung $\approx$ 6,85 Mpc) entdeckt wurde, stellt eine einzigartige Gelegenheit für eine solche Suche dar. Ihre Nähe und die präzise zeitliche Bestimmung der elektromagnetischen Detektion (innerhalb eines Fensters von $\sim$2 Stunden) ermöglichen eine gerichtete Suche. Entscheidend ist, dass der Zeitpunkt des Ereignisses mit dem Engineering Run 15 (ER15) von LIGO überlappt, was einen Datensatz mit koinzidenten Daten aus den Hanford- (H1) und Livingston-Detektoren (L1) unmittelbar vor dem vierten Beobachtungsdurchlauf (O4) bereitstellt.

Methodik
Die Analyse verwendet die GFH-v2-Pipeline, eine hierarchische semikoherente Suchmethode, die auf der verallgemeinerten Frequenz-Hough-Transformation basiert und speziell für Signale mit Potenzgesetz-Frequenzentwicklung entwickelt wurde.

• Signalmodell: Die Suche zielt auf einen starren, nicht-präzedierenden Neutronenstern ab, der primär durch Gravitationswellen-Emission an Rotationsgeschwindigkeit verliert (Spin-down). Dies entspricht einem Bremsindex von $n=5$. Die Gravitationswellenfrequenz $f_{gw}$ entwickelt sich als $\dot{f}_{gw} = -k f_{gw}^n$, und die Schwingungsamplitude $h_0(t)$ zerfällt als $t^{-1/2}$. Die Suche nimmt kanonische Neutronensternparameter an ($M=1,4 M_\odot$, $R=12$ km) und nutzt eine feste Himmelsposition basierend auf der elektromagnetischen Lokalisierung.
• Daten: Die Analyse verwendet kozinzidente wissenschaftliche Segmente aus H1 und L1 während ER15, die insgesamt etwa 48.619 Sekunden (0,56 Tage) an Daten innerhalb eines konservativen On-Source-Fensters umfassen. Die Datenqualität wurde im Bereich von 600–2000 Hz als vergleichbar mit der O4a-Empfindlichkeit verifiziert.
• Suchkonfiguration: Der Parameterraum deckt Anfangsfrequenzen ($f_0$) von 600 bis 2000 Hz (Schrittweite 10 Hz) und Elliptizitäten ($\epsilon$) von $3 \times 10^{-4}$ bis $3 \times 10^{-3}$ (unterteilt in fünf Teilintervalle) ab. Dies ergibt 700 unabhängige Suchkonfigurationen.
• Pipeline-Schritte:
  1. Datenvorbereitung: Die Short Fourier Transform Database (SFDB)-Daten werden bandpassgefiltert und in komplexe analytische Zeitreihen transformiert.
  2. Peakmap-Generierung: Lokale Maxima im egalisierten Leistungsspektrum werden ausgewählt, um eine Zeit-Frequenz-Peakmap zu bilden.
  3. Hough-Transformation: Die Peakmap wird von $(t, f)$ in transformierte Koordinaten $(t, x)$ abgebildet, wobei $x = 1/f^{n-1}$ gilt. In diesem Raum wird die Potenzgesetz-Entwicklung linear. Eine Hough-Transformation bildet dies auf den $(x_0, k)$-Parameterraum (Anfangs-Frequenz und Spin-down-Parameter) ab.
  4. Kandidatenselektion: Kandidaten werden als Bins mit überschüssigen Zählungen relativ zum lokalen Hintergrund identifiziert. Ein kritisches Verhältnis (Critical Ratio, CR) wird unter Verwendung lokal geschätzter Rauschstatistiken (Median und Dispersion innerhalb von 10-Hz-Intervallen) berechnet, um nicht-stationäres Rauschen zu berücksichtigen.
  5. Koinzidenz: Kandidaten aus H1 und L1 werden im $(x_0, k)$-Parameterraum verglichen. Ein Paar gilt als koinzident, wenn die Distanzmetrik $d_{coin} \le 2$ ist.
     6. Follow-up: Koinzidente Kandidaten mit $CR \ge 3,5$ durchlaufen ein Heterodyn-Follow-up. Die erwartete Phasenentwicklung wird aus den Daten entfernt, und das Signal wird erneut mit längeren Kohärenzzeiten ($3\times$ und $6\times$ der ursprünglichen Zeit) analysiert, um zu prüfen, ob eine Zunahme der Signifikanz auftritt, was auf ein echtes astrophysikalisches Signal hindeuten würde.

Wichtigste Ergebnisse

• Kandidatenselektion: Nach Anwendung der Koinzidenzanforderung und des CR-Schwellenwerts ($CR \ge 3,5$) verblieben 131 koinzidente Kandidaten.
• Follow-up-Ergebnis: Keiner der 131 Kandidaten zeigte die erwartete Zunahme des kritischen Verhältnisses während der Heterodyn-Follow-up-Stufen. Folglich wurde kein Kandidat in die zweite Follow-up-Stufe befördert.
• Detektion: Die Analyse ergab keine Evidenz für langanhaltende transiente Gravitationswellen-Emission im Zusammenhang mit SN 2023ixf in den ER15-Daten.
• Obere Grenzwerte: Obere Grenzwerte für die maximal detektierbare Distanz wurden als Funktion der Anfangsfrequenz und Elliptizität festgelegt. Für das höchste Elliptizitätsintervall ($\epsilon \in [2,5 \times 10^{-3}, 3,0 \times 10^{-3}]$) liegen die 90%-Effizienz-Obergrenzen über den Großteil des analysierten Frequenzbandes bei etwa 1 bis 2,5 Mpc. Für geringere Elliptizitäten waren die Grenzwerte entweder nicht bestimmbar oder aufgrund der begrenzten Datendauer und niedrigerer Signal-Rausch-Verhältnisse eingeschränkt.

Bedeutung und Ansprüche
Die Autoren stellen fest, dass, obwohl die abgeleiteten oberen Grenzwerte (1–2,5 Mpc) unter der tatsächlichen Entfernung zu M101 ($\approx$ 6,85 Mpc) liegen, die Arbeit die erste Anwendung der GFH-v2-Pipeline auf eine gerichtete Suche nach einer nahen Kernkollaps-Supernova darstellt. Die Studie charakterisiert erfolgreich die Performance der Pipeline auf realen Detektordaten aus der Kommissionsierungsphase (ER15). Sie zeigt, dass die GFH-v2-Methode effektiv eingesetzt werden kann, um nach tCWs von Magnetar-bildenden Ereignissen zu suchen, und etabliert damit einen Rahmen für zukünftige Suchen unter Verwendung empfindlicherer Datensätze aus kommenden Beobachtungsläufen (O5) und Detektoren der nächsten Generation. Das Paper schließt mit der Feststellung, dass die Suche eine direkte Sensitivitätsschätzung für langanhaltende transiente Suchen in den ER15-Daten liefert und als Benchmark für zukünftige Bemühungen dient.