Abstract
We apply a new galaxy group-finder to the Main Galaxy Sample of the SDSS. This algorithm introduces new freedom to assign halos to galaxies that is self-calibrated by comparing the catalog ...to complementary data. These include galaxy clustering data and measurements of the total satellite luminosity from deep-imaging data. We present constraints on the galaxy-halo connection for star-forming and quiescent populations. The results of the self-calibrated group catalog differ in several key ways from previous group catalogs and halo-occupation analyses. The transition halo mass scale, where half of the halos contain quiescent central galaxies, is at
M
h
∼ 10
12.4
h
−1
M
⊙
, significantly higher than other constraints. Additionally, the width of the transition from predominantly star-forming halos to quiescent halos occurs over a narrower range in halo mass. Quiescent central galaxies in low-mass halos are significantly more massive than star-forming centrals at the same halo mass, but this difference reverses above the transition halo mass. We find that the scatter in
log
M
*
at fixed
M
h
is ∼0.2 dex for massive halos, in agreement with previous estimates, but rises sharply at lower halo masses. The halo masses assigned by the group catalog are in good agreement with weak-lensing estimates for star-forming and quiescent central galaxies. We discuss possible improvements to the algorithm made clear by this first application to data. The group catalog is made publicly available.
In our modern understanding of galaxy formation, every galaxy forms within a dark matter halo. The formation and growth of galaxies over time is connected to the growth of the halos in which they ...form. The advent of large galaxy surveys as well as high-resolution cosmological simulations has provided a new window into the statistical relationship between galaxies and halos and its evolution. Here, we define this galaxy-halo connection as the multivariate distribution of galaxy and halo properties that can be derived from observations and simulations. This galaxy-halo connection provides a key test of physical galaxy-formation models; it also plays an essential role in constraints of cosmological models using galaxy surveys and in elucidating the properties of dark matter using galaxies. We review techniques for inferring the galaxy-halo connection and the insights that have arisen from these approaches. Some things we have learned are that galaxy-formation efficiency is a strong function of halo mass; at its peak in halos around a pivot halo mass of 10
12
M
, less than 20% of the available baryons have turned into stars by the present day; the intrinsic scatter in galaxy stellar mass is small, less than 0.2 dex at a given halo mass above this pivot mass; below this pivot mass galaxy stellar mass is a strong function of halo mass; the majority of stars over cosmic time were formed in a narrow region around this pivot mass. We also highlight key open questions about how galaxies and halos are connected, including understanding the correlations with secondary properties and the connection of these properties to galaxy clustering.
Abstract
We describe an extension of the halo-based galaxy group-finding algorithm. We add freedom to the algorithm in order to more accurately determine which galaxies are central and which are ...satellites, and to provide unbiased estimates of halo masses. We focus on determination of the galaxy-halo relations for star-forming and quiescent galaxies. The added freedom in the group-finding algorithm is self-calibrated using observations of color-dependent galaxy clustering, as well as measurements of the total satellite luminosity in deep imaging data around stacked samples of spectroscopic central galaxies,
L
sat
. We test this approach on a series of mocks that vary the galaxy-halo connection, including one mock constructed from Universe Machine results. Our self-calibrated algorithm shows marked improvement over previous methods in estimating the color-dependent satellite fraction of galaxies. It reduces the error in
log
M
h
for central galaxies by over a factor of two, to ≲0.2 dex. Through the
L
sat
data, it can quantify differences in the luminosity-to-halo mass relations for star-forming and quiescent galaxies, even for groups with only one spectroscopic member. Thus, whereas previous algorithms cannot constrain the scatter in
L
gal
at fixed
M
h
, the self-calibration technique can provide a robust lower limit to this scatter.
Using galaxy group/cluster catalogues created from the Sloan Digital Sky Survey Data Release 7, we examine in detail the specific star formation rate (SSFR) distribution of satellite galaxies and its ...dependence on stellar mass, host halo mass and halo-centric radius. All galaxies, regardless of central satellite designation, exhibit a similar bimodal SSFR distribution, with a strong break at SSFR ≈ 10−11 yr−1 and the same high SSFR peak; in no regime is there ever an excess of galaxies in the 'green valley'. Satellite galaxies are simply more likely to lie on the quenched ('red sequence') side of the SSFR distribution. Furthermore, the satellite quenched fraction excess above the field galaxy value is nearly independent of galaxy stellar mass. An enhanced quenched fraction for satellites persists in groups with halo masses down to 3 × 1011 M⊙ and increases strongly with halo mass and towards halo centre. We find no detectable quenching enhancement for galaxies beyond ∼2 R
vir around massive clusters once these galaxies have been decomposed into centrals and satellites. These trends imply that (1) galaxies experience no significant environmental effects until they cross within ∼R
vir of a more massive host halo; (2) after this, star formation in active satellites continues to evolve in the same manner as active central galaxies for several Gyr; and (3) once begun, satellite star formation quenching occurs rapidly. These results place strong constraints on satellite-specific quenching mechanisms, as we will discuss further in companion papers.
Abstract
We use the scatter in the stellar-to-halo mass relation to constrain galaxy evolution models. If the efficiency of converting accreted baryons into stars varies with time, haloes of the same ...present-day mass but different formation histories will have different z = 0 galaxy stellar mass. This is one of the sources of scatter in stellar mass at fixed halo mass, σlog M
*. For massive haloes that undergo rapid quenching of star formation at z ∼ 2, different mechanisms that trigger this quenching yield different values of σlog M*. We use this framework to test various models in which quenching begins after a galaxy crosses a threshold in one of the following physical quantities: redshift, halo mass, stellar mass and stellar-to-halo mass ratio. Our model is highly idealized, with other sources of scatter likely to arise as more physics is included. Thus, our test is whether a model can produce scatter lower than observational bounds, leaving room for other sources. Recent measurements find σlog M* = 0.16 dex for 1011 M⊙ galaxies. Under the assumption that the threshold is constant with time, such a low value of σlog M* rules out all of these models with the exception of quenching by a stellar mass threshold. Most physical quantities, such as metallicity, will increase scatter if they are uncorrelated with halo formation history. Thus, to decrease the scatter of a given model, galaxy properties would correlate tightly with formation history, creating testable predictions for their clustering. Understanding why σlog M* is so small may be key to understanding the physics of galaxy formation.
Satellite galaxies in groups and clusters are more likely to have low star formation rates (SFRs) and lie on the 'red sequence' than central ('field') galaxies. Using galaxy group/cluster catalogues ...from the Sloan Digital Sky Survey Data Release 7, together with a high-resolution, cosmological N-body simulation to track satellite orbits, we examine the star formation histories and quenching time-scales of satellites of M
star > 5 × 109 M at z 0. We first explore satellite infall histories: group preprocessing and ejected orbits are critical aspects of satellite evolution, and properly accounting for these, satellite infall typically occurred at z ∼ 0.5, or ∼5 Gyr ago. To obtain accurate initial conditions for the SFRs of satellites at their time of first infall, we construct an empirical parametrization for the evolution of central galaxy SFRs and quiescent fractions. With this, we constrain the importance and efficiency of satellite quenching as a function of satellite and host halo mass, finding that satellite quenching is the dominant process for building up all quiescent galaxies at M
star < 1010 M. We then constrain satellite star formation histories, finding a 'delayed-then-rapid' quenching scenario: satellite SFRs evolve unaffected for 2-4 Gyr after infall, after which star formation quenches rapidly, with an e-folding time of <0.8 Gyr. These quenching time-scales are shorter for more massive satellites but do not depend on host halo mass: the observed increase in the satellite quiescent fraction with halo mass arises simply because of satellites quenching in a lower mass group prior to infall (group preprocessing), which is responsible for up to half of quenched satellites in massive clusters. Because of the long time delay before quenching starts, satellites experience significant stellar mass growth after infall, nearly identical to central galaxies. This fact provides key physical insight into the subhalo abundance matching method.
Existing models for the dependence of the halo mass function on cosmological parameters will become a limiting source of systematic uncertainty for cluster cosmology in the near future. We present a ...halo mass function emulator and demonstrate improved accuracy relative to state-of-the-art analytic models. In this work, mass is defined using an overdensity criteria of 200 relative to the mean background density. Our emulator is constructed from the Aemulus simulations, a suite of 40 N-body simulations with snapshots from z = 3 to z = 0. These simulations cover the flat wCDM parameter space allowed by recent cosmic microwave background, baryon acoustic oscillation and SNe Ia results, varying the parameters w, m, b, 8, Neff, ns, and H0. We validate our emulator using five realizations of seven different cosmologies, for a total of 35 test simulations. These test simulations were not used in constructing the emulator, and were run with fully independent initial conditions. We use our test simulations to characterize the modeling uncertainty of the emulator, and introduce a novel way of marginalizing over the associated systematic uncertainty. We confirm nonuniversality in our halo mass function emulator as a function of both cosmological parameters and redshift. Our emulator achieves better than 1% precision over much of the relevant parameter space, and we demonstrate that the systematic uncertainty in our emulator will remain a negligible source of error for cluster abundance studies through at least the LSST Year 1 data set.
We perform the first fit to the anisotropic clustering of Sloan Digital Sky Survey III CMASS data release 10 galaxies on scales of ∼0.8–32 h
−1 Mpc. A standard halo occupation distribution model ...evaluated near the best-fitting Planck Λ cold dark matter (ΛCDM) cosmology provides a good fit to the observed anisotropic clustering, and implies a normalization for the peculiar velocity field of M ∼ 2 × 1013
h
−1 M⊙ haloes of fσ8(z = 0.57) = 0.450 ± 0.011. Since this constraint includes both quasi-linear and non-linear scales, it should severely constrain modified gravity models that enhance pairwise infall velocities on these scales. Though model dependent, our measurement represents a factor of 2.5 improvement in precision over the analysis of DR11 on large scales, fσ8(z = 0.57) = 0.447 ± 0.028, and is the tightest single constraint on the growth rate of cosmic structure to date. Our measurement is consistent with the Planck ΛCDM prediction of 0.480 ± 0.010 at the ∼1.9σ level. Assuming a halo mass function evaluated at the best-fitting Planck
cosmology, we also find that 10 per cent of CMASS galaxies are satellites in haloes of mass M ∼ 6 × 1013
h
−1 M⊙. While none of our tests and model generalizations indicate systematic errors due to an insufficiently detailed model of the galaxy–halo connection, the precision of these first results warrant further investigation into the modelling uncertainties and degeneracies with cosmological parameters.
Using the N-body simulations of the Aemulus Project, we construct an emulator for the nonlinear clustering of galaxies in real and redshift space. We construct our model of galaxy bias using the halo ...occupation framework, accounting for possible velocity bias. The model includes 15 parameters, including both cosmological and galaxy bias parameters. We demonstrate that our emulator achieves ∼1% precision at the scales of interest, 0.1 h−1 Mpc < r < 10 h−1 Mpc, and recovers the true cosmology when tested against independent simulations. Our primary parameters of interest are related to the growth rate of structure, f, and its degenerate combination, f 8. Using this emulator, we show that the constraining power on these parameters monotonically increases as smaller scales are included in the analysis, all the way down to 0.1 h−1 Mpc. For a BOSS-like survey, the constraints on f 8 from r < 30 h−1 Mpc scales alone are nearly a factor of two tighter than those from the fiducial BOSS analysis of redshift-space clustering using perturbation theory at larger scales. The combination of real- and redshift-space clustering allows us to break the degeneracy between f and 8, yielding an 11% constraint on f alone for a BOSS-like analysis. The current Aemulus simulations limit this model to surveys of massive galaxies. Future simulations will allow this framework to be extended to all galaxy target types, including emission-line galaxies.
Central galaxies make up the majority of the galaxy population, including the majority of the quiescent population at . Thus, the mechanism(s) responsible for quenching central galaxies play a ...crucial role in galaxy evolution as whole. We combine a high-resolution cosmological N-body simulation with observed evolutionary trends of the "star formation main sequence," quiescent fraction, and stellar mass function at to construct a model that statistically tracks the star formation histories and quenching of central galaxies. Comparing this model to the distribution of central galaxy star formation rates in a group catalog of the SDSS Data Release 7, we constrain the timescales over which physical processes cease star formation in central galaxies. Over the stellar mass range we infer quenching e-folding times that span 1.5-0.5 Gyr with more massive central galaxies quenching faster. For , this implies a total migration time of from the star formation main sequence to quiescence. Compared to satellites, central galaxies take longer to quench their star formation, suggesting that different mechanisms are responsible for quenching centrals versus satellites. Finally, the central galaxy quenching timescale we infer provides key constraints for proposed star formation quenching mechanisms. Our timescale is generally consistent with gas depletion timescales predicted by quenching through strangulation. However, the exact physical mechanism(s) responsible for this remain unclear.