Resolving Star Cluster Formation in Galaxy Simulations with Cosmic Ray Feedback

Resumen Técnico: Resolución de la formación de cúmulos estelares en simulaciones galácticas con retroalimentación de rayos cósmicos

Planteamiento del problema
Los cúmulos estelares son los sitios primarios de la formación de estrellas masivas, lo que impulera los procesos de retroalimentación que dan forma al medio interestelar (ISM) y regulan la formación estelar futura. Una comprensión exhaustiva de la retroalimentación galáctica requiere resolver la formación de cúmulos a escala de pársecs y el ISM multifase. Se sabe que los rayos cósmicos (CR) contribuyen a la presión y el calentamiento del ISM, pero su impacto específico en las tasas de formación de cúmulos estelares (SFRs), las funciones de masa de los cúmulos y la vinculación gravitacional permanece incierto. Estudios previos han explorado la retroalimentación de CR en el contexto de flujos galácticos o han utilizado modelos de subrejilla para los cúmulos, pero pocos han resuelto simultáneamente partículas estelares individuales, CR acoplados dinámicamente y un ISM multifase regulado por retroalimentación. Este trabajo aborda la brecha en la comprensión de cómo los modelos de transporte de CR, específicamente aquellos basados en la teoría de autoconfinamiento, influyen en la formación y las propiedades de los cúmulos estelares.

Metodología
Los autores analizan simulaciones de alta resolución tipo "tallbox" realizadas con el código de malla móvil Areto. Las simulaciones emplean el modelo de termoquímica fuera de equilibrio y retroalimentación Crisp, que rastrea los procesos de enfriamiento/calentamiento para el hidrógeno y los metales. El estudio utiliza tres casos de física distintos para aislar los efectos del transporte de CR:

1. MHD: Magnetohidrodinámica con CR tratados como un término de calentamiento/ionización uniforme.
2. CR-NL: Magnetohidrodinámica de dos momentos (CRMHD) con amortiguamiento de Landau no lineal (NLLD) como único mecanismo de amortiguamiento, lo que resulta en un acoplamiento CR-gas uniforme.
3. CR-NL-IN (Fiducial): CRMHD con tanto NLLD como amortiguamiento de ion-neutro (IND). En este caso, los CR se desacoplan del gas neutro frío, reduciendo su impacto local en la formación estelar.

Las simulaciones resuelven el gas con una resolución de masa de $10 M_\odot$ y una escala espacial de $\sim 0.5$ pc. La formación estelar se implementa mediante una prescripción de tipo Schmidt donde las celdas de gas que exceden un umbral de densidad ($n_H = 10^3 \text{ cm}^{-3}$) se convierten en partículas estelares colisionales con una eficiencia local del 100% por tiempo de caída libre. La retroalimentación estelar se modela exclusivamente a través de supernovas (SNe); la retroalimentación radiativa y los vientos estelares se omiten para aislar el impacto de los CR y las SNe. Las SNe se inyectan estocásticamente para estrellas con edades de 3 a 38 Myr, con un 5% adicional de energía de SN inyectada como CR en los casos CRMHD.

Los cúmulos estelares se identifican utilizando un algoritmo modificado de Amigos de Amigos (FoF) que incorpora una coordenada temporal, asegurando que los grupos identificados consistan en estrellas coetáneas dentro de una longitud de enlace espacial de 2 pc. El análisis se centra en cúmulos "jóvenes" (edades $\sim 10$ Myr) para caracterizar sus propiedades iniciales antes de una disrupción tidal significativa.

Resultos clave

