Diagnosing Systematic Effects Using the Inferred Initial Power Spectrum

Abstract

The next generation of galaxy surveys has the potential to substantially deepen our understanding of the Universe. This potential hinges on our ability to rigorously address systematic uncertainties. Until now, diagnosing systematic effects prior to inferring cosmological parameters has been out of reach in field-based implicit likelihood cosmological inference frameworks. As a solution, we aim to diagnose a variety of systematic effects in galaxy surveys prior to inferring cosmological parameters, using the inferred initial matter power spectrum. Our approach is built upon a two-step framework. First, we employed the SELFI algorithm to infer the initial matter power spectrum, which we utilised to thoroughly investigate the impact of systematic effects. This investigation relies on a single set of $N$-body simulations. Second, we obtained a posterior on cosmological parameters via implicit likelihood inference, recycling the simulations from the first step for data compression. For demonstration, we relied on a model of large-scale spectroscopic galaxy surveys that incorporates fully non-linear gravitational evolution with COLA and simulates multiple systematic effects encountered in real surveys. We provide a practical guide on how the SELFI posterior can be used to assess the impact of misspecified galaxy bias parameters, selection functions, survey masks, inaccurate redshifts, and approximate gravity models on the inferred initial matter power spectrum. We show that a subtly misspecified model can lead to a bias exceeding $2\sigma$ in the $({\Omega }_{\mathrm{m}},{\sigma }_8)$ plane, which we are able to detect and avoid prior to inferring the cosmological parameters. This framework has the potential to significantly enhance the robustness of physical information extraction from full forward models of large-scale galaxy surveys such as DESI, Euclid, and LSST.

28 pages, 14 figures

References in corpus (87)

