Analysis and implications of the spatio-spectral morphology of the Fermi Bubbles

技术摘要：费米泡（Fermi Bubbles）空间-光谱形态的分析及其意义

问题陈述
费米泡（FBs）——即在银河系平面上下延伸约 50° 的巨型 $\gamma$ 射线结构——其起源和能量来源仍不确定。虽然其形态表明起源于银河系中心（GC），但驱动 $\gamma$ 射线发射的微观物理机制在两种主要通道之间存在争议：强子通道（宇宙线质子与环境气体的非弹性碰撞）和轻子通道（宇宙线电子对星际辐射场光子的逆 컴普顿散射）。以往的分析依赖于粗略的空间表示（例如，纬度平均切片），这可能掩盖了位置相关的光谱变化。此外，粒子输运和冷却所需的时间尺度对于判断这些气泡是由活动星系核（AGN）活动还是恒星形成爆发驱动，提出了显著的约束。本文旨在通过进行高分辨率、逐像素的光谱分析，来约束重建观测到的 $\gamma$ 射线光谱所需的母体宇宙线（CR）种群。

方法论
作者利用了 Platz 等人（2023）提供的十年期全天 Fermi/LAT 数据的无模板重建。具体而言，他们采用了 M2 模型中的“非尘埃”（$I_{nd}$）组分，该组分隔离了与热尘埃图无关的发射，从而在银河系前景中清晰地展现出费米泡。天空使用 HEALPix 方案进行像素化（$N_{side}=128$），产生了约 20,000 个覆盖费米泡的像素。

为了模拟发射，作者：

1. 构建 3D 几何结构： 他们将费米泡近似为通过绕银河系 $z$ 轴旋转眼迹边界而生成的旋转体。这使得计算视线深度以及沿每条射线积分物理量成为可能。
2. 模拟星际辐射场（ISRF）： 他们使用 Popescu 等人（2017）的轴对称辐射传输模型，来确定每条视线上位置相关的 ISRF（包括恒星、尘埃和 CMB 组件）。对于超出模型网格的区域（$z > 10$ kpc），他们使用幂律或 $z^{-1}$ 标度进行外推。
3. 前向建模： 使用 Naima Python 工具包，将每个像素中观测到的 $\gamma$ 射线光谱与三种母体种群模型进行拟合：幂律（PL）、指数截断幂律（EPL）和断裂幂律（BPL）。
4. 强子与轻子： 对于强子模型，他们假设一个典型的均匀氢数密度 $n_H = 10^{-3} \text{ cm}^{-3}$ 并据此归一化结果。对于轻子模型，他们计算针对每个像素平均 ISRF 的逆 컴普顿（IC）散射。
5. 拟合优度： 他们使用降低的 $\chi^2$ 来评估模型，并指出由于数据分解中存在相关不确定性，该指标仅用于比较相对模型的性能，而非绝对统计显著性。

关键结果

• 光谱形状： 对于质子和电子，简单的幂律模型（PL）对数据的拟合效果较差，并已被排除。EPL 和 BPL 模型在强子和轻子情景下均提供了同样良好的拟合。
• 谱指数与硬化： 在低能端，对于宇宙线质子（CRp）和宇宙线电子（CRe），最佳拟合谱指数在气泡南端表现出明显的硬化趋势。这一趋势与之前的发现一致，但涵盖了比以往分析更高的纬度。
• 断裂/截断能量：
  • 对于 EPL 模型，截断能量的数量级为几个 TeV，且随纬度变化很小。
    · 对于 BPL 模型，断裂能量约为 100 GeV，且纬度变化极小。
  • BPL 模型中跨越断裂处的谱指数变化 $\Delta \alpha \lesssim 0.5$。这显著小于由于同步辐射或 IC 损失导致的预期冷却断裂（$\Delta \alpha \approx 1$），表明这些光谱断裂并非冷却特征。
• 辐射场贡献： 在轻子模型中，在所有 Fermi/LAT 能量下，ISRF 的尘埃组分在 $\gamma$ 射线产生中占主导地位，其次是恒星组分（由于克莱因-尼尔森效应在高能端受到抑制）和 CMB。因此，在整个气泡内，IC 损失可以近似为在汤姆逊（Thomson）机制下发生。
• 能量密度与压力：
  • 强子： 总宇宙线能量估计为 $\sim 1.7 \times 10^{57} / n_{H,-3}$ erg。宇宙线压力可能主导热气体压力，但仅高出几倍，且高度依赖于假设的目标气体密度。
  • 轻子： 宇宙线电子能量密度随距离银河系平面的距离增加而增加，并在气泡边缘达到峰值，这在定性上与边缘增亮现象一致。
• 冷却时间尺度：
  • 强子： 冷却时间尺度极长（对于 $n_H = 10^{-3} \text{ cm}^{-3}$，约为 $\sim 100$ Gyr），这意味着气泡可能非常古老，或者远离稳态。
  • 轻子： 高能电子的冷却时间极短（在气泡顶端 $\lesssim 1$ Myr，在靠近银河系中心处降至 $\sim 0.1$ Myr）。

意义与启示
论文得出结论，仅凭光谱形状无法明确排除强子或轻子模型，因为 EPL 和 BPL 模型对两种粒子类型都能提供同样良好的拟合。然而，结果对物理情景施加了严格的约束：

1. 轻子输运 vs. 原位加速： 短冷却时间与观测到的向气泡边缘增加的 CRe 能量密度的结合，对“电子从银河系中心注入并向外输运”的模型提出了严峻挑战。在 $\sim 1$ Myr 的冷却时间内将电子输运到气泡顶端（$\sim 8$ kpc），需要约 8000 km/s 的膨胀速度，这与 X 射线现象学不符。作者认为，**原位再加速（in-situ re-acceleration）**轻子在气泡边缘的模型，是解释观测到的形态和能量密度分布更自然的解释。
2. 强子年龄约束： 长时间的强子冷却时间意味着，如果气泡很年轻（例如 $< 100$ Myr），所需的能量预算将需要 AGN 起源而非恒星形成。相反，如果气泡很古老（数个 Gyr），只要宇宙线种群有足够的时间积累，恒星形成起源仍然是合理的。
3. 光谱硬化： 观测到的向南端硬化的趋势以及缺乏冷却断裂（$\Delta \alpha < 0.5$）表明，光谱特征是由注入或加速机制驱动的，而非辐射冷却损失。

总之，这项工作提供了对费米泡进行首次逐像素光谱约束，证明了虽然简单的幂律是不够的，但带有指数截断或断裂的模型能很好地拟合数据。研究结果强烈支持以下情景：要么是来自核区的相对论性输运，或者更有可能是宇宙线在气泡边缘的分布式原位加速，特别是在轻子通道方面。