Implications of the nanoHertz Gravitational-Wave Background for Galactic Feedback and Massive Black Hole Growth

Resumen Técnico: Implicaciones del Fondo de Ondas Gravitacionales de nanoHertz para la Retroalimentación Galáctica y el Crecimiento de Agujeros Negros Masivos

Enunciado del Problema
Las mediciones recientes realizadas por Arrays de Temporización de Púlsares (PTA) han proporcionado evidencia sólida de un fondo estocástico de ondas gravitacionales (GWB) de nanoHertz, generado probablemente por la superposición incoherente de binarias de agujeros negros supermasivos (SMBHBs). La amplitud de este fondo es sensible a la función de masa de los agujeros negros supermasivos (SMBH). Sin embargo, existe una tensión: las simulaciones hidrodinámicas cosmológicas (como IllustrisTNG, MillenniumTNG y Simba) subestiman consistentemente la amplitud del GWB inferida por los datos de PTA. Esta discrepancia sugiere que los modelos actuales de crecimiento de SMBH, de siembra y de los mecanismos de retroalimentación que los regulan (retroalimentación de AGN y estelar) pueden estar incompletos. Específicamente, no está claro cómo las variaciones en las prescripciones de retroalimentación galáctica alteran el extremo de alta masa de la Función de Masa de Agujeros Negros (BHMF) y, en consecuencia, la señal del GWB predicha.

Metodología
Los autores emplean un marco semianalítico basado en Casey-Clyde et al. (2022) para vincular las BHMFs predichas por simulaciones con la amplitud del GWB. La metodología procede de la siguiente manera:

1. Suites de Simulación: El estudio analiza cuatro suites principales de simulaciones: IllustrisTNG, MillenniumTNG, Simba y el proyecto CAMELS (Cosmología y Astrofísica con Simulaciones de Aprendizaje Automático).
2. Tamaño de Caja y Convergencia: Para abordar la sensibilidad del GWB al extremo de alta masa de la BHMF, los autores realizan estudios de convergencia utilizando múltiples tamaños de caja (desde 25 Mpc/h hasta 500 Mpc/h). Determinan que las cajas menores a ~35 Mpc/h (para TNG) y ~50 Mpc/h (para Simba) sufren de ruido de Poisson y muestreo insuficiente de halos masivos ($M_{BH} > 10^9 M_\odot$).
3. Variaciones de Retroalimentación: El estudio explora el impacto de la retroalimentación galáctica mediante:
   • Comparar modelos fiduciales entre diferentes códigos de simulación.
   • Analizar variaciones "1P" (de un parámetro) en las suites CAMELS, donde se varían individualmente 23 parámetros astrofísicos que controlan la retroalimentación estelar y de AGN.
   • Examinar variantes específicas de retroalimentación de Simba donde los modos de chorro de AGN se eliminan o alteran (por ejemplo, "SW only", "SW + Rad").
4. Cálculo del GWB: Las BHMFs de las simulaciones se convierten en una población de binarias utilizando una Función de Distribución de Razón de Eddington (ERDF) para identificar cuásares activos. Esta población se escala mediante un factor de normalización ($N_{cal} \approx 4$) para tener en cuenta la relación binaria/cuásar, y se asignan masas secundarias mediante una distribución de razón de masas. La población binaria resultante se integra para calcular el espectro de deformación característica y la amplitud del GWB a una frecuencia de referencia de $f = 1 \text{ yr}^{-1}$.

Contribuciones Clave

