Forecast for the detectability of patchy hydrogen reionization in WEAVE-QSO measurements of the Lyman-$\alpha$ forest power spectrum at redshift $z \geq 4$

技術要約：WEAVE-QSO観測におけるパッチ状水素再電離の検出可能性に関する予測

問題提起
宇宙間物質（IGM）は、水素再電離の時代からの「熱的遺物（thermal relics）」を保持している。再電離中、光加熱されたガスは電離バブル内で断熱冷却を起こし、その過程で長期間にわたって、IGMの温度分布に空間的な不均一性を刻み込んだ。これらの不均一性は赤方偏移 $z \sim 4$ まで持続し、ライマン$\alpha$（$Ly\alpha$）フォレストのフラックス・パワースペクトル（$P_{Ly\alpha}$）の大スケール（低波数、$k \lesssim 10^{-3} \, \mathrm{s\,km^{-1}}$）における増幅を予測させる。高解像度データ（KODIAK、SQUAD、XQ-100など）を用いた先行研究では、この増幅に対して暫定的な証拠（$\simeq 2.7\sigma$）が見出されているものの、統計的な精度と赤方偏移の経路長が限られていたため、高い有意性での検出には至らなかった。著者らは、次世代のWEAVE-QSOサーベイが、その膨大なサンプルサイズによって、この大スケールの遺物としての刻印を決定的に検出・特徴付けできるかどうかを予測することを目的としている。

手法
著者らは、WEAVE-QSOサーベイの構成に合わせた模擬$Ly\alpha$フォレスト・スペクトルを生成し、関連する不確実性を予測するための包括的なパイプラインを開発した。

1. 流体シミュレーション： IGMの基礎となる構造と熱履歴は、Sherwood-Relicsシミュレーション・スイートを用いてモデル化された。これらの宇宙論的流体シミュレーションは、再電離時およびその後のIGMを特異的に解像しており、予測される大スケールのパワー増幅を生じさせるパッチ状の再電離効果を組み込んでいる。
2. 模擬スペクトル生成（WEAVEify）： Sherwood-Relicsの出力を現実的なWEAVE-QSOスペクトルへとフォワードモデル化するために、カスタムソフトウェアパッケージであるWEAVEifyが開発された。このプロセスには以下が含まれる：
   • 光学厚さのスキュー（skewers）をQSOの静止系へマッピング。
   • （$R \sim 5000$を想定した）波長依存のガウス型カーネルによる畳み込み。
   • SDSSのQSOから導出された主成分分析（PCA）法による連続成分（continuum）の配置。
   • 検出器の特性、空の明るさ、および露光時間（$3060$ s）に基づき、ポアソン統計を考慮した、検出器の暗電流や読み出しノイズを含む現実的な光子制限ノイズの付加。
3. サーベイ・サンプリング： 著者らは、$z = 4.0, 4.2, 4.4, 4.6$ を中心とする4つの赤方偏移ビンに対して、$10,000$ 個の模擬QSOスペクトルを生成した。サンプルは、QSOの光度関数に基づくモンテカルロ法を用い、$r$バンドの等級カット $m_r < 21.9$ を適用して構築された。これにより、$10,000 \, \mathrm{deg}^2$ のサーベイ面積あたり、各ビンでそれぞれ $2,942, 1,822, 1,385, 1,155$ 個のQSOが期待される。
4. パワースペクトル推定と不確実性モデリング： 1Dパワースペクトルは、高速フーリエ変換（FFT）エスティメータを用いて推定された。著者らは、以下の要素を考慮した完全な共分散行列を構築した：
   • 統計的不確実性： 模擬サンプルサイズのブートストラップ法により導出。
   • 系統的不確実性： 分光分解能の補正、ノイズ減算、相関のない金属線汚染（C IV, Si IV, Mg II）、高カラム密度（HCD）吸収体の不完全性、および連続成分配置の誤差として個別にモデル化。
   • 解析は、分解能補正係数が $0.2$以上に保たれる波数$k < 0.05 , \mathrm{s,km^{-1}}$ に限定された。
