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

Resumen Técnico: Pruebas de la Inferencia de la Cinemática a partir de Observaciones Simuladas de JWST NIRSpec/MSA de Galaxias TNG50 a $z \sim 2-6$

Planteamiento del Problema
Comprender la evolución dinámica de las galaxias requiere medir las velocidades de rotación ($v$) y las dispersiones de velocidad ($\sigma$) del gas ionizado a través del tiempo cósmico. Si bien las simulaciones hidrodinámicas como IllustrisTNG proporcionan datos espaciales y cinemáticos de 6 dimensiones, carecen de soluciones de transferencia de radiación explícitas para configuraciones observacionales específicas. Por el contrario, los estudios observacionales que utilizan instrumentos como el NIRSpec Micro-Shver Array (MSA) de JWST dependen del modelado hacia adelante para inferir la cinemática a partir de espectros 2D. Estos modelos suelen asumir que las galaxias son discos rotatorios, finos y aximétricos, con perfiles de intensidad de Sérsic y curvas de velocidad de arcotangente. Sin embargo, las galaxias de alto desplazamiento al rojo ($z > 2$) a menudo exhiben morfologías complejas (esferoidales, grumosas) y movimientos no circulares (entradas/salidas de flujo) que pueden violar estas suposiciones. Existe una falta crítica de validación estadística sobre si los modelos paramétricos simples pueden recuperar con precisión la cinemática intrínseca para galaxias típicas de alto desplazamiento al rojo, particularmente dada la resolución espacial limitada y la naturaleza basada en rendijas de las observaciones de MSA en comparación con las Unidades de Campo Integral (IFU).

Metodología
Los autores utilizan la simulación cosmológica de alta resolución TNG50 para generar un flujo de trabajo de observación simulada de extremo a extremo para más de 4,000 galaxias con formación estelar ($M_* = 10^8 - 10^{11.5} M_\odot$) en desplazamientos al rojo $z = 2, 3, 4, 6$.

1. Mediciones Intrínsecas: Los autores definen cantidades cinemáticas "intrínsecas" ($v_{max}$ y $\sigma_e$) dentro del radio efectivo ($r_e$) del gas que emite H$\alpha$. Prueban explícitamente la sensibilidad de estos valores intrínsecos a las elecciones de medición, incluyendo el tamaño de la apertura, el agrupamiento espacial (0.5 vs. 5 ckpc) y los esquemas de ponderación (no ponderado vs. ponderado por Tasa de Formación Estelar o SFR).
2. Observaciones Simuladas:
   • Imagen: Se generan imágenes simuladas de NIRCam proyectando celdas de gas con formación estelar, convirtiendo las Tasas de Formación Estelar (SFR) a luminosidad de H$\alpha$ mediante la relación de Kennicutt (1998), y convolucionando con los PSF de NIRCam. Se añade ruido para imitar exposiciones profundas tipo JADES.
   • Espectroscopía: Se crean espectros 2D simulados de NIRSpec/MSA proyectando las celdas de gas a lo largo de un eje aleatorio, asignando velocidades en la línea de visión y convolucionando con la Función de Propagación de la Línea (LSF) y el PSF de NIRSpec. Se aplica ruido gaussiano para lograr una relación señal-ruido (S/N) comparable a los datos más profundos existentes de JWST (S/N $\sim$ 50).
3. Flujo de Inferencia:
   • Morfología: Las imágenes simuladas de NIRCam se ajustan con pysersic para derivar los priors morfológicos (radio efectivo $r_e$, índice de Sérsic $n_s$, relación de ejes $q$, ángulo de posición PA).
   • Cinemática: Los espectros simulados de NIRSpec se ajustan utilizando msafit, asumiendo un modelo de disco delgado y rotatorio con un perfil de velocidad de arcotangente y una dispersión de velocidad constante. Los parámetros morfológicos derivados de la imagen se utilizan como priors fuertes. El ajuste emplea una minimización de $\chi^2$ ponderada, incluyendo un escalado de flujo no paramétrico por fila para dar cuenta de las distribuciones de luz grumosas.

