The remarkable Hubble Space Telescope(HST) data sets from the CANDELS, HUDF09, HUDF12, ERS, and BoRG/HIPPIES programs have allowed us to map the evolution of the rest-frame UV luminosity function (LF) fromz 10 toz 4. We develop new color criteria that more optimally utilize the full wavelength coverage from the optical, near-IR, and mid-IR observations over our search fields, while simultaneously minimizing the incompleteness and eliminating redshift gaps. We have identified 5859, 3001, 857, 481, 217, and 6 galaxy candidates atz 4,z 5,z 6,z 7,z 8, andz 10, respectively, from the ∼1000 arcmin 2 area covered by these data sets. This sample of >10,000 galaxy candidates at ⩾ z 4 is by far the largest assembled to date with HST. The selection ofz 4-8 candidates over the five CANDELS fields allows us to assess the cosmic variance; the largest variations are at ⩾ z 7. Our new LF determinations atz 4 andz 5 span a 6 mag baseline and reach to -16 AB mag. These determinations agree well with previous estimates, but the larger samples and volumes probed here result in a more reliable sampling of >L* galaxies and allow us to reassess the form of the UV LFs. Our new LF results strengthen our earlier findings to s 3.4 significance for a steeper faint-end slope of the UV LF at > z 4, with α evolving from a = - 1.64 0.04 atz 4 to a = - 2.06 0.13 atz 7 (and a = - 2.02 0.23 atz 8), consistent with that expected from the evolution of the halo mass function. We find less evolution in the characteristic magnitude M * fromz 7 toz 4; the observed evolution in the LF is now largely represented by changes in f*. No evidence for a non-Schechter-like form to the z ∼ 4-8 LFs is found. A simple conditional LF model based on halo growth and evolution in the M/L ratio µ +z ( ( 1) ) 1.5 of halos provides a good representation of the observed evolution.

This paper provides an update of our previous scaling relations (Genzel et al. 2015) between galaxy integrated molecular gas masses, stellar masses and star formation rates, in the framework of the star formation main-sequence (MS), with the main goal to test for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and ~1mm dust photometry, in a large sample of 1444 star forming galaxies (SFGs) between z=0 and 4. The sample covers the stellar mass range log(M*/M)=9.0-11.8, and star formation rates relative to that on the MS, δMS=SFR/SFR(MS), from 10 -1.3 to 10 2.2 . Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time tdepl, defined as the ratio of molecular gas mass to star formation rate, scales as (1+z) -0.6 × (δMS) -0.44 , and is only weakly dependent on stellar mass. The ratio of molecular-to-stellar mass μgas depends on (1+z) 2.5 × (δMS) 0.52 × (M*) -0.36 , which tracks the evolution of the specific star formation rate. The redshift dependence of μgas requires a curvature term, as may the mass-dependences of tdepl and μgas. We find no or only weak correlations of tdepl and μgas with optical size R or surface density once one removes the above scalings, but we caution that optical sizes may not be appropriate for the high gas and dust columns at high-z.

We combine molecular gas masses inferred from CO emission in 500 star forming galaxies (SFGs) between z=0 and 3, from the IRAM-COLDGASS, PHIBSS1/2 and other surveys, with gas masses derived from Herschel far-IR dust measurements in 512 galaxy stacks over the same stellar mass/redshift range. We constrain the scaling relations of molecular gas depletion time scale (t depl ) and gas to stellar mass ratio (M molgas /M * ) of SFGs near the star formation 'main-sequence' with redshift, specific star formation rate (sSFR) and stellar mass (M * ). The CO-and dust-based scaling relations agree remarkably well. This suggests that the CO  H 2 mass conversion factor varies little within ±0.6dex of the main sequence (sSFR(ms,z,M * )), and less than 0.3dex throughout this redshift range. This study builds on and strengthens the results of earlier work. We find that t depl scales as (1+z) -0.3  (sSFR/sSFR(ms,z,M * )) -0.5 , with little dependence on M * . The resulting steep redshift dependence of M molgas /M *  (1+z) 3 mirrors that of the sSFR and probably reflects the gas supply rate. The decreasing gas fractions at high M * are driven by the flattening of the SFR-M * relation. Throughout the redshift range probed a larger sSFR at constant M * is due to a combination of an increasing gas fraction and a decreasing depletion time scale. As a result galaxy integrated samples of the M molgas -SFR rate relation exhibit a super-linear slope, which increases with the range of sSFR. With these new relations it is now possible to determine M molgas with an accuracy of ±0.1dex in relative terms, and ±0.2dex including systematic uncertainties.

We extend the phenomenological study of the evolving galaxy population of Peng et al. (2010) to the central/satellite dichotomy in Yang et al. Sloan Digital Sky Survey (SDSS) groups. We find that satellite galaxies are responsible for all the environmental effects in our earlier work. The fraction of centrals that are red does not depend on their environment but only on their stellar masses, whereas that of the satellites depends on both. We define a relative satellite quenching efficiency ε sat , which is the fraction of blue centrals that are quenched upon becoming the satellite of another galaxy. This is shown to be independent of stellar mass, but to depend strongly on local overdensity, δ, ranging between 0.2 and at least 0.8. The red fraction of satellites correlate much better with the local overdensity δ, a measure of location within the group, than with the richness of the group, i.e., dark matter halo mass. This, and the fact that satellite quenching depends on local density and not on either the stellar mass of the galaxy or the dark matter halo mass, gives clues as to the nature of the satellite-quenching process. We furthermore show that the action of mass quenching on satellite galaxies is also independent of the dark matter mass of the parent halo. We then apply the Peng et al. approach to predict the mass functions of central and satellite galaxies, split into passive and active galaxies, and show that these match very well the observed mass functions from SDSS, further strengthening the validity of this phenomenological approach. We highlight the fact that the observed M * is exactly the same for the star-forming centrals and satellites and the observed M * for the star-forming satellites is independent of halo mass above 10 12 M , which emphasizes the universality of the mass-quenching process that we identified in Peng et al. Post-quenching merging modifies the mass function of the central galaxies but can increase the mass of typical centrals by only about 25%.

Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.