Abstract:The non-Gaussisan spatial distribution of galaxies traces the large-scale structure of the Universe and therefore constitutes a prime observable to constrain cosmological parameters. We conduct Bayesian inference of the $\Lambda$CDM parameters $\Omega_m$, $\Omega_b$, $h$, $n_s$, and $\sigma_8$ from the BOSS CMASS galaxy sample by combining the wavelet scattering transform (WST) with a simulation-based inference approach enabled by the ${\rm S{\scriptsize IM}BIG}$ forward model. We design a set of reduced WST statistics that leverage symmetries of redshift-space data. Posterior distributions are estimated with a conditional normalizing flow trained on 20,000 simulated ${\rm S{\scriptsize IM}BIG}$ galaxy catalogs with survey realism. We assess the accuracy of the posterior estimates using simulation-based calibration and quantify generalization and robustness to the change of forward model using a suite of 2,000 test simulations. When probing scales down to $k_{\rm max}=0.5~h/\text{Mpc}$, we are able to derive accurate posterior estimates that are robust to the change of forward model for all parameters, except $\sigma_8$. We mitigate the robustness issues with $\sigma_8$ by removing the WST coefficients that probe scales smaller than $k \sim 0.3~h/\text{Mpc}$. Applied to the BOSS CMASS sample, our WST analysis yields seemingly improved constraints obtained from a standard PT-based power spectrum analysis with $k_{\rm max}=0.25~h/\text{Mpc}$ for all parameters except $h$. However, we still raise concerns on these results. The observational predictions significantly vary across different normalizing flow architectures, which we interpret as a form of model misspecification. This highlights a key challenge for forward modeling approaches when using summary statistics that are sensitive to detailed model-specific or observational imprints on galaxy clustering.

What problem does this paper attempt to address?

The paper aims to address key issues in cosmological parameter estimation, particularly by improving constraints on the Λ Cold Dark Matter (ΛCDM) model parameters through the use of the Wavelet Scattering Transform (WST). Specifically, the research objectives include: 1. **Improving Cosmological Parameter Estimation**: By combining the Simulation-Based Inference (SimBIG) method and the Wavelet Scattering Transform (WST), infer key parameters of the ΛCDM model (such as the matter density parameter $\Omega_m$, baryon density parameter $\Omega_b$, Hubble parameter $h$, scalar index $n_s$, and amplitude parameter $\sigma_8$) from the BOSS CMASS galaxy sample. 2. **Utilizing Non-Gaussian Information**: Use WST to capture non-Gaussian information arising from nonlinear dynamical effects in the galaxy distribution, which is difficult to obtain through traditional power spectrum analysis. 3. **Application of Simulation-Based Inference Methods**: Conduct Simulation-Based Inference (SBI) using a large number of numerical simulations (generated by the SimBIG framework), avoiding strong assumptions about the form of the likelihood function and handling complex nonlinear effects. 4. **Validation and Robustness Assessment**: Validate the proposed inference pipeline's effectiveness and robustness to variations in different forward models through a series of test simulations (including TEST0, TEST1, and TEST2 datasets). In summary, the core objective of this paper is to demonstrate how combining WST with SBI techniques can improve the statistical description of galaxy clustering, leading to more precise cosmological parameter estimates. Additionally, the research explores the performance and potential limitations of these techniques in different simulation environments.

Galaxy Clustering Analysis with SimBIG and the Wavelet Scattering Transform