On gamma rays as predictors of UHECR flux in AGNs

技术摘要：以 $\gamma$ 射线作为 AGN 中 UHECR 通量的预测因子

问题陈述
活动星系核（AGNs）是超高能宇宙线（UHECRs）的主要候选来源。然而，最近一些使用观测到的极高能（VHE）$\gamma$ 射线通量作为 UHECR 光度（$L_{CR}$）直接代理指标的唯象模型，无法重现皮埃尔·奥格天文台（Pierre Auger Observatory）的实验数据。具体而言，这些模型难以同时拟合能量谱、质量组成以及到达方向图。一个关键的失效点在于高估了耀变体 Markarian 421（Mkn 421）的贡献，这产生了一个与观测结果不一致的预测热点和偶极方向。此外，这些模型通常会产生“伪”偶极方向以及较差的光谱形状拟合，特别是在对不同 AGN 类群（例如射电星系和耀变体）进行混合分析时。

方法论
作者认为，这种差异源于先前研究中的一个隐含假设：即 UHECR 的发射与来自相对论性喷流的 $\gamma$ 射线发射具有相同的束缚性（beaming）。论文指出，由于相对论性运动在喷流内的“团块”（blobs）中，$\gamma$ 射线受到强烈的多普勒增益（$D^4$）影响，而 UHECR 在逃逸出源环境之前，很可能经历了各向同性化（去准直化/decollimation）。

本研究通过以下步骤进行：

1. 理论回顾： 作者回顾了 AGN 喷流中的加速机制（激波、重联、剪切加速）以及粒子在喷流瓣（lobes）中的传播过程。他们认为，瓣中的磁场和源包层（cocoon）中的磁场会散射 UHECR，从而有效地使它们的发射趋于各向同性，即使它们起源于相对论性流体中。
2. 代理修正： 作者建议不再使用观测到的 $\gamma$ 射线光度（$L^{obs}_{\gamma}$），而是使用内禀 $\gamma$ 射线光度（$L^{int}_{\gamma}$）或辐射喷流功率作为 $L_{CR}$ 的代理。这涉及通过多普勒因子（$D$）的 $q$ 次幂来修正观测通量，即 $L^{int}_{\gamma} = D^{-q}L^{obs}_{\gamma}$。文中测试了两种情景：$q=4$（假设 $L_{CR} \propto L^{int}_{\gamma}$）和 $q=2$（假设 $L_{CR} \propto$ 辐射喷流功率）。
3. 模拟与拟合： 利用 CRPropa3 框架，作者针对五种主要原子核（H, He, N, Si, Fe），在不同的能量和距离网格上进行了一维传播模拟。他们将模拟的通量与皮埃尔·奥格观测数据（能谱和 $X_{max}$ 矩）及到达方向图进行了对比。
4. 对比分析： 对三种情况进行了比较：
   • 情况 (a)：使用 $L^{obs}_{\gamma}$ 作为代理（标准方法）。
   • 情况 (b)：使用 $D^{-2}L^{obs}_{\gamma}$ 作为代理。
   • 情况 (c)：使用 $D^{-4}L^{obs}_{\gamma}$ 作为代理。

主要贡献与结果

• 抑制耀变体的支配地位： 修正显著降低了高度束缚的耀变体（如 Mkn 421）相对于射电星系（如 Cen A, M87, 和 Fornax A）的权重。在标准模型中，Mkn 421 主导了通量；而在修正后的模型中，其贡献被削减了 2 到 5 个数量级，使其变得微不足道。
• 改善光谱与组成拟合： 当使用内禀光度代理时，能量谱和组成数据的综合拟合效果显著提升。归一化卡方值（$\chi^2/NDF$）从 4.6（使用 $L^{obs}_{\gamma}$）降至 3.1（使用 $D^{-2}L^{obs}_{\gamma}$）。这种改进归功于更真实的近源与远源之间的相对贡献比例，避免了由过度贡献的远距离、高束缚源导致的谱形畸变。
• 解决偶极子张力： 模型与观测偶极方向之间的张力得到了缓解。对于能量 $E > 8$ EeV 的奥格视场，差异从 5.9$\sigma$ 降至 3.5$\sigma$。对于 $E > 32$ EeV，张力从 2.1$\sigma$ 降至 1.1$\sigma$。
• 热点重定向： 预测的到达方向图从由 Mkn 421 主导，转变为围绕 Cen A, Fornax A, M87 以及特定远距离耀变体（Mkn 180, 1ES2344+514）显示热点。这些新的热点与皮埃尔·奥格天文台报告的两个热点以及望远镜阵列（Telescope Array）报告的一个热点更为吻合。

意义与主张
本文声称为先前基于 AGN 的 UHECR 源模型中所存在的矛盾提供了解决方案，其核心在于解耦了 $\gamma$ 射线与 UHECR 的束缚性假设。作者断言：

1. $\gamma$ 射线与 UHECR 之间的关联是微弱但非零的；虽然 $\gamma$ 射线追踪粒子加速区域，但 UHECR 在逃逸前很可能被源瓣中的磁场各向同性化。
2. 使用观测到的 $\gamma$ 射线通量作为代理，隐含地假设了 UHECR 是被束缚的，而作者认为这一假设对于逃逸复杂磁环境的 UHECR 来说在物理上是不太可能的。
3. 通过多普勒增益进行修正（使用 $L^{int}_{\gamma}$ 或辐射喷流功率）提供了一个更稳健的代理，能够调和来自不同 AGN 类别的皮埃尔·奥格天文台的能谱、组成及到达方向数据。

作者总结道，尽管他们的模型是一个简化的综合拟合（未考虑完整的 $X_{max}$ 分布或详细的磁场偏转），但结果强化了 AGNs 作为 UHECR 来源的假设，并提供了一个更一致的框架来解释来自不同 AGN 类别的多信使数据。他们指出，未来的工作必须解决能量依赖性的去准直化以及正则性银河系和河外磁场的具体影响。