ECHO21: a tool for modelling global 21-cm signal from dark ages to reionization

テクニカル・サマリー：「ECHO21：ダークエイジから再電離に至る全宇宙的な21cm信号のモデリングのためのツール」

問題と動機
ダークエイジ（暗黒時代）およびコスミック・ドーン（宇宙の夜明け）は、直接的な観測プローブが不足しているため、宇宙のタイムラインにおいて理解が進んでいない時期である。中性水素の21cm超微細遷移線は、第一世代の星、宇宙間物質（IGM）、および宇宙論パラメータの物理を解明するための潜在的なプローブとなる。21cmFASTやBEORNのような数値的または半数値的なコードは、正確で物理的に豊かな予測を提供するが、計算速度が遅いため、大規模なデータセットからの効率的な天体物理学的または宇宙論的な推論を妨げる。一方、ARESやZeus21のような既存の解析的コードは高速であるが、特定の物理成分、具体的にはLyman-$\alpha$（ライマン$\alpha$）光子による加熱効果が欠落している。著者らは、Ly$\alpha$加熱の欠如が信号モデリングに強いバイアスをもたらす可能性があることを指摘している。迅速で、十分に現実的であり、かつ広大なパラメータ空間を扱うことができる柔軟なツールへのニーズに応えるため、著者らはECHO21（Exploring the Cosmos with Hydrogen Observation）を導入する。

手法
ECHO21は、赤方偏移 $z \approx 1500$ から現在に至るまでの全域的な21cm信号をモデリングするために、解析的なフレームワークを利用したPythonパッケージである。このコードは、バルクのIGMと完全に電離したH II領域を区別する、IGMの二領域モデルに基づいて動作する。

• 物理的フレームワーク: コードは、ガス運動温度 ($T_k$)、バルクIGMの電子分率 ($x_e$)、およびH II領域の体積充填率 ($Q$) に関する結合微分方程式を解く。

  • 熱履歴: $T_k$ の進化は、断熱冷却、コンプトン加熱、X線加熱、そして極めて重要なことに、Ly$\alpha$加熱を考慮して計算される。著者らは、従来の解析的コードはLy$\alpha$加熱を無視していたが、ECHO21はこれを包含しており、星形成率密度（SFRD）から導出されたLy$\alpha$特異強度の $J_{Ly}$ に基づいて加熱率を算出することを述べている。
  • 電離履歴: コードは、X線およびUV光子によって駆動される再結合と電離を追跡する。H II領域の体積充填率は、電離光子の放出率と放射再結合を比較することによってモデル化される。
  • 星形成率密度 (SFRD): ECHO21の斬新な特徴は、3つの異なるSFRDモデルを含んでいることである：
    1. 物理的に動機付けられたモデル: 一定の星形成効率を用いた、ダークマターハローへのバリオンの崩壊に基づくもの（Press-Schechter、Sheth-Tormen、またはTinkerのハロー質量関数を使用）。
    2. 半経験的なモデル: 特徴的な星形成タイムスケールと、最小ビリアル温度における急峻なカットオフを組み込んだ物理モデルの変種。
    3. 経験的に動機付けられたモデル: UV光度関数データ（低 $z$ ではMadau-Dickinson形式、高 $z$ では指数関数的な減少）によって制約されたモデル。

• パラメータ: モデルは、6つの自由な天体物理パラメータ（Ly$\alpha$規格化係数 $f_{Ly}$、ライマン系列スペクトル指数 $s$、X線規格化係数 $f_X$、X線スペクトル指数 $w$、電離光子の脱出率 $f_{esc}$、および最小ビリアル温度 $min(T_{vir})$）によって制御される。さらに、標準的な宇宙論パラメータ（$H_0, \Omega_m, \Omega_b, \sigma_8, n_s, T_{\gamma0}, Y_p$）の変動も可能である。

• パフォーマンス: 解析的なアプローチにより、ECHO21は単一のCPU上で約 $O(1)$ 秒で単一の全域的21cm信号インスタンスを生成できる。このコードはMPI並列化が可能であり、数分間で数千のモデルを生成できるように設計されている。

主要な貢献

