ABSTRACT
Upcoming weak lensing surveys require a detailed theoretical understanding of the matter power spectrum in order to derive accurate and precise cosmological parameter values. While galaxy ...formation is known to play an important role, its precise effects are currently unknown. We present a set of 92 matter power spectra from the OWLS, cosmo-OWLS, and BAryons and HAloes of MAssive Systems simulation suites, including different ΛCDM cosmologies, neutrino masses, subgrid prescriptions, and AGN feedback strengths. We conduct a detailed investigation of the dependence of the relative difference between the total matter power spectra in hydrodynamical and collisionless simulations on the effectiveness of stellar and AGN feedback, cosmology, and redshift. The strength of AGN feedback can greatly affect the power on a range of scales, while a lack of stellar feedback can greatly increase the effectiveness of AGN feedback on large scales. We also examine differences in the initial conditions of hydrodynamic and N-body simulations that can lead to an $\sim 1{{\ \rm per\ cent}}$ discrepancy in the large-scale power, and furthermore show our results to be insensitive to cosmic variance. We present an empirical model capable of predicting the effect of galaxy formation on the matter power spectrum at z = 0 to within $1{{\ \rm per\ cent}}$ for $k\lt 1\, h\, \mathrm{Mpc}^{-1}$, given only the mean baryon fraction in galaxy groups. Differences in group baryon fractions can also explain the quantitative disagreement between predictions from the literature. All total and dark matter only power spectra in this library will be made publicly available at powerlib.strw.leidenuniv.nl.
We report results for the alignments of galaxies in the EAGLE and cosmo-OWLS hydrodynamical cosmological simulations as a function of galaxy separation (−1 ≤ log10(r/ h
−1 Mpc) ≤ 2) and halo mass ...(10.7 ≤ log10(M
200/h
−1 M⊙) ≤ 15). We focus on two classes of alignments: the orientations of galaxies with respect to either the directions to, or the orientations of, surrounding galaxies. We find that the strength of the alignment is a strongly decreasing function of the distance between galaxies. For galaxies hosted by the most massive haloes in our simulations the alignment can remain significant up to ∼100 Mpc. Galaxies hosted by more massive haloes show stronger alignment. At a fixed halo mass, more aspherical or prolate galaxies exhibit stronger alignments. The spatial distribution of satellites is anisotropic and significantly aligned with the major axis of the main host halo. The major axes of satellite galaxies, when all stars are considered, are preferentially aligned towards the centre of the main host halo. The predicted projected direction–orientation alignment, ϵg+(r
p), is in broad agreement with recent observations. We find that the orientation–orientation alignment is weaker than the orientation–direction alignment on all scales. Overall, the strength of galaxy alignments depends strongly on the subset of stars that are used to measure the orientations of galaxies and it is always weaker than the alignment of dark matter haloes. Thus, alignment models that use halo orientation as a direct proxy for galaxy orientation overestimate the impact of intrinsic galaxy alignments.
Halo-based models have been successful in predicting the clustering of matter. However, the validity of the postulate that the clustering is fully determined by matter inside haloes remains largely ...untested, and it is not clear a priori whether non-virialized matter might contribute significantly to the non-linear clustering signal. Here, we investigate the contribution of haloes to the matter power spectrum as a function of both scale and halo mass by combining a set of cosmological N-body simulations to calculate the contributions of different spherical overdensity regions, Friends-of-Friends (FoF) groups and matter outside haloes to the power spectrum. We find that matter inside spherical overdensity regions of size R
200, mean cannot account for all power for 1 ≲ k ≲ 100 h Mpc− 1, regardless of the minimum halo mass. At most, it accounts for 95 per cent of the power (k ≳ 20 h Mpc− 1). For 2 ≲ k ≲ 10 h Mpc− 1, haloes with mass
$M_\mathrm{200,mean} \lesssim 10^{11}\,h^{-1}\,\mathrm{M}_{\odot }$
contribute negligibly to the power spectrum, and our results appear to be converged with decreasing halo mass. When haloes are taken to be regions of size R
200, crit, the amount of power unaccounted for is larger on all scales. Accounting also for matter inside FoF groups but outside R
200, mean increases the contribution of halo matter on most scales probed here by 5–15 per cent. Matter inside FoF groups with M
200, mean > 109 h
−1 M⊙ accounts for essentially all power for 3 < k < 100 h Mpc−1. We therefore expect halo models that ignore the contribution of matter outside R
200, mean to overestimate the contribution of haloes of any mass to the power on small scales (k ≳ 1 h Mpc− 1).
Upcoming weak lensing surveys, such as LSST, EUCLID and WFIRST, aim to measure the matter power spectrum with unprecedented accuracy. In order to fully exploit these observations, models are needed ...that, given a set of cosmological parameters, can predict the non-linear matter power spectrum at the level of 1 per cent or better for scales corresponding to comoving wavenumbers 0.1 ≲k≲ 10 h Mpc−1. We have employed the large suite of simulations from the OverWhelmingly Large Simulations (OWLS) project to investigate the effects of various baryonic processes on the matter power spectrum. In addition, we have examined the distribution of power over different mass components, the back-reaction of the baryons on the cold dark matter and the evolution of the dominant effects on the matter power spectrum. We find that single baryonic processes are capable of changing the power spectrum by up to several tens of per cent. Our simulation that includes AGN feedback, which we consider to be our most realistic simulation as, unlike those used in previous studies, it has been shown to solve the overcooling problem and to reproduce optical and X-ray observations of groups of galaxies, predicts a decrease in power relative to a dark matter only simulation ranging, at z= 0, from 1 per cent at k≈ 0.3 h Mpc−1 to 10 per cent at k≈ 1 h Mpc−1 and to 30 per cent at k≈ 10 h Mpc−1. This contradicts the naive view that baryons raise the power through cooling, which is the dominant effect only for k≳ 70 h Mpc−1. Therefore, baryons, and particularly AGN feedback, cannot be ignored in theoretical power spectra for k≳ 0.3 h Mpc−1. It will thus be necessary to improve our understanding of feedback processes in galaxy formation, or at least to constrain them through auxiliary observations, before we can fulfil the goals of upcoming weak lensing surveys.
We use matter power spectra from cosmological hydrodynamic simulations to quantify the effect of baryon physics on the weak gravitational lensing shear signal. The simulations consider a number of ...processes, such as radiative cooling, star formation, supernovae and feedback from active galactic nuclei (AGN). Van Daalen et al. used the same simulations to show that baryon physics, in particular the strong feedback that is required to solve the overcooling problem, modifies the matter power spectrum on scales relevant for cosmological weak lensing studies. As a result, the use of power spectra from dark matter simulations can lead to significant biases in the inferred cosmological parameters. We show that the typical biases are much larger than the precision with which future missions aim to constrain the dark energy equation of state, w
0. For instance, the simulation with AGN feedback, which reproduces X-ray and optical properties of groups of galaxies, gives rise to a ∼40 per cent bias in w
0. We also explore the effect of baryon physics on constraints on Ωm, σ8, the running of the spectral index, the mass of the neutrinos and models of warm dark matter. We demonstrate that the modification of the power spectrum is dominated by groups and clusters of galaxies, the effect of which can be modelled. We consider an approach based on the popular halo model and show that simple modifications can capture the main features of baryonic feedback. Despite its simplicity, we find that our model, when calibrated on the simulations, is able to reduce the bias in w
0 to a level comparable to the size of the statistical uncertainties for a Euclid-like mission. While observations of the gas and stellar fractions as a function of halo mass can be used to calibrate the model, hydrodynamic simulations will likely still be needed to extend the observed scaling relations down to halo masses of 1012
h
−1 M⊙.
We use cosmological hydrodynamical simulations to investigate how the inclusion of physical processes relevant to galaxy formation (star formation, metal-line cooling, stellar winds, supernovae and ...feedback from active galactic nuclei, AGN) change the properties of haloes, over four orders of magnitude in mass. We find that gas expulsion and the associated dark matter (DM) expansion induced by supernova-driven winds are important for haloes with masses M
200 ≲ 1013 M⊙, lowering their masses by up to 20 per cent relative to a DM-only model. AGN feedback, which is required to prevent overcooling, has a significant impact on halo masses all the way up to cluster scales (M
200 ∼ 1015 M⊙). Baryon physics changes the total mass profiles of haloes out to several times the virial radius, a modification that cannot be captured by a change in the halo concentration. The decrease in the total halo mass causes a decrease in the halo mass function of about 20 per cent. This effect can have important consequences for the abundance matching technique as well as for most semi-analytic models of galaxy formation. We provide analytic fitting formulae, derived from simulations that reproduce the observed baryon fractions, to correct halo masses and mass functions from DM-only simulations. The effect of baryon physics (AGN feedback in particular) on cluster number counts is about as large as changing the cosmology from Wilkinson Microwave Anisotropy Probe 7 to Planck, even when a moderately high-mass limit of M
500 ≈ 1014 M⊙ is adopted. Thus, for precision cosmology the effects of baryons must be accounted for.
ABSTRACT
The clustering of matter, as measured by the matter power spectrum, informs us about cosmology, dark matter, and baryonic effects on the distribution of matter in the universe. Using ...cosmological hydrodynamical simulations from the cosmo-OWLS and BAHAMAS simulation projects, we investigate the contribution of power in haloes with various masses, to the full power spectrum, as well as the power ratio between baryonic and dark matter only (DMO) simulations for a matched (between simulations) and an unmatched set of haloes. We find that the presence of AGN feedback suppresses the power on all scales for haloes of all masses examined (1011.25 ≤ M500, crit ≤ $10^{14.75}\, \mathrm{M_\odot }/h$), by ejecting matter from within $r_{500,\mathrm{c}}\,$ to $r_{200,\mathrm{m}}\,$ and potentially beyond in massive haloes (M500, crit ≳ $10^{{13}}\, \mathrm{M_\odot }/h$), and likely impeding the growth of lower-mass haloes as a consequence. A lower AGN feedback temperature changes the behaviour of high-mass haloes (M500, crit ≥ $10^{{13.25}}\, \mathrm{M_\odot }/h$), damping the effects of AGN feedback at small scales, $k\, {{\gtrsim }}\, {{4}}\, h\mathrm{\, Mpc^{-1}}$. For $k\, {{\lesssim }}\, {{3}}\, h\mathrm{\, Mpc^{-1}}$, group-sized haloes ($10^{{14\pm 0.25}}\, \mathrm{M_\odot }/h$) dominate the power spectrum, while on smaller scales the combined contributions of lower-mass haloes to the full power spectrum rise above that of the group-sized haloes. Finally, we present a model for the power suppression due to feedback, which combines observed mean halo baryon fractions with halo mass fractions and halo-matter cross-spectra extracted from DMO simulations to predict the power suppression to per cent level accuracy down to $k\, {{\approx }}\, {{10}}\, h\mathrm{\, Mpc^{-1}}$ without any free parameters.
The observed clustering of galaxies and the cross-correlation of galaxies and mass provide important constraints on both cosmology and models of galaxy formation. Even though the dissipation and ...feedback processes associated with galaxy formation are thought to affect the distribution of matter, essentially all models used to predict clustering data are based on collisionless simulations. Here, we use large hydrodynamical simulations to investigate how galaxy formation affects the autocorrelation functions of galaxies and subhaloes, as well as their cross-correlation with matter. We show that the changes due to the inclusion of baryons are not limited to small scales and are even present in samples selected by subhalo mass. Samples selected by subhalo mass cluster ∼10 per cent more strongly in a baryonic run on scales r≳1 h
−1Mpc, and this difference increases for smaller separations. While the inclusion of baryons boosts the clustering at fixed subhalo mass on all scales, the sign of the effect on the cross-correlation of subhaloes with matter can vary with radius. We show that the large-scale effects are due to the change in subhalo mass caused by the strong feedback associated with galaxy formation and may therefore not affect samples selected by number density. However, on scales r≲r
vir significant differences remain after accounting for the change in subhalo mass. We conclude that predictions for galaxy–galaxy and galaxy–mass clustering from models based on collisionless simulations will have errors greater than 10 per cent on sub-Mpc scales, unless the simulation results are modified to correctly account for the effects of baryons on the distributions of mass and satellites.
ABSTRACT
To fully take advantage of the data provided by large-scale structure surveys, we need to quantify the potential impact of baryonic effects, such as feedback from active galactic nuclei ...(AGN) and star formation, on cosmological observables. In simulations, feedback processes originate on scales that remain unresolved. Therefore, they need to be sourced via subgrid models that contain free parameters. We use machine learning to calibrate the AGN and stellar feedback models for the FLAMINGO (Fullhydro Large-scale structure simulations with All-sky Mapping for the Interpretation of Next Generation Observations) cosmological hydrodynamical simulations. Using Gaussian process emulators trained on Latin hypercubes of 32 smaller volume simulations, we model how the galaxy stellar mass function (SMF) and cluster gas fractions change as a function of the subgrid parameters. The emulators are then fit to observational data, allowing for the inclusion of potential observational biases. We apply our method to the three different FLAMINGO resolutions, spanning a factor of 64 in particle mass, recovering the observed relations within the respective resolved mass ranges. We also use the emulators, which link changes in subgrid parameters to changes in observables, to find models that skirt or exceed the observationally allowed range for cluster gas fractions and the SMF. Our method enables us to define model variations in terms of the data that they are calibrated to rather than the values of specific subgrid parameters. This approach is useful, because subgrid parameters are typically not directly linked to particular observables, and predictions for a specific observable are influenced by multiple subgrid parameters.