EMU: Cross-correlating EMU Pilot Survey 1 with Dark Energy Survey to validate the radio galaxy bias and redshift distribution

技術要約：ラジオ銀河のバイアスと赤方偏移分布を検証するための、DESを用いたEMUの相互相関

問題提起
電波連続波サーベイは、広大な空領域をカバーできる能力とダスト消散に対する無感度เนื่องจาก、宇宙の大規模構造（LSS）をマッピングするための強力な手法を提供します。しかし、ラジオサーベイを用いた宇宙論的解析における決定的な制限は、ラジオ源に対する直接的な赤方偏移情報の欠如です。支配的な放射メカニズムであるシンクロトロン放射は、特徴のない冪乗則スペクトルを生成するため、赤方偏移の推定は極めて困難です。正確な赤方偏移分布（$dN/dz$）とロバストな銀河バイアスモデルがなければ、$\sigma_8$などの宇宙論パラメータの制約に系統誤差が伝播する可能性があります。これまでの研究では、光学サーベイとのクロスマッチング（ソースのわずかな一部に対してのみフォトメトリック赤方偏移が得られる）や、シミュレーション（TRECS、SKADSなど）に依存してきましたが、個々のソースの特定に頼ることなく、サーベイ領域全体を利用して統計的に赤方偏移分布を推論する手法が求められています。

手法
本研究では、Evolutionary Map of the Universe Pilot Survey 1 (EMU-PS1) の二次元電波連続波マップと、Dark Energy Survey (DES) のフォトメトリック光学銀河サンプルとの間の相互相関技術を採用しています。解析は以下の手順で行われます。

1. データ準備:

   • EMU-PS1: PyBDSFソースファインダーによって生成された、約224,100個のソースを含むカタログを使用します。解析には、保守的なサンプルを作成するために180 $\mu$Jyおよび400 $\mu$Jyのフラックス密度カットを適用します。観測ノイズやソース数の変動を考慮するために、ウェイトマップを構築します。
   • DES: 6つのトモグラフィック赤方偏移ビン（$0.20 \le z < 1.05$）に分割されたマグニチュード制限（MagLim）銀河サンプルを使用します。フォトメトリック赤方偏移は、Directional Neighbourhood Fitting (DNF) アルゴリズムを用いて導出されます。
   • Planck: ラジオバイアスと赤方偏移分布の独立した検証を提供するために、Planck CMB レンジング収束マップ（$\kappa$）を含めます。

2. パワースペクトル推定:

   • 不完全な空の被覆率を扱うために、擬似-$C_\ell$エスティメータ（NaMaster）を使用して角度パワースペクトル（$C_\ell$）を計算します。
   • 解析では、EMU-PS1の自己パワースペクトル（$C^{EMU}_\ell$）、DESビンの自己パワースペクトル（$C^{DES}_\ell$）、EMU-PS1とDESビンの間の相互パワースペクトル（$C^{EMU \times DES}_\ell$）、およびEMU-PS1とPlanck レンジングの間の相互パワースペクトル（$C^{EMU \times \kappa}_\ell$）を測定します。
   • 共分散行列は、サーベイのフットプリントとポアソンノイズを組み込んだ、GLASSによって生成された1,000個の模擬実現体を用いて推定されます。

3. モデリングと推論:

   • 赤方偏移分布 ($n(z)$): 2つのモデルをテストします。第一に、活動銀河核（AGN）と星形成銀河（SFG）の寄与を含む、Tiered Radio Extragalactic Continuum Simulation (TRECS) による赤方偏移分布です。第二に、分布のピークと裾野を物理的に動機付けるための、3つの自由パラメータ（$z_0, r, \gamma$）を持つパラメトリックモデル（Nakoneczny et al. 2024を適応）です。
   • 銀河バイアス: 線形かつスケール独立なバイアスモデルを採用します。パラメトリックな $n(z)$ に対しては、SFGには一定振幅のバイアスモデルを、AGNには多項式進化（または結合された有効バイアス）を使用します。
   • パラメータ推定: 相互パワースペクトルの結合データベクトルにフィットさせるために、最大尤度フレームワーク（PyMCを使用）を用いて、バイアスパラメータと赤方偏移分布のパラメータを制約します。

