The bias of the submillimetre galaxy population: SMGs are poor tracers of the most massive structures in the z ~ 2 Universe

It is often claimed that overdensities of (or even individual bright) submillimetre-selected galaxies (SMGs) trace the assembly of the most-massive dark matter structures in the Universe. We test this claim by performing a counts-in-cells analysis of mock SMG catalogues derived from the Bolshoi cosmological simulation to investigate how well SMG associations trace the underlying dark matter structure. We find that SMGs exhibit a relatively complex bias: some regions of high SMG overdensity are underdense in terms of dark matter mass, and some regions of high dark matter overdensity contain no SMGs. Because of their rarity, Poisson noise causes scatter in the SMG overdensity at fixed dark matter overdensity. Consequently, rich associations of less-luminous, more-abundant galaxies (i.e. Lyman-break galaxy analogues) trace the highest dark matter overdensities much better than SMGs. Even on average, SMG associations are relatively poor tracers of the most significant dark matter overdensities because of 'downsizing': at z < ~2.5, the most-massive galaxies that reside in the highest dark matter overdensities have already had their star formation quenched and are thus no longer SMGs. At a given redshift, of the 10 per cent most-massive overdensities, only ~25 per cent contain at least one SMG, and less than a few per cent contain more than one SMG.Comment: 6 pages, 3 figures, 1 table; accepted for publication in MNRAS; minor revisions from previous version, conclusions unchange

Investigating Overdensities around z > 6 Galaxies through ALMA Observations of [C II]

We present a search for companion [C II] emitters to known luminous sources at 6 < z < 6.5 in deep, archival ALMA observations. The observations are deep enough to detect sources with L_([CII])∼10⁸ at z ∼6. We identify three new robust line detections from a blind search of five deep fields centered on ultraluminous infrared galaxies and QSOs. We calculate the volume density of companions and find a relative overdensity of 6⁺⁴₋₃ and 86⁺⁶⁰₋₃₇ when comparing to current observational constraints and theoretical predictions, respectively. These results suggest that the central sources may be highly biased tracers of mass in the early universe. We find these companion lines to have comparable properties to other known galaxies at the same epoch. All companions lie less than 650 km s⁻¹ and between 25 and 60 kpc (projected) from their central source. To place these discoveries in context, we employ a mock galaxy catalog to estimate the luminosity function for [C II] during reionization and compare to our observations. The simulations support this result by showing a similar level of elevated counts found around such luminous [C II] sources

The relationship between galaxy and dark matter halo size from z ∼ 3 to the present

We explore empirical constraints on the statistical relationship between the radial size of galaxies and the radius of their host dark matter haloes from z similar to 0.1-3 using the Galaxy And Mass Assembly (GAMA) and Cosmic Assembly Near Infrared Deep Extragalactic Legacy Survey (CANDELS) surveys. We map dark matter halo mass to galaxy stellar mass using relationships from abundance matching, applied to the Bolshoi-Planck dissipationless N-body simulation. We define SRHR equivalent to r(e)/R-h as the ratio of galaxy radius to halo virial radius, and SRHR lambda equivalent to r(e)/(lambda R-h) as the ratio of galaxy radius to halo spin parameter times halo radius. At z similar to 0.1, we find an average value of SRHR similar or equal to 0.018 and SRHR. similar or equal to 0.5 with very little dependence on stellar mass. Stellar radius-halo radius (SRHR) and SRHR lambda have a weak dependence on cosmic time since z similar to 3. SRHR shows a mild decrease over cosmic time for low-mass galaxies, but increases slightly or does not evolve formoremassive galaxies. We find hints that at high redshift (z similar to 2-3), SRHR. is lower for more massive galaxies, while it shows no significant dependence on stellar mass at z less than or similar to 0.5. We find that for both the GAMA and CANDELS samples, at all redshifts from z similar to 0.1-3, the observed conditional size distribution in stellar mass bins is remarkably similar to the conditional distribution of lambda R-h. We discuss the physical interpretation and implications of these results

