Improving constraints on primordial non-Gaussianity from Quaia with a new cosmological observable: angular redshift fluctuations

技術要約：角度赤方偏移ゆらぎを用いたQuaiaによる原始非ガウス性の制約の改善

問題と動機
$f_{\text{NL}}$ でパラメータ化される原始非ガウス性（PNG）の検出は、インフレーションモデルに対する重要なテストであり、非ゼロの値はガウス的な初期条件からの逸脱を示す。現在、最もタイトな制約はCMBバイスペクトル（例：Planck 2018）から得られているが、大規模構造（LSS）はスケール依存バイアス効果を通じて補完的なプローブを提供できる。この効果は、Dalal et al. (2008) によって予測されており、局所的なPNGに対して特徴的な $1/k^2$ のシグネチャを局所的な銀河バイアスに導入する。しかし、LSSからこれを測定することは、最大スケールにおける観測的な系統誤差を精密に制御する必要があるため、困難を伴う。

近年の研究では、2次元角度パワースペクトルおよびCMBレンズ効果との相互相関を利用して $f_{\text{NL}}$ を制約している。本研究は、新しい宇宙論的観測量である「角度赤方偏移ゆらぎ（ARF）」を組み込むことで、これらの制約を改善できる可能性に取り組んでいる。ARFは、特定のシェル内における平均的な赤方偏移の2次元的な角度偏差を測定するものである。理論的に、ARFは銀河バイアスに敏感であり、したがってPNGにも敏感であるが、実データへの $f_{\text{NL}}$ 制約への適用はこれまで実証されていなかった。

手法
著者らは、Gaia DR3およびunWISEデータから得られた、広大な宇宙論的体積をカバーする約130万個のクエーサー（QSO）からなる全天サンプルであるQuaiaカタログを分析した。解析パイプラインは以下のステップを含む：

1. データ準備： Quaiaサンプルを2つの赤方偏移ビン（$\bar{z} \approx 0.97$ および $\bar{z} \approx 2.10$）に分割した。銀河数密度（密度）とARFの両方について、HEALPixマップ（$N_{\text{side}}=128$）を生成した。ARF場は、サーベイ固有の赤方偏移分布に相当する広いガウス窓を用いた、ピクセルあたりの加重平均赤方偏移ゆらぎとして定義された。
2. 相互相関： 研究では、Quaiaの密度マップとARFマップをPlanck PR4 CMBレンズ収束マップと組み合わせた。これにより、4つの観測量（密度とARFの自己相関スペクトル $C_{gg}^\ell$ および $C_{zz}^\ell$、およびCMBレンズとの相互相関 $C_{\kappa g}^\ell$ および $C_{\kappa z}^\ell$）が得られた。
3. 理論的モデリング： 理論的なパワースペクトルを計算するために、修正版のCAMB Boltzmannコード（ARFCAMB）を使用した。このコードは、ARFカーネルを組み込み、光度赤方偏移の不確かさ（$\sigma_z \approx 0.06(1+z)$）を考慮している。モデルには、普遍性関係（$p=1$ を仮定）の下で誘起されるスケール依存バイアスと、銀河バイアスの赤方偏移進化に関する柔軟なパラメータ化（$b(z)$）が含まれている。
4. 系統誤差と共分散： 観測的な系統誤差に対処するため、著者らはQuaiaテンプレートに対する線形回帰を用いた系統誤差デプロジェクションを採用した。統計的なゆらぎと矛盾する低$\ell$オフセットを除去するために、データシフトを模擬実現（モック）と比較することで、スケールカット（具体的には密度に対して $\ell_{\text{min}}=14$、ARF自己スペクトルに対して $\ell_{\text{min}}=10$）を適用した。共分散行列は、1,000個の相関したシミュレーションから導出され、実データで見られる低マルチポールでの過剰なパワーに一致するように、モック自己スペクトルに補正が適用された。
5. パラメータ推論： ガウス尤度を構築し、測定されたパワースペクトルと共分散行列を用いた。MCMCサンプリング（emcee）を用いて、$f_{\text{NL}}$ およびバイアスパラメータ（$b_0, \alpha$）を制約した。

主な結果
解析により、局所的な非ガウス性パラメータ $f_{\text{NL}}$（$p=1$ の普遍性関係を仮定）に関する以下の制約が得られた：

• ARF単独： ARF自己相関とARF-CMBレンズ相互相関を組み合わせると、$f_{\text{NL}} = -43^{+50}_{-44}$ となった。
• 密度単独： 密度自己相関と密度-CMBレンズ相互相関を組み合わせると、$f_{\text{NL}} = -25^{+19}_{-17}$ となり、以前のQuaia解析（Fabbian et al. 2026）と一致した。
• 結合解析： 4つの観測量すべて（$C_{gg}^\ell, C_{\kappa g}^\ell, C_{zz}^\ell, C_{\kappa z}^\ell$）を組み合わせた結果、$f_{\text{NL}} = -3 \pm 14$（68%信頼区間）となった。
• 適合度の最適化： 著者らは、ARF-CMB相互相関（$C_{\kappa z}^\ell$）を含めると $\chi^2$ に軽微なテンションが生じることを指摘した。モデルの適合度を改善するためにこの特定の相互相関を除外した結果、$f_{\text{NL}} = -5^{+16}_{-15}$ となった。

この結合制約は、密度とCMBレンズのみを用いた以前のQuaiaの測定値と比較して、約25%の改善を示している。この結果は、LSSの2点統計を用いた $f_{\text{NL}}$ の制約としては、DESI DR1 分光データに次いで2番目にタイトな制約であり、投影された2次元要約統計量を用いて達成されたものとしては最も精密である。

有意性と頑健性
論文は、ARFを宇宙論的観測量として組み込むことが、実データを用いた $f_{\text{NL}}$ 制約の測定可能な改善を提供する、と主張している。著者らは、以下のテストを通じて結果の頑健性を検証した：

• バイアスモデリング： 結果は、基準となる銀河バイアス進化モデルを変えること（Piccirilli et al. 2024 と Laurent et al. 2017 の両方のモデルをテスト）、および密度場とARF場に対して独立したバイアスパラメータを許容することによっても安定していた。
• スケールカット： 制約は最小マルチポール（$\ell_{\text{min}}$）カットの変動に対して安定していることが示されたが、バイアスパラメータについては、ARFカーネルが密度カーネルよりも非線形スケール（$k$モード）に対して高い感度を持つことに起因する、ある程度のテンションが観察された。
• 普遍性関係： 普遍性関係の仮定を $p=1$ から $p=1.6$（一部の文献でQSOに対して示唆されている値）へと緩和すると、制約は $f_{\text{NL}} = 2^{+20}_{-21}$、あるいはモデルに依存しないアプローチでは $b_\phi f_{\text{NL}} = -63^{+79}_{-89}$ へとシフトした。

著者らは、観測量間に内部的なテンションが存在するものの（おそらくARFにおける非線形スケールのモデリング向上の必要性に起因する）、2次元解析にARFを追加することは有望な道であると結論付けている。彼らは、DESI、Euclid、LSSTといった将来の大型サーベイの解析にARFを含めることを推奨しており、ARFが、赤方偏移の不確かさが無視できないフォトメトリックサーベイにおいて、標準的な3次元クラスタリングに代わる堅牢な2次元の選択肢となり得ることを述べている。