Probing globular clusters parameters through gravitational wave lensing with stellar-mass black hole binaries

Resumen Técnico: Sondeo de los Parámetros de Cúmulos Globulares mediante Lentes de Ondas Gravitacionales con Binarias de Agujeros Negros de Masa Estelar

Planteamiento del Problema
Los cúmulos globulares (CG) son sistemas estelares densos que sirven como laboratorios para estudiar la dinámica estelar y la formación temprana de galaxias. Un parámetro crítico para comprender su dinámica interna, formación y evolución es su masa, la cual se infiere típicamente de la dispersión de velocidad estelar central ($\sigma_v$). Los métodos actuales para determinar la masa de los CG dependen de observaciones electromagnéticas (EM), como mediciones de la dispersión de velocidad y perfiles de brillo superficial. Estos métodos enfrentan limitaciones significativas, incluyendo efectos de proyección en la línea de visión, hacinamiento estelar en regiones centrales densas, ruido de disparo (shot noise) de estrellas gigantes luminosas y la incertidumbre en la identificación del centro real del cúmulo. Estos desafíos se ven exacerbados en cúmulos con colapso de núcleo o sistemas que potencialmente albergan agujeros negros de masa intermedia (IMBH). Además, aunque los próximos sondeos de campo amplio (p. ej., LSST, Euclid) ampliarán los catálogos de CG con posiciones y distancias precisas, la medición de la cinemática interna como la dispersión de velocidad sigue siendo difícil y típicamente requiere un seguimiento espectroscópico dedicado.

Los autores proponen una sonda complementaria: utilizar ondas gravitacionales (GW) de coalescencias de binarias compactas (CBC) que son moduladas por lentes de CG. A diferencia del lente fuerte por galaxias, que produce múltiples imágenes, el lensing por CG opera en el régimen de óptica de ondas. En este régimen, el lente no crea múltiples señales, sino que imprime una amplificación dependiente de la frecuencia en la señal de GW. El artículo investiga si estas firmas de óptica de ondas pueden utilizarse para recuperar las propiedades intrínsecas del CG que actúa como lente, específicamente su dispersión de velocidad central.

Metodología
El estudio emplea un marco de simulación y estimación bayesiana de parámetros para probar la recuperabilidad de los parámetros de los CG a partir de señales de GW con efecto de lente.

1. Modelado del Lente: Los CG se modelan como Esferas Isotérmicas Singulares (SIS). La distribución de masa se aproxima mediante un perfil de densidad $\rho(r) = \sigma_v^2 / (2\pi G r^2)$, donde se asume que la dispersión de velocidad $\sigma_v$ es constante. La masa de lensing efectiva ($M_{SIS}$) se deriva como una función de $\sigma_v$ y la distancia del lente ($D_l$).
2. Simulación de la Señal: Los autores simulan señales de GW similares a GW150914 (una fusión de binaria de agujeros negros de masa estelar). Estas señales son moduladas por el factor de amplificación SIS derivado del integral de difracción de Fresnel-Kirchhoff. La amplificación depende de la frecuencia adimensional ($w$) y el parámetro de impacto ($y$).
3. Estudio de Inyección: Las señales con efecto de lente se inyectan en ruido gaussiano generado a partir de las densidades espectrales de potencia (PSD) de la cuarta observación (O4) de una red de tres detectores (LIGO Hanford, LIGO Livingston y Virgo).
4. Estimación de Parámetros: Las señales inyectadas se recuperan utilizando inferencia bayesiana (a través del paquete Bilby y el muestreo anidado dynesty) con plantillas de forma de onda que incluyen tanto parámetros de la fuente (15 parámetros binarios utilizando el aproximante IMRPhenomXPHM) como parámetros del lente (masa de lensing SIS $M_{SIS}$ y parámetro de impacto $y$).
5. Identificación de Candidatos: A partir de las distribuciones posteriores recuperadas, se construye una región de localización en el cielo con un crédito del 90%. Esta región se cruza con un catálogo de CG conocidos de la Vía Láctea (de [37]) para identificar el posible lente.
6. Recuperación de la Dispersión de Velocidad: Una vez identificado un CG candidato, se utiliza su distancia conocida ($D_l$) del catálogo para invertir la posterior de la masa de lensing recuperada ($M_{SIS}$) para estimar la dispersión de velocidad ($\sigma_v$) del cúmulo utilizando la relación $M_{SIS} \propto \sigma_v^4 D_l$.

