Study of Flat Spectrum Radio Quasars and BL Lacertae Objects as Sources of Diffusive Ultra High-Energy Cosmic Rays

技术摘要：平坦谱射电类星体与 BL Lac 天体作为扩散超高能宇宙线源的研究

问题陈述
超高能宇宙线（UHECRs，$E \gtrsim 10^{17}$ eV）的起源仍是现代天体物理学中一个重大的未解之谜。尽管谱特征暗示了其河外起源，但由于银河系和河外磁场的偏转作用，精确识别其源变得十分困难。皮埃尔·奥热天文台（Pierre Auger Observatory）和望远镜阵列（Telescope Array, TA）近期的高精度数据揭示了特定的谱特征（在“踝”处的硬化、变陡以及 GZK 截断）和成分趋势（随着能量增加，从轻核向重核转变）。虽然邻近的射电星系是强有力的候选者，但活动星系核（AGN），特别是平坦谱射电类星体（FSRQs）和 BL Lacertae（BL Lac）天体的潜在贡献，需要针对通量、成分和各向异性的观测约束进行严格检验。

方法论
作者利用标准 $\Lambda$CDM 宇宙学下的综合传播模型，研究 FSRQs 和 BL Lacs 作为扩散 UHECRs 的合理源。

1. 源演化模型：

   • 本研究利用从 Fermi-LAT 数据观测到的伽马射线光度函数（LF）导出的真实光度依赖密度演化（LDDE）函数。
   • FSRQs： 采用 LDDE 函数建模，其中源密度的红移峰值依赖于光度（$z_c \propto L^\alpha$）。
   • BL Lacs： 采用与光度相关的光子指数内禀分布（高斯分布）进行建模，并结合相同的 LDDE 形式。
   • 这些模型用基于观测伽马射线数据的物理驱动演化，取代了简单的唯象演化形式（例如 $1+z$ 的幂律）。

2. 传播物理：

   • 考虑与宇宙微波背景（CMB）和河外背景光（EBL）相互作用导致的能量损失，将 UHECRs 在河外空间中传播。
   • 包含的过程有：Bethe-Heitler 对产生、光致π介子产生（针对质子）和光致裂变（针对原子核）。本研究采用连续能量损失（CEL）近似处理光致裂变，追踪初级原子核（p, He, N, Si, Fe）的存活概率，而不显式求解次级级联。
   • 磁扩散： 传播发生在均匀、湍流的河外磁场（IGMF）中，其均方根强度 $B_{rms} = 1$ nG，相干长度 $l_c = 0.5$ Mpc。基于柯尔莫哥洛夫湍流计算能量依赖的扩散系数 $D(E)$，区分共振和非共振扩散机制。

3. 通量计算与拟合：

   • 通过对红移和光度积分源发射率，并应用传播定理对均匀分布的合成源贡献求和，计算扩散通量。
   • 模型包含混合核成分（p, He, N, Si, Fe），其比例作为自由参数处理。
   • 源发射率归一化通过重子加载因子 $\kappa = L_{CR}/L_\gamma$ 与伽马射线光度相关联。
   • 参数（源分离距离 $d_s$、加速效率 $\eta$ 和核比例 $\epsilon_Z$）根据奥热和 TA 的观测数据（通量和平均簇射深度 $\langle X_{max} \rangle$）进行拟合。
   • 使用简化 $\chi^2$、赤池信息准则（AIC）和贝叶斯信息准则（BIC）进行统计模型比较。
   • 使用 HEALPix 框架生成全天空积分通量（$10-100$ EeV）和偶极各向异性图。

主要结果

1. 通量与成分：

   • FSRQ 和 BL Lac 情景均能利用混合成分重现观测到的 UHECR 能谱和 $\langle X_{max} \rangle$ 数据。
   • 最佳拟合的加速效率为 $\eta \approx 0.1$，与相对论喷流加速的数值模拟一致。
   • 所需的重子加载因子为 $\kappa \approx 5.65$，表明宇宙线功率在能量上是可行的（仅为伽马射线光度的几倍）。
   • FSRQ 情景： 需要源分离距离 $d_s \approx 8 \pm 1$ Mpc 以拟合通量。这意味着源密度显著高于 Fermi-LAT 光度函数观测估计值。
   • BL Lac 情景： 需要 $d_s \approx 160 \pm 21$ Mpc，这与 BL Lac 密度的观测估计一致。

2. 统计比较：

   • 奥热和 TA 联合数据的简化 $\chi^2$ 值分别为 2.37（FSRQ）和 2.55（BL Lac）。
   • 模型选择标准（AIC 和 BIC）倾向于 BL Lac 模型（$\Delta$AIC = 4.17, $\Delta$BIC = 4.83），提供了中等证据表明 BL Lacs 比 FSRQs 能更好地描述数据。

3. 各向异性约束：

   • BL Lacs： 使用通量拟合参数（$z_{max}=0.5$）预测的偶极各向异性与奥热观测结果吻合良好。
   • FSRQs： 使用相同参数时，预测的偶极振幅显著高于观测值。为了匹配观测到的各向异性，必须人为将红移范围限制在 $z_{max} \approx 0.01$，这与通量拟合要求以及 FSRQs 已知的强宇宙学演化相矛盾。

4. 光致裂变：

   • 光深和吸收分数的计算表明，与较轻的原子核（如 He）相比，重原子核（如 Fe）在较低能量下就会发生显著的光致裂变，导致在最高能量处重成分逐渐被抑制。

意义与结论
本文结论认为，虽然 FSRQs 和 BL Lacs 都能重现观测到的 UHECR 通量，但在考虑全套观测约束时，BL Lac 情景受到强烈支持。

• FSRQs 作为扩散 UHECRs 的主导源受到排斥，因为该情景需要不切实际的高源密度（小 $d_s$）来拟合通量，并且预测的偶极各向异性超出了观测限制，除非人为地在低红移处截断源种群。
• BL Lacs 与观测保持一致，提供了通量、成分和偶极各向异性的自洽描述，其源密度和红移演化与 Fermi-LAT 数据相符。

作者指出，虽然 FSRQs 仍可能对多信使通道（中微子和伽马射线）有所贡献，但与 BL Lac 天体相比，它们很可能是扩散 UHECR 通量的次要贡献者。本研究强调了在约束 UHECR 起源时，使用真实的光度依赖源演化模型而非简化的唯象形式的重要性。