[論文レビュー] 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

ひとこと要約

この論文は、Westerlund 1 超泡の終端ショックでのハドロニックCR加速モデルを開発し、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加速の物理モデルを開発する。
• 周囲のガスとのハドロニックなCR相互作用からガンマ線スペクトルと半径方向プロファイルを計算する。
• はハドロニック仮説が観測された HESS J1646-458 のスペクトルと形態を再現できるかどうかを評価する。

提案手法

• Westerlund 1 によって形成される超泡の構造をモデル化し、TS半径 R_TS と殻半径 R_S を導出する。
• CR の定常輸送方程式を拡散 D(p) と対流 u(r) を用いて解く。ソース項を TS に含み、Rankine-Hugoniot 条件を適用する。
• 拡散係数 D(p)=D0β(p/1 GeV/c)^δ を均質と仮定し Kolmogorov (δ=1/3) および IK (δ=1/2)の乱流を試す。
• pp 衝突からの微分断面に基づいてガンマ線放出率 q_pH→γ(E_γ) を計算し、フラックスと強度を得るために積分する。
• 衝撃風腔内の canonical なガス密度 (式 6) と、クランプしたガス量を反映させるため ρ2 を自由に変化させて探索する。
• モデルのガンマ線フラックスと半径方向プロファイルを HESS J1646-458 データにフィットさせ、D0, ρ0, η_inj を変化させて ξ_CR などのパラメータを制約する。

実験結果

リサーチクエスチョン

• RQ1 CR は Westerlund 1 の終端ショックで加速され、観測された非常に高エネルギーガンマ線を説明できるか？
• RQ2 ハドロニック過程は HESS J1646-458 のスペクトルと半径方向プロファイルを説明するのに十分か、それとも電子起源が有利か？
• RQ3 拡散、ガス密度、注入効率は Westerlund 1 の周囲のガンマ線フラックスと空間分布にどう影響するか？
• RQ4 観測を再現するために、Hadronic 仮定の下での CR-to-ram pressure 比 ξ_CR および拡散係数の妥当な値は何か？

主な発見

• ハドロニックモデルは reasonable parameters の下でスペクトルと半径プロファイルの両方を説明できる。
• Best-fit diffusion coefficients are D0 ≈ 3.6–4.0 × 10^25 cm^2 s^-1 for ρ0 = 30–50 m_p cm^-3, with corresponding ρ2 ≈ 1.22–1.31 m_p cm^-3 when fine-tuning is applied.
• With a clumpy, higher average gas density in the shocked wind cavity (ρ2 ≈ 1 m_p cm^-3), a CR-to-ram pressure ratio ξ_CR ≈ 0.1 suffices to match the gamma-ray flux, avoiding ξ_CR > 1.
• The model reproduces the observed angular extent and radial gamma-ray profile, supporting the termination shock around Westerlund 1 as a site of CR acceleration up to sub-PeV energies.
• If ρ2 is constrained by a smooth wind-bubble expectation (Eq. 6), ξ_CR would need to be >1, making the hadronic scenario less favorable unless clumpy gas is present.
• The study reinforces that young massive clusters and their superbubbles can contribute meaningfully to Galactic CRs up to sub-PeV energies