[论文解读] Model analysis of the very high energy detections of the starburst galaxies M82 and NGC 253

E. De Cea del Pozo, Torres, D. F.+2 more|ArXiv.org|Dec 17, 2009

Astronomy and Astrophysical Research被引用 9

一句话总结

本文提出了一种多波长模型，利用改进的Q-Diffuse代码解释星暴星系M82和NGC 253中极高能（VHE）和高能（HE）伽马射线辐射。通过在现有模型中微调宇宙射线注入效率、磁场强度和扩散 timescale 等参数，作者成功再现了VERITAS和H.E.S.S.的VHE探测结果以及Fermi的HE数据，证实了宇宙射线与周围气体的强子相互作用可解释整个能谱范围内的观测伽马射线辐射。

ABSTRACT

In the light of recent results detecting the nearest starburst galaxies M 82 and NGC 253 at very high energy (VHE), we present a multi-wavelength model succesfully explaining their gamma-ray diffuse emission. The detection of M 82 was presented by VERITAS collaboration in the recent ICRC, while the integral flux of the fainter NGC253 has just been published by HESS Collaboration. Also, in the 2009 Fermi Symposium itself, Fermi data coming from the direction of both galaxies was released, thus confirming their detection in the high energy (HE) regime. Slight and reasonable variations in the space parameter of already published models can fully account for the HE and VHE emission coming from M 82 and NGC 253, while agreeing with previous data detected from radio to infrared. We explore these changes and the implications they have for the cosmic-ray distribution in these galaxies.

研究动机与目标

• 在一致的多波长框架内解释VERITAS和H.E.S.S.对M82和NGC 253的近期VHE伽马射线探测结果。
• 将观测到的VHE和HE伽马射线通量与星暴星系中宇宙射线传播和强子相互作用的现有模型相协调。
• 利用Fermi、VERITAS和H.E.S.S.的更新观测数据，约束M82和NGC 253中的宇宙射线分布与能量密度。
• 探讨关键天体物理参数（如超新星爆发率、宇宙射线注入效率和磁场）的微小变化如何再现从射电到VHE伽马射线的全谱能分布。
• 证明单一宇宙射线谱可解释所有观测到的辐射成分：π⁰衰变、轫致辐射和逆康普顿辐射。

提出的方法

• 使用Q-Diffuse代码求解质子和电子的稳态扩散-损失方程，包含能量损失、次级粒子产生和强子相互作用。
• 模型计算中性π⁰衰变、逆康普顿散射和轫致辐射产生的伽马射线，以CMB和红外光子场作为靶场。
• 通过组合同步辐射、自由-自由辐射和尘埃发射组分（冷尘和热尘分别对应45 K和200 K）拟合射电和红外数据。
• 在物理解释范围内调整磁场强度和扩散系数，以匹配VHE和HE通量，M82中为170–210 μG，NGC 253中为200–270 μG。
• 扩散 timescale 在1至10 Myr之间变化，NGC 253中质子注入谱指数调整为2.1，原始模型中为2.3。
• 模型引入宇宙射线以300–600 km/s的速度对流逃逸，并考虑能量依赖的逃逸 timescale，尤其针对E > 1500 GeV的粒子。

实验结果

研究问题

• RQ1 对现有宇宙射线模型进行微小修改，能否解释VERITAS和H.E.S.S.对M82和NGC 253的VHE伽马射线探测？
• RQ2 为再现观测到的VHE和HE伽马射线通量，所需的宇宙射线注入效率、磁场和扩散 timescale 值是多少？
• RQ3 星暴星系中宇宙射线的能量密度与银河系相比如何，特别是在1500 GeV以上能量段？
• RQ4 宇宙射线与周围气体的强子相互作用在多大程度上解释了全谱能分布范围内的观测伽马射线辐射？
• RQ5 距离、气体密度和恒星形成率的不确定性在多大程度上影响模型与多波长观测的一致性？

主要发现

• 通过仅对先前发表参数的微小调整，模型成功再现了VERITAS探测到的M82的VHE伽马射线通量和H.E.S.S.探测到的NGC 253的VHE伽马射线通量。
• 在VHE能量段，主导的伽马射线辐射机制是中性π⁰衰变，而逆康普顿散射和轫致辐射在较低能量段有显著贡献。
• M82和NGC 253中宇宙射线能量密度在1500 GeV以上达到约10 eV cm⁻³，表明其宇宙射线种群相比银河系显著增强。
• M82中所需的磁场强度为170–210 μG，NGC 253中为200–270 μG，与星暴环境的预期一致。
• 两星系中的扩散系数在1–10 GeV能量段估计为~10²⁶–10²⁷ cm² s⁻¹，显著低于银河系中的值（~10²⁸ cm² s⁻¹）。
• NGC 253中质子与电子的初始比例降低至30，且注入谱指数被展平至2.1，以匹配H.E.S.S.的积分通量。