Se seleccionaron tres CG específicos del catálogo de la Vía Láctea para el estudio basados en sus distintas dispersiones de velocidad central: Liller-1 ($\sigma_v \approx 19.6$ km s$^{-1}$), NGC 6093 ($\sigma_v \approx 11.1$ km s$^{-1}$) y NGC 2419 ($\sigma_v \approx 5.1$ km s$^{-1}$). Para cada cúmulo, se realizaron 20 inyecciones con variaciones en el parámetro de impacto ($y \in [0.001, 5]$).

Resultados Clave

• Recuperación de la Masa de Lensing y la Dispersión de Velocidad: Para alineaciones fuente-lente favorables (parámetros de impacto pequeños), los valores inyectados de la dispersión de velocidad fueron recuperados exitosamente dentro de los intervalos de crédito del 68%. Para Liller-1, el valor mediano recuperado fue de 19.7 km s$^{-1}$, coincidiendo estrechamente con el valor inyectado de 19.6 km s$^{-1}$.
• Impacto del Parámetro de Impacto: La precisión de la recuperación depende fuertemente del parámetro de impacto $y$. A medida que $y$ aumenta, la amplificación disminuye, debilitando la impronta de la lente en la señal. Esto conduce a un muestreo más pobre de la masa de lensing y a una mayor incertidumbre en la dispersión de velocidad inferida. Para valores de $y$ más grandes, los valores medianos recuperados pueden desplazarse del valor real debido a las realizaciones del ruido.
• Localización en el Cielo y Cruce de Catálogos: El cúmulo inyectado fue correctamente identificado dentro de la región de localización en el cielo con un crédito del 90% en todos los casos probados. El método distinguió con éxito el cúmulo de lensing de otros en el catálogo, siempre que la fuente estuviera posicionada detrás del cúmulo objetivo específico.
• Degeneración de Parámetros: El estudio observó una anticorrelación entre la masa de lensing recuperada ($M_{SIS}$) y el parámetro de impacto ($y$) para $y \sim 1$, consistente con hallazgos previos para lentes de masa puntual.

Significancia y Reivindicaciones
El artículo demuestra que las GW con efecto de lente pueden servir como una sonda complementaria para la dinámica de los CG, permitiendo específicamente la estimación de la dispersión de velocidad y, por extensión, la masa del cúmulo. Esto ofrece un nuevo método para superar las limitaciones inherentes a la determinación de masa basada en EM, particularmente en regiones densas o de colapso de núcleo.

Los autores señalan que la probabilidad de que un CG de la Vía Láctea actúe como lente para un evento de GW es baja ($\tau \sim 10^{-8}$), lo que destaca el potencial de los detectores de próxima generación para descubrir tales eventos. La metodología se presenta como una herramienta para:

1. Identificar con confianza eventos de GW como señales con efecto de lente si un CG de la Vía Láctea es el lente.
2. Restringir la dispersión de velocidad y la masa del cúmulo de lensing.
3. Indicar potencialmente la presencia de IMBHs en los núcleos de los cúmulos.

El estudio es modesto en su alcance, utilizando el modelo SIS y centrándose en los CG de la Vía Láctea donde las distancias están bien conocidas. Los autores reconocen que, para los CG extragalácticos (donde los catálogos pueden estar incompletos), existiría una degeneración entre la distancia del lente y la dispersión de velocidad, lo que requeriría una estimación de parámetros más extensa. También señalan que, aunque el presente trabajo utiliza el modelo SIS, la metodología podría extenderse a modelos más realistas (p. ej., esferas isotérmicas con núcleo, modelos de King) en trabajos futuros. El artículo concluye que una red global de al menos tres detectores de GW con sensibilidad comparable es esencial para mejorar la precisión de la localización en el cielo necesaria para que esta técnica sea efectiva.