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

Resumen Técnico: Firmas de la Mezcla de $^{56}\text{Ni}$ y Eyecta Rica en Neutrones en Supernovas

Planteamiento del Problema
Las curvas de luz de las supernovas (SN) se interpretan rutinariamente utilizando modelos de decaimiento radiactivo de una sola zona (por ejemplo, Arnett 1980, 1982) para inferir propiedades de la progenitora, tales como la masa de la eyecta ($M_{\text{ej}}$), la masa de níquel ($M_{\text{Ni}}$) y la energía de la explosión ($E_k$). Estos modelos asumen un único tiempo de difusión efectivo e ignoran la distribución espacial del material radiactivo. Sin embargo, las simulaciones hidrodinámicas indican que el $^{56}\text{Ni}$ se sintetiza en lo profundo de la eyecta y posteriormente se mezcla hacia el exterior mediante inestabilidades de Rayleigh–Taylor, convección y turbulencia. Esta distribución espacial es una variable física crítica que los modelos de una sola zona omiten. En consecuencia, las curvas de luz de ejectas altamente mezcladas pueden malinterpretarse como poseedoras de una masa de eyecta baja o como requiriendo potencia adicional de un motor central (por ejemplo, magnetares o acreción de fallback) para explicar ascensos rápidos y luminosidades elevadas. Además, la observabilidad del enriquecimiento por proceso r en colapsares depende sensiblemente de la ubicación y la geometría de la eyecta rica en neutrones respecto al fondo impulsado por el níquel, un factor que suele simplificarse excesivamente en modelos esféricos.

Metodología
Los autores introducen snmix, un modelo semi-analítico de múltiples capas implementado como un plugin basado en JAX para el marco de trabajo Redback. A diferencia de los modelos de una sola zona, snmix discretiza la eyecta en expansión homóloga en $N$ capas con velocidades espaciadas geométricamente.

• Estructura Radiactiva: El modelo permite distribuciones radiales arbitrarias de $^{56}\text{Ni}$ (parametrizadas por una fracción de mezcla $f_{\text{mix}}$ en un perfil de tapa plana o de núcleo más cola) y de material de proceso r.
• Física de Transporte: Cada capa evoluciona de forma independiente, rastreando la energía térmica, las pérdidas adiabáticas y la difusión radiativa. Crucialmente, el modelo calcula la difusión y la deposición de rayos gamma basándose en la densidad de columna acumulada exterior a cada capa específica, en lugar de un promedio global.
• Opacidad y SED: El modelo incorpora opacidades grises y una ley de opacidad dependiente de la temperatura. Para los cálculos de proceso r, incluye opacidades dependientes de lantánidos y tasas de calentamiento analíticas calibradas mediante cálculos de redes nucleares.
• Geometría: Aunque es principalmente esférico, el modelo se extiende a una geometría de "viento de disco" de dos componentes para colapsares, combinando un componente de SN polar con un viento ecuatorial rico en neutrones, ponderado por el ángulo de visión.
• Inferencia: Los autores utilizan este marco para generar curvas de luz bolométricas sintéticas y realizar inferencia bayesiana (usando muestreo NUTS) para comparar los ajustes de una sola zona contra la realidad de múltiples capas, analizando los sesgos de los parámetros.

Contribuciones Clave y Resultados

1. Impacto de la Mezcla de $^{56}\text{Ni}$ en las Curvas de Luz:

   • La mezcla hacia el exterior del $^{56}\text{Ni}$ reduce la columna de difusión suprayacente y la columna de atrapamiento de rayos gamma. Esto resulta en ascensos más rápidos y brillantes, y colas radiactivas más pronunciadas para una $M_{\text{Ni}}$, $M_{\text{ej}}$ y $E_k$ fijas.
   • Un ascenso rápido y brillante no es evidencia exclusiva de una masa de eyecta baja o de un motor; los modelos radiactivos altamente mezclados pueden reproducir estas características sin fuentes de potencia adicionales.

