Accurate Galaxy Cluster Shear and Mass Calibration for LSST with AnaCal

技術要約：AnaCalを用いたLSSTにおける正確な銀河団シアーおよび質量校正

問題提起
銀河団宇宙論は、観測量と銀河団質量の間の正確な関係を確立すること、すなわち質量校正に依存している。弱重力レンズ効果は、動力学的または熱的平衡の仮定に依存せずに全質量（ダークマターを含む）を探索できるため、この校正のための標準的なツールである。しかし、銀河団環境においては重大な課題が生じる。銀河団中心付近のシアー信号は非線形領域（例：$|g| \sim 0.15$）に達することがあり、多くのシアー推定法は明示的な検証を必要とする摂動に基づく仮定に依存している。さらに、銀河団レンズは、標準的な宇宙シアー測定には存在しない特有の複雑さを導入する。これには以下が含まれる：

• 高シアー領域： 中心部（$\sim 0.2$ Mpc）における接線方向シアーは、摂動ベースの校正テストで通常使用される値よりも数倍大きい。
• 不均質な偏向： 銀河団レンズ場は背景銀河の見かけの位置を変化させ、それによって重なりやブレンディングの構成が変化し、検出および形状測定に影響を与える。
• 半径依存性： コンバージェンス場は半径方向に変化し、それが半径方向に減少する平均シアー応答を引き起こす可能性がある。
• 高次効果： 銀河団中心付近では、高次の歪み（例：フレクション）や銀河団内ハローからの異方的な光が、偽の信号を導入する可能性がある。

LSST ダークエネルギー科学コラボレーション（DESC）は、fiducial（基準となる）シアー推定法としてAnaCalおよびmetadetectionを選択したが、それらの性能は、銀河団レンズに特徴的な高シアーおよび不均質な条件下での厳密なテストが行われていなかった。

手法
著者らは、AnaCalシアー推定器をテストするために、10年間のLSSTデータを代表する画像シミュレーションを利用している。その手法は以下の通りである：

1. シミュレーション・フレームワーク： descwl-shear-sims および xlens を用い、現実的な銀河団レンズ画像シミュレーションを生成した。彼らは、$z_l = 0.25$ のレンズ赤方偏移において、NFWハローとして9つのハロー質量（$10^{14} M_\odot$ から $10^{15} M_\odot$）をモデル化した。ソース銀河は $z_s = 1.0$ である。

   • シミュレーションには、質量ごとに2178回の実現数が含まれ、合計で約1700 deg$^2$ をカバーしている。
   • シアー推定の影響を孤立させるため、シミュレーションでは銀河団メンバー銀河、銀河団内光、および赤方偏移校正誤差を除外し、高シアーおよび不均質な偏向場の影響にのみ焦点を当てた。
   • 定数円形Moffat PSF（$FWHM = 0.8''$）を使用し、Rubin天文台の予想される画像品質を表現した。

2. AnaCalの実装： 本研究では、自動微分とシェップレット（shapelet）基底を用いて、フル測定パイプライン（観測量、検出、および選択）のシアー応答を解析的に導出するAnaCalフレームワークを採用している。

   • 観測量： 画像は検出カーネルおよび極座標シェップレット基底関数上に投影される。
   • 応答計算： 制御されたシアー摂動に対して観測量を微分することで、シアー応答を計算する。
   • ノイズ補正： 画像に合成ノイズの実現数を加えることで、ノイズバイアスを補正するための「リノイジング（renoising）」スキームを適用する。
   • 重み付け： 選択および検出の重みは、信号対雑音比（SNR）および銀河のサイズに基づいて適用され、90度回転に対して不変となるよう設計されている。

3. 解析およびフィッティング：

   • シアー推定： 接線方向およびクロス楕円率が計算され、応答行列は接線フレームへと回転される。極めて重要な点として、コンバージェンス場の半径依存性によるバイアスを避けるため、著者らは各角度ビンに対して個別にシアー応答を使用しており、すべてのビンに対して平均化は行っていない。
   • 質量フィッティング： フォワードモデリング手法を用い、10年間のLSST予測から導出された解析的共分散行列を組み込んだ尤度関数を使用して、観測された接線方向縮小シアープロファイル（$g_T$）をNFWモデルにフィットさせる。
   • バイアス定量化： シアーのマルチプリカティブ・バイアス（$m$）を $m = a + b(g_{true})^2$ の形式でフィッティングし、線形および第3次の非線形バイアスをテストする。

主な結果

• シアー回復： AnaCalは、軽微に高いシアーレベル（$|g| \sim 0.15$）においても、最小限のバイアスで入力シアーを回復する。
• 半径方向の応答傾向： 本研究は、コンバージェンス場の半径依存性に起因する、半径方向に減少する平均シアー応答を確認した。著者らは、もしこの効果がモデル化されていない場合（すなわち、グローバルな平均応答が使用された場合）、シアー推論にバイアスが生じることを示している。
• 第3次バイアス： 銀河団中心付近において、縮小シアーの第3次における正のシアー推定バイアスが検出された。
• 質量校正の正確性： 高シアー領域および第3次バイアスの存在にもかかわらず、最終的な質量推定への影響は無視できる。$[10^{14} M_\odot, 10^{15} M_\odot]$ の範囲の銀河団レンズ質量について（$z=0.25$ において $\sim 0.2$ Mpc のスケールカットを採用）、得られた平均銀河団質量バイアスは、理想的な設定下で $0.24 \pm 0.26\%$ であった。この小さなバイアスは、高シアー領域に存在する銀河の割合がわずかであり、それらの測定値が共分散行列によってダウンウェイトされていることに起因する。

意義および主張
本論文は、これらの結果が AnaCalがLSST時代における正確な銀河団質量校正のための堅牢なツールであることを実証している と主張している。具体的には：

• 摂動ベースのシアー推定（AnaCal経由）は、半径方向の応答が適切にモデル化されている限り、銀河団中心に典型的な高シアー領域においても正確であり続けることを検証した。
• 非線形性および高シアー効果が、宇宙論的に関連する銀河団サンプルに対する質量推定の精度を著しく低下させないことを確認した。
  本研究は、選択されたシアー推定法が、重大な系統誤差を導入することなく、銀河団宇宙論の厳格な要件を満たすことができることを保証する、LSST DESCにとって極めて重要な検証ステップを提供するものである。

著者らは、自らのシミュレーションには銀河団メンバー、銀河団内光、またはフレクションのような高次歪みが含まれておらず、高シアーおよび不均質な偏向場に対するシアー推定器の性能を孤立させていることを明記している。