Constraining reionization-era Ly$\alpha$ escape with JELS-MUSE: a highly complete H$\alpha$-selected sample at $z\sim6.1$

技术摘要：利用 JELS-MUSE 约束再电离时期的 Lyα 逃逸

问题与动机
再电离时期（EoR）标志着星际介质（IGM）从中性向电离态的转变过程，这一过程主要由恒星形成星系驱动。理解这一时期的关键量是莱曼连续谱（LyC）逃逸份额（$f_{\text{LyC}}^{\text{esc}}$），它决定了注入 IGM 的电离光子速率。然而，由于 IGM 的不透明度，在 $z \gtrsim 4.5$ 时无法直接观测到 LyC。因此，莱曼阿尔法（Ly$\alpha$）逃逸份额（$f_{\text{Ly}\alpha}^{\text{esc}}$）成为了一个主要的间接示踪指标，因为在关于中性气体柱密度和尘埃几何结构方面，Ly$\alpha$ 光子与 LyC 光子具有相似的逃逸物理机制。

在高红移条件下，表征 $f_{\texttext{Ly}\alpha}^{\text{esc}}$ 的一个重大挑战在于区分星系星际介质（ISM）条件的内在变化与 IGM 引起的衰减。以往的研究通常依赖于 Ly$\alpha$ 选择或 UV 选择的样本，这会引入向具有有利逃逸条件或特定亮度范围的星系倾斜的选择偏差。此外，在 $z \gtrsim 6.5$ 时，随着 IGM 中性比例的增加，Ly$\alpha$ 的可见度受到抑制，使得隔离内在 ISM 驱动因素变得困难。因此，需要一个在过渡红移（$z \sim 6$）处高度完整且无偏的样本，以在 IGM 抑制占据主导地位之前，建立一个内在 Ly$\alpha$ 逃逸的基准。

方法论
本研究利用了 JELS-MUSE 大面积巡天数据，结合了詹姆斯·韦伯空间望远镜（JWST）和甚大望远镜（VLT）/MUSE 的数据。

1. 样本选择： 母样本由 24 个位于 $z \approx 6.1$ 的恒星形成星系组成，这些星系完全通过 **JWST 发射线巡天（JELS）**中的 H$\alpha$ 发射进行选择。这种窄带成像方法（使用 F466N 和 F470N 滤光片）确保了在 $\log(L_{\text{H}\alpha}) \gtrsim 41.5$ erg s$^{-1}$ 以下的恒星形成星系的通量限制完整普查，避免了 Ly$\alpha$ 或 UV 连续谱选择所固有的偏差。
2. 光谱学： Ly$\alpha$ 光谱学通过 VLT/MUSE 在 COSMOS 区域内一个连续的 63 arcmin$^2$ 场内获得，该区域覆盖了 JELS 足迹的相同区域。这为所有 24 个经 H$\alpha$ 选择的源提供了积分场 Ly$\alpha$ 覆盖。
3. 数据处理：
   • 连续谱减除： 对 MUSE 数据立方体应用了光谱中值滤波方法，以在保留谱线通量的同时去除连续谱发射。
   • 提取： 通过生长曲线分析确定了一个 $r = 0.6$ arcsec（在 $z=6.1$ 时约为 $3\text{--}4$ kpc）的最佳空间孔径，以最大化信噪比（SNR），同时最小化来自来源不明的扩展低表面亮度晕的污染。
   • 探测阈值： 通过对 500 个随机空间位置进行假阳性分析，建立了经验探测阈值 $\Sigma_{\text{Ly}\alpha} > 5.09$（第 98 百分位数）用于探测，以及 $\Sigma_{\text{Ly}\alpha} > 6.41$（第 99.5 百分位数）用于强探测。
   • 通量测量： 通过拟合偏斜高斯轮廓来测量 Ly$\alpha$ 通量，以解释共振散射的不对称性。对于未探测到的目标，计算了保守的 $5\sigma$ 上限。
