Decoding the Early-Time Light Curves of Type Ia Supernovae. II. Population Parameters of One Thousand ZTF Supernovae

Resumen Técnico: Decodificación de las curvas de luz en tiempos tempranos de las Supernovas de Tipo Ia. II. Parámetros de la población de mil Supernovas de ZTF

Planteamiento del problema
Las supernovas de tipo Ia (SNe Ia) son indicadores de distancia cosmológica estandarizables, aunque sus mecanismos de explosión y sistemas progenitores subyacentes siguen siendo objeto de debate. Mientras que las curvas de luz en tiempos tardíos cerca del brillo máximo están bien caracterizadas, la emisión en tiempos tempranos (dentro de la primera semana) codifica información crítica sobre los eyecta externos, la distribución de $^{56}\text{Ni}$ radiactivo y las posibles interacciones con material circunestelar (CSM) o compañeras. Históricamente, los estudios de las SNe Ia en tiempos tempranos se han visto limitados por tamaños de muestra pequeños, sesgos de selección (por ejemplo, el sesgo de Malmquist) y la agregación de parámetros ajustados individualmente, lo que puede introducir correlaciones espurias y sistemáticas. Además, la relación entre la morfología del ascenso en tiempos tempranos y el parámetro de estiramiento (stretch) global de la curva de luz ($x_1$) no ha sido cuantificada rigurosamente en una población estadísticamente completa. Este artículo aborda la necesidad de una caracterización robusta, a nivel de población, de las SNe Ia en tiempos tempranos utilizando una muestra completa de volumen para romper las degeneraciones entre los canales de explosión y las propiedades de los progenitores.

Metodología
Los autores aplican un marco bayesiano jerárquico, desarrollado previamente en el "Artículo I" de esta serie, a una muestra completa de volumen de 972 SNe Ia espectroscópicamente normales de los Datos de la Zwicky Transient Facility Data Release 2 (ZTF DR2, $z < 0.06$). Esta muestra es un orden de magnitud más grande que los conjuntos de datos utilizados anteriormente para análisis similares.

Componentes metodológicos clave incluyen:

• Modelo de ascenso de ley de potencia: Las curvas de luz en tiempos tempranos se modelan mediante un ascenso de ley de potencia empírica ($F \propto (t - t_{\text{fl}})^\alpha$) hasta un umbral de truncamiento (30%, 40% o 50% del flujo de pico). El modelo incluye parámetros para el tiempo de la primera luz ($t_{\text{fl}}$), la amplitud y un índice de ascenso específico del filtro ($\alpha$).
• Inferencia jerárquica: En lugar de ajustar SNe individuales y agregar los resultados, el marco infiere simultáneamente los parámetros individuales y los hiperparámetros de la población (medias, dispersiones y covarianzas). Este enfoque pondera naturalmente de forma decreciente los eventos con muestreo escaso o ruidoso, mitigando los sesgos de selección y el sesgo de Malmquist sin requerir cortes de datos estrictos.
• Manejo de la incertidumbre: La función de verosimilitud incorpora un factor de escala ($\beta_d$) para corregir las incertididades fotométricas subestimadas en el pipeline de ZTF y propaga las incertidumbres del tiempo de pico ($t_0$) inferido a partir de los ajustes de la curva de luz completa (SALT2).
• Análisis de correlación: El estudio investiga las correlaciones entre los parámetros de tiempo temprano ($t_{\text{rise}}$, $\alpha$, evolución del color $\alpha_g - \alpha_r$) y el parámetro de estiramiento de SALT2 ($x_1$) utilizando un modelo de regresión de errores en las variables bayesiano con un ajuste lineal fragmentado para identificar rupturas estructurales.
• Comparación teórica: Las correlaciones observadas se comparan con curvas de luz sintéticas de dos marcos teóricos: modelos de masa cercana a la de Chandrasekhar (cerca-$M_{\text{Ch}}$) en 1D con perfiles de densidad y mezcla de $^{56}\text{Ni}$ variables (código TURTLS) y modelos de doble detonación de masa sub-Chandrasekhar ($M_{\text{Ch}}$) en 2D con efectos de ángulo de visión variables.

Resultos clave

1. Parámetros de la población: El estudio proporciona restricciones robustas sobre la población en tiempos tempranos:

   • Tiempo de ascenso medio ($t_{\text{rise}}$): $18.55 \pm 0.08$ días.
   • Índice de ascenso medio en la banda $r$ ($\alpha_r$): $2.10 \pm 0.04$ (más empinado que el modelo de fuego de la batalla o fireball $\alpha=2$).
   • Evolución de color media ($\alpha_g - \alpha_r$): $0.20 \pm 0.02$, lo que indica que el flujo en tiempos tempranos evoluciona más abruptamente en la banda $g$ que en la banda $r$ (evolución más azul).
   • Existe una dispersión intrínseca significativa en todos los parámetros ($\sigma_{t_{\text{rise}}} \approx 1.42$ días, $\sigma_{\alpha_r} \approx 0.48$).

