Stellar Property Statistics of Massive Halos from Cosmological Hydrodynamics Simulations: Common Kernel Shapes

We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generates samples of $1000$ or more halos with $M_{\rm halo}> 10^{13.5} M_{\odot}$ at redshift $z=0$. Applying localized, linear regression (LLR), we extract halo mass-conditioned statistics (normalizations, slopes, and intrinsic covariance) for a three-element stellar property vector consisting of: i) $N_{sat}$, the number of satellite galaxies with stellar mass, $M_{\star, \rm sat} > 10^{10} M_{\odot}$ within radius $R_{200c}$ of the halo; ii) $M_{\star,\rm tot}$, the total stellar mass within that radius, and; iii) $M_{\star,\rm BCG}$, the gravitationally-bound stellar mass of the central galaxy within a $100 \, \rm kpc$ radius. Scaling parameters for the three properties with halo mass show mild differences among the simulations, in part due to numerical resolution, but there is qualitative agreement on property correlations, with halos having smaller than average central galaxies tending to also have smaller total stellar mass and a larger number of satellite galaxies. Marginalizing over total halo mass, we find the satellite galaxy kernel, $p(\ln N_{sat}\,|\,M_{\rm halo},z)$ to be consistently skewed left, with skewness parameter $\gamma = -0.91 \pm 0.02$, while that of $\ln M_{\star,\rm tot}$ is closer to log-normal, in all three simulations. The highest resolution simulations find $\gamma \simeq -0.8$ for the $z=0$ shape of $p(\ln M_{\star,\rm BCG}\,|\,M_{\rm halo},z)$ and also that the fractional scatter in total stellar mass is below $10\%$ in halos more massive than $10^{14.3} M_{\odot}$

Galaxy velocity bias in cosmological simulations: towards per cent-level calibration

Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, bv – defined as the ratio of galaxy and DM velocity dispersions – as a function of redshift, host halo mass, and galaxy stellar mass threshold (M ,sat), for massive haloes (M200c > 1013.5 M ) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of bv increasing with M200c and decreasing with redshift and M ,sat. The ensemble-estimated theoretical uncertainty in bv is 2–3 per cent, but becomes percent-level when considering only the three highest resolution simulations. We update the mass–richness normalization for an SDSS redMaPPer cluster sample, and find our improved bv estimates reduce the normalization uncertainty from 22 to 8 per cent, demonstrating that dynamical mass estimation is competitive with weak lensing mass estimation. We discuss necessary steps for further improving this precision. Our estimates for bv (M200c, M ,sat, z) are made publicly available

The Three Hundred Project: the gizmo-simba run

We introduce gizmo-simba, a new suite of galaxy cluster simulations within The Three Hundred project. The Three Hundred consists of zoom re-simulations of 324 clusters with M 200≳ 1014.8, M ⊙ drawn from the MultiDark-Planck N-body simulation, run using several hydrodynamic and semi-analytical codes. The gizmo-simba suite adds a state-of-the-art galaxy formation model based on the highly successful Simba simulation, mildly re-calibrated to match z = 0 cluster stellar properties. Comparing to The Three Hundred zooms run with gadget-x, we find intrinsic differences in the evolution of the stellar and gas mass fractions, BCG ages, and galaxy colour-magnitude diagrams, with gizmo-simba generally providing a good match to available data at z ≈ 0. gizmo-simba's unique black hole growth and feedback model yields agreement with the observed BH scaling relations at the intermediate-mass range and predicts a slightly different slope at high masses where few observations currently lie. Gizmo-Simba provides a new and novel platform to elucidate the co-evolution of galaxies, gas, and black holes within the densest cosmic environments

Shocks in the stacked Sunyaev-Zel'dovich profiles of clusters II: Measurements from SPT-SZ + Planck Compton-y map

We search for the signature of cosmological shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the latest Compton-y maps from the 2500 deg2 SPT-SZ survey on the locations of clusters identified in that same data set. The sample contains 516 clusters with mean mass 〈 M 200m〉 = 1014.9\, M⊙ and redshift z= 0.55. We analyse in parallel a set of zoom-in hydrodynamical simulations from the three hundred project. The SPT-SZ data show two features: (i) a pressure deficit at R/R200m = 1.08 ± 0.09, measured at 3.1σ significance and not observed in the simulations, and; (ii) a sharp decrease in pressure at R/R200m = 4.58 ± 1.24 at 2.0σ significance. The pressure deficit is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions, and the second feature is consistent with accretion shocks seen in previous studies. We split the cluster sample by redshift and mass, and find both features exist in all cases. There are also no significant differences in features along and across the cluster major axis, whose orientation roughly points towards filamentary structure. As a consistency test, we also analyse clusters from the Planck and Atacama Cosmology Telescope Polarimeter surveys and find quantitatively similar features in the pressure profiles. Finally, we compare the accretion shock radius (R sh,\, acc) with existing measurements of the splashback radius (Rsp) for SPT-SZ and constrain the lower limit of the ratio, R sh,\, acc/R sp\> 2.16\± 0.59

Shocks in the Stacked Sunyaev-Zel'dovich Profiles of Clusters II: Measurements from SPT-SZ + Planck Compton-y Map

