A Radiation-Hydrodynamic Light Curve Grid and Interpolation Framework for Stripped-Envelope Supernovae

Probleemstelling
Stripped-envelope supernovae (SESNe) bieden cruciaal inzicht in de evolutie van massieve sterren en de fysica van kerninstorting. Hoewel tijdreeks-surveys het aantal goed geobserveerde SESNe hebben vergroot, blijft het extraheren van betrouwbare fysische parameters (zoals ejectamassa, explosie-energie en nikkelmassa) uit lichtcurves een uitdaging. De standaardbenadering steunt op analytische Arnett-type modellen, die uitgaan van een constante opaciteit, centraal geconcentreerde verhitting en homologe expansie. De mate waarin deze simplificerende aannames leiden tot systematische biases of spuure correlations in de afgeleide parameters, is echter nog niet rigoureus gekwantificeerd. Daarentegen is de toepassing van stralingshydrodynamische modellen op observationele data gehinderd door computationele kosten en de discrete aard van simulatiegrids, ondanks dat zij een meer fysisch zelfconsistente beschrijving bieden.

Methodologie
De auteurs hebben een uitgebreid grid van 8.148 stralingshydrodynamische lichtcurve-modellen voor SESNe geconstrueerd.

• Progenitoren en Explosie: De evolutie van massieve sterren werd berekend met MESA (niet-roterend, zonne-metalliciteit) voor zero-age-main-sequence massa's ($M_{ZAMS}$) van 11,5 tot 20,0 $M_\odot$. Waterstofenveloppen werden kunstmatig verwijderd om kale heliumsterren te creëren, die werden geëvolueerd tot de uitputting van de kern-zuurstof. Explosies werden fenomenologisch getriggerd door thermische energie te deponeren in de binnenste 0,2 $M_\odot$, waarbij de asymptotische explosie-energie ($E_K$) en de ejectamassa ($M_{ej}$) varieerden.
• Radiatieve Overdracht: De hydrodynamische outputs werden verwerkt met de STELLA-code om bolometrische lichtcurves te berekenen. Het grid beslaat een bereik van $M_{ej}$ (gemotiveerd door late-fase spectrale beperkingen), specifieke kinetische energie ($\eta = E_K/M_{ej}$), radioactieve nikkelmassa ($M_{Ni}$) en een mengparameter ($f_{mix}$) die het fractie van de ejecta vertegenwoordigt die uniform gemengd is.
• Interpolatiekader: Om continue parameterinferentie mogelijk te maken, ontwikkelden de auteurs een vierdimensionaal lineair interpolatieschema. Dit kader schat lichtcurves in op willekeurige punten in de parameterruimte door 16 omliggende gridpunten te combineren, gewogen door hun relatieve afstanden.
• Validatie en Toepassing: De nauwkeurigheid van de interpolatie werd getest via een strikte cross-validatie (injectie-reconstructie methode), waarbij doelmodellen werden gereconstrueerd met behulp van uitsluitend niet-identieke gridpunten. Het kader werd vervolgens toegepast op twee geobserveerde SNe Ib (2007Y en 2009jf) met behulp van Markov Chain Monte Carlo (MCMC) analyse.
• Benchmarking van Arnett-modellen: Het numerieke grid diende als een gecontroleerd laboratorium om Arnett-type modellen te testen. Synthetische lichtcurves uit het grid werden gefit met het Arnett-model (met een vaste opaciteit $\kappa = 0,06$ cm$^2$ g$^{-1}$) om de afgeleide parameters te vergelijken met de bekende "ware" inputs.

Belangrijkste Bijdragen

1. Grootschalig Simulatiegrid: De publicatie van een grid van 8.148 SESN-lichtcurves die een representatief bereik van ejectamassa's, explosie-energieën, nikkelmassa's en menggraden beslaat.
2. Interpolatie-gebaseerde Inferentie: Een nieuw kader dat de kloof overbrugt tussen discrete numerieke simulaties en observationele data, waardoor een continue exploratie van de fysische parameterruimte mogelijk is met een hoge computationele efficiëntie.
3. Kwantificering van Systematische Bias: Een rigoureuze beoordeling van het Arnett-model met het numerieke grid als grondwaarheid, wat specifieke beperkingen in de parameterreconstructie aan het licht bracht.

Resultaten

• Afhankelijkheden van de Lichtcurve: De studie bevestigt dat $L_{peak}$, $t_{peak}$ en de vorm van de lichtcurve complexe functies zijn van $M_{ej}$, $\eta$, $M_{Ni}$ en $f_{mix}$. Opvallend genoeg verhoogt een toenemende $M_{Ni}$ alle tijdschalen, een gedrag dat niet wordt voorspeld door standaard Arnett-modellen. Sterke menging ($f_{mix}$) vermindert $L_{peak}$ en creëert schouderachtige kenmerken in de post-piek evolutie.
• Interpolatienauwkeurigheid: Het interpolatieschema is zeer accuraat, met mediane residuen onder de 0,01 dex en een conservatieve onzekachtschatting van 0,02 dex.
• Toepassing op SNe 2007Y en 2009jf:
  • Voor SN 2007Y reproduceerde het model succesvol de lichtcurve en kleur-evolutie, waarbij een ejectamassa van $\sim 1,86 M_\odot$ werd afgeleid, wat consistent is met zwakke nebulaire [O I] emissie.
  • Voor SN 2009jf onthulde de inferentie een sterke degeneratie tussen ejectamassa ($M_{ej}$) en menging ($f_{mix}$). Twee verschillende oplossingen (lage massa/hoge menging en hoge massa/lage menging) produceerden bijna identieke bolometrische lichtcurves, en de kleur-evolutie bood beperkte aanvullende beperkingen.
• Beperkingen van het Arnett-model:
  • Ejectamassa: Het Arnett-model faalt in het herstellen van de ware ejectamassa; afgeleide waarden vertonen geen significante correlatie met de ware waarden ($\rho = 0,19$) en clusteren in een smal bereik ($\sim 1,7–2,5 M_\odot$), ongeacht de input.
  • Spuure Correlaties: Hoewel $E_K$ en $M_{Ni}$ onafhankelijk zijn in het simulatiegrid, induceert de Arnett-fittingprocedure een spuure positieve correlatie ($M_{Ni} \propto E_K^{4/3}$). Dit komt doordat het Arnett-model uitgaat van een constante opaciteit, terwijl in werkelijkheid een hogere $M_{Ni}$ de ionisatie en opaciteit verhoogt, wat de lichtcurve verbredt. De Arnett-fit compenseert dit door een grotere $M_{ej}$ en $E_K$ af te leiden.

Significantie
Dit artikel demonstreert dat hoewel Arnett-type modellen efficiënte fenomenologische beschrijvingen van SESN-lichtcurves bieden, de fysische parameters die zij afleiden (met name de ejectamassa) vaak onbetrouwbaar zijn en niet de ware onderliggende fysica reflecteren. De schijnbare correlaties die worden waargenomen in populatiestudies (bijv. tussen explosie-energie en nikkelmassa) kunnen artefacten zijn van de inferentiemethode in plaats van intrinsieke fysische eigenschappen. Daarom pleiten de auteurs voor de noodzaak van fysisch gemotiveerde stralingshydrodynamische modellen in populatiestudies om misleidende conclusies over progenitor-evolutie en explosiemechanismen te voorkomen. Het ontwikkelde interpolatiekader biedt een praktische, accurate tool voor de toepassing van deze complexe modellen op grote observationele datasets.