Young massive stellar clusters (YMSCs) have recently regained interest as PeVatron candidates, potentially accounting for the cosmic-ray (CR) knee as alternatives to isolated supernova remnants (SNRs). LHAASO's unique capability to detect photons above 0.1 PeV, hence tracing multi-PeV CRs, can provide critical constraints on galactic acceleration models when combined with H.E.S.S. and Fermi-LAT data. We investigate the transport of particles from YMSCs acceleration sites, namely wind termination shocks (WTS) or embedded SNRs, to nearby dense molecular clouds where proton-proton interactions produce high-energy gamma rays. We determine the necessary conditions, such as the distance between the acceleration site and the target, or the cluster's power and age, for detectable gamma-ray excesses and identify viable systems through parameter space exploration. By comparing with observations, we can constrain key physical parameters including WTS efficiency, diffusion coefficient and injection slope. Our analysis also examines whether some of LHAASO's unidentified sources might correspond to such cluster-cloud systems.
Modeling gamma-ray signatures of particle acceleration in stellar clusters from GeV to PeV 论文ID : 2510.12562标题 : Modeling gamma-ray signatures of particle acceleration in stellar clusters from GeV to PeV作者 : A. Inventar, S. Gabici, E. Peretti分类 : astro-ph.HE (高能天体物理)发表时间/会议 : 39th International Cosmic Ray Conference (ICRC2025), 2025年7月15-24日, 日内瓦论文链接 : https://arxiv.org/abs/2510.12562 年轻大质量恒星团(YMSCs)作为PeVatron候选体重新引起关注,可能作为孤立超新星遗迹(SNRs)的替代解释宇宙线"膝部"现象。LHAASO独特的0.1 PeV以上光子探测能力,能够追踪多PeV宇宙线,结合H.E.S.S.和Fermi-LAT数据可为银河系加速模型提供关键约束。本文研究粒子从YMSCs加速区(恒星风终止激波WTS或嵌入式SNRs)传输到邻近致密分子云的过程,通过质子-质子相互作用产生高能伽马射线。通过参数空间探索确定可探测伽马射线超出所需的必要条件,如加速区与靶区距离、星团功率和年龄等,并识别可行系统。通过与观测比较,约束WTS效率、扩散系数和注入谱指数等关键物理参数。分析还探讨LHAASO未识别源是否对应此类星团-云系统。
宇宙线加速机制 :传统认为超新星遗迹是银河系宇宙线的主要加速源,但在解释PeV能段("膝部"区域)宇宙线时存在困难PeVatron候选体 :需要寻找能够加速粒子到PeV能量的天体物理源观测约束不足 :现有伽马射线观测数据对高能宇宙线加速模型的约束有限能量尺度关键性 :PeV能段是银河系与河外宇宙线成分的分界,理解该能段加速机制对宇宙线起源至关重要多信使天体物理 :伽马射线观测提供宇宙线加速的直接证据新一代探测器 :LHAASO等新设备的超高能探测能力为理论模型提供前所未有的验证机会孤立SNR模型难以解释观测到的超高能伽马射线发射的空间延展性 缺乏系统性的理论框架连接恒星团、宇宙线传输和分子云相互作用 参数空间探索不充分,难以确定可探测信号的条件 建立了完整的理论框架 :从恒星团内宇宙线加速到分子云伽马射线产生的端到端物理模型提供解析解 :针对脉冲式和连续注入两种情况,给出宇宙线传输的解析表达式确定探测条件 :通过参数空间分析,确定产生可探测伽马射线超出的必要条件实际应用验证 :以W43星形成区为例,约束了关键物理参数多能段覆盖 :建立从GeV到PeV的统一模型框架输入 :恒星团物理参数(年龄、质量、功率),宇宙线注入谱,扩散参数,分子云性质
输出 :预期伽马射线流量谱,可探测性评估
约束条件 :恒星泡半径限制,扩散时标与相互作用时标关系
恒星泡半径由Weaver等人的解析解给出:
R b = 260 ( η N ∗ 100 ) 1 / 5 ( n 0 cm − 3 ) − 1 / 5 ( t 10 Myr ) 3 / 5 pc R_b = 260\left(\frac{\eta N_*}{100}\right)^{1/5} \left(\frac{n_0}{\text{cm}^{-3}}\right)^{-1/5} \left(\frac{t}{10 \text{Myr}}\right)^{3/5} \text{pc} R b = 260 ( 100 η N ∗ ) 1/5 ( cm − 3 n 0 ) − 1/5 ( 10 Myr t ) 3/5 pc
其中N ∗ N_* N ∗ 为大质量恒星数目,n 0 n_0 n 0 为星际介质密度,η \eta η 为辐射损失修正因子。
脉冲注入情况 :
宇宙线数密度分布为:
n ( t , E , R ) = Q ( E ) π 3 / 2 R d 3 exp [ − ( R R d ) 2 ] n(t,E,R) = \frac{Q(E)}{\pi^{3/2}R_d^3} \exp\left[-\left(\frac{R}{R_d}\right)^2\right] n ( t , E , R ) = π 3/2 R d 3 Q ( E ) exp [ − ( R d R ) 2 ]
其中R d ( t , E ) = 4 D ( E ) t R_d(t,E) = \sqrt{4D(E)t} R d ( t , E ) = 4 D ( E ) t 为扩散半径。
连续注入情况 :
n ( t , E , R ) = Q ˙ ( E ) 4 π D ( E ) R erfc ( R R d ( t , E ) ) n(t,E,R) = \frac{\dot{Q}(E)}{4\pi D(E)R} \text{erfc}\left(\frac{R}{R_d(t,E)}\right) n ( t , E , R ) = 4 π D ( E ) R Q ˙ ( E ) erfc ( R d ( t , E ) R )
通过质子-质子相互作用产生的伽马射线流量:
ϕ γ = ∫ d σ d E γ ( E , E γ ) n ( E ) d E c 4 π d 2 M cloud m p \phi_\gamma = \int \frac{d\sigma}{dE_\gamma}(E,E_\gamma) n(E) dE \frac{c}{4\pi d^2} \frac{M_{\text{cloud}}}{m_p} ϕ γ = ∫ d E γ d σ ( E , E γ ) n ( E ) d E 4 π d 2 c m p M cloud
统一框架 :首次将恒星风终止激波(WTS)和嵌入式SNR两种加速机制纳入统一分析框架解析方法 :提供简化但准确的解析表达式,便于参数空间快速探索多尺度建模 :从pc尺度的恒星泡到kpc尺度的宇宙线传输观测约束集成 :结合多个伽马射线探测器的敏感度进行可探测性分析LHAASO数据 :超高能伽马射线观测(>0.1 PeV)H.E.S.S.数据 :甚高能伽马射线观测Fermi-LAT数据 :GeV能段伽马射线观测目标源 :W43星形成区及其相关伽马射线源LHAASO J1848-0153u扩散系数 :D ( E ) = D 10 [ E / ( 10 m p c 2 ) ] δ D(E) = D_{10}[E/(10m_pc^2)]^\delta D ( E ) = D 10 [ E / ( 10 m p c 2 ) ] δ ,D 10 ∼ 10 27 − 10 28 D_{10} \sim 10^{27}-10^{28} D 10 ∼ 1 0 27 − 1 0 28 cm²/s注入谱指数 :α = 2.0 − 2.2 \alpha = 2.0-2.2 α = 2.0 − 2.2 扩散指数 :δ = 0.3 − 0.5 \delta = 0.3-0.5 δ = 0.3 − 0.5 最大能量 :E max = 3 E_{\max} = 3 E m a x = 3 PeV恒星团功率 :W ˙ C R = 3 × 10 37 \dot{W}_{CR} = 3 \times 10^{37} W ˙ CR = 3 × 1 0 37 erg/sSNR能量 :W C R = 10 50 W_{CR} = 10^{50} W CR = 1 0 50 erg分子云质量 :10 5 − 10 6 M ⊙ 10^5-10^6 M_\odot 1 0 5 − 1 0 6 M ⊙ 距离范围 :50-300 pc脉冲注入 :
R d ( t max ) ≈ 102 ( α − 2 ) 1 / 3 ( W C R 10 50 erg ) 1 / 3 ( E m p c 2 ) 0.9 − 1 3 α pc R_d(t_{\max}) \approx 102 (\alpha-2)^{1/3} \left(\frac{W_{CR}}{10^{50}\text{erg}}\right)^{1/3} \left(\frac{E}{m_pc^2}\right)^{0.9-\frac{1}{3}\alpha} \text{pc} R d ( t m a x ) ≈ 102 ( α − 2 ) 1/3 ( 1 0 50 erg W CR ) 1/3 ( m p c 2 E ) 0.9 − 3 1 α pc
连续注入 :
R max ≈ 3 × 10 2 ( α − 2 ) ( W ˙ C R 10 37 erg/s ) ( E m p c 2 ) 2.7 − α − δ pc R_{\max} \approx 3 \times 10^2 (\alpha-2) \left(\frac{\dot{W}_{CR}}{10^{37}\text{erg/s}}\right) \left(\frac{E}{m_pc^2}\right)^{2.7-\alpha-\delta} \text{pc} R m a x ≈ 3 × 1 0 2 ( α − 2 ) ( 1 0 37 erg/s W ˙ CR ) ( m p c 2 E ) 2.7 − α − δ pc
最佳拟合参数 :α = 2.0 \alpha = 2.0 α = 2.0 , δ = 0.33 \delta = 0.33 δ = 0.33 , D 10 = 10 27 D_{10} = 10^{27} D 10 = 1 0 27 cm²/s截断能量 :E cut ∼ 0.4 E_{\text{cut}} \sim 0.4 E cut ∼ 0.4 PeVWTS效率约束 :ϵ 10 28 cm 2 / s D 10 ∼ 1 \epsilon \frac{10^{28}\text{cm}^2/\text{s}}{D_{10}} \sim 1 ϵ D 10 1 0 28 cm 2 / s ∼ 1 在300 TeV能段,需要硬注入谱(α ≲ 2.1 \alpha \lesssim 2.1 α ≲ 2.1 )和抑制扩散(D 10 ≲ 10 27 D_{10} \lesssim 10^{27} D 10 ≲ 1 0 27 cm²/s)才能产生可探测信号 WTS模型比单个SNR模型在更大参数空间内可探测 有效最大距离约为300 pc(最优参数下) GeV-TeV能段:主要受宇宙线海背景限制 1 TeV:开始出现探测器敏感度限制的有效最大距离
100 TeV:需要非常特殊的参数组合才可探测
SNR范式 :传统认为超新星遗迹是主要加速源,但在PeV能段面临挑战恒星团模型 :近年来作为PeVatron候选体重新受到关注扩散激波加速 :理论基础来自Fermi加速机制点源与延展源 :LHAASO发现的超高能源多为延展源,支持恒星团起源多信使关联 :伽马射线、宇宙线、中微子的联合分析未识别源问题 :大量伽马射线源缺乏明确对应体相比现有工作,本文首次提供了从加速到探测的完整解析框架,能够快速评估不同参数下的可探测性,为观测规划和数据解释提供理论指导。
WTS机制优势 :恒星风终止激波相比嵌入式SNR在产生可探测伽马射线方面具有更大的参数空间参数约束 :通过W43系统的拟合,约束了注入谱指数(α ∼ 2.0 \alpha \sim 2.0 α ∼ 2.0 )和扩散系数探测条件 :确定了产生可探测超高能伽马射线的必要条件,包括硬谱注入和抑制扩散距离限制 :有效探测距离约为几百pc,限制了可观测的星团-云系统数量简化假设 :采用球对称扩散模型,忽略了磁场结构的复杂性单一目标 :仅考虑单个分子云目标,实际情况可能涉及多个云团参数退化 :某些参数组合存在退化,难以独立约束统计样本 :仅分析了W43一个系统,需要更多案例验证3维MHD模拟 :考虑磁场结构和湍流的复杂影响多目标建模 :同时考虑多个分子云的贡献系统性巡天 :分析更多LHAASO未识别源多信使分析 :结合中微子观测约束强子成分理论完整性 :建立了从微观物理到宏观观测的完整理论链条方法创新性 :提供简洁有效的解析方法,便于参数空间探索实际应用价值 :直接应用于LHAASO观测数据,具有很强的现实意义多尺度统一 :成功连接了不同空间和时间尺度的物理过程模型简化过度 :球对称假设可能过于简单,实际系统具有复杂的3维结构参数退化严重 :多个参数的组合效应使得独立约束困难统计样本不足 :仅基于单一系统的分析,结论的普适性有待验证不确定性分析缺失 :缺乏对模型参数不确定性的系统性分析理论贡献 :为恒星团宇宙线加速理论提供了重要的定量分析工具观测指导 :为LHAASO等新一代探测器的数据解释提供理论框架学科交叉 :连接了恒星形成、宇宙线物理和伽马射线天文学实用价值 :可直接用于评估其他星团-云系统的可探测性LHAASO数据分析 :特别适用于超高能延展伽马射线源的理论解释CTA观测规划 :可用于预测和规划未来观测目标宇宙线起源研究 :为理解银河系宇宙线"膝部"现象提供新视角恒星团研究 :为年轻大质量恒星团的多波段研究提供理论支持论文引用了该领域的关键文献,包括:
Gabici, S. 2023 关于恒星团宇宙线加速的综述 LHAASO Collaboration 2024 关于超高能伽马射线观测 Aharonian & Atoyan 1996 宇宙线传输的经典理论 Bell, A. R. 2004 关于流不稳定性对扩散的影响 该论文在理论建模和观测应用之间建立了很好的桥梁,为理解超高能宇宙线起源提供了新的视角。虽然存在一些简化假设的局限性,但其解析方法的简洁性和实际应用价值使其成为该领域的重要贡献。