Wave-optics imprints of dark matter subhalos on strongly lensed gravitational waves. II. Saddle images and detectability

Resumen Técnico: Improntas de la óptica de ondas de subhalos de materia oscura en ondas gravitacionales fuertemente lentes. II. Imágenes de silla y detectabilidad

Planteamiento del problema
La abundancia y la estructura de densidad de los subhalos de materia oscura de baja masa ($10^4$–$10^7 M_\odot$) sirven como un discriminante primario entre el modelo de Materia Oscura Fría (CDM) y modelos alternativos. Sin embargo, la ausencia de bariones en estos subhalos los hace invisibles para las observaciones electromagnéticas. Mientras que las ondas gravitacionales (GW) ofrecen una sonda de óptica de ondas (WO) de la subestructura mediante distorsiones dependientes de la frecuencia en el factor de amplificación, el trabajo previo (Papel I) se limitó a la paridad de imagen "mínima" en sistemas de lente fuerte. Este artículo aborda dos preguntas abiertas: (1) ¿Cómo impronta los subhalos la paridad de la imagen de "silla" (saddle), la cual es conocida por ser frágil en la óptica geométrica (GO) respecto a las anomalías de la relación de flujo? (2) ¿Es el señal combinado de ambas paridades de imagen detectable por la LISA (Laser Interferometer Space Antenna) al marginalizar sobre parámetros desconocidos de la fuente y la lente?

Metodología
Los autores extienden el formalismo de WO a la imagen de punto de silla para una lente de escala galáctica (halo NFW + esfera isoterma singular) en $z_L=0.5$ lenteando un binario de agujeros negros masivos (MBHB) en $z_S=1.5$.

1. Cálculo del punto de silla: Evaluar el factor de amplificación en una silla es numéricamente delicado porque los contornos de tiempo de llegada iguales son abiertos (hiperbólicos) en lugar de cerrados, y la señal del subhalo es un residuo pequeño sobre un fondo divergente. Los autores desarrollan un método dedicado en el dominio del tiempo:

   • Separan la señal en un template macro de GO analítico (la silla cuadrática pura) y un residuo de WO de subhalo transformado numéricamente.
   • El residuo se calcula utilizando una cuadratura de Filon para integrales altamente oscilatorias en una rejilla no uniforme agrupada cerca del retardo de la silla.
   • Se añade una cola de campo lejano analítica para manejar el potencial logarítmico no cancelable de los subhalos NFW truncados en retardos grandes.
   • Este método se valida contra un árbitro del método de área sin contornos independiente.

2. Ensamble de Monte Carlo: Los autores generan 1,000 realizaciones de poblaciones de subhalos CDM (masas $10^2$–$10^9 M_\odot$) alrededor de imágenes tanto mínimas como de silla. Computan los factores de amplificación por imagen $F(w)$ y la modulación compleja $\eta(w) = F(w)/F_{GO}(w) - 1$.

3. Marco de Detectabilidad: Dado que LISA observa la forma de onda lente total (la suma coherente de la imagen mínima y la de silla), los autores construyen un estadístico de detección que:

   • Desmodula la interferencia de la imagen macro (separando las dos imágenes en el dominio del tiempo o del cepstrum) para recuperar las funciones de transferencia por imagen.
   • Proyecta fuera los parámetros degenerados (amplitud, fase de Morse, tiempo de llegada) utilizando un producto interno de filtro adaptado (matched-filter). Esto aísla la estructura dependiente de la frecuencia no absorbible de la señal del subhalo.
   • Combina la significación de las imágenes mínima y de silla en cuadratura.

Resultados Clave

• Asimetría de Paridad en la Desmagnificación Media: Aunque la amplitud de fluctuación de la impronta del subhalo es ciega a la paridad (afectando a ambas imágenes a un nivel porcentual similar), la modulación media se divide por paridad. La imagen mínima tiene una magnificación neta ($dc \approx +1.3\%$), mientras que la imagen de silla tiene una desmagnificación neta ($dc \approx -0.8\%$). Esta división sistemática se correlaciona con el cizallamiento (shear) de marea local con signos opuestos para las dos paridades, confirmando que la "fragilidad de la silla" de la óptica geométrica persiste en el régimen de óptica de ondas como un efecto de (des)magnificación media.
• Características de la Modulación: Los subhalos inducen modulaciones de amplitud del orden de un porcentaje y desplazajes de fase de $\sim 10^{-2}$ rad a través de la banda de LISA ($10^{-4}$–$10^{-1}$ Hz). La modulación está dominada por subhalos en el rango de $10^4$–$10^7 M_\odot$.
• Significación de la Detección: Para una fuente MBHB fiduciaria de masa total $\sim 10^6 M_\odot$ ubicada cerca de la cáustica de la lente ($y_{src} \lesssim 0.1$), la significación de detección combinada supera los $5\sigma$ en el 62% de las realizaciones de Monte Carlo. La imagen mínima contribuye más a la significación que la de silla, pero incluir la de silla eleva significativamente la tasa de detección combinada (del 40% para la mínima sola al 62% combinado para el caso de $10^6 M_\odot$).
• Dependencia de la Posición de la Fuente: La significación de detección escala estrepitosamente con el parámetro de impacto de la fuente ($\rho_{sub} \propto y_{src}^{-1.6}$). Una detección robusta por evento requiere que la fuente se encuentre muy cerca de la cáustica ($y_{src} \lesssim 0.1$).

Significancia y Reivindicaciones
El artículo afirma que las ondas gravitacionales fuertemente lentes proporcionan una sonda sensible y complementaria de la subestructura de materia oscura en escalas ($10^4$–$10^7 M_\odot$) inaccesibles para las observaciones electromagnéticas.

• Complementariedad: A diferencia de las anomalías de relación de flujo en cuásares que dependen del brillo de la silla, este método detecta la distorsión dependiente de la frecuencia de toda la forma de onda, permitiendo medir incluso la fragilidad de la silla cuando la señal no es anomalamente tenue en un sentido estático.
• Pronósticos de Población: Al plegar las significaciones por evento con pronósticos optimistas de tasas de lente fuerte, se sugieren 10–20 detecciones de subestructura durante la misión de LISA. Sin embargo, los autores señalan que poblaciones de fuentes y lentes conservadoras podrían reducir este rendimiento a bien menos de un evento.
• Poder Diagnóstico: La señal se ve potenciada para subestructuras más densas (por ejemplo, agujeros negros primordiales o materia oscura auto-interactuante de núcleo colapsado), ofreciendo una ruta para restringir no solo la abundancia, sino también la densidad interna de la subestructura oscura.

Los autores concluyen que, si bien el número absoluto de detecciones depende fuertemente de las inciertas tasas astrofísicas, la sensibilidad por evento es robusta: incluso un solo evento fuerte y fuertemente lente cerca de una curva crítica puede servir como un poderoso diagnóstico de la estructura de la materia oscura subgaláctica.