Optical Spectroscopy of Distant Red Galaxies

We present optical spectroscopic follow-up of a sample of Distant Red Galaxies (DRGs) with K 2.3, in the Hubble Deep Field South, the MS 1054-03 field, and the Chandra Deep Field South. Spectroscopic redshifts were obtained for 15 DRGs. Only 2 out of 15 DRGs are located at z < 2, suggesting a high efficiency to select high-redshift sources. From other spectroscopic surveys in the CDFS targeting intermediate to high redshift populations selected with different criteria, we find spectroscopic redshifts for a further 30 DRGs. We use the sample of spectroscopically confirmed DRGs to establish the high quality (scatter in \Delta z/(1+z) of ~ 0.05) of their photometric redshifts in the considered deep fields, as derived with EAZY (Brammer et al. 2008). Combining the spectroscopic and photometric redshifts, we find that 74% of DRGs with K 2. The combined spectroscopic and photometric sample is used to analyze the distinct intrinsic and observed properties of DRGs at z 2. In our photometric sample to K < 22.5, low-redshift DRGs are brighter in K than high-redshift DRGs by 0.7 mag, and more extincted by 1.2 mag in Av. Our analysis shows that the DRG criterion selects galaxies with different properties at different redshifts. Such biases can be largely avoided by selecting galaxies based on their rest-frame properties, which requires very good multi-band photometry and high quality photometric redshifts.Comment: Accepted for publication in the Astrophysical Journal, 13 pages, 8 figures, 5 table

FIREWORKS U38-to-24 micron photometry of the GOODS-CDFS: multi-wavelength catalog and total IR properties of distant Ks-selected galaxies

We present a Ks-selected catalog, dubbed FIREWORKS, for the Chandra Deep Field South (CDFS) containing photometry in U_38, B_435, B, V, V_606, R, i_775, I, z_850, J, H, Ks, [3.6 um], [4.5 um], [5.8 um], [8.0 um], and the MIPS [24 um] band. The imaging has a typical Ks limit of 24.3 mag (5 sigma, AB) and coverage over 113 arcmin^2 in all bands and 138 arcmin^2 in all bands but H. We cross-correlate our catalog with the 1 Ms X-ray catalog by Giacconi et al. (2002) and with all available spectroscopic redshifts to date. We find and explain systematic differences in a comparison with the 'z_850 + Ks'-selected GOODS-MUSIC catalog that covers ~90% of the field. We exploit the U38-to-24 micron photometry to determine which Ks-selected galaxies at 1.5<z<2.5 have the brightest total IR luminosities and which galaxies contribute most to the integrated total IR emission. The answer to both questions is that red galaxies are dominating in the IR. This is true no matter whether color is defined in the rest-frame UV, optical, or optical-to-NIR. We do find however that among the reddest galaxies in the rest-frame optical, there is a population of sources with only little mid-IR emission, suggesting a quiescent nature.Comment: Accepted for publication in the Astrophysical Journal, 20 pages, 10 figures, reference to website correcte

Lyman-alpha Emission from a Luminous z=8.68 Galaxy: Implications for Galaxies as Tracers of Cosmic Reionization

We report the discovery of Lyman-alpha emission (Ly$\alpha$) in the bright galaxy EGSY-2008532660 (hereafter EGSY8p7) using the MOSFIRE spectrograph at the Keck Observatory. First reported by Roberts-Borsani et al. (2015), it was selected for spectroscopic observations because of its photometric redshift ($z_{phot}=8.57^{+0.22}_{-0.43}$), apparent brightness (H$_{160}=25.26\pm0.09$) and red Spitzer/IRAC [3.6]-[4.5] color indicative of contamination by strong oxygen emission in the [4.5] band. With a total integration of $\sim$4.3 hours, our data reveal an emission line at $\simeq$11776 {\AA} which we argue is likely Ly$\alpha$ at a redshift $z_{spec}=8.683^{+0.001}_{-0.004}$, in good agreement with the photometric estimate. The line was detected independently on two nights using different slit orientations and its detection significance is $\sim7.5\sigma$. An overlapping skyline contributes significantly to the uncertainty on the total line flux although the significance of the detected line is robust to a variety of skyline-masking procedures. By direct addition and a Gaussian fit, we estimate a 95\% confidence range of 1.0--2.5$\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$, corresponding to a rest-frame equivalent width of 17--42 {\AA}. EGSY8p7 is the most distant galaxy confirmed spectroscopically to date, and the third luminous source in the EGS field beyond $z_{phot}\gtrsim7.5$ with detectable Ly$\alpha$ emission viewed at a time when the intergalactic medium is believed to be fairly neutral. Although the reionization process was probably patchy, we discuss whether luminous sources with prominent IRAC color excesses may harbor harder ionizing spectra than the dominant fainter population thereby creating earlier ionized bubbles. Further spectroscopic follow-up of such bright sources promises important insight into the early formation of galaxies.Comment: V3: ApJL accepted; 7 pages, 4 figures, 1 tabl

