Accurate modeling for 3$\times$2pt analyses in Roman and Rubin: a study of model approximations

技術要約：RomanおよびRubinにおける3×2pt解析のための正確なモデリング

問題提起
現代の宇宙論は、銀河のクラスタリング、銀河・銀河間重力レンズ、および宇宙シアーを組み合わせることで宇宙論パラメータを制約する「3×2pt」解析に大きく依存している。このフレームワークは、Stage-IIIサーベイ（DES、KiDSなど）では成功裏に適用されてきたが、次世代のStage-IVサーベイ、具体的にはNancy Grace Roman宇宙望遠鏡の高緯度撮像サーベイ（HLIS）およびVera C. Rubin天文台のLegacy Survey of Space and Time（LSST）は、前例のない統計的精度を提供することになる。著者らは、Stage-IIIのデータに対しては許容可能であったモデリングの近似が、Stage-IVの制約において重大なバイアスを導入する可能性があると主張している。本論文は、3×2pt解析で一般的に用いられる3つの特定のモデリング近似を系統的に調査している：

1. 3次元パワースペクトルを角スケールへ投影するためのLimber近似の適用。
2. 銀河過密度の場における**赤方偏移空間歪み（RSD）**の無視。
3. 非線形物質パワースペクトルに対する精度の低いモデルの使用（具体的には、HMCodeとEuclidEmulator2の比較）。

手法
本研究は、Roman HLIS Cosmology PITによって開発されたパイプラインであるCoCoA（CosmoLike-Cobaya Architecture）に実装された、シミュレーションされた尤度解析フレームワークを利用している。堅牢性を確保するため、著者らはまず、Rubin LSST Dark Energy Science Collaboration（DESC）の標準である**Core Cosmology Library (CCL)**に対してCoCoAの検証を行った。

• サーベイ構成： 解析は、5年間のRoman HLISデータセット（ミディアムティア）および10年間のRubin LSSTの制約を予測する。両者はともに10個のレンズ銀河ビンを使用し、ソース銀河の数は異なる（Romanは8、Rubinは5）。
• モデリングの枠組み： 解析は、5つの自由パラメータ（$A_s, n_s, H_0, \Omega_b, \Omega_m$）を持つ平坦な$\Lambda$CDM宇宙論を仮定している。周辺パラメータ（Nuisance parameters）には、photo-zのシフト、シアー校正バイアス、固有整列（NLAモデル）、および線形銀河バイアスが含まれる。
• スケールカット： バリオンフィードバックおよび非線形銀河バイアスの不確実性を軽減するため、著者らは保守的なスケールカットを適用している。実空間では、これは重力レンズに対して最小横半径6 Mpc/h、クラスタリングに対して8 Mpc/hに相当する。調和空間では、これらは特定の多重極範囲（$20 < \ell < 4000$）に変換される。
• 実験設計： 著者らは、近似の影響を分離するために5つのモデルシナリオを定義した：
  1. ベースライン： 非Limberによるクラスタリング、RSDを含む、EuclidEmulator2による非線形パワースペクトル。
  2. Limber： クラスタリングにLimber近似を使用。
  3. noRSD： RSD補正を省略。
  4. HMCode： EuclidEmulator2の代わりにHMCodeを使用。
  5. 結合（Combined）： Limber、RSDの省略、およびHMCodeを使用。

主な結果
本研究は、これらの近似がデータベクトル（$\Delta\chi^2$を介して）および結果としての宇宙論的制約（$\Omega_m - S_8$平面におけるバイアス）に与える影響を定量化している。

1. コードの検証： CoCoAとCCLの比較は、基準設定において$\Delta\chi^2$値が1を大幅に下回る優れた一致を示し、モデリングパイプラインの信頼性を確認した。
2. 個別の近似の影響：
   • Limber近似： 銀河クラスタリングにおける非Limber計算を無視することは、Romanで**0.7–0.8$\sigma$**、Rubinで**1.8–2.0$\sigma$**のバイアスを誘発する。この影響は、スケールカットによって調和空間における影響を受ける大きなスケールのモードが除去されるため、実空間よりも調和空間において顕著である。
   • 赤方偏移空間歪み（RSD）： RSDを無視すると、Limber近似と同等またはそれよりわずかに大きいバイアスが生じ、Romanで**0.7–0.8$\sigma$**、Rubinで**2.0–2.3$\sigma$**に達する。著者らは、Limber近似が本質的にRSDの寄与を抑制するため、RSDがこれら2つのうち支配的な効果であると述べている。
   • 非線形パワースペクトル： EuclidEmulator2からHMCodeへの切り替えは、より小さなバイアス（Romanで**0.6$\sigma$**、Rubinで**1.0–1.2$\sigma$**）をもたらす。これは、これらのモデル間の差異が主に小スケールにあり、それらが基準となるスケールカットによって大部分が除去されているためである。
3. 結合効果： 3つの近似すべてを同時に適用した場合、バイアスは累積する。Rubinの場合、結合されたバイアスは実空間で**2.8$\sigma$**、調和空間で**3.1$\sigma$に達する。Romanの場合、結合されたバイアスは~1.0–1.1$\sigma$**である。
4. サーベイの感度： バイアスは、Rubinの制約力がより高いことから、一貫してRomanよりもRubinの方が大きくなっている。著者らは、Romanについて「ミディアムティア」を想定したが、「ワイドティア」構成であれば、よりタイトな制約が得られ、近似を使用した場合のバイアスもさらに大きくなる可能性があると指摘している。

意義および結論
本論文は、Stage-III解析で使用されている標準的な近似は、Stage-IVサーベイには不十分であると結論付けている。調査された3つの効果のいずれかを無視することは、Romanに対しては1$\sigma$に近く、Rubinに対しては2$\sigma$に達する宇宙論的バイアスを誘発する可能性がある。

著者らは、将来のStage-IV解析に向けて、以下の3つの主要な示唆を強調している：

1. 非Limberアルゴリズム： 高精度サーベイにおいてLimber近似は銀河クラスタリングに対して不十分であるため、非Limberアルゴリズムを最適化し、数値的に安定かつ計算効率的にすることが極めて重要である。
2. RSDモデリング： 線形次数のRSD効果はデータベクトルと制約に大きな影響を与えるため、より現実的なRSDモデルの開発と導入が求められる。
3. 非線形パワースペクトル： 非線形物質パワースペクトルの選択は、特定のスケールカット戦略と組み合わせた場合に重要となる。

本研究は、CCLに対するCoCoAパイプラインの検証として機能し、RomanおよびRubinサーベイの最終データリリースにおいて、より高忠実度なモデリングを必要とする定量的根拠を提供するものである。著者らは、今後の研究として、非線形銀河バイアス、より複雑な固有整列モデル、およびバリオンフィードバックを含めること、ならびにRomanのミディアムティアとワイドティアの組み合わせを調査することを提案している。