Large-scale halo velocity correlations and the impact of finite simulation volumes

技術要約：大規模なハロー速度相関と有限シミュレーション体積の影響

問題提起
ダークマターおよびハローの余剰速度場から測定される速度相関関数は、理論的な線形予測と比較して、大規模スケールにおいて振幅の系統的な抑制を示す。この不一致は、速度統計を用いて成長率パラメータ $f(z)\sigma_8(z)$ を制約し、重力をテストする精密宇宙論にとって特に問題となる。この抑制は、シミュレーションボックスの有限な体積に起因すると考えられている。この有限な体積はフーリエモードにハードカットオフを課し、現実の宇宙には自然に存在する長波長の摂動を排除してしまう。この効果が大規模な相関を変化させることは既知であるが、速度統計への具体的な影響、ハロー質量への依存性、およびシミュレーションベースのモデリング（エミュレータなど）におけるパラメータ推論への結果については、十分に定量化されていない。

手法
著者らは、DARK QUEST プロジェクトのN体シミュレーション・スイートを用いてこの問題を調査しており、具体的には以下の2つのボックスサイズを比較している：

• L1G: 辺の長さ $1 \, h^{-1}$ Gpc
• L2G: 辺の長さ $2 \, h^{-1}$ Gpc

長波長モードの影響とボックスサイズの影響を切り離すために、著者らは制御されたサブボリューム解析を行っている：

1. 細分化: より大きな L2G ボックスを、$1 \, h^{-1}$ Gpc の辺の長さを持つ $L\times0.5_{2G}$ と $0.3 \, h^{-1}$ Gpc の $L\times0.15_{2G}$ の小さな体積に細分化する。
2. 比較: これらのサブボリュームで測定された速度相関関数を、独立した L1G シミュレーションと比較する。重要な点は、これらのサブボリュームは親となる L2G ボックスの長波長モードを保持しているのに対し、独立した L1G シミュレーションは保持していないという点である。
3. モデリング: 理論的な速度自己相関関数は、物質パワースペクトル $P_m(k)$ のハンケル変換を用いてモデル化される。著者らは、欠落している長波長モードを模倣するために、積分の下限を表す自由パラメータ $k_{min}$ を導入する。
4. 解析: モデルをシミュレーションデータに適合させ、$f\sigma_8$ を制約する。その際、（a）$k_{min}$ を自由パラメータとして周辺化するシナリオと、（b）$k_{min} \approx 0$ と固定する（有限体積効果を無視する）シナリオの2つをテストする。また、ハローをグループ質量範囲とクラスター質量範囲に分割することで、これらの効果のハロー質量への依存性を分析する。

主な結果

• 抑制の起源: 速度相関振幅の抑制は、単に解析領域の物理的な大きさによるものではなく、シミュレーション体積における長波長モードの不在によって引き起こされる。より大きな L2G ボックスから切り出されたサブボリューム（$L\times0.5_{2G}$）は、独立した L1G ボックスと同じ体積であるにもかかわらず、元の L2G ボックスと同じレベルの抑制を示す。これは、より大きなボックスのモードの「記憶」がサブボリュームに保持されていることを裏付けている。
• 質量依存性: 抑制は質量に依存する。より質量の大きいハロー（クラスター規模、$M_h > 10^{14} h^{-1} M_\odot$）は、より質量の小さいハロー（グループ規模、$M_h \sim 10^{13} h^{-1} M_\odot$）と比較して、系統的に強い速度相関振幅の抑制を示す。
• パラメータ推論 ($f\sigma_8$):
  • $k_{min}$ の周辺化: $k_{min}$ が自由パラメータとして扱われ、周辺化された場合、回収された $f\sigma_8$ は、ボックスサイズに関わらず、シミュレーションの真の入力値（$f\sigma_8 = 0.4852$）と一致する。
  • 効果の無視: もし $k_{min}$ をゼロに固定した場合（有限体積効果を無視した場合）、$f\sigma_8$ の制約は著しくバイアスを受ける。これは、分離ビン間の相関を無視した場合（すなわち、完全な共分散行列を使用しなかった場合）に特に顕著である。完全な共分散行列を用いた場合でも、最良適合値のシフトは残るが、それは $1\sigma$ 以内に収まっている。
• 共分散の重要性: $f\sigma_8$ のバイアスは、共分散行列の非対角要素を無視した際に増幅される。

意義と主張
本論文は、N体シミュレーションにおける有限体積効果が、速度統計のシミュレーションベースのモデリングにおいて、極めて重要でありながら、しばしば過小評価されている系統誤差であることを主張している。

• モデリングへの示唆: エミュレータのようなシミュレーションベースのアプローチにおいて、欠落している長波長モードを補正せずに有限のボックスの生出力をそのまま使用することは、バイアスのかかった宇宙論的制約をもたらす。著者らは、非バイアスな $f\sigma_8$ の制約を回収するためには、有効な長波長カットオフ（$k_{min}$）を自由パラメータとして導入することが必要なステップであると提案している。
• 観測との区別: 著者らは、シミュレーションのアーティファクトと観測のウィンドウ関数の根本的な違いを強調している。観測のウィンドウ関数は明示的にモデリングしデコンボリューションできるが、シミュレーションにおけるハードカットオフは初期条件に刻み込まれており、非線形進化を通じて持続するため、より強固な制限となる。
• 物理的解釈: 質量依存の抑制は、バイアスのあるトレーサーを分析する際に、有限体積効果と速度バイアスを明確に分離できないことを意味する。これは、シミュレーションで較正されたモデルを用いたアインシュタイン等価原理のテストにおいて、未考慮の有限体積効果が、真の物理的信号を模倣したり、あるいは覆い隠したりする可能性があることを示唆しており、複雑な問題となっている。

結論として、余剰速度および運動学的スニェイ・ゼルドビッチ（kSZ）データの将来の解析においては、シミュレーションを理論的テンプレートとして用いる際、ロバストな宇宙論的解釈を保証するために、欠落している長波長モードを整合的に考慮しなければならない。