The host halo masses of AGNs and quasars at $z \sim 3-7$ with TNG-Cluster, FLAMINGO and other cosmological galaxy simulations

Utilizando simulaciones cosmológicas de gran escala, este estudio revela que, si bien la luminosidad de los AGN generalmente aumenta con la masa del halo anfitrión hasta un umbral, la relación es altamente no lineal con una dispersión significativa, lo que indica que los cuásares más luminosos en $z \sim 3-7$ típicamente habitan halos de masa intermedia ($10^{12-12.5}$ M$_{\odot}$) en lugar de los más masivos, un hallazgo que concuerda bien con las estimaciones observacionales.

Lo esencial

La visión general: Encontrando el "punto ideal" para las superestrellas cósmicas

Imagina el universo como una ciudad gigante en crecimiento. En esta ciudad, los "halos de materia oscura" son los vecindarios, y los "Núcleos Galácticos Activos" (AGN) o "Cuásares" son los faros masivos y cegadoramente brillantes en el centro de los edificios más importantes. Estos faros son alimentados por agujeros negros supermasivos devorando gas.

Durante mucho tiempo, los astrónomos han intentado responder a una pregunta sencilla: ¿Qué tan grande debe ser un vecindario (halo) para albergar el faro más brillante (cuásar)?

La mayoría de las observaciones sugieren que estos cuásares súper brillantes viven en vecindarios que tienen aproximadamente el mismo tamaño, sin importar qué tan atrás en el tiempo (o qué tan lejos) mires. Parecen preferir una zona "Goldilocks" (el punto justo): ni demasiado pequeña, ni demasiado grande.

Este artículo utiliza simulaciones computacionales masivas para probar esa idea. Los investigadores construyeron universos digitales (usando modelos como TNG-Cluster y FLAMINGO) para ver si sus cuásares virtuales se comportan como los reales.

El principal descubrimiento: No es una línea recta

Los investigadores descubrieron que la relación entre el tamaño del vecindario y el brillo del faro no es una línea recta. Es más bien como una colina con una cima.

1. El ascenso: A medida que el vecindario se hace más grande (más masivo), el faro generalmente se vuelve más brillante. Esto tiene sentido; los vecindarios más grandes tienen más gas para alimentar al agujero negro.
2. La cima: Esta tendencia solo sube hasta cierto punto (aproximadamente $10^{12}$ veces la masa de nuestro Sol).
3. La caída: Una vez que el vecindario es demasiado masivo, el faro en realidad se vuelve más tenue o se mantiene igual.

La analogía: Piensa en esto como una fiesta.

• En una casa pequeña (halo pequeño), hay pocas personas y la música es suave.
• En una casa de tamaño medio (el punto ideal), tienes la multitud perfecta y la fiesta es ruidosa y enérgica.
• En un estadio masivo (un halo enorme), la fiesta en realidad se vuelve más silenciosa. ¿Por qué? Porque el "anfitrión" (el agujero negro) se vuelve tan codicioso y poderoso que comienza a expulsar a los invitados de la habitación (esto se llama retroalimentación de AGN). Empuja el gas lejos, matando de hambre a sí mismo y evitando que la fiesta se vuelva más ruidosa.

El factor de "caos": Una enorme dispersión

Uno de los hallazgos más sorprendentes es lo desordenados que son los datos. El artículo compara dos formas de observar los datos:

• Escenario A: Si eliges un tamaño de vecindario específico, ¿qué tan brillante es el faro?
  • Resultado: Locamente impredecible. El brillo puede variar por un factor de 1,000 a 10,000 (3 a 4 "décadas" de diferencia). Una casa del mismo tamaño podría tener una pequeña luz nocturna, mientras que otra tiene un reflector cegador.
• Escenario B: Si eliges un brillo específico (por ejemplo, un cuásar súper brillante), ¿qué tan grande es el vecindario?
  • Resultado: Mucho más predecible. Si ves un faro súper brillante, es casi seguro que está en un rango de tamaño de vecindario específico.

La analogía: Imagina intentar adivinar cuánto pesa una persona solo mirando su altura.

• Si eliges una altura específica (por ejemplo, 1.80 metros), las personas pueden pesar desde 60 kg hasta 135 kg. Hay una gran dispersión.
• Pero si eliges un peso específico (por ejemplo, 100 kg), la persona es casi con seguridad de aproximadamente 1.80 metros de altura. El peso predice la altura mucho mejor de lo que la altura predice el peso.

