Stellar Interpretation of Meteoritic Data and PLotting for Everyone (SIMPLE): Isotope Mixing Lines for Six Sets of Core-Collapse Supernova Models

Technische Samenvatting: SIMPLE – Stellaire Interpretatie van Meteoritische Data en Plotten voor Iedereen

Probleemstelling
Bulkmeteorieten en hun inclusies vertonen isotopische variabiliteit over vele chemische elementen, een signatuur van nucleosynthetische gebeurtenissen in sterren en supernova's voorafgaand aan de vorming van het Zonnestelsel. Hoewel het identificeren van deze variaties inzichten biedt in de geboorteomgeving van de Zon, blijven de exacte astrofysische oorsprong onderwerp van debat. Een primaire uitdaging bij het oplossen hiervan is de vergelijking tussen meteoritische isotoopdata en voorspellingen uit stellaire nucleosynthesemodellen. Specifiek zijn kerninstortende supernova's (CCSNe) van bijzonder belang vanwege hun korte progenitor-levensduur en hun rol bij de productie van presolaire stofdeeltjes. Het vergelijken van modelvoorspellingen met meteoritische data wordt echter bemoeilijkt door niet-lineaire effecten die optreden tijdens interne normalisatie, de complexiteit van datareductie, en het gebrek aan toegankelijke, gedetailleerde abundantieprofielen (als functie van geëjecteerde massa) voor veel moderne CCSN-modellen. Historisch gezien zijn slechts enkele oudere modellen (bijv. Meyer et al. 1995; Rauscher et al. 2002) gebruikt voor dergelijke analyses, wat er toe leidt dat twee decennia aan vooruitgang in het begrip van CCSN niet worden gevangen. Bov�elijk vereist het onderscheiden van nucleosynthetische effecten van chemische massafhankelijkem fractionering en radioactief verval een zorgvuldige behandeling van molaire versus massafracties en een specifieke behandeling van instabiele isotopen.

Methodologie
Om deze uitdagingen aan te pakken, hebben de auteurs SIMPLE (Stellar Interpretation of Meteoritic Data and Plotting for Everyone) ontwikkeld, een open-source Python-toolkit die is ontworpen om de isotopische samenstellingen van CCSN-ejecta te vergelijken met meteoritische data. De methodologie omvat:

1. Data-integratie: De toolkit verwerkt momenteel 18 CCSN-modellen uit zes verschillende gepubliceerde sets (Ra02, Pi16, Si18, Ri18, LC18, La22), die initiële massa's van 15, 20 en 25 $M_\odot$ met een zonne metalliciteit beslaan. Deze modellen maken gebruik van diverse progenitor-codes (Kepler, MESA, GENEC, FRANEC) en explosiemechanismen (piston-gedreven, Sedov blast-wave, kinetische bom, neutrino-gedreven).
2. Dataverwerking: De code leest getabelleerde isotopische abundanties (massafracties) als functie van de geëjecteerde massa. Deze worden omgezet naar molaire fracties ($Y_i$) om consistentie met standaard zonne/aardse abundanties te waarborgen, waardoor kunstmatige anomalieën veroorzaakt door massagetallen ($A_i$) in ratioberekeningen worden voorkomen.
3. Afhandeling van Radioactief Verval: Gebruikers kunnen ervoor kiezen om radioactieve kernen weg te laten of hun volledige abundanties toe te voegen aan stabiele dochterkernen. De code houdt rekening met de korte tijdschalen van CCSN-ejecta (uren na de explosie) versus stofvorming (maanden/jaren), waardoor de inclusie van kortlevende ouders (bijv. $^{60}$Co) mogelijk is, terwijl flexibiliteit wordt geboden voor langere levende ouders zoals $^{60}$Fe.
4. Normalisatie en Verdunning: Om de CCSN-ejecta te vergelijken met meteoritische data, berekent de code menglijnen. Het verdunt stellaire abundanties met zonne abundanties en past interne normalisatie toe met behulp van een referentie-isotoopratio. De code implementeert de standaard exponentiële wet voor normalisatie (het berekenen van $\epsilon$-waarden) en bevat een lineaire benaderingsmethode (Dauphas et al. 2004) als alternatief. Het bepaalt automatisch een verdunningsfactor om de grootste afwijking te schalen naar $\pm 1$ $\epsilon$-eenheid, om zo spookresultaten in regio's met minimale nucleosynthese te vermijden.
5. Visualisatie: De tool genereert plots van isotoopratio's tegen de geëjecteerde massa, waarbij het ejecta wordt onderverdeeld in zones op basis van dominante verbrandingsproducten (bijv. Si, O/Si, O/C), volgens de classificatie van Meyer et al. (1995). Het plot ook correlaties tussen isotoopratio's (bijv. $\epsilon^{64}\text{Ni}$ vs. $\epsilon^{62}\text{Ni}$) om menglijnen te identificeren.