A Public Ks-selected Catalog in the COSMOS/UltraVISTA Field: Photometry, Photometric Redshifts and Stellar Population Parameters

We present a catalog covering 1.62 deg^2 of the COSMOS/UltraVISTA field with PSF-matched photometry in 30 photometric bands. The catalog covers the wavelength range 0.15um - 24um including the available GALEX, Subaru, CFHT, VISTA and Spitzer data. Catalog sources have been selected from the DR1 UltraVISTA Ks band imaging that reaches a depth of K_{s,tot} = 23.4 AB (90% completeness). The PSF-matched catalog is generated using position-dependent PSFs ensuring accurate colors across the entire field. Also included is a catalog of photometric redshifts (z_phot) for all galaxies computed with the EAZY code. Comparison with spectroscopy from the zCOSMOS 10k bright sample shows that up to z ~ 1.5 the z_phot are accurate to dz/(1 + z) = 0.013, with a catastrophic outlier fraction of only 1.6%. The z_phot also show good agreement with the z_phot from the NEWFIRM Medium Band Survey (NMBS) out to z ~ 3. A catalog of stellar masses and stellar population parameters for galaxies determined using the FAST spectral energy distribution fitting code is provided for all galaxies. Also included are rest-frame U-V and V-J colors, L_2800 and L_IR. The UVJ color-color diagram confirms that the galaxy bi-modality is well-established out to z ~ 2. Star-forming galaxies also obey a star forming "main sequence" out to z ~ 2.5, and this sequence evolves in a manner consistent with previous measurements. The COSMOS/UltraVISTA Ks-selected catalog covers a unique parameter space in both depth, area, and multi-wavelength coverage and promises to be a useful tool for studying the growth of the galaxy population out to z ~ 3 - 4.Comment: 20 pages, 14 figures. Accepted to the ApJSS. Catalog data products available for download here: http://www.strw.leidenuniv.nl/galaxyevolution/ULTRAVISTA

Stellar mass functions of galaxies at 4<z<7 from an IRAC-selected sample in COSMOS/UltraVISTA: limits on the abundance of very massive galaxies

We build a Spitzer IRAC complete catalog of objects, obtained by complementing the $K_\mathrm{s}$-band selected UltraVISTA catalog with objects detected in IRAC only. With the aim of identifying massive (i.e., $\log(M_*/M_\odot)>11$) galaxies at $4<z<7$, we consider the systematic effects on the measured photometric redshifts from the introduction of an old and dusty SED template and from the introduction of a bayesian prior taking into account the brightness of the objects, as well as the systematic effects from different star formation histories (SFHs) and from nebular emission lines in the recovery of stellar population parameters. We show that our results are most affected by the bayesian luminosity prior, while nebular emission lines and SFHs only introduce a small dispersion in the measurements. Specifically, the number of $4<z<7$ galaxies ranges from 52 to 382 depending on the adopted configuration. Using these results we investigate, for the first time, the evolution of the massive end of the stellar mass functions (SMFs) at $4<z<7$. Given the rarity of very massive galaxies in the early universe, major contributions to the total error budget come from cosmic variance and poisson noise. The SMF obtained without the introduction of the bayesian luminosity prior does not show any evolution from $z\sim6.5$ to $z\sim 3.5$, implying that massive galaxies could already be present when the Universe was $\sim0.9$~Gyr old. However, the introduction of the bayesian luminosity prior reduces the number of $z>4$ galaxies with best fit masses $\log(M_*/M_\odot)>11$ by 83%, implying a rapid growth of very massive galaxies in the first 1.5 Gyr of cosmic history. From the stellar-mass complete sample, we identify one candidate of a very massive ($\log(M_*/M_\odot)\sim11.5$), quiescent galaxy at $z\sim5.4$, with MIPS $24\mu$m detection suggesting the presence of a powerful obscured AGN.Comment: 23 pages, 18 figures. ApJ accepte

High-precision Photometric Redshifts from Spitzer/IRAC: Extreme [3.6]-[4.5] Colors Identify Galaxies in the Redshift Range z~6.6-6.9

One of the most challenging aspects of studying galaxies in the z>~7 universe is the infrequent confirmation of their redshifts through spectroscopy, a phenomenon thought to occur from the increasing opacity of the intergalactic medium to Lya photons at z>6.5. The resulting redshift uncertainties inhibit the efficient search for [C II] in z~7 galaxies with sub-mm instruments such as ALMA, given their limited scan speed for faint lines. One means by which to improve the precision of the inferred redshifts is to exploit the potential impact of strong nebular emission lines on the colors of z~4-8 galaxies as observed by Spitzer/IRAC. At z~6.8, galaxies exhibit IRAC colors as blue as [3.6]-[4.5] ~-1, likely due to the contribution of [O III]+Hb to the 3.6 mum flux combined with the absence of line contamination in the 4.5 mum band. In this paper we explore the use of extremely blue [3.6]-[4.5] colors to identify galaxies in the narrow redshift window z~6.6-6.9. When combined with an I-dropout criterion, we demonstrate that we can plausibly select a relatively clean sample of z~6.8 galaxies. Through a systematic application of this selection technique to our catalogs from all five CANDELS fields, we identify 20 probable z~6.6-6.9 galaxies. We estimate that our criteria select the ~50% strongest line emitters at z~6.8 and from the IRAC colors we estimate a typical [O III]+Hb rest-frame equivalent width of 1085A for this sample. The small redshift uncertainties on our sample make it particularly well suited for follow-up studies with facilities such as ALMA.Comment: In submission to the Astrophysical Journal, updated in response to the referee report, 13 pages, 11 figures, 1 tabl

