Signatures of $^{56}$Ni Mixing and Neutron-rich Ejecta in Supernovae

Technische Samenvatting: Signaturen van $^{56}\text{Ni}$-menging en Neutronenrijke Ejecta in Supernovae

Probleemstelling
Lichtcurves van supernova's (SN) worden routinematig geïnterpreteerd met behulp van one-zone radioactieve vervalmodellen (bijv. Arnett 1980, 1982) om eigenschappen van de progenitor af te leiden, zoals de ejectamassa ($M_{\text{ej}}$), nikkelmassa ($M_{\text{Ni}}$) en explosie-energie ($E_k$). Deze modellen gaan uit van één effectieve diffusietijd en negeren de ruimtelijke verdeling van radioactief materiaal. Hydrodynamische simulaties geven echter aan dat $^{56}\text{Ni}$ diep in de ejecta wordt gesynthetiseerd en vervolgens naar buiten toe wordt gemengd via Rayleigh–Taylor-instabiliteiten, convectie en turbulentie. Deze ruimtelijke verdeling is een cruciale fysische variabele die one-zone modellen weglaat. Bijgevolg kunnen lichtcurves van sterk gemengde ejecta verkeerd worden geïnterpreteerd als een lage ejectamassa of het vereisen van extra energie van een centrale motor (bijv. magnetars of fallback-accretie) om snelle stijgingen en hoge luminositeiten te verklaren. Bovendien hangt de observeerbaarheid van r-procesverrijking in collapsars gevoelig af van de plaatsing en geometrie van neutronenrijke ejecta ten opzichte van de door nikkel aangedreven achtergrond, een factor die in sferische modellen vaak te eenvoudig wordt weergegeven.

Methodologie
De auteurs introduceren snmix, een semi-analytisch, multi-shell model geïmplementeerd als een JAX-gebaseerde plugin voor het Redback-framework. In tegenstelling tot one-zone modellen discretiseert snmix homogeen expanderende ejecta in $N$ schillen met geometrisch gespatieerde snelheden.

• Radioactieve Structuur: Het model staat willekeurige radiale verdelingen van $^{56}\text{Ni}$ toe (geparameteriseerd door een mengingsfractie $f_{\text{mix}}$ in een top-hat of core-plus-tail profiel) en r-procesmateriaal.
• Transportfysica: Elke schil evolueert onafhankelijk en houdt thermische energie, adiabatische verliezen en radiatieve diffusie bij. Cruciaal is dat het model diffusie en gamma-straal depositie berekent op basis van de cumulatieve kolomdichtheid buiten elke specifieke schil, in plaats van een globaal gemiddelde.
• Opaciteit en SED: Het model incorporeert grijze opaciteiten en een temperatuurafhankelijke opaciteitswet. Voor r-procesberekeningen bevat het lanthanide-afhankelijke opaciteiten en analytische verwarmingssnelheden gekalibreerd op basis van nucleaire netwerkberekeningen.
• Geometrie: Hoewel primair sferisch, breidt het model zich uit naar een twee-componenten "disk-wind" geometrie voor collapsars, waarbij een polaire SN-component wordt gecombineerd met een equatoriale neutronenrijke wind, gewogen door kijkhoek.
• Inference: De auteurs gebruiken dit framework om synthetische bolometrische lichtcurves te genereren en voeren Bayesiaanse inferentie uit (met behulp van NUTS-sampling) om one-zone fits te vergelijken met de multi-shell realiteit, waarbij zij parameter-biases analyseren.

Belangrijkste Bijdragen en Resultaten

1. Impact van $^{56}\text{Ni}$-menging op Lichtcurves:

   • Naar buiten gerichte menging van $^{56}\text{Ni}$ vermindert de bovenliggende diffusiekolom en de gamma-straal trappingkolom. Dit resulteert in snellere, helderdere stijgingen en steilere radioactieve staarten voor een vaste $M_{\text{Ni}}$, $M_{\text{ej}}$ en $E_k$.
   • Een snelle, heldere stijging is geen exclusief bewijs voor een lage ejectamassa of een motor; sterk gemengde radioactieve modellen kunnen deze kenmerken reproduceren zonder extra energiebronnen.

