We present an optimized variant of the halo model, designed to produce accurate matter power spectra well into the non-linear regime for a wide range of cosmological models. To do this, we introduce ...physically motivated free parameters into the halo-model formalism and fit these to data from high-resolution N-body simulations. For a variety of Λ cold dark matter (ΛCDM) and wCDM models, the halo-model power is accurate to ≃ 5 per cent for k ≤ 10h Mpc−1 and z ≤ 2. An advantage of our new halo model is that it can be adapted to account for the effects of baryonic feedback on the power spectrum. We demonstrate this by fitting the halo model to power spectra from the OWLS (OverWhelmingly Large Simulations) hydrodynamical simulation suite via parameters that govern halo internal structure. We are able to fit all feedback models investigated at the 5 per cent level using only two free parameters, and we place limits on the range of these halo parameters for feedback models investigated by the OWLS simulations. Accurate predictions to high k are vital for weak-lensing surveys, and these halo parameters could be considered nuisance parameters to marginalize over in future analyses to mitigate uncertainty regarding the details of feedback. Finally, we investigate how lensing observables predicted by our model compare to those from simulations and from halofit for a range of k-cuts and feedback models and quantify the angular scales at which these effects become important. Code to calculate power spectra from the model presented in this paper can be found at https://github.com/alexander-mead/hmcode.
We present an accurate non-linear matter power spectrum prediction scheme for a variety of extensions to the standard cosmological paradigm, which uses the tuned halo model previously developed in ...Mead et al. We consider dark energy models that are both minimally and non-minimally coupled, massive neutrinos and modified gravitational forces with chameleon and Vainshtein screening mechanisms. In all cases, we compare halo-model power spectra to measurements from high-resolution simulations. We show that the tuned halo-model method can predict the non-linear matter power spectrum measured from simulations of parametrized w(a) dark energy models at the few per cent level for k < 10 h Mpc−1, and we present theoretically motivated extensions to cover non-minimally coupled scalar fields, massive neutrinos and Vainshtein screened modified gravity models that result in few per cent accurate power spectra for k < 10 h Mpc−1. For chameleon screened models, we achieve only 10 per cent accuracy for the same range of scales. Finally, we use our halo model to investigate degeneracies between different extensions to the standard cosmological model, finding that the impact of baryonic feedback on the non-linear matter power spectrum can be considered independently of modified gravity or massive neutrino extensions. In contrast, considering the impact of modified gravity and massive neutrinos independently results in biased estimates of power at the level of 5 per cent at scales k > 0.5 h Mpc−1. An updated version of our publicly available hmcode can be found at https://github.com/alexander-mead/hmcode.
I examine differences in non-linear structure formation between cosmological models that share a z = 0 linear power spectrum in both shape and amplitude, but that differ via their growth history. ...N-body simulations of these models display an approximately identical large-scale-structure skeleton, but reveal deeply non-linear differences in the demographics and properties of haloes. I investigate to what extent the spherical-collapse model can help in understanding these differences, in both real and redshift space. I discuss how this is difficult to do if one attempts to identify haloes directly, because in that case one is subject to the vagaries of halo-finding algorithms. However, I demonstrate that the halo model of structure formation provides an accurate non-linear response in the power spectrum, but only if results from spherical collapse that include formation hysteresis are properly incorporated. I comment on how this fact can be used to provide per cent level accurate matter power-spectrum predictions for dark energy models for ... by using the halo model as a correction to accurate ...CDM simulations. In the Appendix, I provide some fitting functions for the linear-collapse threshold (...) and virialized overdensity (...) that are valid for a wide range of dark energy models. I also make my spherical-collapse code available at https://github.com/alexander-mead/collapse. (ProQuest: ... denotes formulae/symbols omitted.)
ABSTRACT
We derive a simple prescription for including beyond-linear halo bias within the standard, analytical halo-model power spectrum calculation. This results in a corrective term that is added ...to the usual two-halo term. We measure this correction using data from N-body simulations and demonstrate that it can boost power in the two-halo term by a factor of ∼2 at scales $k\sim 0.7\, h\mathrm{Mpc}^{-1}$, with the exact magnitude of the boost determined by the specific pair of fields in the two-point function. How this translates to the full power spectrum depends on the relative strength of the one-halo term, which can mask the importance of this correction to a greater or lesser degree, again depending on the fields. Generally, we find that our correction is more important for signals that arise from lower mass haloes. When comparing our calculation to simulated data, we find that the underprediction of power in the transition region between the two- and one-halo terms, which typically plagues halo-model calculations, is almost completely eliminated when including the full non-linear halo bias. We show improved results for the autospectra and cross-spectra of galaxies, haloes, and matter. In the specific case of matter–matter or matter–halo power, we note that a large fraction of the improvement comes from the non-linear biasing between low- and high-mass haloes. We envisage our model being useful in the analytical modelling of cross-correlation signals. Our non-linear bias halo-model code is available at https://github.com/alexander-mead/BNL.
ABSTRACT
We present an updated version of the hmcode augmented halo model that can be used to make accurate predictions of the non-linear matter power spectrum over a wide range of cosmologies. Major ...improvements include modelling of baryon-acoustic oscillation (BAO) damping in the power spectrum and an updated treatment of massive neutrinos. We fit our model to simulated power spectra and show that we can match the results with an root mean square (RMS) error of 2.5 per cent across a range of cosmologies, scales $k \lt 10\, h\, \mathrm{Mpc}^{-1}$, and redshifts z < 2. The error rarely exceeds 5 per cent and never exceeds 16 per cent. The worst-case errors occur at z ≃ 2, or for cosmologies with unusual dark energy equations of state. This represents a significant improvement over previous versions of hmcode, and over other popular fitting functions, particularly for massive-neutrino cosmologies with high neutrino mass. We also present a simple halo model that can be used to model the impact of baryonic feedback on the power spectrum. This six-parameter physical model includes gas expulsion by active galactic nuclei (AGN) feedback and encapsulates star formation. By comparing this model to data from hydrodynamical simulations, we demonstrate that the power spectrum response to feedback is matched at the <1 per cent level for z < 1 and $k\lt 20\, h\, \mathrm{Mpc}^{-1}$. We also present a single-parameter variant of this model, parametrized in terms of feedback strength, which is only slightly less accurate. We make code available for our non-linear and baryon models at https://github.com/alexander-mead/HMcode and it is also available within camb and soon within class.
We present cosmological parameter constraints from a tomographic weak gravitational lensing analysis of ~450 deg super( 2) of imaging data from the Kilo Degree Survey (KiDS). For a flat ... cold dark ...matter (...CDM) cosmology with a prior on H sub( 0) that encompasses the most recent direct measurements, we find S sub( 8) ... = 0.745 plus or minus 0.039. This result is in good agreement with other low-redshift probes of large-scale structure, including recent cosmic shear results, along with pre-Planck cosmic microwave background constraints. A 2.3... tension in S sub( 8) and 'substantial discordance' in the full parameter space is found with respect to the Planck 2015 results. We use shear measurements for nearly 15 million galaxies, determined with a new improved 'self-calibrating' version of lensfit validated using an extensive suite of image simulations. Four-band ugri photometric redshifts are calibrated directly with deep spectroscopic surveys. The redshift calibration is confirmed using two independent techniques based on angular cross-correlations and the properties of the photometric redshift probability distributions. Our covariance matrix is determined using an analytical approach, verified numerically with large mock galaxy catalogues. We account for uncertainties in the modelling of intrinsic galaxy alignments and the impact of baryon feedback on the shape of the non-linear matter power spectrum, in addition to the small residual uncertainties in the shear and redshift calibration. The cosmology analysis was performed blind. Our high-level data products, including shear correlation functions, covariance matrices, redshift distributions, and Monte Carlo Markov chains are available at http://kids.strw.leidenuniv.nl. (ProQuest: ... denotes formulae/symbols omitted.)
On the scale of galactic haloes, the distribution of matter in the cosmos is affected by energetic, non-gravitational processes, the so-called baryonic feedback. A lack of knowledge about the details ...of how feedback processes redistribute matter is a source of uncertainty for weak-lensing surveys, which accurately probe the clustering of matter in the Universe over a wide range of scales. We developed a cosmology-dependent model for the matter distribution that simultaneously accounts for the clustering of dark matter, gas, and stars. We informed our model by comparing it to power spectra measured from the
BAHAMAS
suite of hydrodynamical simulations. In addition to considering matter power spectra, we also considered spectra involving the electron-pressure field, which directly relates to the thermal Sunyaev-Zel’dovich (tSZ) effect. We fitted parameters in our model so that it can simultaneously model both matter and pressure data and such that the distribution of gas as inferred from tSZ has an influence on the matter spectrum predicted by our model. We present two variants, one that matches the feedback-induced suppression seen in the matter–matter power spectrum at the percent level and a second that matches the matter–matter data to a slightly lesser degree (≃2%). However, the latter is able to simultaneously model the matter–electron pressure spectrum at the ≃15% level. We envisage our models being used to simultaneously learn about cosmological parameters and the strength of baryonic feedback using a combination of tSZ and lensing auto- and cross-correlation data.
ABSTRACT
We present a general method to compute the non-linear matter power spectrum for dark energy (DE) and modified gravity scenarios with per cent-level accuracy. By adopting the halo model and ...non-linear perturbation theory, we predict the reaction of a lambda cold dark matter (ΛCDM) matter power spectrum to the physics of an extended cosmological parameter space. By comparing our predictions to N-body simulations we demonstrate that with no-free parameters we can recover the non-linear matter power spectrum for a wide range of different w0–wa DE models to better than 1 per cent accuracy out to k ≈ 1 $h \,{\rm Mpc}^{-1}$. We obtain a similar performance for both DGP and f(R) gravity, with the non-linear matter power spectrum predicted to better than 3 per cent accuracy over the same range of scales. When including direct measurements of the halo mass function from the simulations, this accuracy improves to 1 per cent. With a single suite of standard ΛCDM N-body simulations, our methodology provides a direct route to constrain a wide range of non-standard extensions to the concordance cosmology in the high signal-to-noise non-linear regime.
We present a tomographic cosmic shear analysis of the Kilo-Degree Survey (KiDS) combined with the VISTA Kilo-Degree Infrared Galaxy Survey. This is the first time that a full optical to near-infrared ...data set has been used for a wide-field cosmological weak lensing experiment. This unprecedented data, spanning 450 deg 2 , allows us to significantly improve the estimation of photometric redshifts, such that we are able to include robustly higher-redshift sources for the lensing measurement, and – most importantly – to solidify our knowledge of the redshift distributions of the sources. Based on a flat ΛCDM model we find S 8 ≡ σ 8 Ω m /0.3 = 0.737 +0.040 −0.036 in a blind analysis from cosmic shear alone. The tension between KiDS cosmic shear and the Planck-Legacy CMB measurements remains in this systematically more robust analysis, with S 8 differing by 2.3 σ . This result is insensitive to changes in the priors on nuisance parameters for intrinsic alignment, baryon feedback, and neutrino mass. KiDS shear measurements are calibrated with a new, more realistic set of image simulations and no significant B-modes are detected in the survey, indicating that systematic errors are under control. When calibrating our redshift distributions by assuming the 30-band COSMOS-2015 photometric redshifts are correct (following the Dark Energy Survey and the Hyper Suprime-Cam Survey), we find the tension with Planck is alleviated. The robust determination of source redshift distributions remains one of the most challenging aspects for future cosmic shear surveys.
We present and test a method for modifying the catalogue of dark matter haloes produced from a given cosmological simulation, so that it resembles the result of a simulation with an entirely ...different set of parameters. This extends the method of Angulo & White, which rescales the full particle distribution from a simulation. Working directly with the halo catalogue offers an advantage in speed, and also allows modifications of the internal structure of the haloes to account for non-linear differences between cosmologies. Our method can be used directly on a halo catalogue in a self-contained manner without any additional information about the overall density field; although the large-scale displacement field is required by the method, this can be inferred from the halo catalogue alone. We show proof of concept of our method by rescaling a matter-only simulation with no baryon acoustic oscillation (BAO) features to a more standard Λ cold dark matter model containing a cosmological constant and a BAO signal. In conjunction with the halo occupation approach, this method provides a basis for the rapid generation of mock galaxy samples spanning a wide range of cosmological parameters.