HETDEX における Lyα 強度マッピング：銀河-Lyα 強度交差パワースペクトル測定の技術的サマリー

本論文は、ホビー・エバーリー望遠鏡ダークエネルギー実験（HETDEX）のデータを用いて、Lyα 強度マッピング（LIM）によるパワースペクトルの測定を行った研究報告です。特に、検出された Lyα 放射銀河（LAEs）と、検出限界以下の未検出源からの Lyα 強度との交差パワースペクトルを初めて高統計的有意性で検出しました。

以下に、問題設定、手法、主要な貢献、結果、および意義について詳細をまとめます。

1. 問題設定と背景

• 線強度マッピング（LIM）の意義: 従来の銀河サーベイが個々の銀河を「ピーク」として検出するのに対し、LIM はノイズを含むデータ全体から、検出限界以下の銀河や拡散ガスの光子を総和して、宇宙の大規模構造を tracer として捉える手法です。
• 既存研究との矛盾: これまでの研究（Croft et al. 2016, 2018; Lin et al. 2022）では、クエーサー（QSO）と Lyα 強度の交差相関を用いて宇宙の Lyα 光度密度を制限してきました。しかし、これらの結果は Lyα 森林からの上限値と矛盾しており、QSO 周辺の放射場が結果を支配している可能性が指摘されていました。
• HETDEX の可能性: HETDEX は事前選択なしに広範囲を観測するため、検出された LAE は信号の約 0.01% しか占めず、残りの大部分は「未検出源」からなる強度マップを提供します。これは、QSO 依存性を排除し、星形成銀河や IGM/CGM からの Lyα 放射を直接探る理想的なデータセットです。
• 課題: 強度マッピングでは、天体物理的信号よりもはるかに大きな天体測光ノイズや大気発光（Sky emission）の残差が支配的です。これらを除去しつつ、大規模スケールでの信号を維持するデータ処理手法が求められていました。

2. 手法とデータ処理

本研究では、HETDEX の内部データリリース（HDR5）およびテキサス・ユークリッド・サーベイ（TESLA）のデータを使用し、以下の厳密な処理ステップを踏みました。

2.1 データ選定とマスキング

• LAE カタログ: 信号対雑音比（SNR）> 5.5 の高信頼度 LAE を選別し、機械学習（Random Forest, CNN）による分類精度を確保。[O III] 線や天体発光線との混同を避けるための波長カットを実施。
• 強度スペクトルの前処理:
  • 天体発光除去: 「フルフレーム」のスカイサブトラクションを採用し、各ファイバーごとの系統的残差を除去。
  • 異常値除去: 連続スペクトルを持つ天体（前景星・銀河）、明るい天体、不良検出器、および [O II] 放射銀河を含むスペクトルをマスキング。
  • PCA による系統除去: 主成分分析（PCA）を用いて、天体発光の残差や検出器の系統的ノイズを除去。この際、信号の損失を最小限に抑えるため、最適な主成分数（$N_{PC}$）を決定（低赤方偏移で 170、中赤方偏移で 130、高赤方偏移で 100）。

2.2 マップ作成とシミュレーション

• マップ生成: 観測データを赤方偏移ビン（$\bar{z} = 2.1, 2.6, 3.2$）および空間領域に分割し、3 次元ボクセル（$2 h^{-1}\text{Mpc}$）に集約して強度マップと銀河過密度マップを生成。
• モックデータ生成: 「Simple」コードを用いて、対数正規分布（Lognormal）に基づくモックデータを生成。これには、HETDEX の感度、選択関数、マスク、赤方偏移空間歪み（RSD）、および PCA 処理による信号減衰効果が含まれています。
• 共分散行列の推定: データのシャッフル（位置のランダム化）による統計的誤差と、モックデータによる宇宙論的変動（Cosmic variance）の両方を考慮した共分散行列を構築。

2.3 パワースペクトル推定