We search for the signature of cosmological shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the latest Compton-y maps from the 2500 deg^2 SPT-SZ survey on the locations of clusters identified in that same dataset. The sample contains 516 clusters with mean mass = 1e14.9 msol and redshift = 0.55. We analyze in parallel a set of zoom-in hydrodynamical simulations from The Three Hundred project. The SPT-SZ data show two features: (i) a pressure deficit at R/R200m = $1.08 \pm 0.09$, measured at $3.1\sigma$ significance and not observed in the simulations, and; (ii) a sharp decrease in pressure at R/R200m = $4.58 \pm 1.24$ at $2.0\sigma$ significance. The pressure deficit is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions, and the second feature is consistent with accretion shocks seen in previous studies. We split the cluster sample by redshift and mass, and find both features exist in all cases. There are also no significant differences in features along and across the cluster major axis, whose orientation roughly points towards filamentary structure. As a consistency test, we also analyze clusters from the Planck and Atacama Cosmology Telescope Polarimeter surveys and find quantitatively similar features in the pressure profiles. Finally, we compare the accretion shock radius (Rsh_acc) with existing measurements of the splashback radius (Rsp) for SPT-SZ and constrain the lower limit of the ratio, Rsh_acc/Rsp > $2.16 \pm 0.59$.Comment: [v1]: 8 Figures, 16 Pages in Main text. [v2]: Added text to discussion. Version accepted in MNRA

Beyond the 3rd moment: A practical study of using lensing convergence CDFs for cosmology with DES Y3

Widefield surveys of the sky probe many clustered scalar fields -- such as galaxy counts, lensing potential, gas pressure, etc. -- that are sensitive to different cosmological and astrophysical processes. Our ability to constrain such processes from these fields depends crucially on the statistics chosen to summarize the field. In this work, we explore the cumulative distribution function (CDF) at multiple scales as a summary of the galaxy lensing convergence field. Using a suite of N-body lightcone simulations, we show the CDFs' constraining power is modestly better than that of the 2nd and 3rd moments of the field, as they approximately capture the information from all moments of the field in a concise data vector. We then study the practical aspects of applying the CDFs to observational data, using the first three years of the Dark Energy Survey (DES Y3) data as an example, and compute the impact of different systematics on the CDFs. The contributions from the point spread function are 2-3 orders of magnitude below the cosmological signal, while those from reduced shear approximation contribute $\lesssim 1\%$ to the signal. Source clustering effects and baryon imprints contribute $1-10\%$. Enforcing scale cuts to limit systematics-driven biases in parameter constraints degrades these constraints a noticeable amount, and this degradation is similar for the CDFs and the moments. We also detect correlations between the observed convergence field and the shape noise field at $13\sigma$. We find that the non-Gaussian correlations in the noise field must be modeled accurately to use the CDFs, or other statistics sensitive to all moments, as a rigorous cosmology tool.Comment: 21 pages, 12 figure

The scatter in the galaxy-halo connection: a machine learning analysis

We apply machine learning (ML), a powerful method for uncovering complex correlations in high-dimensional data, to the galaxy-halo connection of cosmological hydrodynamical simulations. The mapping between galaxy and halo variables is stochastic in the absence of perfect information, but conventional ML models are deterministic and hence cannot capture its intrinsic scatter. To overcome this limitation, we design an ensemble of neural networks with a Gaussian loss function that predict probability distributions, allowing us to model statistical uncertainties in the galaxy-halo connection as well as its best-fitting trends. We extract a number of galaxy and halo variables from the Horizon-AGN and IllustrisTNG100-1 simulations and quantify the extent to which knowledge of some subset of one enables prediction of the other. This allows us to identify the key features of the galaxy-halo connection and investigate the origin of its scatter in various projections. We find that while halo properties beyond mass account for up to 50 per cent of the scatter in the halo-To-stellar mass relation, the prediction of stellar half-mass radius or total gas mass is not substantially improved by adding further halo properties. We also use these results to investigate semi-Analytic models for galaxy size in the two simulations, finding that assumptions relating galaxy size to halo size or spin are not successful

A test of the standard cosmological model with geometry and growth

We perform a general test of the $\Lambda{\rm CDM}$ and $w {\rm CDM}$ cosmological models by comparing constraints on the geometry of the expansion history to those on the growth of structure. Specifically, we split the total matter energy density, $\Omega_M$, and (for $w {\rm CDM}$) dark energy equation of state, $w$, into two parameters each: one that captures the geometry, and another that captures the growth. We constrain our split models using current cosmological data, including type Ia supernovae, baryon acoustic oscillations, redshift space distortions, gravitational lensing, and cosmic microwave background (CMB) anisotropies. We focus on two tasks: (i) constraining deviations from the standard model, captured by the parameters $\Delta\Omega_M \equiv \Omega_M^{\rm grow}-\Omega_M^{\rm geom}$ and $\Delta w \equiv w^{\rm grow}-w^{\rm geom}$, and (ii) investigating whether the $S_8$ tension between the CMB and weak lensing can be translated into a tension between geometry and growth, i.e. $\Delta\Omega_M \neq 0$, $\Delta w \neq 0$. In both the split $\Lambda{\rm CDM}$ and $w {\rm CDM}$ cases, our results from combining all data are consistent with $\Delta\Omega_M = 0$ and $\Delta w = 0$. If we omit BAO/RSD data and constrain the split $w {\rm CDM}$ cosmology, we find the data prefers $\Delta w0$ at $4.2\sigma$ evidence. We also find that for both CMB and weak lensing, $\Delta\Omega_M$ and $S_8$ are correlated, with CMB showing a slightly stronger correlation. The general broadening of the contours in our extended model does alleviate the $S_8$ tension, but the allowed nonzero values of $\Delta\Omega_M$ do not encompass the $S_8$ values that would point toward a mismatch between geometry and growth as the origin of the tension.Comment: 27 pages, 10 figures. References updated; matches version published in JCA