Search

We present the first results from the deep and wide 5 GHz radio observations of the Great Observatories Origins Deep Survey (GOODS)-North ($\sigma=3.5 \; \mu Jy \; beam^{-1}$, synthesized beam size $\theta =$ 1.47 arcsec $\times$ 1.42 arcsec, and 52 sources over 109 arcmin$^{2}$) and GOODS-South ($\sigma=3.0 \; \mu Jy \; beam^{-1}$, $\theta=$0.98 arcsec $\times$ 0.45 arcsec, and 88 sources over 190 arcmin$^{2}$) fields using the Karl G. Jansky Very Large Array. We derive radio spectral indices {\alpha} between 1.4 and 5 GHz using the beam-matched images and show that the overall spectral index distribution is broad even when the measured noise and flux bias are considered. We also find a clustering of faint radio sources around $\alpha=0.8$, but only within $S_{5GHz} < 150 \; \mu Jy$. We demonstrate that the correct radio spectral index is important for deriving accurate rest frame radio power and analyzing the radio-FIR correlation, and adopting a single value of $\alpha=0.8$ leads to a significant scatter and a strong bias in the analysis of the radio-FIR correlation, resulting from the broad and asymmetric spectral index distribution. When characterized by specific star formation rates, the starburst population (58%) dominates the 5 GHz radio source population, and the quiescent galaxy population (30%) follows a distinct trend in spectral index distribution and the radio-FIR correlation. Lastly, we offer suggestions on sensitivity and angular resolution for future ultra-deep surveys designed to trace the cosmic history of star formation and AGN activity using radio continuum as a probe., Comment: 25 pages, 13 figures, Accepted for the publication in the Astrophysical Journal

An understanding of the mass build-up in galaxies over time necessitates tracing the evolution of cold gas (molecular and atomic) in galaxies. To that end, we have conducted a pilot study called CO Observations with the LMT of the Blind Ultra-Deep H I Environment Survey (COOL BUDHIES). We have observed 23 galaxies in and around the two clusters Abell 2192 (z = 0.188) and Abell 963 (z = 0.206), where 12 are cluster members and 11 are slightly in the foreground or background, using about 28 total hours on the Redshift Search Receiver (RSR) on the Large Millimeter Telescope (LMT) to measure the $^{12}$CO J = 1 --> 0 emission line and obtain molecular gas masses. These new observations provide a unique opportunity to probe both the molecular and atomic components of galaxies as a function of environment beyond the local Universe. For our sample of 23 galaxies, nine have reliable detections (S/N$\geq$3.6) of the $^{12}$CO line, and another six have marginal detections (2.0 < S/N < 3.6). For the remaining eight targets we can place upper limits on molecular gas masses roughly between $10^9$ and $10^{10} M_\odot$. Comparing our results to other studies of molecular gas, we find that our sample is significantly more abundant in molecular gas overall, when compared to the stellar and the atomic gas component, and our median molecular gas fraction lies about $1\sigma$ above the upper limits of proposed redshift evolution in earlier studies. We discuss possible reasons for this discrepancy, with the most likely conclusion being target selection and Eddington bias., Comment: MNRAS, submitted

We explore star-formation histories (SFHs) of galaxies based on the evolution of the star-formation rate stellar mass relation (SFR-M*). Using data from the FourStar Galaxy Evolution Survey (ZFOURGE) in combination with far-IR imaging from the Spitzer and Herschel observatories we measure the SFR-M* relation at 0.5 < z < 4. Similar to recent works we find that the average infrared SEDs of galaxies are roughly consistent with a single infrared template across a broad range of redshifts and stellar masses, with evidence for only weak deviations. We find that the SFR-M* relation is not consistent with a single power-law of the form SFR ~ M*^a at any redshift; it has a power-law slope of a~1 at low masses, and becomes shallower above a turnover mass (M_0) that ranges from 10^9.5 - 10^10.8 Msol, with evidence that M_0 increases with redshift. We compare our measurements to results from state-of-the-art cosmological simulations, and find general agreement in the slope of the SFR-M* relation albeit with systematic offsets. We use the evolving SFR-M* sequence to generate SFHs, finding that typical SFRs of individual galaxies rise at early times and decline after reaching a peak. This peak occurs earlier for more massive galaxies. We integrate these SFHs to generate mass-growth histories and compare to the implied mass-growth from the evolution of the stellar mass function. We find that these two estimates are in broad qualitative agreement, but that there is room for improvement at a more detailed level. At early times the SFHs suggest mass-growth rates that are as much as 10x higher than inferred from the stellar mass function. However, at later times the SFHs under-predict the inferred evolution, as is expected in the case of additional growth due to mergers., Comment: Accepted for publication in the ApJ: 16 pages, 11 figures, 1 table

