Impact of initial mass function on the chemical evolution of high-redshift galaxies

Título: Impacto de la función de masa inicial en la evolución química de galaxias de alto corrimiento al rojo

1. Planteamiento del Problema

Las observaciones recientes del Telescopio Espacial James Webb (JWST) han revelado una relación invariante entre la masa estelar, la metalicidad y la tasa de formación estelar (SFR) en galaxias hasta $z \sim 8$, la cual parece romperse a corrimientos al rojo más altos ($z > 8$), donde las galaxias son sistemáticamente más pobres en metales. Comprender la física subyacente que impulsa estas correlaciones es crucial.
El problema central radica en determinar si la Función de Masa Inicial (IMF) de las estrellas, específicamente su extremo de alta masa, ha variado en el universo temprano. Las observaciones de un exceso de galaxias luminosas en UV y firmas químicas peculiares (como abundancias mejoradas de N y C) sugieren la presencia de estrellas muy masivas ($\gtrsim 100 M_\odot$) y supernovas de inestabilidad por pares (PISNe), lo que desafía la suposición de una IMF universal (tipo Kroupa) con un límite superior fijo de $100 M_\odot$.

2. Metodología

Los autores utilizan el código semi-analítico de evolución galáctica a-sloth para modelar el ciclo de bariones a lo largo de árboles de fusión derivados de una simulación cosmológica de alta resolución (Ishiyama et al. 2016) en un universo $\Lambda$CDM.

Actualizaciones y Configuración Clave:

• Prescripción de Formación Estelar: Se actualizó el radio característico de la región de formación estelar ($R_*$). En lugar de usar el radio de escala del halo ($R_s$), se adoptó $R_* = 2.5 R_{50}$ (basado en simulaciones zoom-in como NIHAO y VELA), reduciendo el tamaño efectivo de la región de formación estelar en un factor de $\sim 4$. Esto mejora la concordancia con la dispersión observada en la secuencia principal de formación estelar.
• Poblaciones Estelares y Evolución: Se implementaron tracks de evolución estelar actualizados del código PARSEC (v2.0) que cubren un rango completo de metalicidades ($Z \sim 10^{-11} - 0.03$) y masas estelares ($2 - 600 M_\odot$).
  • Se incluyen rendimientos de metales no solo de supernovas (CCSN, PISNe, PPISNe, SNe Ia), sino también de vientos estelares.
  • Se modelan dos poblaciones: Pop III (extremadamente pobres en metales) y Pop II (enriquecidas).
• Parámetros Libres: Se exploró un espacio de parámetros de 165 modelos variando:
  • El límite superior de la masa de las estrellas Pop II ($m_{max}$): de $100$a$600 M_\odot$.
  • La eficiencia de los flujos galácticos (outflows), controlada por la pendiente ($\alpha_{out}$) y la normalización ($M_{out0}$) en relación con la masa del halo.
• Restricciones Observacionales: Los modelos se compararon con tres aspectos observacionales:
  1. Historia de Formación Estelar Cósmica (CSFH): Densidad de tasa de formación estelar (SFRD) y densidad de masa estelar cósmica (CSMD) de $z \sim 4-16$.
  2. Evolución Química: La relación Masa-Metalicidad-SFR (MZSFR) y la Relación Fundamental de Metalicidad (FMR) para $z \sim 4-10$.
  3. Conexión Galaxia-Halo: La relación Masa Estelar-Masa del Halo (SHMR).

3. Contribuciones Clave

• Actualización del Modelo de Retroalimentación: Integración de rendimientos de metales dependientes de la metalicidad y vientos estelares en un código semi-analítico, permitiendo un tratamiento más realista de la evolución química en el universo temprano.
• Análisis de la IMF Variable: Estudio sistemático de cómo un límite superior de masa variable ($m_{max}$) afecta la evolución química, demostrando que una IMF "top-heavy" (más pesada en el extremo superior) es necesaria para explicar las observaciones de JWST.
• Desacople de Parámetros: Identificación de la tensión entre los requisitos para reproducir la CSFH (que favorece flujos de salida específicos) y la MZSFR (que requiere altos rendimientos de metales), resolviéndola mediante una combinación específica de IMF y eficiencia de flujos.

4. Resultados Principales

• Necesidad de Estrellas Muy Masivas: Solo los modelos con un límite superior de masa $m_{max} \gtrsim 200 M_\odot$ logran reproducir simultáneamente la relación MZSFR observada y la historia de formación estelar cósmica a $z \gtrsim 4$.
  • Esto se debe a que las estrellas con $m_{max} \gtrsim 200 M_\odot$ producen Supernovas de Inestabilidad por Pares (PISNe), que aportan rendimientos de metales significativamente mayores (hasta un factor de 3.4 más que modelos con $m_{max}=100 M_\odot$).
• Mejor Ajuste (Best-Fit): El modelo con mayor verosimilitud global tiene los siguientes parámetros:
  • $\alpha_{out} = 0$ (la eficiencia de los flujos de salida es independiente de la masa del halo en el rango explorado).
  • $M_{out0} = 3 \times 10^9 M_\odot$.
  • $m_{max} = 250 M_\odot$.
  • Este modelo reproduce las observaciones dentro de una desviación de $1\sigma$ para la CSFH y la MZSFR.
• Evolución de la Metalicidad: El modelo predice correctamente la ruptura de la relación invariante MZSFR a $z > 8$, donde las galaxias se vuelven más pobres en metales debido a la dilución por flujos de gas prístino y la menor retención de metales en halos pequeños.
• Galaxias Tenues: El modelo predice una población abundante de galaxias tenues y de baja masa ($M_* \sim 1000 M_\odot$) que, aunque no son detectables actualmente, contribuyen con más del 90% de la SFRD total a $z \gtrsim 12$ y son cruciales para la reionización.

5. Significado e Implicaciones

• Evolución Química del Universo Temprano: La confirmación de que estrellas muy masivas ($\gtrsim 200 M_\odot$) y las PISNe jugaron un papel dominante en el enriquecimiento químico de las galaxias a $z \gtrsim 4$ cambia nuestra comprensión de la nucleosíntesis en el universo temprano.
• Fuentes Transitorias Electromagnéticas: Se predice una tasa de eventos de PISNe mucho más alta de lo esperado si solo existieran estrellas Pop III. El pico de la tasa de eventos es de $\sim 10 \text{ yr}^{-1} \text{ deg}^{-2}$ alrededor de $z \sim 5$, dominado por estrellas Pop II masivas. Esto tiene implicaciones para futuras observaciones con telescopios como el Roman Space Telescope y el Vera C. Rubin Observatory.
• Ondas Gravitacionales: La alta tasa de formación estelar en galaxias pobres en metales (predicha por el modelo) implica una mayor eficiencia en la formación de binarias de agujeros negros. Esto sugiere que las galaxias tenues y pobres en metales son los principales progenitores de las fusiones de agujeros negros de alto corrimiento al rojo, detectables por futuros detectores de ondas gravitacionales de 3ª generación.
• Reionización: El modelo subraya el papel crítico de las galaxias tenues y no detectadas en la reionización del universo, apoyando la evidencia observacional de que estas poblaciones ocultas son las principales impulsoras de este proceso.

En conclusión, el estudio demuestra que asumir una IMF universal con un límite de $100 M_\odot$ es insuficiente para explicar los datos de JWST, y que la inclusión de estrellas extremadamente masivas y sus explosiones asociadas es esencial para modelar correctamente la evolución química y dinámica de las galaxias en el universo temprano.