No Evidence for Superradiant Axions in LIGO-Virgo-KAGRA GWTC-5 Binary Black Hole Spins

Resumen Técnico: Sin evidencia de axiones superradiantes en los espines de agujeros negros binarios de LIGO-Virgo-KAGRA GWTC-5

Planteamiento del problema
El axión de la Cromodinámica Cuántica (QCD) y las partículas similares al axión (ALPs) son candidatos a materia oscura bien motivados que resuelven el problema CP fuerte y surgen genéricamente en las compactaciones de la teoría de cuerdas. Si bien los experimentos terrestres están escalando sus búsquedas para axiones con masas $m_a \lesssim 20$ meV, estos tienen dificultades para detectar masas por debajo de aproximadamente $0.5$ neV debido a fuerzas de interacción ultra-débiles. Por el contrario, la superradiancia (SR) inducida por axiones en agujeros negros (BH) ofrece un mecanismo astrofísico prometedor para sondear la parte más baja del rango de masa de axiones teóricamente viable (potencialmente hasta $m_a \sim 10^{-11}$ eV o superior).

El mecanismo de SR postula que cuando la longitud de onda de Compton de un bosón ultraligero es comparable al radio gravitacional de un agujero negro de Kerr en rotación, un proceso de amplificación de ondas clásica forma una nube de tipo hidrógeno. Esta nube extrae el momento angular del BH, agotando su espín hasta un valor crítico $\chi_c(m_a, M)$ hasta que la inestabilidad se detiene. En consecuencia, existe una "región prohibida" en el plano masa-espín (Regge) de los BH donde los agujeros negros no pueden persistir en escalas de tiempo observacionalmente relevantes; cualquier BH nacido con un espín por encima de la curva crítica debería reducir su espín rápidamente.

Las restricciones previas sobre axiones en este rango de masa se han derivado de mediciones de rayos X de espines de agujeros negros de masa estelar y supermasivos. Sin embargo, estas restricciones están fundamentalmente limitadas por sustanciales sistemáticas de modelado inherentes a los métodos de ajuste de continuo y espectroscopía de reflexión de la línea K$\alpha$ del hierro, que requieren un modelado detallado de los discos de acreción internos. Análisis recientes sugieren que algunos BH de alto espín que impulsan las exclusiones basadas en rayos X pueden sufrir de un modelado sistemático erróneo, lo que podría sesgar las mediciones de espín hacia valores más altos.

Metodología
Este trabajo realiza un análisis bayesiano jerárquico para restringir axiones ultraligeros utilizando las mediciones de espín de agujeros negros binarios (BBH) del catálogo LVK GWTC-5. Este catálogo es el catálogo público de BBH más extenso hasta la fecha, conteniendo $N=257$ eventos de fusión con masas de BH que abarcan aproximadamente $5\text{--}135\,M_\odot$ desde las corridas de observación O1 hasta O4b.

El análisis procede de la siguiente manera:

