[论文解读] HERMES: a Monte Carlo Code for the Propagation of Ultra-High Energy Nuclei

Manlio De Domenico, Haris Lyberis+1 more|arXiv (Cornell University)|Jul 16, 2013

Particle physics theoretical and experimental studies参考文献 14被引用 3

一句话总结

HERMES 是一种新型蒙特卡洛代码，用于模拟超高能（UHE）核素（10^16–10^22 eV）在河外和银河系磁场中的传播，并考虑其与河外背景辐射的相互作用。该代码能够实现通量、GZK 视界和到达方向分布的逼真建模，计划发布稳定公共版本，供社区用于质量组成和能谱分析。

ABSTRACT

Although the recent experimental efforts to improve the observation of Ultra-High Energy Cosmic Rays (UHECRs) above $10^{18}$ eV, the origin and the composition of such particles is still unknown. In this work, we present the novel Monte Carlo code (HERMES) simulating the propagation of UHE nuclei, in the energy range between $10^{16}$ and $10^{22}$ eV, accounting for propagation in the intervening extragalactic and Galactic magnetic fields and nuclear interactions with relic photons of the extragalactic background radiation. In order to show the potential applications of HERMES for astroparticle studies, we estimate the expected flux of UHE nuclei in different astrophysical scenarios, the GZK horizons and we show the expected arrival direction distributions in the presence of turbulent extragalactic magnetic fields. A stable version of HERMES will be released in the next future for public use together with libraries of already propagated nuclei to allow the community to perform mass composition and energy spectrum analysis with our simulator.

研究动机与目标

• 开发一种灵活、模块化的蒙特卡洛代码，用于模拟超高能（UHE）核素在磁化宇宙环境中的传播。
• 模拟河外和银河系磁场（包括规则分量与湍流分量）对 UHECR 轨迹的影响。
• 纳入与河外背景辐射（CMB、CIOB、URB）中遗迹光子的核相互作用，以及次级粒子的产生。
• 在不同天体物理情景下，估算通量、GZK 视界和到达方向分布等可观测量。
• 为社区提供稳定、公开的 HERMES 发布版本，并附带预处理的核素库，用于 UHECR 质量组成与能谱分析。

提出的方法

• HERMES 在 ΛCDM 宇宙学框架下采用弗里德曼-勒梅特-罗伯逊-沃尔克（Friedmann-Robertson-Walker）度规，支持可调的宇宙学参数。
• 通过可自定义的光度、谱指数（γ）和红移演化（如 SFR、GRB、AGN、QSO 或 (1+z)^m 形式）实现源注入模型。
• 代码纳入了 2.725 K 的宇宙微波背景（CMB），通过 Finke 等人的模型模拟红外/光学背景（CIOB），并采用红移演化模型处理通用射电背景（URB）。
• 通过数值积分求解相对论性运动方程计算带电粒子轨迹，对轻粒子包含同步辐射能量损失。
• 湍流磁场通过基于 Giacalone & Jokipii（1999）的局部逐步方法模拟，采用类似柯尔莫哥洛夫的统计方法确定均方根强度与相干长度。
• 核相互作用通过在遗迹光子上的光致分解与光致吸收截面建模，包含次级粒子产生。

实验结果

研究问题

• RQ1 在不同的天体物理源分布与演化模型下，地球处预期的 UHE 核素通量是多少？
• RQ2 不同核素与源红移下，GZK 视界如何变化？
• RQ3 在存在湍流河外磁场的条件下，附近源产生的 UHECR 到达方向分布有何特征？
• RQ4 各向同性背景事件与源星表（如 2MRS、SWIFT-BAT）如何影响模拟的 UHECR 天图？
• RQ5 当与观测曝光函数结合时，HERMES 在多大程度上能再现观测到的 UHECR 能谱与各向异性？

主要发现

• HERMES 成功模拟了 UHE 质子的存活概率及其对应的 GZK 视界，结果表明其依赖于核素质量与源红移。
• 该代码在不同天体物理情景下（包括源演化与不同注入指数）重现了预期的 UHECR 通量与能谱，结果与 HiRes 观测结果相当。
• 来自附近源（4–200 Mpc）的 UHE 质子模拟天图显示，由于河外磁场作用，其轨迹发生显著偏转，2MRS 与 SWIFT-BAT 源星表表现出不同的特征模式。
• 引入基于 Pierre Auger 测量结果的 56% 各向同性事件，改变了到达方向分布，证明了代码对背景污染的敏感性。
• HERMES 能够实现到达方向各向异性的逼真建模，且在天图中正确考虑了 Pierre Auger 观测站的非均匀曝光效应。
• 该代码具有模块化与可扩展性，可根据更新的数据，轻松集成新的背景辐射或核相互作用模型。