Curbing PBHs with PTAs

技術要約：PTAによるPBHの抑制

問題提起
原始ブラックホール（PBH）は、ダークマターの有力な候補であり、超巨大ブラックホール（SMBH）の種となる可能性を持つ。その形成には、宇宙マイクロ波背景放射（CMB）が探査するスケールよりも小さなスケールにおいて、原始曲率パワースペクトル（$P_\zeta$）の大幅な増強が必要である。しかし、このような大きな曲率摂動は、必然的に二次的な信号であるスカラー誘起重力波（SIGW）を生成する。パルサー・タイミング・アレイ（PTA）、具体的にはNANOGrav 15年間のデータセットは、ヘリングス・ダウンズ・パターンを伴う確率的重力波背景放射（SGWB）を検出している。本論文は、ダークマターやSMBHの種を説明するために必要なPBHの存在量と、PTAの観測によって課されるSIGW振幅の上限との間の緊張関係を調査する。この調査における重要な側面は、原始非ガウス性（NG）の役割であり、これはPBHの形成確率と、結果として生じるSIGWスペクトルの両方を著しく変化させる可能性がある。

手法
著者らは、PBH存在量に対する制約を導出するために、多段階の理論的枠組みを用いている：

1. パワースペクトルモデル： 本研究では、増強された原始曲率パワースキクトルのベンチマーク形状として、対数正規（LN）分布と、異なる幅を持つ2つの破断冪乗則（BPL）モデルの計3つを検討する。
2. 非ガウス性（NG）定式化： 原始曲率摂動 $\zeta$ は、二次（$f_{NL}$）および三次（$g_{NL}$）の項を含む局所展開を用いてモデル化される。著者らは、密度ゆらぎの確率密度関数（PDF）へのNGの影響と、誘起GWスペクトルへの影響の両方を分析する。
3. PBH存在量の計算： PBH質量分率（$f_{PBH}$）を推定するために、2つの異なる定式化を使用する：
   • 閾値統計（Threshold Statistics）： コンパクト関数に基づき、密度ゆらぎが臨界閾値（$C_{th}$）を超える確率に基づく。このアプローチは、曲率と密度摂動の間の非線形関係を考慮している。
   • 極値理論（Theory of Peaks）： 曲率場の局所的な極大値の統計に基づく。著者らは、NGをこの定式化に含めることは技術的に困難であるため、比較のための分析は主にガウス限界に限定されていると述べている。
4. SIGW計算： SIGW背景放射は、摂動論の二次のオーダーで計算される。著者らは、$f_{NL}$ がスペクトルの振幅と形状に与える影響を組み込み、GWエネルギー密度スペクトル（$\Omega_{GW}$）を導出する。
5. PTA制約： 理論的なSIGW予測は、14個の周波数ビンにおけるNANOGrav 15年間のデータの95%信頼区間上限と比較される。解析は、異なるウィンドウ関数（Dirac、Sinc、Top-Hat）を用いてスペクトル分解能効果を考慮している。
6. SMBH種としての生存可能性： 本論文は、大きな正のNGが、CMBスペクトル歪み（FIRAS）の制限に抵触することなく、大規模なSMBHの種となるための十分な質量のPBH（$10^3 - 10^6 M_\odot$）の存在量を可能にするかどうかを検証する。

主要な貢献と結果

• 恒星質量PBHへの制約：

  • ガウスゆらぎ（または非線形な密度-曲率関係のみに由来するNG）を仮定した場合、PTA観測は曲率パワースペクトルの振幅に対して厳格な上限を課す。
  • 閾値統計を用いると、著者らは、恒星質量のPBHがダークマターのわずかな割合しか占めることができないことを見出した。狭いスペクトル（LN、BPL1）の場合、$0.1 - 100 M_\odot$ の質量範囲に対して $f_{PBH} \lesssim 10^{-5}$ という制約が適用される。より広いスペクトル（BPL2）の場合、制約はより重い質量（$5 - 500 M_\odot$）に適用される。
  • 極値理論を用いると、閾値統計と比較してPBH存在量への制約が数桁緩和される。これは、PBH存在量の推定における顕著な理論的不確実性を浮き彫りにしている。

