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.

We present deep $J$ and $H$-band images in the extended Great Observatories Origins Deep Survey-North (GOODS-N) field covering an area of 0.22 $\rm{deg}^{2}$. The observations were taken using WIRCam on the 3.6-m Canada France Hawaii Telescope (CFHT). Together with the reprocessed $K_{\rm s}$-band image, the $5\sigma$ limiting AB magnitudes (in 2" diameter apertures) are 24.7, 24.2, and 24.4 AB mag in the $J$, $H$, and $K_{\rm s}$ bands, respectively. We also release a multi-band photometry and photometric redshift catalog containing 93598 sources. For non-X-ray sources, we obtained a photometric redshift accuracy $\sigma_{\mathrm{NMAD}}=0.036$ with an outlier fraction $\eta = 7.3\%$. For X-ray sources, which are mainly active galactic nuclei (AGNs), we cross-matched our catalog with the updated 2M-CDFN X-ray catalog from Xue et al. (2016) and found that 658 out of 683 X-ray sources have counterparts. $GALEX$ UV data are included in the photometric redshift computation for the X-ray sources to give $\sigma_{\mathrm{NMAD}} = 0.040$ with $\eta=10.5\%$. Our approach yields more accurate photometric redshift estimates compared to previous works in this field. In particular, by adopting AGN-galaxy hybrid templates, our approach delivers photometric redshifts for the X-ray counterparts with fewer outliers compared to the 3D-HST catalog, which fit these sources with galaxy-only templates.

We have produced for the first time a detailed velocity map of the giant filamentary nebula surrounding NGC 1275, the Perseus cluster's brightest galaxy, and revealed a previously unknown rich velocity structure across the entire nebula. We present new observations of the low-velocity component of this nebula with the optical imaging Fourier transform spectrometer SITELLE at CFHT. With its wide field of view ($\sim$11'$\times$11'), SITELLE is the only integral field unit spectroscopy instrument able to cover the 80 kpc$\times$55 kpc (3.8'$\times$2.6') large nebula in NGC 1275. Our analysis of these observations shows a smooth radial gradient of the [N II]$\lambda$6583/$\text{H} \alpha$ line ratio, suggesting a change in the ionization mechanism and source across the nebula, while the dispersion profile shows a general decrease with increasing distance from the AGN at up to $\sim 10$ kpc. The velocity map shows no visible general trend or rotation, indicating that filaments are not falling uniformly onto the galaxy, nor being pulled out from it. Comparison between the physical properties of the filaments and Hitomi measurements of the X-ray gas dynamics in Perseus are also explored.

We investigate the relationship between the black hole accretion rate (BHAR) and star-formation rate (SFR) for Milky Way (MW) and Andromeda (M31)-mass progenitors from z = 0.2 - 2.5. We source galaxies from the Ks-band selected ZFOURGE survey, which includes multi-wavelenth data spanning 0.3 - 160um. We use decomposition software to split the observed SEDs of our galaxies into their active galactic nuclei (AGN) and star-forming components, which allows us to estimate BHARs and SFRs from the infrared (IR). We perform tests to check the robustness of these estimates, including a comparison to BHARs and SFRs derived from X-ray stacking and far-IR analysis, respectively. We find as the progenit- ors evolve, their relative black hole-galaxy growth (i.e. their BHAR/SFR ratio) increases from low to high redshift. The MW-mass progenitors exhibit a log-log slope of 0.64 +/- 0.11, while the M31-mass progenitors are 0.39 +/- 0.08. This result contrasts with previous studies that find an almost flat slope when adopting X-ray/AGN-selected or mass-limited samples and is likely due to their use of a broad mixture of galaxies with different evolutionary histories. Our use of progenitor-matched samples highlights the potential importance of carefully selecting progenitors when searching for evolutionary relationships between BHAR/SFRs. Additionally, our finding that BHAR/SFR ratios do not track the rate at which progenitors quench casts doubts over the idea that the suppression of star-formation is predominantly driven by luminous AGN feedback (i.e. high BHARs).

We present observations collected in the CFHTLS-VIPERS region in the ultraviolet (UV) with the GALEX satellite (far and near UV channels) and the near infrared with the CFHT/WIRCam camera ($K_s$-band) over an area of 22 and 27 deg$^2$, respectively. The depth of the photometry was optimized to measure the physical properties (e.g., SFR, stellar masses) of all the galaxies in the VIPERS spectroscopic survey. The large volume explored by VIPERS will enable a unique investigation of the relationship between the galaxy properties and their environment (density field and cosmic web) at high redshift (0.5 < z < 1.2). In this paper, we present the observations, the data reductions and the build-up of the multi-color catalogues. The CFHTLS-T0007 (gri-{\chi}^2) images are used as reference to detect and measure the $K_s$-band photometry, while the T0007 u-selected sources are used as priors to perform the GALEX photometry based on a dedicated software (EMphot). Our final sample reaches $NUV_{AB}$~25 (at 5{\sigma}) and $K_{AB}$~22 (at 3{\sigma}). The large spectroscopic sample (~51,000 spectroscopic redshifts) allows us to highlight the robustness of our star/galaxy separation, and the reliability of our photometric redshifts with a typical accuracy $\sigma_z \le$ 0.04 and a catastrophic failure rate {\eta} < 2% down to i~23. We present various tests on the $K_s$ band completeness and photometric redshift accuracy by comparing with existing, overlapping deep photometric catalogues. Finally, we discuss the BzK sample of passive and active galaxies at high redshift and the evolution of galaxy morphology in the (NUV-r) vs (r-K_s) diagram at low redshift (z < 0.25) thanks to the high image quality of the CFHTLS. The images, catalogues and photometric redshifts for 1.5 million sources (down to $NUV \le$ 25 or $K_s \le$ 22) are released and available at this URL: http://cesam.lam.fr/vipers-mls/

We investigate the evolution of the galaxy stellar mass function (SMF) and stellar mass density from redshift z=0.2 to z=1.5 of a $K_{AB}$<22-selected sample with highly reliable photometric redshifts and over an unprecedentedly large area. Our study is based on NIR observations carried out with WIRCam at CFHT over the footprint of the VIPERS spectroscopic survey and benefits from the high quality optical photometry from the CFHTLS and UV observations with the GALEX satellite. The accuracy of our photometric redshifts is $\sigma_z$ < 0.03 and 0.05 for the bright ($i_{AB}$<22.5) and faint ($i_{AB}$>22.5) samples, respectively. The SMF is measured with ~760,000 galaxies down to $K_s$=22 and over an effective area of ~22.4 deg$^2$, the latter of which drastically reduces the statistical uncertainties (i.e. Poissonian error & cosmic variance). We point out the importance of a careful control of the photometric calibration, whose impact becomes quickly dominant when statistical uncertainties are reduced, which will be a major issue for future generation of cosmological surveys with, e.g. EUCLID or LSST. By exploring the rest-frame (NUV-r) vs (r-$K_s$) color-color diagram separating star-forming and quiescent galaxies, (1) we find that the density of very massive log($M_*/ M_{\odot}$) > 11.5 galaxies is largely dominated by quiescent galaxies and increases by a factor 2 from z~1 to z~0.2, which allows for additional mass assembly via dry mergers, (2) we confirm a scenario where star formation activity is impeded above a stellar mass log($M^*_{SF} / M_{\odot}$) = 10.64$\pm$0.01, a value that is found to be very stable at 0.2 < z < 1.5, (3) we discuss the existence of a main quenching channel that is followed by massive star-forming galaxies, and finally (4) we characterise another quenching mechanism required to explain the clear excess of low-mass quiescent galaxies observed at low redshift.

We build a set of composite galaxy SEDs by de-redshifting and scaling multi-wavelength photometry from galaxies in the ZFOURGE survey, covering the CDFS, COSMOS, and UDS fields. From a sample of ~4000 K_s-band selected galaxies, we define 38 composite galaxy SEDs that yield continuous low-resolution spectra (R~45) over the rest-frame range 0.1-4 um. Additionally, we include far infrared photometry from the Spitzer Space Telescope and the Herschel Space Observatory to characterize the infrared properties of our diverse set of composite SEDs. From these composite SEDs we analyze the rest-frame UVJ colors, as well as the ratio of IR to UV light (IRX) and the UV slope ($\beta$) in the IRX$-\beta$ dust relation at 1<z<3. Blue star-forming composite SEDs show IRX and $\beta$ values consistent with local relations; dusty star-forming galaxies have considerable scatter, as found for local IR bright sources, but on average appear bluer than expected for their IR fluxes. We measure a tight linear relation between rest-frame UVJ colors and dust attenuation for star-forming composites, providing a direct method for estimating dust content from either (U-V) or (V-J) rest-frame colors for star-forming galaxies at intermediate redshifts.

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.

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.

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.

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.

We report the discovery of a grand-design spiral galaxy associated with a double-lobed radio source. J1649+2635 (z = 0.0545) is a red spiral galaxy with a prominent bulge that it is associated with a L$_{1.4{\rm GHz}}\sim$10$^{24}$WHz$^{-1}$ double-lobed radio source that spans almost 100kpc. J1649+2635 has a black hole mass of M$_{\rm BH} \sim$ 3--7 $\times$ 10$^8$M$_{\odot}$ and SFR$\sim$ 0.26 -- 2.6M$_{\odot}$year$^{-1}$. The galaxy hosts a $\sim$96kpc diffuse optical halo, which is unprecedented for spiral galaxies. We find that J1649+2635 resides in an overdense environment with a mass of M$_{dyn} = 7.7^{+7.9}_{-4.3} \times 10^{13}$M$_{\odot}$, likely a galaxy group below the detection threshold of the ROSAT All-Sky Survey. We suggest one possible scenario for the association of double-lobed radio emission from J1649+2635 is that the source may be similar to a Seyfert galaxy, located in a denser-than-normal environment. The study of spiral galaxies that host large-scale radio emission is important because although rare in the local Universe, these sources may be more common at high-redshifts.

Galaxies with stellar masses near M* contain the majority of stellar mass in the universe, and are therefore of special interest in the study of galaxy evolution. The Milky Way (MW) and Andromeda (M31) have present day stellar masses near M*, at 5x10^10 Msol (MW-mass) and 10^11 Msol (M31-mass). We study the typical progenitors of these galaxies using ZFOURGE, a deep medium-band near-IR imaging survey, which is sensitive to the progenitors of these galaxies out to z~3. We use abundance-matching techniques to identify the main progenitors of these galaxies at higher redshifts. We measure the evolution in the stellar mass, rest-frame colors, morphologies, far-IR luminosities, and star-formation rates combining our deep multiwavelength imaging with near-IR HST imaging from CANDELS, and far-IR imaging from GOODS-H and CANDELS-H. The typical MW-mass and M31-mass progenitors passed through the same evolution stages, evolving from blue, star-forming disk galaxies at the earliest stages, to redder dust-obscured IR-luminous galaxies in intermediate stages, and to red, more quiescent galaxies at their latest stages. The progenitors of the MW-mass galaxies reached each evolutionary stage at later times (lower redshifts) and with stellar masses that are a factor of 2-3 lower than the progenitors of the M31-mass galaxies. The process driving this evolution, including the suppression of star-formation in present-day M* galaxies requires an evolving stellar-mass/halo-mass ratio and/or evolving halo-mass threshold for quiescent galaxies. The effective size and star-formation rates imply that the baryonic cold-gas fractions drop as galaxies evolve from high redshift to z~0 and are strongly anticorrelated with an increase in the S\'ersic index. Therefore, the growth of galaxy bulges in M* galaxies corresponds to a rapid decline in the galaxy gas fractions and/or a decrease in the star-formation efficiency.

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.

We seek to derive star formation rates (SFR) and stellar masses (M_star) in distant galaxies and to quantify the main uncertainties affecting their measurement. We explore the impact of the assumptions made in their derivation with standard calibrations or through a fitting process, as well as the impact of the available data, focusing on the role of IR emission originating from dust. We build a sample of galaxies with z>1, all observed from the UV to the IR (rest frame). The data are fitted with the code CIGALE, which is also used to build and analyse a catalogue of mock galaxies. Models with different SFHs are introduced. We define different set of data, with or without a good sampling of the UV range, NIR, and thermal IR data. The impact of these different cases on the determination of M_star and SFR are analysed. Exponentially decreasing models with a redshift formation of the stellar population z ~8 cannot fit the data correctly. The other models fit the data correctly at the price of unrealistically young ages when the age of the single stellar population is taken to be a free parameter. The best fits are obtained with two stellar populations. As long as one measurement of the dust emission continuum is available, SFR are robustly estimated whatever the chosen model is, including standard recipes. M_star measurement is more subject to uncertainty, depending on the chosen model and the presence of NIR data, with an impact on the SFR-M_star scatter plot. Conversely, when thermal IR data from dust emission are missing, the uncertainty on SFR measurements largely exceeds that of stellar mass. Among all physical properties investigated here, the stellar ages are found to be the most difficult to constrain and this uncertainty acts as a second parameter in SFR measurements and as the most important parameter for M_star measurements.

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.

We investigate the relation between total X-ray emission from star-forming galaxies and their star formation activity. Using nearby late-type galaxies and ULIRGs from Paper I and star-forming galaxies from Chandra Deep Fields, we construct a sample of 66 galaxies spanning the redshift range z~0-1.3 and the star-formation rate (SFR) range ~0.1-10^3 M_sun/yr. In agreement with previous results, we find that the Lx-SFR relation is consistent with a linear law both at z=0 and for the z=0.1-1.3 CDF galaxies, within the statistical accuracy of ~0.1 in the slope of the Lx-SFR relation. For the total sample, we find a linear scaling relation Lx/SFR~(4.0\pm 0.4)x10^{39}(erg/s)/(Msun/yr), with a scatter of ~0.4 dex. About ~2/3 of the 0.5-8 keV luminosity generated per unit SFR is expected to be due to HMXBs. We find no statistically significant trends in the mean Lx/SFR ratio with the redshift or star formation rate and constrain the amplitude of its variations by <0.1-0.2 dex. These properties make X-ray observations a powerful tool to measure the star formation rate in normal star-forming galaxies that dominate the source counts at faint fluxes.

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.

We study the mid- to far-IR properties of a 24um-selected flux-limited sample (S24 > 5mJy) of 154 intermediate redshift (<z>~0.15), infrared luminous galaxies, drawn from the 5MUSES survey. By combining existing mid-IR spectroscopy and new Herschel SPIRE submm photometry from the HerMES program, we derived robust total infrared luminosity (LIR) and dust mass (Md) estimates and infered the relative contribution of the AGN to the infrared energy budget of the sources. We found that the total infrared emission of galaxies with weak 6.2um PAH emission (EW<0.2um) is dominated by AGN activity, while for galaxies with EW>0.2um more than 50% of the LIR arises from star formation. We also found that for galaxies detected in the 250-500um Herschel bands an AGN has a statistically insignificant effect on the temperature of the cold dust and the far-IR colours of the host galaxy, which are primarily shaped by star formation activity. For star-forming galaxies we reveal an anti-correlation between the LIR-to-rest-frame 8um luminosity ratio, IR8 = LIR\L8, and the strength of PAH features. We found that this anti-correlation is primarily driven by variations in the PAHs emission, and not by variations in the 5-15um mid-IR continuum emission. Using the [NeIII]/[NeII] line ratio as a tracer of the hardness of the radiation field, we confirm that galaxies with harder radiation fields tend to exhibit weaker PAH features, and found that they have higher IR8 values and higher dust-mass-weighted luminosities (LIR/Md), the latter being a proxy for the dust temperature (Td). We argue that these trends originate either from variations in the environment of the star-forming regions or are caused by variations in the age of the starburst. Finally, we provide scaling relations that will allow estimating LIR, based on single-band observations with the mid-infrared instrument, on board the upcoming JWST.

We present here a new spectral energy distribution (SED) fitting approach that we adopt to select radio-excess sources amongst distant star-forming galaxies in the GOODS-Herschel (North) field and to reveal the presence of hidden, highly obscured AGN. Through extensive SED analysis of 458 galaxies with radio 1.4 GHz and mid-IR 24 um detections using some of the deepest Chandra X-ray, Spitzer and Herschel infrared, and VLA radio data available to date, we have robustly identified a sample of 51 radio-excess AGN (~1300 deg^-2) out to redshift z~3. These radio-excess AGN have a significantly lower far-IR/radio ratio (q<1.68) than the typical relation observed for star-forming galaxies (q~2.2). We find that ~45% of these radio-excess sources have a dominant AGN component in the mid-IR band, while for the remainders the excess radio emission is the only indicator of AGN activity. The fraction of radio-excess AGN increases with X-ray luminosity reaching ~60% at Lx~10^44-10^45 erg/s, making these sources an important part of the total AGN population. However, almost half (24/51) of these radio-excess AGN are not detected in the deep Chandra X-ray data, suggesting that some of these sources might be heavily obscured. We also find that the specific star formation rates (sSFRs) of the radio-excess AGN are on average lower that those observed for X-ray selected AGN hosts, indicating that our sources are forming stars more slowly than typical AGN hosts, and possibly their star formation is progressively quenching.

We explore the effects of active galactic nuclei (AGN) and star formation activity on the infrared (0.3-1000 microns) spectral energy distributions of luminous infrared galaxies from z = 0.5 to 4.0. We have compiled a large sample of 151 galaxies selected at 24 microns (S24 > 100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-IR spectrum into contributions from star formation and AGN activity. A significant portion (~25%) of our sample is dominated by an AGN in the mid-IR. Based on the mid-IR classification, we divide our full sample into four sub-samples: z~1 star-forming (SF) sources; z~2 SF sources; AGN with clear 9.7 micron silicate absorption; and AGN with featureless mid-IR spectra. From our large spectroscopic sample and wealth of multi-wavelength data, including deep Herschel imaging at 100, 160, 250, 350, and 500 microns, we use 95 galaxies with complete spectral coverage to create a composite spectral energy distribution (SED) for each sub-sample. We then fit a two-temperature component modified blackbody to the SEDs. We find that the IR SEDs have similar cold dust temperatures, regardless of the mid-IR power source, but display a marked difference in the warmer dust temperatures. We calculate the average effective temperature of the dust in each sub-sample and find a significant (~20 K) difference between the SF and AGN systems. We compare our composite SEDs to local templates and find that local templates do not accurately reproduce the mid-IR features and dust temperatures of our high redshift systems. High redshift IR luminous galaxies contain significantly more cool dust than their local counterparts. We find that a full suite of photometry spanning the IR peak is necessary to accurately account for the dominant dust temperature components in high redshift IR luminous galaxies.

Dusty galaxies at z ~ 2 span a wide range of relative brightness between rest-frame mid-infrared (8um) and ultraviolet wavelengths. We attempt to determine the physical mechanism responsible for this diversity. Dust-obscured galaxies (DOGs), which have rest-frame mid-IR to UV flux density ratios > 1000, might be abnormally bright in the mid-IR, perhaps due to prominent AGN and/or PAH emission, or abnormally faint in the UV. We use far-infrared data from the GOODS-Herschel survey to show that most DOGs with 10^12 L_Sun < L_IR < 10^13 L_Sun are not abnormally bright in the mid-IR when compared to other dusty galaxies with similar IR (8--1000um) luminosities. We observe a relation between the median IR to UV luminosity ratios and the median UV continuum power-law indices for these galaxies, and we find that only 24% have specific star formation rates which indicate the dominance of compact star-forming regions. This circumstantial evidence supports the idea that the UV- and IR-emitting regions in these galaxies are spatially coincident, which implies a connection between the abnormal UV faintness of DOGs and dust obscuration. We conclude that the range in rest-frame mid-IR to UV flux density ratios spanned by dusty galaxies at z ~ 2 is due to differing amounts of UV obscuration. Of galaxies with these IR luminosities, DOGs are the most obscured. We attribute differences in UV obscuration to either: 1) differences in the degree of alignment between the spatial distributions of dust and massive stars, or 2) differences in the total dust content.

