Implications for Primordial Black Hole Dark Matter from a Single Subsolar Mass Gravitational-wave Detection in LVK O1--O4

Resumen Técnico: Implicaciones para la Materia Oscura de Agujeros Negros Primordiales a partir de una Única Detección de Ondas Gravitacionales de Masa Subsolar en LVK O1–O4

Enunciado del Problema
La detección de fusiones de binarias compactas con masas de componente en el régimen subsolar ($M \lesssim 1 M_\odot$) presenta un desafío significativo para la astrofísica estelar estándar, ya que los mecanismos convencionales de colapso estelar no producen agujeros negros en este rango de masas. Tales detecciones sugerirían fuertemente un origen no estelar, más notablemente Agujeros Negros Primordiales (ANP), los cuales podrían constituir una fracción de la Materia Oscura (MO) del Universo. El 12 de noviembre de 2025, la colaboración LIGO-Virgo-KAGRA (LVK) reportó un evento candidato, S251112cm, con una masa de chirp ($M_c$) en el rango $[0.1, 0.87] M_\odot$ y una probabilidad $>99\%$ de que el componente más ligero sea subsolar. Este artículo investiga si una población de ANP formados durante la época de la Cromodinámica Cuántica (QCD), caracterizada por una función de masa extendida, puede cuantitativamente dar cuenta de la tasa de fusiones observada de tal evento.

Metodología
Los autores modelan la población de ANP utilizando una función de masa extendida ($df_{PBH}/d\ln M$) derivada de G. Hasinger (2020) y modificada por asimetrías de sabor leptónico según lo descrito en D. Bödeker et al. (2021). Este espectro abarca desde $10^{-6} M_\odot$ hasta $4 \times 10^8 M_\odot$, con picos específicos correspondientes a transiciones de fase en el Universo temprano. El modelo asume una fracción de materia oscura de ANP de $f_{PBH} = 0.339$ dentro del rango de masas considerado.

El cálculo de la tasa de fusiones se centra en el canal de formación de "época tardía", donde los ANP son capturados gravitacionalmente dentro de halos de materia oscura después de la formación de estructuras, en lugar del canal de desacoplamiento del Universo temprano. La densidad de tasa de fusión intrínseca ($\tau_L$) se deriva basándose en la densidad local de materia oscura ($\rho_{loc} \approx 0.4 \text{ GeV cm}^{-3}$), la dispersión de velocidad de los ANP (asumiendo que siguen una distribución Maxwell-Boltzmann con $v_{vir} \approx 250 \text{ km s}^{-1}$) y la sección transversal de captura de ondas gravitacionales.

Para determinar la tasa de eventos detectables ($R_{tot}$), los autores integran la tasa intrínseca sobre la función de masa de los ANP, aplicando cortes de selección consistentes con la alerta de S251112cm: una relación de masas $q \in [0.05, 1]$ y una masa de chirp $M_c \in [0.1, 0.87] M_\odot$. La distancia de horizonte del detector se calcula utilizando la densidad espectral de amplitud H1 de O3b, teniendo en cuenta las fases de espiral y fusión de la forma de onda. El análisis asume la cosmología de Planck 2018.

Resultados Clave

1. Tasa de Eventos Predicha: Para la función de masa extendida de ANP adoptada y los parámetros de halo estándar, el modelo predice una tasa de fusiones subsolares detectables de $R_{tot} = 0.8 \text{ yr}^{-1}$ para un detector similar a LVK (sensibilidad O3b). Esto corresponde a una densidad de tasa de fusión local intrínseca de aproximadamente $3.45 \text{ Gpc}^{-3} \text{ yr}^{-1}$ dentro del espacio de parámetros analizado.
2. Consistencia con la Observación: La tasa observada, inferida a partir de la única detección de S251112cm durante un tiempo de observación efectivo de 4.35 años, es $R_{obs} = 0.23^{+0.86}_{-0.218} \text{ yr}^{-1}$ (IC 95%). La tasa predicha de $0.8 \text{ yr}^{-1}$ cae bien dentro de este intervalo de confianza.
3. Restricciones sobre la Abundancia de ANP: Al escalar la tasa predicha con la abundancia de ANP ($R_{det} \propto f_{PBH}^2$), los autores derivan un límite inferior conservador sobre la fracción de ANP. Utilizando el límite inferior de Poisson al 95% para la tasa media real de eventos ($\lambda_{low} = 0.0513$), la detección implica $f_{PBH} > 0.04$ para el modelo adoptado.
4. Consistencia con Masas Estelares: Cuando la misma metodología se aplica al rango de masas estelares ($3 - 200 M_\odot$), el modelo predice una tasa de detección de $\sim 120 \text{ yr}^{-1}$. Esto se alinea con las expectativas actuales de LVK y sugiere que una fracción no despreciable de las fusiones observadas de agujeros negros binarios de masa estelar también podrían originarse de la misma población primordial.

Significado y Afirmaciones
El artículo afirma que si S251112cm se confirma como un evento astrofísico, sirve como un candidato convincente para la primera detección de una binaria de ANP de masa subsolar en fusión. La ausencia de un contraparte electromagnético y la baja masa de los componentes son difíciles de conciliar con los canales de formación estelar estándar (tales como la fragmentación de discos o la fisión de núcleos en supernovas), los cuales son teóricamente posibles pero estadísticamente menos probables dada la falta de transitorios asociados.

Los autores enfatizan que su análisis proporciona un escenario "motivado físicamente" donde los ANP formados durante la época de la QCD pueden dar cuenta de la señal observada sin violar las restricciones existentes de acreción en el Fondo Cósmico de Microondas (CMB), microlente o ondas gravitacionales. Si bien la tasa de captura de época tardía conlleva incertidumbres astrofísicas significativas (dependientes de la densidad del halo y las distribuciones de velocidad), el acuerdo entre las tasas predichas y observadas sugiere que los ANP podrían constituir una fracción significativa de la Materia Oscura ($f_{PBH} \sim 0.339$ en el rango de masas del modelo). El artículo concluye que las futuras campañas de observación de LVK (específicamente O5) proporcionarán las estadísticas necesarias para probar esta hipótesis y potencialmente revelar la naturaleza de la Materia Oscura, transitando la región de masa subsolar de una curiosidad teórica a una frontera observacional.