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

Resumo Técnico: Testando a Inferência de Cinemática de Observações Mock de JWST NIRSpec/MSA de Galáxias TNG50 em $z \sim 2-6$

Enunciado do Problema
Compreender a evolução dinâmica das galáxias requer a medição das velocidades de rotação ($v$) e das dispersões de velocidade ($\sigma$) do gás ionizado ao longo do tempo cósmico. Embora simulações hidrodinâmicas como o IllustrisTNG forneçam dados espaciais e cinemáticos de 6 dimensões, elas carecem de soluções explícitas de transferência de radiação para configurações observacionais específicas. Por outro vez, estudos observacionais utilizando instrumentos como o JWST NIRSpec Micro-Shutter Array (MSA) dependem de modelagem direta para inferir a cinemática a partir de espectros 2D. Esses modelos tipicamente assumem que as galáxias são discos rotativos finos e axissimétricos com perfis de intensidade Sérsic e curvas de velocidade de arcotangente. No entanto, galáxias de alto redshift ($z > 2$) frequentemente exibem morfologias complexas (esferoidais, clumpy) e movimentos não circulares (inflows/outflows) que podem violar essas suposições. Existe uma falta crítica de validação estatística sobre se modelos paramétricos simples podem recuperar com precisão a cinemática intrínseca para galáxias típicas de alto redshift, particularmente dada a resolução espacial limitada e a natureza baseada em fendas (slit-based) das observações de MSA em comparação com Unidades de Campo Integral (IFUs).

Metodologia
Os autores utilizam a simulação cosmológica de alta resolução TNG50 para gerar um pipeline de observação mock de ponta a ponta para mais de 4.000 galáxias formadoras de estrelas ($M_* = 10^8 - 10^{11.5} M_\odot$) em redshifts $z = 2, 3, 4, 6$.

1. Medições Intrínsecas: Os autores definem quantidades cinemáticas "intrínsecas" ($v_{max}$ e $\sigma_e$) dentro do raio efetivo ($r_e$) do gás emissor de H$\alpha$. Eles testam explicitamente a sensibilidade desses valores intrínsecos às escolhas de medição, incluindo tamanho de abertura, agrupamento espacial (0.5 vs. 5 ckpc) e esquemas de ponderação (não ponderado vs. ponderado pela Taxa de Formação Estelar - SFR).
2. Observações Mock:
   • Imagens: Imagens mock de NIRCam são geradas projetando células de gás formadoras de estrelas, convertendo Taxas de Formação Estelar (SFR) em luminosidade de H$\alpha$ via relação de Kennicutt (1998), e convolvendo com PSFs de NIRCam. Ruído é adicionado para mimetizar exposições profundas do tipo JADES.
   • Espectroscopia: Espectros 2D de NIRSpec/MSA mock são criados projetando células de gás ao longo de um eixo aleatório, atribuindo velocidades na linha de visada e convolvendo com a Função de Espalhamento de Linha (LSF) e o PSF do NIRSpec. Ruído gaussiano é aplicado para atingir uma relação sinal-ruído (S/N) comparável aos dados mais profundos existentes do JWST (S/N $\sim$ 50).
3. Pipeline de Inferência:
   • Morfologia: Imagens mock de NIRCam são ajustadas com pysersic para derivar priors morfológicos (raio efetivo $r_e$, índice de Sérsic $n_s$, razão de eixo $q$, ângulo de posição PA).
   • Cinemática: Espectros mock de NIRSpec são ajustados usando msafit, assumindo um modelo de disco fino e rotativo com um perfil de velocidade de arcotangente e dispersão de velocidade constante. Os parâmetros morfológicos derivados da imagem são usados como priors fortes. O ajuste emprega uma minimização de $\chi^2$ ponderada, incluindo uma escala de fluxo por linha não paramétrica para lidar com distribuições de luz clumpy.

Contribuições Principentes e Resultados
O estudo compara os parâmetros cinemáticos de melhor ajuste derivados dos espectros mock contra os valores intrínsecos da simulação para quantificar vieses observacionais.

• Precisão de Recuperação: Para fontes alinhadas e resolvidas em $z=3$, as velocidades rotacionais inferidas ($v$) e as dispersões de velocidade ($\sigma$) geralmente concordam com as medições intrínsecas dentro de um fator de $\sim 2$ e $\sim 1.5$, respectivamente. O desvio (scatter) em $v$ é maior que em $\sigma$.
• Vieses Sistemáticos:
  • Dependência de Morfologia: A recuperação é robusta para sistemas suaves e de aspecto discoide. No entanto, para sistemas intrinsecamente alongados ou esferoidais (grandes razões de eixo intrínsecas $C/A$), $v$ tende a ser superestimado e $\sigma$ subestimado por mais de um fator de 2. Isso ocorre porque a suposição de disco fino força o modelo a interpretar o alargamento não rotacional como efeitos de inclinação.
  • Evolução com Redshift: O desvio na precisão de recuperação aumenta em redshifts mais altos ($z > 3$) à medida que as galáxias se desviam mais fortemente da suposição de disco fino.
  • Individual vs. População: Embora a cinemática de galáxias individuais ($v/\sigma$) possa ser classificada incorretamente por um fator de 2–3 (ex: um sistema suportado por dispersão parecendo suportado por rotação), as tendências populacionais são robustamente recuperadas. A fraca evolução com o redshift de $\sigma$ e $v/\sigma$, bem como a dependência de $\sigma$ com a SFR global, são preservadas nos dados inferidos por mock.
• Relações de Escala: A relação inferida $\sigma$-SFR alinha-se bem com os modelos teóricos de "transporte de massa + feedback" (Krumholz et al. 2018), sugerindo que, apesar das incertezas individuais, a modelagem cinemática pode restringir de forma confiável os drivers físicos da turbulência (feedback vs. gravidade) em populações de galáxias.
• Modos de Falha: Aproximadamente 5–10% das fontes falham em ajustar o modelo de disco fino, frequentemente devido a morfologias irregulares, sistemas de alta massa estelar com gradientes de velocidade centrais íngremes não capturados por perfis de arcotangente, ou fontes onde o $\sigma$ intrínseco está abaixo do limite inferior do modelo ($10$ km s$^{-1}$).

Significância e Alegações
O artigo afirma fornecer o primeiro framework de ponta a ponta para comparar observações de NIRSpec MSA com simulações cosmológicas para quantificar vieses na medição da cinemática de gás ionizado. Os autores concluem que, embora modelos simples de disco fino sejam insuficientes para determinar a "aspecto de disco" ($v/\sigma$) preciso de galáxias individuais de alto redshift devido às grandes incertezas, eles são eficazes para recuperar tendências estatísticas de população.

O trabalho destaca que as inferências observacionais de $v/\sigma$ para sistemas individuais em $z \ge 2$ devem ser tratadas com cautela, pois podem ser enviesadas pela suposição de geometria de disco fino e pela presença de movimentos não circulares. Os autores argumentam que, para melhorar as restrições sobre o estado dinâmico de galáxias precoces, trabalhos futuros devem ir além da espectroscopia de fenda 2D em direção a observações de IFU 3D ou estratégias de degraus de fenda (slit-stepping), e desenvolver modelos cinemáticos dedicados que deem conta da complexidade morfológica 3D e da natureza clumpy inerente às galáxias de alto redshift. O estudo enfatiza que as limitações atuais não são meramente instrumentais, mas decorrem do descompasso entre modelos paramétricos simples e a complexidade intrínseca da física de formação de galáxias.