IllustrisTNG50 angular momentum maps: tracing the morpho-kinematic evolution of galaxies

Resumen Técnico: Mapas de Momento Angular de IllustrisTNG50: Rastreo de la Evolución Morfo-Cinemática de las Galaxias

Problema y Motivación
Aunque las simulaciones cosmológicas modernas reproducen con éxito la relación de escala global entre el momento angular específico estelar total ($j_\star$) y la masa estelar total ($M_\star$), conocida como la relación de Fall, reproducir únicamente esta relación global es insuficiente para comprender plenamente el acoplamiento entre la evolución galáctica y la adquisición del momento angular. Observaciones de galaxias de tipo tardío han revelado recientemente una diversidad morfo-cinemática inesperada dentro de los discos galácticos, caracterizada por distribuciones espaciales distintas de la densidad superficial del momento angular específico estelar (sAMSD). Esta diversidad sugiere que el momento angular se redistribuye en subestructuras específicas ($j_\star$-subestructuras) que no siempre se alinean con las clasificaciones morfológicas tradicionales. El objetivo principal de este estudio es cuantificar y comprender esta diversidad morfo-cinemática dentro de simulaciones hidrodinámicas cosmológicas, probando específicamente si las $j_\star$-subestructuras identificadas en observaciones (Pacheco-Arias et al. 2026, en adelante PA26) se reproducen en la suite de simulaciones IllustrisTNG50.

Metodología
Los autores analizaron una muestra de aproximadamente 8,000 discos estelares del subsample MW/M31 de IllustrisTNG50, abarcando desplazamientos al rojo $0 \le z \le 3.5$ y masas estelares $10^{9.5} \le M_\star \le 10^{11.2} M_\odot$. La muestra se construyó rastreando SubHalos hacia atrás en el tiempo desde el catálogo MW/M31, seleccionando solo aquellos que albergan discos estables y con soporte rotacional ($V/\sigma \ge 0.5$) con resolución espacial suficiente ($R_{max} \ge 7$ kpc) y baja contra-rotación ($f_{CR} \le 0.20$).

La metodología central consistió en:

1. Construcción de Mapas de sAMSD: Cálculo del momento angular específico estelar ($j_{\star,k}$) para partículas individuales y proyección de la densidad superficial ($\iota_\star = dj_\star/dS$) sobre una rejilla regular (píxeles de $500$ pc), suavizada con un núcleo Gaussiano que coincide con la longitud de suavizado estelar ($\epsilon_\star \approx 300$ pc).
2. Métricas Morfo-Cinemáticas: Se definieron cuatro métricas para cuantificar las $j_\star$-subestructuras:
   • $C_{2D}$: Un coeficiente de correlación 2D que mide la semejanza del mapa de sAMSD con un modelo de disco de Freeman (indicativo de $j_\star$-anillos).
   • $a_{max}^1$: La amplitud normalizada máxima del modo de Fourier $m=1$ (cuantifica asimetría/irregularidades de $j_\star$).
   • $a_{max}^2$: La amplitud normalizada máxima del modo de Fourier $m=2$ (cuantifica la bi-simetría).
   • $R_{max}^2$: La fracción de extensión radial donde el modo $m=2$ alcanza su máximo (distingue entre barras centralmente concentradas y espirales extendidas).
3. Clasificación: Se aplicó un Modelo de Mezcla Gaussiana (GMM) con cuatro componentes totalmente covariantes al espacio de parámetros tetradimensional ($C_{2D}, a_{max}^1, a_{max}^2, R_{max}^2$). Este enfoque probabilístico asignó a cada galaxia una probabilidad posterior de pertenecer a uno de los cuatro tipos de $j_\star$ identificados: $j_\star$-irregular, $j_\star$-espiral, $j_\star$-anillo y $j_\star$-barra.

Resultados Clave

• Diversidad Morfo-Cinemática: Las simulaciones TNG50 exhiben una diversidad morfo-cinemática similar a las observaciones. Los discos estelares almacenan y redistribuyen el momento angular predominantemente en las cuatro $j_\star$-subestructuras identificadas.
• Secuencia Evolutiva: Se identificó una vía evolutiva canónica basada en los desplazamientos al rojo medios de cada clase y las probabilidades de transición de Markov:
  • $j_\star$-irregular ($\bar{z} = 0.91$) $\rightarrow$ $j_\star$-espiral ($\bar{z} = 0.76$) $\rightarrow$ $j_\star$-anillo ($\bar{z} = 0.62$) $\rightarrow$ $j_\star$-barra ($\bar{z} = 0.39$).
  • Las $j_\star$-irregulares y las $j_\star$-espirales actúan como estados transicionales, mientras que las $j_\star$-anillos y las $j_\star$-barras actúan como atractores dinámicos; una vez formados, los anillos y las barras tienden a persistir.
• Impulsores Físicos:
  • Fracción de Gas ($f_{gas}$): Es el principal impulsor que separa las clases. Una $f_{gas}$ alta favorece a las $j_\star$-irregulares y $j_\star$-espirales, mientras que una $f_{gas}$ baja favorece a las $j_\star$-anillos y $j_\star$-barras.
  • Soporte de Rotación ($V/\sigma$): Una vez fijada la $f_{gas}$, un $V/\sigma$ alto favorece a las $j_\star$-espirales en sistemas ricos en gas y a las $j_\star$-anillos en sistemas pobres en gas.
  • Entorno: En $z=0$, no se encontraron correlaciones fuertes entre las métricas morfo-cinemáticas y el entorno de la red cósmica local (densidad, distancia a filamentos/nodos) o la alineación del vector de momento angular con los filamentos.
• Diagrama $j_\star$-$M_\star$: Las galaxias se agrupan en regiones distintas en el plano $j_\star$-$M_\star$ dependiendo de su tipo de $j_\star$. La relación de Fall ajustada para la muestra arroja una pendiente $\alpha = 0.80 \pm 0.01$, más pronunciada que en algunos estudios previos de TNG50 debido a la selección cinemática ($V/\sigma \ge 0.5$) que excluye sistemas dominados por dispersión.

Significancia y Conclusiones
El artículo establece una vía canónica para la redistribución del momento angular dentro de los discos que evolucionan secularmente en TNG50. El estudio demuestra que la clasificación de $j_\star$-subestructura proporciona un marco que integra la información dinámica y morfológica de manera más efectiva que la morfología tradicional, ya que las $j_\star$-subestructuras no siempre coinciden con sus contrapartes morfológicas (por ejemplo, una barra morfológica no implica necesariamente una $j_\star$-barra).

Los resultados sugieren que la acreción de gas y la retroalimentación estelar son los procesos físicos primarios que impulsan la redistribución del momento angular. La prevalencia observada de $j_\star$-anillos y $j_\star$-barras en $z=0$ se atribuye a los criterios de selección de la muestra (que favorecen discos estables) y a la falta de mecanismos de regulación eficientes para la destrucción de barras en la física de las simulaciones TNG50. Los autores señalan la ausencia de $j_\star$-clumps (grumos) en la simulación, probablemente debido a límites de resolución, y la sobreabundancia de $j_\star$-barras, lo que sugiere la necesidad de procesos de retroalimentación más refinados en futuras iteraciones de las simulaciones. En última instancia, el análisis de sAMSD vincula las variaciones en la dinámica estelar con la redistribución de masa, ofreciendo una herramienta para reconstruir historias evolutivas galácticas integrales que se extienden más allá de los indicadores cinemáticos estándar.