[論文レビュー] Simulations of PBH formation at the QCD epoch and comparison with the GWTC-3 catalog

Albert Escrivà, E. Bagui+1 more|arXiv (Cornell University)|Sep 13, 2022

Quantum Chromodynamics and Particle Interactions被引用数 5

ひとこと要約

本研究では、時間変化する状態方程式を考慮した初期ブラックホール（PBH）形成の数値相対論的シミュレーションを、クォーク・グルーオン相転移期（QCD期）において初めて実施した。QCD相転移がPBH形成を促進し、太陽質量付近に顕著な質量ピークと20–100 M⊙範囲に二次的なピークを生じさせることを明らかにした。合体レートと重力波イベント分布は、PBHがダークマターの0.1～1の割合を占める範囲でGWTC-3カタログと一致するが、観測された8–15 M⊙の二次的ピークを再現できないため、PBHと恒星起源の混合起源である可能性が示唆される。

ABSTRACT

The probability of primordial black hole (PBH) formation is known to be boosted during the Quantum Chromodynamics (QCD) crossover due to a slight reduction of the equation of state. This induces a high peak and other features in the PBH mass distribution. But the impact of this variation during the PBH formation has been so far neglected. In this work we simulate for the first time the formation of PBHs by taking into account the varying equation of state at the QCD epoch, compute the over-density threshold using different curvature profiles and find that the resulting PBH mass distributions are significantly impacted. The expected merger rate distributions of early and late PBH binaries is comparable to the ones inferred from the GWTC-3 catalog for dark matter fractions in PBHs within $0.1 < f_{ m PBH} <1 $. The distribution of gravitational-wave events estimated from the volume sensitivity could explain mergers around $30-50 M_\odot$, with asymmetric masses like GW190814, or in the pair-instability mass gap like GW190521. However, none of the considered cases leads to a multi-modal distribution with a secondary peak around $8-15 M_\odot$, as suggested by the GWTC-3 catalog, possibly pointing to a mixed population of astrophysical and primordial black holes.

研究の動機と目的

• クォーク・グルーオン相転移期における一時的な状態方程式の低下が、PBH形成ダイナミクスに与える影響を調査すること。
• 変化する曲率プロファイルと時間依存のwをモデル化することで、より正確なPBH形成閾値を計算すること。
• 標準的な定数wモデルと比較して、QCD遷移がPBH質量関数と合体レートに与える影響を特定すること。
• 得られた重力波イベント分布を、GWTC-3カタログと比較し、特にGW190521やGW190814のようなイベントについて検証すること。
• QCDに起因する特徴が、LIGO/Virgoデータにおける観測された多峰性質量分布を説明できるかどうかを評価すること。

提案手法

• ビッグバン後約10⁻⁵ sにおけるQCD相転移を特徴とする時間依存状態方程式w(t)を用いたPBH形成の数値相対論的シミュレーション。
• 動的圧力勾配を考慮した、さまざまな初期曲率プロファイルを用いて、崩壊の臨界過密度閾値δcを計算すること。
• 変化するw(t)に基づくPBH質量関数の計算を行い、QCD遷移期における形成確率の増加に特に注目すること。
• PBHの集団形成と破壊を考慮した修正されたレート形式を用いて、初期および遅延期の二重PBH合体レートを推定すること。
• 予測された重力波イベント分布を、LIGO/Virgoの第3期観測における体積感度とGWTC-3カタログと比較すること。
• 広い質量関数における合体レートの人工的抑制を避けるために、ピーク領域に対して単色質量関数近似を用いること。

実験結果

リサーチクエスチョン

• RQ1 QCD期における時間変化する状態方程式が、PBH形成の臨界過密度閾値にどのように影響を与えるか？
• RQ2 QCD相転移を数値シミュレーションに組み込むと、得られるPBH質量関数はどのように変化するか？
• RQ3 QCDに起因するPBHから予測される合体レートと重力波イベント分布は、GWTC-3カタログの観測結果と一致するか？
• RQ4 なぜモデルはGWTC-3データにおける観測された8–15 M⊙の二次的ピークを再現できないのか？
• RQ5 広いPBH質量関数における集団形成と破壊効果が、合体レートをどの程度抑制するのか？

主な発見

• QCD相転移に起因する一時的なwの低下（1/3から約0.23に）により、PBH形成確率が顕著に増加し、太陽質量付近に顕著なピークを持つPBH質量関数が得られた。
• シミュレーションでは、中間子形成に起因する20–100 M⊙範囲に二次的ピークを持つPBH質量分布が得られ、GW190521の質量範囲と整合的である。
• PBHダークマター割合が0.1 < f_PBH < 1の範囲では、予測された合体レート分布がGWTC-3カタログの観測レートと一致した。
• 体積感度に基づく予測された重力波イベント分布は、GW190814（非対称質量）やGW190521（対安定性ギャップ）のようなイベントを再現した。
• 全体のレートと高質量特徴は一致するが、8–15 M⊙の二次的ピークはいかなるモデルでも再現されず、恒星起源とPBHの混合存在が必要である可能性を示唆した。
• 広い質量関数における標準的なレート抑制処理は人工的であることが判明した。QCDピーク領域に対して単色質量関数近似を用いることで、より信頼性の高い合体レート推定が可能となった。