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

Resumo Técnico: Mapas de Momento Angular de IllustrisTNG50: Rastreando a Evolução Morfo-Cinemática de Galáxias

Problema e Motivação
Embora as simulações cosmológicas modernas reproduzam com sucesso a relação de escala global entre o momento angular específico estelar total ($j_\star$) e a massa estelar total ($M_\star$), conhecida como relação de Fall, reproduzir apenas esta relação global é insuficiente para compreender totalmente o acoplamento entre a evolução galáctica e a aquisição de momento angular. Observações de galáxias do tipo tardio revelaram recentemente uma diversidade morfo-cinemática inesperada dentro de discos galácticos, caracterizada por distribuições espaciais distintas da densidade superficial de momento angular específico estelar (sAMSD). Essa diversidade sugere que o momento angular é redistribuído em subestruturas específicas ($j_\star$-subestruturas) que nem sempre se alinham com as classificações morfológicas tradicionais. O objetivo principal deste estudo é quantificar e compreender essa diversidade morfo-cinemática em simulações hidrodinâmicas cosmológicas, testando especificamente se as $j_\star$-subestruturas identificadas em observações (Pacheco-Arias et al. 2026, doravante PA26) são reproduzidas na suite de simulação IllustrisTNG50.

Metodologia
Os autores analisaram uma amostra de aproximadamente 8.000 discos estelares do subsample MW/M31 do IllustrisTNG50, abrangendo redshifts $0 \le z \le 3.5$ e massas estelares $10^{9.5} \le M_\star \le 10^{11.2} M_\odot$. A amostra foi construída rastreando SubHalos de volta no tempo a partir do catálogo MW/M31, selecionando apenas aqueles que abrigam discos estáveis e rotacionalmente suportados ($V/\sigma \ge 0.5$) com resolução espacial suficiente ($R_{max} \ge 7$ kpc) e baixa contra-rotação ($f_{CR} \le 0.20$).

A metodologia central envolveu:

1. Construção do Mapa de sAMSD: Calculando o momento angular específico estelar ($j_{\star,k}$) para partículas individuais e projetando a densidade superficial ($\iota_\star = dj_\star/dS$) sobre uma grade regular (pixels de $500$ pc), suavizada com um kernel Gaussiano correspondente ao comprimento de suavização estelar ($\epsilon_\star \approx 300$ pc).
2. Métricas Morfo-Cinemáticas: Quatro métricas foram definidas para quantificar as $j_\star$-subestruturas:
   • $C_{2D}$: Um coeficiente de correlação 2D que mede a semelhança do mapa de sAMSD com um modelo de disco de Freeman (indicativo de $j_\star$-anéis).
   • $a_{max}^1$: A amplitude normalizada máxima do modo de Fourier $m=1$ (quantificando assimetria/$j_\star$-irregularidades).
   • $a_{max}^2$: A amplitude normalizada máxima do modo de Fourier $m=2$ (quantificando bi-simetria).
   • $R_{max}^2$: A fração de extensão radial onde o modo $m=2$ atinge o pico (distinguindo entre barras centralizadas e espirais estendidas).
3. Classificação: Um Modelo de Mistura Gaussiana (GMM) com quatro componentes totalmente covariantes foi aplicado ao espaço de parâmetros quadridimensional ($C_{2D}, a_{max}^1, a_{max}^2, R_{max}^2$). Esta abordagem probabilística atribuiu a cada galáxia uma probabilidade posterior de pertencer a um de quatro tipos de $j_\star$ identificados: $j_\star$-irregular, $j_\star$-espiral, $j_\star$-anel e $j_\star$-barra.

Resultos Principais

• Diversidade Morfo-Cinemática: As simulações TNG50 exibem uma diversidade morfo-cinemática semelhante às observações. Os discos estelares armazenam e redistribuem o momento angular predominantemente nas quatro $j_\star$-subestruturas identificadas.
• Sequência Evolutiva: Uma via evolutiva canônica foi identificada com base nos redshifts médios de cada classe e nas probabilidades de transição de Markov:
  • $j_\star$-irregular ($\bar{z} = 0.91$) $\rightarrow$ $j_\star$-espiral ($\bar{z} = 0.76$) $\rightarrow$ $j_\star$-anel ($\bar{z} = 0.62$) $\rightarrow$ $j_\star$-barra ($\bar{z} = 0.39$).
  • $j_\star$-irregulares e $j_\star$-espirais atuam como estados de transição, enquanto $j_\star$-anéis e $j_\star$-barras atuam como atratores dinâmicos; uma vez formados, anéis e barras tendem a persistir.
• Drivers Físicos:
  • Fração de Gás ($f_{gas}$): Este é o principal driver que separa as classes. O alto $f_{gas}$ favorece $j_\star$-irregulares e $j_\star$-espirais, enquanto o baixo $f_{gas}$ favorece $j_\star$-anéis e $j_\star$-barras.
  • Suporte de Rotação ($V/\sigma$): Uma vez fixado o $f_{gas}$, o alto $V/\sigma$ favorece $j_\star$-espirais em sistemas ricos em gás e $j_\star$-anéis em sistemas pobres em gás.
  • Ambiente: Em $z=0$, não foram encontradas correlações fortes entre as métricas morfo-cinemáticas e o ambiente da teia cósmica local (densidade, distância a filamentos/nós) ou o alinhamento do vetor de momento angular com os filamentos.
• Diagrama $j_\star$-$M_\star$: As galáxias agrupam-se em regiões distintas no plano $j_\star$-$M_\star$ dependendo do seu tipo de $j_\star$. A relação de Fall de melhor ajuste para a amostra fornece uma inclinação $\alpha = 0.80 \pm 0.01$, mais íngreme do que alguns estudos anteriores do TNG50 devido à seleção cinemática ($V/\sigma \ge 0.5$) que exclui sistemas dominados por dispersão.

Significado e Conclusões
O artigo estabelece uma via canônica para a redistribuição do momento angular dentro de discos que evoluem secularmente no TNG50. O estudo demonstra que a classificação de $j_\star$-subestrutura fornece um framework que integra informação dinâmica e morfológica de forma mais eficaz do que a morfologia tradicional, visto que as $j_\star$-subestruturas nem sempre coincidem com seus contrapartes morfológicos (por exemplo, uma barra morfológica não implica necessariamente uma $j_\star$-barra).

Os resultados sugerem que a acreção de gás e o feedback estelar são os principais processos físicos que impulsionam a redistribuição do momento angular. A prevalência observada de $j_\star$-anéis e $j_\star$-barras em $z=0$ é atribuída aos critérios de seleção da amostra (favorecendo discos estáveis) e à falta de mecanismos reguladores eficientes para a destruição de barras na física da simulação TNG50. Os autores observam a ausência de $j_\star$-clumps na simulação, provavelmente devido a limites de resolução, e o excesso de $j_\star$-barras, sugerindo a necessidade de processos de feedback refinados em futuras iterações de simulação. Em última análise, a análise de sAMSD vincula variações na dinâmica estelar à redistribuição de massa, oferecendo uma ferramenta para reconstruir histórias evolutivas galácticas abrangentes que se estendem além dos indicadores cinemáticos padrão.