An Increasing Stellar Baryon Fraction in Bright Galaxies at High Redshift

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.Comment: Accepted to ApJ. 15 pages, 5 figures, 6 table

Reproducing the Stellar Mass/Halo Mass Relation in Simulated LCDM Galaxies: Theory vs Observational Estimates

We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses ranging from a few times 10^8 to nearly 10^12 solar masses. At z=0 our simulated galaxies have a baryon content and morphology typical of field galaxies. Over a stellar mass range of 2.2 x 10^3 to 4.5 x 10^10 solar masses, we find extremely good agreement between the SHM ratio in simulations and the present-day predictions from the statistical Abundance Matching Technique presented in Moster et al. (2012). This improvement over past simulations is due to a number systematic factors, each decreasing the SHM ratios: 1) gas outflows that reduce the overall SF efficiency but allow for the formation of a cold gas component 2) estimating the stellar masses of simulated galaxies using artificial observations and photometric techniques similar to those used in observations and 3) accounting for a systematic, up to 30 percent overestimate in total halo masses in DM-only simulations, due to the neglect of baryon loss over cosmic times. Our analysis suggests that stellar mass estimates based on photometric magnitudes can underestimate the contribution of old stellar populations to the total stellar mass, leading to stellar mass errors of up to 50 percent for individual galaxies. These results highlight the importance of using proper techniques to compare simulations with observations and reduce the perceived tension between the star formation efficiency in galaxy formation models and in real galaxies.Comment: Submitted to ApJ 9 pages, 5 figure

An Empirical Mass Function Distribution

The halo mass function, encoding the comoving number density of dark matter halos of a given mass, plays a key role in understanding the formation and evolution of galaxies. As such, it is a key goal of current and future deep optical surveys to constrain the mass function down to mass scales that typically host ${L}_{\star }$ galaxies. Motivated by the proven accuracy of Press–Schechter-type mass functions, we introduce a related but purely empirical form consistent with standard formulae to better than 4% in the medium-mass regime, {10}^{10}\mbox{--}{10}^{13}\,{h}^{-1}M☉. In particular, our form consists of four parameters, each of which has a simple interpretation, and can be directly related to parameters of the galaxy distribution, such as {L}_{\star }$. Using this form within a hierarchical Bayesian likelihood model, we show how individual mass-measurement errors can be successfully included in a typical analysis, while accounting for Eddington bias. We apply our form to a question of survey design in the context of a semi-realistic data model, illustrating how it can be used to obtain optimal balance between survey depth and angular coverage for constraints on mass function parameters. Open-source Python and R codes to apply our new form are provided at http://mrpy.readthedocs.org and https://cran.r-project.org/web/packages/tggd/index.html respectively

Improved Mock Galaxy Catalogs for the DEEP2 Galaxy Redshift Survey from Subhalo Abundance and Environment Matching

We develop empirical methods for modeling the galaxy population and populating cosmological N-body simulations with mock galaxies according to the observed properties of galaxies in survey data. We use these techniques to produce a new set of mock catalogs for the DEEP2 Galaxy Redshift Survey based on the output of the high-resolution Bolshoi simulation, as well as two other simulations with different cosmological parameters, all of which we release for public use. The mock-catalog creation technique uses subhalo abundance matching to assign galaxy luminosities to simulated dark-matter halos. It then adds color information to the resulting mock galaxies in a manner that depends on the local galaxy density, in order to reproduce the measured color-environment relation in the data. In the course of constructing the catalogs, we test various models for including scatter in the relation between halo mass and galaxy luminosity, within the abundance-matching framework. We find that there is no constant-scatter model that can simultaneously reproduce both the luminosity function and the autocorrelation function of DEEP2. This result has implications for galaxy-formation theory, and it restricts the range of contexts in which the mocks can be usefully applied. Nevertheless, careful comparisons show that our new mocks accurately reproduce a wide range of the other properties of the DEEP2 catalog, suggesting that they can be used to gain a detailed understanding of various selection effects in DEEP2.Comment: 24 pages, 13 figures, matches version accepted for publication in ApJS. Catalogs are available for download from the URL referenced in the Appendi