Globular cluster abundance patterns inherited from giant molecular clouds

Resumen Técnico: Patrones de Abundancia de Cúmulos Globulares Heredados de Nubes Moleculares Gigantes

Planteamiento del Problema
Los cúmulos globulares (CG) son sistemas estelares antiguos y densos caracterizados por una anomalía química definitoria: aunque exhiben una dispersión mínima en hierro ([Fe/H]), muestran grandes variaciones de estrella a estrella en elementos ligeros (He, C, N, O, Na) con anticorrelaciones específicas (por ejemplo, [Na/Fe] frente a [O/Fe] y [N/Fe] frente a [O/Fe]). La explicación predominante postula un "autoenriquecimiento intra-cúmulo", donde las estrellas que se forman tempranamente contaminan el gas intra-cúmulo a lo largo de millones de años, dando lugar a una segunda generación de estrellas. Sin embargo, estos modelos enfrentan un severo "problema del presupuesto de masa", que requiere que el cúmulo inicial sea muchas veces más masivo de lo observado hoy para retener suficientes eyectas, o bien invocan poblaciones estelares exóticas y no observadas (por ejemplo, estrellas supermasivas). Además, las recientes observaciones de JWST de galaxias con alto desplazamiento al rojo con un enriquecimiento de nitrógeno extremo han reavivado el interés en los sitios de formación de CG, pero los modelos a menudo dependen de los mismos problemáticos mecanismos de autoenriquecimiento.

Metodología
Los autores utilizan la suite de simulaciones de hidrodinámica de radiación cosmológica THESAN-ZOOM, que modelan la evolución galáctica durante los dos primeros mil millones de años de la historia cósmica. Las simulaciones emplean el código de malla móvil AREPO con una red de enriquecimiento químico "estándar" que rastrea de forma consistente nueve elementos (H, He, C, N, O, Ne, Mg, Si, Fe) utilizando tablas de rendimientos dependientes de la masa y la metalicidad para estrellas AGB, supernovas de colapso de núcleo y supernovas de Tipo Ia.

El estudio identifica Nubes Moleculares Gigantes (GMC, por sus siglas en inglés) utilizando el código CloudPhinder, que localiza estructuras gravitacionalmente ligadas. Para aislar las GMC "similares a CG", los autores aplican criterios de selección química estrictos consistentes con los CG observados en la Vía Láctea:

1. Grandes dispersiones en $\log(\text{N/O})$ ($> 0.5$ dex).
2. $[\text{N/Fe}]$ y $[\text{O/Fe}]$ anticorrelacionados (coeficiente de Pearson $< -0.3$).
3. Dispersión mínima en $[\text{Fe/H}]$ ($\sigma_{[\text{Fe/H}]} < 0.1$ dex).

Para comprender los impulsores físicos de estas nubes, los autores emplean un clasificador de Bosque Aleatorio (Random Forest) para distinguir las GMC similares a CG de la población general de GMC, analizando características como la masa de la nube, el desplazamiento al rojo de formación, las tasas de formación estelar (SFR) durante periodos específicos de tiempo de retroceso y las tasas netas de entrada de gas.

Resultados Clave

1. Identificación de GMC similares a CG: Las simulaciones identifican con éxito una población de GMC masivas que reproducen naturalmente las firmas químicas clave de los CG, incluyendo grandes dispersiones de elementos ligeros y anticorrelaciones N–O con una abundancia de hierro casi constante. Estas nubes son típicamente más masivas que el promedio de las GMC (masa mediana $\sim 0.5$ dex superior).
2. Mecanismo de Formación (El Ciclo de los Bariones): La formación de las GMC similares a CG es impulsada por una fase específica del ciclo de bariones galácticos:
   • Un brote estelar (starburst) expulsa gas rico en oxígeno (proveniente de supernovas de colapso de núcleo) fuera del pozo de potencial poco profundo de una galaxia enana.
   • La galaxia se detiene temporalmente (quench), permitiendo que las estrellas AGB de masa intermedia de poblaciones más antiguas enriquezcan el gas difuso restante con nitrógeno.
   • El gas rico en oxígeno expulsado se enfría en el halo y cae de nuevo hacia la galaxia.
   • La colisión entre esta entrada de flujo de alta velocidad y rica en oxígeno y el gas estacionario y rico en nitrógeno crea una zona de compresión extrema. Esto desencadena el colapso gravitacional rápido de GMC masivas que contienen una mezcla de estos gases químicamente distintos.
