A new overdensity-conditioned emulator trained on small subvolumes from Quijote recovers the global halo mass function via integration over the overdensity distribution at 0.026% of the simulation cost.
Simulation-based cosmological inference from optically-selected galaxy clusters with $\texttt{Capish}$
3 Pith papers cite this work. Polarity classification is still indexing.
abstract
Galaxy clusters are powerful probes of the growth of cosmic structure through measurements of their abundance as a function of mass and redshift. Extracting precise cosmological constraints from cluster surveys is challenging, as we must contend the complex relationship between richness and the underlying halo mass, selection function biases, super-sample covariance, and correlated measurement noise between mass proxies. As upcoming photometric surveys are expected to detect tens to hundreds of thousands of galaxy clusters, controlling these systematics becomes essential. In this paper, we present a forward-modelling approach using Simulation-Based Inference (SBI), which provides a natural framework for jointly modelling cluster abundance and lensing mass observables while capturing systematic uncertainties at higher fidelity than analytic likelihood methods - which rely on simplifying assumptions such as fixed covariances and Gaussianity - without requiring an explicit likelihood formulation. We introduce $\texttt{Capish}$, a Python code for generating forward-modelled galaxy cluster catalogues using halo mass functions and incorporating observational effects. We perform SBI using neural density estimation with normalizing flows, trained on abundance and mean lensing mass measurements in observed redshift-richness bins. Our forward model accounts for realistic noise, redshift uncertainties, selection functions, and correlated scatter between lensing mass and observed richness. We find good agreement with likelihood-based analyses, with broader SBI posteriors reflecting the increased realism of the forward model. We also test $\texttt{Capish}$ on cluster catalogues built from a large cosmological simulation, finding a good fit to cosmological parameters.
years
2026 3representative citing papers
Fixing the covariance at an incorrect cosmology in cluster count analyses leaves Ω_c, σ_8, and w estimates unbiased but distorts their uncertainties, driven by S_8 amplitude effects; a single update at the recovered best-fit cosmology restores correct normalization for LSST-like surveys.
Gaia quasar proper motions give a Solar System acceleration amplitude of 5.72 μas yr^{-1} with vector components (0.40, -5.09, -2.40), consistent with prior work but with 1.5-2.5 times larger credible intervals after marginalizing higher multipoles.
citing papers explorer
-
Efficiently emulating distribution functions in gigaparsec volumes for varying cosmological parameters
A new overdensity-conditioned emulator trained on small subvolumes from Quijote recovers the global halo mass function via integration over the overdensity distribution at 0.026% of the simulation cost.
-
Cosmology-dependent covariance in galaxy cluster number counts: consequences for parameter inference
Fixing the covariance at an incorrect cosmology in cluster count analyses leaves Ω_c, σ_8, and w estimates unbiased but distorts their uncertainties, driven by S_8 amplitude effects; a single update at the recovered best-fit cosmology restores correct normalization for LSST-like surveys.
-
Observing the Peculiar Acceleration of our Solar System with Quasar Proper Motions
Gaia quasar proper motions give a Solar System acceleration amplitude of 5.72 μas yr^{-1} with vector components (0.40, -5.09, -2.40), consistent with prior work but with 1.5-2.5 times larger credible intervals after marginalizing higher multipoles.