The galaxy–halo connection of emission-line galaxies in IllustrisTNG
02 Mar 2021-Monthly Notices of the Royal Astronomical Society (Oxford Academic)-Vol. 502, Iss: 3, pp 3599-3617
TL;DR: In this article, the authors employ the hydrodynamical simulation IllustrisTNG-300-1 to explore the halo occupation distribution and environmental dependence of luminous star-forming emission-line galaxies (ELGs) at $z \sim 1$.
Abstract: We employ the hydrodynamical simulation IllustrisTNG-300-1 to explore the halo occupation distribution (HOD) and environmental dependence of luminous star-forming emission-line galaxies (ELGs) at $z \sim 1$. Such galaxies are key targets for current and upcoming cosmological surveys. We select model galaxies through cuts in colour-colour space allowing for a direct comparison with the Extended Baryon Oscillation Spectroscopic Survey and the Dark Energy Spectroscopic Instrument (DESI) surveys and then compare them with galaxies selected based on specific star-formation rate (sSFR) and stellar mass. We demonstrate that the ELG populations are twice more likely to reside in lower-density regions (sheets) compared with the mass-selected populations and twice less likely to occupy the densest regions of the cosmic web (knots). We also show that the colour-selected and sSFR-selected ELGs exhibit very similar occupation and clustering statistics, finding that the agreement is best for lower redshifts. In contrast with the mass-selected sample, the occupation of haloes by a central ELG peaks at $\sim$20\%. We furthermore explore the dependence of the HOD and the auto-correlation on environment, noticing that at fixed halo mass, galaxies in high-density regions cluster about 10 times more strongly than low-density ones. This result suggests that we should model carefully the galaxy-halo relation and implement assembly bias effects into our models (estimated at $\sim$4\% of the clustering of the DESI colour-selected sample at $z = 0.8$). Finally, we apply a simple mock recipe to recover the clustering on large scales ($r \gtrsim 1 \ {\rm Mpc}/h$) to within 1\% by augmenting the HOD model with an environment dependence, demonstrating the power of adopting flexible population models.
Citations
More filters
01 Jan 2019
TL;DR: The DESI Legacy Imaging Surveys (http://legacysurvey.org/) project is a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image ≈14,000 deg2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory.
Abstract: Author(s): Huang, L; Dey, A; Schlegel, DJ; Lang, D; Blum, R; Burleigh, K; Fan, X; Findlay, JR; Finkbeiner, D; Herrera, D; Juneau, S; Landriau, M; Levi, M; McGreer, I; Meisner, A; Myers, AD; Moustakas, J; Nugent, P; Patej, A; Schlafly, EF; Walker, AR; Valdes, F; Weaver, BA; Yeche, C; Zou, H; Zhou, X; Abareshi, B; Abbott, TMC; Abolfathi, B; Aguilera, C; Alam, S; Allen, L; Alvarez, A; Annis, J; Ansarinejad, B; Aubert, M; Beechert, J; Bell, EF; Benzvi, SY; Beutler, F; Bielby, RM; Bolton, AS; Briceno, C; Buckley-Geer, EJ; Butler, K; Calamida, A; Carlberg, RG; Carter, P; Casas, R; Castander, FJ; Choi, Y; Comparat, J; Cukanovaite, E; Delubac, T; Devries, K; Dey, S; Dhungana, G; Dickinson, M; Ding, Z; Donaldson, JB; Duan, Y; Duckworth, CJ; Eftekharzadeh, S; Eisenstein, DJ; Etourneau, T; Fagrelius, PA; Farihi, J; Fitzpatrick, M; Font-Ribera, A; Fulmer, L; Gansicke, BT; Gaztanaga, E; George, K; Gerdes, DW; A Gontcho, SG; Gorgoni, C; Green, G; Guy, J; Harmer, D; Hernandez, M; Honscheid, K; Huang, LW; James, DJ; Jannuzi, BT; Jiang, L | Abstract: © 2019. The American Astronomical Society. All rights reserved.. The DESI Legacy Imaging Surveys (http://legacysurvey.org/) are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image ≈14,000 deg2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12, and 22 μm) observed by the Wide-field Infrared Survey Explorer satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.
193 citations
TL;DR: In this paper, an analytical forward model based on perturbation theory was developed to predict the redshift-space galaxy overdensity at the field level given a realization of the initial conditions.
Abstract: We develop an analytical forward model based on perturbation theory to predict the redshift-space galaxy overdensity at the field level given a realization of the initial conditions. We find that the residual noise between the model and simulated galaxy density has a power spectrum that is white on large scales, with size comparable to the shot noise. In the mildly nonlinear regime, we see a $k^2\mu^2$ correction to the noise power spectrum, corresponding to larger noise along the line of sight and on smaller scales. The parametric form of this correction has been predicted on theoretical grounds before, and our simulations provide important confirmation of its presence. We have also modeled the galaxy velocity at the field-level and compared it against simulated galaxy velocities, finding that about $10\%$ of the galaxies are responsible for half of the rms velocity residual for our simulated galaxy sample.
30 citations
Additional excerpts
...[11, 12] for two recent studies....
[...]
TL;DR: In this paper, the authors investigate the impact of long-wavelength perturbations on the response of the HODs to longwavelength observations and find that the effect of these perturbation on the bias of the gas distribution is not negligible in general.
Abstract: Halo occupation distribution (HOD) models describe the number of galaxies that reside in different haloes, and are widely used in galaxy-halo connection studies using the halo model (HM). Here, we introduce and study HOD response functions $R_\mathcal{O}^g$ that describe the response of the HODs to long-wavelength perturbations $\mathcal{O}$. The linear galaxy bias parameters $b_\mathcal{O}^g$ are a weighted version of $b_\mathcal{O}^h + R_\mathcal{O}^g$, where $b_\mathcal{O}^h$ is the halo bias, but the contribution from $R_\mathcal{O}^g$ is routinely ignored in the literature. We investigate the impact of this by measuring the $R_\mathcal{O}^g$ in separate universe simulations of the IllustrisTNG model for three types of perturbations: total matter perturbations, $\mathcal{O}=\delta_m$; baryon-CDM compensated isocurvature perturbations, $\mathcal{O}=\sigma$; and potential perturbations with local primordial non-Gaussianity, $\mathcal{O}\propto f_{\rm NL}\phi$. Our main takeaway message is that the $R_\mathcal{O}^g$ are not negligible in general and their size should be estimated on a case-by-case basis. For stellar-mass selected galaxies, the responses $R_\phi^g$ and $R_\sigma^g$ are sizeable and cannot be neglected in HM calculations of the bias parameters $b_\phi^g$ and $b_\sigma^g$; this is relevant to constrain inflation using galaxies. On the other hand, we do not detect a strong impact of the HOD response $R_1^g$ on the linear galaxy bias $b_1^g$. These results can be explained by the impact that the perturbations have on stellar-to-total-mass relations. We also look into the impact on the bias of the gas distribution and find similar conclusions. We show that a single extra parameter describing the overall amplitude of $R_\mathcal{O}^g$ recovers the measured $b_\mathcal{O}^g$ well, which indicates that $R_\mathcal{O}^g$ can be easily added to HM/HOD studies as a new ingredient.
19 citations
TL;DR: In this paper, the authors compare large scale structure observables for stellar mass selected samples at $z = 0, as predicted by two galaxy models, the hydrodynamical simulation IllustrisTNG and the Santa-Cruz semi-analytic model (SC-SAM).
Abstract: In this work, we compare large scale structure observables for stellar mass selected samples at $z=0$, as predicted by two galaxy models, the hydrodynamical simulation IllustrisTNG and the Santa-Cruz semi-analytic model (SC-SAM). Although both models have been independently calibrated to match observations, rather than each other, we find good agreement between the two models for two-point clustering and galaxy assembly bias signatures. The models also show a qualitatively similar response of occupancy and clustering to secondary halo paramaters other than mass, such as formation history and concentration, although with some quantitative differences. Thus, our results demonstrate that the galaxy-halo relationships in SC-SAM and TNG are quite similar to first order. However, we also find areas in which the models differ. For example, we note a strong correlation between halo gas content and environment in TNG, which is lacking in the SC-SAM, as well as differences in the occupancy predictions for low-mass haloes. Moreover, we show that higher-order statistics, such as cumulants of the density field, help to accurately describe the galaxy distribution and discriminate between models that show degenerate behavior for two-point statistics. Our results suggest that SAMs are a promising cost-effective and intuitive method for generating mock catalogues for next generation cosmological surveys.
15 citations
TL;DR: This article examined the impact of baryonic physics on the halo distribution in hydrodynamic simulations (Illustris, IllustrisTNG, and EAGLE), particularly with regards to how it differs from that in dark matter only (DMO) simulations.
Abstract: We examine the impact of baryonic physics on the halo distribution in hydrodynamic simulations (Illustris, IllustrisTNG, and EAGLE), particularly with regards to how it differs from that in dark matter only (DMO) simulations. We find that, in general, DMO simulations produce halo mass functions (HMF) that are shifted to higher halo masses than their hydrodynamic counterparts, due to the lack of baryonic physics. However, the exact nature of this mass shift is a complex function of mass, halo definition, redshift, and larger-scale environment, and it also depends on the specifics of the baryonic physics implemented in the simulation. We present fitting formulae for the corrections one would need to apply to each DMO halo catalogue in order to reproduce the HMF found in its hydrodynamic counterpart. We provide these formulae for all three simulations, for five different halo definitions at redshifts 0, 1, and 2. Additionally, we explore the dependence on environment of this HMF discrepancy, and find that, in most cases, halos in low density environments are slightly more impacted by baryonic physics than halos in high density environments. We thus also provide environment-dependent mass correction formulae that can reproduce the conditional, as well as global, HMF. We show that our mass corrections also repair the large-scale clustering of halos, though the environment-dependent corrections are required to achieve an accuracy better than 2%. Finally, we examine the impact of baryonic physics on the halo mass - concentration relation, and find that its slope in hydrodynamic simulations is consistent with that in DMO simulations. Ultimately, we recommend that any future work relying on DMO halo catalogues incorporate our mass corrections to test the robustness of their results to baryonic effects.
14 citations
References
More filters
TL;DR: In this article, the authors present a cosmological analysis based on full-mission Planck observations of temperature and polarization anisotropies of the cosmic microwave background (CMB) radiation.
Abstract: This paper presents cosmological results based on full-mission Planck observations of temperature and polarization anisotropies of the cosmic microwave background (CMB) radiation. Our results are in very good agreement with the 2013 analysis of the Planck nominal-mission temperature data, but with increased precision. The temperature and polarization power spectra are consistent with the standard spatially-flat 6-parameter ΛCDM cosmology with a power-law spectrum of adiabatic scalar perturbations (denoted “base ΛCDM” in this paper). From the Planck temperature data combined with Planck lensing, for this cosmology we find a Hubble constant, H0 = (67.8 ± 0.9) km s-1Mpc-1, a matter density parameter Ωm = 0.308 ± 0.012, and a tilted scalar spectral index with ns = 0.968 ± 0.006, consistent with the 2013 analysis. Note that in this abstract we quote 68% confidence limits on measured parameters and 95% upper limits on other parameters. We present the first results of polarization measurements with the Low Frequency Instrument at large angular scales. Combined with the Planck temperature and lensing data, these measurements give a reionization optical depth of τ = 0.066 ± 0.016, corresponding to a reionization redshift of . These results are consistent with those from WMAP polarization measurements cleaned for dust emission using 353-GHz polarization maps from the High Frequency Instrument. We find no evidence for any departure from base ΛCDM in the neutrino sector of the theory; for example, combining Planck observations with other astrophysical data we find Neff = 3.15 ± 0.23 for the effective number of relativistic degrees of freedom, consistent with the value Neff = 3.046 of the Standard Model of particle physics. The sum of neutrino masses is constrained to ∑ mν < 0.23 eV. The spatial curvature of our Universe is found to be very close to zero, with | ΩK | < 0.005. Adding a tensor component as a single-parameter extension to base ΛCDM we find an upper limit on the tensor-to-scalar ratio of r0.002< 0.11, consistent with the Planck 2013 results and consistent with the B-mode polarization constraints from a joint analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP B-mode data to our analysis leads to a tighter constraint of r0.002 < 0.09 and disfavours inflationarymodels with a V(φ) ∝ φ2 potential. The addition of Planck polarization data leads to strong constraints on deviations from a purely adiabatic spectrum of fluctuations. We find no evidence for any contribution from isocurvature perturbations or from cosmic defects. Combining Planck data with other astrophysical data, including Type Ia supernovae, the equation of state of dark energy is constrained to w = −1.006 ± 0.045, consistent with the expected value for a cosmological constant. The standard big bang nucleosynthesis predictions for the helium and deuterium abundances for the best-fit Planck base ΛCDM cosmology are in excellent agreement with observations. We also constraints on annihilating dark matter and on possible deviations from the standard recombination history. In neither case do we find no evidence for new physics. The Planck results for base ΛCDM are in good agreement with baryon acoustic oscillation data and with the JLA sample of Type Ia supernovae. However, as in the 2013 analysis, the amplitude of the fluctuation spectrum is found to be higher than inferred from some analyses of rich cluster counts and weak gravitational lensing. We show that these tensions cannot easily be resolved with simple modifications of the base ΛCDM cosmology. Apart from these tensions, the base ΛCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.
10,728 citations
TL;DR: In this paper, the authors present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB, which are consistent with the six-parameter inflationary LCDM cosmology.
Abstract: We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.
9,745 citations
3,633 citations
"The galaxy–halo connection of emiss..." refers background in this paper
...Baryonic matter sinks to the center of their gravitational potential wells and condensation of cold gas allows for galaxies to form and evolve (White & Rees 1978)....
[...]
TL;DR: The Virgo Consortium's EAGLE project as discussed by the authors is a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes, where thermal energy is injected into the gas, allowing winds to develop without predetermined speed or mass loading factors.
Abstract: We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to ≲0.2 dex over the full mass range, 108
2,828 citations
"The galaxy–halo connection of emiss..." refers background in this paper
...…simulations, on the other hand, are an example of an ab initio approach to gaining insight into the formation and evolution of galaxies (e.g. Vogelsberger et al. 2014; Schaye et al. 2015), incorporating baryonic effects such as stellar wind, supernova feedback, gas cooling, and black hole feedback....
[...]
TL;DR: In this article, a robust method to constrain average galaxy star formation rates, star formation histories (SFHs), and the intracluster light (ICL) as a function of halo mass is presented.
Abstract: We present a robust method to constrain average galaxy star formation rates (SFRs), star formation histories (SFHs), and the intracluster light (ICL) as a function of halo mass. Our results are consistent with observed galaxy stellar mass functions, specific star formation rates (SSFRs), and cosmic star formation rates (CSFRs) from z = 0 to z = 8. We consider the effects of a wide range of uncertainties on our results, including those affecting stellar masses, SFRs, and the halo mass function at the heart of our analysis. As they are relevant to our method, we also present new calibrations of the dark matter halo mass function, halo mass accretion histories, and halo-subhalo merger rates out to z = 8. We also provide new compilations of CSFRs and SSFRs; more recent measurements are now consistent with the buildup of the cosmic stellar mass density at all redshifts. Implications of our work include: halos near 1012 M ☉ are the most efficient at forming stars at all redshifts, the baryon conversion efficiency of massive halos drops markedly after z ~ 2.5 (consistent with theories of cold-mode accretion), the ICL for massive galaxies is expected to be significant out to at least z ~ 1-1.5, and dwarf galaxies at low redshifts have higher stellar mass to halo mass ratios than previous expectations and form later than in most theoretical models. Finally, we provide new fitting formulae for SFHs that are more accurate than the standard declining tau model. Our approach places a wide variety of observations relating to the SFH of galaxies into a self-consistent framework based on the modern understanding of structure formation in ΛCDM. Constraints on the stellar mass-halo mass relationship and SFRs are available for download online.
2,085 citations
"The galaxy–halo connection of emiss..." refers background in this paper
...G A ] 1 0 N ov 2 02 halos have been an area of keen interest (Behroozi et al. 2013; Tacchella et al. 2013; Moster et al. 2013, 2018; Tacchella et al. 2018; Behroozi et al. 2019)....
[...]