2. Bias in One-Zone Inference:

   • Het fitten van gemengde lichtcurves met one-zone modellen produceert visueel goede fits, maar leidt tot bevooroordeelde fysische parameters.
   • Bij een vaste opaciteit absorbeert het model de ontbrekende radiale structuur door een aanzienlijk lagere $M_{\text{ej}}$ en een hogere nikkelmassafractie ($f_{\text{Ni}}$) te infereren, terwijl $M_{\text{Ni}}$ relatief stabiel blijft (hersteld binnen $\sim3\%$ met vaste opaciteit).
   • Als de opaciteit variabel wordt toegestaan, verschuift de bias naar het infereren van een veel lagere effectieve opaciteit ($\kappa_{\text{fit}}/\kappa_{\text{input}} \simeq 0.24$ voor volledig gemengde gevallen).
   • Deze biases propageren naar onjuiste conclusies over explosie-energieën en progenitor-mapping.

3. Casestudy: SN 1998bw:

   • Het toepassen van snmix op SN 1998bw levert een matige mengingsfractie op ($f_{\text{mix}} \approx 0.59$). Hoewel de one-zone maximum likelihood estimate overlapt met de posterior, wijkt deze met $\sim13\%$ af in $M_{\text{ej}}$ en $\sim5\%$ in $E_k$, wat illustreert dat bolometrische data alleen niet uniek massa, energie, opaciteit en menging kunnen ontleden.

4. Collapsar r-Proces Signaturen:

   • De signatuur van r-procesverrijking is niet universeel een late-time Near-Infrared (NIR) excess.
   • Plaatsing: Begraven r-procesmateriaal onderdrukt aanvankelijk het optische licht door een hoge opaciteit, waarbij het NIR-excess pas verschijnt nadat de verrijkte lagen transparant worden.
   • Geometrie: In disk-wind modellen, als neutronenrijk materiaal beperkt is tot het equatoriale vlak, vertonen on-axis GRB-SNe (polaire kijkhoeken) een verwaarloosbare r-proces signatuur omdat de polaire SN-component domineert. Het signaal is het sterkst voor off-axis waarnemers of als het materiaal gemengd is naar hogere breedtegraden ($\gtrsim 30^\circ$).
   • Diagnose: De kleur evolutie ($\Delta(r-H)$) is een waardevolle diagnose, maar vereist het ontrafelen van optische onderdrukking van H-band verheldering. Gepaarde optische en NIR monitoring (10–60 dagen) is informatiever dan enkel late-time H-band meetpunten.

Significantie
Dit artikel betoogt dat de "radioactieve structuur" (radiale verdeling van verwarming en opaciteit) een cruciale, vaak verwaarloosde variabele is in de interpretatie van supernova's.

• Voor Nikkel-aangedreven SNe: Het daagt de routinematige aanroeping van centrale motoren uit om snelle stijgende Ic-BL SNe en GRB-SNe te verklaren, en suggereert dat outward mixing een levensvatbaar alternatief is. Het waarschuwt dat one-zone fits "goede" statistische fits kunnen geven terwijl ze fysisch incorrecte parameters leveren.
• Voor Collapsars: Het verfijnt de zoekstrategie naar r-procesverrijking, waarbij verder wordt gekeken dan een eenvoudige "NIR excess" heuristiek. Het benadrukt dat non-detecties in on-axis GRB-SNe niet betekenen dat r-procesproductie is uitgesloten als de ejecta equatoriaal geconcentreerd zijn.
• Observationele Strategie: De auteurs concluderen dat het doorbreken van deze degeneraties vereist dat waarnemingen worden gecombineerd: fotometrische snelheden, kleur-evolutie, spectra en late-time bolometrische dalingen zijn noodzakelijk om onderscheid te maken tussen een lage ejectamassa, hoge menging en motor-energie.

Het werk benadrukt dat hoewel snmix een vereenvoudigd semi-analytisch instrument is, het effectief de impact van de radiale structuur isoleert, wat een noodzakelijke correctie vormt op de standaard one-zone interpretatie van supernova lichtcurves.