[论文解读] Extragalactic Propagation of Ultra-High Energy Cosmic Rays

D. Kuempel|arXiv (Cornell University)|Sep 10, 2014

Astrophysics and Cosmic Phenomena参考文献 5被引用 3

一句话总结

本文综述了超高能宇宙射线（UHECRs）在星系际空间中的传播，重点研究其与宇宙微波背景光子的相互作用以及在星系际磁场中的偏转。利用先进的模拟代码，研究表明，Pierre Auger 和 Telescope Array 的 GZK 抑制效应与成分数据限制了源谱和化学成分，支持硬注入指数并将质子分数限制在约 10%，暗示需要多信使分析以阐明其天体物理起源。

ABSTRACT

More than 100 years after the discovery of cosmic rays and various experimental efforts, the origin of ultra-high energy cosmic rays (E > 100 PeV) remains unclear. The understanding of production and propagation effects of these highest energetic particles in the universe is one of the most intense research fields of high-energy astrophysics. With the advent of advanced simulation engines developed during the last couple of years, and the increase of experimental data, we are now in a unique position to model source and propagation parameters in an unprecedented precision and compare it to measured data from large scale observatories. In this paper we revisit the most important propagation effects of cosmic rays through photon backgrounds and magnetic fields and introduce recent developments of propagation codes. Finally, by comparing the results to experimental data, possible implications on astrophysical parameters are given.

研究动机与目标

• 理解超高能宇宙射线（UHECRs）如何通过宇宙光子背景和磁场传播，影响其能谱和到达方向。
• 评估光致π介子产生及能量损失机制（如 GZK 效应）对 UHECR 通量抑制的影响。
• 评估磁场偏转在改变 UHECR 天空分布及影响成分解释中的作用。
• 将模拟的 UHECR 传播与 Pierre Auger 和 Telescope Array 等大型观测站的数据进行比较，以约束天体物理模型。
• 倡导采用多信使方法，结合能谱、成分、到达方向以及次级信使（中微子、光子）以改善模型约束。

提出的方法

• 使用相互作用长度公式 λ⁻¹(E) = ∫ n(ε)σ_avg(ε) dε 来模拟光子背景引起能量损失，其中 n(ε) 来自 CMB、IRB 和 URB 谱。
• 应用光致π介子产生的阈值能量：E_thres^N,π ≈ 6.8×10¹⁹ eV × (10⁻³ eV / ε)⁻¹，考虑原子核静止参考系中洛伦兹增强的光子能量。
• 通过假设核子相互作用中π介子产生与质子相互作用具有相似动力学，采用核子数比例近似来建模核子相互作用。
• 在计算机集群上进行三维传播模拟，以模拟宇宙射线在磁场和光子背景中的输运。
• 将模拟的能谱、X_max（簇射最大深度）和 σ(X_max)（波动）与 Pierre Auger 观测站的数据进行比较，使用 DPMJET 和 QGSP 等相互作用模型。
• 利用事例形状可观测量和能量-能量相关性，探测湍流和相干磁场，并将次级粒子产生（光子、中微子）与观测限制进行交叉验证。

实验结果

研究问题

• RQ1 宇宙微波背景相互作用对超高能宇宙射线的能谱和成分有何影响？
• RQ2 星系际磁场如何影响 UHECR 的到达方向分布和偏转角？
• RQ3 观测到的 GZK 抑制效应和谱形在多大程度上能约束本征源谱指数和最大能量？
• RQ4 在假设不同源成分和注入模型的前提下，模拟的 X_max 和 σ(X_max) 值能否再现空气簇射数据中的观测结果？
• RQ5 到达方向模式和次级信使（中微子、光子）如何在仅依赖能谱和成分的基础上，进一步改善对 UHECR 天体物理的约束？

主要发现

• 在 ~5×10¹⁸ eV 以上的观测到的 UHECR 能谱需要一个硬的源谱指数 β ≈ -1 才能被再现，表明存在非热、高能加速机制。
• Pierre Auger 数据表明，在 E > 10¹⁸ eV 时，质子组分分数超过 10% 的纯质子模型被各向异性上限所排除，暗示重核有显著贡献。
• 模拟表明，UHECR 能谱的低能部分（低于 ~5×10¹⁸ eV）无法通过单一组分解释，除非引入第二组分，如银河系宇宙射线或轻元素星系际源。
• 只有在使用特定相互作用模型并假设混合成分时，模拟的 X_max 及其弥散 σ(X_max) 才能与 Auger 数据良好吻合，其中铁族核在高能区比质子提供更优拟合。
• 事例形状可观测量和到达方向中的能量-能量相关性对相干和湍流磁场敏感，可实现对磁场强度和结构的间接测量。
• 模拟预测的次级中微子和光子通量与当前观测限制一致，支持 UHECR 的强子起源，并为未来多信使约束提供了可能。