Cosmic variance or galaxy bias? Disentangling finite-volume and galaxy formation effects in cosmological analysis

技術要約：宇宙論的解析における有限体積効果と銀河形成効果の分離

問題提起
現在および今後の大規模銀河赤方偏移サーベイ（DESI、Euclidなど）は、前例のない精度で宇宙論パラメータを決定することを目指している。しかし、偏りのない結果を得るためには、理論的な系統誤差を厳密に制御する必要がある。主要な課題は、これら二つの強く結合した不確実性の源を区別することである：すなわち、有限の体積内における初期ガウス型ランダム場の固有の分散に起因する**宇宙論的分散（cosmic variance）と、銀河形成の複雑な非線形プロセスを記述する銀河バイアス（galaxy bias）**である。

標準的な摂動論的枠組みでは、バイアスパラメータは観測データをアンサンブル平均された物質パワースペクトルに適合させることで推定される。この手法は、調査体積がアンサンブル平均を代表していることを暗黙のうちに仮定しており、観測された特定の宇宙の実現（realization）に伴う宇宙論的分散を無視している。その結果、推定されるバイアスパラメータは退化し、真の銀河形成物理学と、調査された特定の宇宙体積の統計的ゆらぎが混同されてしまう。この退化は、バリオン音響振動（BAO）再構成や、バイアス値に関する仮定が重要となる宇宙論パラメータの推論において、問題を複雑化させる。

手法
著者らは、宇宙論的分散と銀河バイアスを数学的に分離するために、**オイラー摂動論（EPT）**に基づく新しい定式化を開発した。

1. 宇宙論的分散のバイアス化： 著者らは、任意の有限体積の実現（$\delta'$）を、アンサンブル平均された密度場（$\delta$）の「バイアスされた」トレーサーとして扱う手法を提案している。新しい一連のバイアスパラメータ（$b_{1,dm}, b_{2,dm}$）を導入することで、有限体積の場を以下のように表現する：
   $$\delta'(x) = b_{0,dm} + b_{1,dm}\delta(x) + \frac{1}{2}b_{2,dm}\delta^2(x) + \dots$$
   この枠組みにおいて、体積が無限大に近づくにつれ、$b_{1,dm} \to 1$ および $b_{2,dm} \to 0$ となる。この定式化により、「宇宙論的分散によるバイアス」と、物理的な「銀河バイアス」（$b_{1,g}, b_{2,g}$）の分離が可能になる。

2. BAOシフトへの適用： 著者らは、この定式化をBAOスケールの非線形シフトに適用した。［22］によって導出された解析的な式を用いて、全BAOシフト（$\Delta \alpha$）を以下の二つの異なる成分に分解した：

   • ダークマター場の固有の宇宙論的分散によって駆動される成分（$b_{1,dm}, b_{2,dm}$ によって制御される）。
   • 銀河バイアスによって駆動される成分（$b_{1,g}, b_{2,g}$ によって制御される）。
     この分解により、宇宙論的分散単独でも、特に小さな体積においては、銀河バイアスによるシフトに匹敵するBAOシフトを引き起こし得ることが明らかになり、これは非線形重力成長によるシフトと退化している。

3. 積分制約の処理： 本論文では、調査体積の平均密度をゼロに強制することによって生じるバイアス（積分制約）に対処している。著者らは、この効果がバイアス展開の定数項（$b'_0$）に吸収可能であることを示し、それが追加のバイアスパラメータを導入したり、宇宙論的分散の分離を複雑にしたりしないことを証明している。

4. 数値的検証： 解析的な議論を検証するため、著者らはUNITシミュレーション・スイートを利用している。彼らは、4つの異なる$N$体シミュレーションのセット（$(1 h^{-1}\text{Gpc})^3$ 体積内の粒子数 $2 \times 10^3$）を分析した：

   • 100個の標準的なシミュレーション。
   • 100個の固定振幅シミュレーション（振幅の分散を抑制）。
   • 2 $\times$ 50組のペアリングされたシミュレーション（位相の分散を抑制）。
   • 2 $\times$ 50組の固定かつペアリングされたシミュレーション（両方を抑制）。
     彼らは、モデルに依存しない方法（残差相関関数にガウス関数をフィッティングする手法）でBAOピーク位置を測定し、分散減少を定量化した。

主な結果

• 解析的な分離： 新しい定式化は、非線形BAOシフトに対する宇宙論的分散と銀河バイアスの寄与を正常に分離した。解析の結果、宇宙論的分散は銀河バイアスと同等のシフトを引き起こし得ることを確認しており、これは有限体積における摂動論的アプローチを用いる際に細心の注意が必要であることを示唆している。
• 分散の抑制： UNITシミュレーションによる数値結果は、初期条件を制御することでBAOピーク位置の不確実性が大幅に減少することを実証している：
  • 固定（Fixing）（振幅を制御）：不確実性を約 $1.5 - 2.1$ 倍減少させる。
  • ペアリング（Pairing）（位相を制御）：不確実性を約 $1.5 - 1.7$ 倍減少させる。
  • 併用（Fixed-and-Paired）：これら両方の手法を組み合わせることで、標準的なアンサンブルと比較して約 $3.4 - 5.4$ 倍の分散抑制を実現する。
• 収束性： BAOピーク位置は、すべてのシミュレーションタイプにおいてサブパーセントレベルで良好に収束しており、固定またはペアリングの手法によって標準的なアンサンブルに対して明らかな系統誤差が生じることはない。

意義および提案される戦略
本論文は、現在のバイアスパラメータの推定手法が、宇宙論的分散の退化した効果を不当に含んでしまっていると結論付けている。これに対処するため、著者らは将来の宇宙論的解析に向けて二つの実用的な戦略を提案している：

1. 最適なマルチトレーサーの組み合わせ： 単に統計的な精度を最大化するのではなく、将来の解析では、宇宙論的分散を相殺するために補完的なトレーサー（例：銀河と宇宙空洞）を組み合わせることができる。これには、補完的な密度領域をトレースするように事前再構成サンプルを再定義するか、あるいは、ペアリングされた実現の（効果を）模倣するように、結合された相関関数の重みを最適化することなどが含まれる。
2. 分散抑制再構成： 著者らは、観測データ（逆転操作やベイズ的順方向再構成などの再構成手法を介して）および共分散推定用の模擬カタログの両方に、分散抑制技術（固定およびペアリング）を適用することを提案している。これにより、観測された宇宙とその対となる「反宇宙」の「分散抑制された」バージョンを解析することが可能になる。ただし、再構成プロセスで使用される宇宙論によってBAOシフトにバイアスが生じることを避けるため、慎重な実装が必要であると注記している。

著者らは、これらの戦略が有望であることを強調しており、詳細な模擬カタログ（DESIのようなもの）でのさらなるテストと、BAO再構成や有効場理論（EFT）といった最先端の解析手法への統合が必要であるとしている。