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

Résumé technique : SIMPLE – Interprétation stellaire des données météoritiques et traçage pour tous

Énoncé du problème
Les météorites globales et leurs inclusions présentent une variabilité isotopique pour de nombreux éléments chimiques, une signature d'événements de nucléosynthèse dans les étoiles et les supernovae survenus avant la formation du Système solaire. Bien que l'identification de ces variations fournisse des informations sur l'environnement de naissance du Soleil, les origines astrophysiques exactes restent débattues. Un défi majeur pour résoudre cela est la comparaison entre les prédictions des modèles de nucléosynthèse stellaire et les données météoritiques. Plus précisément, les supernovae à effondrement de cœur (CCSN) sont d'un intérêt particulier en raison de la courte durée de vie de leurs progéniteurs et de leur rôle dans la production de poussières présolaires. Cependant, comparer les prédictions des modèles aux données météoritiques est compliqué par les effets non linéaires introduits lors de la normalisation interne, la complexité de la réduction des données, et l'absence de profils d'abondance détaillés et accessibles (en fonction de la masse éjectée) pour de nombreux modèles modernes de CCSN. Historiquement, seuls quelques modèles plus anciens (par exemple, Meyer et al. 1995 ; Rauscher et al. 2002) ont été utilisés pour de telles analyses, ne parvenant pas à capturer deux décennies de progrès dans la compréhension des CCSN. De plus, distinguer les effets de nucléosynthèse de la fractionnement chimique dépendant de la masse et de la désintégration radioactive nécessite un traitement soigneux des fractions molaires par rapport aux fractions de masse et un traitement spécifique des isotopes instables.

Méthodologie
Pour répondre à ces défis, les auteurs ont développé SIMPLE (Stellar Interpretation of Meteoritic Data and Plotting for Everyone), un outil Python open-source conçu pour comparer les compositions isotopiques des éjectas de CCSN avec les données météoritiques. La méthodologie comprend :

