Lightcurve Modelling of 2,205 ZTF DR2 Type~Ia Supernovae: Implications for SN Ia Physics and Cosmology

Resumen Técnico: Modelado de Curvas de Luz de 2,205 Supernovas Tipo Ia de ZTF DR2

Problema y Motivación
Las supernovas de tipo Ia (SNe Ia) sirven como candelas estandarizables para la cosmología, basándose en relaciones empíricas como la relación "más brillante-más lenta" (Phillips) y la relación "más brillante-más azul". Sin embargo, los orígenes físicos de estas relaciones y la diversidad de los mecanismos de explosión (por ejemplo, progenitores de masa de Chandrasekhar frente a sub-Chandrasekhar) siguen sin comprenderse completamente. Estudios sistemáticos previos de parámetros físicos (masa de los eyectos $M_{\rm ej}$, masa de níquel sintetizado $M_{\rm Ni}$) se vieron limitados por tamaños de muestra pequeños (por ejemplo, Scalzo et al. 2014 analizaron 337 SNe) o dependieron de relaciones empíricas en lugar del ajuste directo de las curvas de luz. Además, los modelos estándar de decaimiento radiactivo de una sola zona suelen fallar al reproducir los máximos secundarios observados en las bandas rojas/infrarrojo cercano, los cuales se atribuyen a la física de recombinación Fe III $\to$ Fe II. Este artículo tiene como objetivo realizar la extracción más grande y sistemática de parámetros físicos a partir de curvas de luz de SNe Ia hasta la fecha, aprovechando el Zwicky Transient Facility (ZTF) Data Release 2 (DR2) para restringir la física de la explosión, los canales de los progenitores y la base física de la estandarización cosmológica.

Metodología
Los autores analizaron 2,205 SNe Ia confirmadas espectroscópicamente y de calidad cosmológica de la muestra ZTF DR2. El análisis empleó un marco de inferencia bayesiana utilizando un modelo de decaimiento radiactivo de una sola zona modificado. Las características metodológicas clave incluyen:

• Opacidad Dependiente de la Longitud de Onda: El modelo estándar de opacidad gris se extendió fenomenológicamente para incluir la física de recombinación Fe III $\to$ Fe II. La opacidad óptica $\kappa$ se hizo función de la longitud de onda ($\lambda$) y el tiempo ($t$), permitiendo que el modelo reproduzca los máximos secundarios en las bandas rojas ($i, z$) al reducir la opacidad en longitudes de onda más largas, aumentando así la temperatura y el flujo fotosférico.
• Parámetros Físicos: El modelo constreñyó simultáneamente 10 parámetros libres por evento, incluyendo la masa de los eyectos ($M_{\rm ej}$), la fracción de níquel ($f_{\rm Ni} = M_{\rm Ni}/M_{\rm ej}$), la velocidad de los eyectos ($v_{\rm ej}$), el tiempo de explosión ($t_0$) y los parámetros de opacidad (opacidad base, época de recombinación, anchura y fuerza). La extinción del anfitrión ($A_V$) también fue modelada.
• Marco de Inferencia: Se utilizó una verosimilitud de mezcla gaussiana para manejar valores atípicos y datos dispersos. Los autores utilizaron el muestreador anidado PyMultiNest a través de la interfaz Bilby para muestrear las distribuciones posteriores de todos los 2,205 eventos.
• Modelado Jerárquico: Para abordar los sesgos de selección inherentes a la muestra de calidad cosmológica (que excluye eventos subluminosos como los de tipo 91bg), se realizó un análisis bayesiano jerárquico sobre una submuestra limitada por desplazamiento al rojo ($z \le 0.06$, $N=902$). Este enfoque modeló las distribuciones a nivel de población de $M_{\rm ej}$ y $M_{\rm Ni}$ como gaussianas truncadas mientras se marginalizaba sobre las incertidumbres individuales de los eventos y una función de selección dependiente del brillo.

Resultados Clave

1. Distribuciones de Parámetros Físicos:

   • Para la muestra completa, la masa mediana de los eyectos es $M_{\rm ej} \approx 1.40 M_\odot$ y la masa de níquel mediana es $M_{\rm Ni} \approx 1.19 M_\odot$. La alta $M_{\rm Ni}$ se atribuye al sesgo de selección que favorece a los eventos luminosos.
   • La distribución de la masa de los eyectos es ampliamente unimodal, con un pico cerca de la masa de Chandrasekhar, sin una bimodalidad clara que sugiera distintos canales de progenitores.
   • Subtipos: Los eventos tipo SN 1991T son físicamente distintos, exhibiendo masas de eyectos sistemáticamente más altas (mediana de $1.64 M_\odot$ frente a $1.38 M_\odot$ para las normales) y produciendo un $\sim30\%$ más de níquel. El 84% de los eventos tipo 1991T tienen $M_{\rm ej} > 1.4 M_\odot$.