Contribuciones Clave y Resultados
El estudio compara los parámetros cinemáticos de mejor ajuste derivados de los espectros simulados contra los valores intrínsecos de la simulación para cuantificar los sesgos observacionales.

• Precisión de la Recuperación: Para fuentes alineadas y resueltas en $z=3$, las velocidades de rotación ($v$) y las dispersiones de velocidad ($\sigma$) inferidas generalmente concuerdan con las mediciones intrínsecas dentro de un factor de $\sim 2$ y $\sim 1.5$, respectivamente. La dispersión en $v$ es mayor que en $\sigma$.
• Sesgos Sistemáticos:
  • Dependencia de la Morfología: La recuperación es robusta para sistemas suaves y de tipo disco. Sin embargo, para sistemas intrínsecamente elongados o esferoidales (grandes relaciones de ejes intrínsecos $C/A$), $v$ tiende a sobreestimarse y $\sigma$ a subestimarse por más de un factor de 2. Esto ocurre porque la suposición de disco delgado obliga al modelo a interpretar el ensanchamiento no rotacional como efectos de inclinación.
  • Evolución con el Desplazamiento al Rojo: La dispersión en la precisión de la recuperación aumenta en desplazamientos al rojo más altos ($z > 3$) a medida que las galaxias se desvían más fuertemente de la suposición de disco delgado.
  • Individual vs. Población: Aunque la cinemática de galaxias individuales ($v/\sigma$) puede ser mal clasificada por un factor de 2–3 (por ejemplo, un sistema soportado por dispersión que aparece como soportado por rotación), las tendencias de la población se recuperan de manera robusta. La débil evolución con el desplazamiento al rojo de $\sigma$ y $v/\sigma$, así como la dependencia de $\sigma$ con la SFR global, se preservan en los datos inferidos de la simulación.
• Relaciones de Escala: La relación inferida $\sigma$-SFR se alinea bien con los modelos teóricos de "transporte de masa + retroalimentación" (Krumholz et al. 2018), lo que sugiere que, a pesar de las incertidumbres individuales, el modelado cinemático puede restringir de manera fiable los impulsores físicos de la turbulencia (retroalimentación vs. gravedad) en poblaciones de galaxias.
• Modos de Fallo: Aproximadamente el 5–10% de las fuentes fallan en el ajuste al modelo de disco delgado, a menudo debido a morfologías irregulares, sistemas de alta masa estelar con gradientes de velocidad centrales pronunciados que no son capturados por los perfiles de arcotangente, o fuentes donde la $\sigma$ intrínseca está por debajo del límite inferior del modelo ($10$ km s$^{-1}$).

Significancia y Reivindicaciones
El artículo afirma proporcionar el primer marco de extremo a extremo para comparar las observaciones de NIRSpec MSA con simulaciones cosmológicas para cuantificar los sesgos en la medición de la cinemática del gas ionizado. Los autores concluyen que, si bien los modelos de disco delgado simples son insuficientes para determinar la "forma de disco" ($v/\sigma$) precisa de galaxias individuales de alto desplazamiento al rojo debido a las grandes incertidumbres, son efectivos para recuperar tendencias estadísticas de la población.

El trabajo destaca que las inferencias observacionales de $v/\sigma$ para sistemas individuales en $z \ge 2$ deben tratarse con cautela, ya que pueden estar sesgadas por la suposición de geometría de disco delgado y la presencia de movimientos no circulares. Los autores argumentan que para mejorar las restricciones sobre el estado dinámico de las galaxias tempranas, el trabajo futuro debe ir más allá de la espectroscopía de rendija 2D hacia observaciones de IFU 3D o estrategias de paso de rendija (slit-stepping), y desarrollar modelos cinemáticos dedicados que den cuenta de las morfologías complejas en 3D y la grumosidad inherente a las galaxias de alto desplazamiento al rojo. El estudio enfatiza que las limitaciones actuales no son meramente instrumentales, sino que surgen de la discordancia entre los modelos paramétricos simples y la complejidad intrínseca de la física de la formación de galaxias.