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

Resumen Técnico: SIMPLE – Interpretación Estelar de Datos Meteoríticos y Representación Gráfica para Todos

Planteamiento del Problema
Los meteoritos masivos y sus inclusiones exhiben variabilidad isotópica en muchos elementos químicos, una firma de eventos nucleosintéticos en estrellas y supernovas previos a la formación del Sistema Solar. Si bien la identificación de estas variaciones proporciona información sobre el entorno de nacimiento del Sol, los orígenes astrofísicos exactos siguen siendo objeto de debate. Un desafío principal para resolver esto es la comparación entre los datos de isótopos meteoríticos y las predicciones de los modelos de nucleosíntesis estelar. Específicamente, las supernovas de colapso de núcleo (CCSN, por sus siglas en inglés) son de particular interés debido a sus cortos tiempos de vida de los progenitores y su papel en la producción de polvo presolar. Sin embargo, comparar las predicciones de los modelos con los datos meteoríticos se ve complicado por los efectos no lineales introducidos durante la normalización interna, la complejidad de la reducción de datos y la falta de perfiles de abundancia detallados y accesibles (en función de la masa eyectada) para muchos modelos modernos de CCSN. Históricamente, solo se han utilizado para tales análisis unos pocos modelos antiguos (por ejemplo, Meyer et al. 1995; Rauscher et al. 2002), lo que no logra capturar dos décadas de progreso en la comprensión de las CCSN. Además, distinguir los efectos nucleosintéticos de la fraccionación química dependiente de la masa y la desintegración radiactiva requiere un manejo cuidadoso de las fracciones molares frente a las de masa y un tratamiento específico de los isótopos inestables.

Metodología
Para abordar estos desafíos, los autores desarrollaron SIMPLE (Stellar Interpretation of Meteoritic Data and Plotting for Everyone), un kit de herramientas de Python de código abierto diseñado para comparar las composiciones isotópicas de los eyectos de CCSN con los datos meteoríticos. La metodología consiste en:

1. Integración de Datos: El kit de herramientas actualmente ingiere 18 modelos de CCSN de seis conjuntos publicados distintos (Ra02, Pi16, Si18, Ri18, LC18, La22), que cubren masas iniciales de 15, 20 y 25 $M_\odot$ con metalicidad solar. Estos modelos utilizan varios códigos de progenitores (Kepler, MESA, GENEC, FRANEC) y mecanismos de explosión (pistón, onda de choque de Sedov, bomba cinética, impulsada por neutrinos).
2. Procesamiento de Datos: El código lee abundancias isotópicas tabuladas (fracciones de masa) en función de la masa eyectada. Convierte estas a fracciones molares ($Y_i$) para asegurar la consistencia con las abundancias solares/terrestres estándar, evitando anomalías artificiales causadas por factores de número de masa ($A_i$) en los cálculos de proporciones.
3. Manejo de la Desintegración Radiactiva: Los usuarios pueden elegir omitir los núcleos radiactivos o añadir sus abundancias completas a los núcleos hijas estables. El código tiene en cuenta las escalas de tiempo cortas de los eyectos de CCSN (horas post-explosión) frente a la formación de polvo (meses/años), permitiendo la inclusión de padres de vida corta (por ejemplo, $^{60}$Co) mientras ofrece flexibilidad para padres de vida más larga como el $^{60}$Fe.
4. Normalización y Dilución: Para comparar con los datos meteoríticos, el código calcula líneas de mezcla. Diluye las abundancias estelares con abundancias solares y aplica una normalización interna utilizando una relación de isótopos de referencia. El código implementa la ley exponencial estándar para la normalización (calculando valores $\epsilon$) e incluye un método de aproximación lineal (Dauphas et al. 2004) como alternativa. Determina automáticamente un factor de dilución para escalar el mayor desfase a $\pm 1$ unidad $\epsilon$, evitando resultados espurios en regiones con mínima nucleosíntesis.
5. Visualización: La herramienta genera gráficos de relaciones isotópicas contra la masa eyectada, dividiendo el eyecto en zonas basadas en los productos de combustión dominantes (por ejemplo, Si, O/Si, O/C) siguiendo la clasificación de Meyer et al. (1995). También grafica correlaciones entre relaciones isotópicas (por ejemplo, $\epsilon^{64}\text{Ni}$ vs. $\epsilon^{62}\text{Ni}$) para identificar líneas de mezcla.