• Cuantificación del Impacto de la Retroalimentación: El artículo demuestra que la elección de los modelos de retroalimentación estelar y de AGN altera significativamente el extremo de alta masa de la BHMF. Una retroalimentación de AGN ineficiente o ausente conduce a una sobreabundancia de SMBH masivos, aumentando la amplitud predicha del GWB en factores de 2 (en variaciones fiduciales) a 10 (en variaciones extremas de CAMELS).
• Análisis de Convergencia: Los autores establecen los tamaños de caja mínimos requeridos para predicciones robustas del GWB, identificando que las estadísticas de $M_{BH} > 10^9 M_\odot$ no están convergidas en las cajas estándar de 25 Mpc/h de CAMELS.
• Subestimación Sistemática: El estudio confirma que los modelos fiduciales de IllustrisTNG, MillenniumTNG y Simba subestiman la amplitud del GWB inferida por PTA en un factor de aproximadamente 2. Para coincidir con las observaciones, estos modelos requerirían de 1,5 a 5,5 veces más SMBH en el rango $M_{BH} > 10^8 M_\odot$.

Resultados

• Modelos Fiduciales: Todos los modelos de simulación estándar (fiduciales) subestiman la amplitud del GWB en comparación con los datos de 15 años de NANOGrav. Esto implica que las prescripciones de retroalimentación actuales suprimen el crecimiento de SMBH de manera demasiado eficiente o que la siembra ocurre demasiado tarde.
• Variantes de Retroalimentación: Las simulaciones con retroalimentación de AGN reducida o ausente (específicamente las ejecuciones de Simba sin modos de chorro) producen amplitudes del GWB consistentes con los datos de PTA. Sin embargo, estos modelos no logran reproducir poblaciones galácticas realistas (por ejemplo, no logran apagar la formación estelar de manera efectiva), lo que los hace físicamente poco realistas.
• Sensibilidad a Parámetros: En los conjuntos 1P de CAMELS, parámetros específicos impulsan la amplitud del GWB:
  • Para CAMELS-TNG, el parámetro $\beta$ (que controla el cambio entre estados de AGN de acreción baja y alta) tiene el mayor impacto. Un valor de $\beta=0$ coincide con los datos de PTA pero subestima la función de masa estelar de baja masa.
  • Para CAMELS-Simba, el parámetro de relación de masa de siembra ($A_{BHseed}$) es el más crítico. Disminuir el umbral de siembra permite un crecimiento de SMBH más temprano, coincidiendo con los datos de PTA, pero podría alterar la formación estelar temprana.
• Retroalimentación Estelar: Aunque a menudo se asume que afecta solo a agujeros negros de baja masa, el estudio encuentra que ciertos parámetros de retroalimentación estelar (por ejemplo, la normalización de la velocidad del viento y la fracción térmica) también influyen en la BHMF de alta masa y en la amplitud del GWB al regular el combustible disponible para el crecimiento de los agujeros negros.

Significado y Afirmaciones
El artículo afirma que las mediciones de PTA proporcionan una sonda poderosa e independiente de la física de la retroalimentación galáctica, complementando las restricciones electromagnéticas tradicionales. La discrepancia entre las predicciones de las simulaciones y los datos de PTA sugiere que el crecimiento de SMBH en el universo temprano es más eficiente o ocurre antes de lo capturado por los modelos actuales.

Los autores concluyen que, aunque modificar los parámetros de retroalimentación o las prescripciones de siembra puede alinear artificialmente las simulaciones con los datos de PTA, tales modificaciones a menudo rompen el acuerdo con las estadísticas galácticas observadas (como la función de masa estelar). Por lo tanto, la tensión entre la amplitud del GWB y las predicciones de las simulaciones indica una necesidad fundamental de reevaluar las prescripciones de siembra y crecimiento de SMBH en las simulaciones hidrodinámicas cosmológicas. El estudio enfatiza que los modelos futuros deben reproducir simultáneamente poblaciones galácticas realistas y la amplitud observada del GWB, lo que sugiere que los modelos de retroalimentación actuales podrían estar suprimiendo excesivamente el crecimiento de SMBH. Los autores señalan que su marco se basa en suposiciones simplificadas (por ejemplo, descuidar el endurecimiento ambiental y los detalles del tiempo de fusión), y que se requiere un tratamiento autoconsistente de la evolución binaria para la precisión futura.