Modeling Redshift Uncertainties in Roman Weak Lensing Cosmology

技術要約：Romanの弱重力レンズ宇宙論における赤方偏移不確実性のモデリング

問題提起
Nancy Grace Roman宇宙望遠鏡を用いた弱重力レンズ（宇宙剪断）による宇宙論的制約は、系統誤差、特にソース銀河の赤方偏移分布（$n(z)$）の不確実性に非常に敏感である。銀河サンプルをトモグラフィー・ビンに分割することで宇宙論パラメータの推定精度は向上するが、不完全なフォトメトリック赤方偏移（photo-$z$）の決定は、物質密度（$\Omega_m$）や質量揺らぎの振幅（$S_8$）といったパラメータにバイアスをもたらす。標準的な緩和戦略は、各ビンの赤方偏移分布の平均を調整する「平均シフト」パラメータ（$\Delta z_i$）に依存することが多いが、これでは複雑な形状の不確実性を捉えることはできない。DES Y3などの先行研究で示されているように、全形状の不確実性は宇宙論的推論に影響を与える可能性があり、Romanのような次世代のStage IVサーベイには、より堅牢なモデリング手法が必要となる。

手法
本研究では、Roman High Latitude Imaging Survey (HLIS) Cosmology Project Infrastructure Team (PIT) のパイプラインであるCobaya-Cosmolike Joint Architecture (CoCoA) 内において、赤方偏移分布の全形状の不確実性をモデル化するための、最適化された主成分分析（PCA）フレームワークを実装し、検証を行う。

1. シミュレーション生成： 著者らは、クラスター銀河集団の欠損を修正し、カラー割り当てスキームを洗練させたBuzzardシミュレーションを改良したCardinalシミュレーション・スイートを利用する。彼らは、様々なRoman観測シナリオにわたって100万個の赤方偏移分布の実現体（$N_{sim} = 10^6$）を生成する。これらのシナリオは、面積、深さ、フィルターセットを変化させた3つのワイド・ティア構成（W1, W2, W3）と4つのディープ・ティア構成（D1–D4）を組み合わせ、サーベイ設計がphoto-$z$校正にどのように影響するかを評価する。
2. Photo-$z$推論： 赤方偏移分布を推定するために、Self-Organizing Map Photo-$z$ (SOMPZ) フレームードワークを採用する。この手法は、測光に基づいて銀河を「ワイド・ティア」および「ディープSOM」セルにグループ化し、完全な分光赤方偏移を持つディープフィールド・サンプルを用いてワイドフィールドの分布を校正する。サンプル分散とショットノイズは、3ステップのディリクレ・サンプリング法を用いて伝播される。
3. PCAの実装：
   • 基底の構築： 著者らは、一連の実現体から平均赤方偏移分布を差し引くことで、差分行列 $\Delta$ を構築する。
   • 最適化： [1] の処方に基づき、赤方偏移分布に関する宇宙剪断相関関数（$\xi_{\pm}$）のヤコビアンから導出されたフィッシャー行列（$F$）を用いて、PCA基底を最適化する。この重み付けにより、主成分（PC）が単なる $n(z)$ の統計的分散ではなく、剪断信号に最も関連する変動を捉えるようになる。
   • デコーディング： デコーディング行列 $D$ を用いて、赤方偏移分布を $\hat{n} = \bar{n} + \sum u_i U_i$ として再構成する（ここで $u_i$ はPC振幅、$U_i$ は基底モードである）。
4. 解析： 宇宙論パラメータ（$\Omega_m, S_8$ など）を推論するために、Cobayaを用いたマルコフ連鎖モンテカルロ（MCMC）解析を行う。本研究では、2つの緩和スキームを比較する：
   • 平均シフト： 9つの剰余パラメータ $\Delta z_i$（各トモグラフィー・ビンにつき1つ）を変化させる。
   • PCA： 平均シフトを固定したまま、主要な重み付きPCの振幅 $u_i$ を変化させる。

主な結果

• 整合性チェック： CoCoAにおける重み付きPCAの数値的実装は、[1] で提案された線形近似近似と一致している。重み付きPCAを用いたフル $\chi^2$ 評価は、PCの数が増えるにつれて線形近似と同様に収束し、パイプラインの妥当性が示された。
• バイアス緩和性能：
  • 低/中間程度の誤校正： シミュレーションされたデータベクトルが基準分布に対して強く誤校正されていない場合（初期 $\chi^2_0$ が低い場合）、平均シフト法とPCA法の両方が一貫した宇宙論的制約を与える。
  • 強い誤校正： 顕著な誤校正があるシナリオ（初期 $\chi^2_0$ が高い場合）では、平均シフト法は $\Delta z_i$ への事前分布を大幅に広げない限り（$\sigma=0.003$ から $\sigma=0.01$ へ）、バイアスのない制約を回復することが困難である。これは制約力を低下させる。対照的に、PCA法は少数のPC（例：5つのPC）を含めることで $\Omega_m$ および $S_8$ のバイアスを効果的に緩和し、多くの場合、9つのパラメータを持つ平均シフトモデルよりも少ないパラメータ数で、同等またはよりタイトな制約を達成する。
• シナリオ感受性： 本研究は、PCAの有効性が訓練データに依存することを特定した。「混合」シナリオ（PCがテストデータとは異なるサーベイ構成で訓練された場合）では、性能が変動する。具体的には、W2-D3 シナリオ（40 deg$^2$ の面積を持つディープ・ティア）は困難であることが判明した。W2-D3から導出されたPCAは、W2-D1テストケースに対するバイアス緩和に失敗した。これは、ディープフィールドの面積が大きくなることで宇宙分散が減少し、それによって、基準シナリオの特定の揺らぎに関する情報の有用性が制限されたことに起因する。
• PCの飽和： 解析によれば、20〜30個以上のPCを追加しても、情報量は飽和するため、収穫逓減が見られる。

意義および主張
本論文は、最適化されたPCAアプローチが、Romanの弱重力レンズ宇宙論における赤方偏移分布の不確実性をモデル化するための、柔軟かつ有望なフレームワークを提供すると主張している。この手法は、穏やかなシナリオにおいては標準的な平均シフト・パラメトリゼーションと同等の性能を示すが、強い誤校正条件下において、制約力を低下させるような事前分布の緩和を必要とせずに、バイアスを緩和する優れた能力を示す。著者らは、これらの結果がCoCoAパイプラインを検証するために設計された、ストレス・テスト済みの構成から得られたものであることを強調している。結論として、PCAベースのフレームワークは、photo-$z$ の不確実性をモデル化するために実行可能であり、LSSTやEuclidなどの他のイメージング・サーベイにも転用できる可能性があるとしている。ただし、今後の課題として、動的なダークエネルギーモデルや $3\times2\text{pt}$ 解析の検討が挙げられている。