1. Funciones de masa de cúmulos: Los tres casos de simulación producen funciones de masa de cúmulos que siguen una ley de potencia ($dN/dM \propto M^{\alpha_{PL}}$) con pendientes entre $-2.0$y$-2.4$ para masas $> 10^3 M_\odot$. Esto coincide con las tendencias observacionales donde la pendiente se correlaciona positivamente con la densidad superficial de la tasa de formación estelar ($\Sigma_{SFR}$). La inclusión de CR desplaza los puntos de datos a lo largo de estas relaciones observadas pero no produce resultados no físicos.
2. Eficiencia de formación de cúmulos (CFE): La CFE ($\Gamma$) muestra una correlación positiva con $\Sigma_{SFR}$ cuando se calcula para cúmulos masivos ($M > 5000 M_\odot$), consistente con las predicciones teóricas (p. ej., Kruijssen 2012) y las observaciones. Sin embargo, para umbrales de masa más bajos ($M > 1000 M_\odot$), los casos CRMHD muestran una CFE alta y constante, mientras que el caso MHD muestra valores más bajos. Los autores sugieren que la alta CFE en los casos CRMHD puede ser un artefactio del alto umbral de densidad y la eficiencia utilizada en la prescripción de formación estelar.
3. Agrupamiento de supernovas: Aproximadamente el 50% de las SNe ocurren en gas difuso y precondicionado ($n_H < 10^{-2} \text{ cm}^{-3}$) en los tres casos. Esta fracción es mayormente independiente del modelo de CR, lo que indica que la retroalimentación de CR no altera significativamente el entorno local donde explotan las primeras SNe, a pesar de cambiar la fracción de volumen global del ISM difuso.
4. Propiedades de los cúmulos y parámetros de virialidad:
   • Radios y dispersiones de velocidad: Los cúmulos estelares en el caso MHD exhiben radios de semimasa y dispersiones de velocidad de 1D más grandes en comparación con los casos CRMHD.
   • Vinculación gravitacional: Los cúmulos MHD son, en promedio, no ligados ($\alpha_{vir} > 1$), mientras que los cúmulos CRMHD son, en promedio, ligados ($\alpha_{vir} < 1$).
   • Mecanismo: Los autores atribuyen estas diferencias al presupuesto de energía turbulenta del ISM. El caso MHD tiene la SFR más alta, lo que conduce a una mayor tasa de inyección de energía de SN y un ISM más turbulento. Esta mayor turbulencia ambiental es heredada por las nubes moleculares en formación, resultando en cúmulos estelares con mayores dispersiones de velocidad que están menos ligados gravitacionalmente. La retroalimentación de CR reduce la SFR, disminuyendo así el presupuesto de energía turbulenta y produciendo cúmulos más estrechamente ligados.

Significancia y afirmaciones
El artículo afirma presentar las primeras simulaciones de formación de cúmulos estelares que incluyen simultáneamente CR acoplados dinámicamente, un ISM multifase e partículas estelares individuales representativas de estrellas masivas. La principal significancia radica en demostrar que:

• El modelo de retroalimentación Crisp puede reproducir de manera consistente la formación de cúmulos estelares con pendientes de función de masa y dependencias ambientales observadas sin necesidad de un ajuste fino, incluso en ausencia de retroalimentación radiativa de pre-SN.
• El impacto primario de la retroalimentación de CR en este contexto es la reducción de la SFR global. Esta reducción disminuye el presupuesto de energía turbulenta del ISM, lo que a su vez produce cúmulos estelares que están más ligados gravitacionalmente (menores dispersiones de velocidad) en comparación con las simulaciones sin retroalimentación de CR efectiva.
• Los CR no parecen alterar dinámicamente el colapso de nubes moleculares densas individuales ni alterar significativamente la precondición local de las SNe; su efecto es mediado a través de la regulación global del estado del ISM.

Los autores concluyen que, aunque su configuración idealizada omite la retroalimentación radiativa y los vientos estelares, las diferencias relativas entre los casos respecto al impacto de los CR en la vinculación de los cúmulos y la SFR son robustas. Los resultados proporcionan nuevas restricciones para los modelos de transporte de CR y sugieren que la retroalimentación de CR actúa para mover los resultados de las simulaciones a lo largo de las relaciones de cúmulos estelares observadas, en lugar de crear nuevos regímenes no observados.