2. Correlaciones con Parámetros Empíricos:

   • Se encontró una fuerte correlación positiva ($\rho = 0.57$) entre la masa de níquel sintetizado ($M_{\rm Ni}$) y el parámetro de estiramiento (stretch) de SALT2 ($x_1$). Esto proporciona una base física para la relación "más brillante-más lenta": las SNe más brillantes sintetizan más níquel y, dentro de este modelo, poseen masas de eyectos mayores que conducen a tiempos de difusión de fotones más largos.
   • El estiramiento de SALT2 ($x_1$) se correlaciona con la masa de los eyectos ($M_{\rm ej}$), mientras que el color de SALT2 ($c$) no muestra correlación con $M_{\rm ej}$, lo que sugiere que el color está dominado por la extinción de polvo y los efectos de recombinación secundaria, en lugar de propiedades de la explosión a gran escala.

3. Dependencias de la Galaxia Anfitriona:

   • Las SNe Ia en anfitriones de baja masa ($\log_{10}(M_*/M_\odot) < 10$) sintetizan un $\sim12\%$ más de níquel ($\Delta M_{\rm Ni} = 0.13 M_\odot$) que aquellas en anfitriones de alta masa.
   • Esta diferencia intrínseca en la producción de níquel ofrece una explicación física para el "escalón de masa" observado en los residuos de Hubble. Los autores sugieren que la estandarización basada en parámetros físicos (como $M_{\rm Ni}$) podría potencialmente eliminar este escalón, mientras que la estandarización empírica (SALT2) deja un desfase residual.

4. Inferencia de Población Jerárquica:

   • Corrigiendo los sesgos de selección en la submuestra de bajo desplazamiento al rojo ($z \le 0.06$), los autores infirieron las verdaderas distribuciones de la población.
   • La población está bien descrita por distribuciones gaussianas con $\mu_{\rm ej} = 1.26 \pm 0.01 M_\odot$ y $\mu_{\rm Ni} = 0.64 \pm 0.06 M_\odot$.
   • Las fracciones inferidas son: $43 \pm 2\%$ sub-Chandrasekhar ($M_{\rm ej} < 1.2 M_\odot$), $34 \pm 1\%$ cerca de Chandrasekhar ($1.2 \le M_{\rm ej} \le 1.5 M_\odot$) y $24 \pm 2\%$ super-Chandrasekhar ($M_{\rm ej} > 1.5 M_\odot$).
   • La distribución suave y continua y la división similar entre eventos sub- y cerca de Chandrasekhar sugieren un único mecanismo de explosión dominante (probablemente detonación retardada) que opera a través de un rango de masas de progenitores, en lugar de canales distintos y separados.

Significancia y Reivindicaciones
El artículo afirma proporcionar la caracterización física homogénea más grande de la población local de SNe Ia hasta la fecha. Sus principales contribuciones son:

• Base Física para la Estandarización: Establece un vínculo directo entre la relación empírica de Phillips y la física fundamental de la explosión (masa de níquel y masa de los eyectos), confirmando que la relación "más brillante-más lenta" es una manifestación de la física de difusión.
• Origen del Escalón de Masa: Vincula el escalón de masa de la galaxia anfitriona con diferencias intrínsecas en la eficiencia de producción de níquel, probablemente impulsadas por la edad o la metalicidad del progenitor, en lugar de puramente de sistemáticas observacionales.
• Restricciones de los Progenitores: Al analizar una muestra dos órdenes de magnitud más grande que los esfuerzos previos de modelado físico, el estudio sugiere que la población de SNe Ia útil para la cosmología está dominada por explosiones de masa cercana a la de Chandrasekhar con una distribución continua de masas de progenitores, desafiando la noción de canales de explosión distintos y separados para diferentes subtipos.
• Factibilidad Metodológica: El trabajo demuestra que los modelos físicamente motivados y dependientes de la longitud de onda pueden aplicarse a miles de curvas de luz con una eficiencia computacional comparable al ajuste empírico, allanando el camino para la estandarización física en futuros grandes estudios (por ejemplo, Rubin, Roman).

Los autores se mantienen modestos respecto a las incertidumbres sistemáticas, reconociendo que su modelo de una sola zona y los tratamientos simplificados de opacidad pueden introducir sesgos. Enfatizan que sus resultados caracterizan la población "útil para la cosmología" (sesgada hacia eventos luminosos) en lugar de toda la diversidad de los transitorios térmicos, y que el modelado jerárquico es esencial para recuperar las propiedades reales de la población.