We present an analysis of the deepest Herschel images in four major extragalactic fields GOODS-North, GOODS-South, UDS and COSMOS obtained within the GOODS-Herschel and CANDELS-Herschel key programs. The picture provided by 10497 individual far-infrared detections is supplemented by the stacking analysis of a mass-complete sample of 62361 star-forming galaxies from the CANDELS-HST H band-selected catalogs and from two deep ground-based Ks band-selected catalogs in the GOODS-North and the COSMOS-wide fields, in order to obtain one of the most accurate and unbiased understanding to date of the stellar mass growth over the cosmic history. We show, for the first time, that stacking also provides a powerful tool to determine the dispersion of a physical correlation and describe our method called "scatter stacking" that may be easily generalized to other experiments. We demonstrate that galaxies of all masses from z=4 to 0 follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the logSFR-logM* relation with evidence for a flattening of the main sequence at high masses (log(M*/Msun) > 10.5) that becomes less prominent with increasing redshift and almost vanishes by z~2. This flattening may be due to the parallel stellar growth of quiescent bulges in star-forming galaxies. Within the main sequence, we measure a non varying SFR dispersion of 0.3 dex. The specific SFR (sSFR=SFR/M*) of star-forming galaxies is found to continuously increase from z=0 to 4. Finally we discuss the implications of our findings on the cosmic SFR history and show that more than 2/3 of present-day stars must have formed in a regime dominated by the main sequence mode. As a consequence we conclude that, although omnipresent in the distant Universe, galaxy mergers had little impact in shaping the global star formation history over the last 12.5 Gyr., Comment: 30 pages, 26 figures, accepted for publication in A&A

We use deep panchromatic datasets in the GOODS-N field, from GALEX to the deepest Herschel far-infrared and VLA radio continuum imaging, to explore, using mass-complete samples, the evolution of the star formation activity and dust attenuation of star-forming galaxies to z~4. Our main results can be summarized as follows: i) the slope of the SFR-M correlation is consistent with being constant, and equal to ~0.8 at least up to z~1.5, while its normalization keeps increasing with redshift; ii) for the first time here we are able to explore the FIR-radio correlation for a mass-selected sample of star-forming galaxies: the correlation does not evolve up to z~4; iii) we confirm that galaxy stellar mass is a robust proxy for UV dust attenuation in star-forming galaxies, with more massive galaxies being more dust attenuated, strikingly we find that this attenuation relation evolves very weakly with redshift, the amount of dust attenuation increasing by less than 0.3 magnitudes over the redshift range [0.5-4] for a fixed stellar mass, as opposed to a tenfold increase of star formation rate; iv) the correlation between dust attenuation and the UV spectral slope evolves in redshift, with the median UV spectral slope of star-forming galaxies becoming bluer with redshift. By z~3, typical UV slopes are inconsistent, given the measured dust attenuation, with the predictions of commonly used empirical laws. Finally, building on existing results, we show that gas reddening is marginally larger (by a factor of around 1.3) than stellar reddening at all redshifts probed, and also that the amount of dust attenuation at a fixed ISM metallicity increases with redshift. We speculate that our results support evolving ISM conditions of typical star-forming galaxies such that at z~1.5 Main Sequence galaxies have ISM conditions getting closer to those of local starbursts., Comment: 19 pages, 12 figure, accepted for publication in ApJ