Belangrijkste Bijdragen

• Software Release: De release van SIMPLE als een collaboratieve, open-source Python-package (beschikbaar op GitHub en PyPI) die de verwerking van diverse CCSN-modeloutputs standaardiseert.
• Omvattende Dataset: De integratie van 18 modellen uit zes onafhankelijke onderzoeksgroepen, reikend van modellen uit 2002 tot 2022, wat een breder beeld geeft van systematische onzekerheden en robuuste kenmerken in CCSN-nucleosynthese.
• Methodologische Rigor: Implementatie van consistente conversie tussen massafracties en molaire fracties en flexibele afhandeling van radioactieve vervalketens om artefacten in isotoopratioberekeningen te voorkomen.
• Casusstudie Toepassing: Een gedetailleerde analyse van Nikkel (Ni) isotopen over alle 18 modellen, die het vermogen van de tool demonstreert om eerdere bevindingen te reproduceren en uit te breiden.

Resultaten
Met Ni-isotopen als casusstudie analyseerden de auteurs de 20 $M_\odot$ modellen uit alle zes de sets:

• Zone-signaturen: De analyse onthult duidelijke isotopische signatures in verschillende ejecta-zones. De Si-rijke zone produceert consistent een helling van $\epsilon^{64}\text{Ni}$ vs. $\epsilon^{62}\text{Ni}$ (intern genormaliseerd naar $^{58}\text{Ni}/^{61}\text{Ni}$) van ongeveer 3,0. Dit resultaat is robuust over verschillende progenitor-massa's en modelsets heen, mits een Si-zone wordt geëjecteerd.
• Vergelijking met Meteorieten: De geobserveerde helling in bulkmeteoritische data ($3,003 \pm 0,166$) is compatibel met materiaal dat afkomstig is uit de binnenste Si-rijke regio van CCSN-ejecta.
• Modelvariabiliteit: Terwijl de Si-zone consistente resultaten oplevert, vertonen andere zones (bijv. O/Si, O/C, He/C) aanzienlijke variabiliteit in hellingen, afhankelijk van het specifieke explosiemechanisme, het voorkomen van C-O schil-fusies en de activatie van neutronenvangstprocessen (bijv. de $^{22}\text{Ne}(\alpha,n)^{25}\text{Mg}$ reactie).
• Radioactieve Bijdragen: De tool kwantificeert succesvol de bijdrage van radioactieve ouders (specifiek $^{60}$Fe via $^{60}$Co) aan de abundantie van stabiel $^{60}\text{Ni}$, en illustreert hoe verschillende aannames over de timing van verval de uiteindelijke isotopische samenstelling beïnvloeden.

Betekenis
Het artikel stelt dat SIMPLE de interpretatie van meteoritische data vergemakkelijkt door een gestandaardiseerd, flexibel platform te bieden om een breed scala aan stellaire modellen te vergelijken. Door te bevestigen dat bulkmeteoriet Ni-anomalieën compatibel zijn met materiaal uit de Si-rijke regio van CCSN-ejecta, versterkt dit werk de hypothese dat dergelijke supernova's hebben bijgedragen aan de heterogeniteit van het vroege Zonnestelsel. De auteurs benadrukken dat de tool is ontworpen als een collaboratief platform om de ontwikkeling van stellaire nucleosynthese en meteoritische studies te bevorderen. Zij merken op dat hoewel de huidige versie zich richt op Fe-kerninstortende supernova's, het framework is gebouwd om toekomstige modellen te incorporeren, inclusief die van 3D-simulaties, elektron-capture supernova's en andere stellaire bronnen (bijv. AGB-sterren, nova's), evenals functies voor het mengen van verschillende ejecta-lagen.