Unveiling the population of massive quenched galaxies at $z\ge2$ in the COLIBRE simulations -- I. Galaxy demographics

Resumen Técnico: Revelando la población de galaxias masivas extinguidas a $z \ge 2$ en las simulaciones COLIBRE – I. Demografía galáctica

Planteamiento del Problema
Observaciones recientes del Telescopio Espacial James Webb (JWST) han revelado una población sustancial de Galaxias Masivas y Extinguidas (MQGs; $M_\star \gtrsim 10^{10} M_\odot$) en altos desplazamientos al rojo ($z \gtrsim 2$). Las densidades numéricas inferidas y los rápidos tiempos de formación de estos sistemas presentaron inicialmente tensiones significativas con el modelo jerárquico $\Lambda$CDM y las simulaciones existentes de formación de galaxias, las cuales tuvieron dificultades para reproducir tales sistemas masivos y pasivos tan temprano en la historia cósmica. Si bien análisis observacionales posteriores han revisado las densidades numéricas a la baja, persisten discrepancias entre los modelos actuales y los últimos datos del JWST con respecto a la abundancia, las funciones de masa estelar y las propiedades físicas de las MQGs.

Metodología
Este estudio utiliza las nuevas simulaciones hidrodinámicas cosmológicas COLIBRE, que introducen innovaciones clave en la física de sub-malla en comparación con simulaciones de gran volumen anteriores. El análisis se centra en la simulación de referencia L200m6 (un volumen de $200 \, \text{cMpc}^3$ con masa media de partícula de gas $m_g \sim 1.84 \times 10^6 M_\odot$), junto con ejecuciones comparativas que varían en volumen, resolución y modelos de retroalimentación de Núcleos Galácticos Activos (AGN) (térmica vs. híbrida térmica+chorro/jet).

Características metodológicas clave:

• Física de sub-malla: Seguimiento explícito de la fase molecular fría y granos de polvo (mediante el modelo HYBRID-CHIMES), enriquecimiento químico consistente y prescripciones actualizadas de retroalimentación de AGN (incluyendo acreción super-Eddington y retroalimentación de chorro/jet).
• Selección de la muestra: Las MQGs se definen como galaxias con $M_\star > 10^{10} M_\odot$ y tasas de formación estelar específicas (sSFR) $< 0.2/t_{\text{age}}$ en $z \ge 2$.
• Modelado hacia adelante (Forward Modeling): Para abordar las incertidumbres observacionales, los autores conllevan las masas estelares y las SFR simuladas con una distribución Gaussiana ($\sigma = 0.3$ dex) para imitar los errores inherentes al ajuste de SED (Distribución de Energía Espectral) y los desplazamientos al rojo fotométricos.
• Análisis de propiedades: El estudio examina densidades numéricas, funciones de masa estelar (SMF), historias de formación estelar (SFHs), tiempos de formación y extinción ($t_{50}$ y $t_q$), fracciones de polvo e hidrógeno molecular ($H_2$) y cinemática ($v_\star/\sigma_\star$).

Contribuciones Clave y Resultados

1. Densidades Numéricas y Dependencia de la Resolución:

   • La densidad numérica de las MQGs en COLIBRE depende tanto de la resolución de masa como de la implementación de la retroalimentación de AGN. Las simulaciones de mayor resolución (por ejemplo, L200m7) producen menos galaxias extinguidas porque una mayor resolución conduce a gases más densos alrededor de los agujeros negros, lo que requiere una retroalimentación más fuerte para extinguir, y porque la retroalimentación de AGN se vuelve menos errática (bursty).
   • La retroalimentación de AGN híbrida (térmica + chorros) produce significativamente menos galaxias extinguidas a $z \gtrsim 2$ en comparación con el modelo térmico de referencia, ya que los agujeros negros crecen más lentamente en el universo temprano en el escenario híbrido.
   • Cuando se aplica el modelado hacia adelante de las incertidumbres observacionales (0.3 dex), las densidades numéricas predichas por la simulación L200m6 de referencia están en amplio acuerdo con las últimas restricciones espectroscópicas y fotométricas del JWST, reduciendo significativamente la tensión vista en estudios previos.

2. Funciones de Masa Estelar (SMF):

   • La SMF de las MQGs en COLIBRE tiene un pico alrededor de $M_\star \approx 10^{10.5} M_\odot$, lo cual es consistente con los datos observacionales de Baker et al. (2025c) y Stevenson et al. (2026).
   • A diferencia de otros modelos (por ejemplo, EAGLE, IllustrisTNG, Simba) que luchan por igualar simultáneamente los extremos de baja y alta masa, COLIBRE reproduce la SMF en todo el rango de masa cuando se contabilizan los errores observacionales. Esto se atribuye al modelado explícito de la fase de gas frío, lo que permite densidades de gas más altas, acreción de agujeros negros super-Eddington y, consecuentmente, una extinción temprana más eficiente.

3. Historias de Formación Estelar (SFHs):

   • Tiempos de formación ($t_{50}$): Las MQGs en COLIBRE se forman rápidamente, con tiempos de formación medianos de $t_{50} \approx 0.5 - 1.5$ Gyr. Esto indica una formación extendida seguida de una extinción rápida.
   • Tiempos de extinción ($t_q$): Las simulaciones predicen una extinción rápida con tiempos medianos de $t_q \lesssim 0.6$ Gyr, a menudo impulsada por episodios intensos de brotes estelares (starbursts).
   • Las SFHs frecuentemente exhiben múltiples episodios de brotes estelares antes de la extinción final, una característica que complica la recuperación observacional pero que es cualitativamente consistente con las inferencias de los datos del JWST.

4. Contenido de Polvo y Gas Molecular:

   • El modelo predice que las MQGs poseen fracciones de polvo y $H_2$ más de 1 dex menores que sus contrapartes masivas en formación de estrellas.
   • Estas bajas fracciones son consistentes con las escasas estimaciones observacionales y los límites superiores disponibles, apoyando la visión de que los procesos de retroalimentación agotan eficientemente los reservorios de gas y suprimen el crecimiento de polvo en estos sistemas.

5. Tamaños y Cinemática:

   • Las MQGs y los sistemas masivos en formación en COLIBRE muestran tamaños y cinemáticas ($v_\star/\sigma_\star$) ampliamente similares. Esto sugiere que transformaciones significativas de tamaño o morfología en estas galaxias ocurren después de la extinción, en lugar de ser un prerrequisito para el estado extinto en $z \ge 2$.

Significado
Este artículo establece que las simulaciones COLIBRE, aprovechando la física de sub-malla de vanguardia (particularmente en polvo, gas frío y retroalimentación de AGN), pueden reproducir robustamente la demografía de las galaxias masivas extinguidas en alto desplazamiento al rojo. El estudio demuestra que las aparentes tensiones entre las simulaciones y las observaciones del JWST se alivian en gran medida cuando se modela hacia adelante de forma realista las incertidumbres observacionales. Al proporcionar un marco consistente para la abundancia de las MQGs, el ensamblaje de masa y los tiempos de extinción, este trabajo ofrece una vía viable para comprender los orígenes de estos tempranos sistemas pasivos y sienta las bases para el Papel II, que investigará los mecanismos físicos específicos que impulsan su rápida formación y extinción.