2. Bifurcación con el estiramiento ($x_1$): La relación entre la morfología de tiempo temprano y el estiramiento de la curva de luz exhibe una ruptura estructural distinta cerca de $x_1 \approx 0.2\text{--}0.5$:

   • Población de alto estiramiento ($x_1 > x_{\text{brk}}$): Muestra tendencias claras donde un mayor $x_1$ se correlaciona con ascensos más suaves (menor $\alpha$), colores azules más persistentes (menor $\alpha_g - \alpha_r$) y un aplanamiento de la correlación $t_{\text{rise}}$ vs. $x_1$.
   • Población de bajo estiramiento ($x_1 < x_{\text{brk}}$): Carece de correlaciones significativas entre los parámetros de tiempo temprano y $x_1$.
   • Un "tiempo de ascenso ajustado" (medido desde un epoch de fracción de flujo en lugar de extrapolado a la primera luz) revela una correlación continua y fuerte con $x_1$ sin ruptura, lo que sugiere que la ruptura en el $t_{\text{rise}}$ estándar es impulsada por variaciones en la evolución del flujo temprano ($t_{X\%}$).

3. Anomalías y subclases peculiares:

   • Eventos tipo 03fg: Exhiben tiempos de ascenso significativamente más largos e índices de ascenso más empinados, consistentes con excesos de flujo de corta duración por interacciones de los eyecta con el CSM en fusiones violentas de dobles enanas blancas.
   • Eventos tipo 02cx y 91bg: Muestran tiempos de ascenso más cortos e índices más suaves, consistentes con deflagraciones parciales y bajas masas de eyecta.
   • Anomalías en SNe Ia normales: Los autores identifican SNe Ia normales con tiempos de ascenso inusualmente largos (potencialmente excesos de corta duración) y un subconjunto de eventos tipo 91T de alto estiramiento con ascensos casi lineales ($\alpha \approx 1$). Estos ascendentes lineales parecen ocupar la cola de la distribución de SNe normales luminosas y de evolución lenta, en lugar de formar una subclase distinta, lo que sugiere una mezcla hacia afuera de $^{56}\text{Ni}$ generalizada en la población de alto estiramiento.

4. Comparación de modelos:

   • Modelos cerca-$M_{\text{Ch}}$: Aunque son capaces de producir diversas curvas de luz, los modelos 1D con variaciones de densidad y mezcla de $^{56}\text{Ni}$ no logran reproducir las correlaciones fragmentadas entre los parámetros de tiempo temprano y el estiramiento.
   • Doble detonación sub-$M_{\text{Ch}}$: Los modelos 2D con eyectas asimétricas y efectos de ángulo de visión logran reproducir una correlación fragmentada en la evolución del color ($\alpha_g - \alpha_r$), pero fallan en igualar los tiempos de ascenso observados (prediciendo ascensos más cortos) y no reproducen la ruptura en la relación del índice de ascenso.

Significancia y afirmaciones
El artículo afirma proporcionar la primera caracterización exhaustiva, a nivel de población, de los tiempos tempranos de las SNe Ia utilizando una muestra completa de volumen, yendo más allá de las limitaciones de las muestras pequeñas y sesgadas. Al emplear un marco bayesiano jerárquico, los autores demuestran que los cortes de calidad estrictos pueden introducir sesgos de selección (por ejemplo, sesgando las muestras hacia eventos de alto estiramiento) y que se pueden realizar inferencias de población robustas incluso con datos de muestreo escaso.

La principal significancia reside en el descubrimiento de una dicotomía en el comportamiento de tiempo temprano de las SNe Ia basada en el estiramiento de la curva de luz. La ruptura observada sugiere que los mecanismos físicos que impulsan la emisión temprana difieren fundamentalmente entre las poblaciones de bajo y alto estiramiento. Específicamente, la prevalencia de ascensos suaves y colores azules persistentes en la población de alto estiramiento apunta a una mezcla hacia afuera de $^{56}\text{Ni}$ generalizada en las capas externas de los eyecta.

Los autores concluyen que los modelos teóricos actuales (tanto los de 1D cerca-$M_{\text{Ch}}$ como los de 2D sub-$M_{\text{Ch}}$) son insuficientes para reproducir completamente el espacio de parámetros observado, particularmente las correlaciones específicas con el estiramiento. Argumentan que los modelos de explosión multidimensionales con configuraciones de progenitores más realistas y tratamientos de la distribución de $^{56}\text{Ni}$ y la asimetría de los eyecta son necesarios para explicar la diversidad observada. El artículo posiciona estas restricciones empíricas como un punto de referencia crítico para futuros desarrollos teóricos y destaca el potencial de los próximos grandes sondeos (por ejemplo, el Observatorio Rubin) para refinar aún más estas demografías de población.