scispace - formally typeset
Open accessJournal ArticleDOI: 10.1093/MNRAS/STAB235

Evidence for galaxy assembly bias in BOSS CMASS redshift-space galaxy correlation function

02 Mar 2021-Monthly Notices of the Royal Astronomical Society (Oxford Academic)-Vol. 502, Iss: 3, pp 3582-3598
Abstract: Building accurate and flexible galaxy-halo connection models is crucial in modeling galaxy clustering on non-linear scales. Recent studies have found that halo concentration by itself cannot capture the full galaxy assembly bias effect and that the local environment of the halo can be an excellent indicator of galaxy assembly bias. In this paper, we propose an extended halo occupation distribution model (HOD) that includes both a concentration-based assembly bias term and an environment-based assembly bias term. We use this model to achieve a good fit (chi^2/DoF = 1.35) on the 2D redshift-space 2-point correlation function (2PCF) of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy sample. We find that the inclusion of both assembly bias terms is strongly favored by the data and the standard 5-parameter HOD is strongly rejected. More interestingly, the redshift-space 2PCF drives the assembly bias parameters in a way that preferentially assigns galaxies to lower mass halos. This results in galaxy-galaxy lensing predictions that are within 1sigma agreement with the observation, alleviating the perceived tension between galaxy clustering and lensing. We also showcase a consistent 3-5sigma preference for a positive environment-based assembly bias that persists over variations in the fit. We speculate that the environmental dependence might be driven by underlying processes such as mergers and feedback, but might also be indicative of a larger halo boundaries such as the splashback radius. Regardless, this work highlights the importance of building flexible galaxy-halo connection models and demonstrates the extra constraining power of the redshift-space 2PCF.

... read more

Topics: Halo occupation distribution (62%), Galaxy (55%), Dark matter (52%) ... show more
Citations
  More

20 results found


Open access
01 Jan 2016-
Abstract: The two-point clustering of dark matter halos is influenced by halo properties besides mass, a phenomenon referred to as halo assembly bias. Using the depth of the gravitational potential well, Vmax; as our secondary halo property, in this paper we present the first study of the scale-dependence assembly bias. In the large-scale linear regime, r & 10Mpc=h; our findings are in keeping with previous results. In particular, at the low-mass end (Mvir < Mcoll 10 12:5 M =h), halos with high-Vmax show stronger large-scale clustering relative to halos with low-Vmax of the same mass; this trend weakens and reverses for Mvir & Mcoll: In the nonlinear regime, assembly bias in low-mass halos exhibits a pronounced scale-dependent “bump” at 500kpc=h 5Mpc=h; a new result. This feature weakens and eventually vanishes for halos of higher mass. We show that this scale-dependent signature can primarily be attributed to a special subpopulation of ejected halos, defined as present-day host halos that were previously members of a higher-mass halo at some point in their past history. A corollary of our results is that galaxy clustering on scales of r 1 2Mpc=h can be impacted by up to 15% by the choice of the halo property used in the halo model, even for stellar mass-limited samples.

... read more

Topics: Halo (61%), Dark matter (51%)

37 Citations


Open accessJournal ArticleDOI: 10.1093/MNRAS/STAB243
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.

... read more

Topics: Halo occupation distribution (63%), Galaxy (53%), Halo (50%) ... show more

17 Citations