3. Herencia vs. Autoenriquecimiento: Los patrones de abundancia química se establecen en el gas antes de la formación estelar. Por consiguiente, las estrellas que se forman dentro de estas nubes heredan las dispersiones de abundancia al nacer ("poblaciones múltiples coevales") en lugar de formarse durante épocas extendidas mediante la auto-contaminación.
4. Vínculo con la Formación de Cúmulos: El análisis de los cúmulos estelares ligados formados a partir de estas nubes muestra que los cúmulos masivos ($M_{*,\text{bound}} \gtrsim 10^5 M_\odot$) con altas densidades superficiales estelares ($\Sigma_* \gtrsim 700 M_\odot \text{pc}^{-2}$) exhiben las anticorrelaciones y dispersiones de abundancia requeridas. Los sistemas de menor densidad muestran comportamientos químicos diversos, consistentes con cúmulos en disolución o sistemas débilmente ligados.
5. Escalas de Tiempo: El escenario predice una escala de tiempo de formación total muy corta para los CG ($\lesssim 3-5$ Myr), restringida por los tiempos dinámicos de las GMC y el requisito de que la formación estelar cese antes de que las supernovas de colapso de núcleo inyecten hierro (lo que crearía una dispersión de hierro).

Significancia y Reivindicaciones
El artículo argumenta que los patrones de abundancia de los cúmulos globulares no requieren necesariamente una formación estelar prolongada intra-cúmulo o contaminadores exóticos. En su lugar, estos patrones pueden ser heredados al nacer de un gas interestelar químicamente estructurado y moldeado por el ciclo de bariones galácticos.

• Resolución del Problema del Presupuesto de Masa: Dado que las estrellas de la "primera" y "segunda" generación se forman simultáneamente a partir de gas pre-enriquecido, no es necesario que el cúmulo inicial sea mucho más masivo que el cúmulo observado final para retener las eyectas.
• Conexión con Observaciones de Alto Desplazamiento al Rojo: El mecanismo explica por qué JWST observa galaxias ricas en nitrógeno con alto desplazamiento al rojo; estos sistemas son probablemente sitios de formación activa de CG impulsados por los mismos ciclos de brote-detención-reinicio (burst-quench-restart).
• Registro Fósil: Los CG sirven como un registro fósil del enriquecimiento químico y los flujos de gas en galaxias de alto desplazamiento al rojo, más que simplemente como reliquias de la dinámica interna de los cúmulos.

Limitaciones y Advertencias
Los autores señalan explícitamente que las simulaciones fueron diseñadas para estudiar la formación de galaxias, no específicamente la de cúmulos estelares. Por consiguiente:

• Las simulaciones no rastrean el sodio (Na), por lo que la anticorrelación $[\text{Na/Fe}]$–$[\text{O/Fe}]$ no está verificada directamente (aunque se espera que siga la tendencia N–O mediante los rendimientos de las estrellas AGB).
• La evolución dinámica a largo plazo de los cúmulos estelares densos está por debajo del límite de resolución; el estudio identifica las condiciones de nacimiento de los cúmulos, pero no puede rastrear su supervivencia a lo largo de escalas de tiempo cosmológicas.
• Se requiere trabajo futuro para madurar este escenario al nivel de detalle logrado por décadas de modelado de autoenriquecimiento.