Minkowski Functionals of the 21 cm Signal as a Probe of Primordial Features

技術要約：原始的特徴のプローブとしての21 cm信号のミンコフスキー汎関数

問題提起
宇宙の夜明け（Cosmic Dawn）および宇宙再電離（EoR）における赤方偏移した21 cm信号は、天体物理学的プロセスと原始物理学（特にインフレーション）の両方を解明するためのユニークなプローブである。インフレーションモデルは、標準的なスケール不変の原始パワースペクトルからの偏差を予測しており、これはしばしばインフレーション中の粒子生成に起因する「バンプ状」の特徴として現れるが、これらの特徴の検出は困難である。グローバルに平均化された21 cm信号やパワースペクトルを用いた先行研究では、特定のスケール、すなわち「ターンオーバー・スケール」（$k_{peak} \approx k_{turn}$）において、バンプモデルは標準的なフィデューシャルモデルと区別がつかないことが示されている。さらに、原始的なシグナルを、電離効率やX線加熱などの天体物理学的パラメータの変動から区別することも大きな課題である。21 cm信号は、再電離のパッチ状の性質により本質的に非ガウス的であり、原始的な密度のゆらぎの起源を特定するためには、高次統計量や形態学的記述子が必要となる。

手法
著者らは、21cmFAST v3 を用いた半数値シミュレーションを用いて、赤方偏移 $5 \lesssim z \lesssim 35$ における宇宙論的場の3次元的な実現形を生成している。

• 原始的特徴: 初期条件は、粒子生成から生じるバンプ状の構造を持つ原始パワースペクトルを含むように修正されている。これらの特徴は、振幅（$A_I$）とピーク位置（$k_{peak}$）によってパラメータ化される。本研究では、$k_{peak} \in \{0.4, 0.5, 0.6\} \, \text{Mpc}^{-1}$（ターンオーバー・スケール $k_{turn} \approx 0.5 \, \text{Mpc}^{-1}$ を含む）および振幅 $A_I \in \{10^{-9}, 10^{-8}\}$ に焦点を当てている。
• 天体物理学的モデリング: シミュレーションには、主要なEoRパラメータである電離効率（$\zeta$）、最小ヴィリアル温度（$T_{vir}$）、およびソフトバンドX線光度（$L_X$）が含まれている。
• 形態学的解析: 著者らは、密度場（$\delta_b$）、スピン温度（$T_S$）、中性水素分率（$x_{HI}$）、および輝度温度（$\delta T_b$）に対して**ミンコフスキー汎関数（MFs）**を適用している。4つのMF（$V_0, V_1, V_2, V_3$）は、様々な閾値におけるエクスカーション・セットの体積分率、表面積、平均曲率、およびオイラー標数を定量化する。
• 定量化: モデルを区別するために、著者らはMF曲線における非対称性やピーク構造を捉える派生可能な測定量である、振幅量化器（$A_1, A_2, A_3$）および形状量化器（$E_1, E_2, E_3$）を定義している。統計的有意性は、30個の実現形を用い、これら派生統計量の平均値を標準偏差単位でフィデューシャルモデルと比較することによって評価される。

主な貢献と結果

1. 原始的特徴への感度: 本研究は、MFsが原始的特徴の振幅とスケール（$k_{peak}$）の両方に極めて敏感であることを示している。グローバルな統計量とは異なり、MFsは21 cm場の幾何学的形状やトポロジーに関する豊かな形態学的情報を取り込む。
2. ターンオーバー・スケールでの検出: 重要な知見は、MFsがターンオーバー・スケール（$k_{peak} \approx 0.5 \, \text{Mpc}^{-1}$）においても、バンプモデルをロバストに識別できることである。このスケールでは、グローバルに平均化された21 cm信号や再電離の履歴はフィデューシャルモデルと区別がつかない。これは、バンプモデルが異なる質量スケールにわたってハロー質量関数（HMF）を変化させ、平均化された量が変化しない場合でも、崩壊する構造の分布を修正するためである。
3. 赤方偏移依存のシグネチャ:
   • 高赤方偏移 ($z > 20$): 信号は密度ゆらぎによって支配されている。バンプモデルは、この領域においてフィデューシャルモデルおよびEoRパラメータの変動の両方からロバストに区別可能であり、原始物理学を制約するためのクリーンな窓を提供している。
   • 中間赤方偏移 ($10 < z < 20$): 区別可能性は微妙であり、特定の統計量や $k_{peak}$ の値に依存する。ターンオーバー・スケールのモデルは、特定の狭い赤方偏移ウィンドウ（例：$z \sim 16\text{--}17$ 付近）においてのみ識別可能である。
   • 低赤方偏移 ($z < 10$): 原始的特徴とEoRパラメータ（特に $\zeta$ と $L_X$）との間に顕著な縮退が存在し、分離は困難であるが、特定の統計量については限定的な識別可能ウィンドウが存在する。
4. 場特有の形態: 解析により、原始的特徴が構成要素となる場を通じてどのように伝播するかが明らかになった。より低い $k_{peak}$ 値は大規模な過密領域を強化してMFの振幅を増大させる一方で、より高い $k_{peak}$ 値は構造形成と再電離を加速させ、電離バブルのタイミングと形態を変化させる。

意義
本論文は、ミンコフスキー汎関数が、迫り来る21 cmサーベイにおいてインフレーション物理学を探索するための、パワースペクトルやバイスペクトルに補完的な強力なツールであることを主張している。非ガウス的な形態学的シグネチャを捉えることで、MFsは、グローバルな統計量が機能しない領域において、原始的特徴を天体物理学的な不確実性から分離することができる。著者らは、現在の解析は観測ノイズのない理想化されたシミュレーションを使用しているものの、その枠組みは一般的であり、スクエア・キロメートル・アレイ（SKA）からの将来のデータにも適用可能であることを強調している。本研究は、形態学的統計量を用いてインフレーションモデルや原始的な密度ゆらぎの起源を制約するための基礎を築いており、これらのシグネチャをロバストに特定するためには、複数の赤方偏移にわたるMFを組み合わせることが不可欠であることを示唆している。