ENGRAVE follow-up of a type IIb supernova spatially coincident with the sub-threshold gravitational wave trigger S250818k

Resumen Técnico: Seguimiento ENGRAVE de una supernova de tipo IIb espacialmente coincidente con el disparador de ondas gravitacionales subumbral S250818k

Enunciado del Problema
Durante la cuarta campaña de observación (O4) de la red LIGO-Virgo-KAGRA (LVK), se identificó un candidato de ondas gravitacionales (GW) de baja significancia, S250818k, con una probabilidad de clasificación como fusión de estrellas de neutrones binarias (BNS) del 29% y una región de localización extensa (949 deg² al 90% de credibilidad). El Zwicky Transient Facility (ZTF) descubrió una transición óptica, SN 2025ulz, dentro de esta región aproximadamente 2.7 horas después del disparador. La transición exhibió un decaimiento inicial rápido y un espectro sin características, rasgos que inicialmente imitaban una kilonova (KN) o un escenario de "superkilonova" que involucra estrellas de neutrones con masa subsolar. El desafío científico principal fue clasificar definitivamente SN 2025ulz para determinar si era un contraparte electromagnético genuino del evento de GW o una coincidencia casual, distinguiéndola específicamente de la población contaminante de supernovas jóvenes que pueden imitar las firmas de KN.

Metodología
La Colaboración ENGRAVE (Contrapartes electromagnéticas de fuentes de ondas gravitacionales en el Very Large Telescope) llevó a cabo una extensa campaña observacional multibanda utilizando el Very Large Telescope (VLT) del ESO y otras instalaciones.

• Imagen: Se realizó imagen óptica e infrarroja cercana (NIR) utilizando VLT/FORS2 (bandas R, I, z), VLT/HAWK-I (banda K) y NOT/ALFOSC (banda r). Estas se complementaron con observaciones del Telescopio Espacial Hubble (HST) (F336W, F606W, F110W, F160W) para construir plantillas de referencia profundas para la sustracción de la galaxia anfitriona. La imagen de diferencia se ejecutó utilizando los pipelines ECSNOoPY y HOTPANTS para aislar el flujo de la transición.
• Espectroscopía: Se obtuvo seguimiento espectral utilizando VLT/X-shooter (UVB, VIS, NIR), VLT/FORS2 (rejilla 300I) y VLT/MUSE (Unidad de Campo Integral). Los datos de MUSE fueron cruciales para resolver espacialmente la transición de su galaxia anfitriona, permitiendo la extracción de espectros limpios mediante el modelado y la sustracción del fondo de la anfitriona utilizando múltiples métodos de aperturas y anillos.
• Observaciones de Radio: La fuente fue monitoreada en longitudes de onda de radio utilizando VLA, MeerKAT, uGMRT y e-MERLIN para buscar firmas de resplandor posterior o restricciones sobre la inclinación del chorro.
• Análisis: Se construyeron curvas de luz y se analizó la evolución del color. Se realizó la clasificación espectral utilizando SNID-SAGE. El modelado de la curva de luz comparó dos escenarios: un modelo de kilonova de dos componentes (incluyendo el enfriamiento del capullo) y un modelo de supernova de tipo IIb (enfriamiento por choque más desintegración radiactiva). Se calcularon tasas volumétricas para evaluar la probabilidad de alineación casual.

Resultados Clave

• Clasificación Espectroscópica: Los espectros tempranos (20 de agosto) fueron en gran parte sin características. Sin embargo, los espectros de MUSE obtenidos el 24 y 28 de agosto revelaron un perfil P-Cygni prominente y ancho de Hα y características de Hβ. El análisis SNID-SAGE proporcionó una clasificación de alta confianza de SN 2025ulz como una supernova de tipo IIb (que coincide específicamente con las fases tempranas de SN 2011dh y SN 2016gkg). El corrimiento al rojo de la galaxia anfitriona se refinó a $z = 0.0849 \pm 0.0004$.
• Evolución de la Curva de Luz: La transición mostró un desvanecimiento inicial rápido seguido de un rebrillanteo a un segundo pico a $\sim$15–20 días. Este comportamiento de "cola de enfriamiento por choque" es característico de las SNe de tipo IIb con envolturas extendidas, pero no es predicho por los modelos estándar de kilonova, que típicamente muestran un desvanecimiento monótono.
• Restricciones de Radio: Las observaciones de radio detectaron una densidad de flujo constante consistente con la formación estelar en la galaxia anfitriona (TFS $\sim 1.3 M_\odot$ yr$^{-1}$), sin evidencia de variabilidad o un resplandor posterior transitorio. Los límites superiores descartaron un resplandor posterior similar a GW170817 dentro de ángulos de visión de $\theta \lesssim 12.5^\circ$.
• Modelado:
  • Ajuste de Kilonova: Un modelo de KN de dos componentes ajustó solo los primeros 4 días de datos, pero requirió parámetros no físicos (por ejemplo, masas de eyección y opacidades inconsistentes con las expectativas de fusión BNS) y no logró explicar el rebrillanteo posterior.
  • Ajuste de Supernova: Un modelo de SN de tipo IIb (enfriamiento por choque + desintegración radiactiva) describió con éxito toda la curva de luz hasta 60 días. Los parámetros derivados incluyeron una masa de eyección de $2.2 M_\odot$, una masa de $^{56}$Ni de $0.14 M_\odot$ y un tiempo de explosión aproximadamente 1.6 días antes del disparador de GW.
• Galaxia Anfitriona: La anfitriona es una galaxia con metalicidad moderadamente subsolar ($12 + \log(\text{O/H}) \approx 8.51$) con una tasa de formación estelar de $1.6 \pm 0.5 M_\odot$ yr$^{-1}$.

Conclusiones y Significancia
El artículo concluye que SN 2025ulz es una supernova de tipo IIb que ocurrió en coincidencia temporal y espacial casual con el candidato de GW de baja significancia S250818k. La clasificación es robusta basándose en características espectroscópicas, morfología de la curva de luz y restricciones de tiempo de explosión que son incompatibles con el disparador de GW. El escenario de "superkilonova" está fuertemente desfavorecido por los datos y la navaja de Occam.

Los autores destacan una implicación significativa para las búsquedas futuras de contrapartes de GW: las supernovas de tipo IIb jóvenes son contaminantes principales en las búsquedas de kilonovas. Aproximadamente el 30–50% de las SNe de tipo IIb exhiben colas de enfriamiento por choque que imitan el declive inicial rápido y los espectros sin características de las kilonovas. Los autores estiman una probabilidad de aproximadamente el 30% de que tal intruso aparezca en una región de localización típica de GW (como S250818k) puramente por casualidad. Esto subraya la necesidad de un seguimiento espectroscópico rápido para distinguir contaminantes de kilonovas genuinas, particularmente en la era de encuestas de alta cadencia como Rubin/LSST, donde tales alineaciones casuales se volverán frecuentes. El estudio refuerza que las colas de enfriamiento por choque están entre los contaminantes más prominentes en las búsquedas de kilonovas.