We present deep IRAM Plateau de Bure Interferometer (PdBI) observations, searching for CO-emission toward two massive, non-lensed Lyman break galaxies (LBGs) at z=3.216 and 4.058. With one low significance CO detection (3.5 sigma) and one sensitive upper limit, we find that the CO lines are >~ 3-4 times weaker than expected based on the relation between IR and CO luminosities followed by similarly, massive galaxies at z=0-2.5. This is consistent with a scenario in which these galaxies have low metallicity, causing an increased CO-to-H_2 conversion factor, i.e., weaker CO-emission for a given molecular (H_2) mass. The required metallicities at z>3 are lower than predicted by the fundamental metallicity relation (FMR) at these redshifts, consistent with independent evidence. Unless our galaxies are atypical in this respect, detecting molecular gas in normal galaxies at z>3 may thus remain challenging even with ALMA., Comment: 9 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

We report the detection of 6.2um polycyclic aromatic hydrocarbon (PAH) and rest-frame 4-7um continuum emission in the z=4.055 submillimeter galaxy GN20, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. This represents the first detection of PAH emission at z>4. The strength of the PAH emission feature is consistent with a very high star formation rate of ~1600Msun/yr. We find that this intense starburst powers at least ~1/3 of the faint underlying 6um continuum emission, with an additional, significant (and perhaps dominant) contribution due to a power-law-like hot dust source, which we interpret to likely be a faint, dust-obscured active galactic nucleus (AGN). The inferred 6um AGN continuum luminosity is consistent with a sensitive upper limit on the hard X-ray emission as measured by the Chandra X-Ray Observatory if the previously undetected AGN is Compton-thick. This is in agreement with the finding at optical/infrared wavelengths that the galaxy and its nucleus are heavily dust-obscured. Despite the strong power-law component enhancing the mid-infrared continuum emission, the intense starburst associated with the photon-dominated regions that give rise to the PAH emission appears to dominate the total energy output in the infrared. GN20 is one of the most luminous starburst galaxies known at any redshift, embedded in a rich protocluster of star-forming galaxies. This investigation provides an improved understanding of the energy sources that power such exceptional systems, which represent the extreme end of massive galaxy formation at early cosmic times., Comment: 7 pages, 4 figures, to appear in ApJ (accepted March 10, 2014) [v3: title updated to be an exact match to the print version, optimized figure placement, one grammatical mistake fixed; no content changes relative to v2]

Measuring star formation rates (SFRs) in high-z galaxies with their rest-frame ultraviolet (UV) continuum can be uncertain because of dust obscuration. Prior studies had used the submillimeter emission at 850 um to determine the intrinsic SFRs of rest-frame UV selected galaxies, but the results suffered from the low sensitivity and poor resolution (~15''). Here, we use ultradeep Very Large Array 1.4 GHz images with ~1''-2'' resolutions to measure the intrinsic SFRs. We perform stacking analyses in the radio images centered on ~3500 Lyman break galaxies (LBGs) at z~4 in the Great Observatories Origins Deep Survey-North and South fields selected with Hubble Space Telescope/Advanced Camera for Surveys data. The stacked radio flux is very low, 0.08+/-0.15 uJy, implying a mean SFR of 6+/-11 M/yr. This is comparable to the uncorrected mean UV SFRs of 5 M/yr, implying that the z~4 LBGs have little dust extinction. The low SFR and dust extinction support the previous results that z~4 LBGs are in general not submillimeter galaxies. We further show that there is no statistically significant excess of dust-hidden star-forming components within ~22 kpc from the LBGs., Comment: 6 pages, 3 figures, and 2 tables. Accepted to the Astrophysical Journal (ApJ)

