Measuring local primordial non-Gaussianity from the clustering of DESI DR1 LRGs and QSOs

技術要約：DESI DR1 クラスターリングによる局所的な原始非ガウス性の測定

問題提起
初期宇宙におけるインフレーション期の性質は、依然として宇宙論における主要な未解決問題である。宇宙マイクロ波背景放射（CMB）は、これまでで最も高感度な原始非ガウス性（PNG）の制約を提供してきたが、大規模構造（LSS）サーベイは、補完的かつ潜在的に競争力のあるプローブを提供する。パラメータ $f_{\text{NL}}$ で表される局所的な非線形性（PNG）は、トレーサーのクラスター化にスケール依存のバイアスを誘起する。しかし、構成空間（二点相関関数、すなわち2pcf）における $f_{\text{NL}}$ の測定は、イメージング・システマティクスによる汚染のため、極めて困難であることが知られている。サーベイの深さの変動、ダスト減光、および恒星密度の変化は、PNGの広域的なクラスター化信号を模倣する角密度ゆらぎを誘起し得る。本論文は、DESI DR1（Dark Energy Spectroscopic Instrument Data Release 1）から堅牢な $f_{\text{NL}}$ 測定値を抽出する際の課題に対処し、これらの系統的な影響を厳密に軽減する方法に取り組むものである。

手法
著者らは、DESI DR1 サンプルからの赤色銀河（LRG）およびクエーサー（QSO）の2pcf モノポールをモデル化するために、シミュレーションベースの統計フレームワークを採用している。

• データ処理: 解析には、2pcf の推定に ConKer アルゴリズムを使用する。トレーサーは、イメージング・システマティクス（$w_{\text{sys}}$）、カタログの不完全性（$w_{\text{comp}}$）、赤偏移失敗（$w_{\text{zfail}}$）、および FKP 重み（$w_{\text{FKP}}$）を考慮した重み（$w_{\text{tot}}$）を持つ3次元デカルト格子にマッピングされる。解析は、LRG（$0.4 < z < 1.1$）および QSO（$0.8 < z < 3.1$）について、北半球および南半球の銀河キャップ（NGC/SGC）を対象としている。
• 理論モデル: モデルは完全にシミュレーションベースであり、標準的なケース（$f_{\text{NL}}=0$）については AbacusSummit 高解像度 N体シミュレーションに基づき、PNG 依存性については「FastPM」シミュレーションに基づいている。2pcf は、$f_{\text{NL}}=0$ と $f_{\text{NL}}=100$ の実現例の間の線形補間としてモデル化される。
• バイアスと普遍性: モデルには線形銀河バイアス（$\tilde{b}_1$）と PNG バイアス（$b_\phi$）が含まれる。両者の関係は、普遍性関係 $b_\phi = 2\delta_c(b_1 - p)$ によって支配される。ここで $p$ はトレーサー依存のパラメータである。解析では、$p_{\text{LRG}} = 1.0 \pm 0.1$ および $p_{\text{QSO}} = 1.6 \pm 0.1$ というガウス事前分布を課している。
• システマティクス緩和: 手法の重要な構成要素は、イメージング・システマティクスの明示的なモデリングである。著者らは、特定のサーベイ領域（MzLS, (N)DEC, (S)DEC）においてイメージング重みの過剰または過少適用を表すニューストランス・パラメータ（$K^R_{\text{sys}}$）を導入している。これらのパラメータは、重みの変化が 2pcf に与える影響を事前に計算し、二次関数にフィットさせるシミュレーションベースのアプローチによって制約されている。
• 統計的推論: 1000 個の EZ（有効ゼルドビッチ）モックから導出された共分散行列を含む $\chi^2$ 統計量を用いた尤度関数を使用する。関心対象のパラメータ（$f_{\text{NL}}$ および線形バイアス $b_1$）と 7 つのニューストランス・パラメータは、マルコフ連鎖モンテカルロ法（MCMC）を用いてサンプリングされる。非線形モデリングのバイアスを避けるため、バリオン音響振動（BAO）スケール（$100 < s < 120 \, h^{-1}\text{Mpc}$）はマスクされている。

主要な結果
本解析は、DESI データを用いた構成空間 2pcf における初の $f_{\text{NL}}$ 測定結果をもたらした。

• LRG サンプル: 測定値は $f_{\text{NL}} = -3^{+22}_{-21}$（68% 信頼区間 CL）であった。
• QSO サンプル: 測定値は $f_{\text{NL}} = 0^{+17}_{-16}$（68% CL）であった。
• 結合制約: 両方のトレーサーサンプルを組み合わせた結合制約は、$f_{\text{NL}} = -3^{+12}_{-12}$（68% CL）となった。
• バイアス・パラメータ: 線形バイアスは、LRG に対して $b_1 = 2.01 \pm 0.06$、QSO に対して $b_1 = 2.57 \pm 0.10$ と制約された。
• システマティクス: イメージング・システマティクスを制御するニューストランス・パラメータ（$K_{\text{sys}}$）および形成パラメータ $p$ は、概ね事前分布と一致しており、データが標準的なシステマティクス・モデルからの偏差を強く支持していないことを示している。

意義と主張
本論文は、主に 3 つの結論を主張している。

1. 構成空間解析の堅牢性: 本研究は、堅牢な PNG 研究にはフーリエ空間と構成空間の両方での解析が必要であることを実証している。PNG の影響はこれら両方の領域で異なるが、構成空間の結果が以前のフーリエ空間解析（Chaussidon et al. 2025）と一致していることは、このアプローチの妥当性を裏付けている。
2. システマティクス緩和の重要性: 本研究は、LSS による PNG 測定にはイメージング・システマティクスの制御が極めて重要であることを強調している。この構成空間の結果と、異なるシステマティクス緩和戦略を用いる独立したフーリエ空間の結果との一致は、サーベイの不完全性や角密度変動を処理するパイプラインの能力に対する強力な検証となっている。
3. LSS の制約能力: 最も重要な主張は、LSS サーベイが CMB プローブの感度に急速に近づいていることである。わずか 1 年目の DESI データのみで、結合制約（$\sigma_{f_{\text{NL}}} \approx 12$）は、従来の LSS 制約を大幅に上回り、Planck CMB の感度（$\sigma_{f_{\text{NL}}} \approx 5$）に接近している。著者らは、この分野が、後期宇宙の構造特性化が原始物理学の主要なプローブとなる時代に入ったと断言している。

結果は標準的なインフレーションの「スローロール」モデル（$f_{\text{NL}} \approx 0$）と一致しており、現在の精度限界内ではマルチフィールド・インフレーションのシナリオを示す証拠は見られない。