Exploiting the full statistical power of future cosmic shear surveys will necessitate improvements to the accuracy with which the gravitational lensing signal is measured. We present a framework for ...calibrating shear with image simulations that demonstrates the importance of including realistic correlations between galaxy morphology, size, and more importantly, photometric redshifts. This realism is essential to ensure that selection and shape measurement biases can be calibrated accurately for a tomographic cosmic shear analysis. We emulate Kilo-Degree Survey (KiDS) observations of the COSMOS field using morphological information from Hubble Space Telescope imaging, faithfully reproducing the measured galaxy properties from KiDS observations of the same field. We calibrate our shear measurements from lensfit, and find through a range of sensitivity tests that lensfit is robust and unbiased within the allowed two per cent tolerance of our study. Our results show that the calibration has to be performed by selecting the tomographic samples in the simulations, consistent with the actual cosmic shear analysis, because the joint distributions of galaxy properties are found to vary with redshift. Ignoring this redshift variation could result in misestimating the shear bias by an amount that exceeds the allowed tolerance. To improve the calibration for future cosmic shear analyses, it will also be essential to correctly account for the measurement of photometric redshifts, which requires simulating multi-band observations.
We constrain the luminosity and redshift dependence of the intrinsic alignment (IA) of a nearly volume-limited sample of luminous red galaxies selected from the fourth public data release of the ...Kilo-Degree Survey (KiDS-1000). To measure the shapes of the galaxies, we used two complementary algorithms, finding consistent IA measurements for the overlapping galaxy sample. The global significance of IA detection across our two independent luminous red galaxy samples, with our favoured method of shape estimation, is ∼10.7
σ
. We find no significant dependence with redshift of the IA signal in the range 0.2 <
z
< 0.8, nor a dependence with luminosity below
L
r
≲ 2.9 × 10
10
h
−2
L
r
, ⊙
. Above this luminosity, however, we find that the IA signal increases as a power law, although our results are also compatible with linear growth within the current uncertainties. This behaviour motivates the use of a broken power law model when accounting for the luminosity dependence of IA contamination in cosmic shear studies.
Intrinsic galaxy alignments are a source of bias for weak lensing measurements as well as a tool for understanding galaxy formation and evolution. In this work, we measure the alignment of shapes of ...satellite galaxies, in galaxy groups, with respect to the brightest group galaxy (BGG), as well as alignments of the BGG shape with the satellite positions, using the highly complete Galaxy And Mass Assembly (GAMA) spectroscopic survey and deep imaging from the Kilo Degree Survey. We control systematic errors with dedicated image simulations and measure accurate shapes using the DEIMOS shape measurement method. We find a significant satellite radial alignment signal, which vanishes at large separations from the BGG. We do not identify any strong trends of the signal with galaxy absolute magnitude or group mass. The alignment signal is dominated by red satellites. We also find that the outer regions of galaxies are aligned more strongly than their inner regions, by varying the radial weight employed during the shape measurement process. This behaviour is evident for both red and blue satellites. BGGs are also found to be aligned with satellite positions, with this alignment being stronger when considering the innermost satellites, using red BGGs and the shape of the outer region of the BGG. Lastly, we measure the global intrinsic alignment signal in the GAMA sample for two different radial weight functions and find no significant difference.
ABSTRACT
We use KiDS weak lensing data to measure variations in mean halo mass as a function of several key galaxy properties (namely stellar colour, specific star formation rate, Sérsic index, and ...effective radius) for a volume-limited sample of GAMA galaxies in a narrow stellar mass range M* ∼ (2–5) × 1010 M⊙. This mass range is particularly interesting, inasmuch as it is where bimodalities in galaxy properties are most pronounced, and near to the break in both the galaxy stellar mass function and the stellar-to-halo mass relation (SHMR). In this narrow mass range, we find that both size and Sérsic index are better predictors of halo mass than either colour or SSFR, with the data showing a slight preference for Sérsic index. In other words, we find that mean halo mass is more tightly correlated with galaxy structure than either past star formation history or current star formation rate. Our results lead to an approximate lower bound on the dispersion in halo masses among log M* ≈ 10.5 galaxies: We find that the dispersion is ≳0.3 dex. This would imply either that offsets from the mean SHMR are closely coupled to size/structure or that the dispersion in the SHMR is larger than what past results have suggested. Our results thus provide new empirical constraints on the relationship between stellar and halo mass assembly at this particularly interesting mass range.
ABSTRACT
We present a comparison between approximated methods for the construction of mock catalogues based on the halo-bias mapping technique. To this end, we use as reference a high-resolution ...N-body simulation of 38403 dark matter particles on a 400 h−1 Mpc cube box from the Multidark suite. In particular, we explore parametric versus non-parametric bias mapping approaches and compare them at reproducing the halo distribution in terms of the two- and three-point statistics down to $\sim 10^8\, {\rm M}_{\odot }\, h^{-1}$ halo masses. Our findings demonstrate that the parametric approach remains inaccurate even including complex deterministic and stochastic components. On the contrary, the non-parametric one is indistinguishable from the reference N-body calculation in the power spectrum beyond $k=1\, h\, {\rm Mpc}^{-1}$, and in the bispectrum for typical configurations relevant to baryon acoustic oscillation analysis. We conclude that approaches which extract the full bias information from N-body simulations in a non-parametric fashion are ready for the analysis of the new generation of large-scale structure surveys.
Clustering of dark matter halos has been shown to depend on halo properties beyond mass, such as halo concentration, a phenomenon referred to as halo assembly bias. Standard halo-occupation models ...(HOD) in large-scale structure studies assume that halo mass alone is sufficient in characterizing the connection between galaxies and halos. Modeling of galaxy clustering can face systematic effects if the number of galaxies within a halo is correlated with other halo properties. Using the Small MultiDark-Planck high-resolution N-body simulation and the clustering measurements of the Sloan Digital Sky Survey DR7 main galaxy sample, we investigate the extent to which the concentration dependence of halo occupation can be constrained. Furthermore, we study how allowing for the concentration dependence can improve our modeling of galaxy clustering. Our constraints on HOD with assembly bias suggest that satellite population is not correlated with halo concentration at fixed halo mass. At fixed halo mass, our constraints favor lack of correlation between the occupation of centrals and halo concentration in the most luminous samples (Mr < −21.5, −21), and modest correlation in the Mr < −20.5, −20, −19.5 samples. We show that, in comparison with abundance-matching mock catalogs, our findings suggest qualitatively similar but modest levels of the impact of halo assembly bias on galaxy clustering. The effect is only present in the central occupation and becomes less significant in brighter galaxy samples. Furthermore, by performing model comparison based on information criteria, we find that in most cases, the standard mass-only HOD model is still favored by the observations.
Abstract
Standard approaches to Bayesian parameter inference in large-scale structure assume a Gaussian functional form (chi-squared form) for the likelihood. This assumption, in detail, cannot be ...correct. Likelihood free inferences such as approximate Bayesian computation (ABC) relax these restrictions and make inference possible without making any assumptions on the likelihood. Instead ABC relies on a forward generative model of the data and a metric for measuring the distance between the model and data. In this work, we demonstrate that ABC is feasible for LSS parameter inference by using it to constrain parameters of the halo occupation distribution (HOD) model for populating dark matter haloes with galaxies. Using specific implementation of ABC supplemented with population Monte Carlo importance sampling, a generative forward model using HOD and a distance metric based on galaxy number density, two-point correlation function and galaxy group multiplicity function, we constrain the HOD parameters of mock observation generated from selected ‘true' HOD parameters. The parameter constraints we obtain from ABC are consistent with the ‘true' HOD parameters, demonstrating that ABC can be reliably used for parameter inference in LSS. Furthermore, we compare our ABC constraints to constraints we obtain using a pseudo-likelihood function of Gaussian form with MCMC and find consistent HOD parameter constraints. Ultimately, our results suggest that ABC can and should be applied in parameter inference for LSS analyses.
This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the ...covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ice-cola, peak patch, and pinocchio), methods that require a calibration with N-body simulations (patchy and halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (lognormal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5 per cent for our lower mass threshold and 10 per cent for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20 per cent depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others.
ABSTRACT
We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences ...in covariances into parameter constraints related to growth rate of structure, Alcock–Paczynski distortions, and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE–COLA, pinocchio, and peakpatch), methods that rely on halo assignment schemes on to dark matter overdensities calibrated with a target N-body run (halogen, patchy), and two standard assumptions about the full density probability distribution function (Gaussian and lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5{{\ \rm per\ cent}}$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10{{\ \rm per\ cent}}$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10{{\ \rm per\ cent}}$ of those from the N-body-based covariance.
We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogues to ...the same quantities obtained from full N-body simulations. To this purpose we employ a large set of 300 realizations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and non-linear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10 per cent of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to 2-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation by comparing with a larger set of several thousands of realizations performed with one approximate method.