Contribuciones Clave

• Lanzamiento de Software: El lanzamiento de SIMPLE como un paquete de Python colaborativo y de código abierto (disponible en GitHub y PyPI) que estandariza el procesamiento de diversos resultados de modelos de CCSN.
• Conjunto de Datos Exhaustivo: La integración de 18 modelos de seis grupos de investigación independientes, que abarcan modelos desde 2002 hasta 2022, proporcionando una visión más amplia de las incertidumbres sistemáticas y las características robustas en la nucleosíntesis de CCSN.
• Rigor Metodológico: Implementación de una conversión consistente entre fracciones de masa y molares y un manejo flexible de las cadenas de desintegración radiactiva para prevenir artefactos en los cálculos de relaciones isotópicas.
• Aplicación de Caso de Estudio: Un análisis detallado de los isótopos de Níquel (Ni) en todos los 18 modelos, demostrando la capacidad de la herramienta para reproducir y extender hallazgos previos.

Resultos
Utilizando los isótopos de Ni como caso de estudio, los autores analizaron los modelos de 20 $M_\odot$ de los seis conjuntos:

• Firmas de Zona: El análisis revela firmas isotópicas distintas en diferentes zonas de eyecta. La zona rica en Si produce consistentemente una pendiente de $\epsilon^{64}\text{Ni}$ vs. $\epsilon^{62}\text{Ni}$ (internamente normalizada a $^{58}\text{Ni}/^{61}\text{Ni}$) de aproximadamente 3.0. Este resultado es robusto a través de diferentes masas de progenitores y conjuntos de modelos, siempre que se eyecte una zona de Si.
• Comparación con Meteoritos: La pendiente observada en los datos meteoríticos masivos ($3.003 \pm 0.166$) es compatible con material originado en la región más interna rica en Si de los eyectos de CCSN.
• Variabilidad de los Modelos: Mientras que la zona de Si arroja resultados consistentes, otras zonas (por ejemplo, O/Si, O/C, He/C) muestran una variabilidad significativa en las pendientes dependiendo del mecanismo de explosión específico, la ocurrencia de fusiones de capas C-O y la activación de procesos de captura de neutrones (por ejemplo, la reacción $^{22}\text{Ne}(\alpha,n)^{25}\text{Mg}$).
• Contribuciones Radiactivas: La herramienta cuantifica con éxito la contribución de los padres radiactivos (específicamente $^{60}\text{Fe}$ vía $^{60}\text{Co}$) a la abundancia del $^{60}\text{Ni}$ estable, ilustrando cómo diferentes supuestos sobre el tiempo de desintegración afectan la composición isotópica final.

Significancia
El artículo afirma que SIMPLE facilita la interpretación de los datos meteoríticos al proporcionar una plataforma estandarizada y flexible para comparar una amplia gama de modelos estelares. Al confirmar que las anomalías de Ni en meteoritos masivos son compatibles con el material de la región rica en Si de los eyectos de CCSN, el trabajo refuerza la hipótesis de que tales supernovas contribuyeron a la heterogeneidad del sistema solar temprano. Los autores enfatizan que la herramienta está diseñada como una plataforma colaborativa para mejorar el desarrollo de los estudios de nucleosíntesis estelar y meteorítica. Señalan que, si bien la versión actual se centra en las supernovas de colapso de núcleo de Fe, el marco está construido para incorporar futuros modelos, incluyendo aquellos de simulaciones 3D, supernovas de captura de electrones y otras fuentes estelares (por ejemplo, estrellas AGB, novae), así como funciones para mezclar diferentes capas de eyecta.