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

Sintesi Tecnica: Mappe del Momento Angolare in IllustrisTNG50: Tracciare l'Evoluzione Morfo-Cinematica delle Galassie

Problema e Motivazione
Sebbene le moderne simulazioni cosmologiche riescano a riprodurre la relazione di scala globale tra il momento angolare specifico stellare totale ($j_\star$) e la massa stellare totale ($M_\star$), nota come relazione di Fall, riprodurre solo questa relazione globale è insufficiente per comprendere appieno l'accoppiamento tra l'evoluzione galattica e l'acquisizione del momento angolare. Le osservazioni delle galassie di tipo tardivo hanno recentemente rivelato una diversità morfo-cinematica inaspettata all'interno dei dischi galattici, caratterizzata da distinte distribuzioni spaziali della densità superficiale del momento angolare specifico stellare (sAMSD). Questa diversità suggerisce che il momento angolare venga ridistribuito in specifiche sottostrutture ($j_\star$-sottostrutture) che non sempre si allineano con le tradizionali classificazioni morfologiche. L'obiettivo primario di questo studio è quantificare e comprendere questa diversità morfo-cinematica all'interno di simulazioni idrodinamiche cosmologiche, testando specificamente se le $j_\star$-sottostrutture identificate nelle osservazioni (Pacheco-Arias et al. 2026, hereafter PA26) siano riprodotte nella suite di simulazioni IllustrisTNG50.

Metodologia
Gli autori hanno analizzato un campione di circa 8.000 dischi stellari dal subsampolo MW/M31 di IllustrisTNG50, coprendo redshift $0 \le z \le 3.5$ e masse stellari $10^{9.5} \le M_\star \le 10^{11.2} M_\odot$. Il campione è stato costruito tracciando i SubHalo nel catalogo MW/M31 a ritroso nel tempo, selezionando solo quelli che ospitano dischi stabili e rotazionalmente supportati ($V/\sigma \ge 0.5$) con risoluzione spaziale sufficiente ($R_{max} \ge 7$ kpc) e bassa contro-rotazione ($f_{CR} \le 0.20$).

Il nucleo della metodologia ha previsto:

1. Costruzione delle mappe sAMSD: Calcolo del momento angolare specifico stellare ($j_{\star,k}$) per singole particelle e proiezione della densità superficiale ($\iota_\star = dj_\star/dS$) su una griglia regolare (pixel di $500$ pc), mediata con un kernel gaussiano che corrisponde alla lunghezza di smoothing stellare ($\epsilon_\star \approx 300$ pc).
2. Metriche Morfo-Cinematiche: Sono state definite quattro metriche per quantificare le $j_\star$-sottostrutture:
   • $C_{2D}$: Un coefficiente di correlazione 2D che misura la somiglianza della mappa sAMSD con un modello di disco di Freeman (indicativo di $j_\star$-anelli).
   • $a_{max}^1$: L'ampiezza massima normalizzata del modo di Fourier $m=1$ (quantifica l'asimmetria/$j_\star$-irregolarità).
   • $a_{max}^2$: L'ampiezza massima normalizzata del modo di Fourier $m=2$ (quantifica la bi-simmetria).
   • $R_{max}^2$: L'estensione radiale frazionaria dove il modo $m=2$ raggiunge il picco (distinguendo tra barre centralmente concentrate e spirali estese).
3. Classificazione: Un Modello a Miscela Gaussiana (GMM) con quattro componenti completamente covarianti è stato applicato allo spazio dei parametri quadridimensionale ($C_{2D}, a_{max}^1, a_{max}^2, R_{max}^2$). Questo approccio probabilistico ha assegnato a ciascuna galassia una probabilità a posteriori di appartenenza a uno dei quattro tipi di $j_\star$ identificati: $j_\star$-irregolare, $j_\star$-spirale, $j_\star$-anello e $j_\star$-barra.