ZFIRE: Using Hα\alpha equivalent widths to investigate the in situ initial mass function at z~2

We use the ZFIRE survey (http://zfire.swinburne.edu.au) to investigate the high mass slope of the initial mass function (IMF) for a mass-complete (log10(M$_*$/M$_\odot$)~9.3) sample of 102 star-forming galaxies at z~2 using their H$\alpha$ equivalent widths (H$\alpha$-EW) and rest-frame optical colours. We compare dust-corrected H$\alpha$-EW distributions with predictions of star-formation histories (SFH) from PEGASE.2 and Starburst99 synthetic stellar population models. We find an excess of high H$\alpha$-EW galaxies that are up to 0.3--0.5 dex above the model-predicted Salpeter IMF locus and the H$\alpha$-EW distribution is much broader (10--500 \AA) than can easily be explained by a simple monotonic SFH with a standard Salpeter-slope IMF. Though this discrepancy is somewhat alleviated when it is assumed that there is no relative attenuation difference between stars and nebular lines, the result is robust against observational biases, and no single IMF (i.e. non-Salpeter slope) can reproduce the data. We show using both spectral stacking and Monte Carlo simulations that starbursts cannot explain the EW distribution. We investigate other physical mechanisms including models with variations in stellar rotation, binary star evolution, metallicity, and the IMF upper-mass cutoff. IMF variations and/or highly rotating extreme metal poor stars (Z~0.1Z$_\odot$) with binary interactions are the most plausible explanations for our data. If the IMF varies, then the highest H$\alpha$-EWs would require very shallow slopes ($\Gamma$>-1.0) with no one slope able to reproduce the data. Thus, the IMF would have to vary stochastically. We conclude that the stellar populations at z~2 show distinct differences from local populations and there is no simple physical model to explain the large variation in H$\alpha$-EWs at z~2.Comment: Accepted to MNRAS. 43 pages, 27 Figures. Survey website: http://zfire.swinburne.edu.au

Discovery of a dark, massive, ALMA-only galaxy at z~5-6 in a tiny 3-millimeter survey

We report the serendipitous detection of two 3 mm continuum sources found in deep ALMA Band 3 observations to study intermediate redshift galaxies in the COSMOS field. One is near a foreground galaxy at 1.3", but is a previously unknown dust-obscured star-forming galaxy (DSFG) at probable $z_{CO}=3.329$, illustrating the risk of misidentifying shorter wavelength counterparts. The optical-to-mm spectral energy distribution (SED) favors a grey $\lambda^{-0.4}$ attenuation curve and results in significantly larger stellar mass and SFR compared to a Calzetti starburst law, suggesting caution when relating progenitors and descendants based on these quantities. The other source is missing from all previous optical/near-infrared/sub-mm/radio catalogs ("ALMA-only"), and remains undetected even in stacked ultradeep optical ($>29.6$ AB) and near-infrared ($>27.9$ AB) images. Using the ALMA position as a prior reveals faint $SNR\sim3$ measurements in stacked IRAC 3.6+4.5, ultradeep SCUBA2 850$\mu$m, and VLA 3GHz, indicating the source is real. The SED is robustly reproduced by a massive $M^*=10^{10.8}$M$_\odot$ and $M_{gas}=10^{11}$M$_\odot$, highly obscured $A_V\sim4$, star forming $SFR\sim300$ M$_{\odot}$yr$^{-1}$ galaxy at redshift $z=5.5\pm$1.1. The ultrasmall 8 arcmin$^{2}$ survey area implies a large yet uncertain contribution to the cosmic star formation rate density CSFRD(z=5) $\sim0.9\times10^{-2}$ M$_{\odot}$ yr$^{-1}$ Mpc$^{-3}$, comparable to all ultraviolet-selected galaxies combined. These results indicate the existence of a prominent population of DSFGs at $z>4$, below the typical detection limit of bright galaxies found in single-dish sub-mm surveys, but with larger space densities $\sim3 \times 10^{-5}$ Mpc$^{-3}$, higher duty cycles $50-100\%$, contributing more to the CSFRD, and potentially dominating the high-mass galaxy stellar mass function.Comment: Accepted for publication in ApJ. 2 galaxies, too many pages, 8 figures, 2 table