Improved recipes for peculiar velocity power spectra using Evolution Mapping

技術的サマリー：進化マッピングを用いた特異速度パワースペクトルの改良されたレシピ

問題提起
現在および将来の銀河サーベイ（例：DES、DESI、Euclid、LSST）は、赤方偏移空間歪み（RSD）から宇宙論情報を抽出するために、非線形構造形成の精密な理論モデルを必要とする。RSD モデリングの重要な構成要素は、速度発散場（$\theta = -\nabla \cdot \mathbf{v}$）であり、特にその自己パワースペクトル $P_{\theta\theta}(k)$ と密度場との交差パワースペクトル $P_{\delta\theta}(k)$ である。

既存のフィッティング処方箋、特に Bel ら (2019、以下 B19) によるものは、以下の 3 つの主要な制限に苦しんでいる：

1. 数値収束性： 以前の較正は、離散粒子分布からの速度場の推定に起因するサンプリングアーティファクトが、$k \sim 0.2\text{--}0.3 \, h\text{Mpc}^{-1}$ といったスケールでの測定を信頼できなくする解像度のシミュレーションに依存していた。
2. 人工的な $h$ 依存性： B19 は、非線形進化を $\sigma_8$（$8 \, h^{-1}\text{Mpc}$ で平滑化された密度揺らぎの RMS）を用いてパラメータ化している。平滑化スケールがハッブルパラメータ $h$ に依存するため、これは $h$ に対する非物理的な依存性を導入する。その結果、線形物理学は同一だが $h$ の値が異なるモデルは発散する予測を生み出し、宇宙論パラメータ空間全体への外挿時の精度を損なう。
3. 関数形： 既存のモデルは、追加の自由パラメータなしには、モード結合効果やバリオン音響振動（BAO）の非線形減衰を完全に捉えていない。

手法
著者らは、進化マッピングフレームワークを用いて、重力のみを扱う $N$ 体シミュレーション（Aletheia、AletheiaMass、AletheiaEmu）のスイート上で較正された新しいフィッティング関数を提示する。このフレームワークは、形状パラメータ（$\Theta_s$）が固定されている場合、進化パラメータ（$\Theta_e$）と赤方偏移の影響は、$h$ に依存しない単一のクラスターリング振幅パラメータ $\sigma_{12}$（12 Mpc で平滑化された密度揺らぎの RMS）に吸収され得ると仮定している。

主要な手法のステップには以下が含まれる：

• 収束戦略： 著者らは、マルチ解像度シミュレーション（異なるボックスサイズと粒子数）を特徴とする AletheiaMass スイートを利用する。異なる体積からの測定値を継ぎ接ぎする「パッチング」戦略を採用する。各ボックスに対して、高解像度と低解像度のランを比較することで収束スケール $k_{\text{crop}}$ を定義し、最終的なパワースペクトルは $k < k_{\text{crop}}$ を満たす最大の体積からの測定値を使用する。これにより、$k \simeq 0.56 \, \text{Mpc}^{-1}$ までの収束した測定値が得られる。
• 速度推定： サンプリングアーティファクトを軽減するため、著者らは MC-Voronoi 法を使用する。これは、最も近い $N$ 体粒子に基づき、体積重み付けされた一様なガラス状の点分布 onto 速度場を再サンプリングし、その後 Cloud-in-Cell (CIC) 割当てとデコンボリューションを行うものである。
• 関数形：
  • $P_{\theta\theta}(k)$： 著者らは、非線形 BAO 減衰を考慮するためにデ・ウィグド線形パワースペクトル $P_{dw}$ を用いた減衰線形項と、モード結合を捉えるための $\ln k$ における加法的ガウス項を組み合わせたアンサッツを提案する。パラメータは $\sigma_{12}$ の関数である。
  • $P_{\delta\theta}(k)$： これは、非線形密度スペクトル $P_{\delta\delta}(k)$（Aletheia エミュレータによって予測される）と新しい $P_{\theta\theta}(k)$ の幾何平均を、$\sigma_{12}$ に依存する指数減衰項を乗じたものとしてモデル化される。
• 較正： パラメータは、広い範囲の $\sigma_{12}$（0.2 から 1.0）と赤方偏移（$z \approx -0.4$ から 4.2）をカバーする AletheiaMass スイートに対して、ベイズ的フレームワーク（nautilus）を用いてフィッティングされる。

主要な貢献

1. 新しいフィッティング関数： 本論文は、$\sigma_{12}$ に明示的に依存し、物理単位（Mpc$^{-1}$）を使用する $P_{\theta\theta}(k)$ と $P_{\delta\theta}(k)$ の較正済み式を提供し、人工的な $h$ 依存性を排除する。
2. 堅牢な収束性： マルチ解像度シミュレーションを用いた保守的なパッチング戦略を実装することで、著者らは $k \simeq 0.56 \, \text{Mpc}^{-1}$ までパーセントレベルの精度を達成し、以前の研究と比較して信頼できる範囲を大幅に拡大した。
3. 進化マッピングの検証： 本研究は、速度統計を $\sigma_{12}$ の観点で表現することが進化マッピング関係を維持することを確認し、線形物理学と $\sigma_{12}$ が同一であれば、その成長履歴や $h$ の値に関わらず、一貫した非線形予測が得られることを保証する。

結果

• 精度： 新しい処方は、較正セット内で広い範囲の $\sigma_{12}$ と $k \simeq 0.56 \, \text{Mpc}^{-1}$ までのスケールにおいて、$P_{\theta\theta}(k)$ と $P_{\delta\theta}(k)$ 双方で 1–2% の精度 を達成する。
• B19 との比較： 新しいモデルは B19 を体系的に上回る。B19 の予測は、解像度の制限により小スケール（$k > 0.2 \, \text{Mpc}^{-1}$）で大きく逸脱し、異なる $h$ の値を持つ宇宙論に対して非物理的なオフセットを示す。新しいモデルは、異なるハッブルパラメータにわたる速度統計を正しく予測する。
• 成長履歴の影響： 形状パラメータを固定し進化パラメータを変化させることで成長履歴の影響を分離する Aletheia スイートを用いて、著者らは完全な進化マッピングからの逸脱を定量化する。実用上の関心を持つほとんどの宇宙論において、成長履歴の差異を無視することによる誤差は 1% 未満 であることを発見する。有意な逸脱（>1%）は、極端な観測的に排除された宇宙論（例：アインシュタイン・ド・ジッター宇宙や極端な $h$ の値）でのみ発生する。
• 交差スペクトル： $P_{\delta\theta}(k)$ モデルは、$P_{\delta\delta}(k)$ に対する Aletheia エミュレータの恩恵を受け、これは暗黙的に成長履歴の依存性の一部を捉えており、誤差をさらに低減する。

重要性
本論文は、パラメータ化と単位的选择が速度統計モデルの堅牢性にとって決定的であると主張する。$h$ 依存単位（$\sigma_8$ など）が系統的バイアスを導入することを実証することで、著者らは物理単位と主要なスケーリングパラメータとしての $\sigma_{12}$ の使用に対する強力な実証的支援を提供する。これらのフィッティング関数は、現在および将来の銀河赤方偏移サーベイにおける RSD モデリングに直接適用可能な、堅牢でパーセントレベルの速度パワースペクトルの記述を提供する。著者らは、現在のフィッティング関数が実用的な用途には十分であるが、将来の研究には、残留する成長履歴の逸脱を明示的にモデル化し、宇宙論パラメータ空間全体でより高い精度を達成するための専用速度エミュレータの構築が含まれる可能性があることに言及している。