Search

Large-scale structure surveys measure the shape and position of millions of galaxies in order to constrain the cosmological model with high precision. The resulting large data volume poses a challenge for the analysis of the data, from the estimation of photometric redshifts to the calibration of shape measurements. We present GalSBI, a model for the galaxy population, to address these challenges. This phenomenological model is constrained by observational data using simulation-based inference (SBI). The $\texttt{galsbi}$ Python package provides an easy interface to generate catalogs of galaxies based on the GalSBI model, including their photometric properties, and to simulate realistic images of these galaxies using the $\texttt{UFig}$ package., Comment: 10 pages, 5 figures, submitted to JOSS

With the rise of simulation-based inference (SBI) methods, simulations need to be fast as well as realistic. $\texttt{UFig v1}$ is a public Python package that generates simulated astronomical images with exceptional speed - taking approximately the same time as source extraction. This makes it particularly well-suited for simulation-based inference (SBI) methods where computational efficiency is crucial. To render an image, $\texttt{UFig}$ requires a galaxy catalog and a description of the point spread function (PSF). It can also add background noise, sample stars using the Besan\c{c}on model of the Milky Way, and run $\texttt{SExtractor}$ to extract sources from the rendered image. The extracted sources can be matched to the intrinsic catalog, flagged based on $\texttt{SExtractor}$ output and survey masks, and emulators can be used to bypass the image simulation and extraction steps. A first version of $\texttt{UFig}$ was presented in Berg\'e et al. (2013) and the software has since been used and further developed in a variety of forward modelling applications., Comment: 8 pages, 4 figures, submitted to JOSS

We present GalSBI, a phenomenological model of the galaxy population for cosmological applications using simulation-based inference. The model is based on analytical parametrizations of galaxy luminosity functions, morphologies and spectral energy distributions. Model constraints are derived through iterative Approximate Bayesian Computation, by comparing Hyper Suprime-Cam deep field images with simulations which include a forward model of instrumental, observational and source extraction effects. We developed an emulator trained on image simulations using a normalizing flow. We use it to accelerate the inference by predicting detection probabilities, including blending effects and photometric properties of each object, while accounting for background and PSF variations. This enables robustness tests for all elements of the forward model and the inference. The model demonstrates excellent performance when comparing photometric properties from simulations with observed imaging data for key parameters such as magnitudes, colors and sizes. The redshift distribution of simulated galaxies agrees well with high-precision photometric redshifts in the COSMOS field within $1.5\sigma$ for all magnitude cuts. Additionally, we demonstrate how GalSBI's redshifts can be utilized for splitting galaxy catalogs into tomographic bins, highlighting its potential for current and upcoming surveys. GalSBI is fully open-source, with the accompanying Python package, $\texttt{galsbi}$, offering an easy interface to quickly generate realistic, survey-independent galaxy catalogs., Comment: 29 pages, 14 figures, submitted to Phys. Rev. D

Accurate redshift calibration is required to obtain unbiased cosmological information from large-scale galaxy surveys. In a forward modelling approach, the redshift distribution n(z) of a galaxy sample is measured using a parametric galaxy population model constrained by observations. We use a model that captures the redshift evolution of the galaxy luminosity functions, colours, and morphology, for red and blue samples. We constrain this model via simulation-based inference, using factorized Approximate Bayesian Computation (ABC) at the image level. We apply this framework to HSC deep field images, complemented with photometric redshifts from COSMOS2020. The simulated telescope images include realistic observational and instrumental effects. By applying the same processing and selection to real data and simulations, we obtain a sample of n(z) distributions from the ABC posterior. The photometric properties of the simulated galaxies are in good agreement with those from the real data, including magnitude, colour and redshift joint distributions. We compare the posterior n(z) from our simulations to the COSMOS2020 redshift distributions obtained via template fitting photometric data spanning the wavelength range from UV to IR. We mitigate sample variance in COSMOS by applying a reweighting technique. We thus obtain a good agreement between the simulated and observed redshift distributions, with a difference in the mean at the 1$\sigma$ level up to a magnitude of 24 in the i band. We discuss how our forward model can be applied to current and future surveys and be further extended. The ABC posterior and further material will be made publicly available at https://cosmology.ethz.ch/research/software-lab/ufig.html., Comment: Accepted for publication in JCAP. 35 pages, 9 figures

Intrinsic alignment (IA) modelling and photometric redshift estimation are two of the main sources of systematic uncertainty in weak lensing surveys. We investigate the impact of redshift errors and their interplay with different IA models. Generally, errors on the mean $\delta_z$ and on the width $\sigma_z$ of the redshift bins can both lead to biases in cosmological constraints. We find that such biases can, however, only be partially resolved by marginalizing over $\delta_z$ and $\sigma_z$. For Stage-III surveys, $\delta_z$ and $\sigma_z$ cannot be well constrained due to limited statistics. The resulting biases are thus sensitive to prior volume effects. For Stage-IV surveys, we observe that marginalizing over the redshift parameters has an impact and reduces the bias. We derive requirements on the uncertainty of $\sigma_z$ and $\delta_z$ for both Stage-III and Stage-IV surveys. We assume that the redshift systematic errors on $S_8$ should be less than half of the statistical errors, and the median bias should be smaller than $0.25\sigma$. We find that the uncertainty on $\delta_z$ has to be $\lesssim0.025$ for the NLA IA model with a Stage-III survey. For $\sigma_z$, the requirement is met even for large uncertainties $\leq0.3$. For the TATT IA model, the uncertainty on $\delta_z$ has to be $\lesssim0.02$ and the uncertainty on $\sigma_z$ has to be $\lesssim0.2$. For Stage-IV surveys, the uncertainty on $\delta_z$ has to be $\lesssim0.005$ and the uncertainty on $\sigma_z$ should be $\lesssim0.1$, with no significant dependence on the IA model. This required high precision will be a challenge for the redshift calibration of these future surveys. Finally, we investigate whether the interplay between redshift systematics and IA modelling can explain the $S_8$-tension between cosmic shear results and CMB measurements. We find that this is unlikely to explain the current $S_8$-tension., Comment: 37 pages including 15 figures and 9 tables, accepted for publication in JCAP

The next generation of weak lensing surveys will measure the matter distribution of the local Universe with unprecedented precision, allowing the resolution of non-Gaussian features of the convergence field. This encourages the use of higher-order mass-map statistics for cosmological parameter inference. We extend the forward-modelling based methodology introduced in a previous forecast paper to match these new requirements. We provide multiple forecasts for the wCDM parameter constraints that can be expected from stage 3 and 4 weak lensing surveys. We consider different survey setups, summary statistics and mass map filters including wavelets. We take into account the shear bias, photometric redshift uncertainties and intrinsic alignment. The impact of baryons is investigated and the necessary scale cuts are applied. We compare the angular power spectrum analysis to peak and minima counts as well as Minkowski functionals of the mass maps. We find a preference for Starlet over Gaussian filters. Our results suggest that using a survey setup with 10 instead of 5 tomographic redshift bins is beneficial. Adding cross-tomographic information improves the constraints on cosmology and especially on galaxy intrinsic alignment for all statistics. In terms of constraining power, we find the angular power spectrum and the peak counts to be equally matched for stage 4 surveys, followed by minima counts and the Minkowski functionals. Combining different summary statistics significantly improves the constraints and compensates the stringent scale cuts. We identify the most `cost-effective' combination to be the angular power spectrum, peak counts and Minkowski functionals following Starlet filtering., Comment: Monthly Notices of the Royal Astronomical Society, 28 July 2023