4. 検証:

   • 真の $n(z)$ とバイアスが既知である模擬カタログを用いて、解析パイプラインを検証します。これにより、赤方偏移分布のピーク位置を回収する本手法の能力を確認します。
   • ラジオカタログの残留する不完全性の影響を評価し、自己パワースペクトルは不完全性に敏感である一方、推論に使用される相互パワースペクトルは頑健であることを示します。

主な結果

• 検出有意性: EMU-PS1と6つのDESトモグラフィックビンの間の相互相関は、マルチポール範囲 $50 < \ell < 250$ において、$3.5\sigma$ から $5.5\sigma$ の有意性で検出されました。
• DESバイアスの検証: クラスタリング測定のみを用いて推定されたDESビンの銀河バイアスは、以前のDESの結果（例：Abbott et al. 2022）と$1\sigma$以内で一致しており、解析パイプラインの妥当性を検証しています。
• 赤方偏移分布の回収:
  • TRECSモデル: TRECSシミュレーションによる $n(z)$ は、銀河バイアスの混合モデルと組み合わせることで、相互パワースペクトル、EMU-PS1の自己パワースペクトル、およびCMB レンジングとの相互相関によく適合します。総 $\chi^2$ 値は、20自由度に対して、180 $\mu$Jyで7.78、400 $\mu$Jyで3.15です。
  • パラメトリックモデル: パラメトリックモデル（式11）も、AGN/SFGの混合に関する事前知識を必要とせずに、相互パワースペクトルへの良好なフィットを提供します。最良フィットの分布は $z \sim 0.85$ にピークを持ちます。このピークは、TRECSやLOFARディープフィールドの予測よりも高い値ですが、COSMOS 3GHzおよびMIGHTEEサーベイの結果とは一致しています。
• バイアス進化: 回収されたEMU-PS1の銀河バイアス進化は、ソース集団を明示的に分離していないパラメトリックモデルを使用した場合でも、TRECSの予測および文献からの独立した推定値の両方と一致しています。
• 系統誤差: パラメトリックモデルが180 $\mu$JyカットにおいてEMU-PS1の自己パワースペクトルへの適合に失敗しているという不一致が認められ、これは（潜在的に不完全性を示す）ラジオデータの残留系統誤差を示唆していますが、相互相関の結果には影響しません。

意義と主張
著者らは、本研究が、個々のソースの赤方偏移やシミュレーションのみに依存することなく、大規模領域の電波連続波サーベイから赤方偏移分布と銀河バイアスを直接推定できる可能性を示したものであると主張しています。

• 手法の検証: 相互相関技術は、赤方偏移分布の形状を正常に回収し、バイアスモデルを検証しました。これは、直接的な光学マッチングを補完するアプローチとなります。
• 宇宙論的有用性: 本結果は、将来の宇宙論的解析において、パラメトリック $n(z)$ モデルおよびTRECSシミュレーションの使用を正当化するものです。これは、フルEMUサーベイやスクエア・キロメートル・アレイ観測装置（SKAO）のような、今後の大規模サーベイにおいて重要となります。なぜなら、すべてのソースに対して直接的な赤方偏移測定を行うことは不可能だからです。
• 頑健性: 本手法は、相互パワースペクトルを使用する場合、ラジオカタログの内部系統誤差や残留不完全性に対して頑健であることが示されました。ただし、著者は、パイロット調査の小さなサーベイ面積（$\sim$270 deg$^2$）が、メインサーベイで達成可能となる精度と比較して、制約の精度を制限していることに注意を促しています。

結論として、光学サーベイとの相互相関は、精密宇宙論時代における電波天文学の実行可能かつ必要なツールであり、現在ラジオ連続波データの科学的リターンを制限している赤方偏移分布の不確実性を緩和することを可能にします。