Balancing bias, baryons, and scale cuts in LSST 3x2pt analysis

技術要約：LSST $3 \times 2$pt 解析におけるバイアス、バリオン、およびスケールカットのバランス調整

問題提起
LSST（ルビン天文台 レガシー・サーベイ・オブ・スペース・アンド・タイム）のようなステージIVの宇宙論的サーベイは、銀河バイアスやバリオンフィードバックによる系統誤差が支配的となる非線形構造形成の領域を調査することになる。$3 \times 2$pt解析（銀河クラスタリング、宇宙剪断、およびそれらの相互相関を組み合わせたもの）において、これらの効果は、しばしば制約が弱く、宇宙論パラメータとの縮退を示す可能性のある摂動パラメータを導入する。具体的には、非線形銀河バイアスと物質パワースペクトルのバリオン抑制との相互作用が課題となる。純粋な線形バイアスモデルは小スケール（$k \gtrsim 0.1 \, h/\text{Mpc}$）では不十分であるが、より複雑な摂動モデルは、バリオン物理学と縮退する可能性のある追加の自由度をもたらす。対処すべき中心的な問題は、宇宙論パラメータ（$\Omega_m, \sigma_8$）および拡張パラメータ（ニュートリノ質量 $M_\nu$）に対する制約力を最大化しつつ、重大なバイアスを導入しないように、モデルの複雑さとスケールカットをどのようにバランスさせるかである。

手法
著者らは、MGL（Modified Gravity Lensing）と呼ばれる新しいオープンソースのベイズ推論パイプラインを用いた包括的な解析を提示している。この研究では、宇宙剪断と銀河クラスタリングのスペクトルを含む、LSSTの第1年（Y1）および第10年（Y10）の仕様に類似したモックデータベクトルを生成し、20個の対数間隔のビン $\ell \in [20, 5000]$ にわたって解析を行っている。

主な手法構成要素は以下の通りである：

• モデリング・フレームワーク： 解析では、非線形物質パワースペクトルと「バリオニフィケーション（baryonification）」メカニズムを介したバリオンフィードバックの両方をモデル化するために、BACCOエミュレータを採用している。銀河バイアスについては、以下の3つのモデルを比較している：
  1. 線形バイアス（Linear Bias）： 大スケールでのみ有効。
  2. ハイブリッド有効場理論（HEFT）： ラグランジュ的なバイアス展開とN体シミュレーションの進化を組み合わせた摂動的アプローチであり、$k_{\text{max}} \approx 0.7 \, h/\text{Mpc}$ まで有効。
  3. 最小バイアス（Minimal Bias）： 高次のバイスターム（$b_{s^2}$ および $b_{\nabla^2}$）をゼロに固定した簡略化されたHEFTの変種。
• バリオンの取り扱い： バリオンフィードバックは、物質パワースペクトルに適用される抑制因子 $S_{mm}$ を用いてモデル化される。銀河・銀河重力レンズ（$P_{gm}$）については、Zennaro et al. (2024) の近似を採用し、$P_{gm}$ は $\sqrt{S_{mm}}$ によって補正される一方、$P_{gg}$ はバイアス展開の柔軟性に依存する。
• 解析戦略： 著者らは、フルHEFTモデルとフィダル（基準）パラメータを用いて生成されたモックデータに対して尤度解析を行う。これには、不一致モデル（例：HEFTで生成されたデータを線形または最小バイアスモデルで解析する場合）の採用や、バリオンパラメータの包含の有無を変えるなど、様々な解析上の仮定のテストが含まれる。
• 指標： パフォーマンスは、制約力としてのFigure of Merit (FoM) と、フィダル値からの偏差を測定するFigure of Bias (FoB) を用いて定量化される。

主要な貢献

1. MGLパイプラインの開発： 本論文は、パワースペクトルの効率的な計算のためにエミュレータを利用し、nautilus サンプラーを用いたMCMC解析をサポートする、ステージIVサーベイ向けに設計された独立したオープンソース・パイプラインであるMGLを紹介している。
2. バイアスモデルの系統的な比較： 本研究は、LSSTのY1およびY10の文脈において、線形、HEFT、および最小バイアスモデルを直接比較し、特にバリオンフィードバックとスケールカットに対する頑健性を評価している。
3. 縮退解析： 高次のバイアスパラメータとバリオンフィードバックの間の縮退、およびバイアスモデリングと質量ニュートリノ検出との相互作用を明示的に調査している。

