Disentangling the Halo: Joint Model for Measurements of the Kinetic Sunyaev-Zeldovich Effect and Galaxy-Galaxy Lensing

技術要約：kSZとGGLの結合モデリングによるハローの分離

問題提起
現代の弱重力レンズ（WL）およびクラスタリング解析は、非線形スケールにおいて物質パワースペクトルを抑制するバリオンフィードバックをモデリングするという大きな課題に直面している。現在の流体解析シミュレーションは、この抑制の振幅とスケールに関する予測において幅広い差異を示しており、これが宇宙論的推論における重大な不確eties（不確実性）につながっている。運動学的スニヤエフ・ゼルドビッチ（kSZ）効果は、銀河ハローの外縁部における電子密度（したがってバリオン）の直接的なプローブとなる一方、銀河・銀河重力レンズ（GGL）は全物質分布をプローブするが、これらのプローブを個別にフィッティングすることはモデルの縮退という問題に直面する。具体的には、kSZのみの解析はしばло（しばしば）ハロー質量を推定するために不確実な恒星・ハロー質量関係（SHMR）に依存し、一方でGGLのみの解析は、外部からの制約なしにはバリオン密度プロファイルとダークマター分布を分離することに苦慮する。

手法
著者らは、BOSS（Baryon Oscillation Spectroscopic Survey）におけるCMASS銀河のkSZ効果測定とGGL測定の初の結合解析を提示している。本研究では以下を用いる：

• データ: 再構成された特異速度によって重み付けされた、Atacama Cosmology Telescope (ACT) DR5およびPlanckマップからのスタックされたkSZ温度プロファイル。GGLの過剰表面質量密度（$\Delta\Sigma$）プロファイルは、DES Y3（Dark Energy Survey Year 3）およびKiDS-1000の背景ソース銀河との相互相関から導出される。
• モデリング・フレームワーク: 軸対称な密度プロファイルの仮定に基づいた、高速で解析的なハローモデル。全物質密度は、バリオン成分（$\rho_b$）とダークマター成分（$\rho_{dm}$）に分けられ、それぞれが1ハロー項（中心ハロー）と2ハロー項（周囲の近傍）で構成される。
  • ダークマター: Navarro-Frank-White (NFW) プロファイルとしてモデル化される。
  • バリオン: 振幅、コアスケール分率、および内側・中間・外側のべき指数（power-law slopes）によってパラメータ化されたGeneralized NFW (GNNF) プロファイルとしてモデル化される。モデルは、バリオンが完全に電離したガスであることを想定し、バリオン質量分率が大きな切断半径（$r_b \approx 10 \times r_{200}$）において宇宙バリオン分率に達するようにプロファイルを正規化する。
  • 観測量: モデルは、ビームとの畳み込みを伴うガス密度の視線方向への積分を介してkSZ温度ゆらぎを予測し、投影された表面質量密度を介してGGL $\Delta\Sigma$ を予測する。
• 統計的手法: アフィン不変アンサンブル・マルコフ連鎖モンテカルロ（MCMC）アルゴリズムを用いた結合尤度解析が行われる。解析では、ハロー質量（$M_{halo}$）、バリオンプロファイルパラメータ（特に外側の傾き $\beta$）、および2ハロー振幅の摂動パラメータ（$A_{2h}$）についてフィッティングを行い、他のパラメータ（例：内側の傾き $\gamma$、中間傾き $\alpha$）は縮退を緩和するために固定する。

主な貢献

1. 結合モデリング・フレームワーク: 本論文は、単一のハローモデルの枠組み内でkSZとGGLの両方の信号を一貫してモデリングし、ハロー質量とバリオン密度パラメータの同時推論を可能にする公開コード glasz を導入している。
2. 成分の分離: 本研究は、kSZとGGLのデータを組み合わせることが、いずれか一方のデータセットのみをフィッティングする場合に存在する縮退を効果的に打破することを実証している。kSZ信号はハロー外縁部のバリオン分布を制約し、GGLは全ハロー質量と2ハロー項の振幅を厳密に制約する。
3. 直接的なシミュレーション比較: バリオン密度プロファイルを良好に制約することで、著者らは中間的なSHMRの推定に頼ることなく、流体解析シミュレーションと直接比較することを可能にしている。

結果

• パラメータ制約: 結合フィッティングにより、ハロー質量（$\log_{10}(M_{halo}/M_\odot) \approx 13.44 \pm 0.05$）およびバリオン分布の外側のべき指数（$\beta \approx 2.45^{+0.30}_{-0.23}$）に対するタイトな制約が得られた。バリオン密度プロファイルは、$0.3$から$50 \ h^{-1}$ Mpc のスケールにわたって制約されている。
• ハロー質量推定への影響: バリオンを無視したGGLのみのモデルは、結合kSZ+GGLモデルと比較して、推定されるハロー質量を約20%シフトさせることを著者らは発見した。これは、質量推定におけるバイアスを避けるためにバリオン物理学を含める必要性を強調している。
• バリオンフィードバックの強さ: 推定された外側の傾き $\beta$ は、Battaglia (2016) の流体解析シミュレーションにおけるAGNフィードバックモデルの予測よりも緩やかである。より緩やかな傾きは、バリオンがより広がっていることを意味し、観測されたフィードバックがこれらの特定のシミュレーションが予測するものよりも強いことを示唆している。
• シミュレーションとの比較: $\Delta\Sigma$ の小スケールにおける抑制は、Illustrisシミュレーション（より極端なフィードバックを特徴とする）には近いが、IllustrisTNG-300シミュレーションよりもIllustrisシミュレーションに近い結果となっている。

意義と主張
本論文は、この結合解析が銀河ハローにおけるバリオンとダークマターを「分離」するための堅牢な手法を提供すると主張している。モデルの縮退を緩和することにより、このアプローチは以下を可能にする：

• バリオン含有量を考慮した、より正確なハロー質量の推定。
• バリオンフィードバックの強さに対する直接的な経験的制約。これは、標準的な流体解析シミュレーションのいくつかが予測するものよりも強いフィードバックを明らかにしている。
• 現在、バリオンモデリングの不確実性のためにしばしば破棄されている小スケールの弱重力レンズ測定を、将来の宇宙論的解析で利用するための経路。

著者らは、現在のモデルがサテライト銀河やミスセンタリング効果を無視しており、これらが（質量と $\beta$ に対して $\sim 10-20\%$ レベルの）系統的なバイアスをもたらす可能性があると述べている。彼らは、精密宇宙論には将来的にハロ占有分布（HOD）フレームワークを組み込むことが必要になると考えているが、現在の解析的なフレームワークは、kSZとGGLのプローブの相補性を成功裏に実証している。