Tracing Large-scale Structure with the MeerKLASS On-the-Fly Survey: Angular Clustering of Radio Sources at 816 MHz

技術要約：MeerKLASS On-the-Fly サーベイによる大規模構造の追跡

問題と動機
電波連続波サーベイは、宇宙の大規模構造（LSS）に対する、塵による減光の影響を受けない独自の窓を提供しており、電波放射の強い活動銀河核（AGN）や星形成銀河（SFG）を通じて物質密度場を追跡する。しかし、ラジオカタログからのロバストな角相関測定は、位置依存のノイズによる空間的に変化する完備性、複雑なソース形態（単一の物理的銀河が複数のカタログコンポーネントに分解される現象）、および赤方偏移情報の直接的な欠如により、容易ではない。FIRSTやNVSSなどのサーベイを用いたこれまでの測定では、明るいAGN主体の集団のクラスタリングが確立されているが、特性がよく把握された深さ変動を持ち、現代の光学分光プログラムと重複する現代的な干渉計を用いた、中間周波数（$\sim$GHz）における高フィデリティなクラスタリング測定へのニーズが依然として存在する。本論文は、MeerKAT Large Area Synoptic Survey (MeerKLASS) UHF on-the-fly (OTF) continuum Data Release 1 (DR1) の最初の角相関解析を提示することで、これらの課題に対処する。

手法
解析には、参照周波数816 MHzにおいて約800 deg$^2$の範囲をカバーするMeerKLASS UHF DR1データセットを使用する。本サーベイは、交差リンクされた昇角・降角スキャンを用いる一定高度OTF戦略を採用しており、その結果として生じる非一様な深さパターンを明示的にモデル化する必要がある。

1. サンプル定義: DR1ソースカタログ（SRL）に含まれる95,483個のユニークなソースから、保守的なクラスタリングサンプルを構築した。標準的な選択基準は以下の通りである：

   • フラックス密度閾値：$S_{816} > 2$ mJy。
   • 深さ制御マスク：局所的なRMSノイズ $\sigma_{local} \le 250$ $\mu$Jy beam$^{-1}$の領域に解析を制限。
   • 「コンパクトのみ」基準（$S\_Code = 'S'$ かつ $N\_Gaus = 1$）：分割されたマルチコンポーネント電波銀河による小スケールの過剰パワーを軽減するため。
   • 保守的な赤経（RA）エッジカット（$RA \le 184.07^\circ$）：データとランダムカタログ間の一様性を確保するため。
     これにより、最終的なクラスタリングカタログは39,483個のソースとなった。

2. マスクとランダムカタログ: 幾何学的フットプリントとRMS深さマップからバイナリサーベイマスクを構築した。ランダムカタログは、このマスク内で一様に生成され、データの対ランダムの表面密度比を確認することによって、角選択関数を追跡するように検証された。

3. 相関測定: 角二点相関関数 $w(\theta)$ は、20個の対数間隔ビン（$0.02^\circ \le \theta \le 10^\circ$）を用いて、Landy–Szalay推定量により測定された。不確実性とビン間の共分散は、30個のパッチを用いた空間的ジャックナイフ再サンプリング手法を用いて推定された。また、有限なサーベイフットプリントを考慮するために、積分制約（IC）補正が適用された。

4. モデリングとバイアス推論:

   • 現象論的フィット: 信号は、固定された傾き $\gamma=1.8$ を持つべき乗則モデル $w(\theta) \approx A(\theta/1^\circ)^{1-\gamma}$ を用いて特徴付けられた。
   • 物理的解釈: 測定されたクラスタリングは、CAMBとリンビ近似（Limber approximation）を用いて計算された、$\Lambda$CDM投影物質テンプレート（$w_{DM}(\theta)$）へのフィッティングによって解釈された。個々の赤方偏移は利用できないため、解析ではT-RECS半経験的シミュレーションから導出された2つの事前分布（AGN主体の事前分布と、TOTAL（AGN + SFG）事前分布）を用いて、主要な系統的不確定性（赤方偏移分布 $N(z)$）を挟み撃ちにした。これにより、テンプレートをデータにスケーリングすることで、有効な大規模バイアス（$b_{eff}$）が推論された。

主な結果

• 検出: $0.02^\circ \lesssim \theta \lesssim 10^\circ$ の範囲で正のクラスタリング信号が検出された。測定値は、中間スケールおよび大規模スケールにおいて、深さマスクおよびフラックス閾値の変化に対して安定している。
• べき乗則振幅: 固定された傾き $\gamma=1.8$ で中間スケールの信号（$0.112^\circ \le \theta \le 1.36^\circ$）をフィッティングすると、振幅は $A(1^\circ) = (1.434 \pm 0.475) \times 10^{-3}$ となり、これは $\log_{10} A = -2.843^{+0.124}_{-0.175}$ に相当する。
• 有効バイアス: 推論された有効バイアスは、仮定された赤方偏移分布に大きく依存する：
  • AGN事前分布：$b_{eff} = 1.998 \pm 0.350$。
  • TOTAL事前分布：$b_{eff} = 1.530 \pm 0.265$。
    これらの値の差（$\Delta b_{eff} \approx 0.47$）は、未知の集団混合に起因する主要なモデリングの不確定性として特定されている。
• 相関長: べき乗則の振幅をリンビ反転（Limber-inverting）することにより、空間相関長は $r_0 = 6.18 \pm 1.13$ $h^{-1}$Mpc（AGN事前分布）および $r_0 = 5.59 \pm 1.02$ $h^{-1}$Mpc（TOTAL事前分布）となる。
• 系統誤差: コンパクト性の基準を緩和すると、小スケール（$\theta < 0.1^\circ$）で $w(\theta)$ の顕著な上昇が見られた。これは、ロバストな大規模クラスタリング測定のために「コンパクトのみ」の選択が必要であることを裏付けている。

意義と主張
本論文は、MeerKLASS UHF DR1が、電波連続波のクラスタリングにおける新しい、広域かつ中間周波数のデータセットであることを確立している。著者らは、自らの結果が、古典的な全天GHzサーベイ（例：NVSS）と新世代の低周波サーベイ（例：LoTSS）の間のギャップを埋める、精密な大規模構造測定のためのOTF干渉計連続波イメージングの実行可能性を示していると主張している。

本研究は、現在の結果が赤方偏移分布（$N(z)$）の不確実性に制限されているものの、DESIフットプリントとの重複が明確な道筋を提供していることを強調している。著者らは、改良された赤方偏移キャリブレーションのために重複する光学/赤外分光を取り入れた将来のリリースによって、集団分離クラスタリングおよびバイアス進化の測定が可能になり、モデリングの不確定性が大幅に減少すると述べている。本研究は基礎的なステップとして機能し、宇宙論的応用に対するサーベイ戦略とデータプロダクトを検証する、保守的で高度に制御された測定を提供している。