• Detection of the Baryon Acoustic Peak in the Large-Scale Correlation Function of SDSS Luminous Red Galaxies
• Cosmological parameters from SDSS and WMAP
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• The 2dF Galaxy Redshift Survey: Spectra and redshifts
• The Cosmic Linear Anisotropy Solving System (CLASS) II: Approximation schemes
• Baryonic Features in the Matter Transfer Function
• The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 10 and 11 galaxy samples
• Spectroscopic Target Selection for the Sloan Digital Sky Survey: The Luminous Red Galaxy Sample
• Why your model parameter confidences might be too optimistic -- unbiased estimation of the inverse covariance matrix
• Large-Scale Galaxy Bias
• Adaptive approximate Bayesian computation
• Euclid. I. Overview of the Euclid mission
• The 6dF Galaxy Survey: z \approx 0 measurement of the growth rate and sigma_8
• Solving Large Scale Structure in Ten Easy Steps with COLA
• Bayesian physical reconstruction of initial conditions from large scale structure surveys
• The Clustering of the SDSS Main Galaxy Sample II: Mock galaxy catalogues and a measurement of the growth of structure from Redshift Space Distortions at $z=0.15$
• Massive Lossless Data Compression and Multiple Parameter Estimation from Galaxy Spectra
• On Galactic density modeling in the presence of dust extinction
• Baryonic Signatures in Large-Scale Structure
• Target Selection and Validation of DESI Luminous Red Galaxies
• The power spectrum and bispectrum of SDSS DR11 BOSS galaxies I: bias and gravity
• Global, Exact Cosmic Microwave Background Data Analysis Using Gibbs Sampling
• COSMOPOWER: emulating cosmological power spectra for accelerated Bayesian inference from next-generation surveys
• Massive optimal data compression and density estimation for scalable, likelihood-free inference in cosmology
• ELUCID - Exploring the Local Universe with reConstructed Initial Density field I: Hamiltonian Markov Chain Monte Carlo Method with Particle Mesh Dynamics
• Physical Bayesian modelling of the non-linear matter distribution: new insights into the Nearby Universe
• The Clustering of Galaxies in SDSS-III DR9 Baryon Oscillation Spectroscopic Survey: Constraints on Primordial Non-Gaussianity
• Effects of Systematic Uncertainties on the Supernova Determination of Cosmological Parameters
• ELUCID - Exploring the Local Universe with reConstructed Initial Density field III: Constrained Simulation in the SDSS Volume
• The effect of peculiar velocities on supernova cosmology
• Joint analysis of 6dFGS and SDSS peculiar velocities for the growth rate of cosmic structure and tests of gravity
• Nuisance hardened data compression for fast likelihood-free inference
• Bayesian power-spectrum inference for Large Scale Structure data
• Automatic physical inference with information maximising neural networks
• Likelihood-Free Cosmological Inference with Type Ia Supernovae: Approximate Bayesian Computation for a Complete Treatment of Uncertainty
• Generalized massive optimal data compression
• Lensing, reddening and extinction effects of MgII absorbers from z=0.4 to z=2
• cosmoabc: Likelihood-free inference via Population Monte Carlo Approximate Bayesian Computation
• Systematic effects on dark energy from 3D weak shear
• Hierarchical Cosmic Shear Power Spectrum Inference
• Bayesian optimisation for likelihood-free cosmological inference
• Cosmic transparency: A test with the baryon acoustic feature and type Ia supernovae
• Calibration errors unleashed: effects on cosmological parameters and requirements for large-scale structure surveys
• On the impact of galaxy bias uncertainties on primordial non-Gaussianity constraints
• The Impact of Cosmic Dust on Supernova Cosmology
• Neutrino masses and beyond-$\Lambda$CDM cosmology with LSST and future CMB experiments
• Sloan Digital Sky Survey III Photometric Quasar Clustering: Probing the Initial Conditions of the Universe using the Largest Volume
• Exploiting the full potential of photometric quasar surveys: Optimal power spectra through blind mitigation of systematics
• Bayesian analysis of the dynamic cosmic web in the SDSS galaxy survey
• Galaxy bias from forward models: linear and second-order bias of IllustrisTNG galaxies
• Cosmological inference from Bayesian forward modelling of deep galaxy redshift surveys
• JAX-COSMO: An End-to-End Differentiable and GPU Accelerated Cosmology Library
• How much information can be extracted from galaxy clustering at the field level?
• FlowPM: Distributed TensorFlow Implementation of the FastPM Cosmological N-body Solver
• Discrepancies between CFHTLenS cosmic shear & Planck: new physics or systematic effects?
• Cosmic Shear: Inference from Forward Models
• LtU-ILI: An All-in-One Framework for Implicit Inference in Astrophysics and Cosmology
• Consistency tests of field level inference with the EFT likelihood
• Kernel-Based Emulator for the 3D Matter Power Spectrum from CLASS
• Interloper bias in future large-scale structure surveys
• Galaxy Clustering Analysis with SimBIG and the Wavelet Scattering Transform
• Impact of systematics on cosmological parameters from future Galaxy Clusters surveys
• Impacts of the physical data model on the forward inference of initial conditions from biased tracers
• Galaxy Bias and ${\sigma }_8$ from Counts in Cells from the SDSS Main Sample
• SED-dependent Galactic Extinction Prescription for Euclid and Future Cosmological Surveys
• EFTofLSS meets simulation-based inference: ${\sigma }_8$ from biased tracers
• Bayesian Inference of Initial Conditions from Non-Linear Cosmic Structures using Field-Level Emulators
• Explicit Bayesian treatment of unknown foreground contaminations in galaxy surveys
• The Impact of Line Misidentification on Cosmological Constraints from Euclid and other Spectroscopic Galaxy Surveys
• Primordial power spectrum and cosmology from black-box galaxy surveys
• Growth-rate measurement with type-Ia supernovae using ZTF survey simulations
• Bayesian power-spectrum inference with foreground and target contamination treatment
• Almanac: Weak Lensing power spectra and map inference on the masked sphere
• The Effect of Systematic Redshift Biases in BAO Cosmology
• Compensated Isocurvature Perturbations in the Galaxy Power Spectrum
• Field-Based Physical Inference From Peculiar Velocity Tracers
• Line confusion in spectroscopic surveys and its possible effects: Shifts in Baryon Acoustic Oscillations position
• Accurate field-level weak lensing inference for precision cosmology
• UNIONS: The impact of systematic errors on weak-lensing peak counts
• DISCO-DJ I: a differentiable Einstein-Boltzmann solver for cosmology
• Simulation-Based Inference Benchmark for Weak Lensing Cosmology
• $\text{PineTree}$: A generative, fast, and differentiable halo model for wide-field galaxy surveys
• SIDE-real: Supernova Ia Dust Extinction with truncated marginal neural ratio estimation applied to real data
• Almanac: MCMC-based signal extraction of power spectra and maps on the sphere
• Beyond dark energy Fisher forecasts: how DESI will constrain LCDM and quintessence models
• Joint constraints on galaxy bias and ${\sigma }_8$ through the N-pdf of the galaxy number density
• Bye-bye, Local-in-matter-density Bias: The Statistics of the Halo Field Are Poorly Determined by the Local Mass Density