Ultra-slow-roll Inflation with Non-perturbative Non-Gaussianity and Scalar Induced Gravitational Waves

技術要約：非摂動的非ガウス性とスカラー誘起重力波を伴う超スローロール・インフレーション

問題提起
原始ブラックホール（PBH）は、特にアステロイド質量範囲（$\sim 10^{-15} M_\odot$）において、ダークマター（DM）の全密度を構成し得る有力な候補であり続けている。このようなPBHの形成には、宇宙マイクロ波背景放射（CMB）で観測されるものよりも著しく大きい原始曲率揺らぎが必要である。超スローロール（USR）インフレーションは、必要なパワースペクトルの増幅を生み出すことができるが、結果として生じるPBHの存在量を予測することは容易ではない。従来の研究の多くは、曲率揺らぎと密度ゆらぎの間の線形近似（$\delta \rho \propto \zeta$）や、非ガウス性（NG）に対する摂動論的な扱いへの依存度が高かった。しかし、USRフェーズでは$O(1)$の非ガウス性が生成されることがあり、これらの近似では不十分となる。さらに、大きな密度ゆらぎを介したPBHの形成は、必然的に二次オーダーのスカラ誘起重力波（SIGW）を生成する。したがって、CMBスケールからインフレーションの終焉までをモデル化し、非摂動的なNG効果によるPBH存在量を計算し、LISAのような将来の観測装置に向けたGW信号を決定するための、一貫したフレームワークが必要である。

手法
著者らは、2つのスローロール（SR）フェーズに挟まれたUSRフェーズを特徴とする、現象論的なインフレーションモデルを構築している。

• インフレーション・ポテンシャル: インフラトン・ポテンシャル $V(\phi)$ は、CMBスケールから終焉に至るインフレーションの全履歴を記述するようにパラメータ化されている。これは、二次領域（USR前）、完全に平坦なプラトー（USRフェーズ）、およびSRへのスムーズな遷移と有効な極小値を保証する三次領域（USR後）で構成される。モデルは、USRの継続期間（$\Delta N$）、遷移の滑らかさ（$h$）、およびUSR後の継続期間（$N_-$）を含む自由パラメータによって定義される。
• 曲率揺らぎ: ムクアノフ・サスキ方程式を数値的に解き、パワースペクトル $\mathcal{P}_\zeta(k)$ を得る。USRフェーズ中のモードの進化を追跡し、超地平線モードは凍結されるのではなく、指数関数的に増大し得ることに注意する。
• 非ガウス性（$\delta N$ 定式化）: NGに対する摂動展開を用いる代わりに、著者らは$\delta N$定式化を用いて、全曲率揺らぎ $\zeta$ とそのガウス成分 $\zeta_G$ との間の正確な非線形関係を導出する。これにより、$f_{NL} \sim O(1)$ の場合に極めて重要となる、揺らぎの完全な確率密度関数（PDF）を考慮することができる。
• PBH存在量: 密度ゆらぎ $\delta \propto \zeta$ という線形近似を回避する。代わりに、$\zeta$ と非線形に関係するコンパクション関数 $C$ を用いる。PBHの存在量は、最終的な質量に関するクリティカル・スケーリング則を組み込みながら、コンパクション関数のPDFを崩壊閾値にわたって積分することで計算される。
• 重力波: 曲率パワースペクトルから、二次オーダーのSIGWスペクトルを計算する。著者らは、非ガウス性がPBHに関連するPDFの裾（テイル）には影響を与えるものの、より広い範囲のモードを積分するSIGW信号への影響は、本モデルで見出された $f_{NL}$ の値においては無視できる程度であることを指摘している。

主要な貢献

