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

Technische Samenvatting: IllustrisTNG50 Angular Momentum Kaarten: Het traceren van de Morfo-Kinematische Evolutie van Sterrenstelsels

Probleem en Motivatie
Hoewel moderne kosmologische simulaties succesvol de globale schaleringsrelatie tussen het totale stellaire specifieke impulsmoment ($j_\star$) en de totale stellaire massa ($M_\star$), bekend als de Fall-relatie, reproduceren, is het reproduceren van enkel deze globale relatie onvoldoende om de koppeling tussen de evolutie van sterrenstelsels en de acquisitie van impulsmoment volledig te begrijpen. Observaties van type-late sterrenstelsels hebben onlangs een onverwachte morfo-kinematische diversiteit binnen galactische schijven onthuld, gekenmerkt door onderscheidende ruimtelijke distributies van de stellaire specifieke impulsmoment-oppervlakte-dichtheid (sAMSD). Deze diversiteit suggereert dat impulsmoment wordt herverdeeld naar specifieke substructuren ($j_\star$-substructuren) die niet altijd overeenkomen met traditionele morfologische classificaties. Het primaire doel van dit onderzoek is om deze morfo-kinematische diversiteit binnen kosmologische hydrodynamische simulaties te kwantificeren en te begrijpen, specifiek om te testen of de $j_\star$-substructuren geïdentificeerd in observaties (Pacheco-Arias et al. 2026, hierna PA26) worden gereproduceerd in de IllustrisTNG50-simulatiesuite.

Methodologie
De auteurs analyseerden een steekproef van ongeveer 8.000 stellaire schijven uit de IllustrisTNG50 MW/M31-subsample, reikend over roodverschuivingen $0 \le z \le 3.5$ en stellaire massa's $10^{9.5} \le M_\star \le 10^{11.2} M_\odot$. De steekproef werd samengesteld door SubHalos in de MW/M31-catalogus terug in de tijd te volgen, waarbij alleen die werden geselecteerd die stabiele, rotationeel ondersteunde schijven huisvesten ($V/\sigma \ge 0.5$) met voldoende ruimtelijke resolutie ($R_{max} \ge 7$ kpc) en lage tegenrotatie ($f_{CR} \le 0.20$).

De kernmethodologie omvatte:

1. Constructie van sAMSD-kaarten: Het berekenen van het stellaire specifieke impulsmoment ($j_{\star,k}$) voor individuele deeltjes en het projecteren van de oppervlakte-dichtheid ($\iota_\star = dj_\star/dS$) op een regelmatig rooster (500 pc pixels), gesmooth met een Gaussische kernel die overeenkomt met de stellaire softening-lengte ($\epsilon_\star \approx 300$ pc).
2. Morfo-kinematische metrieken: Vier metrieken werden gedefinieerd om $j_\star$-substructuren te kwantificeren:
   • $C_{2D}$: Een 2D correlatiecoëfficiënt die de gelijkenis meet tussen de sAMSD-kaart en een Freeman-schijfmodel (indicatief voor $j_\star$-ringen).
   • $a_{max}^1$: De maximale genormaliseerde amplitude van de $m=1$ Fourier-modus (het kwantificeren van asymmetrie/$j_\star$-onregelmatigheden).
   • $a_{max}^2$: De maximale genormaliseerde amplitude van de $m=2$ Fourier-modus (het kwantificeren van bi-symmetrie).
   • $R_{max}^2$: Het radiale extensie-fractie waar de $m=2$ modus piekt (het onderscheiden tussen centraal geconcentreerde staven en uitgebreide spiralen).
3. Classificatie: Een Gaussian Mixture Model (GMM) met vier volledig covariante componenten werd toegepast op de vierdimensionale parameterruimte ($C_{2D}, a_{max}^1, a_{max}^2, R_{max}^2$). Deze probabilistische benadering wijsde elk sterrenstelsel een posterieure waarschijnlijkheid toe dat het behoort tot een van de vier geïdentificeerde $j_\star$-typen: $j_\star$-onregelmatig, $j_\star$-spiraal, $j_\star$-ring en $j_\star$-staaf.

Belangrijkste Resultaten