• 原始非ガウス性の影響：

  • 正のNG ($f_{NL} > 0$): これらは二重の効果を持つ。これらはSIGW振幅を増加させ（パワースペクトル振幅 $A$ の制約を強める）、より重要なことに、大きな密度ゆらぎの確率を指数関数的に増大させる。その結果、同じPBH存在量を生成するためにより小さなパワースペクトル振幅が必要となる。総体的な結果として、正のNGはPBH存在量に対する制約を緩和する。$f_{NL} \gtrsim 5$ の場合、広いスペクトルにおける恒星質量範囲のPBH存在量に対するPTAの制約は事実上消失する。
  • 負のNG ($f_{NL} < 0$): これらはPDFの裾を抑制し、特定のPBH存在量を達成するためにより大きなパワースペクトル振幅を必要とする。これは、$f_{PBH}$ に対するより厳しい制約につながる。最も厳格な制約は $f_{NL} \approx -2$ 付近で発生する。
  • 超スローロール（USR）モデル: スペクトルの傾きによってNGが生成される特定のUSRインフレーションシナリオでは、誘起されるNGは比較的小さい（$f_{NL} \sim 1$）。これらの場合、PTAの制約は依然として有効であり、特にNGによる抑制が弱い広いパワースペクトルの場合に有効である。

• SMBH種と大きなNG:

  • 本論文は、「SMBH種問題」に対処する。これは、標準的なガウスモデルでは、CMBスペクトル歪みの制限に抵触することなく、十分な質量のSMBHの種となるPBHを生成できないという問題である。
  • 著者らは、大きな三次非ガウス性（$g_{NL} \sim 10^5$）がこの障害を克服できることを示し、FIRAS $\mu$-歪み制限と整合性を保ちながら、SMBHの種となる十分な存在量の $10^3 - 10^6 M_\odot$ のPBHを可能にすることを実証した。ただし、そのような大きな $g_{NL}$ を生成する具体的なモデルは、現在文献には欠けていると注記している。

• 理論的不確実性：

  • 存在量計算における定式化の選択（閾値統計 vs 極値理論）は、導出される制約にオーダーの違いをもたらす。
  • 強いNGが存在する場合の崩壊閾値 $C_{th}$ の扱いは不確実性を伴うが、著者らは摂動による $C_{th}$ への補正が主要な結論を劇的に変えることはないと推定している。

意義と主張
本論文は、最新のPTAデータを用い、原始非ガウス性の影響を明示的に組み込んだ、PBH存在量に関する更新された保守的な境界を提供すると主張している。その主な意義は以下の点にある：

1. PTAデータは初期宇宙物理学の強力なプローブであることであり、標準的なガウス仮定の下では、恒星質量範囲におけるPBHをダークマターの主要成分から排除できる。
2. 非ガウス性は決定的な縮退要因である： 正のNGは制約を「抑制（curb）」し、ガウス仮定の下で考えられていたよりも高いPBH存在量を可能にする一方で、負のNGは制約を強める。
3. 理論的な堅牢性は限定的である： 閾値統計と極値理論から導かれる制約の間の大きな乖離は、PBH形成の閾値および非ガウス領域における存在量のより精密な計算の必要性を強調している。
4. SMBHの種となることは可能だがモデル依存である： 大きなNGは、PBHをSMBHの種とするための理論的な経路を提供するが、そのような大きなNGを生み出す具体的なモデルの欠如が、このシナリオを確定させる上での障壁となっている。

著者らは、PTA観測はPBHダークマターのパラメータ空間を著しく制約するものの、PBH存在量の最終的な判断は、崩壊閾値および原始非ガウス性の具体的な性質に関する理論的不確実性に左右されると結論付けている。