Using deep 100-160 micron observations in GOODS-S from the GOODS-H survey, combined with HST/WFC3 NIR imaging from CANDELS, we present the first morphological analysis of a complete, FIR selected sample of 52 ULIRGs at z~2. We also make use of a comparison sample of galaxies without Herschel detections but with the same z and magnitude distribution. Our visual classifications of these two samples indicate that the fraction of objects with disk and spheroid morphologies is roughly the same but that there are significantly more mergers, interactions, and irregular galaxies among the ULIRGs. The combination of disk and irregular/interacting morphologies suggests that early stage interactions and minor mergers could play an important role in ULIRGs at z~2. We compare these fractions with those of a z~1 sample across a wide luminosity range and find that the fraction of disks decreases systematically with L_IR while the fraction of mergers and interactions increases, as has been observed locally. At comparable luminosities, the fraction of ULIRGs with various morphological classifications is similar at z~2 and z~1. We investigate the position of the ULIRGs, along with 70 LIRGs, on the specific star formation rate versus redshift plane, and find 52 systems to be starbursts (lie more than a factor of 3 above the main sequence relation). The morphologies of starbursts are dominated by interacting and merging systems (50%). If irregular disks are included as potential minor mergers, then we find that up to 73% of starbursts are involved in a merger or interaction at some level. Although the final coalescence of a major merger may not be required for the high luminosities of ULIRGs at z~2 as is the case locally, the large fraction of interactions at all stages and potential minor mergers suggest that the high star formation rates of ULIRGs are still largely externally triggered at z~2.

Models of galaxy evolution assume some connection between the AGN and star formation activity in galaxies. We use the multi-wavelength information of the CDFS to assess this issue. We select the AGNs from the 3Ms XMM-Newton survey and measure the star-formation rates of their hosts using data that probe rest-frame wavelengths longward of 20 um. Star-formation rates are obtained from spectral energy distribution fits, identifying and subtracting an AGN component. We divide the star-formation rates by the stellar masses of the hosts to derive specific star-formation rates (sSFR) and find evidence for a positive correlation between the AGN activity (proxied by the X-ray luminosity) and the sSFR for the most active systems with X-ray luminosities exceeding Lx=10^43 erg/s and redshifts z~1. We do not find evidence for such a correlation for lower luminosity systems or those at lower redshifts. We do not find any correlation between the SFR (or the sSFR) and the X-ray absorption derived from high-quality XMM-Newton spectra either, showing that the absorption is likely to be linked to the nuclear region rather than the host, while the star-formation is not nuclear. Comparing the sSFR of the hosts to the characteristic sSFR of star-forming galaxies at the same redshift we find that the AGNs reside mostly in main-sequence and starburst hosts, reflecting the AGN - sSFR connection. Limiting our analysis to the highest X-ray luminosity AGNs (X-ray QSOs with Lx>10^44 erg/s), we find that the highest-redshift QSOs (with z>2) reside predominantly in starburst hosts, with an average sSFR more than double that of the "main sequence", and we find a few cases of QSOs at z~1.5 with specific star-formation rates compatible with the main-sequence, or even in the "quiescent" region. (abridged)

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