1. Intégration des données : L'outil ingère actuellement 18 modèles de CCSN provenant de six ensembles publiés distincts (Ra02, Pi16, Si18, Ri18, LC18, La22), couvrant des masses initiales de 15, 20 et 25 $M_\odot$ avec une métallicité solaire. Ces modèles utilisent divers codes de progéniteurs (Kepler, MESA, GENEC, FRANEC) et divers mécanismes d'explosion (piston, onde de choc de type Sedov, bombe cinétique, piloté par les neutrinos).
2. Traitement des données : Le code lit les abondances isotopiques tabulées (fractions de masse) en fonction de la masse éjectée. Il les convertit en fractions molaires ($Y_i$) pour assurer la cohérence avec les abondances solaires/terrestres standards, évitant ainsi des anomalies artificielles causées par les facteurs de nombre de masse ($A_i$) dans les calculs de rapports.
3. Gestion de la désintégration radioactive : Les utilisateurs peuvent choisir d'omettre les noyaux radioactifs ou d'ajouter leurs abondances complètes aux noyaux filles stables. Le code tient compte des échelles de temps courtes des éjectas de CCSN (heures après l'explosion) par rapport à la formation de la poussière (mois/années), permettant l'inclusion de parents à courte durée de vie (ex: $^{60}$Co) tout en offrant la flexibilité pour des parents à plus longue durée de vie comme le $^{60}$Fe.
4. Normalisation et dilution : Pour comparer avec les données météoritiques, le code calcule des lignes de mélange. Il dilue les abondances stellaires avec les abondances solaires et applique une normalisation interne en utilisant un rapport d'isotope de référence. Le code implémente la loi exponentielle standard pour la normalisation (calcul des valeurs $\epsilon$) et inclut une méthode d'approximation linéaire (Dauphas et al. 2004) comme alternative. Il détermine automatiquement un facteur de dilution pour mettre à l'échelle le plus grand écart à $\pm 1$ unité $\epsilon$, évitant ainsi des résultats spécieux dans les régions présentant une nucléosynthèse minimale.
5. Visualisation : L'outil génère des graphiques de rapports isotopiques en fonction de la masse éjectée, divisant les éjectas en zones basées sur les produits de combustion dominants (ex: Si, O/Si, O/C) suivant la classification de Meyer et al. (1995). Il trace également des corrélations entre rapports isotopiques (ex: $\epsilon^{64}\text{Ni}$ vs $\epsilon^{62}\text{Ni}$) pour identifier les lignes de mélange.

Contributions clés

• Publication du logiciel : La publication de SIMPLE en tant que package Python collaboratif et open-source (disponible sur GitHub et PyPI) qui standardise le traitement des sorties de divers modèles de CCSN.
• Jeu de données complet : L'intégration de 18 modèles provenant de six groupes de recherche indépendants, couvrant des modèles de 2002 à 2022, offrant une vue plus large des incertitudes systématiques et des caractéristiques robustes de la nucléosynthèse des CCSN.
• Rigueur méthodologique : Mise en œuvre d'une conversion cohérente entre fractions de masse et fractions molaires et gestion flexible des chaînes de désintégration radioactive pour prévenir les artefacts dans les calculs de rapports isotopiques.
• Application par étude de cas : Une analyse détaillée des isotopes du Nickel (Ni) à travers les 18 modèles, démontrant la capacité de l'outil à reproduire et étendre les découvertes précédentes.

Résultats
En utilisant les isotopes du Ni comme étude de cas, les auteurs ont analysé les modèles de 20 $M_\odot$ de l'ensemble des six séries :

• Signatures de zone : L'analyse révèle des signatures isotopiques distinctes dans différentes zones d'éjectas. La zone riche en Si produit systématiquement une pente de $\epsilon^{64}\text{Ni}$ vs $\epsilon^{62}\text{Ni}$ (normalisée en interne à $^{58}\text{Ni}/^{61}\text{Ni}$) d'environ 3,0. Ce résultat est robuste à travers différentes masses de progéniteurs et séries de modèles, à condition qu'une zone Si soit éjectée.
• Comparaison avec les météorites : La pente observée dans les données météoritiques globales ($3,003 \pm 0,166$) est compatible avec un matériau provenant de la région la plus interne, riche en Si, des éjectas de CCSN.
• Variabilité des modèles : Alors que la zone Si produit des résultats cohérents, d'autres zones (ex: O/Si, O/C, He/C) montrent une variabilité significative des pentes selon le mécanisme d'explosion spécifique, l'occurrence de fusions de couches C-O, et l'activation des processus de capture neutronique (ex: la réaction $^{22}\text{Ne}(\alpha,n)^{25}\text{Mg}$).
• Contributions radioactives : L'outil quantifie avec succès la contribution des parents radioactifs (spécifiquement le $^{60}$Fe via le $^{60}$Co) à l'abondance du $^{60}\text{Ni}$ stable, illustrant comment différentes hypothèses sur le timing de désintégration affectent la composition isotopique finale.

Signification
L'article affirme que SIMPLE facilite l'interprétation des données météoritiques en fournissant une plateforme standardisée et flexible pour comparer un large éventail de modèles stellaires. En confirmant que les anomalies de Ni dans les météorites globales sont compatibles avec le matériel provenant de la région riche en Si des éjectas de CCSN, ce travail renforce l'hypothèse selon laquelle ces supernovae ont contribué à l'hétérogénéité du système solaire primitif. Les auteurs soulignent que l'outil est conçu comme une plateforme collaborative pour améliorer le développement de la nucléosynthèse stellaire et des études météoritiques. Ils notent que bien que la version actuelle se concentre sur les supernovae à cœur de fer, le cadre est construit pour incorporer de futurs modèles, y compris ceux issus de simulations 3D, des supernovae à capture électronique, et d'autres sources stellaires (ex: étoiles AGB, novae), ainsi que des fonctionnalités pour mélanger différentes couches d'éjectas.