2. Sesgo en la Inferencia de Una Sola Zona:

   • Ajustar curvas de luz mezcladas con modelos de una sola zona produce ajustes visualmente buenos pero con parámetros físicos sesgados.
   • A una opacidad fija, el modelo absorbe la estructura radial faltante infiriendo una $M_{\text{ej}}$ significativamente menor y una fracción de níquel ($f_{\text{Ni}}$) más alta, mientras que la $M_{\text{Ni}}$ permanece relativamente estable (recuperada dentro de un $\sim3\%$ con opacidad fija).
   • Si se permite que la opacidad varíe, el sesgo se desplaza hacia la inferencia de una opacidad efectiva mucho menor ($\kappa_{\text{fit}}/\kappa_{\text{input}} \simeq 0.24$ para casos totalmente mezclados).
   • Estos sesgos se propagan a inferencias incorrectas sobre las energías de explosión y los mapeos de progenitores.

3. Caso de Estudio: SN 1998bw:

   • La aplicación de snmix a la SN 1998bw arroja una fracción de mezcla moderada ($f_{\text{mix}} \approx 0.59$). Aunque la estimación de máxima verosimilitud de una sola zona se solapa con la posterior, difiere en un $\sim13\%$ en $M_{\text{ej}}$ y un $\sim5\%$ en $E_k$, ilustrando cómo los datos bolométricos por sí solos no pueden descomponer de forma única la masa, la energía, la opacidad y la mezcla.

4. Firmas de Proceso r en Colapsares:

   • La firma del enriquecimiento por proceso r no es universalmente un exceso en el infrarrojo cercano (NIR) en tiempos tardíos.
   • Ubicación: El material de proceso r enterrado suprime inicialmente la luz óptica debido a la alta opacidad, y el exceso en el NIR aparece solo después de que las capas enriquecidas se vuelven transparentes.
   • Geometría: En los modelos de viento de disco, si el material rico en neutrones está confinado al plano ecuatorial, las SN-GRB en el eje (vistas polares) muestran una firma de proceso r insignificante porque el componente de la SN polar domina. La señal es más fuerte para observadores fuera del eje o si el material se mezcla hacia latitudes más altas ($\gtrsim 30^\circ$).
   • Diagnóstico: La evolución del color ($\Delta(r-H)$) es un diagnóstico valioso, pero requiere desentrañar la supresión óptica del aumento de brillo en la banda H. El monitoreo emparejado de óptico y NIR (10–60 días) es más informativo que los puntos de banda H en tiempos tardíos por sí solos.

Significancia
El artículo argumenta que la "estructura radiactiva" (distribución radial del calentamiento y la opacidad) es una variable crítica y a menudo descuidada en la interpretación de las supernovas.

• Para las SN impulsadas por el Níquel: Desafía la invocación rutinaria de motores centrales para explicar las SN Ic-BL y las SN-GRB de ascenso rápido, sugiriendo que la mezcla hacia el exterior es una alternativa viable. Advierte que los ajustes de una sola zona pueden proporcionar ajustes estadísticos "buenos" mientras ofrecen parámetros físicamente incorrectos.
• Para los Colapsares: Refina la estrategia de búsqueda de enriquecimiento por proceso r, yendo más allá de la heurística simple de "exceso en el NIR". Destaca que las no detecciones en SN-GRB en el eje no descartan la producción de proceso r si la eyecta está confinada ecuatorialmente.
• Estrategia Observacional: Los autores concluyen que romper estas degeneraciones requiere combinar observables: las velocidades fotosféricas, la evolución del color, los espectros y los descensos bolométricos en tiempos tardíos son necesarios para distinguir entre una masa de eyecta baja, una mezcla alta, o la potencia de un motor.

El trabajo enfatiza que, aunque snmix es una herramienta semi-analítica simplificada, efectivamente aísla el impacto de la estructura radial, proporcionando una corrección necesaria a la interpretación estándar de una sola zona de las curvas de luz de las supernovas.