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

技术摘要：测试从模拟的 JWST NIRSpec/MSA 观测中推断运动学的方法

问题陈述
理解星系的动力学演化需要测量电离气体在宇宙时标上的旋转速度 ($v$) 和速度弥散 ($\sigma$)。虽然像 IllustrisTNG 这样的流体动力学模拟提供了六维空间和运动学数据，但它们缺乏针对特定观测设置的显式辐射传输解。相反，使用诸如 JWST 的 NIRSpec 微分光计（MSA）等仪器的观测研究依赖于前向建模来从二维光谱中推断运动学。这些模型通常假设星系是轴对称、薄的旋转盘，并具有 Sérsic 强度分布和反正切速度曲线。然而，高红移（$z > 2$）星系往往表现出复杂的形态（类球状、团块状）和非圆周运动（流入/流出），这可能违反这些假设。目前极度缺乏统计学验证，以确定简单的参数化模型是否能够准确恢复典型高红移星系的内在运动学，特别是考虑到与积分场单元（IFU）相比，MSA 观测具有有限的空间分辨率和狭缝式观测的特性。

方法论
作者利用高分辨率 TNG50 宇宙学模拟生成了一个端到端的模拟观测流水线，涵盖了超过 4,000 个处于 $z = 2, 3, 4, 6$ 的恒星形成星系（$M_* = 10^8 - 10^{11.5} M_\odot$）。

1. 内在测量： 作者在 H$\alpha$ 发射气体有效半径 ($r_e$) 内定义了“内在”运动学量（$v_{max}$ 和 $\sigma_e$）。他们明确测试了这些内在值对测量选择（包括孔径大小、空间分箱 [0.5 vs. 5 ckpc] 以及权重方案 [无权重 vs. 星形成率加权]）的敏感性。
2. 模拟观测：
   • 成像： 通过投影恒星形成气体单元，通过 Kennicutt (1998) 关系将恒星形成率 (SFR) 转换为 H$\alpha$ 光度，并与 NIRCam PSF 进行卷积，生成模拟 NIRCam 图像。加入噪声以模拟深度 JADES 级的曝光。
   • 光谱： 通过沿随机轴投影气体单元，分配视线速度，并与 NIRSpec 线扩散函数 (LSF) 和 PSF 进行卷积，创建模拟 NIRSpec/MSA 二维光谱。应用高斯噪声以达到与现有的最深 JWST 数据相当的信噪比 (S/N $\sim$ 50)。
3. 推断流水线：
   • 形态： 使用 pysersic 对模拟 NIRCam 图像进行拟合，以导出形态先验（有效半径 $r_e$、Sérsic 指数 $n_s$、轴比 $q$、位置角 PA）。
   • 运动学： 使用 msafit 对模拟 NIRSpec 光谱进行拟合，假设一个具有反正切速度剖面和恒定速度弥散的薄旋转盘模型。从成像中导出的形态参数被用作强先验。拟合过程采用加权 $\chi^2$ 最小化，包括一个非参数化的行向通量缩放，以解释团块状的光分布。

主要贡献与结果
本研究将从模拟光谱中推导出的最佳拟合运动学参数与内在模拟值进行比较，以量化观测偏差。

• 恢复准确度： 对于在 $z=3$ 时对齐且解析的源，推断的旋转速度 ($v$) 和速度弥散 ($\sigma$) 通常在内在测量值的 $\sim 2$ 倍和 $\sim 1.5$ 倍因子范围内一致。$v$ 的散射比 $\sigma$ 更大。
• 系统性偏差：
  • 形态依赖性： 对于平滑、盘状系统，恢复是稳健的。然而，对于本质上呈拉长状或类球状的系统（较大的内在轴比 $C/A$），$v$ 往往被高估，而 $\sigma$ 被低估超过 2 倍。这是因为薄盘假设迫使模型将非旋转展宽解释为倾角效应。
  • 红移演化： 随着星系更强烈地偏离薄盘假设，恢复准确度的散射在更高红移（$z > 3$）处增加。
  • 个体 vs. 总体： 虽然单个星系的运动学 ($v/\sigma$) 可能会被错误分类（例如，一个弥散支持系统被误认为旋转支持系统），但总体趋势是稳健恢复的。$\sigma$ 和 $v/\sigma$ 的微弱红移演化，以及 $\sigma$ 对全局 SFR 的依赖性，在模拟推断的数据中得到了保留。
• 标度关系： 推断的 $\sigma$-SFR 关系与“质量-输运 + 反馈”理论模型（Krumholz 等人 2018）吻合良好，这表明尽管存在个体不确定性，运动学建模可以可靠地约束星系族群中湍流的物理驱动因素（反馈 vs. 引力）。
• 失效模式： 大约 5–10% 的源无法拟合薄盘模型，这通常是由于不规则形态、具有无法被反正切剖面捕捉到的陡峭中心速度梯度的重质量星系，或者是内在 $\sigma$ 低于模型下限（$10$ km s$^{-1}$）的源。

意义与主张
本文声称提供了第一个将 NIRSpec MSA 观测与宇宙学模拟进行对比以量化测量电离气体运动学偏差的端到端框架。作者得出结论，虽然简单的薄盘模型不足以确定单个高红移星系的精确“盘状程度”（$v/\sigma$），但它们对于恢复统计学上的族群趋势是有效的。

这项工作强调，对于 $z \ge 2$ 的单个系统，其 $v/\sigma$ 的观测推断应谨慎对待，因为它们可能受到薄盘几何假设和非圆周运动存在的偏差影响。作者认为，为了改进对早期星系动力学状态的约束，未来的工作必须超越二维狭缝光谱，转向三维 IFU 观测或狭缝步进策略，并开发能够处理高红极星系中固有的复杂三维形态和团块性的专用运动学模型。该研究强调，目前的局限性不仅在于仪器层面，更源于简单参数化模型与复杂的星系形成物理过程之间的不匹配。