[论文解读] On the hadronic origin of the very high energy $γ$-ray emission surrounding the young massive stellar cluster Westerlund 1

Zhaodong Shi, Rui-zhi Yang|arXiv (Cornell University)|Jan 24, 2026

Astrophysics and Cosmic Phenomena被引用 0

一句话总结

该论文提出了一个 CR 加速的 hadronic 模型，研究 Westerlund 1 超泡的终止冲击处的粒子加速，并显示若发射源为 hadronic，尤其是在冲击风腔内的团聚气体存在时，可以解释来自 HESS J1646-458 的观测到的超高能伽马射线。

ABSTRACT

The Westerlund 1 (Wd 1) is the most massive known young star cluster in the Galaxy, and an extended $γ$-ray source HESS J1646-458 surrounding it has been detected up to 80 TeV in the very high energy, implying that cosmic rays (CRs) are accelerated effectively in the region. However, the dominant radiation process contributing to the $γ$-ray emission is not well constrained. In the present work, we develop a model of CR acceleration at the termination shock in the superbubble inflated by the interaction of the cluster wind from the Wd 1 with the surrounding interstellar medium. We then calculate the flux and radial profile of $γ$ rays produced by the inelastic collisions of the hadronic CRs with the ambient gas. Our results with reasonable parameters can explain well the spectrum and radial profile of the $γ$-ray emission of HESS J1646-458, and consequently the $γ$-ray emission of HESS J1646-458 is likely to be of hadronic origin.

研究动机与目标

• 激发在年轻大质量星团中寻找能够超越超新星残骸（SNR）加速粒子的 CR 来源。
• 为 Westerlund 1 星团风的终止冲击处的 CR 加速建立物理模型。
• 计算 hadronic CR 与环境气体相互作用产生的伽马射线光谱和径向分布。
• 评估 hadronic 场景是否能够再现观测到的 HESS J1646-458 光谱和形态。

提出的方法

• 建模 Westerlund 1 形成的超泡结构并推导 TS 半径 R_TS 和壳层半径 R_S。
• 求解带扩散系数 D(p) 与对流速度 u(r) 的稳态传输方程，其中包含在 TS 的源项以及 Rankine-Hugoniot 边界条件。
• 假设均匀扩散系数 D(p)=D0β(p/1 GeV/c)^δ 并测试 Kolmogorov (δ=1/3) 与 IK (δ=1/2) 湍动。
• 利用 aafragpy 的微分截面，计算 pp 碰撞的伽玛射线发现率 q_pH→γ(E_γ) 并积分得到通量与强度。
• 在冲击风腔内的 Canonical 气体密度（式 (Eq. 6)）与反映可能的团聚气体含量的自由变化 ρ2 之间进行探索。
• 在拟合 HESS J1646-458 数据时，改变 D0、ρ0 与 η_inj，以约束 ξ_CR 与其他参数，从而拟合伽马射线通量和径向轮廓。

实验结果

研究问题

• RQ1 CR 是否能在 Westerlund 1 的终止冲击处被加速，以解释观测到的极高能伽马射线？
• RQ2 hadronic 过程是否足以解释 HESS J1646-458 的光谱和径向分布，还是偏向 leptonic 来源？
• RQ3 扩散、气体密度与注入效率如何影响 Westerlund 1 周围的伽马射线通量与空间分布？
• RQ4 在 hadronic 假设下，可再现观测的 ξ_CR 与扩散系数的合理取值范围有什么？

主要发现

• hadronic 模型在合理参数下可以同时解释 HESS J1646-458 的光谱和径向分布。
• 最佳拟合的扩散系数为 D0 ≈ 3.6–4.0 × 10^25 cm^2 s^-1，当 ρ0 = 30–50 m_p cm^-3 时，相应的 ρ2 ≈ 1.22–1.31 m_p cm^-3（经过微调后）。
• 在冲击风腔内具有团聚且平均气体密度较高且 ρ2 ≈ 1 m_p cm^-3 时，CR 与行星压比 ξ_CR 约为 0.1 即足以匹配伽马射线通量，避免 ξ_CR > 1。
• 模型能够再现观测到的角度范围和伽马射线的径向轮廓，支持 Westerlund 1 周围的终止冲击位点作为 CR 加速至亚 PeV 能级的场所。
• 如果 ρ2 受到光滑风泡预期（式 6）的约束，则 ξ_CR 需 >1，使得只有在存在团聚气体的情况下才使 hadronic 场景更具可行性。
• 该研究强化了年轻大质量星簇及其超泡对银河 CR 在亚 PeV 能量水平的贡献潜力。