El artículo concluye que el tamaño del vecindario es un mal predictor de qué tan brillante es el agujero negro en este momento. El brillo del agujero negro es caótico y depende de muchas otras cosas (como cuánto gas hay disponible actualmente), no solo del tamaño del vecindario.

El "Punto Ideal" confirmado

A pesar del caos, las simulaciones confirmaron la idea observacional: Los cuásares prefieren un tamaño de vecindario específico.

• Ya sea que el universo tenga 3 mil millones de años o 7 mil millones de años, los cuásares más brillantes tienden a vivir en halos con una masa de aproximadamente $10^{12}$ a $10^{12.5}$ masas solares.
• Rara vez viven en los halos más masivos (las superciudades).
• Rara vez viven en los halos más pequeños (los pueblos diminutos).

Esto sugiere que existe un "punto ideal" para el crecimiento de los agujeros negros. Si el vecindario es demasiado pequeño, el agujero negro no puede obtener suficiente comida. Si es demasiado grande, la propia retroalimentación del agujero negro (expulsar el gas) detiene su crecimiento.

Por qué esto es importante

El artículo argumenta que no podemos simplemente mirar el tamaño del vecindario de una galaxia para adivinar qué tan brillante es su agujero negro. La conexión es demasiado laxa y caótica. Sin embargo, si vemos un cuásar súper brillante, podemos estar bastante seguros del tamaño de su hogar.

Las simulaciones también muestran que este "punto ideal" mantiene aproximadamente el mismo tamaño a lo largo de la historia cósmica, a pesar de que el universo se está expandiendo y cambiando. Esto implica que las reglas que gobiernan cómo se alimentan los agujeros negros y cómo se apagan son consistentes a través de miles de millones de años.

Resumen en una frase

El artículo encuentra que, si bien los agujeros negros supermasivos viven en un tamaño de vecindario cósmico "Goldilocks" específico, el brillo del agujero negro es tan caótico que conocer el tamaño del vecindario te dice muy poco sobre qué tan brillante es el agujero negro en un momento dado.

Resumen técnico

Planteamiento del Problema
Los estudios observacionales y los análisis de agrupamiento han sugerido durante mucho tiempo que los cuásares habitan un estrecho rango de masas de halos de materia oscura ($\sim 10^{12} - 10^{13} \, M_\odot$) a través del tiempo cósmico ($z \lesssim 7$). Sin embargo, estas inferencias son a menudo indirectas, ya que dependen de muestras limitadas por flujo y de supuestos sobre la monotonicidad de la relación entre la luminosidad de los Núcleos Galácticos Activos (AGN) y la masa del halo anfitrión. Existe una brecha crítica en la comprensión de la dispersión intrínseca y la naturaleza precisa de esta relación antes del "Mediodía Cósmico" ($z \gtrsim 3$), particularmente en lo que respecta a si los AGNs más luminosos residen en los halos más masivos y cómo los mecanismos de retroalimentación regulan esta conexión. Las simulaciones previas a menudo estaban limitadas por volumen (incapaces de muestrear halos masivos y raros) o resolución (incapaces de resolver la dinámica del gas a pequeña escala que impulsa las luminosidades de los AGN).

Metodología
Este trabajo utiliza una serie de simulaciones hidrodinámicas cosmológicas para comparar directamente las imágenes observacionales con modelos de formación de galaxias autoconsistentes. El estudio se centra en AGNs con luminosidades bolométricas $L_{\text{AGN}}^{\text{bol}} \geq 10^{42} \, \text{erg s}^{-1}$ en los desplazamientos al rojo $z = 3 - 7$.

El análisis principal se basa en dos grandes series de simulaciones de gran volumen:

1. TNG300 + TNG-Cluster: Parte del proyecto IllustrisTNG, que ofrece una alta resolución (masa objetivo bariónica $\sim 1.1 \times 10^7 \, M_\odot$) a través de un volumen de $300^3 \, \text{cMpc}^3$ y una extensión de zoom dirigida a halos masivos de hasta $10^{15} \, M_\odot$.
2. FLAMINGO: Que abarca volúmenes de $1 - 3 \, \text{Gpc}^3$ por lado (corridas L1_m8 y L2p8_m9), proporcionando un muestreo estadístico superior del extremo de alta masa de la función de masa de halos, aunque con una resolución menor que TNG.

Para la comparación de referencia (benchmarking), los autores incluyen simulaciones de volúmenes más pequeños ($\sim 100 \, \text{cMpc}^3$) de generaciones anteriores: Illustris, EAGLE, TNG100 y Simba.