1. Modelado de la población: La distribución natal de espín de los BH se modela como una distribución Beta, $p(\chi | \alpha, \beta) \propto \chi^{\alpha-1}(1-\chi)^{\beta-1}$, donde $\alpha$ y $\beta$ son hiperparámetros. Esto permite la flexibilidad para capturar poblaciones que van desde picos de alto espín hasta picos de bajo espín (esto último favorecido por las inferencias actuales de LVK).
2. Física de la superradiancia: Para una masa de axión dada $m_a$, el modelo incorpora el efecto de reducción de espín por SR. Los BH nacidos con espines por encima de la curva crítica $\chi_{\text{sat}}(m_a, M, \tau_m)$ son reducidos a la curva de saturación, mientras que aquellos por debajo retienen sus espines natales. El tiempo entre la formación del BH y la fusión, $\tau_m$, se trata como un parámetro de molestia extraído de una distribución de tiempo de retraso (DTD) log-uniforme $p(\tau_m) \propto \tau_m^{-1}$ sobre $[10^7, \tau_H]$ años, contabilizando tanto escenarios de formación aislada como de ensamblaje dinámico.
3. Verosimilitud jerárquica: La posterior conjunta sobre los hiperparámetros de la población $\Lambda = (m_a, \alpha, \beta, \tau_{m,1}, \dots)$ se computa utilizando las muestras de la estimación de parámetros (PE) de LVK para cada evento. El análisis preserva las correlaciones por evento entre las masas y los espines de los componentes primario y secundario.
4. Prueba de hipótesis: El estudio compara la hipótesis de señal de axión ($H_B$) contra la hipótesis nula ($H_A$, sin axión). El análisis incluye modos hasta $m=6$ y utiliza soluciones numéricas para las ecuaciones de SR en acoplamientos más altos ($\alpha \gtrsim 0.5$) para asegurar la robustez, aunque las soluciones perturbativas dominan el análisis fiduciario.
5. Selección de datos: Los eventos se seleccionan con masas componentes medianas $M_{1,2}^{\text{med}} \geq 3\,M_\odot$ para excluir candidatos a estrellas de neutrones, resultando en la muestra final de 257 BBHs.

Resultados clave

• Límites de exclusión: El análisis no encuentra evidencia de axiones a través de más de dos décadas de masa. Con un 95% de confianza, el estudio excluye masas de axiones en el rango:
  $$1.7 \times 10^{-14} \, \text{eV} \lesssim m_a \lesssim 3.3 \times 10^{-12} \, \text{eV}$$
  Esto corresponde a un rango de constante de decaimiento de aproximadamente $10^{17} \text{--} 10^{19}$ GeV para el axión de la QCD.
• Evidencia Bayesiana: El factor de Bayes entre la hipótesis de señal de axión y la hipótesis nula es $\log B^B_A = 0.22 \pm 0.06$, lo cual es consistente con la ausencia de evidencia para la hipótesis del axión.
• Inferencia de la población: Simultáneamente, el análisis infiere la población de espín natal de los BH, obteniendo medianas posteriores de $\alpha = 0.69^{+0.21}_{-0.12}$ y $\beta = 2.28^{+0.62}_{-0.43}$. Esto favorece una población de espín natal bajo ($\alpha < \beta$), consistente con las inferencias de población de LVK en general.
• Robustez: La región de exclusión es robusta a las variaciones en el prior sobre $m_a$ (incluyendo un prior de Jeffreys no informativo), la distribución de tiempo de retraso y modelos alternativos de parametrización de espín natal (por ejemplo, modelos de mezcla Beta o Gaussianos truncados).

Significancia y afirmaciones
Los autores afirman que este resultado representa uno de los límites inferiores robustos más fuertes para la masa del axión de la QCD hasta la fecha. A diferencia de las restricciones previas basadas en rayos X, este análisis depende de datos de ondas gravitacionales donde la inferencia de espín se deriva únicamente de las formas de onda de inspiración-fusión-amortiguamiento (inspiral-merger-ringdown), evitando las sistemáticas de modelado de discos de acreción que plagan las mediciones de rayos X.

Al excluir el rango de masa $1.7 \times 10^{-14} \text{ eV} \lesssim m_a \lesssim 3.3 \times 10^{-12} \text{ eV}$, este trabajo elimina efectivamente la parte de menor masa del espacio de parámetros del axión de la QCD que es objetivo de experimentos terrestres como CASPEr. Además, los resultados restringen las partículas similares al axión y los escalares ultraligeros en este rango de masa, impactando las construcciones de la "axiverse" de cuerdas que predicen masas de partículas distribuidas casi log-uniformemente.

El artículo señala que se espera que estas restricciones mejoren sustancialmente con las futuras corridas de observación de LVK (O4c, O5) y los detectores de tercera generación (Cosmic Explorer, Einstein Telescope), que proporcionarán inferencias de espín natal más estrictas, mayores relaciones señal-ruido y acceso a un rango de masa más amplio, potencialmente extendiendo la sensibilidad a masas de axiones por encima de $10^{-11}$ eV.