Simulation of PBH formation in a matter-dominated universe

技術要約：物質優勢宇宙におけるPBH形成のシミュレーション

問題提起
原始ブラックホール（PBH）の形成は、支配的な宇宙成分の状態方程式に強く依存する。放射優勢期におけるPBH形成は比較的よく理解されているが、初期の物質優勢期におけるプロセスについては依然として不確実性が高い。物質優勢期では、圧力勾配が無視できるため、過密領域がより容易に崩壊する。しかし、形成基準は非球対称のダイナミクス、不均一性、および楕円性に大きく影響される。従来の解析的推定は、摂動近似やゼルドビッチ近似に依拠することが多いが、これらは高度な非線形領域においては成立しない可能性がある。この領域を調査する上での大きな課題は、数値シミュレーションにおいて非球対称性を実装することであり、特にダスト流体を扱う場合、ブラックホール・ホライゾンが完全に解像される前にシェルクロッシング（殻の交差）特異点によって数値計算がクラッシュしてしまう問題がある。

手法
著者らは、膨張する物質優勢背景下でのPBH形成をシミュレートするために、バウムガルテ・シャピロ・シバタ・ナカムラ（BSSN）形式を用いた完全非線形数値相対論を採用している。

• 初期条件： シミュレーションでは、楕円性を含む曲率摂動 $\zeta$ の関数形式を用い、三軸構成を作成している。初期データは、振幅パラメータ $\mu$ と楕円性パラメータ ($e, p$) によって特徴付けられる、ダスト流体の長波長勾配展開解を用いて構築されている。
• ゲージと座標： 新しい時刻スライス条件を採用しており、e-folding数に基づくゲージから、ホライゾン進入に近づくにつれて特異点回避型の「1+log」型スライシングへと遷移させている。また、シミュレーションボックスの中心領域の解像度を高めるため、非デカルト座標系を使用している。
• 物質の記述： 2つの異なる物質モデルをテストしている：
  1. ダスト流体： 無圧流体の記述。著者らは、この記述が必然的に、流体密度が発散するシェルクロッシング特異点を引き起こし、数値コードのクラッシュを招くことを指摘している。
  2. 衝突のない粒子： 流体近似の破綻を回避するため、著者らはフル数値相対論内に粒子インセル（PIC）法を実装している。粒子には所望の密度不均一性を再現するように質量が割り当てられ、その運動は測地線方程式に従う。これにより、流体要素が交差する地点を過ぎてもシミュレーションを継続することが可能となる。

主な結果

• 球対称性の検証： コードは、球対称なダスト崩壊のレマイトル・トルマン・ボンド（LTB）解に対して検証された。数値結果は、ホライゾン形成の時点まで解析解を正確に再現し、局所的な裸特異点（PBH形成）と全域的な裸特異点の間の閾値を正しく特定した。
• ダスト流体シミュレーション（非球対称）：
  • 三軸崩壊の場合、初期振幅 $\mu$ が十分に大きくない限り、ダスト流体シミュレーションは、見かけのホライゾンが形成される前にシェルクロッシング特異点によってクラッシュする。
  • ホライゾン形成がコードのクラッシュよりも先に観察されるために必要な閾値振幅は、ゼルドビッチ近似やホップ予想に基づく従来の解析的推定値の約2倍であることが判明した。
  • これは、流体近似が有効な範囲内では、PBH形成は以前考えられていたよりも困難であることを示唆している。
• 衝突のない粒子シミュレーション：
  • 同じシステムを衝突のない粒子でモデル化した場合、シミュレーションはシェルクロッシングによるクラッシュを起こさない。
  • 初期振幅 $\mu \approx 0.045$ （$e=0.2$ の場合）という低い値でも、ブラックホール形成（見かけのホライゾン形成によって示される）が観察された。
  • この閾値は、従来の解析的推定よりも桁違いに小さく、またダスト流体シミュレーションから導出された閾値よりも大幅に低い。
  • しかし、粒子シミュレーションは、粒子同士の交差が発生した後にハミルトニアン制約の著しい違反を示す。これは、崩壊自体は進行するものの、粒子記述と幾何学との間のミスマッチにより、後期進化の定量的精度が制限されていることを示している。

意義と主張
本論文は、PBH形成の基準が、物質場の物理的記述に依存して大きな理論的不確定性を有していることを示すと主張している：

1. 流体近似の限界： もし流体記述におけるシェルクロッシングに伴う特異点の形成が、物理的な抑制メカニズム（あるいは流体近似の破綻）として解釈されるならば、PBH形成の閾値は従来よりも大幅に高くなり、PBH形成はより困難になる。
2. 粒子ダイナミクス： もし衝突のない粒子記述が、単一ストリーム流体近似を超えた有効な拡張として受け入れられるならば、PBH形成は予想よりもはるかに容易（より低い閾値）になる。
3. 理論的不確定性： 流体と粒子の結果の相違は、標準的な解析基準（ゼルドビッチ近似）が三軸崩壊の非線形ダイナミクスを正確に捉えきれていない可能性を浮き彫りにしている。著者らは、閾値の値は高密度における物質場の取り扱いに極めて敏感であると結論付けている。

本論文は、PBHの存在量を決定的に解決することを目的としているのではなく、むしろ、物質モデルの扱いおよび特異点の数値的処理に対する感度を明らかにすることを目的としている。高密度領域の取り扱いの精緻化、粒子シミュレーションにおける制約保存の改善、および結果としてのPBH存在量と重力波シグネチャの定量化が必要な今後の課題として挙げられている。