We report the detection of the CO J=1-0 emission line in three near-infrared selected star-forming galaxies at z~1.5 with the Very Large Array (VLA) and the Green Bank telescope (GBT). These observations directly trace the bulk of molecular gas in these galaxies. We find H_2 gas masses of 8.3 \pm 1.9 x 10^{10} M_sun, 5.6 \pm 1.4 x 10^{10} M_sun and 1.23 \pm 0.34 x 10^{11} M_sun for BzK-4171, BzK-21000 and BzK-16000, respectively, assuming a conversion alpha_CO=3.6 M_sun (K km s^{-1} pc^{2})^{-1}. We combined our observations with previous CO 2-1 detections of these galaxies to study the properties of their molecular gas. We find brightness temperature ratios between the CO 2-1 and CO 1-0 emission lines of 0.80_{-0.22}^{+0.35}, 1.22_{-0.36}^{+0.61} and 0.41_{-0.13}^{+0.23} for BzK-4171, BzK-21000 and BzK-16000, respectively. At the depth of our observations it is not possible to discern between thermodynamic equilibrium or sub-thermal excitation of the molecular gas at J=2. However, the low temperature ratio found for BzK-16000 suggests sub-thermal excitation of CO already at J=2. For BzK-21000, a Large Velocity Gradient model of its CO emission confirms previous results of the low-excitation of the molecular gas at J=3. From a stacked map of the CO 1-0 images, we measure a CO 2-1 to CO 1-0 brightness temperature ratio of 0.92_{-0.19}^{+0.28}. This suggests that, on average, the gas in these galaxies is thermalized up to J=2, has star-formation efficiencies of ~100 L_sun (K km s^{-1} pc^2)^{-1} and gas consumption timescales of ~0.4 Gyr, unlike SMGs and QSOs at high redshifts., Comment: Accepted for publication in ApJ

We describe deep, new, wide-field radio continuum observations of the Great Observatories Origins Deep Survey -- North (GOODS-N) field. The resulting map has a synthesized beamsize of ~1.7" and an r.m.s. noise level of ~3.9uJy/bm near its center and ~8uJy/bm at 15', from phase center. We have cataloged 1,230 discrete radio emitters, within a 40' x 40' region, above a 5-sigma detection threshold of ~20uJy at the field center. New techniques, pioneered by Owen & Morrison (2008), have enabled us to achieve a dynamic range of 6800:1 in a field that has significantly strong confusing sources. We compare the 1.4-GHz (20-cm) source counts with those from other published radio surveys. Our differential counts are nearly Euclidean below 100uJy with a median source diameter of ~1.2". This adds to the evidence presented by Owen & Morrison (2008) that the natural confusion limit may lie near ~1uJy. If the Euclidean slope of the counts continues down to the natural confusion limit as an extrapolation of our log N - log S, this indicates that the cutoff must be fairly sharp below 1uJy else the cosmic microwave background temperature would increase above 2.7K at 1.4 GHz., Comment: Accepted for publication in ApJS. 16 pages, 19 figures. Radio data and source list can be found at http://www.ifa.hawaii.edu/~morrison/GOODSN/

