Helium superluminous SN 2021bnw : an explosion of a massive star with a pre-outburst

Technische Zusammenfassung: Helium-Superluminöse Supernova 2021bnw

Problemstellung
Superluminöse Supernovae (SLSNe) sind transiente Ereignisse mit peak-Absoluthelligkeiten, die −20 mag überschreiten. Während viele durch den Spin-down von Magnetaren oder radioaktiven Zerfall angetrieben werden, bleiben die spezifischen Mechanismen, die heliumreiche SLSNe (He-SLSNe) antreiben, umstritten. SN 2021bnw ist eine nahe He-SLSN, die durch eine breite Lichtkurve (LC) von etwa 140 Tagen und eine Peak-Helligkeit von −20,5 mag gekennzeichnet ist. Eine reine radioaktive Antriebskraft durch den Zerfall von $^{56}$Ni würde eine unphysikalisch große Nickelmasse (~14 $M_\odot$) erfordern, was auf einen Ursprung als Paarinstabilitäts-Supernova (PISN) hindeutet. Die beobachtete LC fällt jedoch nach dem Peak schnell ab und fehlt der langsam evolvierende, hochleuchtende Schwanz, der für PISNe charakteristisch ist. Umgekehrt wird ein Ursprung als pulsierende Paarinstabilitäts-Supernova (PPISN) durch die spezifischen Energie- und Masseneinschränkungen, die aus den Daten abgeleitet wurden, in Frage gestellt. Die Studie zielt darauf ab, den Antriebsmechanismus und die Natur des Vorläufersterns von SN 2021bnw durch die Modellierung ihrer pseudo-bolometrischen Lichtkurve zu bestimmen.

Methodik
Die Autoren verwendeten hydrodynamische Strahlungstransport-Simulationen mit dem Code STELLA, um die Wechselwirkung zwischen Supernova-Ejekta und zirkumstellarem Material (CSM) zu modellieren, ergänzt durch radioaktive Heizung.

• Ejekta-Modelle: Die Ejekta-Profile wurden auf eine hydrodynamische Vorlage eines 90 $M_\odot$ schweren Heliumsterns skaliert (He90-Modell von Heger & Woosley 2002). Die Autoren testeten Ejekta-Massen im Bereich von 15 bis 90 $M_\odot$. Die Zusammensetzung bestand hauptsächlich aus Sauerstoff, Helium, Kohlenstoff, Silizium und Schwefel.
• CSM-Modelle: Die Ejekta waren von einem windähnlichen CSM ($\rho \propto r^{-2}$) mit heliumreicher Zusammensetzung (He:C:O-Verhältnis von 5:3:2) umgeben, motiviert durch NIR-Spektrallinien, die 87 Tage nach dem Peak beobachtet wurden. Die CSM-Massen variierten zwischen 1 und 22,2 $M_\odot$.
• Energiequellen: Die Simulationen beinhalteten kinetische Energie aus der Explosion, Wechselwirkungsenergie zwischen Ejekta und CSM sowie radioaktive Heizung durch den Zerfall von $^{56}$Ni und $^{56}$Co. Während der Hauptexplosionsphase wurde keine zusätzliche Energieinjektion simuliert.
• Frühe Epochen-Analyse: Um frühe Datenpunkte (ca. 38 Tage vor dem Peak) zu erklären, die die Leuchtkraft des Vorläufersterns überschritten, testeten die Autoren Modelle, die ein kompaktes CSM und eine expandierende Hülle beinhalten, die einen prä-supernova-Ausbruch darstellt.

Hauptergebnisse

1. Best-Fit-Explosionsparameter: Die beobachtete Lichtkurve wird am besten durch eine Kombination aus CSM-Wechselwirkung und radioaktivem Zerfall reproduziert. Die optimalen Modelle erfordern:

   • Ejekta-Masse: 15–22,5 $M_\odot$.
   • CSM-Masse: ~7 $M_\odot$.
   • Kinetische Energie: 4–5 foe ($10^{51}$ erg).
   • Radioaktives Nickel: 1,7 $M_\odot$ $^{56}$Ni.
   • Wechselwirkung: Die Kollision der 15–22,5 $M_\odot$ schweren Ejekta mit der 7 $M_\odot$ schweren CSM treibt die ersten ~80 Tage der Lichtkurve an. Die 1,7 $M_\odot$ $^{56}$Ni sind erforderlich, um die Leuchtkraft aufrechtzuerhalten und die Lichtkurve bis zum Tag 140 zu verlängern.