• 検出された LAE の過密度（$\delta_g$）と、検出された LAE をマスキングした Lyα 強度の揺らぎ（$\delta_I$）との交差パワースペクトルを FFT を用いて計算。
• モデルとデータの比較により、パラメータ $A = b_g b_I \langle I \rangle \bar{F}_{RSD} / \bar{F}^{fid}_{RSD}$ をフィッティング。

3. 主要な貢献

1. HETDEX における初の高統計的検出: 検出された LAE と未検出源からの Lyα 強度との交差パワースペクトルを、$k \simeq 0.08 - 1 \, h\text{Mpc}^{-1}$ の範囲で高有意性で検出しました。
2. QSO 依存性の排除: 従来の QSO-Lyα 相関とは異なり、星形成銀河（LAE）を tracer として使用することで、QSO 周辺の局所的な放射場の影響を受けにくい、より普遍的な宇宙の Lyα 放射場を制限しました。
3. 高度な系統誤差制御: 大規模な観測データにおける天体発光残差や検出器ノイズを、PCA と厳密なマスキングを組み合わせることで効果的に除去し、信号を抽出する手法を実証しました。
4. 物理的解釈の枠組み: 測定された強度を、銀河形成シミュレーション（Byrohl & Nelson 2023）や光度関数の外挿と比較可能な形に変換し、Lyα 放射の起源（星形成、CGM、IGM）についての制約を提示しました。

4. 結果

• 測定値: 3 つの赤方偏移ビン（$\bar{z} = 2.1, 2.6, 3.2$）において、以下の積を制限しました：
  $$b_g b_I \langle I \rangle \frac{\bar{F}_{RSD}}{\bar{F}^{fid}_{RSD}} = (6.7 \pm 3.1), (11.7 \pm 1.4), (8.3 \pm 1.5) \times 10^{-22} \, \text{erg s}^{-1} \text{cm}^{-2} \text{arcsec}^{-2} \text{\AA}^{-1}$$
  （ここで $b_g$ は LAE バイアス、$b_I$ は強度バイアス、$\langle I \rangle$ は未検出源の平均強度）。
• 他研究との比較:
  • 従来の QSO-Lyα 交差相関の結果（Croft et al., Lin et al.）よりも有意に小さい値を示しました。これは、QSO 周辺の過剰な Lyα 放射が以前の測定を過大評価させていた可能性を示唆しています。
  • Lyα 森林からの上限値とは概ね整合的ですが、中赤方偏移ビンではわずかに高い値を示しています。
• 物理的解釈:
  • 仮定（$b_g=2, b_I=2$）のもと、測定された光度密度は、光度関数の外挿による予測値と同等かやや大きい値です。
  • 宇宙論的シミュレーション（Byrohl & Nelson 2023）と比較すると、測定値は CGM や IGM からの散乱光子を含む総強度のオーダーと一致します。
  • 星形成率密度から推定される Lyα 光度と比較すると、銀河からの Lyα 光子の逃げ率（escape fraction）$f_{esc} < 1$ であることが支持されます。

5. 意義と将来展望

• 宇宙論的意義: この結果は、HETDEX が Lyα 強度マッピングにおいて強力な統計的パワーを持つことを実証しました。また、現代の銀河形成・進化シミュレーションにおける Lyα 放射モデルを制約する重要な観測データとなります。
• 理論的整合性: 測定値は、Lyα 放射が主に星形成銀河に由来し、CGM や IGM での散乱を経て観測されるというシミュレーションの予測と矛盾しません。
• 今後の展望:
  • データ処理の改善（より高度なスカイサブトラクション、重み付けの最適化）により、信号対雑音比（SNR）をさらに向上させる余地があります。
  • 四極子パワースペクトルの解析や、Lyα 吸収（吸収トローグ）の影響をより詳細にモデル化することで、赤方偏移空間歪みや Lyα 放射輸送（RT）効果のより精密な測定が可能になると期待されます。

総じて、本論文は HETDEX データを用いた Lyα 強度マッピングの最初の成功例であり、高赤方偏移宇宙における未検出銀河と拡散ガスの性質を解明するための重要な基盤を築きました。