Pressure-regulated feedback-modulated star formation as a subgrid model for galaxy formation simulations

Resumen Técnico: Modelo de Subrejilla de Formación Estelar Modulada por Retroalimentación y Regulada por Presión

Planteamiento del Problema
Las simulaciones cosmológicas de gran volumen de formación de galaxias (por ejemplo, IllustrisTNG, EAGLE, SIMBA) reproducen con éxito las propiedades estadísticas de las poblaciones de galaxias, pero no pueden resolver las escalas físicas de las nubes moleculares gigantes (GMCs) o el medio interestelar (ISM) multifásico. En consecuencia, estas simulaciones dependen de prescripciones de rejilla subcelular (subgrid) calibradas empíricamente para las tasas de formación estelar (SFR) y la ecuación de estado (EoS) efectiva. Los modelos tradicionales suelen asumir que la SFR específica escala con la densidad del gas como $\dot{m}_*/m_g \propto \rho^{1/2}$ (tiempo de caída libre inverso), con una normalización establecida por observaciones del universo local.

Los autores identifican dos limitaciones primarias en este enfoque:

1. Falta de Conexión Física: El escalamiento empírico carece de un vínculo directo con los procesos físicos que regulan la formación estelar, específicamente el equilibrio entre la gravedad y la retroalimentación estelar.
2. Sensibilidad a la Resolución y al Entorno: Estas prescripciones a menudo fallan en regímenes donde las densidades de gas y presión son significativamente más altas que en el universo local (por ejemplo, galaxias de alto desplazamiento al rojo). Además, los modelos que dependen de la densidad volumétrica local son sensibles a la resolución numérica; si la escala de altura del disco de gas no está resuelta, la densidad volumétrica medida se vuelve inexacta, lo que conduce a predicciones de formación estelar inconsistentes a través de diferentes resoluciones.

Metodología
El artículo introduce un nuevo marco de rejilla subcelular basado en la teoría de la formación estelar Regulada por Presión y Modulada por Retroalimentación (PRFM, por sus siglas en inglés) (Ostriker & Kim 2022). En este marco, la SFR se ajusta para mantener el equilibrio térmico y dinámico vertical en el ISM, donde la presión total del plano medio ($P_{tot}$) equilibra el peso del ISM ($W$) y es proporcional a la SFR por unidad de área ($\Sigma_{SFR}$) mediante un rendimiento de retroalimentación ($\Upsilon_{tot}$).

Los autores implementan dos versiones complementarias de este modelo en el código de malla móvil Arepo, calibradas utilizando simulaciones magnetohidrodinámicas de alta resolución de la suite TIGRESS:

1. PRFM-vol (Volumétrica): Diseñada para simulaciones donde el espesor vertical del disco de gas está resuelto numéricamente (típicamente resolución de masa $\sim 10^5 M_\odot$).

   • Calcula la presión efectiva ($P_{eff}$) en cada celda utilizando una EoS calibrada por TIGRESS ($P_{tot} \propto n_H^{1.8}$).
   • El tiempo dinámico ($t_{dyn}$) se deriva de las densidades volumétricas locales de gas, estrellas y materia oscura.
   • El tiempo de agotamiento ($t_{dep}$) se establece basándose en la presión local y un factor de renormalización ($R_f \approx 2$) para contabilizar la integración vertical de la ley volumétrica.

2. PRFM-int (Integrada): Diseñada para resoluciones más gruesas donde la escala de altura del disco no está resuelta (típicamente resolución de masa $\gtrsim 10^6 M_\odot$).

   • En lugar de depender de la densidad volumétrica local, reconstruye la densidad y la presión del plano medio a partir de consideraciones de equilibrio vertical utilizando densidades superficiales ($\Sigma_g$, $\Sigma_*$) y el peso del ISM ($W$).
   • Resuelve la escala de altura del gas en equilibrio ($H_g$) y deriva $t_{dep}$ directamente de estas cantidades integradas, evitando el uso de densidades volumétricas sensibles a la resolución en la prescripción de la SFR.

Contribuciones Clave