結果

• スケールカットとバイアスモデル： 線形バイアスモデルは、$k_{\text{max}} \approx 0.1 \, h/\text{Mpc}$ までは $\Omega_m$ および $\sigma_8$ に関して偏りのない制約を与える（FoB $< 0.5\sigma$）。より小さなスケールへの拡張には摂動的アプローチが必要となる。
• 最小バイアスの有効性： 「最小バイアス」モデル（$b_{s^2}$ と $b_{\nabla^2}$ をゼロに固定）は、フルHEFTモデルと同等の性能を示す。フィダルなバイアスパラメータが標準的な範囲内であれば、$k_{\text{max}} = 0.7 \, h/\text{Mpc}$ においても $\Omega_m$ および $\sigma_8$ に対して偏りなく機能する。これは、標準的な宇宙論的制約のためには、フル形式のHEFTの複雑さが必ずしも必要ではないことを示唆しており、パラメータ空間を削減する道筋を示している。
• バリオンとバイアスの相互作用： 高次のバイアスパラメータは、銀河・銀河パワースペクトルにおけるバリオン抑制の効果を模倣することができる。しかし、バリオンフィードバックは高次のバイアス項によって生成される挙動の範囲を完全には再現できない。銀河・銀河重力レンズ成分においてバリオン補正を無視しても、$\Omega_m$ や $\sigma_8$ に重大なバイアスを与えることはないが、推論されるバリオンパラメータ（$\log_{10} M_c$）をシフトさせる可能性がある。
• ニュートリノ質量の検出： 非線形スケール（$k \geq 0.3 \, h/\text{Mpc}$）まで拡張した場合、Y1およびY10の両方で全ニュートリノ質量（$M_\nu$）の検出が可能である。しかし、推論される $M_\nu$ の具体的な値は、あらゆる妥当なシナリオにおいてロバストではない。高次のバイアスパラメータが増大したモックデータを最小バイアスモデルで解析すると、有意な $M_\nu$ の偏り（Y10で最大 $\sim 2\sigma$）が生じる。これは、ニュートリノ質量、高次バイアス項、および投影効果の間の強い縮退を示している。
• Y1 vs. Y10： Y1からY10へ移行することで、$k_{\text{max}} = 0.7 \, h/\text{Mpc}$ においてFoMは約2.3倍向上する。$\Omega_m$ と $\sigma_8$ は頑健であるが、簡略化されたバイアスモデル下での $M_\nu$ 制約の脆弱性は依然として存在する。

意義と主張
本論文は、標準的な $\Lambda$CDMパラメータ（$\Omega_m, \sigma_8$）については、最小のバイアス・パラメトリゼーションを用いることで、HEFTの有効限界（$k_{\text{max}} = 0.7 \, h/\text{Mpc}$）までのLSSTデータから偏りのない制約を引き出すことが可能であると主張している。これにより、15項におよぶフル形式のバイアス展開に伴う計算的・統計的なペナルティなしに、準非線形スケールの包含が可能となる。

しかし、著者らは、この結論がバイアスパラメータの特定のフィダル値に依存していることに慎重な注意を払っている。LSSTは非ゼロのニュートリノ質量の明確な検出が可能であるはずだが、その検出の具体的な有意性と値は、バイアスモデリングの正確性に強く依存する。本研究は、簡略化されたモデリング手法（最小バイアスなど）が標準的な宇宙論に対しては頑健な結果をもたらす可能性がある一方で、基礎となる銀河バイアスの物理が想定されたフィダル値から逸脱している場合、ニュートリノ質量のような拡張パラメータに対しては重大なバイアスを導入する可能性があることを強調している。したがって、著者らは、採用前に特定の科学的ケースに対して簡略化されたモデルを注意深くテストすることを推奨している。