Inference on inner galaxy structure via gravitational waves from supermassive binaries

Resumen Técnico: Inferencia sobre la estructura del centro galáctico mediante ondas gravitacionales de binarias supermasivas

Planteamiento del problema
La detección de un fondo estocástico de ondas gravitacionales (SGWB) por parte de las redes de temporización de púlsares (PTAs), específicamente el conjunto de datos de 15 años de NANOGrav, indica una población de binarias de agujeros negros supermasivos (SMBHB). Aunque el espectro observado generalmente se alinea con las predicciones para binarias circulares impulsadas por la emisión de ondas gravitacionales ($h_c \propto f^{-2/3}$), los datos exhiben una ligera preferencia por un giro espectral en las frecuencias más bajas observadas (alrededor de 4 nHz). Esta característica sugiere una aceleración en el endurecimiento orbital antes del régimen dominado por ondas gravitacionales, lo que potencialmente resuelve el "problema del parsec final". El origen físico de este giro se hipotetiza como interacciones ambientales —específicamente, la expulsión de estrellas y partículas de materia oscura mediante lanzamientos gravitacionales de tres cuerpos— que extraen la energía orbital de manera más eficiente que la fricción dinámica de dos cuerpos. Sin embargo, las propiedades específicas de la distribución de materia en el centro galáctico (densidad y pendiente) que impulsan este proceso, y su interacción con la excentricidad orbital inicial, permanecen sin constreñir.

Metodología
Los autores modelan la coevolución de los parámetros orbitales de las SMBHB (semieje mayor $a$ y excentricidad $e$) y el perfil de densidad de la materia circundante. El análisis integra dos mecanismos primarios:

1. Emisión de ondas gravitacionales: Gobierna el decaimiento orbital y la circularización, particularmente a altas frecuencias.
2. Dispersión de tres cuerpos: Describe el endurecimiento de la órbita de la binaria debido a las interacciones con el fondo de estrellas y materia oscura. La tasa de endurecimiento se parametriza mediante la densidad de materia local $\rho$ y la dispersión de velocidades $\sigma$.

El estudio asume un perfil de densidad de ley de potencia en el centro galáctico previo al barrido (scouring): $\rho(r) = \rho_{pc} (r/1 \text{ pc})^{-\gamma}$, donde $\rho_{pc}$ es la normalización de la densidad y $\gamma$ es la pendiente radial. El radio de influencia $r_i$ se define donde la masa encerrada es igual a dos veces la masa de la binaria. Los autores calculan la frecuencia de giro $f_t$ donde la dispersión de tres cuerpos transiciona al dominio de las ondas gravitacionales, la cual depende de $\rho_i$, $\sigma_i$, la masa de la binaria $M$ y la relación de masas $q$.

Para inferir las restricciones, los autores:

• Adoptaron el conjunto de datos de 15 años de NANOGrav, centrándose en los cinco intervalos de frecuencia más bajos (2–10 nHz) donde la relación señal-ruido es robusta.
• Modelaron la población de SMBHB utilizando una distribución astrofísica fiduciaria (masa, corrimiento al rojo, relación de masas) derivada de análisis previos de NANOGrav.
• Realizaron una inferencia bayesiana para constreñir el espacio de parámetros conjunto de la excentricidad inicial ($e_0$), la normalización de la densidad ($\rho_{pc}$) y la pendiente ($\gamma$).
• Probaron parametrizaciones alternativas, incluyendo distribuciones de ley de potencia para la excentricidad inicial y distribuciones gaussianas para la varianza de la densidad, para asegurar la robustez.

Resultos clave

• Restricciones de densidad: El análisis favorece una densidad de materia en la escala de un pársec de aproximadamente $\rho_{pc} \sim 10^6 M_\odot \text{pc}^{-3}$ en la mayor parte del espacio de parámetros. La distribución posterior para $\log_{10}(\rho_{pc}/10^5 M_\odot \text{pc}^{-3})$ se estima en $0.6^{+1.4}_{-2.6}$ (1$\sigma$).
• Degeneración excentricidad-densidad: Existe una degeneración entre la excentricidad inicial ($e_0$) y la densidad ($\rho_{pc}$). Una densidad menor requiere una excentricidad inicial extremadamente alta ($e_0 > 0.999$) para reproducir el giro observado en bajas frecuencias, ya que las densidades más bajas implican un radio de influencia mayor donde la emisión de ondas gravitacionales circularizaría la órbita antes de alcanzar las frecuencias observadas.
• Pendiente del perfil: Se prefieren perfiles de densidad más planos (menor $\gamma$) frente a otros más empinados. Los perfiles más empinados conducen a una distribución más amplia de frecuencias de giro en la población de SMBHB, lo que requiere un factor de normalización inconsistente con el modelo fiduciario.
• Comparaciones de referencia: El rango de densidad favorecido incluye naturalmente distribuciones estelares de tipo núcleo (por ejemplo, el núcleo de la Vía Láctea y la distribución estelar de M87), pero desfavorece los picos de materia oscura empinados (por ejemplo, picos adiabáticos en M87), ya que estos últimos resultan en densidades más bajas en el radio de influencia y, por lo tanto, no logran producir el giro observado sin excentricidades extremas.
• Robustez: Las parametrizaciones alternativas que permiten distribuciones en la excentricidad y la densidad no alteran la conclusión principal de que la densidad a la escala del pársec es del orden de $10^6 M_\odot \text{pc}^{-3}$.

Significado y afirmaciones
El artículo afirma que los datos actuales de las PTA, a pesar de su sensibilidad marginal, ya codifican información medible sobre la distribución de materia en escalas de pársec en los centros galácticos. Al modelar conjuntamente la excentricidad y el endurecimiento ambiental, los autores demuestran que el giro espectral en bajas frecuencias en el espectro del SGWB puede utilizarse para inferir el perfil de densidad previo al barrido de los núcleos galácticos.

El estudio destaca que:

1. La dispersión de tres cuerpos es un mecanismo dominante para el endurecimiento orbital en la banda de frecuencia de las PTA, particularmente para binarias de masas comparables.
2. Las características espectrales observadas proporcionan una restricción sobre la densidad de materia combinada (estrellas y materia oscura) cerca de los centros galácticos, sugiriendo densidades mayores que algunas expectativas previas para ciertos entornos.
3. Si la materia oscura es el componente dominante, estos hallazgos ofrecen posibles perspectivas sobre la naturaleza de partícula de la materia oscura y su perfil de densidad en los centros galácticos, aunque los modelos específicos de materia oscura no fría requieren una evaluación adicional.
4. Se espera que las observaciones futuras de las próximas generaciones de PTA (por ejemplo, FAST, DSA-2000, SKA) y las misiones de astrometría proporcionen mediciones espectrales más precisas, resolviendo potencialmente la degeneración entre la excentricidad y la densidad mediante la resolución de binarias individuales o el análisis de correlación a través de los intervalos de frecuencia.

Los autores enfatizan que, si bien establecen un límite superior estricto para la distribución de densidad en el centro galáctico, la determinación precisa del perfil de densidad requiere resolver la degeneración con la excentricidad inicial, potencialmente a través de la resolución de binarias individuales o el análisis de correlación entre los intervalos de frecuencia.