Impact of 21-cm foreground mitigation strategies on reionization power spectrum constraints

技術要約：21cmフォアグラウンド除去戦略が再電離パワースペクトル制約に与える影響

問題提起
宇宙再電離（EoR）期における中性水素（HI）からの赤方偏移した21cm信号は、高赤方偏移の構造形成および初期銀河天文学を知るためのユニークなプローブを提供する。しかし、この宇宙論的信号は、数桁も強力な天体物理学的フォアグラウンドによって圧倒されてしまう。フォアグラウンドはスペクトル的に滑らかであるが、観測装置の特性（例：ビームの波長依存性、校正残差）がそのスペクトル構造を損ない、ラインオブサイト方向のフーリエモードへのパワーの漏出を引き起こす。これにより、$(k_\perp, k_\parallel)$ フーリエ平面内に「フォアグラウンド・ウェッジ」が生じ、検出に不可欠な「EoRウィンドウ」が制限される。本研究が取り組む中心的な課題は、これらフォアグラウンドの軽減戦略（具体的にはフォアグラウンド回避（Foreground Avoidance）とフォアグラウンド除去（Foreground Removal））が、天体物理パラメータの推論および全域的な再電離史に対して、どのように系統的なバイアスを伝播させるかを定量化することである。

手法
著者らは、未知のグラウンドトゥルース（真値）に対して解析チームが未知のまま進めるブラインド解析である、SKAO（Square Kilometre Array Observatory）の科学データチャレンジ3a（SDC3a）のデータを利用している。これは、シミュレーションされた1000時間のSKA-Lowデータセットである。本研究では、DOTSSチームによって提出された2つの具体的な軽減パイプラインを評価している：

1. フォアグラウンド回避（Foreground Avoidance）： 解析をフーリエ空間における汚染されていないEoRウィンドウに限定し、統計的なモデリングを行わずにデータの忠実性を確保するために大規模モードを犠牲にする。
2. フォアグラウンド除去（Foreground Removal）： ガウス過程回帰（GPR）と機械学習で訓練されたカーネルを用いて、スペクトル的に滑らかなフォアグラウンドをモデル化および減算し、ウェッジ内のモードの回収を試みる。

著者らは、ガウス尤度とnautilusダイナミック・ネステッドサンプラーを用いて、回収されたパワースペクトルのベイズ推論を行う。また、パワースペクトルを4つの自由な天体物理パラメータ（最小星形成ハロー質量 $M_{turn}$、星形成効率の正規化 $f_{\star,10}$、その質量依存性 $\alpha_\star$、および電離脱出率の正規化 $f_{esc,10}$）へとマッピングするために、21cmFASTで訓練されたエミュレータ21cmEMUを採用している。導出された量には、特定の赤方偏移における中性分率（$x_{HI}$）が含まれる。解析では、様々なシナリオ（全赤方偏移範囲、選択的な赤方偏移ビン、および変動するモデリング不確実性を含む）において、回収されたパワースペクトルを「フォアグラウンドフリー」のグラウンドトゥルースと比較する。

主な結果

• 天体物理パラメータにおける系統的バイアス： 両方の軽減戦略は、天体物理パラメータにおいて最大 $\approx 1\sigma$ の系統的バイアスをもたらす。
  • **回避（Avoidance）**は、大規模モード（$k \lesssim 0.2 \, h\text{Mpc}^{-1}$）の喪失により、一貫して事後分布が広くなる。また、最小ハロー質量（$M_{turn}$）をより高い値へとバイアスさせる。
  • **除去（Removal）**は、よりタイトな制約を与えるが、特に $M_{turn}$ を下限の事前分布境界（$7.0$）へと押し込み、$f_{esc,10}$をより低い値へとバイアスさせるなど、より強いバイアスをもたらす。これらのバイアスは、回収されたスペクトル、特に信号が複雑な構造を示す赤方偏移（例：$z \approx 8.90$）において、大規模なパワーが抑制されることによって引き起こされる。
• 赤方偏移依存の性能： 両手法とも、$z \approx 8.90$ において最も苦戦する。ここでは、X線加熱の開始と初期の再電離が同時に起こることで、信号の複雑さがフォアグラウンドの残差と区別することが困難になる。また、$z \approx 6.54$（再電離の後期）において、両手法は $x_{HI}$ に関して低い値への持続的なバイアスを示す。これは、電離バブルの重なりによる非線形性が、パワースペクトルの振幅と中性分率の結合をデカップリングしているためと考えられる。
• 全域的な履歴のレジリエンス（強靭性）： 源パラメータへのバイアスにもかかわらず、全域的な再電離史（$x_{HI}$）は、すべてのシナリオにおいて95%信頼区間内で概ね回収されている。このレジリエンスは、$x_{HI}$ が詳細な形状よりも、パワースペクトルの全体的な振幅に主に符号化されていることに起因する。
• 赤方偏移選択の影響： 回収が不十分なビン（具体的には $z=8.90$ と $z=6.54$）を除外すると、「最適な赤方偏移（Best redshifts）」シナリオにおいてパラメータのバイアスが大幅に減少し、ほとんどのパラメータでバイアスが $\approx 0.5\sigma$ まで低減する。ただし、実際の観測においてこれらのビンを特定するには、独立した品質指標が必要となる。
• 同一スケールへの収束： 除去（Removal）戦略を回避（Avoidance）と同じ $k$ モードに制限した場合、得られる事後分布は同様の領域に収束する。これは、両戦略間の相違が、主に除去によって回収された追加のモードによって駆動されていることを裏付けている。

意義と主張
本論文は、SKA時代におけるフォアグラウンド軽減の情報のトレードオフに関する明示的な評価を提供することを目的としている。その主要な貢献は、全域的な再電離史は頑健に回収される一方で、特定の天体物理源パラメータの推論は、選択された軽減戦略および回収されたパワースペクトルの品質に対して非常に敏感であることを示した点にある。

著者らは以下の点を強調している：

1. 現在のパワースペクトル解析において、軽減によるバイアスは避けられないものであり、$M_{turn}$ や $f_{esc}$ といったパラメータに最大 $1\sigma$ の影響を与える。
2. マルチ赤方偏移推論は、適切に軽減されていないビンからの汚染に対して脆弱であるため、実際の観測においては、それらのビンを特定し除外するための独立した品質指標（例：Lyman-$\alpha$ フォレストやCMBレンシングとの相互相関）が必要である。
3. 戦略の選択は観測対象に依存する： 回避と除去は同一スケールでは同様の結果をもたらすが、天球面へのアクセスを必要とする画像ベースの統計には除去が必要であり、一方で回避は感度の犠牲となるものの、モデルの誤設定に対する堅牢性を提供する。

本研究は、21cmパワースペクトルが依然として初期のEoR科学における主要なベンチマークであるが、将来の戦略としては、フォアグラウンドモデルを推論のフォワードモデルに直接組み込むか、標準的な軽減手法で失われた情報を回収するために高次統計量を利用する必要がある可能性があると結論付けている。