1. 非摂動的なNGの取り扱い: 本論文は、$O(1)$ の非ガウス性が（USR-SR遷移から自然に生じ）ガウス的な予測と比較してPBHの存在量を著しく増大させることを示している。この増幅は、切り捨てられた展開ではなく、完全な非線形 $\zeta(\zeta_G)$ 関係を用いて導出される。
2. 現実的なポテンシャル・モデリング: 鋭いカットオフやUSR後のスケール不変なスペクトルを仮定するモデルとは異なり、本研究では、CMBスケールから終焉までインフレーションを追跡する連続的なポテンシャルを使用しており、パワースペクトルの特定のスケール依存性や、ピーク前の「ディップ」を捉えている。
3. PBHとGWの相関: 本研究は、PBH存在量（$f_{PBH}$）とSIGW信号強度の間の定量的な結びつきを確立している。それは、直感に反する結果を明らかにしている。すなわち、固定されたPBH存在量に対して、非ガウス性を高める（より鋭い遷移パラメータ $|h|$ を用いる）と、必要なパワースペクトルの振幅が減少し、それによって関連するGW信号が抑制される。
4. 観測的制約: 著者らは、将来の検出器（LISA, BBO, Einstein Telescope）に対する信号対雑音比（SNR）を計算し、PBHがダークマターである場合の除外領域をマッピングしている。

結果

• PBH存在量: 本モデルは、アステロイド質量窓（$\sim 10^{-15} M_\odot$）においてダークマターの100%を説明できるPBHを生成することに成功している。PDFの非ガウス的な裾は、同じパワースペクトルに対して、ガウス的な推定値よりも存在量を数桁増大させる。
• GW信号: 関連するSIGW信号は $O(0.1)$ Hz の周波数でピークに達する。
  • もしこの質量範囲のPBHが全ダークマター密度の重要な成分を占める場合、GW信号は、非ガウス性による抑制があるにもかかわらず、LISAの感度を十分に上回ると予測される。
  • $O(1)$ の非ガウス性の存在は、固定されたPBH存在量に対してGW信号の強度を減少させる。具体的には、$\Omega_{GW} \propto (f_{PBH})^{0.059}$ であることを見出しており、これは、PBHの存在量がモデルパラメータに対して高い感度を持つのに対し、GW信号はPBH存在量の変化に対して比較的鈍感であることを意味する。
• 検出可能性: 非ガウス性による抑制がある状況でも、PBHがダークマターの重要な成分を構成する場合、LISAおよびBBOのSNRは $>1$ となる。逆に、LISAによるノー結果（null result）は、ポテンシャルの局所的な特徴によって予測されるものよりも質量スペクトルが大幅に広い場合を除き、アステロイド質量のPBHがダークマターの重要な割合を占めるというシナリオを事実上閉鎖することになる。
• パラメータ依存性: 滑らかさパラメータ $h$ は遷移速度を制御する。より大きな $|h|$（より鋭い遷移）は、より大きな非ガウス性と、高 $k$ におけるパワースペクトルの急激な低下をもたらし、その結果、高周波数におけるGW信号の急激な減衰をもたらす。

意義および主張
著者らは、非ガウス性を非摂動的に扱い、現実的なインフレーション・ポテンシャルを使用することで、PBH存在量と関連するGW信号を予測するための、より堅牢で一貫したフレームワークを提供していると主張している。彼らは、非ガウス性がPBHの形成を促進する一方で、固定されたDM分率に対してGW信号を同時に抑制するという、単純なモデルでは見落とされがちなニュアンスを強調している。

本論文は、将来の重力波観測装置、特にLISAが、アステロイド質量のPBH形成からの確率的背景放射を検出するか、あるいはこのシナリオをダークマターの重要な成分としての可能性を排除できる可能性があると結論付けている。著者らは、自身のポテンシャルおよび非ガウス性の具体的な取り扱いが、スケール不変なスペクトルを仮定した先行研究よりも低いGW信号をもたらすことを述べているが、PBHがダークマターの重要な部分を構成する場合、その検出の挑戦はより厳しくなるものの、依然としてLISAの予測範囲内にあるとしている。