Modelling Galactic neutrino emission: contributions from massive star clusters and interstellar cosmic rays

技術要約：大質量星団および銀河間宇宙線からの銀河ニュートリノ放射のモデリング

問題提起
宇宙における最も高エネルギーな宇宙線（特にPeVエネルギーに達するもの、すなわちPeVatron）の起源は、依然として不透明である。超新星残骸（SNR）は伝統的に主要な加速器と考えられてきたが、$\gamma$線観測によれば、それらはPeVエネルギーに達するには不十分である可能性が示唆されており、代替となるソースクラスの必要性が高まっている。若い大質量星団（YMSC）は、衝突する星風や複数のSNRによって生成される強化された磁気乱流により、有望なPeVatron候補として提案されている。同時に、IceCube観測所は銀河面からの高エネルギーニュートリノフラックスを検出している。しかし、現在の解析は、宇宙線の海（CR sea）と銀河間ガスとの相互作用によって生成される純粋に拡散的なニュートリノ放射のみを考慮したテンプレートに依存しており、これらのモデルは、YMSCのような未分解のハドロン加速器による寄与を明示的に解明できておらず、これらソース固有のシグネチャーを覆い隠している可能性がある。

手法
著者らは、以下の2つの異なる成分を統合して、銀河からの期待ニュートリノ放射を計算するための新しい理論的枠組みを提示する。

1. 合成YMSC集団： 太陽近傍（$\lesssim 2$ kpc）のYMSCのみを観測できるという制限を克服するため、著者らは銀河のYMSC集団に対して100個の独立した実現化（realizations）を生成した。この合成集団は、観測的制約に基づき、星団質量（$10^3 - 6.3 \times 10^4 M_\odot$）、年齢（$<30$ Myr）、および空間分布（螺旋腕、銀河棒、およびローカルスパーに従う）をサンプリングしている。
2. 粒子加速と輸送： 各合成星団について、著者らは恒星集団（主系列星、ウォルフ・ライエ星、および超新年代）をシミュレートし、風と超新星からの集団的な機械的光度を計算する。その後、風によって吹き飛ばされたバブル内のガス密度と磁場プロファイルをモデル化する。
   • 加速メカニズム： 粒子加速は、集団的な風の終端衝撃波（WTS）と超新星残骸（SNR）の衝撃波という2つのサイトでモデル化される。
   • 拡散シナリオ： 星団バブル内での宇宙線の伝播は、3つの異なる乱流仮定（コルモゴロフ、クライナックス、およびボーム拡散）の下で計算される。
   • ニュートリノ生成： 加速された陽子が標的ガス（主にバブルの密な掃き寄せられたシェル内）と衝突するハドロン衝突から生じるニュートリノフラックスを計算する。
3. 拡散銀河成分： モデルは、銀河宇宙線の海と銀河間ガスとの相互作用から生じる拡散的なニュートリノ放射も計算する。これには、宇宙線伝播モデル（銀河中心部におけるスペクトル硬化を組み込んだもの）とガスマップ（GALPROPおよび代替の3Dモデル）を用いる。
4. テンプレート生成： 予測される全フラックスは、YMSCの寄与とCR誘起の拡散成分の和である。これらは、IceCubeのデータと比較するために、空間的およびスペクトル的なテンプレートへと合成される。

主な貢献

• 新規ニュートリノテンプレート： 本論文は、標準的なCR誘起の拡散成分とともに、YMSC集団の放射率を明示的に組み込んだ、初の包括的なニュートリノテンプレートを提供する。
• 合成集団モデリング： 近傍の解明された少数の星団のみに依存するのではなく、銀河全体のYMSC集団からの累積的なニュートリノ出力を統計的に堅牢な推定手法で提供する。
• マルチコンポーネント解析： 星団内での「in-situ」のニュートリノ生成と、背景となる拡散放射を区別することで、ソースの寄与の分離を可能にする。
• パラメータ感度： 粒子拡散レジーム（コルモゴロフ、クライナックス、ボーム）および銀河中心部におけるCRスペクトルの硬化への依存性を調査している。

結果

• 個々のソースの検出可能性： 個々のYMSCからの予測フラックスは、立方キロメートル規模の検出器（例：IceCube, KM3NeT）であっても検出限界に近いことが示された。したがって、YMSC集団からの累積的な放射は、未分解の拡散成分として現れることが予想される。
• フラックスの大きさ： モデルの正規化は、既存の銀河テンプレートを用いたIceCubeのベストフィットと一致している。著者らは、クラスターからの放射による未分解の寄与が無視できない可能性があることを見出した。
• 拡散への依存性：
  • ボーム拡散シナリオ（最も効率的な閉じ込め）では、YMSCは$\sim 50$ TeV以下でニュートリノフラックスを支配し、スペクトル硬化を仮定しない場合、KraichnanおよびKolmogorovシナリオと比較してIceCubeのデータに対して過大予測する可能性がある。
  • コルモゴロフおよびクライナックスシナリオでは、一般にCR誘起の拡散成分がYMSC成分と同等、あるいはそれを上回る寄与を示す。
• 空間分布： YMSCの寄与は、内銀河経度範囲（$|l| \le 30^\circ$、銀河リッジ）に集中している一方、CR誘起の放射はより大きな経度において支配的となる。
• スペクトル硬化： 銀河中心部におけるCRのスペクトル硬化を含めることで、コルモゴロフおよびクライナックスシナリオにおけるYMSCの相対的な寄与が減少し、全予測スペクトルが$\sim 100$ TeVにおけるIceCubeの観測結果とより良く一致するようになる。
• 整合性： ハドロン衝突の$\gamma$線カウンターパートのモデル予測は、現在のTeV-PeV領域における拡散$\gamma$線観測と矛盾しない。

意義と主張
著者らは、自らの研究が銀河ニュートリノ放射の理論的モデリングにおける必要な進化を提供すると主張している。純粋に拡散的なテンプレートを超えて、YMSC集団の特定の放射率を含めることで、銀河ニュートリノ放射の2つの成分を「分離」するためのツールを提供している。本論文は、これらの新しいテンプレートを将来のニュートリノデータセットに適用することで、極限的な宇宙線加速とニュートリノ生成における星団の役割を制約できることを示唆している。著者らは、クラスターからの未分解の寄与はおそらく無視できないものであることを示唆しており、ソースの寄与と拡散成分の同時フィッティングが、銀河加速器の特定と宇宙線の起源の解決に不可リストであると強調している。本研究は、IceCubeの信号源としてYMSCを決定的に特定したと主張するものではなく、この仮説を観測データに対して検証するために必要な理論的枠組みを提供するものである。