Forest without Trees is still Fruitful: Constraints on the thermal state of the neutral IGM at $z\approx5.6$ with the 21-cm forest power spectrum

技術要約：21cmフォレストのパワースペクトルによる $z \approx 5.6$ における中性IGMの熱状態への制約

問題と動機
再電離の最終段階（$z \gtrsim 5.5$）における、主に中性である銀河間物質（IGM）の熱的および電離状態は、依然として十分に制約されていない。Ly$\alpha$ フォレストのような従来のプローブは、急速な飽和のために、相当量の中性領域に対して感度が低い。また、全天的な信号やトモグラフィー的な21cm放射実験は、明るい天体前景放射による深刻な汚染に直面している。21cmフォレストは、明るい無線力強いクエーサー（RLQSO）のスペクトル上に、中性IGMによって刻み込まれる狭い吸収特徴であり、前景放射の影響をほとんど受けない微分吸収プローブである。しかし、理論的な予測にもかかわらず、これまでターゲットとされたソースの赤方偏移が低く、感度も不十分であったため、21cmフォレストはまだ統計的に検出されていない。本研究は、直接検出がない状況下で、統計的な観測量としてフォレストの1次元（1D）パワースペクトルを利用することで、この課題に対処するものである。

手法
著者らは、$z > 5.5$ ($z=5.82$) において最も明るい既知のRLQSOである J352–15 (PSO J352.4034–15.3373) を対象とした、アップグレードされたGiant Metrewave Radio Telescope (uGMRT) のアーカイブデータを分析した。

• 観測: バンド2レシーバーを用いた合計17.5時間のオンソース時間を、RFI（無線周波数干渉）フラッギングとキャリブレーションを経て、実効積分時間11.73時間として減少させた。データは、サイドローブのアーティファクトを軽減するために高度な方向依存キャリブレーション（DDFacetおよびCUBICAL）を用いて処理され、その後、連続波減算とウェーブレットベースの残留バンドパス除去が行われた。
• データプロダクト: 209.0–222.5 MHz ($5.4 \le z \le 5.82$) にわたる無線スペクトルが得られ、チャンネルあたりのRMSノイズは 3.62 mJy beam$^{-1}$ (分解能 6.1 kHz) であった。主要な観測量は、正規化されたフラックスゆらぎ ($\delta F_{21}$) の1次元（1D）パワースペクトルであり、RFI除去による不規則なサンプリングを扱うためにLomb–Scargleピリオドグラムを用いて計算された。
• フォワードモデリング: 半数値的なコード 21cmFAST (v3.3.1) を用いて、534個のIGMモデルのグリッドにわたる合成21cmフォレストスペクトルを生成した。これらのモデルは、平均中性水素分率 ($\langle x_{\text{HI}} \rangle$) と、中性IGMの平均温度 ($\langle T_{\text{HI}} \rangle$) を決定するX線加熱効率 ($f_X$) を変化させている。シミュレーションでは、スピン温度 ($T_S$) と運動学的温度 ($T_K$) の間の飽和した結合を仮定している。
• 推論: ガウス尤度を用いたベイズ・マルコフ連鎖モンテカルロ（MCMC）解析を実施した。観測されたパワースペクトルは、インストルメンタルノイズ（オフソースの視線から推定）とサンプル分散の両方から導出された共分散行列を考慮した上で、フォワードモデリングされた合成スペクトルと比較された。解析では、予測される信号が観測された残差を超える領域を除外することで、非検出を制約として扱った。

主な結果

• 検出なし（Null Detection）: 測定されたJ352–15の1Dパワースペクトルは、インストルメンタルノイズと一致している。ターゲットのスペクトル分布とオフソースのノイズ分布を比較する統計テスト（コルモゴロフ–スミルノフ検定およびアンダーソン–ダリング検定）の結果、21cmフォレストの検出を示す有意な証拠は得られなかった（$p$-値 = 0.353）。
• 熱的制約: 検出なしであるにもかかわらず、本解析は $z \approx 5.6$ におけるIGMの状態に対して共同の制約を課している。68%信頼区間において、データは低温で相当量の中性IGMを持つモデルを否定している。具体的には：
  • $\langle x_{\text{HI}} \rangle = 0.1$ の場合、$\langle T_{\text{HI}} \rangle \lesssim 26.9$ K のモデルが除外される。
  • $\langle x_{\text{HI}} \rangle = 0.2, 0.5, 0.8$ の場合、温度の下限はそれぞれ $\langle T_{\text{HI}} \rangle > 39.8, 64.7, 79.9$ K である。
  • $\langle x_{\text{HI}} \rangle \gtrsim 0.05$ かつ $\langle T_{\text{HI}} \lesssim 18$ K のモデルは除外される。
• 他のプローブとの比較: これらの制約は、Ly$\alpha$ フォレスト（電離状態は制約できるが、中性領域の温度は制約できない）および既存の21cmトモグラフィーの制限（現在はより高い赤方偏移に限定されているか、多くのパラメータを周辺化している）を補完するものである。21cmフォレストの結果は、中性IGMが断熱冷却の底よりも高いレベルまで、実質的な予熱を受けていることを示唆している。
• MWA解析: アーカイブされたMurchison Widefield Array (MWA) データの独立した解析も行われた。感度が低く、視野が広いこと（サイドローブ汚染の増大につながる）により、MWAのデータは著しく弱い制約しか得られず、最終的な推論には含まれなかった。

意義と主張
論文は、21cmフォレストが「観測的に情報価値のある統計量」の領域に入ったと主張している。主な意義は、統計的な検出が得られない場合でも、パワースペクトルを観測量として用いることで、中性IGMの熱的および電離状態に関する有意義な制約を導き出せることを示した点にある。

• 本研究は、EoR期における21cmフォレストの1Dパワースペクトルの最初の観測的推定値を提供する。
• 現世代の無線施設（uGMRT）が、既存のLy$\alpha$ および21cmトモグラフィーの観測によって許容されているパラメータ空間を探索可能であり、特に極めて低温な中性IGMのシナリオを否定できることを確立した。
• 著者らは、将来のより深い観測（例：仮想的な100時間のキャンペーン）によって、統計的な検出がなされるか、あるいは現在の制限よりも温度が30〜100%高いモデルを排除するなど、これらの制約が大幅に強化されるであろうことを強調している。

限界
著者らは、モデリングにおけるいくつかの限界を認めている：

1. 近傍効果: クエーサーからの放射フィードバックや、潜在的なホスト銀河による吸収はモデル化されていないが、クエーサーに近い領域は解析から除外されている。
2. 小スケール構造: 解像度の制約により、未分解のミニヘイローや高密度構造は含まれていない。
3. 赤方偏移の進化: 解析は共存するIGMを仮定しており、観測された赤方偏移区間（$5.4 < z < 5.8$）にわたるIGMの状態の進化を無視している。
4. 電離モデリング: 21cmFASTの離脱集合形式（excursion-set formalism）は、主に完全に電離したセルまたは中性のセルを生成するため、電離フロントにおける部分的に電離した構造の詳細な取り扱いを欠いている。
5. 系統誤差: ノイズモデルは、将来のより深い観測において残存する系統誤差やキャリブレーション誤差が支配的になる可能性を考慮せず、相関のないガウスノイズを想定している。

結論として、本研究は、21cmフォレストが中性IGMの熱史を探るための実行可能なプローブであることを検証し、この分野を理論的予測から観測的制約へと移行させた。