We present a multi-wavelength, UV-to-radio analysis for a sample of massive (M$_{\ast}$ $\sim$ 10$^{10}$ M$_\odot$) IRAC- and MIPS 24$\mu$m-detected Lyman Break Galaxies (LBGs) with spectroscopic redshifts z$\sim$3 in the GOODS-North field (L$_{\rm UV}$$>1.8\times$L$^{\ast}_{z=3}$). For LBGs without individual 24$\mu$m detections, we employ stacking techniques at 24$\mu$m, 1.1mm and 1.4GHz, to construct the average UV-to-radio spectral energy distribution and find it to be consistent with that of a Luminous Infrared Galaxy (LIRG) with L$\rm_{IR}$=4.5$^{+1.1}_{-2.3}$$\times 10^{11}$ L$_{\odot}$ and a specific star formation rate (SSFR) of 4.3 Gyr$^{-1}$ that corresponds to a mass doubling time $\sim$230 Myrs. On the other hand, when considering the 24$\mu$m-detected LBGs we find among them galaxies with L$\rm_{IR}> 10^{12}$ L$_{\odot}$, indicating that the space density of $z\sim$3 UV-selected Ultra-luminous Infrared Galaxies (ULIRGs) is $\sim$(1.5$\pm$0.5)$\times 10^{-5}$ Mpc$^{-3}$. We compare measurements of star formation rates (SFRs) from data at different wavelengths and find that there is tight correlation (Kendall's $\tau >$ 99.7%) and excellent agreement between the values derived from dust-corrected UV, mid-IR, mm and radio data for the whole range of L$\rm_{IR}$ up to L$\rm_{IR}$ $\sim$ 10$^{13}$ L$_{\odot}$. This range is greater than that for which the correlation is known to hold at z$\sim$2, possibly due to the lack of significant contribution from PAHs to the 24$\mu$m flux at $z\sim$3. The fact that this agreement is observed for galaxies with L$\rm_{IR}$ $>$ 10$^{12}$ L$_{\odot}$ suggests that star-formation in UV-selected ULIRGs, as well as the bulk of star-formation activity at this redshift, is not embedded in optically thick regions as seen in local ULIRGs and submillimeter-selected galaxies at $z=2$., Comment: Accepted for publication in ApJ

We present a study of a 20cm selected sample in the Deep SWIRE VLA Field, reaching a limiting flux density of ~13.5 uJy at the image center. In a 0.6x0.6 square degrees field, we are able to assign an optical/IR counterpart to 97% of the radio sources. Up to 11 passbands from the NUV to 4.5um are then used to sample the spectral energy distribution (SED) of these counterparts in order to investigate the nature of the host galaxies. By means of an SED template library and stellar population synthesis models we estimate photometric redshifts, stellar masses, and stellar population properties, dividing the sample in three sub-classes of quiescent, intermediate and star-forming galaxies. We focus on the radio sample in the redshift range 0.3

We present results from deep radio observations taken with the VLA at a center frequency of 1400 MHz cover a region of the SWIRE Spitzer Legacy survey, centered at 10 46 00, 59 01 00 (J2000). The reduction and cataloging of the radio sources are described. The survey presented is the deepest so far in terms of the radio source density on the sky. Perhaps surprisingly, the sources down to the bottom of the catalog have median angular sizes greater than 1 arcsecond, like their cousins 10-100 times stronger. If the log N - log S normalization remains constant at the lowest flux densities, there are about 6 sources per square arcminute down to 15 microJy at 20cm. Given the finite source sizes this implies we may reach the natural confusion limit near 1 microJy., Comment: 29 pages, 11 figures, accepted AJ

(Abridged) We report on the properties of radio-selected galaxies within 30 very-rich Abell clusters with z < 0.25. These radio data sample the ultra-faint (L(1.4) > 2E22 W/Hz) radio galaxy population with M_R < -21 using the well-known FIR/radio correlation to link the radio with ongoing star formation within individual cluster galaxies. These radio-selected galaxies reveal the `active' galaxy population (starburst and active galactic nuclei) within these rich cluster environments that can be identified regardless of their level of dust obscuration. For all clusters in our sample, we are sensitive to star formation rates (M > 5 M_sun) > 5 M_sun/yr. We have found that the excess number of low-luminosity `starburst' radio-selected galaxies (SBRG) found by Owen et. al. 1999 in Abell 2125 is not indicative of other rich clusters in our sample. The average fraction of SBRG is = 0.022+/-0.003. The A2125 fraction is f(SBRG) = 0.09+/-10.03 which is significantly different from the sample average at a >99.99% confidence level. The bimodal structure of both the x-ray brightness distribution and optical adaptively smoothed images of A1278 and A2125 suggests that ongoing cluster-cluster mergers may be enhancing this SBRG population. There is a significant spatial distribution difference between the low and high-luminosity (HLRG) radio-selected populations. The HLRGs seem to be a subclass of the cluster's massive red elliptical population, while the SBRGs have a projected radial distribution more like the blue spiral population. Our results indicate that most of the SBRGs are probably gas-rich disk galaxies undergoing > 5 M_sun/yr of star-formation., Comment: 18 pages, 6 figures, accepted by AJ

We report optical imaging and spectroscopy of 41 galaxies in a 22 arcmin square region surrounding Cygnus A. The results show that there is an extensive rich cluster associated with Cygnus A of Abell richness at least 1 and possibly as high as 4. The velocity histogram has two peaks, one centered on Cygnus A, and a more significant peak redshifted by about 2060 km/s from the velocity of Cygnus A. The dynamical centroid of the spatial distribution is also shifted somewhat to the NW. However, statistical tests show only weak evidence that there are two distinct clusters. The entire system has a velocity dispersion of 1581 km/s which is slightly larger than other, well studied, examples of rich clusters., Comment: 11 pages LaTeX, 1 table, 3 postscript figures. Accepted for publication in Ap.J. Letters