ABSTRACT
Luminous matter produces very energetic events, such as active galactic nuclei and supernova explosions, that significantly affect the internal regions of galaxy clusters. Although the ...current uncertainty in the effect of baryonic physics on cluster statistics is subdominant as compared to other systematics, the picture is likely to change soon as the amount of high-quality data is growing fast, urging the community to keep theoretical systematic uncertainties below the ever-growing statistical precision. In this paper, we study the effect of baryons on galaxy clusters, and their impact on the cosmological applications of clusters, using the magneticum suite of cosmological hydrodynamical simulations. We show that the impact of baryons on the halo mass function can be recast in terms on a variation of the mass of the haloes simulated with pure N-body, when baryonic effects are included. The halo mass function and halo bias are only indirectly affected. Finally, we demonstrate that neglecting baryonic effects on haloes mass function and bias would significantly alter the inference of cosmological parameters from high-sensitivity next-generations surveys of galaxy clusters.
Abstract
Using stellar population synthesis models to infer star formation histories (SFHs), we analyze photometry and spectroscopy of a large sample of quiescent galaxies that are members of ...Sunyaev–Zel’dovich (SZ)-selected galaxy clusters across a wide range of redshifts. We calculate stellar masses and mass-weighted ages for 837 quiescent cluster members at 0.3 <
z
< 1.4 using rest-frame optical spectra and the Python-based
Prospector
framework, from 61 clusters in the SPT-GMOS Spectroscopic Survey (0.3 <
z
< 0.9) and three clusters in the SPT Hi-z cluster sample (1.25 <
z
< 1.4). We analyze spectra of subpopulations divided into bins of redshift, stellar mass, cluster mass, and velocity-radius phase-space location, as well as by creating composite spectra of quiescent member galaxies. We find that quiescent galaxies in our data set sample a diversity of SFHs, with a median formation redshift (corresponding to the lookback time from the redshift of observation to when a galaxy forms 50% of its mass,
t
50
) of
z
= 2.8 ± 0.5, which is similar to or marginally higher than that of massive quiescent field and cluster galaxy studies. We also report median age–stellar mass relations for the full sample (age of the universe at
t
50
(Gyr) = 2.52 (±0.04)–1.66 (±0.12) log
10
(
M
/10
11
M
⊙
)) and recover downsizing trends across stellar mass; we find that massive galaxies in our cluster sample form on aggregate ∼0.75 Gyr earlier than lower-mass galaxies. We also find marginally steeper age–mass relations at high redshifts, and report a bigger difference in formation redshifts across stellar mass for fixed environment, relative to formation redshifts across environment for fixed stellar mass.
We have simulated the formation of a galaxy cluster in a Λ cold dark matter universe using 13 different codes modelling only gravity and non-radiative hydrodynamics (ramses, ART, arepo, hydra and ...nine incarnations of gadget). This range of codes includes particle-based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various gadget implementations span classic and modern smoothed particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at z = 0, global properties such as mass and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ramses, art and arepo form extended entropy cores in the gas with rising central gas temperatures. Those codes employing classic SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid-based methods.
Abstract
The environments where galaxies reside crucially shape their star formation histories. We investigate a large sample of 1626 cluster galaxies located within 105 galaxy clusters spanning a ...large range in redshift (0.26 <
z
< 1.13). The galaxy clusters are massive (
M
500
≳ 2 × 10
14
M
⊙
) and uniformly selected from the SPT and ACT Sunyaev–Zel’dovich surveys. With spectra in hand for thousands of cluster members, we use the galaxies’ position in projected phase space as a proxy for their infall times, which provides a more robust measurement of environment than quantities such as projected clustercentric radius. We find clear evidence for a gradual age increase of the galaxy’s mean stellar populations (∼0.71 ± 0.4 Gyr based on a 4000 Å break, D
n
4000) with the time spent in the cluster environment. This environmental quenching effect is found regardless of galaxy luminosity (faint or bright) and redshift (low or high-
z
), although the exact stellar age of galaxies depends on both parameters at fixed environmental effects. Such a systematic increase of D
n
4000 with infall proxy would suggest that galaxies that were accreted into hosts earlier were quenched earlier due to longer exposure to environmental effects such as ram pressure stripping and starvation. Compared to the typical dynamical timescales of 1–3 Gyr of cluster galaxies, the relatively small age increase (∼0.71 ± 0.4 Gyr) found in our sample galaxies seems to suggest that a slow environmental process such as starvation is the dominant quenching pathway. Our results provide new insights into environmental quenching effects spanning a large range in cosmic time (∼5.2 Gyr,
z
= 0.26–1.13) and demonstrate the power of using a kinematically derived infall time proxy.
Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 ...state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern smoothed particle hydrodynamics (SPH) approaches. For each code represented we have a dark-matter-only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency ⪅20 per cent) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R
200. However, we see larger differences for quantities within R
2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include active galactic nucleus feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.
ABSTRACT
We employ a set of Magneticum cosmological hydrodynamic simulations that span over 15 different cosmologies, and extract masses and concentrations of all well-resolved haloes between z = ...0 and 1 for critical overdensities $\Delta _\textrm {vir}, \Delta _{200c}, \Delta _{500c}, \Delta _{2500c}$ and mean overdensity Δ200m. We provide the first mass–concentration (Mc) relation and sparsity relation (i.e. MΔ1 − MΔ2 mass conversion) of hydrodynamic simulations that is modelled by mass, redshift, and cosmological parameters Ωm, Ωb, σ8, h0 as a tool for observational studies. We also quantify the impact that the Mc relation scatter and the assumption of Navarro–Frank–White (NFW) density profiles have on the uncertainty of the sparsity relation. We find that converting masses with the aid of an Mc relation carries an additional fractional scatter ($\approx 4{{\ \rm per\ cent}}$) originated from deviations from the assumed NFW density profile. For this reason, we provide a direct mass–mass conversion relation fit that depends on redshift and cosmological parameters. We release the package hydro_mc, a python tool that perform all kind of conversions presented in this paper.
In this study, we present a detailed, statistical analysis of black hole growth and the evolution of active galactic nuclei (AGN) using cosmological hydrodynamic simulations run down to z = 0. The ...simulations self-consistently follow radiative cooling, star formation, metal enrichment, black hole growth and associated feedback processes from both Type II/Ia supernovae and AGN. We consider two simulation runs, one with a large comoving volume of (500 Mpc)3 and one with a smaller volume of (68 Mpc)3 but with a factor of almost 20 higher mass resolution. We compare the predicted statistical properties of AGN with results from large observational surveys. Consistently with previous results, our simulations can widely match observed black hole properties of the local Universe. Furthermore, our simulations can successfully reproduce the evolution of the bolometric AGN luminosity function for both the low-luminosity and the high-luminosity end up to z = 3.0, only at z = 1.5–2.5, the low-luminosity end is overestimated by up to 1 dex. In addition, the smaller but higher resolution run is able to match the observational data of the low bolometric luminosity end at higher redshifts z = 3–4. We also perform a direct comparison with the observed soft and hard X-ray luminosity functions of AGN, including an empirical correction for a torus-level obscuration, and find a similarly good agreement. These results nicely demonstrate that the observed ‘antihierarchical’ trend in the AGN number density evolution (i.e. the number densities of luminous AGN peak at higher redshifts than those of faint AGN) is self-consistently predicted by our simulations. Implications of this downsizing behaviour on active black holes, their masses and Eddington ratios are discussed. Overall, the downsizing behaviour in the AGN number density as a function of redshift can be mainly attributed to the evolution of the gas density in the resolved vicinity of a (massive) black hole (which is depleted with evolving time as a consequence of star formation and AGN feedback).
ABSTRACT
We study the potential of the kinematic Sunyaev–Zel’dovich (kSZ) effect as a probe for cosmology, focusing on the pairwise method. The main challenge is disentangling the cosmologically ...interesting mean pairwise velocity from the cluster optical depth and the associated uncertainties on the baryonic physics in clusters. Furthermore, the pairwise kSZ signal might be affected by internal cluster motions or correlations between velocity and optical depth. We investigate these effects using the Magneticum cosmological hydrodynamical simulations, one of the largest simulations of this kind performed to date. We produce thermal SZ and kSZ maps with an area of ≃ 1600 deg2, and the corresponding cluster catalogues with M500c ≳ 3 × 1013 h−1 M⊙ and z ≲ 2. From these data sets, we calibrate a scaling relation between the average Compton-y parameter and optical depth. We show that this relation can be used to recover an accurate estimate of the mean pairwise velocity from the kSZ effect, and that this effect can be used as an important probe of cosmology. We discuss the impact of theoretical and observational systematic effects, and find that further work on feedback models is required to interpret future high-precision measurements of the kSZ effect.
ABSTRACT
We present an analysis aimed at combining cosmological constraints from number counts of galaxy clusters identified through the Sunyaev–Zeldovich effect, obtained with the South Pole ...Telescope (SPT), and from Lyman α spectra obtained with the MIKE/HIRES and X-shooter spectrographs. The SPT cluster analysis relies on mass calibration based on weak lensing measurements, while the Lyman α analysis is built over mock spectra extracted from hydrodynamical simulations. The resulting constraints exhibit a tension (∼3.3σ) between the low σ8 values preferred by the low-redshift cluster data, $\sigma _8=0.74 ^{+0.03}_{-0.04}$, and the higher one preferred by the high-redshift Lyman α data, $\sigma _8=0.91 ^{+0.03}_{-0.03}$. We present a detailed analysis to understand the origin of this tension and to establish whether it arises from systematic uncertainties related to the assumptions underlying the analyses of cluster counts and/or Lyman α forest. We found this tension to be robust with respect to the choice of modelling of the IGM, even when including possible systematics from unaccounted sub-Damped Lyman α (DLA) and Lyman-limit systems (LLS) in the Lyman α data. We conclude that to solve this tension would require a large bias on the cluster mass estimate, or large unaccounted errors on the Lyman α mean fluxes. Our results have important implications for future analyses based on cluster number counts from future large photometric surveys (e.g. Euclid and LSST) and on larger samples of high-redshift quasar spectra (e.g. DESI and WEAVE surveys). If confirmed at the much higher statistical significance reachable by such surveys, this tension could represent a significant challenge for the standard ΛCDM paradigm.
We employ a set of Magneticum cosmological hydrodynamic simulations that span over \(15\) different cosmologies, and extract masses and concentrations of all well-resolved haloes between \(z=0-1\) ...for critical over-densities \(\Delta_\texttt{vir}, \Delta_{200c}, \Delta_{500c}, \Delta_{2500c}\) and mean overdensity \(\Delta_{200m}.\) We provide the first mass-concentration (Mc) relation and sparsity relation (i.e. \(M_{\Delta1} - M_{\Delta2}\) mass conversion) of hydrodynamic simulations that is modelled by mass, redshift and cosmological parameters \(\Omega_m, \Omega_b, \sigma_8, h_0\) as a tool for observational studies. We also quantify the impact that the Mc relation scatter and the assumption of NFW density profiles have on the uncertainty of the sparsity relation. We find that converting masses with the aid of a Mc relation carries an additional fractional scatter (\(\approx 4\%\)) originated from deviations from the assumed NFW density profile. For this reason we provide a direct mass-mass conversion relation fit that depends on redshift and cosmological parameters. We release the package hydro\_mc, a python tool that perform all kind of conversions presented in this paper.