We report the results of the counterpart identification and a detailed analysis of the physical properties of the 48 sources discovered in our deep 1.1mm wavelength imaging survey of the GOODS-South field using the AzTEC instrument on the Atacama Submillimeter Telescope Experiment (ASTE). One or more robust or tentative counterpart candidate is found for 27 and 14 AzTEC sources, respectively, by employing deep radio continuum, Spitzer MIPS & IRAC, and LABOCA 870 micron data. Five of the sources (10%) have two robust counterparts each, supporting the idea that these galaxies are strongly clustered and/or heavily confused. Photometric redshifts and star formation rates (SFRs) are derived by analyzing UV-to-optical and IR-to-radio SEDs. The median redshift of z~2.6 is similar to other earlier estimates, but we show that 80% of the AzTEC-GOODS sources are at z>2, with a significant high redshift tail (20% at z>3.3). Rest-frame UV and optical properties of AzTEC sources are extremely diverse, spanning 10 magnitude in the i- and K-band photometry with median values of i=25.3 and K=22.6 and a broad range of red colour (i-K=0-6). These AzTEC sources are some of the most luminous galaxies in the rest-frame optical bands at z>2, with inferred stellar masses of (1-30) x 10^{10} solar masses and UV-derived star formation rates of SFR(UV) > 10-1000 solar masses per year. The IR-derived SFR, 200-2000 solar masses per year, is independent of redshift or stellar mass. The resulting specific star formation rates, SSFR = 1-100 per Gyr, are 10-100 times higher than similar mass galaxies at z=0, and they extend the previously observed rapid rise in the SSFR with redshift to z=2-5. These galaxies have a SFR high enough to have built up their entire stellar mass within their Hubble time. We find only marginal evidence for an AGN contribution to the near-IR and mid-IR SEDs. (abridged)

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

Using extremely deep PACS 100- and 160um Herschel data from the GOODS-Herschel program, we identify 21 infrared bright galaxies previously missed in the deepest 24um surveys performed by MIPS. These MIPS dropouts are predominantly found in two redshift bins, centred at z ~0.4 and ~1.3. Their S_100/S_24 flux density ratios are similar to those of local LIRGs and ULIRGs, whose silicate absorption features at 18um (at z ~ 0.4) and 9.7um (at z ~ 1.3) are shifted into the 24um MIPS band at these redshifts. The high-z sub-sample consists of 11 infrared luminous sources, accounting for ~2% of the whole GOODS-Herschel sample and putting strong upper limits on the fraction of LIRGs/ULIRGs at 1.0<z<1.7 that are missed by the 24um surveys. We find that a S_100/S_24 > 43 colour cut selects galaxies with a redshift distribution similar to that of the MIPS dropouts and when combined with a second colour cut, S_16/S_8 > 4, isolates sources at 1.0 < z < 1.7. We show that these sources have elevated specific star formation rates (sSFR) compared to main sequence galaxies at these redshifts and are likely to be compact starbursts with moderate/strong 9.7um silicate absorption features in their mid-IR spectra. Herschel data reveal that their infrared luminosities extrapolated from the 24um flux density are underestimated, on average, by a factor of ~3. These silicate break galaxies account for 16% (8%) of the ULIRG (LIRG) population in the GOODS fields, indicating a lower limit in their space density of 2.0 \times 10^(-5) Mpc^(-3). Finally, we provide estimates of the fraction of z < 2 MIPS dropout sources as a function of the 24-, 100-, 160-, 250- and 350um sensitivity limits, and conclude that previous predictions of a population of silicate break galaxies missed by the major 24um extragalactic surveys have been overestimated.

We take advantage of the sensitivity and resolution of Herschel at 100 and 160 micron to directly image the thermal dust emission and investigate the infrared luminosities, L(IR), and dust obscuration of typical star-forming (L*) galaxies at high redshift. Our sample consists of 146 UV-selected galaxies with spectroscopic redshifts 1.5<z<2.6 in the GOODS-North field. Supplemented with deep Very Large Array (VLA) and Spitzer imaging, we construct median stacks at the positions of these galaxies at 24, 100, and 160 micron, and 1.4 GHz. The comparison between these stacked fluxes and a variety of dust templates and calibrations implies that typical star-forming galaxies with UV luminosities L(UV)>1e10 Lsun at z~2 are luminous infrared galaxies (LIRGs) with a median L(IR)=(2.2+/-0.3)e11 Lsun. Typical galaxies at 1.5<z<2.6 have a median dust obscuration L(IR)/L(UV) = 7.1+/-1.1, which corresponds to a dust correction factor, required to recover the bolometric star formation rate (SFR) from the unobscured UV SFR, of 5.2+/-0.6. This result is similar to that inferred from previous investigations of the UV, H-alpha, 24 micron, radio, and X-ray properties of the same galaxies studied here. Stacking in bins of UV slope implies that L* galaxies with redder spectral slopes are also dustier, and that the correlation between UV slope and dustiness is similar to that found for local starburst galaxies. Hence, the rest-frame 30 and 50 micron fluxes validate on average the use of the local UV attenuation curve to recover the dust attenuation of typical star-forming galaxies at high redshift. In the simplest interpretation, the agreement between the local and high redshift UV attenuation curves suggests a similarity in the dust production and stellar and dust geometries of starburst galaxies over the last 10 billion years.

Dust attenuation curves in external galaxies are useful to study their dust properties as well as to interpret their intrinsic spectral energy distributions. In particular the presence or absence of a UV bump at 2175 A remains an open issue which has consequences on the interpretation of broad band colours of distant galaxies. We study the dust attenuation curve in the UV range at z >1. In particular we search for the presence of a UV bump. We use deep photometric data of the CDFS obtained with intermediate and broad band filters by the MUSYC project to sample the UV rest-frame of galaxies with 1<z <2. Herschel/PACS and Spitzer/MIPS data are used to measure the dust emission. 30 galaxies were selected with high S/N in all bands. Their SEDs from the UV to the far-IR are fitted using the CIGALE code and the characteristics of the dust attenuation curve are obtained. The mean dust attenuation curve we derive exhibits a UV bump at 2175A whose amplitude corresponds to 35 % (76%) that of the MW (LMC2 supershell) extinction curve. An analytical expression of the average attenuation curve is given, it is found slightly steeper than the Calzetti et al. one, although at a 1 sigma level. Our galaxy sample is used to study the derivation of the slopes of the UV continuum from broad band colours, including the GALEX FUV-NUV colour. Systematic errors induced by the presence of the bump are quantified. We compare dust attenuation factors measured with CIGALE to the slope of the UV continuum and find that there is a large scatter around the relation valid for local starbursts (0.7 mag). The uncertainties on the determination of the UV slope lead to an extra systematic error of the order of 0.3 to 0.7 mag on dust attenuation when a filter overlaps the UV bump.

We present observations of the molecular gas in the GN20 proto-cluster of galaxies at $z =4.05$ using the Expanded Very Large Array (EVLA). This group of galaxies is the ideal laboratory for studying the formation of massive galaxies via luminous, gas-rich starbursts within 1.6 Gyr of the Big Bang. We detect three galaxies in the proto-cluster in CO 2-1 emission, with gas masses (H$_2$) between $10^{10}$ and $10^{11} \times (\alpha/0.8)$ M$_\odot$. The emission from the brightest source, GN20, is resolved with a size $\sim 2"$, and has a clear north-south velocity gradient, possibly indicating ordered rotation. The gas mass in GN20 is comparable to the stellar mass ($1.3\times 10^{11} \times (\alpha/0.8)$ M$_\odot$ and $2.3\times 10^{11}$ M$_\odot$, respectively), and the sum of gas plus stellar mass is comparable to the dynamical mass of the system ($\sim 3.4\times 10^{11} [sin(i)/sin(45^o)]^{-2}$ M$_\odot$), within a 5kpc radius. There is also evidence for a tidal tail extending another $2"$ north of the galaxy with a narrow velocity dispersion. GN20 may be a massive, gas rich disk that is gravitationally disturbed, but not completely disrupted. There is one Lyman-break galaxy (BD29079) in the GN20 proto-cluster with an optical spectroscopic redshift within our search volume, and we set a 3$\sigma$ limit to the molecular gas mass of this galaxy of $1.1\times 10^{10} \times (\alpha/0.8)$ M$_\odot$.

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24{\mu}m data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24{\mu}m and 15/45 (all MIPS-24{\mu}m detections) are detected at SPIRE-250{\mu}m, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).

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.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

In this contribution, we present some preliminary observational results from the completed ultra-deep survey of 21cm emission from neutral hydrogen at redshifts z=0.164-0.224 with the Westerbork Synthesis Radio Telescope. In two separate fields, a total of 160 individual galaxies has been detected in neutral hydrogen, with HI masses varying from 1.1x10^9 to 4.0x10^10 Msun. The largest galaxies are spatially resolved by the synthesized beam of 23x37 arcsec^2 while the velocity resolution of 19 km/s allowed the HI emission lines to be well resolved. The large scale structure in the surveyed volume is traced well in HI, apart from the highest density regions like the cores of galaxy clusters. All significant HI detections have obvious or plausible optical counterparts which are usually blue late-type galaxies that are UV-bright. One of the observed fields contains a massive Butcher-Oemler cluster but none of the associated blue galaxies has been detected in HI. The data suggest that the lower-luminosity galaxies at z=0.2 are more gas-rich than galaxies of similar luminosities at z=0, pending a careful analysis of the completeness near the detection limit. Optical counterparts of the HI detected galaxies are mostly located in the 'blue cloud' of the galaxy population although several galaxies on the 'red sequence' are also detected in HI. These results hold great promise for future deep 21cm surveys of neutral hydrogen with MeerKAT, APERTIF, ASKAP, and ultimately the Square Kilometre Array.

We present first insights into the far-IR properties for a sample of IRAC and MIPS-24um detected Lyman Break Galaxies (LBGs) at z ~ 3, as derived from observations in the northern field of the Great Observatories Origins Survey (GOODS-N) carried out with the PACS instrument on board the Herschel Space Observatory. Although none of our galaxies are detected by Herschel, we employ a stacking technique to construct, for the first time, the average spectral energy distribution of infrared luminous LBGs from UV to radio wavelengths. We derive a median IR luminosity of L_{IR} = 1.6 x 10^12 Lo, placing the population in the class of ultra luminous infrared galaxies (ULIRGs). Complementing our study with existing multi-wavelength data, we put constraints on the dust temperature of the population and find that for their L_{IR}, MIPS-LBGs are warmer than submm-luminous galaxies while they fall in the locus of the L_{IR}-T_{d} relation of the local ULIRGs. This, along with estimates based on the average SED, explains the marginal detection of LBGs in current sub-mm surveys and suggests that these latter studies introduce a bias towards the detection of colder ULIRGs in the high-z universe, while missing high-z ULIRGS with warmer dust.

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .

