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.
The rapidly growing statistical precision of galaxy surveys has led to a need for ever more precise predictions of the observables used to constrain cosmological and galaxy formation models. The ...primary avenue through which such predictions will be obtained is suites of numerical simulations. These simulations must span the relevant model parameter spaces, be large enough to obtain the precision demanded by upcoming data, and be thoroughly validated in order to ensure accuracy. In this paper, we present one such suite of simulations, forming the basis for the Aemulus Project, a collaboration devoted to precision emulation of galaxy survey observables. We have run a set of 75 (1.05 h−1 Gpc)3 simulations with mass resolution and force softening of and 20 h−1 kpc, respectively, in 47 different wCDM cosmologies spanning the range of parameter space allowed by the combination of recent cosmic microwave background, baryon acoustic oscillation, and Type Ia supernova results. We present convergence tests of several observables including spherical overdensity halo mass functions, galaxy projected correlation functions, galaxy clustering in redshift space, and matter and halo correlation functions and power spectra. We show that these statistics are converged to 1% (2%) or to the sample variance of the statistic, whichever is larger, for halos with more than 500 (200) particles, respectively, and scales of r > 200 h−1 kpc in real space or k ∼ 3 h Mpc−1 in harmonic space for z ≤ 1. We find that the dominant source of uncertainty comes from varying the particle loading of the simulations. This leads to large systematic errors for statistics using halos with fewer than 200 particles and scales smaller than k ∼ 4 h Mpc−1. We provide the halo catalogs and snapshots detailed in this work to the community at https://AemulusProject.github.io.
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.
ABSTRACT Empirical methods for connecting galaxies to their dark matter halos have become essential for interpreting measurements of the spatial statistics of galaxies. In this work, we present a ...novel approach for parameterizing the degree of concentration dependence in the abundance matching method. This new parameterization provides a smooth interpolation between two commonly used matching proxies: the peak halo mass and the peak halo maximal circular velocity. This parameterization controls the amount of dependence of galaxy luminosity on halo concentration at a fixed halo mass. Effectively this interpolation scheme enables abundance matching models to have adjustable assembly bias in the resulting galaxy catalogs. With the new DarkSky Simulation, whose larger volume provides lower sample variance, we further show that low-redshift two-point clustering and satellite fraction measurements from SDSS can already provide a joint constraint on this concentration dependence and the scatter within the abundance matching framework.
We present the survey strategy and early results of the "Satellites Around Galactic Analogs" (SAGA) Survey. The SAGA Survey's goal is to measure the distribution of satellite galaxies around 100 ...systems analogous to the Milky Way down to the luminosity of the Leo I dwarf galaxy ( ). We define a Milky Way analog based on K-band luminosity and local environment. Here, we present satellite luminosity functions for eight Milky-Way-analog galaxies between 20 and 40 Mpc. These systems have nearly complete spectroscopic coverage of candidate satellites within the projected host virial radius down to using low-redshift gri color criteria. We have discovered a total of 25 new satellite galaxies: 14 new satellite galaxies meet our formal criteria around our complete host systems, plus 11 additional satellites in either incompletely surveyed hosts or below our formal magnitude limit. Combined with 13 previously known satellites, there are a total of 27 satellites around 8 complete Milky-Way-analog hosts. We find a wide distribution in the number of satellites per host, from 1 to 9, in the luminosity range for which there are 5 Milky Way satellites. Standard abundance matching extrapolated from higher luminosities predicts less scatter between hosts and a steeper luminosity function slope than observed. We find that the majority of satellites (26 of 27) are star-forming. These early results indicate that the Milky Way has a different satellite population than typical in our sample, potentially changing the physical interpretation of measurements based only on the Milky Way's satellite galaxies.
We develop a comprehensive and flexible model for the connection between satellite galaxies and dark matter subhalos in dark matter-only zoom-in simulations of Milky Way (MW)-mass host halos. We ...systematically identify the physical and numerical uncertainties in the galaxy-halo connection and simulations underlying our method, including (i) the influence of host halo properties; (ii) the relationship between satellite luminosities and subhalo properties, including the effects of reionization; (iii) the relationship between satellite and subhalo locations; (iv) the relationship between satellite sizes and subhalo properties, including the effects of tidal stripping; (v) satellite and subhalo disruption due to baryonic effects; and (vi) artificial subhalo disruption and orphan satellites. To illustrate our approach, we fit this model to the luminosity distribution of both classical MW satellites and those discovered in the Sloan Digital Sky Survey by performing realistic mock observations that depend on the luminosity, size, and distance of our predicted satellites, and we infer the total satellite population that will be probed by upcoming surveys. We argue that galaxy size and surface brightness modeling will play a key role in interpreting current and future observations, as the expected number of observable satellites depends sensitively on their surface brightness distribution.
Abstract
We present the Stage II results from the ongoing Satellites Around Galactic Analogs (SAGA) Survey. Upon completion, the SAGA Survey will spectroscopically identify satellite galaxies ...brighter than
M
r
,
o
= −12.3 around 100 Milky Way (MW) analogs at
z
∼ 0.01. In Stage II, we have more than quadrupled the sample size of Stage I, delivering results from 127 satellites around 36 MW analogs with an improved target selection strategy and deep photometric imaging catalogs from the Dark Energy Survey and the Legacy Surveys. We have obtained 25,372 galaxy redshifts, peaking around
z
= 0.2. These data significantly increase spectroscopic coverage for very low redshift objects in 17 <
r
o
< 20.75 around SAGA hosts, creating a unique data set that places the Local Group in a wider context. The number of confirmed satellites per system ranges from zero to nine and correlates with host galaxy and brightest satellite luminosities. We find that the number and luminosities of MW satellites are consistent with being drawn from the same underlying distribution as SAGA systems. The majority of confirmed SAGA satellites are star-forming, and the quenched fraction increases as satellite stellar mass and projected radius from the host galaxy decrease. Overall, the satellite quenched fraction among SAGA systems is lower than that in the Local Group. We compare the luminosity functions and radial distributions of SAGA satellites with theoretical predictions based on cold dark matter simulations and an empirical galaxy–halo connection model and find that the results are broadly in agreement.
We explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)-mass systems. We consider a generic ...self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, w, sets the velocity scale above which momentum transfer due to DM self-interactions becomes inefficient. We perform high-resolution zoom-in simulations of an MW-mass halo for values of w that span scenarios in which self-interactions either between the host and its subhalos or only within subhalos efficiently transfer momentum, and we study the effects of self-interactions on the host halo and on the abundance, radial distribution, orbital dynamics, and density profiles of subhalos in each case. The abundance and properties of surviving subhalos are consistent with being determined primarily by subhalo-host halo interactions. In particular, subhalos on radial orbits in models with larger values of the cross section at the host halo velocity scale are more susceptible to tidal disruption owing to mass loss from ram pressure stripping caused by self-interactions with the host. This mechanism suppresses the abundance of surviving subhalos relative to collisionless DM simulations, with stronger suppression for larger values of w. Thus, probes of subhalo abundance around MW-mass hosts can be used to place upper limits on the self-interaction cross section at velocity scales of , and combining these measurements with the orbital properties and internal dynamics of subhalos may break degeneracies among velocity-dependent SIDM models.
Abstract
Secondary halo bias, commonly known as ‘assembly bias’, is the dependence of halo clustering on a halo property other than mass. This prediction of the Λ Cold Dark Matter cosmology is ...essential to modelling the galaxy distribution to high precision and interpreting clustering measurements. As the name suggests, different manifestations of secondary halo bias have been thought to originate from halo assembly histories. We show conclusively that this is incorrect for cluster-size haloes. We present an up-to-date summary of secondary halo biases of high-mass haloes due to various halo properties including concentration, spin, several proxies of assembly history, and subhalo properties. While concentration, spin, and the abundance and radial distribution of subhaloes exhibit significant secondary biases, properties that directly quantify halo assembly history do not. In fact, the entire assembly histories of haloes in pairs are nearly identical to those of isolated haloes. In general, a global correlation between two halo properties does not predict whether or not these two properties exhibit similar secondary biases. For example, assembly history and concentration (or subhalo abundance) are correlated for both paired and isolated haloes, but follow slightly different conditional distributions in these two cases. This results in a secondary halo bias due to concentration (or subhalo abundance), despite the lack of assembly bias in the strict sense for cluster-size haloes. Due to this complexity, caution must be exercised in using any one halo property as a proxy to study the secondary bias due to another property.
We present the first stable release of Halotools (v0.2), a community-driven Python package designed to build and test models of the galaxy-halo connection. Halotools provides a modular platform for ...creating mock universes of galaxies starting from a catalog of dark matter halos obtained from a cosmological simulation. The package supports many of the common forms used to describe galaxy-halo models: the halo occupation distribution, the conditional luminosity function, abundance matching, and alternatives to these models that include effects such as environmental quenching or variable galaxy assembly bias. Satellite galaxies can be modeled to live in subhalos or to follow custom number density profiles within their halos, including spatial and/or velocity bias with respect to the dark matter profile. The package has an optimized toolkit to make mock observations on a synthetic galaxy population-including galaxy clustering, galaxy-galaxy lensing, galaxy group identification, RSD multipoles, void statistics, pairwise velocities and others-allowing direct comparison to observations. Halotools is object-oriented, enabling complex models to be built from a set of simple, interchangeable components, including those of your own creation. Halotools has an automated testing suite and is exhaustively documented on http://halotools.readthedocs.io, which includes quickstart guides, source code notes and a large collection of tutorials. The documentation is effectively an online textbook on how to build and study empirical models of galaxy formation with Python.