Constraining Galaxy Cluster Triaxiality via Weak Lensing -- I. Preparation for the Rubin Data Beyond Leading Order

技術的サマリー：弱い重力レンズによる銀河団の三軸性の制約 — I. ルビン天文台データにおける先頭次を超える準備

問題提起
銀河団は宇宙における最大の重力崩壊構造であり、球対称ではなく三軸性の 3 次元質量分布を有している。この幾何学は、球対称密度モデルを仮定した場合、個々の銀河団の弱い重力レンズ（WL）質量測定において、約 20% の大きなばらつきをもたらす。さらに、三軸性と投影効果の組み合わせは、検出および選択バイアスを誘発する可能性がある。すなわち、視線方向に整列した銀河団はより豊かで検出しやすく、垂直方向のものはより丸く検出しにくく見える。銀河団集合体の投影形状（楕円率）を正確に制約することは、質量測定のばらつきを低減し、選択バイアスを補正し、宇宙論パラメータやダークマター物理学（例えば、自己相互作用ダークマターの制約など）を探求する上で不可欠である。過去の研究では銀河団の形状が測定されてきたが、高度な多重極モデルを用いたダークエネルギー調査（DES）第 3 年（Y3）データの包括的な集合分析は行われていなかった。

手法
著者らは、ダークエネルギー科学コラボレーション（DESC）によって開発された**クラスターレンズ質量モデリング（CLMM）**ソフトウェアを用いて、DES Y3 の公開データリリースを分析する。分析パイプラインは以下の通りである：

1. データ選択：

   • 銀河団： redMaPPer カタログから 9,962 個の銀河団（豊かさ $20 < \lambda < 222$、赤方偏移 $0.2 < z < 0.65$）を選択。中心確率が高い（$p_{cen} > 0.9$）もののみを用いて、緩和された系と正確な中心定位を確保する。
   • 源銀河： DES Y3 メタ較正カタログからの形状と、方向性近傍適合（DNF）カタログからの測光赤方偏移（photo-z）を使用する。
   • 検証： 楕円 Navarro–Frenk–White（NFW）ハローをシミュレートしたモックカタログと、DES Y1 redMaPPer 銀河団の一部を交差検証用として使用する。

2. 方位決定：

   • 各銀河団の主軸方位は、中心銀河のみではなく、衛星メンバー銀河の分布（重み付き 2 次モーメントを使用）によって決定され、中心定位誤差の影響を最小化する。
   • 弱い重力レンズのせん断推定値は、スタッキング前にこの主軸に揃うように回転される。

3. 理論的モデリング（先頭次を超える）：

   • 表面質量密度は、二次（四重極）まで多重極展開を用いて展開される。
   • 単極子（$\Sigma_0$）： 円対称成分を表す。著者らは、楕円率による質量再分配を考慮するため、単極子に二次の補正項（$\epsilon^2$）を含める。
   • 四重極（$\Sigma_2$）： 楕円変形を表す。
   • 分析では、質量、濃度、および投影楕円率（$\epsilon$）を制約するために、単極子と四重極の信号を同時にフィットする。

4. 系統的補正：

   • 希釈因子： 方位定義に用いられた衛星銀河分布と真のダークマターハロー方位の間の不一致を補正する。これは銀河団の豊かさに基づいたモンテカルロサンプリングによって推定される。
   • ブースト因子： 測光赤方偏移誤差による、レンズ効果を受けない銀河団メンバーの背景源銀河への混入を補正する。これは半径方向ビン内の銀河の赤方偏移分布を用いて推定される。
   • 半径カット： 300 kpc 以内の源は、ブレンディングと強い重力レンズ効果の影響を軽減するために除外される。

主要な貢献

• 初の DES Y3 集合形状測定： 本論文は、弱い重力レンズを用いて DES Y3 redMaPPer カタログからの大規模な銀河団集合体の有効投影楕円率を測定した最初の研究である。
• 先頭次を超える分析： 本分析では、単極子項に対して楕円率の二次まで多重極展開を拡張した。これにより質量プロファイルに対してパーセントレベルの補正が得られ、先頭次近似と比較して質量測定の精度が向上する。
• 堅牢な系統誤差の扱い： 本研究では、希釈（方位不一致）とブースト（混入）因子を明示的に補正するパイプラインを実装・検証し、将来の大規模サーベイに向けて CLMM ツールが準備できていることを示した。
• 検証： 本手法は、既知の投入形状を持つモックカタログおよび DES Y1 データに対して検証され、測定手法（個体別対 HEALPix 別応答）の変動に対する一貫性と堅牢性が示された。

結果

• 投影楕円率： ベースラインサンプルにおいて、測定された投影ハロー楕円率は $\epsilon = 0.310^{+0.017}_{-0.016}$ であり、軸比 $q = 0.527^{+0.018}_{-0.019}$ に対応する。これは非ゼロの楕円率の $\sim 19\sigma$ 検出を表す。
• 性質への依存性： 分析の結果、現在のサンプルサイズと統計的精度の範囲内では、投影楕円率が銀河団の質量（豊かさ）や赤方偏移に統計的に有意に依存することは見出されなかった。高/低豊かさおよび高/低赤方偏移ビン間の差は約 $0.4\sigma$ である。
• 質量バイアス： 三軸モデルフィットを球対称モデルと比較すると、同じ単極子信号に対して三軸モデルは約 5% 低い質量推定値を与える。これは楕円 NFW プロファイルに対する理論的期待と一致する。
• モック検証： 測定アルゴリズムは、真の軸比 $q > 0.35$ に対して、投入された軸比を 5% 未満のバイアスで回復する。また、ハロー濃度の変化（$c=2, 4, 6$ でテスト）に対して分析が感度を持たないことが確認された。

意義と主張
本論文は、この研究を**ヴェラ・C・ルビン天文台の「宇宙と時間の遺産サーベイ（LSST）」**に向けた重要な準備段階として位置づけている。著者らは以下を主張する：

1. 三軸性モジュールを備えた CLMM ツールは、大規模な銀河団の弱い重力レンズデータを効率的かつ一貫して分析する能力を有している。
2. 実証された手法は、LSST、ユークリッド、ローマン宇宙望遠鏡といった今後の広域サーベイに直接適用可能である。
3. より大規模なサンプルサイズと深層観測を有する将来のデータセットは、銀河団の楕円率に対する制約を大幅に強化する（DES Y3 と比較して LSST 第 10 年（Y10）では誤差範囲が 82% 削減されると推定される）ため、ハロー形状の質量および赤方偏移への依存性に関するより精密なテストが可能となり、銀河団宇宙論における系統誤差の制御が向上する。
4. 本論文はシリーズの第一弾であり、将来の研究では宇宙論シミュレーション（SkySim5000）を活用して、選択バイアス、中心定位誤差、およびハローの三軸性と周囲の大規模構造との関係をさらに定量化する予定である。