On the Contribution of Local Sources to the Galactic Cosmic-Ray Spectrum: An Exact Series Solution for Two-Zone Diffusion

技術要約：二層拡散モデルによる銀河宇宙線スペクトルへの局所的ソースの寄与

問題提起
「ニー（knee）」以下の銀河宇宙線（CR）の陽子およびヘリウムのスペクトルの観測結果は、単純な冪乗則からの逸脱を示しており、これにはマルチTeV領域における構造が含まれる。一つの、あるいは少数の近傍ソースによる支配的な寄与が潜在的な説明となり得るが、均一拡散モデル（HDM）を用いた先行研究では、そのようなシナリオは統計的に可能性が低いと結論付けられている。具体的には、モンテカルロ・シミュレーションによれば、近傍のソースが局所的なスペクトルを著しく変化させる確率は非常に低い（約0.4%）とされる。しかし、これらの結論は、宇宙線がソースから離脱した直後に効率的に拡散するという仮定に基づいている。近年の銀河加速器（例：Cygnus X、パルサー・ヘイロー）周囲における拡張ガンマ線放射の観測は、ソース近傍での輸送が抑制され、拡散係数が減少している可能性を示唆している。本論文では、宇宙線がソース周辺の局所的な領域において非効率な輸送を経験する場合、局所的な支配的ソースが存在する統計的確率がどのように変化するかを調査する。

手法
局所的な低拡散領域を伴う宇宙線伝播のモデリングにおける計算上の課題に対処するため、著者らは球状の二層拡散問題に対する厳密な解析的級数解（グリーン関数）を導出した。

• モデルの幾何学： 本モデルは、拡散係数 $D_1$ を持つ内側領域（$0 < r < r_1$）と、より大きな係数 $D_2$ を持つ外側領域（$r_1 < r < r_2$）で構成される。界面 $r_1$ において、密度と拡散フラックスの連続性が強制される。
• 解析解： 著者らは、バースト的注入および連続注入の両方のシナリオに対して、固有モード展開を導出した。この解析的なカーネルにより、粒子分布の高速な計算が可能となり、数値コードにおいて小スケールと大スケールの空間スケールを同時に解像するという非効率性を回避しながら、大規模なモンテカルロ・シミュレーションを実行することが可能となった。
• モンテカルロ・シミュレーション： 導出されたグリーン関数を用い、Evoliらによる「Jelly」モデルを採用した $10^5$ 個の銀河超新星残骸（SNR）集団の実現例を生成した。10 TeVにおいて、最も強力な局所ソース（距離 $d < 0.5$ kpc かつ 年齢 $t < 1$ Myr と定義）のフラックス寄与を、他のすべてのソースからの背景フラックスと比較した。
• 候補解析： 本研究では、特定のカタログ化された近傍候補（パルサーおよびLoop I）についても検討し、異なる輸送仮定（HDM、二層拡散モデル［2ZDM］、空間依存伝播［SDP］）の下で、既知のソースが観測されたスペクトル特徴を再現できるかどうかを判定した。

主な貢献

1. 厳密な級数解： 局所的な輸送抑制を伴うソース集団の効率的なモンテカルロ研究を可能にする、二層拡散問題のための高速かつ解析的なグリーン関数の導出。
2. 定量的確率の変化： 近傍での輸送が抑制された場合、局所ソースが宇宙線スペクトルを支配する統計的な尤度がどのように変化するかについての厳密な再評価。
3. モデル依存性の実証： 仮定された輸送モデルおよび注入スペクトルのパラメータに応じて、特定の近傍ソース（例：Vela、Geminga、Monogem）の予測される寄与が劇的に変化することを示す系統的な比較。

結果

• 確率の増大： 標準的なHDMでは、最も強力な局所ソースが10 TeVにおいて背景と同等のフラックスに寄与する確率は約0.4%である。2ZDMの下では、この確率はバースト的注入で1.7%、連続注入で2.2%に上昇する。これは4倍から5倍の増加ではあるが、絶対的な確率は依然として低い（パーセントレベル）。
• 物理的メカニズム： 内側の低拡散領域は一時的なリザーバーとして機能し、粒子の脱出を遅らせ、時間およびエネルギーの再分配を行う。この効果は、内側領域の拡散時間（$\tau \sim 1$）に相当するタイムスケールにおいて、内側領域の外側にいる観測者（$\xi > 1$）に対して最も顕著となる。
• 候補ソース解析： パルサーをSNRのトレーサーとして用いて既知の近傍候補に2ZDMを適用した場合、10 TeVのスペクトル特徴を再現するには、以下の追加的な仮定が同時に満たされる必要がある：
  • より硬い注入スペクトル： 局所ソースは、背景集団（$\gamma \approx 2.34$）と比較して、著しく硬い注入指数（$\gamma \approx 2.08$）を持つ必要がある。もし局所的な注入指数が背景と一致する場合、その寄与は無視できるものとなる。
  • 輸送モデルへの感度： 特定のソースからの予測フラックスは、HDM、2ZDM、およびSDPモデルの間で強く変動する。例えば、SDPの下では、GemingaとMonogemの相対的な寄与は2ZDMと比較して異なっており、「ソースの重要性」は固有の性質ではなく、仮定された磁気環境や輸送特性に依存することを浮き彫りにしている。

意義と主張
本論文は、近傍での輸送抑制は、局所ソースの支配に対する統計的な反論を大幅に弱めるものの、その困難さを排除するものではないと結論付けている。局所ソースが支配的となる確率は依然として低い（約1〜2%）。したがって、マルチTeVのスペクトル構造を特定の局所ソースに由来すると解釈することは、極めてモデル依存性が高い。

著者らは、近傍のソースがスペクトル特徴を説明できるのは、有利な条件が同時に満たされる場合に限られると主張している：すなわち、特定のソースの年齢／距離、硬い局所注入スペクトル、および特定の輸送環境である。したがって、スペクトルのフィッティングのみに基づいて、あるスペクトル特徴の起源を局所ソースと特定することは堅牢ではない。本論文は、スペクトル、組成、異方性、およびガンマ線データの結合モデリングを通じて、ソースの注入履歴、粒子輸送メカニズム、および局所的な星間乱流に対する独立した制約を確立することが、決定的な局所起源の特定には必要であると論じている。