1. Ly$\alpha$加熱の包含: ECHO21は、IGMの熱進化にLy$\alpha$加熱を明示的に含めた最初の解析的コードである。著者らは、この項を無視すると、特にX線放射が弱くLy$\alpha$放射が強いシナリオにおいて、顕著なバイアスが生じることを示している。
2. SFRDモデルの多様性: コードは、SFRDのモデリングに対して3つの異なるアプローチ（物理的、半経験的、経験的）を独自に提供しており、ユーザーが星形成の記述の選択が結果としての21cm信号にどのように影響するかを評価できるようにしている。
3. 宇宙論における柔軟性: 多くの専門的なコードとは異なり、ECHO21は宇宙論パラメータを天体物理パラメータと同様に容易に扱うことができ、21cm信号に対する宇宙論的な不確実性（例：ハッブルテンション）の研究を容易にする。
4. 速度とスケーラビリティ: コードは迅速な推論のために最適化されており、50個のCPUを使用して、6次元の天体物理パラメータ空間における15,625個のモデルを約17分間で生成できる。

結果

• 基準モデル: 物理的に動機付けられたSFRDとPress-Schechterハロー質量関数を用いた基準モデルは、$z \approx 5.3$ まで完了する再電離史と、$z_{re} = 7.42$ という50%電離赤方偏移を予測する。電子散乱光学的厚さは $\tau_e \approx 0.0592$ と計算され、Planck 2018の制約と一致している。
• 信号特性: 基準となる全域的21cm信号は、$z=87.5$（ダークエイジ）において $-41.1$ mK の第1吸収ピーク、および $z=19.2$（コスミック・ドーン）において $-142.4$ mK のより深い第2吸収ピークを示し、続いて $z=11.2$ で $22.2$ mK の放出ピークを示す。
• Ly$\alpha$加熱の影響: Ly$\alpha$加熱を除外したモデルとの比較により、この項の欠如がより強い（深い）21cm吸収信号をもたらすことが示された。著者らは、Ly$\alpha$加熱は冷却効果を相殺し、X線パラメータに応じて非単調に信号強度を変化させ得るため、正確なモデリングにおける極めて重要な構成要素であると強調している。
• SFRDモデルの影響: SFRDモデルの選択は信号を大きく変化させる。物理的に動機付けられたモデルと比較して、星形成の開始を遅らせる経験的に動機付けられたモデルは、信号を約8.8の赤方偏移分遅延させ、約39.4 mK 浅くする。同様に、ハロー質量関数（HMF）の選択も、信号の特徴のタイミングと深さに影響を与える。
• パラメータ感度: 信号は $f_{Ly}, f_X, w,$ および $min(T_{vir})$ に対して最も敏感である。現在の観測的不確定性範囲内での宇宙論パラメータの変動は、信号振幅（例：ピーク吸収深さにおける $\sim 18$ mK の範囲）に対して比較的小さな変化しか与えない。

意義と主張
著者らは、ECHO21を、次世代の21cmデータセット（REACH、SARAS、または将来の宇宙ベースの実験など）から「迅速かつ十分に現実的な天体物理学的推論」を行うために設計された、柔軟で拡張性のあるオープンソースツールとして位置づけている。彼らは、本コードが以下の用途に特に適していると主張している：

• 大規模なパラメータ探索を通じた、堅牢な天体物理学的および宇宙論的推論の実行。
• 純粋に経験的なフィッティングではなく、物理的に動機付けられた信号を用いたデータ解析パイプライン（特にREACHコラボレーション向け）のテスト。
• 21cm観測と、BBN（ビッグバン元素合成）予測、CMB（宇宙マイクロ波背景放射）異方性、ハッブルテンションなどの他の宇宙論的プローブとの相乗効果の構築。

論文では、ECHO21が（複数のLy$\alpha$散乱を無視するなど、解析的な近似 due に）半数値的モデルよりも精度が低いことは謙虚に認めているが、その速度とモジュール性は、初期の推論や将来のバージョンにおける新しい物理の組み込みのための理想的なツールであることを強調している。コードはPyPIおよびGitHubを通じて公開されている。