Improving low-latency multi-messenger follow-up of neutron star-black hole mergers with mode-by-mode filtering

Resumen Técnico: Mejora del seguimiento multi-mensajero de baja latencia de fusiones de estrella de neutrones–agujero negro mediante el filtrado modo a modo

Planteamiento del problema
La estimación rápida de parámetros (PE, por sus siglas en inglés) para las fusiones de estrella de neutrones–agujero negro (NSBH) es crítica para guiar las observaciones de seguimiento electromagnético (EM). Los pipelines actuales de baja latencia, como BAYESTAR, típicamente modelan solo el armónico cuadrupolar dominante de ondas gravitacionales (GW), $(\ell, |m|) = (2, 2)$. Para binarias asimétricas como las NSBH (donde las razones de masa son $q \gtrsim 5\text{--}10$), este enfoque deja fuertes degeneraciones entre la distancia de luminosidad ($d_L$), la inclinación ($\iota$) y los parámetros intrínsecos de la binaria. Específicamente, en las formas de onda de solo cuadrupolo, $d_L$ se correlaciona fuertemente con $\iota$, y la razón de masa se correlaciona con el espín efectivo. Estas degeneraciones ensanchan los volúmenes de localización, debilitan las restricciones sobre la naturaleza del objeto secundario (NS vs. BH) y degradan las estimaciones de inclinación, lo que dificulta las decisiones sobre la detectabilidad EM y la asociación con galaxias anfitrionas. Si bien los modos de orden superior (HMs) pueden romper estas degeneraciones, el muestreo estocástico completo incluyendo HMs es computacionalmente prohibitivo para aplicaciones de baja latencia, requiriendo a menudo horas o días.

Metodología
Los autores proponen un método de PE de baja latencia que realiza un filtrado modo a modo de la serie temporal de la relación señal-ruido (SNR) para los armónicos $(2, 2)$, $(3, 3)$ y $(4, 4)$. El enfoque utiliza la librería flywheel, una adaptación de baja latencia del código cogwheel, para marginalizar analítica o semi-analíticamente sobre los parámetros que entran de forma simple en la verosimilitud (por ejemplo, $d_L$ y la fase de referencia) mientras utiliza el muestreo de importancia adaptativo para los parámetros extrínsecos restantes.

El algoritmo toma tres entradas principales:

1. $\langle d | h_{\ell m} \rangle$: La serie temporal de la SNR para cada plantilla de modo.
2. $\langle h_{\ell m} | h_{\ell' m'} \rangle$: La matriz de covarianza de los modos.
3. Muestras de la distribución previa para las razones de SNR de los modos ($R_{\ell m}$) derivadas de una distribución astrofísica asumida.

El método marginaliza sobre los parámetros extrínsecos de la binaria y el subconjunto de parámetros intrínsecos que afectan la amplitud de los modos. La integral de la distancia de luminosidad se evalúa mediante la interpolación de una tabla precomputada, la integral de la fase de referencia mediante la cuadratura del trapecio, y las integrales restantes (inclinación, ubicación en el cielo, polarización) mediante el muestreo de importancia adaptativo. Las muestras de la posterior se generan mediante el re-muestreo de importancia de las muestras de cuasi-Monte Carlo utilizadas en la integral.

Resultados clave
El método fue validado utilizando dos conjuntos de datos: detecciones simuladas de NSBH en una red de LIGO-Virgo con la sensibilidad de diseño de O5 y datos públicos de eventos NSBH detectados previamente.

• Detecciones simuladas de O5: Aplicado a $\sim 1.6 \times 10^4$ eventos simulados con SNR $\ge 8$, la inclusión de HMs mejoró significativamente las restricciones en comparación con los análisis de solo $(2, 2)$:

  • Distancia y Volumen: El método redujo el ancho del intervalo de credibilidad del 90% para la distancia de luminosidad y el volumen de localización comóvil por factores que a menudo superan 2. Para un evento similar a GW190814, el volumen comóvil se redujo por un factor de $\sim 4.2$.
  • Inclinación y Ángulo de Visión: Los HMs rompieron la degeneración $d_L$–$\iota$, mejorando las restricciones del ángulo de visión hasta en un factor de 5 para sistemas de cara al frente/atrás (face-on/off). Esto impacta directamente en la probabilidad de detectar brotes gamma de rayos X (GRBs) de chorros relativistas.
  • Masa Secundaria: Las restricciones sobre la masa secundaria en el marco de la fuente ($m_{src,2}$) mejoraron, particularmente para sistemas más pesados. Para sistemas de alta masa, la probabilidad de que la secundaria sea una estrella de neutrones ($m_{src,2} \le 2.7 M_\odot$) aumentó hasta en un 50%, ayudando en las estrategias de seguimiento de kilonovas.
  • Costo Computacional: La marginalización de modos superiores toma aproximadamente $0.13 \pm 0.03$ segundos por evento en un solo núcleo de CPU, con $0.11$ segundos adicionales para extraer 1000 muestras de la posterior. Esto es comparable a los pipelines actuales de solo cuadrupolo.

• Validación con Datos Reales:

  • GW190814: Como un evento de alta SNR ($\sim 25$) y alta razón de masa ($q \sim 9$), GW190814 mostró una firma clara de modos de orden superior. La aplicación del método a datos reales estrechó la posterior de $d_L$ por un factor de 1.6, la de la inclinación por un factor de 2, y la del volumen de localización por un factor de 2.1, acercando la posterior de baja latencia a los resultados de la PE bayesiana completa.
  • Otras Eventos NSBH: Para las otras tres detecciones confirmadas de NSBH (GW200105, GW200115, GW230529), que tienen SNRs más bajas ($\sim 11\text{--}14$) y masas menores, no se observó una mejora significativa, lo cual es consistente con su contenido de modos subdominantes más débil.

Significancia y Afirmaciones
El artículo afirma que incluir modos de orden superior en la PE rápida produce mejoras sustanciales en la localización y en las mediciones de la masa secundaria con un costo computacional comparable a los métodos actuales de solo $(2, 2)$. Los autores argumentan que estas mejoras son cruciales para las estrategias de seguimiento EM, permitiendo:

1. Una identificación más eficiente de las galaxias anfitrionas a partir de catálogos.
2. Una mejor determinación de los ángulos de visión para evaluar la detectabilidad de la emisión de chorro (GRBs/Rayos X).
3. Una mejor evaluación de la detectabilidad de kilonovas basada en el potencial de disrupción de marea.
4. Una clasificación de la fuente (NS vs. BH) más confiable para activar los protocolos de seguimiento apropiados.

Los autores concluyen que su método está listo para su despliegue en las próximas corridas de observación (específicamente O5) para generalizar las búsquedas actuales de solo cuadrupolo de baja latencia. Observan que, aunque el método actualmente asume binarias no precesantes, la incorporación de la precesión es una extensión natural para trabajos futuros. El enfoque está diseñado para ser agnóstico al pipeline, requiriendo solo las series temporales de la SNR modo a modo y su covarianza, que pueden ser generadas por cualquier pipeline de búsqueda.