Testing the inference of kinematics from mock JWST NIRSpec/MSA observations of TNG50 galaxies at $z\sim2-6$

Technische Samenvatting: Het testen van de inferentie van kinematica uit gesimuleerde JWST NIRSpec/MSA-observaties van TNG50-stelsels bij $z \sim 2-6$

Probleemstelling
Het begrijpen van de dynamische evolutie van sterrenstelsels vereist het meten van rotatiesnelheden ($v$) en snelheidsdispersies ($\sigma$) van geïoniseerd gas door de tijd heen. Hoewel hydrodynamische simulaties zoals IllustrisTNG beschikken over 6-dimensionale ruimtelijke en kinematische gegevens, missen zij expliciete stralingsoverdrachtoplossingen voor specifieke observationele opstellingen. Daarentegen vertrouwen observationele studies met instrumenten zoals de JWST NIRSpec Micro-Shutter Array (MSA) op forward modelling om kinematica te infereren uit 2D-spectra. Deze modellen gaan doorgaans uit van axiaal symmetrische, dunne, roterende schijven met Sérsic-intensiteitsprofielen en arctangens-snelheidscurven. Echter, hoog-roodverschuvelingsstelsels ($z > 2$) vertonen vaak complexe morfologieën (sferoïdaal, klonterig) en niet-circulaire bewegingen (inflows/outflows) die deze aannames schenden. Er is een kritiek gebrek aan statistische validatie met betrekking tot de vraag of eenvoudige parametrische modellen de intrinsieke kinematica van typische hoog-roodverschuvelingsstelsels accuraat kunnen herstellen, vooral gezien de beperkte ruimtelijke resolutie en de slit-gebaseerde aard van MSA-observaties vergeleken met Integral Field Units (IFUs).

Methodologie
De auteurs maken gebruik van de hoog-resolutie TNG50 kosmologische simulatie om een end-to-end mock observation pipeline te genereren voor meer dan 4.000 stervormende sterrenstelsels ($M_* = 10^8 - 10^{11.5} M_\odot$) bij roodverschuvelingen $z = 2, 3, 4, 6$.

1. Intrinsieke Metingen: De auteurs definiëren "intrinsieke" kinematische grootheden ($v_{max}$ en $\sigma_e$) binnen de effectieve straal ($r_e$) van het H$\alpha$-emitterende gas. Ze testen expliciet de gevoeligheid van deze intrinsieke waarden voor meetkeuzes, inclusief apertuurgrootte, ruimtelijke binning (0.5 vs. 5 ckpc) en weegschema's (ongewogen vs. Star Formation Rate gewogen).
2. Mock Observaties:
   • Imaging: Mock NIRCam-beelden worden gegenereerd door stervormende gascellen te projecteren, het Star Formation Rate (SFR) om te zetten naar H$\alpha$-luminositeit via de Kennicutt (1998) relatie, en te convolueren met NIRCam PSFs. Ruis wordt toegevoegd om diepe JADES-achtige blootstellingen na te bootsen.
   • Spectroscopie: Mock NIRSpec/MSA 2D-spectra worden gecreëerd door gascellen langs een willekeurige as te projecteren, lijn-van-zicht snelheden toe te wijzen, en te convolueren met de NIRSpec Line Spread Function (LSF) en PSF. Gaussische ruis wordt toegepast om een signaal-ruisverhouding (S/N) te bereiken die vergelijkbaar is met de diepste bestaande JWST-data (S/N $\sim$ 50).
3. Inferentie Pipeline:
   • Morfologie: Mock NIRCam-beelden worden gefit met pysersic om morfologische priors af te leiden (effectieve straal $r_e$, Sérsic-index $n_s$, assenverhouding $q$, positiehoek PA).
   • Kinematica: Mock NIRSpec-spectra worden gefit met behulp van msafit, waarbij wordt uitgegaan van een dun, roterend schijfmodel met een arctangens-snelheidsprofiel en een constante snelheidsdispersie. De uit de imaging afgeleide morfologische parameters worden gebruikt als sterke priors. De fitting maakt gebruik van een gewogen $\chi^2$-minimalisatie, inclusief een niet-parametrische rij-gewijze fluxschaling om rekening te houden met klonterige lichtverdelingen.

