Beyond power spectrum to unveil systematics on HI intensity maps

技術的概要：Hi 強度マップにおける系統誤差の解明に向けたパワースペクトルを超えたアプローチ

問題提起
中性水素（Hi）強度マッピングは、特に非線形重力進化により物質分布が強く非ガウス的となる中間赤方偏移（$z < 1$）領域における大規模構造の探査における主要な手法である。従来、これらのサーベイからの宇宙論的情報は、ガウス的特性を効果的に捉える一方で高次相関を棄却する 2 点統計量である角度パワースペクトル（$C_\ell$）を用いて抽出されてきた。しかし、$z \approx 0.425$（MeerKLASS L バンド観測のターゲット）のような赤方偏移では、Hi 信号は $C_\ell$ では完全に特徴付けられない非線形構造（フィラメント、クラスター、ボイド）に支配されている。さらに、望遠鏡ビームの平滑化、部分的な天域カバレッジ（マスキング）、機器雑音、前景残差（銀河シンクロトロン、銀河系外点源）といった観測上の系統誤差は、標準的な解析にバイアスを及ぼしうる。本論文は、非ガウス情報を捉え、かつこれらの系統誤差に対して頑健であり、宇宙論パラメータ推論を改善するための要約統計量の必要性に取り組む。

手法
著者らは、JaxILI ライブラリを用いたシミュレーションベース推論（SBI）フレームワークを採用し、尤度なしベイズ推論を実行した。

• シミュレーション: 赤方偏移 $z = 0.425$ において、GLASS および CAMB フレームワークを用いて、10,000 枚の全天 Hi 輝度温度マップを生成した。これらのマップは、赤方偏移依存の Hi バイアスと平均輝度温度で重み付けされた対数正規物質場をモデル化している。
• 観測効果: シミュレーションには現実的な観測条件が組み込まれている。
  • ビーム: ガウスビームとの畳み込み（単一パラボラアンテナ型の MeerKLASS 類似で 1.34°、将来の干渉計構成で 0.5°）。
  • マスキング: 天域割合（$f_{sky}$）が 0.60、0.19、および 0.014 を表す二値マスクの適用。
  • 雑音: システム温度と積分時間に基づいた無相関ガウス熱雑音の付加。
  • 系統誤差: 乗法的および加法的なストライプパターン（走査アーティファクト）、銀河シンクロトロン残差、銀河系外点源残差、および X 字型走査パターンの制御された玩具モデル。
• 要約統計量: 2 つの統計量が抽出され比較された。
  1. 角度パワースペクトル（$C_\ell$）: 球面調和関数係数から計算され、系統誤差を緩和するため低次多重極（$\ell < 10$）を除外。
  2. スターレット $\ell_1$ ノルム: 等方性アンデシメートウェーブレット変換（スターレット）に基づくマルチスケール、高次統計量。マップは複数のスケール（$J=5$ および $J=8$）にわたるウェーブレット係数に分解され、正規化され、信号対雑音比空間でビン分けされる。これらの係数の $\ell_1$ ノルムは、非ガウス的特徴の疎性と強度を定量化する。
• 推論: 条件付きニューラル密度推定量（マスク自己回帰フロー）が、シミュレーションペア（パラメータ $\theta = \{\Omega_c, h, A_s\}$、要約統計量 $x$）に基づいて訓練され、事後分布 $p(\theta|x)$ を近似した。

主要な貢献

1. スターレット $\ell_1$ ノルムの拡張: 以前は弱い重力レンズで成功を収めていたスターレット $\ell_1$ ノルムを Hi 強度マッピングに適応させ、21cm 輝度温度場における非ガウス情報を捉えるその有効性を実証した。
2. 系統誤差に対する頑健性分析: 本研究は、ビーム、マスキング、雑音、ストライプ、および前景残差といった様々な系統誤差下において、$C_\ell$ とスターレット $\ell_1$ ノルムの両方がどのように機能するかを包括的に評価した。
3. 診断能力: 著者らは、スターレット $\ell_1$ ノルムが異なる種類の前景残差（銀河シンクロトロン対銀河系外点源）を区別し、残留汚染の存在を検出する診断ツールとして機能しうることを実証した。
4. 効率性: 従来の MCMC 法に対する、コンパクトな要約統計量（スターレット $\ell_1$ ノルムは完全な $C_\ell$ ベクトルよりもはるかに少ないデータポイントを使用する）を用いた SBI の計算的および環境的な利点を強調した。

結果

• 宇宙論的制約: スターレット $\ell_1$ ノルムは、宇宙論パラメータ（$\Omega_c$、$h$、$A_s$）の制約において、一貫して角度パワースペクトルを上回る性能を示した。
  • 雑音のない全天シナリオでは、スターレット $\ell_1$ ノルムは $C_\ell$ に対して図の指標（FoM）を 1.59 倍改善した。
  • 0.5° ビームと機器雑音を伴う現実的なシナリオでは、スターレットのスケール数を 5 から 8 に増加させることで、$C_\ell$ に対して FoM が 2.81 倍改善された。
• 系統誤差に対する頑健性:
  • ビームとマスキング: ビームの平滑化は小スケールのパワーを抑制することで $C_\ell$ を著しく劣化させるが、スターレット $\ell_1$ ノルムはスケールにわたる構造コントラストに対する感度を維持し、頑健である。また、$C_\ell$ がモード結合に苦しむ強いマスキング（$f_{sky}=0.014$）下でも性能を維持する。
  • ストライプと走査: スターレット $\ell_1$ ノルムは、乗法的および加法的なストライプパターンや X 字型走査アーティファクトに対してほとんど感度を示さないのに対し、$C_\ell$ は強いバイアスと事後分布の人工的な引き締めを示す。
  • 前景残差: スターレット $\ell_1$ ノルムは、$C_\ell$ よりも残留汚染に対して高い感度を示す。これは高汚染レベルで軽度のバイアスを導入するが、シンクロトロンと点源残差を区別することを可能にし、$C_\ell$ が失敗するタスクを達成する。
• パラメータの縮退: スターレットアプローチは、非ガウス情報を利用することで、パラメータ間の縮退（特に $C_\ell$ によって弱く制約される $h$）を解消する。

意義
本論文は、マルチスケールの高次統計量であるスターレット $\ell_1$ ノルムが、Hi 強度マッピングにおける従来の 2 点解析にとって不可欠な補完であることを結論付けている。パワースペクトルによって見逃される非ガウス的特徴を捉えることで、スターレット $\ell_1$ ノルムは、特に非線形領域（$z < 1$）において宇宙論的推論を著しく強化する。複雑な観測系統誤差に対する実証された頑健性と、前景残差の診断ツールとして機能する能力は、現在の MeerKLASS などのサーベイや、将来の SKA 観測所（SKAO）のような大規模サーベイにとって有望なツールである。この研究は、今後の Hi 強度マッピングデータの科学的成果を最大化するために、そのような統計量を用いたシミュレーションベース推論の採用を提唱している。