KiDS-Legacy: Joint analysis of second- and third-order cosmic shear

技術要約：KiDS-Legacyにおける2次および3次宇宙剪断の結合解析

問題と背景
大規模構造による弱重力レンズ効果は、主要な宇宙論パラメータのプローブである。しかし、2次統計量（2点相関関数）に依存する標準的な解析は、主にパラメータの組み合わせである $S_8 = \sigma_8(\Omega_m/0.3)^{0.5}$ に対して感度を持つ。この感度は、物質密度パラメータ（$\Omega_m$）とクラスターリング振幅（$\sigma_8$）の間の縮退を生み出し、これらを独立して制約することを制限している。さらに、観測された遅い時刻の弱重力レンズ測定（しばしばより低い $S_8$ 値を示す）と、Planckによる初期宇宙の宇宙マイクロ波背景放射（CMB）観測との間に存在する、$S_8$ テンションとして知られる持続的な不一致が存在する。高次統計量（HOS）は、密度場の非ガウス的特徴を捉えることでこれらの縮退を打破できるが、その適用は計算の複雑さとモデリングの課題によって限定されてきた。本論文は、第3次統計量を通じて非ガウス情報を組み込むことが、宇宙論的制約を改善し、$S_8$ テンションを解決するかどうかを検証する必要性に取り組むものである（Kilo-Degree Surveyの最終データリリースであるKiDS-Legacyにおいて）。

手法
著者らは、967.4 deg²をカバーし、有効銀河数密度が $n_{eff} = 8.79$ arcmin$^{-2}$ であるKiDS-Legacyカタログを用いて、2次および3次の宇宙剪断の結合解析を行う。

• 統計量:
  • 2次: 2′から300′のスケールにおける、Eモード/Bモード積分に関する完全直交集合（COSEBIs）を使用する。
  • 3次: 4′から32′のスケールにおける、3次アパチャスマスモーメント（$\langle M_{ap}^3 \rangle$）を利用する。解析には、等スケールおよび不等スケールの構成の両方が含まれる。
• モデリングの進展: 前回のKiDS解析（例：Burger et al. 2024）と比較して、本研究ではいくつかの手法的な改善を実施している：
  • 固有配置（Intrinsic Alignments, IA）: 赤方偏移と恒星質量への明示的な依存性を持つIAモデル（NLA-Mモデル）を採用し、スケールに依存しない、あるいは純粋に赤方偏移に依存するモデルを超えたものとしている。
  • バリオンフィードバック: 新しいバリオン補正モデル（Burger et al. 2026）を採用し、流体計算によるFLAMINGOシミュレーションに対して検証を行っている。このモデルは複数のフィードバックシナリオ間を補間し、高温ガスの質量（$M_c$）と放出されたガスの密度（$\eta$）という2つの主要なパラメータに焦点を当てている。
  • 高次補正: シミュレーションされた剪断および収束マップを用いて、縮退剪断（reduced shear）およびソース・クラスタリングの補正を適用する。
  • エミュレーション: 推論を加速させるため、ニューラルネットワーク・エミュレータを開発した。具体的には、ソースの赤方偏移分布をエミュレータから切り離し、任意のサーベイへの柔軟な適用を可能にする、赤方偏移分解された $\langle M_{ap}^3 \rangle_z$ 用のエミュレータを訓練した。特定のIAモデルに合わせたCOSEBIモード（$E_n$）用の別個のエミュレータも訓練されている。
• 推論:
  • 1,944個のTakahashiシミュレーションによって推定された共分散行列のノイズを考慮した、パーシバルら（2014）の尤度定式化を用いたベイズ・フレームワークを使用する。
  • 効率的な事後分布サンプリングのために、ハミルトニアン・モンテカルロ法（NUTS）を採用する。
  • 第3次データの解析中の確認バイアスを防ぐため、厳格なブラインディング手順が実施された。
  • データベクトルの次元数は、シミュレーションに基づく共分散推定の限界と宇宙論的情報の獲得とのバランスをとるために最適化された（291エントリーに削減）。

主な結果

• 宇宙論的制約: 結合解析により、$\Omega_m = 0.297^{+0.056}_{-0.040}$ および $S_8 = 0.806^{+0.025}_{-0.023}$ が得られた。
• 制約の改善: 第3次統計量を第2次COSEBIsと組み合わせることで、$\Omega_m$ の制約が大幅に強化され（周辺誤差を約36%減少）、2次のみの解析と比較して $S_8$ の制らも約24%改善された。
• 性能指数（Figure of Merit, FoM）: $\Omega_m$–$S_8$ 平面におけるFoMは、2次のみの429から、結合解析による895へと2倍以上に増加した。これは、第3次統計量が $\Omega_ m$–$\sigma_8$ の縮退を効果的に打破していることを示している。
• 一貫性:
  • 第3次測定値は、様々なデータ分割（南北、スケールカット、トモグラフィービン）にわたる厳格な内部一貫性テスト（Tier 2およびTier 3）を通過しており、有意な偏差は見られない。
  • 結合制約は、KiDS-Legacyの2次解析の結果と完全に互換性があり、Planck CMBの測定値とも $1\sigma$ 以内で一致している。
  • 本結果は、KiDS-1000、DES、HSCによる他の2+3点解析の結果とも一致しているが、KiDS-Legacyの結果は、Planckにより近い、わずかに高い $S_8$ 値を支持している。
• モデリングの感度: 固有配置またはバリオンフィードバックを無視すると、$\Omega_m$ と $S_8$ の両方が低く偏った値になることを解析は確認している。しかし、特定のIAモデル（NLA, NLA-z, NLA-M）の選択は、最終的な宇宙論的制約を $1\sigma$ 以内で有意に偏らせることはない。

意義と主張
本論文は、KiDS-Legacyの結合解析が、採用されたモデリングの枠組み内で、他の弱重力レンズの結果やPlanckと一致する競争力のある宇宙論的制約を提供し、$S_8$ テンションの証拠は認められないことを主張している。本研究は、2次解析を制限するパラメータの縮退を打破するために、非ガウス的情報を抽出できる堅牢なプローブとして、第3次宇宙剪断が成熟した手法であることを実証している。

著者らは、不等アパチャースケールの導入が、特にバリオンフィードバックパラメータおよび $\Omega_m$ に対して、第3次統計量の制約力を大幅に高めることを強調している。さらに、赤方偏移分解されたエミュレータおよび柔軟なバリオン補正モデルの開発は、EuclidやLSSTといった次世代のStage-IVサーベイへの適用に向けたパイプラインを確立するものである。本研究は、縮退剪断やソース・クラスタリングといった系統誤差の補正を含む、第2次および第3次統計量の堅牢な結合モデリングが、現在および将来の弱重力レンズサーベイの宇宙論的可能性を最大化するために、実現可能であり不可欠であることを結論づけている。