Probing Memory-Burdened Primordial Black Holes with High-Energy Neutrinos

技術要約：高エネルギーニュートリノによるメモリ負担を伴う原始ブラックホールの探索

問題と動機
原始ブラックホール（PBH）は、初期宇宙において形成された仮説上の天体であり、ダークマターの一部を構成している可能性がある。標準的な半古典的描像では、初期質量が $\sim 10^{15}$ g 未満の PBH は、ホーキング放射を通じて現在までに完全に蒸発しているはずである。しかし、「メモリ負担効果（memory-burden effect）」は、この蒸発プロセスに対する修正を提案している。この枠組みでは、PBH に蓄えられた情報が成長するメモリ負担を生み出し、放出率を抑制することで、より軽い PBH が今日まで生存することを可能にする可能性がある。これらの生存している軽い PBH は、より高温であり、高エネルギー（TeV）から超高エネルギー（EeV）のニュートリノを放出する。

メモリ負担を伴う PBH に関するこれまでの研究は、主に単色質量関数（単一の初期質量）に依存していた。しかし、原始的な密度揺らぎの崩壊といった現実的な形成メカニズムは、一般に拡張された質量関数を予測する。本論文は、拡張された質量分布（具体的には対数正規分布）が、単色仮定と比較してメモリ負担のある PBH に対する制約をどのように変化させるのか、また将来のニュートリノ望遠鏡がこれらの仮説を区別できるのかという理解における空白を埋めるものである。

手法
著者らは、2段階の蒸発モデルを用いて、メモリ負担のある PBH からの拡散ニュートリノフラックスを分析している。このモデルは、初期の標準的なホーキング相と、PBH の質量が初期値の $q=1/2$ に減少した時に始まる抑制されたメモリ負担相で構成される。抑制の強さはパラメータ $k$ によって制御され、$k$ が大きいほど抑制が強く、寿命が長くなる。

本研究では、以下の2つの初期質量関数を用いる：

1. 単色（Monochromatic）: $\psi(M_i) = \delta(M_i - M_0)$。
2. 対数正規（Log-normal）: $\psi(M_i) \propto \frac{1}{M_i} \exp\left[-\frac{\ln^2(M_i/M_c)}{2\sigma^2}\right]$。これは中央質量 $M_c$ と幅 $\sigma$ によって特徴付けられる。

分析には以下の3つのデータセットを利用する：

• 現在のデータ: IceCube HESE（高エネルギー開始イベント）、MESE（中エネルギー開始イベント）、および EHE（極高エネルギー）データ。
• 将来の予測: 10 年間の露光とバックグラウンドのないゼロカウント・シナリオを想定した、IceCube-Gen2 ラジオおよび GRAND200k。

統計的手法には以下が含まれる：

• 頻度論的解析: バックグラウンドを考慮しない尤度を用いた、現在における PBH 存在量（$f_{\text{PBH}}$）の 90% 信頼区間（CL）上限の導出。
• ベイズ的解析: シミュレーションされた 30 イベントの信号を含むデータセットを用いたモデル比較（単色 vs 対数正規）およびパラメータ再構成。これには、ベイズ因子（$B_{\text{LN,mono}}$）およびパラメータ（$M_0/M_c, k, \sigma, f_{\text{PBH}}$）の事後分布の計算が含まれる。

主要な貢献

1. MESE データの導入: 本研究は、IceCube MESE データを用いたメモリ負担のある PBH に対する初の制約を提示している。
2. 拡張質量関数の分析: 単色質量近似を超え、単色質量関数と対数正規質量関数を比較することで、制約がどのように変化するかを体系的に分析した。
3. ベイズ的識別研究: 将来のラジオ検出器が、質量関数仮説を区別し、PBH パラメータを再構成する能力について、これまで定量化されていなかった予測を提供した。

結果

• 存在量のリミット:

  • 低いメモリ負担（$k=1$）の場合、現在の IceCube データ（主に HESE）が最も強い制約を与える。将来のラジオ検出器は、生存している PBH が重すぎて、すでに現在のデータでカバーされているエネルギー範囲のニュートリノを放出するため、これらの制約を改善することはない。
  • 高いメモリ負担（$k=2, 4$）の場合、生存している PBH はより軽く、より高温となり、ニュートリノスペクトルが高エネルギー側にシフトする。この領域では、予測される IceCube-Gen2 ラジオおよび GRAND200k の感度は、現在の制約よりも大幅に強力である。
  • 質量関数の影響: 対数正規分布は、同じ中央質量を持つ単色分布よりも大幅に強い制約を与える。これは、対数正規分布の低質量側（low-mass tail）に、より効率的にニュートリノを放出する、より軽く高温の PBH が含まれているためである。$\sigma=1$ かつ $k=2, 4$ の場合、対数正規分布による制約は、同じ特性質量における単色分布の制約よりも数桁（$10^{-6}$ から $10^{-9}$）厳しくなる。

• モデル識別:

  • シミュレートされた 30 イベントの信号を用いた研究では、質量関数の識別は困難であることが示された。
  • $k=2$ の場合、真の分布が対数正規分布であれば、結合データセットは対数正規仮説に対して「非常に強い」証拠（$\ln B \approx 6.6$）を与える。しかし、真の分布が単色分布である場合、正しいモデルに対する証拠は「肯定的」ではあるが強くはなく（$\ln B \approx -1.88$）、対数正規モデルは強く否定されることもない。
  • $k=4$ の場合、識別は失敗する。強い抑制を伴う単色集団と、拡張された質量関数のスペクトル的特徴が退化（degenerate）するため、誤った、あるいは決定不可能なモデル選択につながる。

• パラメータ再構成:

  • $k=2$ の場合、単色モデルは質量（$M_0$）とメモリ負担（$k$）の正確な再構成を可能にする。
  • 対数正規モデルは、顕著な退化をもたらす。幅のパラメータ $\sigma$ は弱くしか制約されず、$M_c$ は $\sigma$ および $k$ と相関している。その結果、特性質量は $\sim 0.7$ dex の範囲内でしか決定できず、$f_{\text{PBH}}$ も両方のケースで十分に再構成されない。
  • $k=4$ では、パラメータの退化がより強いため、再構成能力はさらに低下する。

意義と結論
本論文は、初期質量分布が、メモリ負担効果による PBH からのニュートリノ信号を解釈する上で極めて重要な要因であると結論付けている。対数正規集団の中央質量において単色分布を仮定することは、制約を（数桁にわたって）過小評価することにつながる。

高エネルギーニュートリノは、メモリ負担効果によって可能となった軽い PBH を直接探査する手段を提供するが、著者らは、質量分布（単色か広範な対数正規か）を区別することは、特に高いメモリ負担（$k=4$）の値において、スペクトルの退化により困難であると警告している。30 イベントの信号があったとしても、データは基礎となる質量関数を確実に特定したり、物理パラメータを精密に再構成したりするには不十分な可能性がある。本研究は、拡張された質量関数が高エネルギーの裾野を通じて発見の可能性を高める一方で、パラメータの退化を増大させることで、根本的な PBH の特性抽出を複雑にするという側面を浮き彫りにしている。