4. 逃逸份额计算： $f_{\text{Ly}\alpha}^{\text{esc}}$ 被计算为观测到的 Ly$\alpha$ 光度与经过尘埃校正后的 H$\alpha$ 光度的比值，假设采用 Case B 再结合比例 8.7。基于 SED 推导的 $A_V$ 对 H$\alpha$ 进行尘埃校正，而 Ly$\alpha$ 则保持不校正，以反映光子逃逸星系的总体份额（即损失于尘埃或被散射出视线的比例）。
5. 统计分析： 利用逆卡普兰-梅耶（reverse Kaplan-Meier）生存分析（以纳入非探测值/上限）和堆叠通量法得出群体平均逃逸份额。使用 Spearman 等级和受限 Kendall $\tau$ 统计量评估了其与星系属性（UV 星等、斜率、恒星质量、尘埃消光）的相关性。

关键结果

• 探测比例： 在 24 个源中有 12 个探测到 Ly$\alpha$ 发射（$50 \pm 10\%$）。在亮度范围 $-20.25 < M_{\text{UV}} < -18.75$ 内，使用经典的 Ly$\alpha$ 辐射体定义（EW(Ly$\alpha$) > 25 Å），Ly$\alpha$ 辐射体比例为 $X_{\text{Ly}\alpha} = 33 \pm 12\%$。
• 逃逸份额分布： 被探测到的源表现出广泛的逃逸份额范围，$f_{\text{Ly}\alpha}^{\text{esc}} = 0.04\text{--}0.83$。
• 群体平均逃逸份额： 通过逆卡普兰-梅耶分析纳入非探测值后，得到的其中值逃逸份额为 $\langle f_{\text{Ly}\alpha}^{\text{esc}} \rangle = 0.07^{+0.04}_{-0.03}$。这与独立的堆叠通量估计值 $0.08^{+0.02}_{-0.02}$ 一致。
• IGM 均匀性： 与巡天区域和红移匹配的再电离子模拟表明，该区域内的 IGM 透射率大致均匀且高度趋向于 1。因此，观察到的 $f_{\text{Ly}\alpha}^{\text{esc}}$ 离散度归因于内在的 ISM 和星系周介质（CGM）的变化，而非 IGM 的不均匀性。
• 驱动因素：
  • 全局属性： 使用 Kolmogorov-Smirnov 统计量测试发现，在 UV 星等（$M_{\text{UV}}$）、恒星质量（$M_*$）、恒星形成率（SFR）或尘埃消光（$A_V$）方面，Ly$\alpha$ 探测与非探测之间没有显著差异。
  • 相关性： 定性趋势表明，具有较高 $f_{\text{Ly}\alpha}^{\text{esc}}$ 的星系往往具有较低的星际尘埃消光（$E(B-V)$）、较蓝的 UV 斜率（$\beta$）和较低的恒星质量。然而，由于样本量和内在离散度，这些相关性在单独测试时并不具备高置信度的显著性。
  • 离散度： 在具有相似全局属性的星系中，$f_{\text{Ly}\alpha}^{\text{esc}}$ 的巨大离散度意味着 Ly$\alpha$ 的逃逸受控于小尺度的、各向异性的中性气体几何结构（例如“栅栏式”通道）和动力学，而这些特征是集成宽带光度诊断无法捕捉到的。

意义与主张
本文声称，这项工作为再电离子结束阶段（$z \approx 6.1$）的内在 Ly$\alpha$ 逃逸提供了一个稳健的经验基准。通过使用 H$\alpha$ 选择的样本，本研究避免了困扰 Ly$\alpha$ 选择类巡天的选择偏差，从而提供了对恒星形成星系种群更具代表性的观察。

作者认为，在这一单场、窄红移样本中观察到的显著 $f_{\text{Ly}\alpha}^{\text{esc}}$ 离散度是星系 ISM 条件的内在特性。这一发现具有重要意义，因为它确立了在 $z \gtrsim 6.5$ 时观察到的 Ly$\alpha$ 可见度下降不能仅仅归因于 IGM 的抑制；很大一部分变率是由 ISM 逃逸通道的随机性质驱动的。因此，解释更高红移处的 Ly$\alpha$ 观测时，需要考虑这种内在的变率。

测得的约 7% 的群体平均逃逸份额与近期基于 JWST 的 H$\alpha$ 锚定研究一致，但低于基于光度函数的估计值或 Ly$\alpha$ 选择类样本，这表明以往的方法可能由于选择偏差而高估了典型的逃逸份额。结果支持这样一种情景：一小部分具有有利、瞬态 ISM 条件的星系（可能由反馈调节的低柱密度通道驱动）主导了电离光子的预算，而非一个均匀的高效漏光种群。