[论文解读] Cosmology with the Galaxy Bispectrum Multipoles: Optimal Estimation and Application to BOSS Data

Mikhail M. Ivanov, Oliver H. E. Philcox+4 more|arXiv (Cornell University)|Feb 9, 2023

Cosmology and Gravitation Theories被引用 4

一句话总结

本文提出了一种新型的无窗最大似然估计器，用于测量红移空间大尺度结构中星系三谱多极矩（ℓ=2 和 ℓ=4），并在高保真度模拟中进行了验证，应用于BOSS DR12数据。结果表明，包含这些高阶多极矩可使一维后验分布的宇宙学约束提升5–10%，获得紧致的约束结果：$H_0 = 68.2 \pm 0.8~\mathrm{km\,s^{-1}\,Mpc^{-1}}$，$\Omega_m = 0.33 \pm 0.01$，$\sigma_8 = 0.736 \pm 0.033$，且 $S_8 = 0.77 \pm 0.04$ 与CMB和弱引力透镜观测结果一致。

ABSTRACT

We present a framework for self-consistent cosmological analyses of the full-shape anisotropic bispectrum, including the quadrupole $(\ell=2)$ and hexadecapole $(\ell=4)$ moments. This features a novel window-free algorithm for extracting the latter quantities from data, derived using a maximum-likelihood prescription. Furthermore, we introduce a theoretical model for the bispectrum multipoles (which does not introduce new free parameters), and test both aspects of the pipeline on several high-fidelity mocks, including the PT Challenge suite of gigantic cumulative volume. This establishes that the systematic error is significantly below the statistical threshold, both for the measurement and modeling. As a realistic example, we extract the large-scale bispectrum multipoles from BOSS DR12 and analyze them in combination with the power spectrum data. Assuming a minimal $Λ$CDM model, with a BBN prior on the baryon density and a extit{Planck} prior on $n_s$, we can extract the remaining cosmological parameters directly from the clustering data. The inclusion of the unwindowed higher-order $(\ell>0)$ large-scale bispectrum multipoles is found to moderately improve one-dimensional cosmological parameter posteriors (at the $5\%-10\%$ level), though these multipoles are detected only in three out of four BOSS data segments at $\approx 5σ$. Combining information from the power spectrum and bispectrum multipoles, the real space power spectrum, and the post-reconstructed BAO data, we find $H_0 = 68.2\pm 0.8~\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, $Ω_m =0.33\pm 0.01$ and $σ_8 = 0.736\pm 0.033$ (the tightest yet found in perturbative full-shape analyses). Our estimate of the growth parameter $S_8=0.77\pm 0.04$ agrees with both weak lensing and CMB results.

研究动机与目标

• 开发一种稳健、自洽的宇宙学分析框架，用于完整形状的各向异性星系三谱多极矩分析，包括四极矩（ℓ=2）和十六极矩（ℓ=4）多极矩。
• 通过引入一种新型无窗最大似然估计器，解决由观测窗函数和多极矩中模式混叠引起的系统误差。
• 在高保真度模拟上验证三谱多极矩的理论建模，包括566 $h^{-3}$Gpc$^3$ 体积的PT挑战模拟套件。
• 将该分析流程应用于真实BOSS DR12数据，评估除单极矩外高阶多极矩带来的宇宙学信息增益。
• 公开发布估计器和数据，以支持未来大尺度结构巡天的可重复性和可扩展分析。

提出的方法

• 采用最大似然方法推导出无窗的三谱多极矩估计器，无需对观测窗函数进行正向建模。
• 在N系列模拟和PT挑战模拟上测试该方法，验证其对系统误差的鲁棒性，并确保多极矩估计的偏差极低。
• 采用基于现有理论框架的微扰理论模型来描述三谱多极矩，扩展至包含ℓ=2和ℓ=4矩，且未引入新的自由参数。
• 在最小ΛCDM模型下，将功率谱、BAO和三谱多极矩（包括单极矩及更高阶矩）联合纳入似然分析框架。
• 通过BBN对重子密度施加先验，通过Planck对谱指数 $n_s$ 施加先验，以最小化对外部先验的依赖。
• 所有估计器和数据均通过GitHub公开发布，以确保可重复性，并促进未来巡天中的应用。

实验结果

研究问题

• RQ1 无窗最大似然估计器能否在不进行窗函数正向建模的情况下，可靠地从真实巡天数据中提取星系三谱四极矩和十六极矩？
• RQ2 与仅使用单极矩的分析相比，将ℓ=2和ℓ=4多极矩纳入分析是否能显著提升宇宙学参数的约束精度？
• RQ3 在高保真度模拟中测试时，测量和建模各向异性三谱多极矩的系统误差水平如何？
• RQ4 在最小ΛCDM模型下，与功率谱和BAO数据相比，完整形状三谱多极矩的宇宙学约束结果如何？
• RQ5 高阶三谱多极矩的信息是否与CMB和弱引力透镜等独立探针一致，特别是在$S_8$方面？

主要发现

• 无窗最大似然估计器在高保真度模拟中成功测量了星系三谱多极矩（ℓ=2 和 ℓ=4），系统误差低于统计阈值。
• 与仅使用单极矩的分析相比，包含ℓ>0多极矩可使一维宇宙学参数后验分布的约束提升5–10%，其中$H_0$、$\Omega_m$和$\sigma_8$的增益最为显著。
• 对BOSS DR12的分析得到：$H_0 = 68.2 \pm 0.8~\mathrm{km\,s^{-1}\,Mpc^{-1}}$，$\Omega_m = 0.33 \pm 0.01$，$\sigma_8 = 0.736 \pm 0.033$，这是迄今为止通过微扰理论完整形状分析获得的最紧约束。
• $S_8 = 0.77 \pm 0.04$ 与CMB和弱引力透镜结果一致，在1σ水平上缓解了$S_8$张力。
• 在四个BOSS数据段中的三个，三谱多极矩的检测显著性达到≈5σ，证实了其稳健性和信息含量。
• 整个分析流程，包括估计器和数据，均已公开发布，以支持可重复性，并为未来巡天应用提供支持。