Risultati Chiave

• Diversità Morfo-Cinematica: Le simulazioni TNG50 mostrano una diversità morfo-cinematica simile alle osservazioni. I dischi stellari immagazzinano e ridistribuiscono il momento angolare prevalentemente nelle quattro $j_\star$-sottostrutture identificate.
• Sequenza Evolutiva: È stata identificata una via evolutiva canonica basata sui redshift medi di ciascuna classe e sulle probabilità di transizione di Markov:
  • $j_\star$-irregolare ($\bar{z} = 0.91$) $\rightarrow$ $j_\star$-spirale ($\bar{z} = 0.76$) $\rightarrow$ $j_\star$-anello ($\bar{z} = 0.62$) $\rightarrow$ $j_\star$-barra ($\bar{z} = 0.39$).
  • Le $j_\star$-irregolari e le $j_\star$-spirali agiscono come stati transitori, mentre le $j_\star$-anelli e le $j_\star$-barre agiscono come attrattori dinamici; una volta formate, gli anelli e le barre tendono a persistere.
• Driver Fisici:
  • Frazione di Gas ($f_{gas}$): È il driver principale che separa le classi. Un alto $f_{gas}$ favorisce $j_\star$-irregolari e $j_\star$-spirali, mentre un basso $f_{gas}$ favorisce $j_\star$-anelli e $j_\star$-barre.
  • Supporto Rotazionale ($V/\sigma$): Una volta fissato $f_{gas}$, un alto $V/\sigma$ favorisce le $j_\star$-spirali nei sistemi ricchi di gas e le $j_\star$-anelli nei sistemi poveri di gas.
  • Ambiente: A $z=0$, non sono state trovate correlazioni forti tra le metriche morfo-cinematiche e l'ambiente locale della ragnatela cosmica (densità, distanza da filamenti/nodi) o l'allineamento del vettore del momento angolare con i filamenti.
• Diagramma $j_\star$-$M_\star$: Le galassie si raggruppano in regioni distinte nel piano $j_\star$-$M_\star$ a seconda del loro tipo di $j_\star$. La relazione di Fall con il miglior fit per il campione fornisce una pendenza $\alpha = 0.80 \pm 0.01$, più ripida rispetto ad alcuni studi precedenti di TNG50 a causa della selezione cinematica ($V/\sigma \ge 0.5$) che esclude i sistemi dominati dalla dispersione.

Significato e Conclusioni
L'articolo stabilisce un percorso canonico per la ridistribuzione del momento angolare all'interno dei dischi in evoluzione secolare in TNG50. Lo studio dimostra che la classificazione delle $j_\star$-sottostrutture fornisce un quadro che integra le informazioni dinamiche e morfologiche in modo più efficace rispetto alla morfologia tradizionale, poiché le $j_\star$-sottostrutture non sempre coincidono con i corrispettivi morfologici (ad esempio, una barra morfologica non implica necessariamente una $j_\star$-barra).

I risultati suggeriscono che l'accrezione di gas e il feedback stellare siano i processi fisici primari che guidano la ridistribuzione del momento angolare. La prevalenza osservata di $j_\star$-anelli e $j_\star$-barre a $z=0$ è attribuita ai criteri di selezione del campione (che favoriscono dischi stabili) e alla mancanza di meccanismi regolatori efficienti per la distruzione delle barre nella fisica delle simulazioni TNG50. Gli autori notano l'assenza di $j_\star$-clump (ammassi) nella simulazione, probabilmente dovuta ai limiti di risoluzione, e l'eccesso di $j_\star$-barre, suggerendo la necessità di processi di feedback più raffinati nelle future iterazioni delle simulazioni. In definitiva, l'analisi sAMSD collega le variazioni nella dinamica stellare alla ridistribuzione della massa, offrendo uno strumento per ricostruire storie evolutive galattiche complete che vanno oltre gli indicatori cinematici standard.