5. パラメータ推論： 模擬データと理論モデルを比較するために、ベイズ的マルコフ連鎖モンテカルロ（MCMC）フレームワークが採用された。モデルには標準的なIGMパラメータ（有効光学厚さ $\tau_{\mathrm{eff}}$、温度 $T_0$、傾き $\gamma$、熱エネルギー $u_0$）に加え、パッチ状再電離の補正を表す追加の振幅パラメータ $A_p$ が含まれる。2つのシナリオがテストされた：「フィデューシャル（基準）」フィット（$A_p$ を変化させるもの）と、「no-$A_p$」フィット（パッチ状補正を無効にしたもの）である。

主な貢献

• WEAVE-QSOに関する初の詳細な予測： 本研究は、$z \ge 4$ における WEAVE-QSO サーベイによる $P_{Ly\alpha}$ 測定の精度に関する最初の具体的な予測を提供する。
• 包括的な系統誤差モデリング： 低分解能かつ大規模なサーベイに特有の系統的不確実性、特にHCDのダンピング翼（damping wings）や連続成分の配置ミスが大スケールのパワーに与える影響に関する厳密な扱いを導入した。
• 検出能力の検証： Sherwood-Relicsの物理学を現実的な計器モデリングと統合することで、パッチ状再電離の特異なシグネチャが、標準的なIGMの熱的変動とは明確に区別できることを本研究は示している。

結果

• 不確実性予算： $k_{\mathrm{max}}$ 以下の全スケールおよび全赤方偏移ビンで平均すると、予測される $P_{Ly\alpha}$ の総相対不確実性は約 $5.2\%$ である。大スケール（$k \lesssim 0.01 \, \mathrm{s\,km^{-1}}$）において、不確実性は $10\%$ 未満であり、中間領域では統計誤差が、極めて大きなスケールでは系統誤差（特にHCDと連続成分配置）が支配的である。
• 検出有意性： パッチ状再電離の信号を含む模擬データに対してMCMC解析を適用した結果：
  • フィデューシャルモデル（$A_p$ あり）は、モデルなし（no-$A_p$）よりも有意に優れた適合を示した（$\chi^2/\mathrm{dof}$ が $65.28/32$ から $36.07/28$ へ改善）。
  • この改善は、4つの追加パラメータに対して $\Delta\chi^2 = 29.21$ に相当し、これは両側ガウス等価有意性で $\simeq 4.5\sigma$ を意味する。
  • $A_p$ の事後分布は、すべての赤方偏移ビンにおいて入力された振幅を正確に回収している。
• スケールへの感度： 検出の有意性は最大スケールに強く依存している。最も低い $k$ ビン（$k \sim 10^{-3} \, \mathrm{s\,km^{-1}}$）を除去すると、有意性は $\simeq 1.2\sigma$ に低下し、このことは検出が熱的遺物の刻印の影響を最も受けるモードに依存していることを示している。

意義と主張
本論文は、WEAVE-QSOサーベイが、その予測されるサンプルサイズと共分散特性を備えていることにより、$z \ge 4$ における $Ly\alpha$ フォレストのパワースペクトルにおけるパッチ状水素再電離の大スケールな遺物としての刻印を、高い統計的有意性（$\simeq 4.5\sigma$）をもって直接検出できる可能性があると主張している。これは、より小規模で高解像度のデータセットによるこれまでの暫定的な検出（$\simeq 2.7\sigma$）から大幅な進歩である。著者らは、開発されたフレームワークにより、大スケールのパワースペクトル測定を通じて再電離の時期を評価することが可能であり、WEAVE-QSOの構成は、系統的不確実性が予測されたレベルで制御されている限り、暫定的な証拠を確実な検出へと変えるのに十分であると結論付けている。なお、本研究は再電離の正確な赤方偏移を測定することを主張しているのではなく、一様なUV背景モデルに対するパッチ状の刻印の振幅を測定することを目的としている。