Inferring the photometric and size evolution of galaxies from image simulations

Abstract

Current constraints on models of galaxy evolution rely on morphometric catalogs extracted from multi-band photometric surveys. However, these catalogs are altered by selection effects that are difficult to model, that correlate in non trivial ways, and that can lead to contradictory predictions if not taken into account carefully. To address this issue, we have developed a new approach combining parametric Bayesian indirect likelihood (pBIL) techniques and empirical modeling with realistic image simulations that reproduce a large fraction of these selection effects. This allows us to perform a direct comparison between observed and simulated images and to infer robust constraints on model parameters. We use a semi-empirical forward model to generate a distribution of mock galaxies from a set of physical parameters. These galaxies are passed through an image simulator reproducing the instrumental characteristics of any survey and are then extracted in the same way as the observed data. The discrepancy between the simulated and observed data is quantified, and minimized with a custom sampling process based on adaptive Monte Carlo Markov Chain methods. Using synthetic data matching most of the properties of a CFHTLS Deep field, we demonstrate the robustness and internal consistency of our approach by inferring the parameters governing the size and luminosity functions and their evolutions for different realistic populations of galaxies. We also compare the results of our approach with those obtained from the classical spectral energy distribution fitting and photometric redshift approach.Our pipeline infers efficiently the luminosity and size distribution and evolution parameters with a very limited number of observables (3 photometric bands). When compared to SED fitting based on the same set of observables, our method yields results that are more accurate and free from systematic biases.

24 pages, 12 figures, accepted for publication in A&A

References in corpus (15)

• Gaia Data Release 1. Summary of the astrometric, photometric, and survey properties
• Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems
• Galaxy Zoo: the dependence of morphology and colour on environment
• Optimal whitening and decorrelation
• Constraining the thick disc formation scenario of the Milky Way
• The Kormendy relation of massive elliptical galaxies at z~1.5. Evidence for size evolution ?
• Bayesian Indirect Inference Using a Parametric Auxiliary Model
• astroABC: An Approximate Bayesian Computation Sequential Monte Carlo sampler for cosmological parameter estimation
• The number density of superdense early-type galaxies at 1<z<2 and the local cluster galaxies
• Approximate Bayesian Computation in Large Scale Structure: constraining the galaxy-halo connection
• Towards a fully consistent Milky Way disc model - III. Constraining the initial mass function
• On estimating redshift and luminosity distributions in photometric redshift surveys
• Fanaroff-Riley dichotomy of radio galaxies and the Malmquist bias
• The First Source Counts at 18 microns from the AKARI NEP Survey
• Simulating Deep Hubble Images With Semi-empirical Models of Galaxy Formation

Cited by in corpus (12)

• Massive optimal data compression and density estimation for scalable, likelihood-free inference in cosmology
• Likelihood-free inference with neural compression of DES SV weak lensing map statistics
• Nuisance hardened data compression for fast likelihood-free inference
• Automatic physical inference with information maximising neural networks
• Generalized massive optimal data compression
• Bayesian optimisation for likelihood-free cosmological inference
• Simulation-Based Inference of Reionization Parameters From 3D Tomographic 21 cm Lightcone Images
• Inference of the optical depth to reionization from low multipole temperature and polarisation Planck data
• Measurement of the B-band Galaxy Luminosity Function with Approximate Bayesian Computation
• Constraining the intrinsic population of Long Gamma-Ray Bursts: implications for spectral correlations, cosmic evolution and their use as tracers of star formation
• Catalog-free modeling of galaxy types in deep images: Massive dimensional reduction with neural networks
• GalSBI: Phenomenological galaxy population model for cosmology using simulation-based inference