• Calibración con Base Física: El modelo deriva la SFR y la EoS de la condición de equilibrio PRFM ($W = P_{tot} = \Upsilon_{tot}\Sigma_{SFR}$), calibrada contra simulaciones TIGRESS en lugar de ajustes empíricos a galaxias locales.
• Robustez de la Resolución: El desarrollo del enfoque PRFM-int aborda específicamente el fallo de los modelos volumétricos a baja resolución al utilizar densidades superficiales y restricciones de equilibrio, que son menos sensibles a la incapacidad de resolver la estructura vertical del disco.
• Eficiencia Variable: A diferencia de los modelos tradicionales que asumen una eficiencia constante por tiempo de caída libre, el modelo PRFM predice que la eficiencia de formación estelar por tiempo dinámico ($\epsilon_{dyn}$) aumenta con la presión ($\epsilon_{dyn} \propto P_{tot}^{0.43}$).

Resultos
Los autores probaron ambas implementaciones en simulaciones de discos aislados similares a la Vía Láctea a través de resoluciones de masa de $10^5$, $10^6$ y $10^7 M_\odot$, comparándolas con la prescripción estándar de IllustrisTNG.

• Acuerdo de Alta Resolución: A $10^5 M_\odot$, tanto PRFM-vol como PRFM-int concuerdan estrechamente entre sí y con los resultados de TIGRESS para las SFR y los tiempos de agotamiento. Ambos recuperan con éxito la relación Ostriker-Kim calibrada por TIGRESS ($\Sigma_{SFR} \propto P_{tot}^{1.2}$).
• Sensibilidad a la Resolución de los Modelos Volumétricos: El modelo PRFM-vol y el modelo TNG muestran una sensibilidad significativa a la resolución. A medida que la resolución se vuelve más gruesa a $10^6$ y $10^7 M_\odot$, la incapacidad de resolver la escala de altura del disco conduce a mediciones de densidad volumétrica inexactas. Esto causa que el modelo PRFM-vol sobreestime los tiempos de agotamiento y falle en igualar la historia de la SFR de alta resolución.
• Robustez de PRFM-int: La implementación PRFM-int permanece robusta en todas las resoluciones probadas. Rastrea con éxito la historia global de la formación estelar y las distribuciones de densidad superficial de la simulación PRFM-vol de alta resolución ($10^5 M_\odot$), incluso a $10^7 M_\odot$.
• Comparación con TNG: Comparado con el modelo TNG, los modelos PRFM predicen tiempos de agotamiento de gas más cortos, particularmente en regiones de alta presión y alta densidad. Los modelos PRFM exhiben un rango dinámico mucho más amplio en los tiempos de agotamiento y una dependencia más pronunciada de la presión, mientras que el modelo TNG mantiene una eficiencia constante por tiempo de caída libre.

Significancia
El artículo afirma que las prescripciones de rejilla subcelular derivadas de PRFM proporcionan un marco físicamente fundamentado y numéricamente estable para la formación estelar a través del rango dinámico de las simulaciones de formación de galaxias.

• Para Simulaciones de Alta Resolución: El modelo PRFM-vol ofrece un tratamiento más realista de la física del ISM que las leyes empíricas tradicionales, vinculando la SFR directamente con el equilibrio dinámico de la galaxia.
• Para Simulaciones Cosmológicas de Gran Escala: El modelo PRFM-int resuelve el problema crítico de la dependencia de la resolución. Permite que las simulaciones cosmológicas de cajas grandes (que típicamente no pueden resolver las escalas de altura de los discos) utilicen leyes de formación estelar motivadas físicamente que permanezcan consistentes con la física de alta resolución.
• Implicaciones de Alto Desplazamiento al Rojo: Los autores señalan que la eficiencia dependiente de la presión del modelo PRFM sugiere que la formación estelar en galaxias de alto desplazamiento al rojo (donde las presiones son más altas) puede ser significativamente más eficiente de lo predicho por los modelos actuales tipo TNG, impactando potencialmente la interpretación de las galaxias luminosas observadas por JWST.

El trabajo se presenta como un paso hacia el objetivo de la colaboración "Learning the Universe" de desarrollar modelos de rejilla subcelular predictivos, en lugar de empíricamente ajustados, que puedan actualizarse a medida que las simulaciones del ISM de alta resolución (como TIGRESS) se expandan para cubrir espacios de parámetros más amplios (por ejemplo, variando la metalicidad).