Belangrijkste Bijdragen en Resultaten
De studie vergelijkt de best-fit kinematische parameters afgeleid van de mock-spectra met de intrinsieke simulatiegegevens om observationele biases te kwantificeren.

• Herstelnauwkeurigheid: Voor uitgelijnde en opgeloste bronnen bij $z=3$ komen de geïnfereerde rotatiesnelheden ($v$) en snelheidsdispersies ($\sigma$) over het algemeen overeen met de intrinsieke metingen binnen een factor van $\sim 2$ en $\sim 1.5$, respectievelijk. De spreiding in $v$ is groter dan in $\sigma$.
• Systematische Biases:
  • Afhankelijkheid van Morfologie: De herwinning is robuust voor gladde, schijfvormige systemen. Echter, voor intrinsiek verlengde of sferoïdale systemen (grote intrinsieke assenverhoudingen $C/A$), wordt $v$ geneigd overschat en $\sigma$ onderschat met meer dan een factor 2. Dit gebeurt omdat de dunne-schijf aanname het model dwingt om niet-rotationele verbreding te interpreteren als inclinatie-effecten.
  • Roodverschuvelingsevolutie: De spreiding in de nauwkeurigheid van het herstel neemt toe bij hogere roodverschuvelingen ($z > 3$) naarmate sterrenstelsels sterker afwijken van de dunne-schijf aanname.
  • Individueel vs. Populatie: Hoewel de kinematica van individuele sterrenstelsels ($v/\sigma$) met een factor 2–3 verkeerd geclassificeerd kan worden (bijv. een door dispersie ondersteund systeem dat als rotationeel ondersteund verschijnt), worden de populatietrends robuust hersteld. De zwakke evolutie van $\sigma$ en $v/\sigma$ met de roodverschuveling, evenals de afhankelijkheid van $\sigma$ van de globale SFR, blijven behouden in de mock-geïnfereerde data.
• Schaalrelaties: De geïnfereerde $\sigma$-SFR relatie komt goed overeen met de "mass-transport + feedback" theoretische modellen (Krumholz et al. 2018), wat suggereert dat ondanks individuele onzekerheden, kinematische modellering betrouwbaar de fysieke drijfveren van turbulentie (feedback vs. zwaartekracht) in populaties van sterrenstelsels kan vastleggen.
• Foutmodi: Ongeveer 5–10% van de bronnen faalt om het dunne-schijf model te fitten, vaak door onregelmatige morfologieën, sterrenstelsels met hoge stellaire massa en steile centrale snelheidsgradiënten die niet door arctangens-profielen worden gevangen, of bronnen waarbij de intrinsieke $\sigma$ onder de ondergrens van het model ($10$ km s$^{-1}$) ligt.

Significantie en Claims
Het artikel claimt het eerste end-to-end framework te bieden voor het vergelijken van NIRSpec MSA-observaties met kosmologische simulaties om biases bij het meten van geïoniseerde gaskinematica te kwantificeren. De auteurs concluderen dat hoewel eenvoudige dunne-schijf modellen onvoldoende zijn voor het bepalen van de precieze "schijfvormigheid" ($v/\sigma$) van individuele hoog-roodverschuvelingsstelsels vanwege grote onzekerheden, ze wel effectief zijn voor het herstellen van statistische populatietrends.

Het werk benadrukt dat observationele inferenties van $v/\sigma$ voor individuele systemen bij $z \ge 2$ met voorzichtigheid moeten worden behandeld, aangezien ze beïnvloed kunnen worden door de aanname van een dunne-schijf geometrie en de aanwezigheid van niet-circulaire bewegingen. De auteurs stellen dat om de beperkingen van de dynamische staat van vroege sterrenstelsels te verbeteren, toekomstig werk voorbij 2D-slit spectroscopie moet bewegen naar 3D IFU-observaties of slit-stepping strategieën, en toegewijde kinematische modellen moet ontwikkelen die rekening houden met de complexe 3D-morfologieën en klonterigheid die inherent zijn aan hoog-roodverschuvelingsstelsels. De studie benadrukt dat de huidige beperkingen niet louter instrumenteel zijn, maar voortkomen uit de mismatch tussen eenvoudige parametrische modellen en de intrinsieke complexiteit van de fysica van sterrenstelselformatie.