Neutrino Constraints on memory-burdened Primordial Black Holes from Dwarf Spheroidal Galaxies

技術要約：矮小楕円銀河からのメモリ負担を伴う原始ブラックホールのニュートリノ制約

問題提起
初期宇宙の過密領域から形成された原始ブラックホール（PBH）は、依然として有力なダークマター（DM）候補である。しかし、標準的な半古典的ホーキング蒸発は、質量が約 $5 \times 10^{14}$ g 未満のPBHが現在までに蒸発してしまったことを示唆している。ブラックホールの量子情報パラドックスに由来する「メモリ負担（Memory Burden, MB）」効果に関する近年の理論的発展は、PBHが初期質量の約半分を失った後、その蒸発率が大幅に抑制される可能性を示唆している。このメカニズムにより、PBHの生存可能な質量窓は $\sim 10^4$ g まで拡大される。これらのメモリ負担を伴うPBHは高エネルギー粒子（ニュートリノを含む）を放出すると予想されるが、その特有のシグネチャーについては、現在のニュートリノ観測器のデータを用いた包括的な探索が行われていない。ダークマター含有量が高く、天体背景放射が低い矮小楕円銀河（dSph）は、このような間接検出探索にとって理想的なターゲットとなる。

手法
著者らは、IceCubeニュートリノ観測器による10年間の公開ミューオン・トラック・データ（2008年4月〜2018年7月）を用いて、メモリ負担を伴うPBHからのニュートリノ信号の探索を行った。解析は、4つの独立した文献ソースから提供されたダークマター（D因子）測定値に基づく、14のdSphに焦点を当てている：

1. GSセット： 一様なジーンズ・モデリングを用いたGeringer-Samethらによるもの。
2. Eセット： 単純な解析公式を用いたEvansらによるもの。
3. Pセット： 経験的なスケーリング関係を用いたPace & Stigariによるもの。
4. Cセット： 更新された制約を提供するLHAASOコラボレーションによるもの。

本研究では、PBHのニュートリノ放出の有意性を評価するために、アンビン・最大尤度解析を採用している。尤度関数には以下が含まれる：

• 信号モデル： BlackHawk 2.30コードを用いて計算された、メモリ負担PBHの蒸発スペクトルから導出されたニュートリノ・フラックス。このモデルは、抑制因子 $k$（$k=1$ および $k=2$ をテスト）と、パラメータ $q$（通常 $q=0.5$ である、MB相への遷移を定義するもの）を考慮している。
• フレーバー振動： 解析はフレーバーに偏りのない初期フラックス（$1:1:1$）を仮定し、地球におけるフラックスを決定するためにポンテコルボ・マキ・ナカガワ・サカタ（PMNS）混合を適用しているが、全イベント率は特定の混合角に依存しないことが記されている。
• D因子処理： ダークマター分布の不確実性を扱うため、14のdSphに対して独立した誤差処理を可能にする非対称ガウス事前分布を用いたプロファイル尤度フレームワークを利用している。
• 統計的枠組み： 背景事象のみの仮説とPBH信号仮説を比較するために、尤度比に基づくテスト統計量（TS）を使用している。

主な貢献

• 包括的なデータセット： 4つの独立したD因子データセットを相互参照することにより、天体物理学的モデリングの変動に対処した、堅牢な14のdSphのサンプルを構築した。
• メモリ負担の適用： メモリ負担効果を、IceCubeデータを用いたdSphからのニュートリノ探索に適用した初の事例であり、特に $k=1$ および $k=2$ の抑制シナリオをテストしている。
• 結合解析： 14のすべてのdSphの結合解析を実施し、独立した誤差処理が単一ソースの解析と比較して統計的不確実性を効果的に低減することを実証した。

結果

• 検出なし： 10年間のIceCubeデータセットにおいて、背景事象のみの仮説に対する有意なニュートリノイベントの超過は検出されなかった。
• 存在量（$f_{PBH}$）の上限：
  • $k=1$ の場合： 本解析は、先行研究と比較して大幅に改善されたPBH存在量の絶対上限値を得た。$k=1$ の場合、結合解析により、質量約 $10^{10}$ g までのPBHがダークマターの主要成分であることを排除している。制約は、ニュートリノ・スペクトルがIceCubeの感度とよく一致する $10^7$ から $10^{10}$ g の間のPBH質量において特に厳しい。
  • $k=2$ の場合： エントロピー抑制がIceCubeが最も感度の高いエネルギー領域におけるニュートリノ・フラックスを強く減少させるため、制約は緩和される。
• データセット比較： Cセット（LHAASO）は、高いD因子により一般に最も強い制約を与えるが、GSおよびEセットも保守的な仮定の下で先行文献よりも厳しい制限を与えている。
• 個別 vs 結合： 全14のdSphの結合解析は、個別のdSph解析よりも、パラメータ空間全体にわたって実質的に厳格な排除限界を生み出す。

意義
本論文は、ニュートリノによるdSphからの探索が、メモリ負担を伴うPBHに対して、従来のガンマ線や宇宙線の研究と比較して競争力のある、そして（$k=1$ のシナリオにおいては）より優れた制約を提供することを確立している。著者らは、 $k \lesssim 1.3$ （または結合データセットを用いた場合は $k \lesssim 1.5$）の場合、PBHが $\sim 10^{10}$ g までの全ダークマターを構成することはできないという結果を強固に示していると主張している。より低質量のPBHに対する制約は既存の制限よりも緩いままであるが、本解析は、MBフレームワーク下でのより質量の大きいPBHに関連する硬いスペクトルに対する、ニュートリノ検出器の感度の高まりを浮き彫りにしている。本研究は、IceCube-Gen2のような次世代の観測装置を用いた将来の探索のための方法論的基礎を提供しており、これらは $k=2$ のシナリオをさらに深く調査し、質量到達範囲を広げることが期待されている。