[論文レビュー] First-order Phase Transition interpretation of PTA signal produces solar-mass Black Holes

Yann Gouttenoire|arXiv (Cornell University)|Jul 9, 2023

Pulsars and Gravitational Waves Research参考文献 139被引用数 9

ひとこと要約

この論文はNG15とIPTA2 PTAデータのベイズ解析を行い、確率的GW背景が初期宇宙の強い1次相転移に起因する可能性を示し、太陽質量範囲の原始ブラックホールおよびマルチメッセージャー検出可能性への影響を示唆する。

ABSTRACT

We perform a Bayesian analysis of NANOGrav 15yr and IPTA DR2 pulsar timing residuals and show that the recently detected stochastic gravitational-wave background (SGWB) is compatible with a SGWB produced by bubble dynamics during a cosmological first-order phase transition. The timing data suggests that the phase transition would occur around QCD confinement temperature and would have a slow rate of completion. This scenario can naturally lead to the abundant production of primordial black holes (PBHs) with solar masses. These PBHs can potentially be detected by current and advanced gravitational wave detectors LIGO-Virgo-Kagra, Einstein Telescope, Cosmic Explorer, by astrometry with GAIA and by 21-cm survey.

研究の動機と目的

• PTA SGWBの起源を SMBH バイナリ以外にも探る動機づけ。
• NG15およびIPTA2 PTA 残差を適合させる能力を強い1次相転移 (1stOPT) が持つかを検証。
• 1stOPT パラメータの事前としてBBN境界とPBH過剰生成プリオリを組み込む。
• ベイズモデル比較を用いて1stOPTシナリオがSMBHバイナリと比較してどれだけ競合するかを定量化。
• PBHおよび将来のGW/天文学観測機器に対する観測的結果を予測。

提案手法

• アルファ>>1の性格を持つ強く過冷却された1stOPT からのSGWBをバルクフロー GWスペクトル（式6）でモデル化。
• enterprise/enterprise_extensionsとPTMCMC/PTArcadeを用いてNG15とIPTA2のタイミング残差に対してベイズ推論を実行。
• 1stOPTパラメータ（T_reh, beta/H, alpha）に対するBBN境界とPBH過剰生成プリオリを制約として組み込み。
• 積集合サンプリングに従ってベイズ因子（BF）によりモデルを比較（1stOPT単独、SMBHバイナリ、または両方）。
• 検出分析と完了率beta/Hの排除（下限）分析の両方を探索。

Figure 1: The violin diagrams depict the posterior probability distribution of the SGWB energy density in each frequency bins of NG15 and IPTA2 data sets. We overlay with solid lines the SWGB from 1stOPT, obtained using Eq. ( 6 ), using mean posterior value for the PT parameters. The dotted lines il

実験結果

リサーチクエスチョン

• RQ1 NG15とIPTA2で観測されたPTA SGWBを過冷却型の1次相転移が説明できるか？
• RQ2 BBNとPBHプリオリの下でPTAデータと一致する推定1stOPTパラメータ（T_reh, beta/H, alpha）はどれか？
• RQ3 PBH過剰生成制約を含めることで1stOPTとSMBHバイナリの解釈の優位性は変わるか？
• RQ4 太陽質量のPBHは1stOPTシナリオの実現可能な副産物であり、他の探査（GW検出器、GAIA、21-cm）で検出可能か？
• RQ5 PTA SGWBが1stOPT以外のものである場合の完了率beta/Hの保守的な下限はどれくらいか？

主な発見

• NG15データはSMBHバイナリよりも強い1stOPTからのGW信号を支持し、ベイズ因子は24（BBNのみプリオリ）。
• PBH過剰生成プリオリを含めると1stOPT vs SMBHのベイズ因子はNG15で15へ低下、IPTA2では未確定。
• NG15（BBN）下の平均後方分布の1stOPTパラメータ：log10(T_reh) = -0.80^{+0.23}_{-0.23}, beta/H = 9.8^{+4.0}_{-2.5}; PBHプリオリを含めると：log10(T_reh) = -0.86^{+0.37}_{-0.26}, beta/H = 10.7^{+2.8}_{-3.3}。
• IPTA2結果（BBN）：log10(T_reh) = -1.34^{+0.3}_{-0.3}, beta/H = 13.8^{+6.6}_{-4.2}。
• 1stOPTシナリオはNG15の後方領域で太陽質量PBHを生み出せる可能性があり、マルチメッセージャー窓（LIGO-Virgo-Kagra, ET/CE, GAIAレンズ効果, 21-cm加熱）を提供。
• 排除境界はβ/H ≳ 10–20 を示し、1stOPT由来でない場合、太陽質量1個から10^3太陽質量までPBHが検出可能なシナリオを制限。

Figure 2: Left : Colored regions are posterior distributions in term of the reheating temperature $T_{\rm reh}$ and rate of completion $\beta/H$ of a strong 1stOPT ( $\alpha\gg 1$ ). They are obtained after performing a Bayesian analysis of PTA dataset. We overlay the CMB, LIGO/Virgo and microlensin