Intrinsic galaxy alignments in CAMELS simulations and the significant impact of baryon model

技術要約：CAMELSシミュレーションにおける銀河の固有整列とバリオンモデルの重大な影響

問題提起
銀河の固有整列（Intrinsic Alignments; IA）は、銀河形成のプローブであると同時に、弱重力レンズ解析における主要な系統誤差の要因でもある。フォトメトリック・サーベイによる現在の宇宙論的推論は、バリオンフィードバック過程（超新星および活動銀河核：AGN）や宇宙論パラメータが整列信号にどのように影響するかに関する、IAモデリングの不確実性によってボトルネックとなっている。理論的な枠組み（例：初期型銀河に対する潮汐整列や、後期型銀河に対する潮汐トルク）は存在するが、広範な天体物理学的および宇宙論的パラメータに対するIAの具体的な依存性は、特に異なる銀河集団（星形成銀河 vs 擬似静止銀河）を区別する場合において、依然として十分に制約されていない。

手法
著者らは、IllustrisTNGコードを用いた、1000個の流体解析シミュレーションからなるCAMELS（Cosmology and Astrophysics with MachinE Learning Simulations）スイート、特にラテン超格子の（LH）セットを利用している。これらのシミュレーションでは、6つの主要なパラメータ（宇宙論パラメータ $\Omega_m, \sigma_8$ および天体物理学的フィードバックパラメータ $A_{SN1}, A_{SN2}$ [超新星]; $A_{AGN1}, A_{AGN2}$ [AGN]）を変化させている。

• データ構築: 銀河の形状は、$z=0$ においてSUBFINDアルゴリズムで特定されたサブハロから測定されている。選択基準は、恒星質量（$M_{2r} > 10^8 M_\odot$）および粒子数（$N_s > 10$）である。楕円率は、3つのデカルト座標投影における恒星密度の第2モーメント・テンソルから算出され、計3,000個のカタログが得られた。
• 統計的尺度: 銀河の楕円率と背景のダークマター密度場との間の射影相互相関関数 $w_{m+}(r_\perp)$ を測定した。無効性テストとして、交差成分相関 $w_{m\times}$ を用いた。
• サブサンプル: 銀河は、特定の比星形成率（sSFR）の閾値 $10^{-10.5} \text{ yr}^{-1}$ に基づき、「星形成銀河」と「擬似静止銀河」の集団に分割された。
• 方位のみの解析: 方位（orientation）の影響を楕円率の大きさから分離するため、著者らはすべての楕円率の大きさを平均値 $\langle|\varepsilon|\rangle$ に正規化した「方位のみ」のカタログを構築し、対応する相関関数 $\tilde{w}_{m+}$ を算出した。
• パラメータ依存性: 整列振幅は $w_{m+}$ の加重平均として定義される。振幅とシミュレーションパラメータとの線形依存性をテストするために、ピアソンの相関係数とp値が計算された。また、線形整列モデル（NLA/TATT）との整合性をテストするために、比 $w_{m+}/w_{mm}$ も分析された。

主な結果

1. 検出と集団の違い: CAMELSシミュレーションにおいて、有意な固有整列信号が検出された。擬似静止銀河は、星形成銀河よりも約1桁大きい整列振幅を示す。
2. パラメータ依存性:
   • 宇宙論: 整列振幅 $w_{m+}$ は $\Omega_m$ および $\sigma_8$ と統計的に有意な相関を示す。しかし、信号を物質密度自己相関（$w_{mm}$）で正規化すると、$\sigma_8$ への依存性は消失し、これは線形整列モデルの予測と一致する。
   • フィードバック: 超新星フィードバックパラメータ（$A_{SN1}$ および $A_{SN2}$）と整列振幅との間に有意な相関が見られた。AGNフィードバックパラメータとの有意な依存性は観察されなかったが、これはシミュレーションの体積（$(25 h^{-1} \text{ Mpc})^3$）が限定的であり、大規模なモードが制限されていることに起因すると著者らは考えている。
3. 星形成 vs 擬似静止のメカニズム:
   • 整列振幅と $\sigma_8$ の相関は、星形成銀河では有意ではないが、擬似静止銀河および全サンプルにおいては強い。著者らは、全サンプルにおける依存性は、高い $\sigma_8$ 値において高度に整列した擬似静河の割合が増加することによって駆動されていると示唆している。
   • 超新星フィードバック（$A_{SN1}$）は、擬似静止銀河よりも星形成銀河において整列との強い相関を示す。逆に、擬似静止銀河においては、$A_{SN1}$ は整列ではなく特定の比星形成率を制御している。
4. 楕円率 vs 方位: 楕円率の大きさを正規化する（方位のみの解析）ことで、超新星フィードバックパラメータへの感度が保持される。これは、超新星フィードバックが銀河の楕円率の大きさだけでなく、ダークマター分布に対する銀河形状の方位にも影響を与えることを示している。
5. 比率分析: 比 $w_{m+}/w_{mm}$ は $\sigma_8$ への依存性を除去するが、$\Omega_m$ および超新星フィードバックパラメータとの有意な相関は保持する。著者らは、この比において単純なNLA/TATTの予測に反する驚くべき負の相関が $\Omega_m$ と見られることを指摘しており、さらなる調査が必要であるとしている。

意義と主張
本論文は、大規模な流体解析シミュレーション・スイートを用いて、宇宙論的および天体物理学的パラメータに対する銀河の固有整列の依存性を直接測定した最初の成果であると主張している。その結果は以下の点において重要な知見を提供する：

• モデリング: バリオンフィードバック（特に超新星）が整列信号に大きな影響を与えることを示しており、宇宙論パラメータ推定におけるバイアスを避けるために、これらのプロセスを考慮したモデルが必要であることを示唆している。
• 銀河進化: 銀河のタイプによって異なる整列メカニズムが働いていることを強調しており、超新星フィードバックが星形成銀河では方位の変化を駆動する一方で、擬似静河では星形成率を制御していることを明らかにしている。
• 将来のサーベイ: 次世代の弱重力レンズ・サーベイに向けたパイプライン検証およびシミュレーションに基づく推論に必要な制約を提供しており、仮定されたバリオンフィードバックモデルとIA処方との間の不一致が、重大なバイアスを導入する可能性があることを強調している。

著者らは、AGNフィードバックの結果については慎重な姿勢を保っており、検出された影響の欠如は小さなシミュレーション体積によるものである可能性が高いと明記しており、正規化された比における予期せぬ $\Omega_m$ 依存性の背後にある具体的なメカニズムを解明したとは主張していない。