2. Anomalie der frühen Lichtkurve: Die initialen Datenpunkte (vor dem Peak) können nicht durch die Haupt-Ejekta-CSM-Wechselwirkungsmodelle erklärt werden. Die Autoren schlagen vor, dass diese Punkte einen prä-Ausbruch darstellen, der etwa 18 Tage vor der terminalen Explosion stattfand. Dies wird als eine expandierende Hülle mit einer Masse von 0,5–1 $M_\odot$ und einer kinetischen Energie von 0,7–1 foe modelliert, die einer adiabatischen Abkühlung unterliegt.

3. Vorläufer-Einschränkungen:

   • Die Gesamtmasse aus Ejekta plus CSM wird auf 22–30 $M_\odot$ geschätzt.
   • Unter Verwendung von Sternentwicklungsbeziehungen (Heger & Woosley 2002) und unter Berücksichtigung des Massenverlusts wird die initiale Null-Alter-Hauptreihen-Masse (ZAMS) des Vorläufers auf nicht weniger als 61 $M_\odot$ geschätzt (und wahrscheinlich unter 120 $M_\odot$).
   • Die Explosionsenergie (4–5 foe) und die hohe $^{56}$Ni-Ausbeute (1,7 $M_\odot$) sind konsistent mit einem magnetorotationellen Kernkollaps (MRCC)-Szenario, das potenziell einen Kollapsar-Mechanismus beinhaltet, anstatt eines standardmäßigen neutrino-getriebenen Kollapses.

Ausschluss von PPISN
Die Autoren argumentieren gegen einen PPISN-Ursprung für SN 2021bnw basierend auf mehreren Faktoren:

• Massenlimit: Die abgeleitete Gesamtmasse (Ejekta + CSM) liegt unter dem typischen unteren Limit für PPISN-Vorläufer (34–48,5 $M_\odot$ Heliumkerne).
• Restkörper-Physik: Standard-PPISN-Modelle im relevanten Massenbereich hinterlassen typischerweise ein massives Schwarzes Loch (30–40 $M_\odot$) ohne signifikante Massenauswurf oder $^{56}$Ni-Produktion.
• Energieanforderungen: Die Erreichung der erforderlichen kinetischen Energie (4–5 foe) und Ejekta-Masse in einem realistischen Kernkollaps-Szenario mit einem massiven Restkörper ist für Standard-PPISN-Modelle problematisch.
• Nickel-Produktion: Die Notwendigkeit von 1,7 $M_\odot$ $^{56}$Ni widerspricht PPISN-Szenarien, die im Allgemeinen keine derart großen Mengen an radioaktivem Nickel produzieren.

Bedeutung und Behauptungen
Die Studie schließt, dass SN 2021bnw höchstwahrscheinlich das Ergebnis einer Kernkollaps-Explosion eines massereichen Sterns (ZAMS $\ge$ 61 $M_\odot$) ist, unterstützt durch magnetorotationelle Effekte. Das Ereignis wurde durch die Wechselwirkung der Ejekta mit einer massiven, heliumreichen CSM-Hülle und einer beträchtlichen Menge synthetisiertem $^{56}$Ni angetrieben. Die Studie hebt hervor, dass, obwohl aktuelle MRCC-Simulationen typischerweise niedrigere $^{56}$Ni-Massen (0,5–1,14 $M_\odot$) ergeben, zukünftige hochauflösende 3D-Simulationen und verbesserte Neutrino-Transport-Schemata diese Lücke schließen könnten. Die Autoren stellen ausdrücklich fest, dass, obwohl das MRCC-Szenario die plausibelste Erklärung ist, die spezifischen Dynamiken, die erforderlich sind, um 1,7 $M_\odot$ $^{56}$Ni in einem so massereichen Vorläufer zu produzieren, eine offene Frage für zukünftige Studien zur explosiven Nukleosynthese bleiben. Die Studie behauptet nicht, den MRCC-Mechanismus endgültig zu beweisen, sondern präsentiert ihn als bevorzugtes Modell gegenüber PISN- oder PPISN-Alternativen angesichts der aktuellen Daten und Modellierungseinschränkungen.