Open accessJournal ArticleDOI: 10.1093/MNRAS/STAA3955
Abstract: We develop a series of N-body data challenges, functional to the final analysis of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 (DR16) galaxy sample. The challenges are primarily based on high-fidelity catalogs constructed from the Outer Rim simulation - a large box size realization (3 Gpc/h) characterized by an unprecedented mass resolution, down to 1.85 x 10^9 M_sun/h. We generate synthetic galaxy mocks by populating Outer Rim halos with a variety of halo occupation distribution (HOD) schemes of increasing complexity, spanning different redshift intervals. We then assess the performance of three complementary redshift space distortion (RSD) models in configuration and Fourier space, adopted for the analysis of the complete DR16 eBOSS sample of Luminous Red Galaxies (LRGs). We find that all the methods are mutually consistent, with comparable systematic errors on the Alcock-Paczynski parameters and the growth of structure, and robust to different HOD prescriptions - thus validating the robustness of the models and the pipelines used for the baryon acoustic oscillation (BAO) and full shape clustering analysis. In particular, all the techniques are able to recover alpha_par and alpha_perp to within 0.9%, and fsig8 to within 1.5%. As a by-product of our work, we are also able to gain interesting insights on the galaxy-halo connection. Our study is relevant for the final eBOSS DR16 `consensus cosmology', as the systematic error budget is informed by testing the results of analyses against these high-resolution mocks. In addition, it is also useful for future large-volume surveys, since similar mock-making techniques and systematic corrections can be readily extended to model for instance the DESI galaxy sample.

... read more

13 Citations


Open accessJournal ArticleDOI: 10.1093/MNRAS/STAB1613
Rafael C. Nunes1, Sunny Vagnozzi2Institutions (2)
Abstract: Within the $\Lambda$CDM model, measurements from recent Cosmic Microwave Background (CMB) and weak lensing (WL) surveys have uncovered a $\sim 3\sigma$ disagreement in the inferred value of the parameter $S_8 \equiv \sigma_8\sqrt{\Omega_m/0.3}$, quantifying the amplitude of late-time matter fluctuations. Before questioning whether the $S_8$ discrepancy calls for new physics, it is important to assess the view of measurements other than CMB and WL ones on the discrepancy. Here, we examine the role of measurements of the growth rate $f(z)$ in arbitrating the $S_8$ discrepancy, considering measurements of $f\sigma_8(z)$ from Redshift-Space Distortions (RSD). Our baseline analysis combines RSD measurements with geometrical measurements from Baryon Acoustic Oscillations (BAO) and Type Ia Supernovae (SNeIa), given the key role of the latter in constraining $\Omega_m$. From this combination and within the $\Lambda$CDM model we find $S_8 = 0.762^{+0.030}_{-0.025}$, and quantify the agreement between RSD+BAO+SNeIa and \textit{Planck} to be at the $2.2\sigma$ level: the mild disagreement is therefore compatible with a statistical fluctuation. We discuss combinations of RSD measurements with other datasets, including the $E_G$ statistic. This combination increases the discrepancy with \textit{Planck}, but we deem it significantly less robust. Our earlier results are stable against an extension where we allow the dark energy equation of state $w$ to vary. We conclude that, from the point of view of combined growth rate and geometrical measurements, there are hints, but no strong evidence yet, for the \textit{Planck}$\Lambda$CDM cosmology over-predicting the amplitude of matter fluctuations at redshifts $z \lesssim 1$. From this perspective, it might therefore still be premature to claim the need for new physics from the $S_8$ discrepancy.

... read more

10 Citations


Open accessJournal ArticleDOI: 10.1093/MNRAS/STAB189
Abstract: The canonical $\Lambda$CDM cosmological model makes precise predictions for the clustering and lensing properties of galaxies. It has been shown that the lensing amplitude of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) is lower than expected given their clustering properties. We present new measurements and modelling of galaxies in the BOSS LOWZ sample. We focus on the radial and stellar mass dependence of the lensing amplitude mis-match. We find an amplitude mis-match of around $35\%$ when assuming $\Lambda$CDM with Planck Cosmological Microwave Background (CMB) constraints. This offset is independent of halo mass and radial scale in the range $M_{\rm halo}\sim 10^{13.3} - 10^{13.9} h^{-1} M_\odot$ and $r=0.1 - 60 \, h^{-1} \mathrm{Mpc}$ ($k \approx 0.05 - 20 \, h \, {\rm Mpc}^{-1}$). The observation that the offset is both mass and scale independent places important constraints on the degree to which astrophysical processes (baryonic effects, assembly bias) can fully explain the effect. This scale independence also suggests that the "lensing is low" effect on small and large radial scales probably have the same physical origin. Resolutions based on new physics require a nearly uniform suppression, relative to $\Lambda$CDM predictions, of the amplitude of matter fluctuations on these scales. The possible causes of this are tightly constrained by measurements of the CMB and of the low-redshift expansion history.

... read more

Topics: Dark matter (52%), Halo (50%)

10 Citations


References
  More

123 results found


Open accessJournal ArticleDOI: 10.1051/0004-6361/201525830
Peter A. R. Ade1, Nabila Aghanim2, Monique Arnaud3, M. Ashdown4  +334 moreInstitutions (82)
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.

... read more

Topics: Planck energy (70%), Planck particle (68%), Planck temperature (68%) ... show more

10,334 Citations


Open accessJournal ArticleDOI: 10.1051/0004-6361/201525830
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.

... read more

Topics: Planck temperature (67%), Planck (60%), Cosmic microwave background (60%) ... show more

9,349 Citations


Open accessJournal ArticleDOI: 10.1086/304888
Abstract: We use high-resolution N-body simulations to study the equilibrium density profiles of dark matter halos in hierarchically clustering universes. We find that all such profiles have the same shape, independent of the halo mass, the initial density fluctuation spectrum, and the values of the cosmological parameters. Spherically averaged equilibrium profiles are well fitted over two decades in radius by a simple formula originally proposed to describe the structure of galaxy clusters in a cold dark matter universe. In any particular cosmology, the two scale parameters of the fit, the halo mass and its characteristic density, are strongly correlated. Low-mass halos are significantly denser than more massive systems, a correlation that reflects the higher collapse redshift of small halos. The characteristic density of an equilibrium halo is proportional to the density of the universe at the time it was assembled. A suitable definition of this assembly time allows the same proportionality constant to be used for all the cosmologies that we have tested. We compare our results with previous work on halo density profiles and show that there is good agreement. We also provide a step-by-step analytic procedure, based on the Press-Schechter formalism, that allows accurate equilibrium profiles to be calculated as a function of mass in any hierarchical model.

... read more

Topics: Halo mass function (72%), Halo occupation distribution (66%), Cuspy halo problem (64%) ... show more

9,149 Citations


Open accessJournal ArticleDOI: 10.1086/670067
Daniel Foreman-Mackey1, David W. Hogg, Dustin Lang2, Dustin Lang3  +1 moreInstitutions (3)
Abstract: We introduce a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010). The code is open source and has already been used in several published projects in the astrophysics literature. The algorithm behind emcee has several advantages over traditional MCMC sampling methods and it has excellent performance as measured by the autocorrelation time (or function calls per independent sample). One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to ~N2 for a traditional algorithm in an N-dimensional parameter space. In this document, we describe the algorithm and the details of our implementation. Exploiting the parallelism of the ensemble method, emcee permits any user to take advantage of multiple CPU cores without extra effort. The code is available online at http://dan.iel.fm/emcee under the GNU General Public License v2.

... read more

Topics: Markov chain Monte Carlo (52%)

6,914 Citations


Open accessJournal ArticleDOI: 10.1086/670067
Daniel Foreman-Mackey1, David W. Hogg, Dustin Lang2, Dustin Lang3  +1 moreInstitutions (3)
Abstract: We introduce a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010). The code is open source and has already been used in several published projects in the astrophysics literature. The algorithm behind emcee has several advantages over traditional MCMC sampling methods and it has excellent performance as measured by the autocorrelation time (or function calls per independent sample). One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to $\sim N^2$ for a traditional algorithm in an N-dimensional parameter space. In this document, we describe the algorithm and the details of our implementation and API. Exploiting the parallelism of the ensemble method, emcee permits any user to take advantage of multiple CPU cores without extra effort. The code is available online at this http URL under the MIT License.

... read more

Topics: Markov chain Monte Carlo (53%)

3,846 Citations


Performance
Metrics
No. of citations received by the Paper in previous years
YearCitations
202117
20202
20161