We present a search for Herschel-PACS counterparts of dust-obscured, high-redshift objects previously selected at submillimeter and millimeter wavelengths in the Great Observatories Origins Deep Survey North field. We detect 22 of 56 submillimeter galaxies (39%) with a SNR of >=3 at 100 micron down to 3.0 mJy, and/or at 160 micron down to 5.7 mJy. The fraction of SMGs seen at 160 micron is higher than that at 100 micron. About 50% of radio-identified SMGs are associated with PACS sources. We find a trend between the SCUBA/PACS flux ratio and redshift, suggesting that these flux ratios could be used as a coarse redshift indicator. PACS undetected submm/mm selected sources tend to lie at higher redshifts than the PACS detected ones. A total of 12 sources (21% of our SMG sample) remain unidentified and the fact that they are blank fields at Herschel-PACS and VLA 20 cm wavelength may imply higher redshifts for them than for the average SMG population (e.g., z>3-4). The Herschel-PACS imaging of these dust-obscured starbursts at high-redshifts suggests that their far-infrared spectral energy distributions have significantly different shapes than template libraries of local infrared galaxies.

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.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

We present deep {\it Spitzer} mid-infrared spectroscopy, along with 16, 24, 70, and 850\,$\micron$\ photometry, for 22 galaxies located in the Great Observatories Origins Deep Survey-North (GOODS-N) field. The sample spans a redshift range of $0.6\la z \la 2.6$, 24~$\mu$m flux densities between $\sim$0.2$-$1.2 mJy, and consists of submillimeter galaxies (SMGs), X-ray or optically selected active galactic nuclei (AGN), and optically faint ($z_{AB}>25$\,mag) sources. We find that infrared (IR; $8-1000~\micron$) luminosities derived by fitting local spectral energy distributions (SEDs) with 24~$\micron$ photometry alone are well matched to those when additional mid-infrared spectroscopic and longer wavelength photometric data is used for galaxies having $z\la1.4$ and 24~$\micron$-derived IR luminosities typically $\la 3\times 10^{12}~L_{\sun}$. However, for galaxies in the redshift range between $1.4\la z \la 2.6$, typically having 24~$\micron$-derived IR luminosities $\ga 3\times 10^{12}~L_{\sun}$, IR luminosities are overestimated by an average factor of $\sim$5 when SED fitting with 24~$\micron$ photometry alone. This result arises partly due to the fact that high redshift galaxies exhibit aromatic feature equivalent widths that are large compared to local galaxies of similar luminosities. Through a spectral decomposition of mid-infrared spectroscopic data, we are able to isolate the fraction of IR luminosity arising from an AGN as opposed to star formation activity. This fraction is only able to account for $\sim$30\% of the total IR luminosity among the entire sample.

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.

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$.

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<z<1.3 where we estimate to have a redshift completeness higher than 90%, and study the properties and redshift evolution of these sub-populations. We find that, as expected, the relative contributions of AGN and star-forming galaxies to the uJy population depend on the flux density limit of the sample. At all flux levels a significant population of "green-valley" galaxies is observed. While the actual nature of these sources is not definitely understood, the results of this work may suggest that a significant fraction of faint radio sources might be composite (and possibly transition) objects, thus a simple "AGN vs star-forming" classification might not be appropriate to fully understand what faint radio populations really are.

We present the largest spectroscopic follow-up performed in SWIRE ELAIS-N1. We were able to determine redshifts for 289 extragalactic sources. The values of spectroscopic redshifts of the latter have been compared with the estimated values from our photometric redshift code with very good agreement between the two for both galaxies and quasars. Six of the quasars are hyperluminous infrared galaxies all of which are broad line AGN. We have performed emission line diagnostics for 30 sources in order to classify them into star-forming, Seyferts, composite and LINER and compare the results to the predictions from our SED template fitting methods and mid-IR selection methods.