• Morfo-kinematische Diversiteit: De TNG50-simulaties vertonen een morfo-kinematische diversiteit die vergelijkbaar is met observaties. De stellaire schijven slaan impulsmoment op en herverdelen dit voornamelijk in de vier geïdentificeerde $j_\star$-substructuren.
• Evolutionair Verloop: Een canoniek evolutionair pad werd geïdentificeerd op basis van de gemiddelde roodverschuivingen van elke klasse en de Markov-transitie-waarschijnlijkheden:
  • $j_\star$-onregelmatig ($\bar{z} = 0.91$) $\rightarrow$ $j_\star$-spiraal ($\bar{z} = 0.76$) $\rightarrow$ $j_\star$-ring ($\bar{z} = 0.62$) $\rightarrow$ $j_\star$-staaf ($\bar{z} = 0.39$).
  • $j_\star$-onregelmatige en $j_\star$-spiraalvormige systemen fungeren als transitietoestanden, terwijl $j_\star$-ringen en $j_\star$-staven fungeren als dynamische attractoren; eenmaal gevormd, hebben ringen en staven de neiging om te persistenter te zijn.
• Fysieke Drijvers:
  • Gasfractie ($f_{gas}$): Dit is de primaire drijfveer die de klassen scheidt. Een hoge $f_{gas}$ bevoordeelt $j_\star$-onregelmatige en $j_\star$-spiraalvormige systemen, terwijl een lage $f_{gas}$ $j_\star$-ringen en $j_\star$-staven bevoordeelt.
  • Rotatie-ondersteuning ($V/\sigma$): Zodra $f_{gas}$ is vastgesteld, bevoordeelt een hoge $V/\sigma$ $j_\star$-spiralen in gasrijke systemen en $j_\star$-ringen in gasarme systemen.
  • Omgeving: Bij $z=0$ werden geen sterke correlaties gevonden tussen de morfo-kinematische metrieken en de lokale kosmische web-omgeving (dichtheid, afstand tot filamenten/knopen) of de uitlijning van de impulsmomentvector met filamenten.
• $j_\star$-$M_\star$ Diagram: Sterrenstelsels clusteren in afzonderlijke regio's in het $j_\star$-$M_\star$-vlak afhankelijk van hun $j_\star$-type. De best-fit Fall-relatie voor de steekproef levert een helling $\alpha = 0.80 \pm 0.01$, die steiler is dan sommige eerdere TNG50-studies vanwege de kinematische selectie ($V/\sigma \ge 0.5$) die dispersie-gedomineerde systemen uitsluit.

Betekenis en Conclusies
Dit artikel vestigt een canoniek pad voor de herverdeling van impulsmoment binnen seculair evoluerende schijven in TNG50. De studie toont aan dat de $j_\star$-substructuur-classificatie een kader biedt dat meer effectief is in het integreren van dynamische en morfologische informatie dan traditionele morfologie alleen, aangezien $j_\star$-substructuren niet altijd overeenkomen met hun morfologische tegenhangers (bijv. een morfologische staaf impliceert niet noodzakelijkerwijs een $j_\star$-staaf).

De resultaten suggereren dat gasaccretie en stellaire feedback de primaire fysieke processen zijn die de herverdeling van impulsmoment aansturen. De waargenomen prevalentie van $j_\star$-ringen en $j_\star$-staven bij $z=0$ wordt toegeschreven aan de selectiecriteria van de steekproef (die stabiele schijven bevoordeelt) en het gebrek aan efficiënte regulerende mechanismen voor de destructie van staven in de TNG50-simulatie-fysica. De auteurs merken de afwezigheid van $j_\star$-klompen in de simulatie op, waarschijnlijk door resolutiebeperkingen, evenals de overvloed aan $j_\star$-staven, wat wijst op een behoefte aan verfijnde feedbackprocessen in toekomstige simulatie-iteraties. Uiteindelijk koppelt de sAMSD-analyse variaties in stellaire dynamica aan massherverdeling, en biedt het een instrument om uitgebreide galactische evolutiegeschiedenissen te reconstrueren die verder gaan dan standaard kinematische indicatoren.