CHEX-MATE: Cluster Multi-Probes in Three Dimensions (CLUMP-3D) II. Combined Gas and Dark Matter Analysis from X-ray, SZE, and WL

技術要約：CHEX-MATE: CLUMP-3D II. X線、SZE、およびWLを用いたガスとダークマターの結合解析

問題提起
階層的構造形成の標準モデルの下では、銀河団は球対称ではなく、本質的に三軸性を持つ。これらの構造に対して球対称モデルを適用することは、導出される特性、特に弱重力レンズ（WL）による質量推定において、クラスターの向きや視線方向（LOS）への伸長に対して非常に敏感であるため、重大なバイアスを導入する。三軸モデルはこれらのバイアスを軽減できる可能性があるが、三次元的な幾何学構造を制約するためには、深いマルチプローブ・データと物理的に動機付けられたモデルが必要となる。先行研究（Paper I, Kim et al. 2024, 以下K24）では、ガスの幾何学構造を制約するためにX線とSunyaev-Zel'dovich（SZ）データのみを用いた三軸フィッティングの定式化を確立した。しかし、このガスのみのアプローチでは、ガスと全質量分布の関係に関する仮定なしには、クラスターの全質量や、銀河団内ガス（ICM）における非熱的圧力支持を直接測定することはできなかった。

手法
本論文では、既存のX線およびSZ解析に弱重力レンズ（WL）データを統合することにより、K24の三軸フィッティング・パイプラインの拡張を提示する。手法は以下の通りである：

1. 三軸幾何学： モデルは、ICMガスの分布および重力ポテンシャルに対して三軸楕円体であることを仮定する。固有の幾何学構造は、観測者に対する軸比（$q_1, q_2$）およびオイラー角（$\theta, \phi, \psi$）によって定義される。
2. ICMモデリング： 電子密度は修正された$\beta$プロファイル（式3）に従い、電子圧力は一般化Navarro-Frenk-White（gNFW）プロファイル（式4）に従う。これらのプロファイルは視線方向に投影され、X線表面輝度（SB）、SZコンプトン$y$パラメータ、および分光温度のモデルマップを生成する。
3. 質量再構成： 全密度プロファイルに対して楕円対称性を仮定していた従来のアプローチとは異なり、本研究では、重力ポテンシャルがICMの三軸幾何学に従う（$q_{pot} = q_{ICM}$）と仮定する。これは、ポテンシャルが全密度よりも滑らかであるというシミュレーションの結果に基づいている。全質量プロファイルは、球対称基底におけるNFW形式によって制約され、その後、ポアソン方程式を用いて三軸格子上に数値的に反転され、楕円状のポテンシャルを生成する。このポテンシャルは投影され、WL観測量（縮退せん断マップ）を生成する。
4. フィッティング・フレームワーク： 結合データセットを適合させるために、emceeを用いたマルコフ連鎖モンテカルロ（MCMC）法が採用されている。尤度は、SZマップ、X線SB（クランピング・バイアスを軽減するための二重の2D/1Dアプローチを使用）、X線温度プロファイル、および2D WLせん断マップを組み合わせた総$\chi^2$統計量によって定義される。
5. 非熱的圧力： 全質量密度およびガス密度プロファイルが制約された後、一般化静水圧平衡の状態を仮定することで、非熱的圧力分率（$P_{nt}/P_{tot}$）が計算される。重力に対してクラスターを支えるために必要な全圧力は、重力ポテンシャルから導出され、非熱的成分は、ガス温度から導出される熱圧力とこの全圧力との差である。

主な貢献

• パイプラインの統合： 本論文は、WLデータをK24のガスのみの三軸パイプラインに統合することに成功し、球対称性の仮定に依存することなく、クラスターの幾何学、全質量、および集中度を同時に測定することを可能にした。
• ポテンシャルに基づくモデリング： 著者らは、全密度プロファイルではなく重力ポテンシャルに三軸性を適用するという新しいアプローチを導入しており、これにより、ICM等密度面と一致するより滑らかなポтенシャルを得ることができる。
• 検証： 本手法は、既知の入力パラメータ（$M_{200c}, c_{200c}$、および幾何学）を持つ模擬観測を用いて検証されている。テストの結果、パイプラインは入力値が事前分布から大きく外れる場合のわずかな集中度のバイアスを除き、入力質量と幾何学を正確に回収できることが示された。

結果：Abell 1689への適用
本パイプラインは、XMM-NewtonのX線データ、ACT/PlanckのSZデータ、およびSubaruのWLデータを用いて、CHEX-MATEクラスターPSZ2 G313.33+61.13（Abell 1689）に適用された。

• 幾何学： クラスターは天球面に対して視線方向に伸長しており、伸長係数は $R_{LP} = 1.27 \pm 0.02$ である。
• 質量： この伸長により、三軸フィットによる全質量は $M_{200c} = (13.69^{+1.56}_{-1.41}) \times 10^{14} M_\odot$ となる。これは、同じ手法を用いた球対称フィットによる $(17.77^{+2.00}_{-1.75}) \times 10^{14} M_\odot$ という値よりも有意に低くなっており、球対称の仮定によって導入されるバイアスを浮き彫りにしている。
• 集中度： 三軸フィットによる集中度は $c_{200c} = 8.55^{+2.20}_{-1.61}$ である。これは球対称の値である $9.99^{+2.26}_{-1.78}$ と一致しており、Abell 1689で見られる高い集中度が三軸性や方位の効果によるアーティファクトではないことを示唆している。
• 非熱的圧力： 非熱的圧力分率は、0.18–1.37 Mpcの半径範囲で測定された。この分率は中間半径で約20%の最小値を示し、最小半径および最大半径の両方で30%近くまで増加しており、典型的な測定精度は $\pm 5\%$ である。

意義
本論文は、個々の銀河団に対して偏りのない質量および集中度の測定を得るためには、X線、SZ、およびWLデータを三軸フレームワーク内で組み合わせることが不可欠であることを示している。三次元的な幾何学を明示的にモデリングすることで、球形解析を悩ませる質量と伸長の間の縮退を解消している。さらに、三軸基底で全質量プロファイルを再構成する能力により、クラスターの外縁部における非熱的圧力支持のロバストな測定が可能となり、降着物質の熱化効率について流体力学シミュレーションと比較可能な観測的制約を提供できる。Abell 1689に関する結果は、その高い集中度が投影効果ではなく固有の性質であることを示しており、また、球対称を仮定した場合よりも質量が低いことを示している。