Las luminosidades bolométricas de los AGN no son salidas directas de la simulación, sino que se calculan en el post-procesamiento a partir de las tasas de acreción de SMBH instantáneas utilizando un modelo de conversión bimodal (Churazov et al. 2005) que distingue entre regímenes de acreción radiativamente eficientes e ineficientes. Se adopta una eficiencia radiativa uniforme de $\epsilon_r = 0.1$ en todos los modelos para aislar las diferencias en la física de acreción subyacente.

Contribuciones Clave y Resultados

• Relaciones No Lineales y No Monotónicas: En todos los modelos, la relación mediana entre la luminosidad bolométrica del AGN y la masa del halo anfitrión es altamente no lineal. Si bien los halos más masivos tienden a albergar AGNs más luminosos en promedio, esta tendencia persiste solo hasta una masa de halo característica de $M_{200,\text{crit}} \sim 10^{12} \, M_\odot$.

  • En TNG300+TNG-Cluster, la relación mediana se aplana e incluso presenta un giro (caída) en $M_{200,\text{crit}} \gtrsim 10^{12} \, M_\odot$ para $z < 5-6$.
  • En FLAMINGO, la relación se aplana en el extremo de alta masa pero no muestra la misma caída pronunciada.
  • Esta extinción (quenching) en altas masas se atribuye a la retroalimentación de los AGN (cinética en TNG, térmica en FLAMINGO) que suprime la acreción de gas hacia los SMBH en halos masivos.

• Gran Dispersión Intrínseca: El artículo cuantifica una variación significativa de objeto a objeto. La dispersión en la luminosidad bolométrica del AGN a una masa de halo fija es mucho mayor (hasta 3 dex para los percentiles 5º–95º) que la dispersión en la masa del halo anfitrión a una luminosidad de AGN fija ($\lesssim 1$ dex). Esto implica un acoplamiento débil entre el crecimiento del halo y la acreción instantánea del SMBH, sugiriendo que la masa del halo por sí sola es un mal predictor de la luminosidad instantánea del AGN.

• Masas de los Halos Anfitriones de Cuásares:

  • TNG300+TNG-Cluster: Los cuásares ($L_{\text{AGN}}^{\text{bol}} \sim 10^{45} - 10^{47} \, \text{erg s}^{-1}$) residen típicamente en halos de $10^{12} - 10^{12.5} \, M_\odot$ en $z=3-6$.
  • FLAMINGO: Los cuásares se extienden a masas medianas de $\sim 10^{12.8} \, M_\odot$ en $z \sim 3-4$.
  • Crucialmente, los AGNs más luminosos no habitan los halos más masivos disponibles en cualquier época dada hasta $z \approx 7$. En su lugar, residen en halos progresivamente más raros en épocas tempranas, pero típicamente evitan los sistemas de mayor masa ($M_{200,\text{crit}} \gtrsim 10^{13} \, M_\odot$) debido a la extinción por retroalimentación.

• Comparación con Observaciones: Las masas de los halos anfitriones simulados para los cuásares coinciden ampliamente con las estimaciones observacionales derivadas de mediciones de agrupamiento y correlación cruzada, las cuales generalmente infieren masas características de $10^{12} - 10^{13} \, M_\odot$. La gran dispersión predicha por las simulaciones ayuda a reconciliar las amplias restricciones observacionales.

Significado y Reivindicaciones
El artículo sostiene que las simulaciones cosmológicas modernas, a pesar de diferir en resolución, volumen y física de sub-malla (semilla, acreción, retroalimentación), convergen en una imagen consistente: la relación entre la luminosidad del AGN y la masa del halo se caracteriza por un "punto ideal" de masa de halo característica ($\sim 10^{12} \, M_\odot$) donde los cuásares son más prevalentes, flanqueado por la extinción en halos más masivos y la supresión en halos menos masivos.

Los autores enfatizan que la gran dispersión intrínseca en la luminosidad del AGN a una masa de halo fija desafía los supuestos de muchos modelos semiempíricos de ocupación de halos, que a menudo asumen una relación más estrecha. Los resultados sugieren que, si bien el cuásar típico reside en un rango de masa específico, los sistemas individuales exhiben una vasta diversidad, y los halos más masivos no son necesariamente los anfitriones de los cuásares más luminosos en un momento dado. Este trabajo proporciona una base cuantitativa para interpretar las observaciones de cuásares de alto desplazamiento al rojo y destaca el papel de la retroalimentación de los AGN en la configuración de la demografía de los agujeros negros más luminosos antes del Mediodía Cósmico.