On the Origin of the Ly$\alpha$ Damping Wing in Galaxies at $8\le z \le 10$: Explorations using the NINJA Simulations

技术摘要：关于 $8 \le z \le 10$ 星系中 Ly$\alpha$ 阻尼翼起源的研究

问题陈述
再电离时期（EoR）的特征是星际介质（IGM）从中性向电离态的转变。虽然通过詹姆斯·韦伯空间望远镜（JWST）在观测高红移（$7 \le z \le 14$）星系光谱时发现的 Ly$\alpha$ 阻尼翼为这一时期提供了强有力的探测手段，但解释这些特征非常复杂。阻尼翼的形状具有简并性，是由弥散 IGM、星系周介质（CGM）以及源星系自身的星际介质（ISM）的共同吸收作用引起的。目前的观测显示，存在高发生率的具有高中性氢柱密度（$N_{\text{H I}}$）的阻尼 Ly$\alpha$（DLA）系统，但目前尚不清楚这些特征在多大程度上是由全球 IGM 中性分率（$x_{\text{HI}}$）驱动的，还是由星系周围的局部气体库驱动的。现有的模型往往难以在缺乏对多个尺度下气体电离状态进行物理动机约束的情况下，同时重现观测到的强吸收器发生率和中值柱密度。

方法论
作者利用了 Ninja (NISER) 系列大规模宇宙学流体动力学模拟，该模拟包含了恒星形成、恒星反馈和黑洞生长，但不包含自洽的辐射传输（RT）。为了模拟 EoR 环境，作者采用了一种理想化的后处理方法，将这些模拟与电离泡的玩具模型相结合。

研究生成了 $z=8$ 和 $z=10$ 时刻的模拟 Ly$\alpha$ 吸收光谱，采用了三种不同的情景来分离不同气体组分的贡献：

1. 情况 I： 吸收仅源于均匀电离（或部分电离）的弥散 IGM，具有全局中性分率 $x_{\text{HI}}$。
2. 情况 II： 星系嵌入在半径为 $R_{\text{H II}}$（范围从 0 到 400 pkpc）的球形完全电离 H II 区内。吸收包括来自 IGM 以及泡外星系周围过密区域的气体贡献，但假设泡内没有中性气体。
3. 情况 III： 情况 II 的一种修正版本，在哈洛维里尔半径（$R_{\text{vir}}$）内加入了具有特定中性分率（$x^{\text{V}}_{\text{HI}}$）的残余中性气体，而 $R_{\text{vir}}$ 与 $R_{\text{H II}}$ 之间的气体保持完全电离。

模拟使用了三个不同分辨率和体积的盒子（$L50N640$、$L50N1008$ 和 $L140N1008$）来评估数值收敛性和哈洛质量依赖性。模拟光谱使用高斯轮廓进行卷积，以匹配 JWST/NIRSpec PRISM 模式的低光谱分辨率（$R \sim 100-1000$），并包含泊松噪声（SNR=30）。吸收强度通过参数 $x_p$（10% 透射率处的速度间隔）进行量化，该参数被映射为等效中性氢柱密度（$N_{\text{H I}}$）。

主要贡献与结果
本文系统地量化了 IGM、CGM 和 ISM 的贡献如何塑造 Ly$\alpha$ 阻尼翼：

• 气体过密度与哈洛质量： 模拟表明，星系被延伸至 $\sim 60-100$ pkpc 的过密气体所包围。这一过密区域的范围（$R_{\text{eq}}$）随哈洛质量而变化（对于低质量哈洛，其规模从 $\sim 66$ pkpc 增加到高质量哈洛的 $\sim 105$ pkpc），但在 $z=8$ 到 $z=10$ 之间几乎没有演化。
• IGM 贡献（情况 I）： 仅具有均匀电离 IGM 或低中性分率（$x_{\text{HI}} = 0.1$）的模型无法重现观测到的强 DLA 系统。即使在完全中性的 IGM（$x_{\text{HI}} = 1.0$）情况下，模拟的 IGM 视线也极少超过 $N_{\text{H I}} > 10^{22} \text{ cm}^{-2}$。这表明仅靠 IGM 无法解释在 JWST 数据中观测到的最强吸收体。
• H II 区的作用（情况 II）： 电离泡的存在显著改变了吸收统计特性。
  • 如果 $R_{\text{H II}} \approx 0$（无电离泡），模型预测超过 90% 的星系具有 $N_{\text{H I}} > 10^{21} \text{ cm}^{-2}$，这超过了观测约束，意味着在相当一部分星系周围必须存在有限大小的电离区。
  • 相反，如果 $R_{\text{H II}} > 50$ pkpc，即使在存在中性 IGM 的情况下，模型也无法重现最强的吸收体（$N_{\text{H I}} > 10^{22} \text{ cm}^{-2}$）。这表明大型电离泡过度抑制了阻尼翼，导致无法匹配观测到的强 DLA。
• 残余中性气体（情况 III）： 在维里半径内引入残余中性气体（$x^{\text{V}}_{\text{HI}} > 0$）对于产生最强的阻尼 Ly$\alpha$ 吸收器是必要的。然而，即使在维里半径内是完全中性的气体，预测的 $N_{\text{H I}} > 10^{22} \text{ cm}^{-2}$ 系统的发生率仍低于模拟中所解析的哈洛质量的观测估计值。
• 质量依赖性： 推断的 $N_{\text{H I}}$ 与恒星质量和哈洛质量均呈正相关。由于拥有更大的气体库和更密的物理环境，较强的吸收优先出现在更高质量的哈洛中。

意义与结论
本文得出结论，在 $z \gtrsim 8$ 的星系中观测到的 Ly$\alpha$ 阻尼翼是 IGM、CGM 和 ISM 吸收共同作用的复合特征。作者断言：

1. IGM 的局限性： 单纯依靠弥散 IGM 的中性分率（$x_{\text{HI}}$）无法解释 JWST 观测到的高发生率强阻尼 Ly$\alpha$ 吸收体。
2. 局部气体的必要性： 要重现最强的吸收体，需要来自星系维里半径内及其紧邻过密环境的中性气体贡献。
3. 缺失的物理机制： 模拟中强吸收器发生率（即使考虑了维里半径内的中性气体）与观测值之间的差异表明，观测到的星系可能居住在比当前模拟所解析的更高质量的哈洛中，或者额外的未解析中性气体（例如在恒星诞生云中）对吸收贡献显著。
4. 方法论要求： 对 IGM 中性分率（$x_{\text{HI}}$）的稳健约束需要能够同时考虑 ISM、CGM 和 IGM 电离状态的统一模型。本研究强调了将小尺度辐射传输计算与大规模再电离模拟相结合的必要性，以充分解释低分辨率的 JWST 光谱。

作者强调，其结果在不同的模拟分辨率和体积下在数值上是稳健的，为未来整合多尺度辐射传输以解决 EoR 观测中简并问题的研究提供了基准。