We present a high resolution (down to 0.18"), multi-transition imaging study of the molecular gas in the z = 4.05 submillimeter galaxy GN20. GN20 is one of the most luminous starburst galaxy known at z > 4, and is a member of a rich proto-cluster of galaxies at z = 4.05 in GOODS-North. We have observed the CO 1-0 and 2-1 emission with the VLA, the CO 6-5 emission with the PdBI Interferometer, and the 5-4 emission with CARMA. The H_2 mass derived from the CO 1-0 emission is 1.3 \times 10^{11} (\alpha/0.8) Mo. High resolution imaging of CO 2-1 shows emission distributed over a large area, appearing as partial ring, or disk, of ~ 10kpc diameter. The integrated CO excitation is higher than found in the inner disk of the Milky Way, but lower than that seen in high redshift quasar host galaxies and low redshift starburst nuclei. The VLA CO 2-1 image at 0.2" resolution shows resolved, clumpy structure, with a few brighter clumps with intrinsic sizes ~ 2 kpc. The velocity field determined from the CO 6-5 emission is consistent with a rotating disk with a rotation velocity of ~ 570 km s^{-1} (using an inclination angle of 45^o), from which we derive a dynamical mass of 3 \times 10^{11} \msun within about 4 kpc radius. The star formation distribution, as derived from imaging of the radio synchrotron and dust continuum, is on a similar scale as the molecular gas distribution. The molecular gas and star formation are offset by ~ 1" from the HST I-band emission, implying that the regions of most intense star formation are highly dust-obscured on a scale of ~ 10 kpc. The large spatial extent and ordered rotation of this object suggests that this is not a major merger, but rather a clumpy disk accreting gas rapidly in minor mergers or smoothly from the proto-intracluster medium. ABSTRACT TRUNCATED

We present evidence for very high gas fractions and extended molecular gas reservoirs in normal, near-infrared selected (BzK) galaxies at z~1.5, based on multi-configuration CO[2-1] observations obtained at the IRAM PdBI. Six of the six galaxies observed were securely detected. High resolution observations resolve the CO emission in four of them, implying sizes of order of 6-11 kpc and suggesting the presence of rotation. The UV morphologies are consistent with clumpy, unstable disks, and the UV sizes are consistent with the CO sizes. The star formation efficiencies are homogeneously low and similar to local spirals - the resulting gas depletion times are ~0.5 Gyr, much higher than what is seen in high-z submm galaxies and quasars. The CO luminosities can be predicted to within 0.15 dex from the star formation rates and stellar masses, implying a tight correlation of the gas mass with these quantities. We use dynamical models of clumpy disk galaxies to derive dynamical masses. These models are able to reproduce the peculiar spectral line shapes of the CO emission. After accounting for the stellar and dark matter masses we derive gas masses of 0.4-1.2x10^11 Msun. The conversion factor is very high: alpha_CO=3.6+-0.8, consistent with the Galaxy but four times higher than that of local ultra-luminous IR galaxies. The gas accounts for an impressive 50-65% of the baryons within the galaxies' half light radii. We are witnessing truly gas-dominated galaxies at z~1.5, a finding that explains the high specific SFRs observed for z>1 galaxies. The BzK galaxies can be viewed as scaled-up versions of local disk galaxies, with low efficiency star formation taking place inside extended, low excitation gas disks. They are markedly different than local ULIRGs and high-z submm galaxies, which have more excited and compact gas.

We present results from a multi-wavelength study of 29 sources (false detection probabilities <5%) from a survey of the Great Observatories Origins Deep Survey-North field at 1.1mm using the AzTEC camera. Comparing with existing 850um SCUBA studies in the field, we examine differences in the source populations selected at the two wavelengths. The AzTEC observations uniformly cover the entire survey field to a 1-sigma depth of ~1mJy. Searching deep 1.4GHz VLA, and Spitzer 3--24um catalogues, we identify robust counterparts for 21 1.1mm sources, and tentative associations for the remaining objects. The redshift distribution of AzTEC sources is inferred from available spectroscopic and photometric redshifts. We find a median redshift of z=2.7, somewhat higher than z=2.0 for 850um-selected sources in the same field, and our lowest redshift identification lies at a spectroscopic redshift z=1.1460. We measure the 850um to 1.1mm colour of our sources and do not find evidence for `850um dropouts', which can be explained by the low-SNR of the observations. We also combine these observed colours with spectroscopic redshifts to derive the range of dust temperatures T, and dust emissivity indices $\beta$ for the sample, concluding that existing estimates T~30K and $\beta$~1.75 are consistent with these new data.

We present evidence for Milky-Way-like, low-excitation molecular gas reservoirs in near-IR selected massive galaxies at z~1.5, based on IRAM Plateau de Bure Interferometer CO[3-2] and NRAO Very Large Array CO[1-0] line observations for two galaxies that had been previously detected in CO[2-1] emission. The CO[3-2] flux of BzK-21000 at z=1.522 is comparable within the errors to its CO[2-1] flux, implying that the CO[3-2] transition is significantly sub-thermally excited. The combined CO[1-0] observations of the two sources result in a detection at the 3 sigma level that is consistent with a higher CO[1-0] luminosity than that of CO[2-1]. Contrary to what is observed in submillimeter galaxies and QSOs, in which the CO transitions are thermally excited up to J>=3, these galaxies have low-excitation molecular gas, similar to that in the Milky Way and local spirals. This is the first time that such conditions have been observed at high redshift. A Large Velocity Gradient analysis suggests that molecular clouds with density and kinetic temperature comparable to local spirals can reproduce our observations. The similarity in the CO excitation properties suggests that a high, Milky-Way-like, CO to H_2 conversion factor could be appropriate for these systems. If such low-excitation properties are representative of ordinary galaxies at high redshift, centimeter telescopes such as the Expanded Very Large Array and the longest wavelength Atacama Large Millimeter Array bands will be the best tools for studying the molecular gas content in these systems through the observations of CO emission lines.

We present the results of spectroscopy using HYDRA on the WIYN 3.5m telescope of objects in the deep SWIRE radio field. The goal of the project was to determine spectroscopic redshifts for as many of the brighter objects in the field as possible, especially those detected in the radio and at 24 microns. These redshifts are primarily being used in studies of galaxy evolution and the connection of that evolution to AGN and star-formation. Redshifts measured for 365 individual objects are reported. The redshifts range from 0.03 to 2.5, mostly with z < 0.9. The sources were selected to be within the WIYN HYDRA field of approximately 30' in radius from the center of the SWIRE deep field, 10h46m00s, 59d 01'00" (J2000). Optical sources for spectroscopic observation were selected from a r-band image of the field. A priority list of spectroscopic targets was established in the following order: 20cm detections, 24 micron detections, galaxies with r < 20 and the balance made up of fainter galaxies in the field. We provide a table listing the galaxy positions, measured redshift and error, and note any emission lines that were visible in the spectrum. In practice almost all the galaxies with r < 19 were observed including all of the radio sources and most of the 24 microns sources with r < 20 and a sample of radio sources which had fainter optical counterparts on the r-band image.

(abridged) We present deep Spitzer mid-infrared spectroscopy, along with 16, 24, 70, and 850um photometry, for 22 galaxies located in GOODS-N. The sample spans a redshift range of 0.6 < z < 2.6, 24um flux densities between ~0.2-1.2 mJy, and consists of SMGs, AGN, and optically faint (z_AB > 25) sources. We find that IR luminosities derived by fitting local SEDs with 24um photometry alone are well matched to those when additional mid-infrared spectroscopic and longer wavelength photometric data is used for galaxies having z < 1.4 and 24um-derived IR luminosities typically > 3x10^12 L_sun. However, for galaxies in the redshift range between 1.4 < z < 2.6, typically having 24um-derived IR luminosities > 3x10^12 L_sun, IR luminosities are overestimated by an average factor of ~5 when SED fitting with 24um photometry alone. This result arises partly due to the fact that high redshift galaxies exhibit PAH EQWs that are large compared to local galaxies of similar luminosities. Using improved estimates for the IR luminosities of these sources, we investigate whether their IR emission is found to be in excess relative to that expected based on extinction corrected UV SFRs, possibly suggesting the presence of an obscured AGN. Through a spectral decomposition of IRS spectroscopic data, we are able to isolate the fraction of IR luminosity arising from an AGN as opposed to star formation activity. This fraction is only ~30% of the total IR luminosity among the entire sample, on average. Of the sources identified as having mid-infrared excesses, half are accounted for by using proper bolometric corrections while half show the presence of obscured AGN. We do not find evidence for evolution in the FIR-radio correlation over this redshift range, although the SMGs have IR/radio ratios which are, on average, ~3 times lower than the nominal value.

We present one of the deepest radio continuum surveys to date at a wavelength ~1 meter, at 324.5 MHz. The data reduction and analysis are described and an electronic catalog of the sources detected above 5 sigma is presented. We also discuss the observed angular size distribution for the sample. Using our deeper 20cm survey of the same field, we calculate spectral indices for sources detected in both surveys. The spectral indices for 90cm-selected sources, defined as S ~nu^(-alpha}, shows a peak near 0.7 and only a few sources with very steep spectra. Thus no large population of very steep spectrum microJy sources seems to exist down to the limit of our survey. For 20cm-selected sources, we find similar mean spectral indices for sources with S_20>1 mJy. For weaker sources, below the detection limit for individual sources at 90cm, we use stacking to study the radio spectra. We find that the spectral indices of small (<3") 20cm-selected sources with S_20< 10 mJy have mean and median alpha(90,20)~0.3-0.5. This is flatter than the spectral indices of the stronger source population. We report log N-log S counts at 90cm which show a flattening below 5 mJy. Given the median redshift of the population, z~1, the spectral flattening and the flattening of the log N-log S counts occurs at radio luminosities normally associated with AGN rather than with galaxies dominated by star-formation.

We report the detection of a CO emission line from the submillimiter galaxy (SMG) GN10 in the GOODS-N field. GN10 lacks any counterpart in extremely deep optical and near-IR imaging obtained with the Hubble Space Telescope and ground-based facilities. This is a prototypical case of a source that is extremely obscured by dust, for which it is practically impossible to derive a spectroscopic redshift in the optical/near-IR. Under the hypothesis that GN10 is part of a proto-cluster structure previously identified at z~4.05 in the same field, we searched for CO[4-3] at 91.4 GHz with the IRAM Plateau de Bure Interferometer, and successfully detected a line. We find that the most likely redshift identification is z=4.0424+-0.0013, based on: 1) the very low chance that the CO line is actually serendipitous from a different redshift; 2) a radio-IR photometric redshift analysis; 3) the identical radio-IR SED, within a scaling factor, of two other SMGs at the same redshift. The faintness at optical/near-IR wavelengths requires an attenuation of A_V~5-7.5 mag. This result supports the case that a substantial population of very high-z SMGs exists that had been missed by previous spectroscopic surveys. This is the first time that a CO emission line has been detected for a galaxy that is invisible in the optical and near-IR. Our work demonstrates the power of existing and planned facilities for completing the census of star formation and stellar mass in the distant Universe by measuring redshifts of the most obscured galaxies through millimeter spectroscopy.

We use extensive new observations of the very rich z ~ 0.4 cluster of galaxies A851 to examine the nature and origin of starburst galaxies in intermediate-redshift clusters. New HST observations, Spitzer photometry and ground-based spectroscopy cover most of a region of the cluster about 10 arcmin across, corresponding to a clustercentric radial distance of about 1.6 Mpc. This spatial coverage allows us to confirm the existence of a morphology-density relation within this cluster, and to identify several large, presumably infalling, subsystems. We confirm our previous conclusion that a very large fraction of the starforming galaxies in A851 have recently undergone starbursts. We argue that starbursts are mostly confined to two kinds of sites: infalling groups and the cluster center. At the cluster center it appears that infalling galaxies are undergoing major mergers, resulting in starbursts whose optical emission lines are completely buried beneath dust. The aftermath of this process appears to be proto-S0 galaxies devoid of star formation. In contrast, major mergers do not appear to be the cause of most of the starbursts in infalling groups, and fewer of these events result in the transformation of the galaxy into an S0. Some recent theoretical work provides possible explanations for these two distinct processes, but it is not clear whether they can operate with the very high efficiency needed to account for the very large starburst rate observed.

We present the serendipitous discovery of z=4.05 molecular gas CO emission lines with the IRAM Plateau de Bure interferometer coincident with GN20 and GN20.2, two luminous submillimeter galaxies (SMGs) in the Great Observatories Origins Deep Survey North field (GOODS-N). These are among the most distant submillimeter-selected galaxies reliably identified through CO emission and also some of the most luminous known. In terms of CO to bolometric luminosities, stellar mass and star formation rates (SFRs), these newly discovered z>4 SMGs are similar to z~1.5-3 SMGs studied to date. These z~4 SMGs have much higher specific SFRs than typical B-band dropout Lyman break galaxies at the same redshift. The stellar mass-SFR correlation for normal galaxies does not seem to evolve much further, between z~2 and z~4. A significant z=4.05 spectroscopic redshift spike is observed in GOODS-N, and a strong spatial overdensity of B-band dropouts and IRAC selected z>3.5 galaxies appears to be centered on the GN20 and GN20.2 galaxies. This suggests a proto-cluster structure with total mass ~10^14 Msun. Using photometry at mid-IR, submm and radio wavelengths, we show that reliable photometric redshifts (Dz/(1+z)~0.1) can be derived for SMGs over 1<z<4. This new photometric redshift technique has been used to provide a first estimate of the space density of 3.5<z<6 hyper-luminous starburst galaxies, and to show that they contribute substantially to the SFR density at early epochs. Many of these high-redshift starbursts will be within reach of Herschel. We find that the radio to mid-IR flux density ratio can be used to select z>3.5 starbursts, regardless of their submm/mm emission [abridged].

We present detailed observations of a z~1.99 cluster of submillimeter galaxies (SMGs), discovered as the strongest redshift spike in our entire survey of ~100 SMGs across 800 square arcmin. It is the largest blank-field SMG concentration currently known and has <0.01% chance of being drawn from the underlying selection function for SMGs. We have compared UV observations of galaxies at this redshift, where we find a much less dramatic overdensity, having an 11% chance of being drawn from its selection function. We use this z~1.99 overdensity to compare the biasing of UV- and submm-selected galaxies, and test whether SMGs could reside in less overdense environments, with their apparent clustering signal being dominated by highly active merger periods in modest mass structures. This impressively active formation phase in a low mass cluster is not something seen in simulations, although we propose a toy model using merger bias which could account for the bias seen in the SMGs. While enhanced buildup of stellar mass appears characteristic of other high-z galaxy clusters, neither the UV- nor submm-galaxies in this structure exhibit larger stellar masses than their field galaxy counterparts (although the excess of SMGs in the structure represents a larger volume-averaged stellar mass than the field). Our findings have strong implications for future surveys for high-z galaxies at long wavelengths such as SCUBA2 and Herschel. We suggest that since these surveys will select galaxies during their episodes of peak starbursts, they could probe a much wider range of environments than just the progenitors of rich clusters, revealing more completely the key events and stages in galaxy formation and assembly.

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.

We use deep far-IR, submm, radio and X-ray imaging and mid-IR spectroscopy to explore the nature of a sample of Spitzer-selected dust-obscured galaxies (DOGs) in GOODS-N. A sample of 79 galaxies satisfy the criteria R-[24]>14 (Vega) down to S24>100 microJy. Twelve of these galaxies have IRS spectra available which we use to measure redshifts and classify these objects as being dominated by star formation or active galactic nuclei (AGN) activity in the mid-IR. The IRS spectra and Spitzer photometric redshifts confirm that the DOGs lie in a tight redshift distribution around z~2. Based on mid-IR colors, 80% of DOGs are likely dominated by star formation; the stacked X-ray emission from this sub-sample of DOGs is also consistent with star formation. Since only a small number of DOGs are individually detected at far-IR and submm wavelengths, we use a stacking analysis to determine the average flux from these objects and plot a composite IR (8-1000 microns) spectral energy distribution (SED). The average luminosity of these star forming DOGs is LIR~1e12 Lsun. We compare the average star forming DOG to the average bright (S850>5 mJy) submillimeter galaxy (SMG); the S24>100 microJy DOGs are 3 times more numerous but 8 times less luminous in the IR. The far-IR SED shape of DOGs is similar to that of SMGs (average dust temperature of around 30 K) but DOGs have a higher mid-IR to far-IR flux ratio. The average star formation-dominated DOG has a star formation rate of 200 Msun/yr which, given their space density, amounts to a contribution of 0.01 Msun/yr/Mpc3 (or 5-10%) to the star formation rate density at z~2. We use the composite SED to predict the average flux of DOGs in future Herschel/PACS 100 micron and SCUBA-2 450 micron surveys and show that the majority of them will be detected.

High-resolution submm imaging of the HzRG, 4C60.07, at z=3.8, has revealed two dusty components. Spitzer imaging shows that one of these components (B) is coincident with an extremely red AGN, offset by ~4" (~30 kpc) from the HzRG core. The other submm component (A) - resolved by our beam and devoid of emission at 3.6-8.0um - lies between B and the HzRG core. Since the HzRG was discovered via its young, steep-spectrum lobes and their creation was likely triggered by the interaction, we argue that we are witnessing an early-stage merger, prior to its eventual equilibrium state. The interaction is between the host galaxy of an actively-fueled BH, and a gas-rich starburst/AGN (B) marked by the compact submm component and coincident with broad CO emission. `A' is a plume of cold, dusty gas, associated with a narrow (~150 km/s) CO feature, and may represent a short-lived tidal structure. It has been claimed that HzRGs and SMGs differ only in the activity of their AGNs, but such complex submm morphologies are seen only rarely amongst SMGs. Our study has important implications: where a galaxy's gas is not aligned with its central BH, CO may be an unreliable probe of dynamical mass, affecting work on the co-assembly of BHs and spheroids. Our data support the picture wherein close binary AGN are induced by mergers. They also raise the possibility that some supposedly jet-induced starbursts may have formed co-evally with (yet independently of) the radio jets, both triggered by the same interaction. We note that the HzRG host would have gone unnoticed without its jets/companion, so there may be many other unseen BHs at high redshift, lost in the sea of ~5 x 10^8 similarly bright IRAC sources - sufficiently massive to drive a >10^27-W radio source, yet practically invisible unless actively fueled (abridged).

We have conducted a deep and uniform 1.1 mm survey of the GOODS-N field with AzTEC on the James Clerk Maxwell Telescope (JCMT). Here we present the first results from this survey including maps, the source catalogue, and 1.1 mm number-counts. The results presented here were obtained from a 245 sq-arcmin region with near uniform coverage to a depth of 0.96-1.16 mJy/beam. Our robust catalogue contains 28 source candidates detected with S/N >= 3.75, only 1-2 of which are expected to be spurious detections. Of these source candidates, 8 are also detected by SCUBA at 850 um in regions where there is good overlap between the two surveys. The major advantage of our survey over that with SCUBA is the uniformity of coverage. We calculate number counts using two different techniques: the first using a frequentist parameter estimation, and the second using a Bayesian method. The two sets of results are in good agreement. We find that the 1.1 mm differential number counts are well described in the 2-6 mJy range by the functional form dN/dS = N' (S'/S) exp(-S/S') with fitted parameters S' = 1.25 +/-0.38 mJy and dN/dS = 300 +/- 90 per mJy per sq-deg at 3 mJy.

We present evidence for the existence of an IRAC excess in the spectral energy distribution (SED) of 5 galaxies at 0.6<z<0.9 and 1 galaxy at z=1.7. These 6 galaxies, located in the Great Observatories Origins Deep Survey field (GOODS-N), are star forming since they present strong 6.2, 7.7, and 11.3 um polycyclic aromatic hydrocarbon (PAH) lines in their Spitzer IRS mid-infrared spectra. We use a library of templates computed with PEGASE.2 to fit their multiwavelength photometry and derive their stellar continuum. Subtraction of the stellar continuum enables us to detect in 5 galaxies a significant excess in the IRAC band pass where the 3.3 um PAH is expected. We then assess if the physical origin of the IRAC excess is due to an obscured active galactic nucleus (AGN) or warm dust emission. For one galaxy evidence of an obscured AGN is found, while the remaining four do not exhibit any significant AGN activity. Possible contamination by warm dust continuum of unknown origin as found in the Galactic diffuse emission is discussed. The properties of such a continuum would have to be different from the local Universe to explain the measured IRAC excess, but we cannot definitively rule out this possibility until its origin is understood. Assuming that the IRAC excess is dominated by the 3.3 um PAH feature, we find good agreement with the observed 11.3 um PAH line flux arising from the same C-H bending and stretching modes, consistent with model expectations. Finally, the IRAC excess appears to be correlated with the star-formation rate in the galaxies. Hence it could provide a powerful diagnostic for measuring dusty star formation in z>3 galaxies once the mid-infrared spectroscopic capabilities of the James Webb Space Telescope become available.

We have obtained a position (at sub-arcsecond accuracy) of the submillimeter bright source GOODS 850-5 (also known as GN10) in the GOODS North field using the IRAM Plateau de Bure interferometer at 1.25 mm wavelengths (MM J123633+6214.1, flux density: S(1.25 mm)=5.0+-1.0 mJy). This source has no optical counterpart in deep ACS imaging down to a limiting magnitude of i(775)=28.4 mag and its position is coincident with the position found in recent sub-millimeter mapping obtained at the SMA (Wang et al. 2007). Using deep VLA imaging at 20 cm, we find a radio source (S(20 cm)=32.7+-4.3 microJy) at the same position that is significantly brighter than reported in Wang et al. The source is detected by Spitzer in IRAC as well as at 24 microns. We apply different photometric redshift estimators using measurements of the dusty, mid/far-infrared part of the SED and derive a redshift z~4. Given our detection in the millimeter and radio we consider a significantly higher redshift (e.g., z~6 Wang et al. 2007) unlikely. MM J123633+6214.1 alias GOODS 850-5 nevertheless constitutes a bright representative of the high-redshift tail of the submillimeter galaxy population that may contribute a significant fraction to the (sub)millimeter background.

We present the first detection of molecular gas cooling CO emission lines from ordinary massive galaxies at z=1.5. Two sources were observed with the IRAM Plateau de Bure Interferometer, selected to lie in the mass-star formation rate correlation at their redshift, thus being representative of massive high-z galaxies. Both sources were detected with high confidence, yielding L'_CO~10^{10}K km/s pc^2. For one of the sources we find evidence for velocity shear, implying CO sizes of ~10 kpc. With an infrared luminosity of L_FIR~10^{12}L_sun, these disk-like galaxies are borderline ULIRGs but with star formation efficiency similar to that of local spirals, and an order of magnitude lower than that in submm galaxies. This suggests a CO to total gas conversion factor similar to local spirals, gas consumption timescales approaching 1 Gyr or longer and molecular gas masses reaching ~10^11 M_sun, comparable to or larger than the estimated stellar masses. These results support a major role of 'in situ' gas consumption over cosmological timescales and with relatively low star formation efficiency, analogous to that of local spiral disks, for the formation of today's most massive galaxies and their central black holes. Given the high space density of similar galaxies, ~10^{-4}/Mpc^3, this implies a widespread presence of gas rich galaxies in the early Universe, many of which might be within reach of detailed investigations of current and planned facilities.

We present the IR luminosity function derived from ultra-deep 70 micron imaging of the GOODS-North field. The 70 micron observations are longward of the PAH and silicate features which complicate work in the MIR. We derive far-infrared luminosities for the 143 sources with S_{70} > 2 mJy (S/N > 3 \sigma). The majority (81%) of the sources have spectroscopic redshifts, and photometric redshifts are calculated for the remainder. The IR luminosity function at four redshifts (z ~ 0.28, 0.48, 0.78, and 0.97) is derived and compared to the local one. There is considerable degeneracy between luminosity and density evolution. If the evolving luminosity function is described as \rho(L, z) = (1 + z)^q \rho(L/(1 + z)^p, 0), we find q = -2.19p + 6.09. In the case of pure luminosity evolution, we find a best fit of p = 2.78^{+0.34}_{-0.32}. This is consistent with the results from 24 micron and 1.4 GHz studies. Our results confirm the emerging picture of strong evolution in LIRGs and ULIRGs at 0.4 < z < 1.1, but we find no evidence of significant evolution in the sub-LIRG (L < 10^{11} L_{\odot}) population for z < 0.4.

Approximately 20-30% of 1.4<z<2.5 galaxies with K<22 (Vega) detected with Spitzer MIPS at 24um show excess mid-IR emission relative to that expected based on the rates of star formation measured from other multiwavelength data.These galaxies also display some near-IR excess in Spitzer IRAC data, with a spectral energy distribution peaking longward of 1.6um in the rest frame, indicating the presence of warm-dust emission usually absent in star forming galaxies. Stacking Chandra data for the mid-IR excess galaxies yields a significant hard X-ray detection at rest-frame energies >6.2 keV. The stacked X-ray spectrum rises steeply at >10 keV, suggesting that these sources host Compton-thick Active Galactic Nuclei (AGNs) with column densities N_H~10^{24} cm^-2 and an average, unobscured X-ray luminosity L_{2-8 keV}~(1-4)x10^43 erg/s. Their sky density(~3200 deg^-2) and space density (~2.6x10^-4 Mpc^-3) are twice those of X-ray detected AGNs at z~2, and much larger than those of previously-known Compton thick sources at similar redshifts. The mid-IR excess galaxies are part of the long sought-after population of distant heavily obscured AGNs predicted by synthesis models of the X-ray background. The fraction of mid-IR excess objects increases with galaxy mass, reaching ~50-60% for M~10^11 M_sun, an effect likely connected with downsizing in galaxy formation. The ratio of theinferred black hole growth rate from these Compton-thick sources to the global star formation rate at z=2 is similar to the mass ratio of black holes to stars inlocal spheroids, implying concurrent growth of both within the precursors oftoday's massive galaxies.

Examining a sample of massive galaxies at 1.4<z<2.5 with K_{Vega}<22 from the Great Observatories Origins Deep Survey, we compare photometry from Spitzer at mid- and far-IR, to submillimeter, radio and rest-frame ultraviolet wavelengths, to test the agreement between different tracers of star formation rates (SFRs) and to explore the implications for galaxy assembly. For z~2 galaxies with moderate luminosities(L_{8um}<10^{11}L_sun), we find that the SFR can be estimated consistently from the multiwavelength data based on local luminosity correlations. However,20--30% of massive galaxies, and nearly all those with L_{8um}>10^{11}L_sun, show a mid-IR excess which is likely due to the presence of obscured active nuclei, as shown in a companion paper. There is a tight and roughly linear correlation between stellar mass and SFR for 24um-detected galaxies. For a given mass, the SFR at z=2 was larger by a factor of ~4 and ~30 relative to that in star forming galaxies at z=1 and z=0, respectively. Typical ultraluminous infrared galaxies (ULIRGs) at z=2 are relatively 'transparent' to ultraviolet light, and their activity is long lived (~400 Myr), unlike that in local ULIRGs and high redshift submillimeter-selected galaxies. ULIRGs are the common mode of star formation in massive galaxies at z=2, and the high duty cycle suggests that major mergers are not the dominant trigger for this activity.Current galaxy formation models underpredict the normalization of the mass-SFR correlation by about a factor of 4, and the space density of ULIRGs by an orderof magnitude, but give better agreement for z>1.4 quiescent galaxies.

We present SEDs, Spitzer colours, and IR luminosities for 850 micron selected galaxies in the GOODS-N field. Using the deep Spitzer Legacy images and new data and reductions of the VLA-HDF radio data, we find statistically secure counterparts for 60 per cent (21/35) of our submm sample, and identify tentative counterparts for another 12 objects. This is the largest sample of submm galaxies with statistically secure counterparts detected in the radio and with Spitzer. We find that in most cases the 850 micron emission is dominated by a single 24 micron source. A composite rest-frame SED shows that the submm sources peak at longer wavelengths than those of local ULIRGs of the same luminosity and therefore appear to be cooler. The SEDs of submm galaxies are also different from those of their high redshift neighbours, the near-IR selected BzK galaxies, whose mid-IR to radio SEDs are more like those of local ULIRGs. Using 24 micron. 850 micron and 1.4 GHz observations, we fit templates that span the mid-IR through radio to derived the integrated IR luminosity of the submm galaxies. By themselves, 24 micron and radio fluxes are able to predict LIR reasonably well because they are relatively insensitive to temperature. However, the submm flux by itself consistently overpredicts LIR when using spectral templates which obey the local ULIRG temperature-luminosity relation. The shorter Spitzer wavelengths sample the stellar bump at the redshifts of the submm sources, and we find that the Spitzer photometry alone provides a model independent estimate of the redshift. Using X-ray and mid-IR data, only 5 per cent of our secure counterparts show strong evidence for an active galactic nucleus dominating the IR luminosity.

We report on a preliminary study concerning the origin of radio emission within radio galaxies at L(1.4GHz)>1E24 W/Hz in the GOODS-N field. In the local universe, Condon et al. (2002) and Yun et al. (2001) have shown that in galaxies with radio luminosities greater than 1E23 W/Hz the majority of the radio emission originates from a `monster' i.e., an AGN. Using the Chandra 2Msec X-ray image centered on the GOODS-N field and a reprocessed VLA HDF A-array data plus newly acquired VLA B-array data (rms=5.3microJy), we find that radio galaxies (with spectroscopic redshifts; all have z>1) with L(1.4GHz)>1E24 W/Hz typically have an X-ray detection rate of 72% (60% emit hard X-rays suggesting an AGN origin for the radio emission) in contrast to 25% for radio galaxies with L < 1E23 W/Hz. The ACS images of these L(1.4 GHz) > 1E24 W/Hz galaxies typically show compact rather than extended galaxy morphology which is generally found for the less luminous radio emitting galaxies but a few appear to be ongoing galaxy mergers. We also present SED fitting for these luminous radio galaxies including Spitzer IRAC & MIPS 24um photometry and 60% show distinct power-law SED indicative of an AGN. Initial results tell us that the X-ray emitting radio galaxy population are generally not submm sources but the few (~10%) that are SCUBA sources appear to be the small AGN population found by Pope et al. and others.

How does the submm galaxy population detected by SCUBA fit into galaxy evolution? How do these rare starbursting systems, which contribute significantly to high redshift star-formation, relate to other galaxy populations? Deep radio observations have been most useful for studying these systems, but still leave a significant fraction of the population unidentified. Now with the deep IRAC and MIPS images from the GOODS Spitzer Legacy program and a re-analysis of the deep radio data, we are able to identify counterparts for a large fraction of SCUBA galaxies in GOODS-N. All of these counterparts are detected by Spitzer. Given the vast multi-wavelength data in this field, we can study the spectral energy distributions (SEDs) of these systems and determine what is fueling their intense infrared luminosities. A rest-frame composite optical-to-radio SED for all spectroscopically identified submm sources shows that the average SCUBA galaxy is consistent with models of ultraluminous starburst galaxies, although cooler than those observed locally. Because of this, the submm flux alone consistently overestimates IR when using spectral templates which obey the local ULIRG temperature-luminosity relation. The wide range of 24/850 micron flux ratios as a function of redshift indicates the presence of strong mid-IR features, to be confirmed with deep IRS spectroscopy. The IRAC colours of the submm systems provide useful redshift constraints, since, at these redshifts, IRAC samples the stellar bump. The Spitzer photometry of this large sample of submm galaxies has allowed us to put constraints on many of the outstanding issues in submm astronomy.

Using the large multi-wavelength data set in the chandra/SWIRE Survey (0.6 square degrees in the Lockman Hole), we show evidence for the existence of highly obscured (Compton-thick) AGN, estimate a lower limit to their surface density and characterize their multi-wavelength properties. Two independent selection methods based on the X-ray and infrared spectral properties are presented. The two selected samples contain 1) 5 X-ray sources with hard X-ray spectra and column densities > 10^24 cm-2, and 2) 120 infrared sources with red and AGN-dominated infrared spectral energy distributions (SEDs). We estimate a surface density of at least 25 Compton-thick AGN per square degree detected in the infrared in the chandra/SWIRE field of which ~40% show distinct AGN signatures in their optical/near-infrared SEDs, the remainings being dominated by the host-galaxy emission. Only ~33% of all Compton-thick AGN are detected in the X-rays at our depth (F(0.3-8 keV)>10^-15 erg/cm2/s. We report the discovery of two sources in our sample of Compton-thick AGN, SWIRE_J104409.95+585224.8 (z=2.54) and SWIRE_J104406.30+583954.1 (z=2.43), which are the most luminous Compton-thick AGN at high-z currently known. The properties of these two sources are discussed in detail with an analysis of their spectra, SEDs, luminosities and black-hole masses.

We present evidence that the mid infrared (MIR) is a good tracer of the total infrared luminosity, L(IR), and star formation rate (SFR), of galaxies up to z 1.3. We use deep MIR images from the Infrared Space Observatory (ISO) and the Spitzer Space Telescope in the Northern field of the Great Observatories Origins Deep Survey (GOODS-N) together with VLA radio data to compute three independant estimates of L(IR). The L(IR,MIR) derived from the observed 15 and/or 24 um flux densities using a library of template SEDs, and L(IR,radio), derived from the radio (1.4 and/or 8.5 GHz) using the radio-far infrared correlation, agree with a 1-sigma dispersion of 40 %. We use the k-correction as a tool to probe different parts of the MIR spectral energy distribution (SED) of galaxies as a function of their redshift and find that on average distant galaxies present MIR SEDs very similar to local ones. However, in the redshift range z= 0.4-1.2, L(IR,24um) is in better agreement with L(IR,radio) than L(IR,15 um) by 20 %, suggesting that the warm dust continuum is a better tracer of the SFR than the broad emission features due to polycyclic aromatic hydrocarbons (PAHs). We find marginal evidence for an evolution with redshift of the MIR SEDs: two thirds of the distant galaxies exhibit rest-frame MIR colors (L(12 um)/L(7 um) and L(10 um)/L(15 um) luminosity ratios) below the median value measured for local galaxies. Possible explanations are examined but these results are not sufficient to constrain the physics of the emitting regions. We compare three commonly used SED libraries which reproduce the color-luminosity correlations of local galaxies with our data and discuss possible refinements to the relative intensities of PAHs, warm dust continuum and silicate absorption.

We use radio, near-IR, optical, and X-ray observations to examine dynamic processes in the central region of Abell 2125. In addition to the central triple, including members of both major dynamical subsystems identified from a redshift survey, this region features a galaxy showing strong evidence for ongoing gas stripping during a high-velocity passage through the gas in the cluster core. The disk galaxy C153 exhibits a plume stretching toward the cluster center seen in soft X-rays by Chandra, parts of which are also seen in [O II] emission and near-UV continuum light. HST imaging shows a distorted disk, with star-forming knots asymmetrically distributed and remnant spiral structure possibly defined by dust lanes. The stars and ionized gas in its disk are kinematically decoupled, demonstrating that pressure stripping must be important, and that tidal disruption is not the only mechanism at work. Comparison of the gas properties seen in the X-ray and optical data on the plume highlight significant features of the history of stripped gas in the intracluster medium. The nucleus of C153 also hosts an AGN, shown by the weak and distorted extended radio emission and a radio compact core. The unusual strength of the stripping signatures in this instance is likely related to the high relative velocity of the galaxy with respect to the intracluster medium, during a cluster/cluster merger, and its passage very near the core of the cluster. Another sign of recent dynamical events is diffuse starlight asymmetrically placed about the central triple in a cD envelope. Transient and extreme dynamical events as seen in Abell 2125 may be important drivers of galaxy evolution in the cores of rich clusters.

We investigate the multi-wavelength emission of BzK selected star forming galaxies at z~2 in the Great Observatories Origins Deep Survey (GOODS) North region. Most (82%) of the sources are individually detected at 24um in the Spitzer MIPS imaging, and one fourth (26%) in the VLA radio data. Significant detections of the individually undetected objects are obtained through stacking in the radio, submm and X-ray domains. The typical star forming galaxy with stellar mass ~10^{11}Mo at z=2 is an Ultra-luminous Infrared Galaxy (ULIRG), with L_IR ~ 1-2x10^{12}Lo and star formation rate SFR 200-300Mo/yr, implying a comoving density of ULIRGs at z=2 at least 3 orders of magnitude above the local one. SFRs derived from the reddening corrected UV luminosities agree well, on average, with the longer wavelength estimates. The high 24um detection rate suggests a relatively large duty cycle for the BzK star forming phase, consistently with the available independent measurements of the space density of passively evolving galaxies at z>1.4. If the IMF at z=2 is similar to the local one, and in particular is not a top-heavy IMF, this suggests that a substantial fraction of the high mass tail (>10^{11}Mo) of the galaxy stellar mass function was completed by z~1.4.

A spectroscopic follow-up of SWIRE sources from the Lockman Validation Field has allowed the determination of the SWIRE 24 micron Luminosity Function (LF). The spectroscopic sample was chosen above a 24 micron flux limit at 260 microJy and an r-band optical limit of r<21. A spectroscopic completeness of 82.5% was achieved. We found the median redshift for the sample to be z_med=0.29. Markov-Chain Monte-Carlo (MCMC) and Maximum Likelihood Estimator (MLE) techniques were employed to fit a parametric LF. Our result of the local LF (LLF) is consistent with the local 25 micron determination from Shupe et al. 1998. We split the sample at a redshift of z_split=0.36 and find strong evidence for galaxy evolution.

We discuss optical associations, spectral energy distributions and photometric redshifts for SWIRE sources in the ELAIS-N1 area and the Lockman Validation Field. The band-merged IRAC (3.6, 4.5, 5.8 and 8.0 mu) and MIPS (24, 70, 160 mu) data have been associated with optical UgriZ data from the INT Wide Field Survey in ELAIS-N1, and with our own optical Ugri data in Lockman-VF. The spectral energy distributions of selected ELAIS sources in N1 detected by SWIRE, most with spectroscopic redshifts, are modelled in terms of a simple set of galaxy and quasar templates in the optical and near infrared, and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst, and AGN dust torus) in the mid infrared. The optical data, together with the IRAC 3.6 and 4.5 mu data, have been used to determine photometric redshifts. For galaxies with known spectroscopic redshifts there is a notable improvement in the photometric redshift when the IRAC data are used, with a reduction in the rms scatter from 10% in (1+z) to 7%. The photometric redshifts are used to derive the 3.6 and 24 mu redshift distribution and to compare this with the predictions of models. For those sources with a clear mid infrared excess, relative to the galaxy starlight model used for the optical and near infrared, the mid and far infrared data are modelled in terms of the same dust emission templates. The proportions found of each template type are: cirrus 31%, M82 29%, Arp 220 10%, AGN dust tori 29%. The distribution of the different infrared sed types in the L_{ir}/L_{opt} versus L_{ir} plane, where L_{ir} and L_{opt} are the infrared and optical bolometric luminosities, is discussed.

We present a description of deep radio, optical, and near IR observations taken with the VLA, the KPNO 2m and the KPNO 4m of the region containing the rich cluster of galaxies Abell 2125. The reduction of each dataset is described. A catalog of radio sources apparently not associated with members of Abell 2125 and the associated R-band magnitudes is presented.

Using our extensive radio, optical, near-IR and X-ray imaging and spectroscopy, we consider the reason for the unusually large number of radio detected galaxies, mostly found outside the cluster core, in Abell 2125 (z=0.2465, richness class 4). With 20-cm VLA data, we detect continuum emission from 90 cluster members. The multiwavelength properties of these galaxies suggest that most of the radio emission is due to an enhanced star-formation rate. The dynamical study of Miller et al (2004) suggests that Abell 2125 is undergoing a major cluster-cluster merger, with our view within 30 degrees of the merger axis and within 0.2 Gyr of core passage. The combination of projection effects and the physical processes at work during this special time in the cluster's evolution seem likely to be responsible for the unusual level of activity we see in the cluster. We argue that tidal effects on individual cluster members, often far from the cluster core, are responsible for the increased star formation. Our results are consistent with the idea that disk galaxies during this phase of a cluster's evolution undergo rapid evolution, through a burst of star formation, on their way to becoming S0's.

We characterize the SWIRE galaxy populations in the SWIRE validation field within the Lockman Hole, based on the 3.6-24$\mu$ Spitzer data and deep U,g',r',r' optical imaging within an area ~1/3 sq. deg for ~16,000 Spitzer-SWIRE sources. The entire SWIRE survey will discover over 2.3 million galaxies at 3.6$\mu$m and almost 350,000 at 24$\mu$m; ~70,000 of these will be 5-band 3.6-24$\mu$ detections. The colors cover a broad range, generally well represented by redshifted spectral energy distributions of known galaxy populations, however significant samples of unusually blue objects in the [3.6-4.5]$\mu$m color are found, as well as many objects very red in the 3.6-24$\mu$m mid-IR. Nine of these are investigated and are interpreted as star-forming systems, starbursts and AGN from z=0.37 to 2.8, with luminosities from L$_{IR}$=10$^{10.3}$ to 10$^{13.7}$ L$_{\odot}$

The SIRTF Wide-area InfraRed Extragalactic survey (SWIRE), the largest SIRTF Legacy program, is a wide-area, imaging survey to trace the evolution of dusty, star-forming galaxies, evolved stellar populations, and AGN as a function of environment, from redshifts z~3 to the current epoch. SWIRE will survey 7 high-latitude fields, totaling 60 - 65 sq. deg. in all 7 SIRTF bands: IRAC 3.6, 4.5, 5.6, 8 microns and MIPS 24, 70, 160 microns. The Legacy Extragalactic Catalog may contain in excess of 2 million IR-selected galaxies, dominated by (1) ~150,000 luminous infrared galaxies (LIRGs: L{FIR}>10^11 L_sun), ~7000 of these with z>2; (2) 1 million early-type galaxies, ~10,000 with z>2; and (3) \~20,000 classical AGN, plus significantly more dust-obscured QSO/AGN among the LIRGs. SWIRE will provide an unprecedented view of the evolution of galaxies, structure, and AGN. The key scientific goals of SWIRE are: (1) to determine the evolution of actively star-forming and passively evolving galaxies in order to understand the history of galaxy formation in the context of cosmic structure formation; (2) to determine the evolution of the spatial distribution and clustering of evolved galaxies, starbursts and AGN in the key redshift range, 0.5<z<3, over which much of cosmic evolution has occurred; (3) to determine the evolutionary relationship between ``normal galaxies'' and AGN, and the contribution of AGN accretion energy vs stellar nucleosynthesis to the cosmic backgrounds. SWIRE's large area is important to create statistically significant population samples over enough volume cells that we can resolve the star formation history as a function of epoch and environment The large volume is also optimised for finding rare objects.

(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 <f(SBRG)> = 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.

Radio observations were used to detect the `active' galaxy population within rich clusters of galaxies in a non-biased manner that is not plagued by dust extinction or the K-correction. We present wide-field radio, optical (imaging and spectroscopy), and ROSAT All-Sky Survey (RASS) X-ray data for a sample of 30 very rich Abell (R > 2) cluster with z < 0.25. The VLA radio data samples the ultra-faint radio (L(1.4 GHz) > 2E22 W/Hz) galaxy population within these extremely rich clusters for galaxies with M_R < -21. This is the largest sample of low luminosity 20 cm radio galaxies within rich Abell clusters collected to date. The radio-selected galaxy sample represents the starburst (Star formation rate > 5 M_sun/yr) and active galactic nuclei (AGN) populations contained within each cluster. Archival and newly acquired redshifts were used to verify cluster membership for most (~95%) of the optical identifications. Thus we can identify all the starbursting galaxies within these clusters, regardless of the level of dust obscuration that would affect these galaxies being identified from their optical signature. Cluster sample selection, observations, and data reduction techniques for all wavelengths are discussed.

We present a multi-wavelength study of Extremely Red Objects (EROs) employing deep RIzJHK photometry of a 8.5'x8.5' region to identify 68 EROs with R-K>5.3 and K<20.5 (5-sigma). This is combined with an extremely deep 1.4-GHz radio map (sigma=3.5uJy), sensitive enough to detect an active galaxy with L_1.4>10^23 W/Hz at z>1 or a SFR of >25Mo/yr. We detect radio emission from 21 EROs at >12.6uJy and resolve a third of these with our 1.6'' FWHM beam. The SEDs of most of these radio EROs are consistent with dust-reddened starbursts at z~1. At z~1 the radio luminosities of these EROs indicate far-infrared luminosities of L_FIR>10^12 Lo, meaning half are ultraluminous infrared galaxies. We conclude that >16+/-5% of EROs with K<20.5 are luminous infrared galaxies at z>1. We also photometrically classify the EROs to investigate the mix of dusty/active and evolved/passive systems in the radio-undetected EROs. We suggest that at least 30%, and perhaps up to ~60%, of all EROs with R-K>5.3 and K<20.5 are dusty, star-forming systems at z>1. The SFD in this optically faint (R>26) population is rho^* (0.1-100Mo)=0.11+/-0.03 Mo/yr/Mpc^3, comparable to that in H-alpha emitting galaxies at z~1, and greater than that in UV-selected samples at these epochs. This support the claim of a strong increase in obscured star formation at high redshifts. Using the observed counts of the radio-detected EROs we model the break in the K-band number counts of all EROs at K~19-20 and propose that the passive ERO class dominates the total population in a narrow range around K~20, with dusty EROs dominating at fainter magnitudes. [Abridged]

We have identified counterparts to two submillimeter (submm) sources, SMM J09429+4659 and SMM J09431+4700, seen through the core of the z=0.41 cluster Abell 851. We employ deep 1.4-GHz observations and the far-infrared/radio correlation to refine the submm positions and then optical and near-infrared imaging to locate their counterparts. We identify an extremely red counterpart to SMM J09429+4659, while GMOS spectroscopy with Gemini-North shows that the R=23.8 radio source identified with SMM J09431+4700 is a hyperluminous infrared galaxy (L_FIR~1.5x10^13 L_sun) at z=3.35, the highest spectroscopic redshift so far for a galaxy discovered in the submm. The emission line properties of this galaxy are characteristic of a narrow-line Seyfert-1, although the lack of detected X-ray emission in a deep XMM-Newton observation suggests that the bulk of the luminosity of this galaxy is derived from massive star formation. We suggest that active nuclei, and the outflows they engender, may be an important part of the evolution of the brightest submm galaxies at high redshifts.

We present interplanetary network localization, spectral, and time history information for four short-duration, hard spectrum gamma-ray bursts, GRB000607, 001025B, 001204, and 010119. All of these events were followed up with sensitive radio and optical observations (the first and only such bursts to be followed up in the radio to date), but no detections were made, demonstrating that the short bursts do not have anomalously intense afterglows. We discuss the upper limits, and show that the lack of observable counterparts is consistent both with the hypothesis that the afterglow behavior of the short bursts is like that of the long duration bursts, many of which similarly have no detectable afterglows, as well as with the hypothesis that the short bursts have no detectable afterglows at all. Small number statistics do not allow a clear choice between these alternatives, but given the present detection rates of various missions, we show that progress can be expected in the near future.

Of all the well localized gamma-ray bursts, GRB 000911 has the longest duration (T_90 ~ 500 s), and ranks in the top 1% of BATSE bursts for fluence. Here, we report the discovery of the afterglow of this unique burst. In order to simultaneously fit our radio and optical observations, we are required to invoke a model involving an hard electron distribution, p ~ 1.5 and a jet-break time less than 1.5 day. A spectrum of the host galaxy taken 111 days after the burst reveals a single emission line, interpreted as [OII] at a redshift z = 1.0585, and a continuum break which we interpret as the Balmer limit at this redshift. Despite the long T_90, the afterglow of GRB 000911 is not unusual in any other way when compared to the set of afterglows studied to date. We conclude that the duration of the GRB plays little part in determining the physics of the afterglow.

We present multi-color light-curves of the optical afterglow of GRB 000926. Beginning ~1.5 days after the burst, the light-curves of this GRB steepen measurably. The existence of such achromatic breaks are usually taken to be an important observational signature that the ejecta are not expanding isotropically, but rather have a collimated jet-like geometry. If we interpret the data in this context, we derive an opening angle of 5 degrees, which reduces the energy release compared to an isotropic model by a factor of ~275, to 1.7 x 10^51 erg. To fit the data with a simple jet model requires extinction along the line of sight. The derived A_V is in the range 0.11 - 0.82 mag, depending on the adopted extinction law and whether the electrons giving rise to the optical emission are undergoing synchrotron cooling or not. Since this is in excess of the expected extinction from our Galaxy, we attribute this to the GRB host. We note that this extinction is typical of a galactic disk, and therefore the event likely took place in the disk of its host.

We use deep ground-based imaging in the near-IR to search for counterparts to the luminous submm sources in the catalog of Smail et al (1998). For the majority of the submm sources the near-IR imaging supports the counterparts originally selected from deep optical images. However, in two cases (10% of the sample) we find a relatively bright near-IR source close to the submm position, sources that were unidentified in the deep HST and ground-based R-band images used in Smail et al (1998). We place limits on colours of these sources from deep high-resolution Keck II imaging and find they have 2-sigma limits of (I-K)>6.8 and (I-K)>6.0 respectively. Both sources thus class as EROs. Using the spectral properties of the submm source in the radio and submm we argue that these EROs are probably the source of the submm emission, rather than the bright spiral galaxies previously identified by Smail et al. (1998). From the surface density of these submm-bright EROs we suggest that this class accounts for the majority of the reddest members of the ERO population, in good agreement with the preliminary conclusions of pointed submm observations of individual EROs. We conclude that the most extreme EROs represent a population of dusty, ultraluminous galaxies at high redshifts; further study of these will provide insights into the nature of star formation in obscured galaxies in the early Universe. The identification of similar counterparts in blank field submm surveys will be extremely difficult owing to their faintness (K~20.5, I>26.5). Finally, we discuss the radio and submm properties of the two submm-bright EROs discovered here and suggest that both galaxies lie at z>2.

We present deep Hubble Space Telescope NICMOS near-infrared and WFPC2 optical imaging of a small region in the core of the distant rich cluster Cl0939+4713 (z=0.41). We compare the optical and near-infrared morphologies of cluster members and find apparent small-scale optical structures within the galaxies which are absent in the near-infrared. We conclude that strong dust obscuration is a common feature in the late-type galaxies in distant clusters. We then concentrate on a sample of ten faint radio galaxies lying within our NICMOS field and selected from a very deep 1.4-GHz VLA map of the cluster with a 1sigma flux limit of 9uJy. Using published data we focus on the spectral properties of the eight radio-selected cluster members and show that these comprise a large fraction of the post-starburst population in the cluster. The simplest interpretation of the radio emission from these galaxies is that they are currently forming massive stars, contradicting their classification as post-starburst systems based on the optical spectra. We suggest that this star formation is hidden from view in the optical by the same obscuring dust which is apparent in our comparison on the optical and near-infrared morphologies of these galaxies. We caution that even in the restframe optical the effects of dust cannot be ignored when comparing samples of distant galaxies to low-redshift systems, particularly if dust is as prevelant in distant galaxies as appears to be the case in our study.

We report a detailed optical study of the clusters Abell 2125 and 2645. Both clusters have z=0.25 and richness class 4, yet contrast strongly in blue fraction and radio galaxy population. We find 27 spectroscopically confirmed radio galaxies in the blue cluster, A2125 and only four in A2645. The excess radio population in A2125 occurs entirely at L(20cm) < 10^23 W/Hz, where one expects star-formation to be responsible for the radio emission. Most of the radio galaxies have optical properties consistent with types later than E/S0, but emission lines weaker than would be expected for the SFR's implied by the radio emission. Thus we suspect dust obscuration is important. The cluster- cluster merger in A2125 seems likely to play a part in these phenomena.

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.