Assembly bias and the redshift evolution of intrinsic alignments for LRGs

技術要約：アセンブリー・バイアスとLRGにおける固有整列の赤方偏移進化

問題提起
銀河の固有整列（Intrinsic Alignment; IA）は、弱重力レンズ観測における主要な天体物理学的系統誤差であり、正しくモデル化されない場合、宇宙論パラメータの推論にバイアスをもたらす可能性がある。非線形整列（NLA）モデルなどの現在のモデリング手法は、多くの場合、整列振幅がハロー質量のみに依存すると仮定している。しかし、EuclidやLSSTのような高統計観測の到来に伴い、観測的な不確実性によってこれまで隠されていた二次的な依存性が重要になる可能性がある。具体的には、アセンブリー・バイアス（ハローの性質が質量以外の形成履歴にも依存すること）の影響と、整列信号の赤方偏移進化について、将来の宇宙シア解析におけるロバストなモデリングを確保するために、大規模な流体動力学シミュレーションを用いて調査する必要がある。

手法
著者らは、バリオン物理学と標準的なフィードバック変異を含むFLAMINGOシミュレーション・スイート（Schaye et al. 2023）の2.8 Gpcボックスサイズ・ランを利用している。本研究は、ハロー質量 $M_h > 10^{12} M_\odot$ および300個以上のスター粒子（$\sim 2.22 \times 10^{11} M_\odot$）を持つ銀河によって定義される、Luminous Red Galaxy（LRG）に類似した中心銀河のサンプルに焦点を当てている。

主な手法ステップは以下の通りである：

1. 形成赤方偏移の定義： 形成赤方偏移（$z_f$）は、ハローが現在の質量の半分（$0.5 M_{200\text{mean}}$）を含んでいた時の赤方偏移として定義される。この定義は、ハローの集中度（$c_{200\text{crit}}$）との強い相関によって検証されている。
2. 相関関数の推定： 位置－位置（$gg$）および位置－形状（$g+$）相関関数を、一般化されたLandy-Szalay推定法を用いて計算し、IATheoryコードを用いたNLAフレームワークでモデル化する。
3. アセンブリー・バイアス解析： 固定質量における $z_f$ に対する整列振幅（$A_1$）の依存性を分離するため、サンプルをハロー質量と形成赤方偏移によってビン分割する。
4. 赤方偏移進化のモデリング： さまざまなシミュレーション・スナップショット（$z=0$ から $z=1.5$）にわたって、整列振幅とハロー質量の間のべき乗則の関係（$A_1 \propto M^{\beta_M}$）をフィッティングする。
5. 銀河のトラッキング： 赤方偏移進化を単純な質量成長から区別するため、HBT-HERONSサブハローファインダーのTrack_idを用いて、中心銀河をスナップショット間で追跡する。サンプルは初期赤方偏移（$z_i$）で選択され、$z=0$ まで追跡される。合体した場合、またはスター粒子の閾値を下回った場合はオブジェクトを除外する。

主な結果

• アセンブリー・バイアスの依存性： 整列振幅はハロー質量だけでなく、形成赤方偏移にも依存する。固定されたハロー質量において、より早期に形成された（高い $z_f$ を持つ）ハローとその中心銀河は、系統的に高い整列振幅を示す。この効果は、$A_1$–質量関係における垂直方向のシフトとして現れるが、べき乗則の傾き自体はほとんど変化しない。
• パラメータの赤方偏移進化： 質量－整列のべき乗則の振幅（$\alpha_M$）は、赤方偏移とともに着実に増加する。傾き（$\beta_M$）は $z \sim 1$ まで赤方偏移とともに減少し、その後平坦化する。著者らは、$\alpha_M(z)$ および $\beta_M(z)$ に対する経験的な二次関数近似を提供している（式13および14）。
• 質量－赤方偏移モデル： 与えられた赤方偏移と平均ハロー質量に対してLRGサンプルの整列振幅を予測するための、結合質量－赤方偏移モデルが提案されている（式12–14）。
• トラッキング解析： 高い初期赤方偏移（$z_i$）から銀河を追跡すると、質量降着と一致して、赤方偏移が減少するにつれて整列振幅が増加する。しかし、質量－赤方偏移モデルの予測値で正規化した後でも、系統的なオフセットが残る。すなわち、高い $z_i$ から追跡された銀河は、低い $z_i$ から追跡されたものよりも高い整列振幅を保持している。これは、アセンブリー・バイアスが、与えられた質量と赤方偏移において整列振幅に不可避な散らばりを導入することを裏付けている。
• べき乗則からの逸脱： べき乗則のフィッティングは、低質量（$< 10^{13} M_\odot$）および高質量（$> 3 \times 10^{14} M_\odot$）において逸脱する。低質量域では、著者らはべき乗則からの減少を観測しており（これは一部の重力のみのシミュレーション結果とは対照的である）、高質量域では、進行中のハロー形成に起因する減少が観測される。

意義と主張
本論文は、アセンブリー・バイアスが固定質量における振幅の根本的な広がりをもたらすため、弱重力レンズ観測におけるIAモデリングにおいて、整列振幅に対して任意に小さな事前分布の範囲を設定することはできないと主張している。既存のモデル（NLAなど）の関数形式自体は、アセンブリー・バイアスを考慮するために修正を必要としないかもしれないが、この二次的な依存性を収容するために、振幅パラメータの事前分布の制約を広げる必要がある。

さらに、本結果は、重力のみのシミュレーションを用いたハロー占有分布（HOD）に基づく整列モデルの実装を困難にするものである。なぜなら、HODモデルは通常、質量のみに基づいて特性を割り当てるため、明示的に組み込まれない限り、形成時間の依存性を捉えることができないからである。本研究は、IAの赤方偏移進化モデルに関するシミュレーションに基づくガイドラインを提供しており、LRG様サンプルに対する質量－赤方偏移依存性の整列振幅に関する経験的な近似式を提供している。

限界
著者らは、FLAMINGOシミュレーションの解像度が銀河ディスクの進化を解像するには不十分であること（ディスク成分の割合は $z \sim 1.5$ まで 1% 未満である）を指摘している。したがって、結果は厳密には純粋なLRGサンプル（楕円銀河）にのみ適用され、楕円銀河からディスク銀河への進化や、ディスク銀河の整列信号については考慮していない。また、質量－赤方偏移モデルは、ハロー質量が $10^{13}$ から $10^{14.5} M_\odot$ の範囲であり、赤方偏移が $z \sim 1.5$ までの場合に有効である。