Passive early-type galaxies dominate cluster cores at z $\lesssim$1.5. At higher redshift, cluster core galaxies are observed to have still on-going star-formation, fuelled by cold molecular gas. We measure the molecular gas reservoir of the central region around the radio-loud AGN in the cluster CARLA J1103+3449 at z=1.44 with NOEMA. The AGN synchrotron emission dominates the continuum emission at 94.48 GHz, and we measure its flux at the AGN position and at the position of two radio jets. Combining our measurements with published results over the range 4.71 GHz-94.5 GHz, we obtain a flat spectral index $\alpha = 0.14 \pm 0.03$ for the AGN core emission, and a steeper index $\alpha = 1.43 \pm 0.04$ and $\alpha = 1.15 \pm 0.04$ at positions close to the western and eastern lobe, respectively. The total spectral index is $\alpha = 0.92 \pm 0.02$ over the range 73.8 MHz-94.5 GHz. We detect two CO(2-1) emission lines, both blue-shifted with respect to the AGN. Their emission corresponds to two regions, ~17 kpc south-east and ~14 kpc south-west of the AGN, not associated with galaxies. In these two regions, we find a total massive molecular gas reservoir of $M_{gas}$ = 3.9 $\pm$ 0.4 $10^{10} M_{\odot}$, which dominates (~ 60%) the central total molecular gas reservoir. These results can be explained by massive cool gas flows in the center of the cluster. The AGN early-type host is not yet quenched; its star formation rate is consistent with being on the main sequence of star-forming galaxies in the field (SFR~30-140 $M_{\odot}$/yr), and the cluster core molecular gas reservoir is expected to feed the AGN and the host star-formation before quiescence. The other cluster confirmed members show star formation rates at ~2 $\sigma$ below the field main sequence at similar redshifts and do not have molecular gas masses larger than galaxies of similar stellar mass in the field.

The entire southern sky (Declination, $\delta <$ 30 deg) has been observed using the Murchison Widefield Array (MWA), which provides radio imaging of $\sim$2-arcmin resolution at low frequencies (72-231 MHz). This is the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, and we have previously used a combination of visual inspection, cross-checks against the literature, and internal matching to identify the 'brightest' radio-sources ($S_{\mathrm{151MHz}} >$ 4 Jy) in the extragalactic catalogue (Galactic latitude, $|b| >$ 10 deg). We refer to these 1,863 sources as the GLEAM 4-Jy (G4Jy) Sample, and use radio images (of $\leq$ 45-arcsec resolution), and multi-wavelength information, to assess their morphology and identify the galaxy that is hosting the radio emission (where appropriate). Details of how to access all of the overlays used for this work are available at https://github.com/svw26/G4Jy. Alongside this we conduct further checks against the literature, which we document in this paper for individual sources. Whilst the vast majority of the G4Jy Sample are active galactic nuclei with powerful radio-jets, we highlight that it also contains a nebula, two nearby, star-forming galaxies, a cluster relic, and a cluster halo. There are also three extended sources for which we are unable to infer the mechanism that gives rise to the low-frequency emission. In the G4Jy catalogue we provide mid-infrared identifications for 86% of the sources, and flag the remainder as: having an uncertain identification (129 sources), having a faint/uncharacterised mid-infrared host (126 sources), or it being inappropriate to specify a host (2 sources). For the subset of 129 sources, there is ambiguity concerning candidate host-galaxies, and this includes four sources (B0424$-$728, B0703$-$451, 3C 198, and 3C 403.1) where we question the existing identification.

The Murchison Widefield Array (MWA) has observed the entire southern sky (Declination, $\delta <$ 30 deg) at low radio-frequencies, over the range 72-231 MHz. These observations constitute the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, and we use the extragalactic catalogue (Galactic latitude, $|b| >$ 10 deg) to define the GLEAM 4-Jy (G4Jy) Sample. This is a complete sample of the 'brightest' radio-sources ($S_{\mathrm{151MHz}} >$ 4 Jy), the majority of which are active galactic nuclei with powerful radio-jets. Crucially, low-frequency observations allow the selection of such sources in an orientation-independent way (i.e. minimising the bias caused by Doppler boosting, inherent in high-frequency surveys). We then use higher-resolution radio images, and information at other wavelengths, to morphologically classify the brightest components in GLEAM. We also conduct cross-checks against the literature, and perform internal matching, in order to improve sample completeness (which is estimated to be $>$ 95.5%). This results in a catalogue of 1,863 sources, making the G4Jy Sample over 10 times larger than that of the revised Third Cambridge Catalogue of Radio Sources (3CRR; $S_{\mathrm{178MHz}} >$ 10.9 Jy). Of these G4Jy sources, 78 are resolved by the MWA (Phase-I) synthesised beam ($\sim$2 arcmin at 200 MHz), and we label 67% of the sample as 'single', 26% as 'double', 4% as 'triple', and 3% as having 'complex' morphology at $\sim$1 GHz (45-arcsec resolution). Alongside this, our value-added catalogue provides mid-infrared source associations (subject to 6-arcsec resolution at 3.4 micron) for the radio emission, as identified through visual inspection and thorough checks against the literature. As such, the G4Jy Sample can be used as a reliable training set for cross-identification via machine-learning algorithms. [Abstract abridged for arXiv submission.]

The search for emission from weakly interacting massive particle (WIMP) dark matter annihilation and decay has become a multi-pronged area of research not only targeting a diverse selection of astrophysical objects, but also taking advantage of the entire electromagnetic spectrum. The decay of WIMP particles into standard model particles has been suggested as a possible channel for synchrotron emission to be detected at low radio frequencies. Here, we present the stacking analysis of a sample of 33 dwarf spheroidal (dSph) galaxies with low-frequency (72 - 231 MHz) radio images from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. We produce radial surface brightness profiles of images centred upon each dSph galaxy with background radio sources masked. We remove ten fields from the stacking due to contamination from either poorly subtracted, bright radio sources or strong background gradients across the field. The remaining 23 dSph galaxies are stacked in an attempt to obtain a statistical detection of any WIMP-induced synchrotron emission in these systems. We find that the stacked radial brightness profile does not exhibit a statistically significant detection above the 95% confidence level of $\sim$1.5 mJy beam$^{-1}$. This novel technique shows the potential of using low-frequency radio images to constrain fundamental properties of particle dark matter.

We present a detailed analysis of the radio galaxy PKS 2250-351, a giant of 1.2 Mpc projected size, its host galaxy, and its environment. We use radio data from the Murchison Widefield Array, the upgraded Giant Metre-wavelength Radio Telescope, the Australian Square Kilometre Array Pathfinder, and the Australia Telescope Compact Array to model the jet power and age. Optical and infra-red data come from the Galaxy And Mass Assembly (GAMA) survey and provide information on the host galaxy and environment. GAMA spectroscopy confirms that PKS 2250-351 lies at z=0.2115 in the irregular, and likely unrelaxed, cluster Abell 3936. We find its host is a massive, `red and dead' elliptical galaxy with negligible star formation but with a highly obscured active galactic nucleus dominating the mid-infrared emission. Assuming it lies on the local M-sigma relation it has an Eddington accretion rate of lambda_EDD~0.014. We find that the lobe-derived jet power (a time-averaged measure) is an order of magnitude greater than the hotspot-derived jet power (an instantaneous measure). We propose that over the lifetime of the observed radio emission (~300 Myr) the accretion has switched from an inefficient advection dominated mode to a thin-disc efficient mode, consistent with the decrease in jet power. We also suggest that the asymmetric radio morphology is due to its environment, with the host of PKS 2250-351 lying to the west of the densest concentration of galaxies in Abell 3936.

We present a new image of the 9.0 GHz radio emission from the extended Chandra Deep Field South. A total of 181 hours of integration with the Australia Telescope Compact Array has resulted in a 0.276 square degree image with a median sensitivity of $\sim$20 $\mu$Jy/beam rms, for a synthesised beam of 4.0 $\times$ 1.3 arcsec. We present a catalogue of the 9.0 GHz radio sources, identifying 70 source components and 55 individual radio galaxies. Source counts derived from this sample are consistent with those reported in the literature. The observed source counts are also generally consistent with the source counts from simulations of the faint radio population. Using the wealth of multiwavelength data available for this region, we classify the faint 9 GHz population and find that 91% are radio loud AGN, 7% are radio quiet AGN and 2% are star forming galaxies. The 9.0 GHz radio sources were matched to 5.5 and 1.4 GHz sources in the literature and we find a significant fraction of flat or inverted spectrum sources, with 36% of the 9 GHz sources having $\alpha_{5.5GHz}^{9.0GHz}$ $>$ -0.3 (for $S \propto \nu^\alpha$). This flat or inverted population is not well reproduced by current simulations of radio source populations.

The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories.

We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observed at 960 MHz (4489 sources) and 1320 MHz (5954 sources) with a bandwidth of 192 MHz and beam sizes of 30.0"x30.0" and 16.3"x15.1", respectively. The median Root Mean Squared (RMS) noise values are 186$\mu$Jy beam$^{-1}$ (960 MHz) and 165$\mu$Jy beam$^{-1}$ (1320 MHz). To create point source catalogues, we use these two source lists, together with the previously published Molonglo Observatory Synthesis Telescope (MOST) and the Australia Telescope Compact Array (ATCA) point source catalogues to estimate spectral indices for the whole population of radio point sources found in the survey region. Combining our ASKAP catalogues with these radio continuum surveys, we found 7736 point-like sources in common over an area of 30 deg$^2$. In addition, we report the detection of two new, low surface brightness supernova remnant candidates in the SMC. The high sensitivity of the new ASKAP ESP survey also enabled us to detect the bright end of the SMC planetary nebula sample, with 22 out of 102 optically known planetary nebulae showing point-like radio continuum emission. Lastly, we present several morphologically interesting background radio galaxies.

We have observed the G23 field of the Galaxy And Mass Assembly (GAMA) survey using the Australian Square Kilometre Array Pathfinder (ASKAP) in its commissioning phase, to validate the performance of the telescope and to characterize the detected galaxy populations. This observation covers $\sim$48 deg$^2$ with synthesized beam of 32.7$^{\prime\prime}$ by 17.8$^{\prime\prime}$ at 936 MHz, and $\sim$39 deg$^2$ with synthesized beam of 15.8$^{\prime\prime}$ by 12.0$^{\prime\prime}$ at 1320 MHz. At both frequencies, the r.m.s. (root-mean-square) noise is $\sim$0.1 mJy/beam. We combine these radio observations with the GAMA galaxy data, which includes spectroscopy of galaxies that are i-band selected with a magnitude limit of 19.2. Wide-field Infrared Survey Explorer (WISE) infrared (IR) photometry is used to determine which galaxies host an active galactic nucleus (AGN). In properties including source counts, mass distributions, and IR vs. radio luminosity relation, the ASKAP detected radio sources behave as expected. Radio galaxies have higher stellar mass and luminosity in IR, optical and UV than other galaxies. We apply optical and IR AGN diagnostics and find that they disagree for $\sim$30% of the galaxies in our sample. We suggest possible causes for the disagreement. Some cases can be explained by optical extinction of the AGN, but for more than half of the cases we do not find a clear explanation. Radio sources are more likely ($\sim$6%) to have an AGN than radio quiet galaxies ($\sim$1%), but the majority of AGN are not detected in radio at this sensitivity.

Future radio surveys will generate catalogues of tens of millions of radio sources, for which redshift estimates will be essential to achieve many of the science goals. However, spectroscopic data will be available for only a small fraction of these sources, and in most cases even the optical and infrared photometry will be of limited quality. Furthermore, radio sources tend to be at higher redshift than most optical sources and so a significant fraction of radio sources hosts differ from those for which most photometric redshift templates are designed. We therefore need to develop new techniques for estimating the redshifts of radio sources. As a starting point in this process, we evaluate a number of machine-learning techniques for estimating redshift, together with a conventional template-fitting technique. We pay special attention to how the performance is affected by the incompleteness of the training sample and by sparseness of the parameter space or by limited availability of ancillary multi-wavelength data. As expected, we find that the quality of the photometric-redshift degrades as the quality of the photometry decreases, but that even with the limited quality of photometry available for all sky-surveys, useful redshift information is available for the majority of sources, particularly at low redshift. We find that a template-fitting technique performs best with high-quality and almost complete multi-band photometry, especially if radio sources that are also X-ray emitting are treated separately. When we reduced the quality of photometry to match that available for the EMU all-sky radio survey, the quality of the template-fitting degraded and became comparable to some of the machine learning methods. Machine learning techniques currently perform better at low redshift than at high redshift, because of incompleteness of the currently available training data at high redshifts.

The GaLactic and Extragalactic All-sky MWA survey (GLEAM) is a radio continuum survey at 72-231 MHz of the whole sky south of declination +30 deg, carried out with the Murchison Widefield Array (MWA). In this paper, we derive source counts from the GLEAM data at 200, 154, 118 and 88 MHz, to a flux density limit of 50, 80, 120 and 290 mJy respectively, correcting for ionospheric smearing, incompleteness and source blending. These counts are more accurate than other counts in the literature at similar frequencies as a result of the large area of sky covered and this survey's sensitivity to extended emission missed by other surveys. At S_154MHz > 0.5 Jy, there is no evidence of flattening in the average spectral index (alpha approx. -0.8 where S proportional to nu^alpha) towards the lower frequencies. We demonstrate that the SKA Design Study (SKADS) model by Wilman et al. (2008) significantly underpredicts the observed 154 MHz GLEAM counts, particularly at the bright end. Using deeper LOFAR counts and the SKADS model, we find that sidelobe confusion dominates the thermal noise and classical confusion at nu >~ 100 MHz due to both the limited CLEANing depth and undeconvolved sources outside the field-of-view. We show that we can approach the theoretical noise limit using a more efficient and automated CLEAN algorithm.

Powerful radio-galaxies feature heavily in our understanding of galaxy evolution. However, when it comes to studying their properties as a function of redshift and/or environment, the most-detailed studies tend to be limited by small-number statistics. During Focus Meeting 3, on "Radio Galaxies: Resolving the AGN phenomenon", SVW presented a new sample of nearly 2,000 of the brightest radio-sources in the southern sky (Dec. $<$ 30 deg). These were observed at low radio-frequencies as part of the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, which is a continuum survey conducted using the Murchison Widefield Array (MWA). This instrument is the precursor telescope for the low-frequency component of the Square Kilometre Array, and allows us to select radio galaxies in an orientation-independent way (i.e. minimising the bias caused by Doppler boosting, inherent in high-frequency surveys). Being brighter than 4 Jy at 151 MHz, we refer to these objects as the GLEAM 4-Jy (G4Jy) Sample. The G4Jy catalogue is close to being finalised, with SVW describing how multi-wavelength data have been used to determine the morphology of the radio emission, and identify the host galaxy. In addition, the MWA's excellent spectral-coverage and sensitivity to extended/diffuse emission were highlighted. Both of these aspects are important for understanding the physical mechanisms that take place within active galaxies, and how they interact with their environment.

We present 0.3" (band 6) and 1.5" (band 3) ALMA observations of the (sub)millimeter dust continuum emission for 25 radio galaxies at 1<z<5.2. Our survey reaches a rms flux density of ~50$\mu$Jy in band 6 and ~20$\mu$Jy in band 3. This is an order of magnitude deeper than single-dish 850 $\mu$m observations, and reaches fluxes where synchrotron and thermal dust emission are expected to be of the same order of magnitude. Combining our sensitive ALMA observations with radio data from ATCA, VLA, and IR photometry from Herschel and Spitzer, we have disentangled the synchrotron and thermal dust emission. We determine the star-formation rates (SFR) and AGN IR luminosities using our newly developed spectral energy distribution fitting code MrMoose. We find that synchrotron emission contributes substantially at ~1 mm. Through our sensitive flux limits and accounting for a contribution from synchrotron emission in the mm, we revise downward the median SFR by a factor of 7 compared to previous estimates based solely on Herschel and Spitzer data. The hosts of these radio-loud AGN appear predominantly below the main sequence of star-forming galaxies, indicating that the star formation in many of the host galaxies has been quenched. Future growth of the host galaxies without substantial black hole mass growth will be needed to bring these objects on the local relation between the supermassive black holes and their host galaxies. Given the mismatch in the timescales of any star formation that took place in the host galaxies and lifetime of the AGN, we hypothesize that a key role is played by star formation in depleting the gas before the action of the powerful radio jets quickly drives out the remaining gas. This positive feedback loop of efficient star formation rapidly consuming the gas coupled to the action of the radio jets in removing the residual gas is how massive galaxies are rapidly quenched.

We report spectroscopic results from our 40-orbit $Hubble~Space~Telescope$ slitless grism spectroscopy program observing the 20 densest CARLA candidate galaxy clusters at $1.4 < z < 2.8$. These candidate rich structures, among the richest and most distant known, were identified on the basis of $[3.6]-[4.5]$ color from a $408~$hr multi-cycle $Spitzer$ program targeting $420$ distant radio-loud AGN. We report the spectroscopic confirmation of $16$ distant structures at $1.4 < z < 2.8$ associated with the targeted powerful high-redshift radio-loud AGN. We also report the serendipitous discovery and spectroscopic confirmation of seven additional structures at $0.87 < z < 2.12$ not associated with the targeted radio-loud AGN. We find that $10^{10} - 10^{11}\,M_{\odot}$ member galaxies of our confirmed CARLA structures form significantly fewer stars than their field counterparts at all redshifts within $1.4\leq z\leq 2$. We also observe higher star-forming activity in the structure cores up to $z = 2$, finding similar trends as cluster surveys at slightly lower redshifts ($1.0 < z < 1.5$). By design, our efficient strategy of obtaining just two grism orbits per field only obtains spectroscopic confirmation of emission-line galaxies. Deeper spectroscopy will be required to study the population of evolved, massive galaxies in these (forming) clusters. Lacking multi-band coverage of the fields, we adopt a very conservative approach of calling all confirmations "structures", although we note that a number of features are consistent with some of them being bona fide galaxy clusters. Together, this survey represents a unique and large homogenous sample of spectroscopically confirmed structures at high redshifts, potentially more than doubling the census of confirmed, massive clusters at $z > 1.4$.

We present a new broadband radio frequency catalogue across 0.12 GHz $\le \nu \le$ 20 GHz created by combining data from the Murchison Widefield Array Commissioning Survey, the Australia Telescope 20 Gigahertz survey, and the literature. Our catalogue consists of 1285 sources limited by $S_{20\,\rm{GHz}} >$ 40 mJy at $5\sigma$, and contains flux density measurements (or estimates) and uncertainties at 0.074, 0.080, 0.119, 0.150, 0.180, 0.408, 0.843, 1.4, 4.8, 8.6, and 20 GHz. We fit a second-order polynomial in log-log space to the spectral energy distributions of all these sources in order to characterise their broadband emission. For the 994 sources which are well-described by a linear or quadratic model we present a new diagnostic plot arranging sources by the linear and curvature terms. We demonstrate the advantages of such a plot over the traditional radio colour-colour diagram. We also present astrophysical descriptions of the sources found in each segment of this new parameter space and discuss the utility of these plots in the upcoming era of large area, deep, broadband radio surveys.

Hybrid morphology radio sources are a rare type of radio galaxy that display different Fanaroff-Riley classes on opposite sides of their nuclei. To enhance the statistical analysis of hybrid morphology radio sources, we embarked on a large-scale search of these sources within the international citizen science project, Radio Galaxy Zoo (RGZ). Here, we present 25 new candidate hybrid morphology radio galaxies. Our selected candidates are moderate power radio galaxies (L_median = 4.7x10^{24} W/(Hz sr) at redshifts 0.14<z<1.0. Hosts of nine candidates have spectroscopic observations, of which six are classified as quasars, one as high- and two as low-excitation galaxies. Two candidate HyMoRS are giant (>1Mpc) radio galaxies, one resides at a centre of a galaxy cluster, and one is hosted by a rare green bean galaxy. Although the origin of the hybrid morphology radio galaxies is still unclear, this type of radio source starts depicting itself as a rather diverse class. We discuss hybrid radio morphology formation in terms of the radio source environment (nurture) and intrinsically occurring phenomena (nature; activity cessation and amplification), showing that these peculiar radio galaxies can be formed by both mechanisms. While high angular resolution follow-up observations are still necessary to confirm our candidates, we demonstrate the efficacy of the Radio Galaxy Zoo in the pre-selection of these sources from all-sky radio surveys, and report the reliability of citizen scientists in identifying and classifying complex radio sources.

Cross-identifying complex radio sources with optical or infra red (IR) counterparts in surveys such as the Australia Telescope Large Area Survey (ATLAS) has traditionally been performed manually. However, with new surveys from the Australian Square Kilometre Array Pathfinder (ASKAP) detecting many tens of million of radio sources such an approach is no longer feasible. This paper presents new software (LRPY - Likelihood Ratio in PYthon) to automate the process of cross-identifying radio sources with catalogues at other wavelengths. LRPY implements the Likelihood Ratio (LR) technique with a modification to account for two galaxies contributing to a sole measured radio component. We demonstrate LRPY by applying it to ATLAS DR3 and a {\it Spitzer}-based multi-wavelength fusion catalogue, identifying 3,848 matched sources via our LR-based selection criteria. A subset of 1987 sources have flux density values for all IRAC bands which allow us to use criteria to distinguish between active galactic nuclei (AGN) and star-forming galaxies (SFG). We find that 936 radio sources ($\approx47\,\%$) meet both of the Lacy and Stern AGN selection criteria. Of the matched sources, 295 have spectroscopic redshifts and we examine the radio to IR flux ratio vs redshift, proposing an AGN selection criterion below the Elvis radio-loud (RL) AGN limit for this dataset. Taking the union of all three AGN selection criteria we identify 956 as AGN ($\approx48\,\%$). From this dataset, we find a decreasing fraction of AGN with lower radio flux densities consistent with other results in the literature.

We have acquired radio continuum data between 70\,MHz and 48\,GHz for a sample of 19 southern starburst galaxies at moderate redshifts ($0.067 < z < 0.227$) with the aim of separating synchrotron and free-free emission components. Using a Bayesian framework we find the radio continuum is rarely characterised well by a single power law, instead often exhibiting low frequency turnovers below 500\,MHz, steepening at mid-to-high frequencies, and a flattening at high frequencies where free-free emission begins to dominate over the synchrotron emission. These higher order curvature components may be attributed to free-free absorption across multiple regions of star formation with varying optical depths. The decomposed synchrotron and free-free emission components in our sample of galaxies form strong correlations with the total-infrared bolometric luminosities. Finally, we find that without accounting for free-free absorption with turnovers between 90 to 500\,MHz the radio-continuum at low frequency ($\nu < 200$\,MHz) could be overestimated by upwards of a factor of twelve if a simple power law extrapolation is used from higher frequencies. The mean synchrotron spectral index of our sample is constrained to be $\alpha=-1.06$, which is steeper then the canonical value of $-0.8$ for normal galaxies. We suggest this may be caused by an intrinsically steeper cosmic ray distribution.

The MIGHTEE large survey project will survey four of the most well-studied extragalactic deep fields, totalling 20 square degrees to $\mu$Jy sensitivity at Giga-Hertz frequencies, as well as an ultra-deep image of a single ~1 square degree MeerKAT pointing. The observations will provide radio continuum, spectral line and polarisation information. As such, MIGHTEE, along with the excellent multi-wavelength data already available in these deep fields, will allow a range of science to be achieved. Specifically, MIGHTEE is designed to significantly enhance our understanding of, (i) the evolution of AGN and star-formation activity over cosmic time, as a function of stellar mass and environment, free of dust obscuration; (ii) the evolution of neutral hydrogen in the Universe and how this neutral gas eventually turns into stars after moving through the molecular phase, and how efficiently this can fuel AGN activity; (iii) the properties of cosmic magnetic fields and how they evolve in clusters, filaments and galaxies. MIGHTEE will reach similar depth to the planned SKA all-sky survey, and thus will provide a pilot to the cosmology experiments that will be carried out by the SKA over a much larger survey volume.

We present results for the first three years of OzDES, a six-year programme to obtain redshifts for objects in the Dark Energy Survey (DES) supernova fields using the 2dF fibre positioner and AAOmega spectrograph on the Anglo-Australian Telescope. OzDES is a multi-object spectroscopic survey targeting multiple types of targets at multiple epochs over a multi-year baseline, and is one of the first multi-object spectroscopic surveys to dynamically include transients into the target list soon after their discovery. At the end of three years, OzDES has spectroscopically confirmed almost 100 supernovae, and has measured redshifts for 17,000 objects, including the redshifts of 2,566 supernova hosts. We examine how our ability to measure redshifts for targets of various types depends on signal-to-noise, magnitude, and exposure time, finding that our redshift success rate increases significantly at a signal-to-noise of 2 to 3 per 1-Angstrom bin. We also find that the change in signal-to-noise with exposure time closely matches the Poisson limit for stacked exposures as long as 10 hours. We use these results to predict the redshift yield of the full OzDES survey, as well as the potential yields of future surveys on other facilities such as the 4m Multi-Object Spectroscopic Telescope (4MOST), the Subaru Prime Focus Spectrograph (PFS), and the Maunakea Spectroscopic Explorer (MSE). This work marks the first OzDES data release, comprising 14,693 redshifts. OzDES is on target to obtain over a yield of approximately 5,700 supernova host-galaxy redshifts.

It is not yet known if the properties of molecular gas in distant protocluster galaxies are significantly affected by their environment as galaxies are in local clusters. Through a deep, 64 hours of effective on-source integration with the Australian Telescope Compact Array (ATCA), we discovered a massive, M_mol=2.0+-0.2x10^11 M_sun, extended, ~40 kpc, CO(1-0)-emitting disk in the protocluster surrounding the radio galaxy, MRC1138-262. The galaxy, at z_CO=2.1478, is a clumpy, massive disk galaxy, M_star~5x10^11 M_sun, which lies 250 kpc in projection from MRC1138-262 and is a known H-alpha emitter, HAE229. HAE229 has a molecular gas fraction of ~30%. The CO emission has a kinematic gradient along its major axis, centered on the highest surface brightness rest-frame optical emission, consistent with HAE229 being a rotating disk. Surprisingly, a significant fraction of the CO emission lies outside of the UV/optical emission. In spite of this, HAE229 follows the same relation between star-formation rate and molecular gas mass as normal field galaxies. HAE229 is the first CO(1-0) detection of an ordinary, star-forming galaxy in a protocluster. We compare a sample of cluster members at z>0.4 that are detected in low-order CO transitions with a similar sample of sources drawn from the field. We confirm findings that the CO-luminosity and FWHM are correlated in starbursts and show that this relation is valid for normal high-z galaxies as well as those in overdensities. We do not find a clear dichotomy in the integrated Schmidt-Kennicutt relation for protocluster and field galaxies. Not finding any environmental dependence in the "star-formation efficiency" or the molecular gas content, especially for such an extended CO disk, suggests that environmentally-specific processes such as ram pressure stripping are not operating efficiently in (proto)clusters. (abridged)

We examine the behaviour of the infrared-radio correlation (IRRC) over the range $0<z<6$ using new, highly sensitive 3GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared data from the Herschel Space Observatory in the 2deg$^{2}$ COSMOS field. We distinguish between objects where emission is believed to arise solely from star-formation, and those where an active galactic nucleus (AGN) is thought to be present. We account for non-detections in the radio or in the infrared using a doubly-censored survival analysis. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4GHz radio luminosity ratio ($q_{\rm TIR}$), decreases with increasing redshift: $q_{\rm TIR}(z)=(2.88\pm0.03)(1+z)^{-0.19\pm0.01}$. Moderate-to-high radiative luminosity AGN do not follow the same $q_{\rm TIR}$$(z)$ trend, having a lower normalisation and steeper decrease with redshift. We cannot rule out the possibility that unidentified AGN contributions only to the radio regime may be steepening the observed $q_{\rm TIR}(z)$ trend of the star-forming population. An increasing fractional contribution to the observed 3GHz flux by free-free emission of star-forming galaxies may also affect the derived evolution. However, we find that the standard (M82-based) assumption of the typical radio spectral energy distribution (SED) for star-forming galaxies is inconsistent with our results. This suggests a more complex shape of the typical radio SED for star-forming galaxies, and that imperfect $K$ corrections in the radio may govern the derived redshift trend of $q_{\rm TIR}$. Lastly, we present a redshift-dependent relation between rest-frame 1.4GHz radio luminosity and star formation rate taking the derived redshift trend into account.

We present new radio continuum observations of NGC253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of central starburst and extended emission. The central component, corresponding to the inner 500pc of the starburst region of the galaxy, is best modelled as an internally free-free absorbed synchrotron plasma, with a turnover frequency around 230 MHz. The extended emission component of the NGC253 spectrum is best described as a synchrotron emission flattening at low radio frequencies. We find that 34% of the extended emission (outside the central starburst region) at 1 GHz becomes partially absorbed at low radio frequencies. Most of this flattening occurs in the western region of the SE halo, and may be indicative of synchrotron self-absorption of shock re-accelerated electrons or an intrinsic low-energy cut off of the electron distribution. Furthermore, we detect the large-scale synchrotron radio halo of NGC253 in our radio images. At 154 - 231 MHz the halo displays the well known X-shaped/horn-like structure, and extends out to ~8kpc in z-direction (from major axis).

We present a robust calibration of the 1.4GHz radio continuum star formation rate (SFR) using a combination of the Galaxy And Mass Assembly (GAMA) survey and the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. We identify individually detected 1.4GHz GAMA-FIRST sources and use a late-type, non-AGN, volume-limited sample from GAMA to produce stellar mass-selected samples. The latter are then combined to produce FIRST-stacked images. This extends the robust parametrisation of the 1.4GHz-SFR relation to faint luminosities. For both the individually detected galaxies and our stacked samples, we compare 1.4GHz luminosity to SFRs derived from GAMA to determine a new 1.4GHz luminosity-to-SFR relation with well constrained slope and normalisation. For the first time, we produce the radio SFR-M* relation over 2 decades in stellar mass, and find that our new calibration is robust, and produces a SFR-M* relation which is consistent with all other GAMA SFR methods. Finally, using our new 1.4GHz luminosity-to-SFR calibration we make predictions for the number of star-forming GAMA sources which are likely to be detected in the upcoming ASKAP surveys, EMU and DINGO.

We present the first results from our on-going Australia Telescope Compact Array survey of CO(1-0) in ALMA-identified submillimetre galaxies in the Extended Chandra Deep Field South. Strong detections of CO(1-0) emission from two submillimetre galaxies, ALESS 122.1 (z = 2.0232) and ALESS 67.1 (z = 2.1230), were obtained. We estimate gas masses of M_gas ~ 1.3 \times 10^{11} M_odot and M_gas ~ 1.0 \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively, adopting alpha_CO = 1.0. Dynamical mass estimates from the kinematics of the CO(1-0) line yields M_dyn (sin i)^2 = 2.1 +- 1.1 \times 10^{11} M_odot and (3.2 +- 0.9) \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively. This is consistent with the total baryonic mass estimates of these two systems. We examine star formation efficiency using the L_FIR versus L'_CO(1-0) relation for samples of local ULIRGs and LIRGs, and more distant star-forming galaxies, with CO(1-0) detections. We find some evidence of a shallower slope for ULIRGs and SMGs compared to less luminous systems, but a larger sample is required for definite conclusions. We determine gas-to-dust ratios of 170 +- 30 and 140 +- 30 for ALESS 122.1 and ALESS 67.1, respectively, showing ALESS 122.1 has an unusually large gas reservoir. By combining the 38.1 GHz continuum detection of ALESS 122.1 with 1.4 and 5.5 GHz data, we estimate that the free-free contribution to radio emission at 38.1 GHz is 34 +- 17 microJy, yielding a star formation rate (1400 +- 700 M_\odot yr^{-1}) consistent with that from the infrared luminosity.

The largest galaxies in the Universe reside in galaxy clusters. Using sensitive observations of carbon-monoxide, we show that the Spiderweb Galaxy -a massive galaxy in a distant protocluster- is forming from a large reservoir of molecular gas. Most of this molecular gas lies between the protocluster galaxies and has low velocity dispersion, indicating that it is part of an enriched inter-galactic medium. This may constitute the reservoir of gas that fuels the widespread star formation seen in earlier ultraviolet observations of the Spiderweb Galaxy. Our results support the notion that giant galaxies in clusters formed from extended regions of recycled gas at high redshift.

We compare the kinetic energy and momentum injection rates from intense star formation, bolometric AGN radiation, and radio jets with the kinetic energy and momentum observed in the warm ionized gas in 24 powerful radio galaxies at z~2. These galaxies are amongst our best candidates for being massive galaxies near the end of their active formation period, when intense star formation, quasar activity, and powerful radio jets all co-exist. All galaxies have VLT/SINFONI imaging spectroscopy of the rest-frame optical line emission, showing emission-line regions with large velocity offsets (up to 1500 km/s) and line widths (typically 800-1000 km/s) consistent with very turbulent, often outflowing gas. As part of the HeRGE sample, they also have FIR estimates of the star formation and quasar activity obtained with Herschel/PACS and SPIRE, which enables us to measure the relative energy and momentum release from each of the three main sources of feedback in massive, star-forming AGN host galaxies during their most rapid formation phase. We find that star formation falls short by factors 10-1000 of providing the energy and momentum necessary to power the observed gas kinematics. The obscured quasars in the nuclei of these galaxies provide enough energy and momentum in about half of the sample, however, only if these are transfered to the gas relatively efficiently. We compare with theoretical and observational constraints on the efficiency of the energy and momentum transfer from jet and AGN radiation, which advocates that the radio jet is the main driver of the gas kinematics.

Infrared-faint radio sources (IFRS) are a class of radio-loud (RL) active galactic nuclei (AGN) at high redshifts (z > 1.7) that are characterised by their relative infrared faintness, resulting in enormous radio-to-infrared flux density ratios of up to several thousand. We aim to test the hypothesis that IFRS are young AGN, particularly GHz peaked-spectrum (GPS) and compact steep-spectrum (CSS) sources that have a low frequency turnover. We use the rich radio data set available for the Australia Telescope Large Area Survey fields, covering the frequency range between 150 MHz and 34 GHz with up to 19 wavebands from different telescopes, and build radio spectral energy distributions (SEDs) for 34 IFRS. We then study the radio properties of this class of object with respect to turnover, spectral index, and behaviour towards higher frequencies. We also present the highest-frequency radio observations of an IFRS, observed with the Plateau de Bure Interferometer at 105 GHz, and model the multi-wavelength and radio-far-infrared SED of this source. We find IFRS usually follow single power laws down to observed frequencies of around 150 MHz. Mostly, the radio SEDs are steep, but we also find ultra-steep SEDs. In particular, IFRS show statistically significantly steeper radio SEDs than the broader RL AGN population. Our analysis reveals that the fractions of GPS and CSS sources in the population of IFRS are consistent with the fractions in the broader RL AGN population. We find that at least 18% of IFRS contain young AGN, although the fraction might be significantly higher as suggested by the steep SEDs and the compact morphology of IFRS. The detailed multi-wavelength SED modelling of one IFRS shows that it is different from ordinary AGN, although it is consistent with a composite starburst-AGN model with a star formation rate of 170 solar masses per year.

High-redshift radio galaxies present signs of both star formation and AGN activity, making them ideal candidates to investigate the connection and coevolution of AGN and star formation in the progenitors of present-day massive galaxies. We make use of a sample of 11 powerful radio galaxies spanning 1<z<4 which have complete coverage of their spectral energy distribution (SED) from UV to FIR wavelengths. Using Herschel data, we disentangle the relative contribution of the AGN and star formation by combining the galaxyevolutioncodePEGASE.3 with an AGN torus model. We find that three components are necessary to reproduce the observed SEDs: an evolved and massive stellar component, a submm bright young starburst, and an AGN torus. We find that powerful radio galaxies form at very high-redshift, but experience episodic and important growth at 1<z<4 as the mass of the associated starburst varies from 5 to 50% of the total mass of the system. The properties of star formation differ from source to source, indicating no general trend of the star formation properties in the most infrared luminous high-redshift radio galaxies and no correlation with the AGN bolometric luminosity. Moreover, we find that AGN scattered light have a very limited impact on broad-band SED fitting on our sample. Finally, our analysis also suggests a wide range in origins for the observed star formation,which we partially constrain for some sources.

We have discovered a previously unreported poor cluster of galaxies (RGZ-CL J0823.2+0333) through an unusual giant wide-angle tail radio galaxy found in the Radio Galaxy Zoo project. We obtained a spectroscopic redshift of $z=0.0897$ for the E0-type host galaxy, 2MASX J08231289+0333016, leading to M$_r = -22.6$ and a $1.4\,$GHz radio luminosity density of $L_{\rm 1.4} = 5.5\times10^{24}$ W Hz$^{-1}$. These radio and optical luminosities are typical for wide-angle tailed radio galaxies near the borderline between Fanaroff-Riley (FR) classes I and II. The projected largest angular size of $\approx8\,$arcmin corresponds to $800\,$kpc and the full length of the source along the curved jets/trails is $1.1\,$Mpc in projection. X-ray data from the XMM-Newton archive yield an upper limit on the X-ray luminosity of the thermal emission surrounding RGZ J082312.9+033301,at $1.2-2.6\times10^{43}$ erg s$^{-1}$ for assumed intra-cluster medium temperatures of $1.0-5.0\,$keV. Our analysis of the environment surrounding RGZ J082312.9+033301 indicates that RGZ J082312.9+033301 lies within a poor cluster. The observed radio morphology suggests that (a) the host galaxy is moving at a significant velocity with respect to an ambient medium like that of at least a poor cluster, and that (b) the source may have had two ignition events of the active galactic nucleus with $10^7\,$yrs in between. This reinforces the idea that an association between RGZ J082312.9+033301, and the newly discovered poor cluster exists.

The GaLactic and Extragalactic All-sky MWA survey (GLEAM) is a new relatively low resolution, contiguous 72-231 MHz survey of the entire sky south of declination +25 deg. In this paper, we outline one approach to determine the relative contribution of system noise, classical confusion and sidelobe confusion in GLEAM images. An understanding of the noise and confusion properties of GLEAM is essential if we are to fully exploit GLEAM data and improve the design of future low-frequency surveys. Our early results indicate that sidelobe confusion dominates over the entire frequency range, implying that enhancements in data processing have the potential to further reduce the noise.

This paper outlines how the new GaLactic and Extragalactic All-sky MWA Survey (GLEAM, Wayth et al. 2015), observed by the Murchison Widefield Array covering the frequency range 72 - 231 MHz, allows identification of a new large, complete, sample of more than 2000 bright extragalactic radio sources selected at 151 MHz. With a flux density limit of 4 Jy this sample is significantly larger than the canonical fully-complete sample, 3CRR (Laing, Riley & Longair 1983). In analysing this small bright subset of the GLEAM survey we are also providing a first user check of the GLEAM catalogue ahead of its public release (Hurley-Walker et al. in prep). Whilst significant work remains to fully characterise our new bright source sample, in time it will provide important constraints to evolutionary behaviour, across a wide redshift and intrinsic radio power range, as well as being highly complementary to results from targeted, small area surveys.

We analyse a 154 MHz image made from a 12 h observation with the Murchison Widefield Array (MWA) to determine the noise contribution and behaviour of the source counts down to 30 mJy. The MWA image has a bandwidth of 30.72 MHz, a field-of-view within the half-power contour of the primary beam of 570 deg^2, a resolution of 2.3 arcmin and contains 13,458 sources above 5 sigma. The rms noise in the centre of the image is 4-5 mJy/beam. The MWA counts are in excellent agreement with counts from other instruments and are the most precise ever derived in the flux density range 30-200 mJy due to the sky area covered. Using the deepest available source count data, we find that the MWA image is affected by sidelobe confusion noise at the ~3.5 mJy/beam level, due to incompletely-peeled and out-of-image sources, and classical confusion becomes apparent at ~1.7 mJy/beam. This work highlights that (i) further improvements in ionospheric calibration and deconvolution imaging techniques would be required to probe to the classical confusion limit and (ii) the shape of low-frequency source counts, including any flattening towards lower flux densities, must be determined from deeper ~150 MHz surveys as it cannot be directly inferred from higher frequency data.

We present an intrinsic AGN SED extending from the optical to the submm, derived with a sample of unobscured, optically luminous (vLv(5100)>10^43.5 erg/s) QSOs at z<0.18 from the Palomar Green survey. The intrinsic AGN SED was computed by removing the contribution from stars using the 11.3um polycyclic aromatic hydrocarbon (PAH) feature in the QSOs' mid-IR spectra; the 1sigma uncertainty on the SED ranges between 12 and 45 per cent as a function of wavelength and is a combination of PAH flux measurement errors and the uncertainties related to the conversion between PAH luminosity and star-forming luminosity. Longwards of 20um the shape of the intrinsic AGN SED is independent of the AGN power indicating that our template should be applicable to all systems hosting luminous AGN (vLv(5100) or L_X(2-10keV) > 10^43.5 erg/s). We note that for our sample of luminous QSOs, the average AGN emission is at least as high as, and mostly higher than, the total stellar powered emission at all wavelengths from the optical to the submm. This implies that in many galaxies hosting powerful AGN, there is no `safe' broadband photometric observation (at lambda<1000um) which can be used in calculating star-formation rates without subtracting the AGN contribution. Roughly, the AGN contribution may be ignored only if the intrinsic AGN luminosity at 5100 Ang is at least a factor of 4 smaller than the total infrared luminosity (L_IR; 8-1000um) of the galaxy. Finally, we examine the implication of our work in statistical studies of star-formation in AGN host galaxies.

We present 0."5 resolution ALMA detections of the observed 246GHz continuum, [CI]^3P_2-^3P_1 fine structure line ([CI]2-1), CO(7-6) and H2O lines in the z=2.161 radio galaxy MRC1138-262, the 'Spiderweb Galaxy'. We detect strong [CI]2-1 emission both at the position of the radio core, and in a second component ~4kpc away from it. The 1100km/s broad [CI]2-1 line in this latter component, combined with its H2 mass of 1.6x10^10Msun implies this emission must come from a compact region <60pc, possibly containing a second AGN. The combined H2 mass derived for both objects using the [CI]2-1 emission is 3.3x10^10Msun. The total CO(7-6)/[CI]2-1 line flux ratio of 0.2 suggests a low excitation molecular gas reservoir and/or enhanced atomic carbon in cosmic-ray dominated regions. We detect spatially-resolved H2O 2_{11}-2_{02} emission - for the first time in a high-z un-lensed galaxy - near the outer radio lobe to the east, and near the bend of the radio jet to the west of the radio galaxy. No underlying 246GHz continuum emission is seen at either position. We suggest that the H2O emission is excited in the cooling region behind slow (10-40km/s) shocks in dense molecular gas (10^{3-5} cm^-3). The extended water emission is likely evidence of the radio jet's impact in cooling and forming molecules in the post-shocked gas in the halo and inter-cluster gas similar to what is seen in low-z clusters and other high-z radio galaxies. These observations imply that the passage of the radio jet in the interstellar and inter-cluster medium not only heats gas to high temperatures as is commonly assumed or found in simulations, but also induces cooling and dissipation which can lead to substantial amounts of cold dense molecular gas. The formation of molecules and strong dissipation in the halo gas of MRC1138-262 may explain both the extended diffuse molecular gas and young stars observed around MRC1138-262.

We used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02<z<0.5$, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500$\,\mu$m. SPIRE flux densities were also combined with Spitzer photometry and multi-wavelength archival data to perform a complete SED fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (8-1000$\,\mu$m) luminosity functions and their low-$z$ evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density (LLD) and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The luminosity functions show significant and rapid luminosity evolution already at low redshifts, $0.02<z<0.2$, with L$_{IR}^* \propto (1+z)^{6.0\pm0.4}$ and $\Phi_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $\Phi_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,\mu$m LFs respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter IMF and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate a SFRD scaling of SFRD$_0+0.08 z$, where SFRD$_0\simeq (1.9\pm 0.03)\times 10^{-2} [\mathrm{M}_\odot\,\mathrm{Mpc}^{-3}]$ is our total SFRD estimate at $z\sim0.02$.

We present a pilot study of the z=2.923 radio galaxy MRC0943-242, where we for the first time combine information from ALMA and MUSE data cubes. Even with modest integration times, we disentangle an AGN and a starburst dominated set of components. These data reveal a highly complex morphology, as the AGN, starburst, and molecular gas components show up as widely separated sources in dust continuum, optical continuum and CO line emission observations. CO(1-0) and CO(8-7) line emission suggest that there is a molecular gas reservoir offset from both the dust and the optical continuum that is located ~90kpc from the AGN. The UV line emission has a complex structure in emission and absorption. The line emission is mostly due to i) a large scale ionisation cone energised by the AGN, ii) a Ly-alpha emitting bridge of gas between the radio galaxy and a heavily star-forming set of components. Strangely, the ionisation cone has no Ly-alpha emission. We find this is due to an optically thick layer of neutral gas with unity covering fraction spread out over a region of at least ~100kpc from the AGN. Other, less thick absorption components are associated with Ly-alpha emitting gas within a few tens of kpc from the radio galaxy and are connected by a bridge of emission. We speculate that this linear structure of dust, Ly-alpha and CO emission, and the redshifted absorption seen in the circum-nuclear region may represent an accretion flow feeding gas into this massive AGN host galaxy.

The Dragonfly Galaxy (MRC0152-209), at redshift z~2, is one of the most vigorously star-forming radio galaxies in the Universe. What triggered its activity? We present ALMA Cycle 2 observations of cold molecular CO(6-5) gas and dust, which reveal that this is likely a gas-rich triple merger. It consists of a close double nucleus (separation ~4 kpc) and a weak CO-emitter at ~10 kpc distance, all of which have counterparts in HST/NICMOS imagery. The hyper-luminous starburst and powerful radio-AGN were triggered at this precoalescent stage of the merger. The CO(6-5) traces dense molecular gas in the central region, and complements existing CO(1-0) data, which revealed more widespread tidal debris of cold gas. We also find ~10$^{10}$ M(sun) of molecular gas with enhanced excitation at the highest velocities. At least 20-50% of this high-excitation, high-velocity gas shows kinematics that suggests it is being displaced and redistributed within the merger, although with line-of-sight velocities of |v| < 500 km/s, this gas will probably not escape the system. The processes that drive the redistribution of cold gas are likely related to either the gravitational interaction between two kpc-scale discs, or starburst/AGN-driven outflows. We estimate that the rate at which the molecular gas is redistributed is at least ~1200 +- 500 M(sun)/yr, and could perhaps even approach the star formation rate of ~3000 +- 800 M(sun)/yr. The fact that the gas depletion and gas redistribution timescales are similar implies that dynamical processes can be important in the evolution of massive high-z galaxies.

The Evolutionary Map of the Universe (EMU) is a proposed radio continuum survey of the Southern Hemisphere up to declination +30 deg., with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU will use an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. As part of the process of achieving this aim, a "Data Challenge" has been conducted, providing international teams the opportunity to test a variety of source finders on a set of simulated images. The aim is to quantify the accuracy of existing automated source finding and measurement approaches, and to identify potential limitations. The Challenge attracted nine independent teams, who tested eleven different source finding tools. In addition, the Challenge initiators also tested the current ASKAPsoft source-finding tool to establish how it could benefit from incorporating successful features of the other tools. Here we present the results of the Data Challenge, identifying the successes and limitations for this broad variety of the current generation of radio source finding tools. As expected, most finders demonstrate completeness levels close to 100% at 10sigma dropping to levels around 10% by 5sigma. The reliability is typically close to 100% at 10sigma, with performance to lower sensitivities varying greatly between finders. All finders demonstrate the usual trade-off between completeness and reliability, whereby maintaining a high completeness at low signal-to-noise comes at the expense of reduced reliability, and vice-versa. We conclude with a series of recommendations for improving the performance of the ASKAPsoft source-finding tool.

We present a new image of the 5.5 GHz radio emission from the extended Chandra Deep Field South. Deep radio observations at 5.5 GHz were obtained in 2010 and presented in the first data release. A further 76 hours of integration has since been obtained, nearly doubling the integration time. This paper presents a new analysis of all the data. The new image reaches 8.6 microJy rms, an improvement of about 40% in sensitivity. We present a new catalogue of 5.5 GHz sources, identifying 212 source components, roughly 50% more than were detected in the first data release. Source counts derived from this sample are consistent with those reported in the literature for S_{5.5GHz} > 0.1 mJy but significantly lower than published values in the lowest flux density bins (S_{5.5GHz} < 0.1 mJy), where we have more detected sources and improved statistical reliability. The 5.5 GHz radio sources were matched to 1.4 GHz sources in the literature and we find a mean spectral index of -0.35 +- 0.10 for S_{5.5GHz} > 0.5 mJy, consistent with the flattening of the spectral index observed in 5 GHz sub-mJy samples. The median spectral index of the whole sample is \alpha_{med} = -0.58, indicating that these observations may be starting to probe the star forming population. However, even at the faintest levels (0.05 < S_{5.5GHz} < 0.1 mJy), 39% of the 5.5 GHz sources have flat or inverted radio spectra. Four flux density measurements from our data, across the full 4.5 to 6.5 GHz bandwidth, are combined with those from literature and we find 10% of sources (S_{5.5GHz} >~ 0.1 mJy) show significant curvature in their radio spectral energy distribution spanning 1.4 to 9 GHz.

We present the third data release from the Australia Telescope Large Area Survey (ATLAS). These data combine the observations at 1.4 GHz before and after upgrades to the Australia Telescope Compact Array reaching a sensitivity of 14 microJy/beam in 3.6 deg^2 over the Chandra Deep Field South (CDFS) and of 17 microJy/beam in 2.7 deg^2 over the European Large Area ISO Survey South 1 (ELAIS-S1). We used a variety of array configurations to maximise the uv coverage resulting in a resolution of 16 by 7 arcsec in CDFS and of 12 by 8 arcsec in ELAIS-S1. After correcting for peak bias and bandwidth smearing, we find a total of 3034 radio source components above 5 sigma in CDFS, of which 514 (17 per cent) are considered to be extended. The number of components detected above 5 sigma in ELAIS-S1 is 2084, of which 392 (19 per cent) are classified as extended. The catalogues include reliable spectral indices (delta alpha < 0.2) between 1.40 and 1.71 GHz for ~350 of the brightest components.

We present results from the first twelve months of operation of Radio Galaxy Zoo, which upon completion will enable visual inspection of over 170,000 radio sources to determine the host galaxy of the radio emission and the radio morphology. Radio Galaxy Zoo uses $1.4\,$GHz radio images from both the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) and the Australia Telescope Large Area Survey (ATLAS) in combination with mid-infrared images at $3.4\,\mu$m from the {\it Wide-field Infrared Survey Explorer} (WISE) and at $3.6\,\mu$m from the {\it Spitzer Space Telescope}. We present the early analysis of the WISE mid-infrared colours of the host galaxies. For images in which there is $>\,75\%$ consensus among the Radio Galaxy Zoo cross-identifications, the project participants are as effective as the science experts at identifying the host galaxies. The majority of the identified host galaxies reside in the mid-infrared colour space dominated by elliptical galaxies, quasi-stellar objects (QSOs), and luminous infrared radio galaxies (LIRGs). We also find a distinct population of Radio Galaxy Zoo host galaxies residing in a redder mid-infrared colour space consisting of star-forming galaxies and/or dust-enhanced non star-forming galaxies consistent with a scenario of merger-driven active galactic nuclei (AGN) formation. The completion of the full Radio Galaxy Zoo project will measure the relative populations of these hosts as a function of radio morphology and power while providing an avenue for the identification of rare and extreme radio structures. Currently, we are investigating candidates for radio galaxies with extreme morphologies, such as giant radio galaxies, late-type host galaxies with extended radio emission, and hybrid morphology radio sources.

We separate the extragalactic radio source population above ~50 uJy into active galactic nuclei (AGN) and star-forming sources. The primary method of our approach is to fit the infrared spectral energy distributions (SEDs), constructed using Spitzer/IRAC and MIPS and Herschel/SPIRE photometry, of 380 radio sources in the Extended Chandra Deep Field-South. From the fitted SEDs, we determine the relative AGN and star-forming contributions to their infrared emission. With the inclusion of other AGN diagnostics such as X-ray luminosity, Spitzer/IRAC colours, radio spectral index and the ratio of star-forming total infrared flux to k-corrected 1.4 GHz flux density, qIR, we determine whether the radio emission in these sources is powered by star formation or by an AGN. The majority of these radio sources (60 per cent) show the signature of an AGN at some wavelength. Of the sources with AGN signatures, 58 per cent are hybrid systems for which the radio emission is being powered by star formation. This implies that radio sources which have likely been selected on their star formation have a high AGN fraction. Below a 1.4 GHz flux density of 1 mJy, along with finding a strong contribution to the source counts from pure star-forming sources, we find that hybrid sources constitute 20-65 per cent of the sources. This result suggests that hybrid sources have a significant contribution, along with sources that do not host a detectable AGN, to the observed flattening of the source counts at ~1mJy for the extragalactic radio source population.

We use a sample of 37 of the densest clusters and protoclusters across $1.3 \le z \le 3.2$ from the Clusters Around Radio-Loud AGN (CARLA) survey to study the formation of massive cluster galaxies. We use optical $i'$-band and infrared 3.6$\mu$m and 4.5$\mu$m images to statistically select sources within these protoclusters and measure their median observed colours; $\langle i'-[3.6] \rangle$. We find the abundance of massive galaxies within the protoclusters increases with decreasing redshift, suggesting these objects may form an evolutionary sequence, with the lower redshift clusters in the sample having similar properties to the descendants of the high redshift protoclusters. We find that the protocluster galaxies have an approximately unevolving observed-frame $i'-[3.6]$ colour across the examined redshift range. We compare the evolution of the $\langle i'-[3.6] \rangle$ colour of massive cluster galaxies with simplistic galaxy formation models. Taking the full cluster population into account, we show that the formation of stars within the majority of massive cluster galaxies occurs over at least 2Gyr, and peaks at $z \sim 2$-3. From the median $i'-[3.6]$ colours we cannot determine the star formation histories of individual galaxies, but their star formation must have been rapidly terminated to produce the observed red colours. Finally, we show that massive galaxies at $z>2$ must have assembled within 0.5Gyr of them forming a significant fraction of their stars. This means that few massive galaxies in $z>2$ protoclusters could have formed via dry mergers.

In the low-redshift Universe, the most powerful radio sources are often associated with gas-rich galaxy mergers or interactions. We here present evidence for an advanced, gas-rich (`wet') merger associated with a powerful radio galaxy at a redshift of z~2. This radio galaxy, MRC 0152-209, is the most infrared-luminous high-redshift radio galaxy known in the southern hemisphere. Using the Australia Telescope Compact Array, we obtained high-resolution CO(1-0) data of cold molecular gas, which we complement with HST/WFPC2 imaging and WHT long-slit spectroscopy. We find that, while roughly M(H2) ~ 2 x 10$^{10}$ M$_{\odot}$ of molecular gas coincides with the central host galaxy, another M(H2) ~ 3 x 10$^{10}$ M$_{\odot}$ is spread across a total extent of ~60 kpc. Most of this widespread CO(1-0) appears to follow prominent tidal features visible in the rest-frame near-UV HST/WFPC2 imaging. Ly$\alpha$ emission shows an excess over HeII, but a deficiency over L(IR), which is likely the result of photo-ionisation by enhanced but very obscured star formation that was triggered by the merger. In terms of feedback, the radio source is aligned with widespread CO(1-0) emission, which suggests that there is a physical link between the propagating radio jets and the presence of cold molecular gas on scales of the galaxy's halo. Its optical appearance, combined with the transformational stage at which we witness the evolution of MRC 0152-209, leads us to adopt the name `Dragonfly Galaxy'.

This is the collection of Chapters on 'Continuum Science' presented at the Conference 'Advancing Astrophysics with the SKA' (AASKA14), held in Giardini Naxos (Italy), June 9-13, 2014

In this chapter we provide an overview of the science enabled by radio continuum surveys in the SKA era, focusing on galaxy/galaxy cluster physics and evolution studies, and other relevant continuum science in the >2020 scientific framework. We outline a number of 'reference' radio-continuum surveys for SKA1 that can address such topics, and comprehensively discuss the most critical science requirements that we have identified. We highlight what should be achieved by SKA1, to guarantee a major leap forwards with respect to the pre-SKA era, considering the science advances expected in the coming years with existing and upcoming telescopes (JVLA, LOFAR, eMERLIN, and the three SKA precursors: MWA, ASKAP and MeerKAT). In this exercise we take in due account also the other waveband facilities coming online at the same time (e.g. Euclid, LSST, etc.), which tackle overlapping scientific goals, but in a different manner. In this respect particular attention has been payed to ensure that the proposed reference surveys are able to exploit the existing synergies with such facilities, so as to generate strong involvement from all astronomical communities, and leave a lasting legacy value. It is clear that a certain degree of freedom is allowed to some of the observational parameters. We believe it is very important to best fine-tune such parameters taking into proper account existing commensalities with SKA1 surveys addressing other science areas (HI galaxy science, magnetism, cosmology).

OzDES is a five-year, 100-night, spectroscopic survey on the Anglo-Australian Telescope, whose primary aim is to measure redshifts of approximately 2,500 Type Ia supernovae host galaxies over the redshift range 0.1 < z < 1.2, and derive reverberation-mapped black hole masses for approximately 500 active galactic nuclei and quasars over 0.3 < z < 4.5. This treasure trove of data forms a major part of the spectroscopic follow-up for the Dark Energy Survey for which we are also targeting cluster galaxies, radio galaxies, strong lenses, and unidentified transients, as well as measuring luminous red galaxies and emission line galaxies to help calibrate photometric redshifts. Here we present an overview of the OzDES program and our first-year results. Between Dec 2012 and Dec 2013, we observed over 10,000 objects and measured more than 6,000 redshifts. Our strategy of retargeting faint objects across many observing runs has allowed us to measure redshifts for galaxies as faint as m_r=25 mag. We outline our target selection and observing strategy, quantify the redshift success rate for different types of targets, and discuss the implications for our main science goals. Finally, we highlight a few interesting objects as examples of the fortuitous yet not totally unexpected discoveries that can come from such a large spectroscopic survey.

In this Chapter we present the motivation for undertaking both a wide and deep survey with the SKA in the context of studying AGN activity across cosmic time. With an rms down to 1 $\mu$Jy/beam at 1 GHz over 1,000 - 5,000 deg$^2$ in 1 year (wide tier band 1/2) and an rms down to 200 nJy/beam over 10 - 30 deg$^2$ in 2000 hours (deep tier band 1/2), these surveys will directly detect faint radio-loud and radio-quiet AGN (down to a 1 GHz radio luminosity of about $2\times10^{23}$ W/Hz at $z=6$). For the first time, this will enable us to conduct detailed studies of the cosmic evolution of radio AGN activity to the cosmic dawn ($z\gtrsim6$), covering all environmental densities.

Faraday rotation of polarised background sources is a unique probe of astrophysical magnetic fields in a diverse range of foreground objects. However, to understand the properties of the polarised sources themselves and of depolarising phenomena along the line of sight, we need to complement Faraday rotation data with polarisation observations over very broad bandwidths. Just as it is impossible to properly image a complex source with limited u-v coverage, we can only meaningfully understand the magneto-ionic properties of polarised sources if we have excellent coverage in $\lambda^2$-space. We here propose a set of broadband polarisation surveys with the Square Kilometre Array, which will provide a singular set of scientific insights on the ways in which galaxies and their environments have evolved over cosmic time.

In this chapter, we highlight a number of science investigations that are enabled by the inclusion of Band~5 ($4.6-13.8$ GHz) for SKA1-MID science operations, while focusing on the astrophysics of star formation over cosmic time. For studying the detailed astrophysics of star formation at high-redshift, surveys at frequencies $\gtrsim$10 GHz have the distinct advantage over traditional $\sim$1.4 GHz surveys as they are able to yield higher angular resolution imaging while probing higher rest frame frequencies of galaxies with increasing redshift, where emission of star-forming galaxies becomes dominated by thermal (free-free) radiation. In doing so, surveys carried out at $\gtrsim$10 GHz provide a robust, dust-unbiased measurement of the massive star formation rate by being highly sensitive to the number of ionizing photons that are produced. To access this powerful star formation rate diagnostic requires that Band~5 be available for SKA1-MID. We additionally present a detailed science case for frequency coverage extending up to 30 GHz during full SKA2 operations, as this allows for highly diverse science while additionally providing contiguous frequency coverage between the SKA and ALMA, which will likely be the two most powerful interferometers for the coming decades. To enable this synergy, it is crucial that the dish design of the SKA be flexible enough to include the possibility of being fit with receivers operating up to 30 GHz.

We show that, in addition to specific science goals, there is a strong case for conducting an all-sky (i.e. the visible 3-pi steradians) SKA continuum survey which does not fit neatly into conventional science cases. History shows that the greatest scientific impact of most major telescopes (e.g., HST, VLA) lies beyond the original goals used to justify the telescope. The design of the telescope therefore needs to maximise the ultimate scientific productivity, in addition to achieving the specific science goals. In this chapter, we show that an all-sky continuum survey is likely to achieve transformational science in two specific respects: (1) Discovering the unexpected (2) Transforming radio-astronomy from niche to mainstream

Astronomy at or below the 'survey threshold' has expanded significantly since the publication of the original 'Science with the Square Kilometer Array' in 1999 and its update in 2004. The techniques in this regime may be broadly (but far from exclusively) defined as 'confusion' or 'P(D)' analyses (analyses of one-point statistics), and 'stacking', accounting for the flux-density distribution of noise-limited images co-added at the positions of objects detected/isolated in a different waveband. Here we discuss the relevant issues, present some examples of recent analyses, and consider some of the consequences for the design and use of surveys with the SKA and its pathfinders.

Radio wavelengths offer the unique possibility of tracing the total star-formation rate in galaxies, both obscured and unobscured. As such, they may provide the most robust measurement of the star-formation history of the Universe. In this chapter we highlight the constraints that the SKA can place on the evolution of the star-formation history of the Universe, the survey area required to overcome sample variance, the spatial resolution requirements, along with the multi-wavelength ancillary data that will play a major role in maximising the scientific promise of the SKA. The required combination of depth and resolution means that a survey to trace the star formation in the Universe should be carried out with a facility that has a resolution of at least ~0.5arcsec, with high sensitivity at < 1 GHz. We also suggest a strategy that will enable new parameter space to be explored as the SKA expands over the coming decade.

It has become apparent that active galactic nuclei (AGN) may have a significant impact on the growth and evolution of their host galaxies and vice versa but a detailed understanding of the interplay between these processes remains elusive. Deep radio surveys provide a powerful, obscuration-independent tool for measuring both star formation and AGN activity in high-redshift galaxies. Multiwavelength studies of deep radio fields show a composite population of star-forming galaxies and AGN, with the former dominating at the lowest flux densities (S$_{1.4\mathrm{GHz}}<$100~$\mu$Jy). The sensitivity and resolution of the SKA will allow us to identify, and separately trace, the total star formation in the bulges of individual high-redshift galaxies, the related nuclear activity and any star formation occurring on larger scales within a disc. We will therefore gain a detailed picture of the apparently simultaneous development of stellar populations and black holes in the redshift range where both star-formation and AGN activity peak (1$\leq$z$\leq$4). In this chapter we discuss the role of the SKA in studying the connection between AGN activity and galaxy evolution, and the most critical technical requirements for such of studies

One of the most challenging and exciting subjects in modern astrophysics is that of galaxy formation at the epoch of reionisation. The SKA, with its revolutionary capabilities in terms of frequency range, resolution and sensitivity, will allow to explore the first Gyr of structure formation in the Universe, in particular, with the detection and study of the earliest manifestations of the AGN phenomenon. The tens of QSOs that are currently known out to the highest redshifts (z~7), many of them exhibiting powerful radio emission, imply that super-massive black holes can be grown on a very short timescale and support the existence of very high redshift (z > 7) radio loud sources - sources that have so far escaped detection. Not only would such detections be paramount to the understanding of the earliest stages of galaxy evolution, they are necessary for the direct study of neutral hydrogen in the Epoch of Reionisation, through observations of the HI 21cm forest against such background sources. In order to understand how SKA and SKA1 observations can be optimised to reveal these earliest AGN, we have examined the effect of a hot CMB on the emission of powerful and young radio galaxies. By looking at the SKA1 capabilities, in particular in terms of wavelength coverage and resolution, we determine how the effects of "CMB-muting" of a radio loud source can be observationally minimised and how to identify the best highest-redshift radio candidates. Considering different predictions for the space density of radio loud AGN at such redshifts, we identify the survey characteristics necessary to optimize the detection and identification of the very first generation of radio loud AGN in the Universe.

We present APEX LABOCA 870 micron observations of the field around the high-redshift radio galaxy MRC1138-262 at z=2.16. We detect 16 submillimeter galaxies in this ~140 square arcmin bolometer map with flux densities in the range 3-11 mJy. The raw number counts indicate a density of submillimeter galaxies (SMGs) that is up to four times that of blank field surveys. Based on an exquisite multiwavelength database, including VLA 1.4 GHz radio and infrared observations, we investigate whether these sources are members of the protocluster structure at z=2.2. Using Herschel PACS+SPIRE and Spitzer MIPS photometry, we derive reliable far-infrared photometric redshifts for all sources. Follow-up VLT ISAAC and SINFONI near-infrared spectra confirm that four of these SMGs have redshifts of z=2.2. We also present evidence that another SMG in this field, detected earlier at 850 micron, has a counterpart that exhibits Halpha and CO(1-0) emission at z=2.15. Including the radio galaxy and two SMGs with far-IR photometric redshifts at z=2.2, we conclude that at least eight submm sources are part of the protocluster at z=2.16 associated with the radio galaxy MRC1138-262. We measure a star formation rate density SFRD ~1500 Msun yr^-1 Mpc^-3, four magnitudes higher than the global SFRD of blank fields at this redshift. Strikingly, these eight sources are concentrated within a region of 2 Mpc (the typical size of clusters in the local universe) and are distributed within the filaments traced by the Halpha emitters at z=2.2. This concentration of massive, dusty starbursts is not centered on the submillimeter-bright radio galaxy which could support the infalling of these sources into the cluster center. Approximately half (6/11) of the SMGs that are covered by the Halpha imaging data are associated with Halpha emitters, demonstrating the potential of tracing SMG counterparts with this population (abridged).

Distant powerful radio-loud active galactic nuclei (RLAGN) tend to reside in dense environments and are commonly found in proto-clusters at z > 1.3. We examine whether this occurs because RLAGN are hosted by massive galaxies, which preferentially reside in rich environments. We compare the environments of powerful RLAGN at 1.3 < z < 3.2 from the CARLA survey to a sample of radio-quiet galaxies matched in mass and redshift. We find the environments of RLAGN are significantly denser than those of radio-quiet galaxies, implying that not more than 50% of massive galaxies in this epoch can host powerful radio-loud jets. This is not an observational selection effect as we find no evidence to suggest it is easier to observe the radio emission when the galaxy resides in a dense environment. We therefore suggest that the dense Mpc-scale environment fosters the formation of a radio-jet from an AGN. We show that the number density of potential RLAGN host galaxies is consistent with every > 10^14 solar mass cluster having experienced powerful radio-loud feedback of duration ~60 Myr during 1.3 < z < 3.2. This feedback could heat the intracluster medium to the extent of 0.5-1 keV per gas particle, which could limit the amount of gas available for further star formation in the proto-cluster galaxies.

We re-examine a series of archived centimetre radio-continuum observations (lambda=16, 6 and 3 cm) focusing on NGC7793 using the Australia Telescope Compact Array. These new images are both very sensitive (rms=0.011 mJy/beam) and feature reasonably high angular resolution (down to < 3"). Using these images, a total of 76 discrete radio sources are identified, of which 57 have been classified. We also studied the radio component of the micro-quasar NGC7793-S26 which shows two distinct regions of somewhat steep spectral index between -0.3 and -0.7.

We present Herschel/PACS observations of the nearby (z=0.1055) dwarf galaxy that has hosted the long gamma ray burst (LGRB) 031203. Using the PACS data we have been able to place constraints on the dust temperature, dust mass, total infrared luminosity and infrared-derived star-formation rate (SFR) for this object. We find that the GRB host galaxy (GRBH) 031203 has a total infrared luminosity of 3x10^10 L_sun placing it in the regime of the IR-luminous galaxy population. Its dust temperature and specific SFR are comparable to that of many high-redshift (z=0.3-2.5) infrared (IR)-detected GRB hosts (T_dust>40K ; sSFR>10 Gyr^-1), however its dust-to-stellar mass ratio is lower than what is commonly seen in IR-luminous galaxies. Our results suggest that GRBH 031203 is undergoing a strong starburst episode and its dust properties are different to those of local dwarf galaxies within the same metallicity and stellar mass range. Furthermore, our measurements place it in a distinct class to the well studied nearby host of GRB 980425 (z=0.0085), confirming the notion that GRB host galaxies can span a large range in properties even at similar cosmological epochs, making LGRBs an ideal tool in selecting samples of star-forming galaxies up to high redshift.

We present results from a survey of 70 radio galaxies (RGs) at redshifts 1<z<5.2 using the PACS and SPIRE on-board Herschel. Combined with existing mid-IR photometry from Spitzer and observations obtained with LABOCA, the SEDs of galaxies in our sample are continuously covered across 3.6-870um. The total infrared luminosities of these RGs are such that they almost all are either ultra-or hyper-luminous infrared galaxies. We fit the infrared SEDs with a set of empirical templates which represent dust heated (1) by a variety of SB and (2) by a AGN. We find that the SEDs of RGs require the dust to be heated by both AGN and SB, but the luminosities of these two components are not strongly correlated. Assuming empirical relations and simple physical assumptions, we calculate the SFR, the black hole mass accretion rate (MdotBH), and the black hole mass (MBH) for each RG. We find that the host galaxies and their BHs are growing extremely rapidly, having SFR~100-5000 Msun/yr and MdotBH~1-100 Msun/yr. The mean sSFR of RGs at z>2.5 are higher than the sSFR of typical star-forming galaxies over the same redshift range but are similar or perhaps lower than the galaxy population for RGs at z<2.5. By comparing the sSFR and the specific black hole mass accretion rate, we conclude that BHs in radio loud AGN are already, or soon will be, overly massive compared to their host galaxies in terms of expectations from the local MBH-MGal relation. In order to ``catch up'' with the BH, the galaxies require about an order-of magnitude more time to grow in mass, at the observed SFRs, compared to the time the BH is actively accreting. However, during the current cycle of activity, we argue that this catching-up is likely to be difficult due to the short gas depletion times. Finally, we speculate on how the host galaxies might grow sufficiently in stellar mass to ultimately fall onto the local MBH-MGal relation.

We present 4.5 {\mu}m luminosity functions for galaxies identified in 178 candidate galaxy clusters at 1.3 < z < 3.2. The clusters were identified as Spitzer/IRAC color-selected overdensities in the Clusters Around Radio-Loud AGN (CARLA) project, which imaged 421 powerful radio-loud AGN at z > 1.3. The luminosity functions are derived for different redshift and richness bins, and the IRAC imaging reaches depths of m*+2, allowing us to measure the faint end slopes of the luminosity functions. We find that {\alpha} = -1 describes the luminosity function very well in all redshifts bins and does not evolve significantly. This provides evidence that the rate at which the low mass galaxy population grows through star formation, gets quenched and is replenished by in-falling field galaxies does not have a major net effect on the shape of the luminosity function. Our measurements for m* are consistent with passive evolution models and high formation redshifts z_f ~ 3. We find a slight trend towards fainter m* for the richest clusters, implying that the most massive clusters in our sample could contain older stellar populations, yet another example of cosmic downsizing. Modelling shows that a contribution of a star-forming population of up to 40% cannot be ruled out. This value, found from our targeted survey, is significantly lower than the values found for slightly lower redshift, z ~ 1, clusters found in wide-field surveys. The results are consistent with cosmic downsizing, as the clusters studied here were all found in the vicinity of radio-loud AGNs -- which have proven to be preferentially located in massive dark matter halos in the richest environments at high redshift -- they may therefore be older and more evolved systems than the general protocluster population.

We present a CO(1-0) survey for cold molecular gas in a representative sample of 13 high-z radio galaxies (HzRGs) at 1.4<z<2.8, using the Australia Telescope Compact Array. We detect CO(1-0) emission associated with five sources: MRC 0114-211, MRC 0152-209, MRC 0156-252, MRC 1138-262 and MRC 2048-272. The CO(1-0) luminosities are in the range $L'_{\rm CO} \sim (5 - 9) \times 10^{10}$ K km/s pc$^{2}$. For MRC 0152-209 and MRC 1138-262 part of the CO(1-0) emission coincides with the radio galaxy, while part is spread on scales of tens of kpc and likely associated with galaxy mergers. The molecular gas mass derived for these two systems is M$_{\rm H2} \sim 6 \times 10^{10}\, {\rm M}_{\odot}$ (M$_{\rm H2}$/$L'_{\rm CO}$=0.8). For the remaining three CO-detected sources, the CO(1-0) emission is located in the halo (~50-kpc) environment. These three HzRGs are among the fainter far-IR emitters in our sample, suggesting that similar reservoirs of cold molecular halo gas may have been missed in earlier studies due to pre-selection of IR-bright sources. In all three cases the CO(1-0) is aligned along the radio axis and found beyond the brightest radio hot-spot, in a region devoid of 4.5$\mu$m emission in Spitzer imaging. The CO(1-0) profiles are broad, with velocity widths of ~ 1000 - 3600 km/s. We discuss several possible scenarios to explain these halo reservoirs of CO(1-0). Following these results, we complement our CO(1-0) study with detections of extended CO from the literature and find at marginal statistical significance (95% level) that CO in HzRGs is preferentially aligned towards the radio jet axis. For the eight sources in which we do not detect CO(1-0), we set realistic upper limits of $L'_{\rm CO} \sim 3-4 \times 10^{10}$ K km/s pc$^{2}$. Our survey reveals a CO(1-0) detection rate of 38%, allowing us to compare the CO(1-0) content of HzRGs with that of other types of high-z galaxies.

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.

This paper presents the first results of a far-infrared search for protocluster-associated galaxy overdensities using the SPIRE instrument on-board the {\it Herschel} Space Observatory. Large ($\sim$400 arcmin$^{2}$) fields surrounding 26 powerful high-redshift radio galaxies ($2.0 < z < 4.1$; $L_{\rm 500 MHz} > 10^{28.5}$ WHz$^{-1}$) are mapped at 250, 350 and 500\mic to give a unique wide-field sample. On average the fields have a higher than expected, compared to blank fields, surface density of 500\mic sources within 6 comoving Mpc of the radio galaxy. The analysis is then restricted to potential protocluster members only, which are identified using a far-infrared colour selection; this reveals significant overdensities of galaxies in 2 fields, neither of which are previously known protoclusters. The probability of finding 2 overdensities of this size by chance, given the number of fields observed is $5 \times 10^{-4}$. Overdensities here exist around radio galaxies with $L_{\rm 500 MHz} \gtrsim 10^{29}$ WHz$^{-1}$ and $z < 3$. The radial extent of the average far-infrared overdensity is found to be $\sim$6 comoving Mpc. Comparison with predictions from numerical simulations shows that the overdensities are consistent with having masses $> 10^{14}$Msolar. However, the large uncertainty in the redshift estimation means that it is possible that these far-infrared overdensities consist of several structures across the redshift range searched.

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

We present the spectroscopic confirmation of a structure of galaxies surrounding the radio galaxy MRC0156-252 at z = 2.02. The structure was initially discovered as an overdensity of both near-infrared selected z > 1.6 and mid-infrared selected z > 1.2 galaxy candidates. We used the VLT/FORS2 multi-object spectrograph to target ~80 high-redshift galaxy candidates, and obtain robust spectroscopic redshifts for more than half the targets. The majority of the confirmed sources are star-forming galaxies at z > 1.5. In addition to the radio galaxy, two of its close-by companions (< 6'') also show AGN signatures. Ten sources, including the radio galaxy, lie within |z - 2.020 | < 0.015 (i.e., velocity offsets < 1500 km/s) and within projected 2 Mpc comoving of the radio galaxy. Additional evidence suggests not only that the galaxy structure associated with MRC0156-252 is a forming galaxy cluster but also that this structure is most probably embedded in a larger scale structure.

We present the results of a program of optical and near-infrared spectroscopic follow-up of candidate Active Galactic Nuclei (AGN) selected in the mid-infrared. This survey selects both normal and obscured AGN closely matched in luminosity across a wide range, from Seyfert galaxies with bolometric luminosities L_bol~10^10L_sun, to highly luminous quasars (L_bol~10^14L_sun), and with redshifts from 0-4.3. Samples of candidate AGN were selected through mid-infrared color cuts at several different 24 micron flux density limits to ensure a range of luminosities at a given redshift. The survey consists of 786 candidate AGN and quasars, of which 672 have spectroscopic redshifts and classifications. Of these, 137 (20%) are type-1 AGN with blue continua, 294 (44%) are type-2 objects with extinctions A_V>~5 towards their AGN, 96 (14%) are AGN with lower extinctions (A_V~1) and 145 (22%) have redshifts, but no clear signs of AGN activity in their spectra. 50% of the survey objects have L_bol >10^12L_sun, in the quasar regime. We present composite spectra for type-2 quasars and for objects with no signs of AGN activity in their spectra. We also discuss the mid-infrared - emission-line luminosity correlation and present the results of cross-correlations with serendipitous X-ray and radio sources. The results show that: (1) obscured objects dominate the overall AGN population, (2) there exist mid-infrared selected AGN candidates which lack AGN signatures in their optical spectra, but have AGN-like X-ray or radio counterparts, and (3) X-ray and optical classifications of obscured and unobscured AGN often differ.

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

We describe the far-infrared (FIR; rest-frame 8--1000\mu m) properties of a sample of 443 H\alpha-selected star-forming galaxies in the COSMOS and UDS fields detected by the HiZELS imaging survey. Sources are identified using narrow-band filters in combination with broad-band photometry to uniformly select H\alpha\ (and [OII] if available) emitters in a narrow redshift slice at z = 1.47+/-0.02. We use a stacking approach in Spitzer, Herschel (from PEP and HerMES surveys) and AzTEC images to describe their typical FIR properties. We find that HiZELS galaxies with observed H\alpha\ luminosities of ~ 10^{8.1-9.1} Lo have bolometric FIR luminosities of typical LIRGs, L_FIR ~ 10^{11.48+/-0.05} Lo. Combining the H\alpha\ and FIR luminosities, we derive median SFR = 32+/-5 Mo/yr and H\alpha\ extinctions of A(H\alpha) = 1.0+/-0.2 mag. Perhaps surprisingly, little difference is seen in typical HiZELS extinction levels compared to local star-forming galaxies. We confirm previous empirical stellar mass (M*) to A(H\alpha) relations and the little or no evolution up to z = 1.47. For HiZELS galaxies, we provide an empirical parametrisation of the SFR as a function of (u-z)_rest colours and 3.6\mu m photometry. We find that the observed H\alpha\ luminosity is a dominant SFR tracer when (u-z)_rest ~< 0.9 mag or when 3.6\mu m photometry > 22 mag (Vega) or when M* < 10^9.7 Mo. We do not find any correlation between the [OII]/H\alpha\ and FIR luminosity, suggesting that this emission line ratio does not trace the extinction of the most obscured star-forming regions. The luminosity-limited HiZELS sample tends to lie above of the so-called `main sequence' for star-forming galaxies, especially at low M*. This work suggests that obscured star formation is linked to the assembly of M*, with deeper potential wells in massive galaxies providing dense, heavily obscured environments in which stars can form rapidly.

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.

We report the first results from the Clusters Around Radio-Loud AGN (CARLA) program, a Cycle 7 and 8 Spitzer Space Telescope snapshot program to investigate the environments of a large sample of obscured and unobscured luminous radio-loud AGN at 1.2 < z < 3.2. These data, obtained for 387 fields, reach 3.6 and 4.5 um depths of [3.6] (AB) = 22.6 and [4.5] (AB) = 22.9 at the 95% completeness level, which is two to three times fainter than L* in this redshift range. By using the color cut [3.6]-[4.5] > -0.1 (AB), which efficiently selects high-redshift (z > 1.3) galaxies of all types, we identify galaxy cluster member candidates in the fields of the radio-loud AGN. The local density of these IRAC-selected sources is compared to the density of similarly selected sources in blank fields. We find that 92% of the radio-loud AGN reside in environments richer than average. The majority (55%) of the radio-loud AGN fields are found to be overdense at a > 2 {\sigma} level; 10% are overdense at a > 5 {\sigma} level. A clear rise in surface density of IRAC-selected sources towards the position of the radio-loud AGN strongly supports an association of the majority of the IRAC-selected sources with the radio-loud AGN. Our results provide solid statistical evidence that radio-loud AGN are likely beacons for finding high-redshift galaxy (proto-)clusters. We investigate how environment depends on AGN type (unobscured radio-loud quasars vs. obscured radio galaxies), radio luminosity and redshift, finding no correlation with either AGN type or radio luminosity. We find a decrease in density with redshift, consistent with galaxy evolution for this uniform, flux-limited survey. These results are consistent with expectations from the orientation-driven AGN unification model, at least for the high radio luminosity regimes considered in this sample.

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.1<z<2 and determine how these evolve with cosmic time. The unique angle of this work is the rigorous analysis of survey selection effects, making this the first study of the star-formation-dominated, IR-luminous population within a framework almost entirely free of selection biases. We find that IR-luminous galaxies have SEDs with broad far-IR peaks characterised by cool/extended dust emission and average dust temperatures in the 25-45K range. Hot (T>45K) SEDs and cold (T<25K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity-temperature (L-T) relation, where the average dust temperature of log[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase.

We present CO(2-1) observations of the submillimetre galaxy ALESS65.1 performed with the Australia Telescope Compact Array at 42.3 GHz. A previous ALMA study of submillimetre galaxies in the Extended Chandra Deep Field South detected [CII] 157.74 micron emission from this galaxy at a redshift of z = 4.44. No CO(2-1) emission was detected but we derive a firm upper limit to the cold gas mass in ALESS65.1 of M_gas < 1.7 x 10^10 M_odot. The estimated gas depletion timescale is <50 Myr, which is similar to other high redshift SMGs, and consistent with z > 4 SMGs being the likely progenitors of massive red-and-dead galaxies at z > 2. The ratio of the [CII], CO and far-infrared luminosities implies a strong far-ultraviolet field of G_0 > 10^3, as seen in Galactic star forming regions or local ULIRGs. The observed L_[CII]/L_FIR = 2.3 x 10^{-3} is high compared to local ULIRGs and, combined with L_[CII]/L_CO > 2700, it is consistent with ALESS65.1 either having an extended (several kpc) [CII] emitting region or lower than solar metallicity.

The high-redshift radio galaxy MRC 1138-262 (`Spiderweb Galaxy'; z = 2.16), is one of the most massive systems in the early Universe and surrounded by a dense `web' of proto-cluster galaxies. Using the Australia Telescope Compact Array, we detected CO(1-0) emission from cold molecular gas -- the raw ingredient for star formation -- across the Spiderweb Galaxy. We infer a molecular gas mass of M(H2) = 6x10^10 M(sun) (for M(H2)/L'(CO)=0.8). While the bulk of the molecular gas coincides with the central radio galaxy, there are indications that a substantial fraction of this gas is associated with satellite galaxies or spread across the inter-galactic medium on scales of tens of kpc. In addition, we tentatively detect CO(1-0) in the star-forming proto-cluster galaxy HAE 229, 250 kpc to the west. Our observations are consistent with the fact that the Spiderweb Galaxy is building up its stellar mass through a massive burst of widespread star formation. At maximum star formation efficiency, the molecular gas will be able to sustain the current star formation rate (SFR ~ 1400 M(sun)/yr, as traced by Seymour et al.) for about 40 Myr. This is similar to the estimated typical lifetime of a major starburst event in infra-red luminous merger systems.

The new evolutionary code Pegase.3 is used for the spectral synthesis of continuous optical-Spitzer-Herschel-submm spectral energy distributions (SEDs) of two z=3.8 radio galaxies 4C 41.17 and TN J2007-1316. Both galaxies are selected from the Herschel Radio Galaxy Evolution Project for their faint AGN contribution and because they show evidence of a large stellar contribution to their bolometric luminosity. Pegase.3 coherently models the reprocessing of the stellar luminosity to dust emission, allowing to build UV to IR-submm libraries by types. These templates allow to predict SEDs at redshifts z in the observer's frame by using z=0 templates corrected for cosmology and evolution. Our best fits are a sum of two distinct evolving stellar populations: i) a very massive starburst observed ~30 Myrs after formation from optical and far-IR peaks ii) a significantly older (~1Gyr old) massive early-type population revealed from mid-IR Spitzer data. The AGN emission fits a roughly thermal model. These results confirm that many of the stellar populations in high-z radio galaxies were already formed by massive (a few thousands of billion solar masses) starbursts at z=4, confirming the previous K-z diagram interpretation. Gas-rich mergers and/or jet-cloud interactions are favored for triggering the intense star formation, possibly in relation with massive super black holes of these radio galaxies. These multiple stellar populations may be a generic feature of the luminous infrared radio galaxy population as a whole and strongly constrain mass accumulation and galaxy formation models.

We present the mid-infrared spectra of seven of the most powerful radio-galaxies known to exist at 1.5 < z < 2.6. The radio emission of these sources is dominated by the AGN with 500 MHz luminosities in the range 10^27.8 - 10^29.1 W/Hz. The AGN signature is clearly evident in the mid-infrared spectra, however, we also detect polycyclic aromatic hydrocarbons emission, indicative of prodigious star formation at a rate of up to ~1000 Msun/yr. Interestingly, we observe no significant correlation between AGN power and star formation in the host galaxy. We also find most of these radio galaxies to have weak 9.7 um silicate absorption features (tau_{9.7 um} < 0.8) which implies that their mid-IR obscuration is predominantly due to the dusty torus that surrounds the central engine, rather than the host galaxy. The tori are likely to have an inhomogeneous distribution with the obscuring structure consisting of individual clouds. We estimate that these radio galaxies have already formed the bulk of their stellar mass and appear to lie at a stage in their evolution where the obscured AGN dominates the energy output of the system but star formation is also prevalent.

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

In the lead-up to the Square Kilometre Array (SKA) project, several next-generation radio telescopes and upgrades are already being built around the world. These include APERTIF (The Netherlands), ASKAP (Australia), eMERLIN (UK), VLA (USA), e-EVN (based in Europe), LOFAR (The Netherlands), Meerkat (South Africa), and the Murchison Widefield Array (MWA). Each of these new instruments has different strengths, and coordination of surveys between them can help maximise the science from each of them. A radio continuum survey is being planned on each of them with the primary science objective of understanding the formation and evolution of galaxies over cosmic time, and the cosmological parameters and large-scale structures which drive it. In pursuit of this objective, the different teams are developing a variety of new techniques, and refining existing ones. Here we describe these projects, their science goals, and the technical challenges which are being addressed to maximise the science return.

We present Keck spectroscopic observations and redshifts for a sample of 767 Herschel-SPIRE selected galaxies (HSGs) at 250, 350, and 500um, taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey (HerMES) peaks at z=0.85, with 731 sources at z<2 and a tail of sources out to z~5. We measure more significant disagreement between photometric and spectroscopic redshifts (<delta_z>/(1+z)>=0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. We estimate that the vast majority (72-83%) of z<2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1mm. We estimate the luminosity function and implied star-formation rate density contribution of HSGs at z<1.6 and find overall agreement with work based on 24um extrapolations of the LIRG, ULIRG and total infrared contributions. This work significantly increased the number of spectroscopically confirmed infrared-luminous galaxies at z>>0 and demonstrates the growing importance of dusty starbursts for galaxy evolution studies and the build-up of stellar mass throughout cosmic time. [abridged]

We present spectroscopic observations for a sample of 36 Herschel-SPIRE 250-500um selected galaxies (HSGs) at 2<z<5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Redshifts are confirmed as part of a large redshift survey of Herschel-SPIRE-selected sources covering ~0.93deg^2 in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts at 24um or 1.4GHz using a cross-identification likelihood matching method. Individual source luminosities range from log(L_IR/Lsun)=12.5-13.6 (corresponding to star formation rates 500-9000Msun/yr, assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well characterized across large areas of sky in contrast to most z>2 HyLIRGs in the literature which are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at (7+-2)x10^(-3)Msun/yr h^3Mpc^(-3) at z~2.5, which is >10% of the estimated total star formation rate density (SFRD) of the Universe from optical surveys. The contribution at z~4 has a lower limit of 3x10^(-3)Msun/yr h^3 Mpc^(-3), ~>20% of the estimated total SFRD. This highlights the importance of extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, even at the earliest epochs.

We present Herschel observations at 70, 160, 250, 350 and 500 micron of the environment of the radio galaxy 4C+41.17 at z = 3.792. About 65% of the extracted sources are securely identified with mid-IR sources observed with the Spitzer Space Telescope at 3.6, 4.5, 5.8, 8 and 24 micron. We derive simple photometric redshifts, also including existing 850 micron and 1200 micron data, using templates of AGN, starburst-dominated systems and evolved stellar populations. We find that most of the Herschel sources are foreground to the radio galaxy and therefore do not belong to a structure associated with 4C+41.17. We do, however, find that the SED of the closest (~ 25" offset) source to the radio galaxy is fully consistent with being at the same redshift as 4C+41.17. We show that finding such a bright source that close to the radio galaxy at the same redshift is a very unlikely event, making the environment of 4C+41.17 a special case. We demonstrate that multi-wavelength data, in particular on the Rayleigh-Jeans side of the spectral energy distribution, allow us to confirm or rule out the presence of protocluster candidates that were previously selected by single wavelength data sets.

Lyman Break Galaxies (LBGs) are widely thought to be prototypical young galaxies in the early universe, particularly representative of those undergoing massive events of star formation. Therefore, LBGs should produce significant amounts of X-ray emission. We aim to trace the X-ray luminosity of Lyman Break Galaxies across cosmic time and from that derive constraints on their star formation history. We utilize the newly released 4 Ms mosaic obtained with the Chandra X-ray Observatory, the deepest X-ray image to date, alongside with the superb spectroscopic data sets available in the CDF-S survey region to construct large but nearly uncontaminated samples of LBGs across a wide range of redshift (0.5 < z < 4.5) which can be used as input samples for stacking experiments. This approach allows us to trace the X-ray emission of Lyman Break Galaxies to even lower, previously unreachable, flux density limits (~10^-18 mW m^-2) and therefore to larger redshifts. We reliably detect soft-band X-ray emission from all our input redshift bins except for the highest redshift (z~4) one. From that we derive rest-frame 2-10 keV luminosities and infer star formation rates and stellar masses. We find that star formation in LBGs peaks at a redshift of z_peak~3.5 and then decreases quickly. We also see a characteristic peak in the specific star formation rate (sSFR=SFR/M_*) at this redshift. Furthermore, we calculate the contribution of LBGs to the total cosmic star formation rate density (SFRD) and find that the contribution of LBGs is negligible. Therefore, we conclude that most of the star formation in the early universe takes place in lower luminosity galaxies as suggested by hierarchical structure formation models.

We examine the relative orientation of radio jets and dusty tori surrounding the AGN in powerful radio galaxies at z > 1. The radio core dominance R = P(20 GHz) /P(500 MHz) serves as an orientation indicator, measuring the ratio between the anisotropic Doppler-beamed core extended core emission and the isotropic lobe emission. Assuming a fixed cylindrical geometry for the hot, dusty torus, we derive its inclination i by fitting optically-thick radiative transfer models to spectral energy distributions obtained with the Spitzer Space Telescope. We find a highly significant anti-correlation (p < 0.0001) between R and i in our sample of 35 type 2 AGN combined with a sample of 18 z = 1 3CR sources containing both type 1 and 2 AGN. This analysis provides observational evidence both for the Unified scheme of AGN and for the common assumption that radio jets are in general perpendicular to the plane of the torus. The use of inclinations derived from mid-infrared photometry breaks several degeneracies which have been problematic in earlier analyses. We illustrate this by deriving the core Lorentz factor Gamma from the R-i anti-correlation, finding Gamma > 1.3.

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

We present a pan-chromatic analysis of an unprecedented sample of 1402 250 micron-selected galaxies at z < 0.5 (mean z = 0.24) from the Herschel-ATLAS survey. We complement our Herschel 100-500 micron data with UV-K-band photometry from the Galaxy And Mass Assembly (GAMA) survey and apply the MAGPHYS energy-balance technique to produce pan-chromatic SEDs for a representative sample of 250 micron selected galaxies spanning the most recent 5 Gyr of cosmic history. We derive estimates of physical parameters, including star formation rates, stellar masses, dust masses and infrared luminosities. The typical H-ATLAS galaxy at z < 0.5 has a far-infrared luminosity in the range 10^10 - 10^12 Lsolar (SFR: 1-50 Msolar/yr) thus is broadly representative of normal star forming galaxies over this redshift range. We show that 250 micron-selected galaxies contain a larger mass of dust at a given infra-red luminosity or star formation rate than previous samples selected at 60 micron from IRAS. We derive typical SEDs for H-ATLAS galaxies, and show that the emergent SED shape is most sensitive to specific star formation rate. The optical-UV SEDs also become more reddened due to dust at higher redshifts. Our template SEDs are significantly cooler than existing infra-red templates. They may therefore be most appropriate for inferring total IR luminosities from moderate redshift submillimetre selected samples and for inclusion in models of the lower redshift submillimetre galaxy populations.

The Australia Telescope Large Area Survey (ATLAS) has surveyed seven square degrees of sky around the Chandra Deep Field South (CDFS) and the European Large Area ISO Survey - South 1 (ELAIS-S1) fields at 1.4 GHz. ATLAS aims to reach a uniform sensitivity of $10 \mu$Jy beam$^{-1}$ rms over the entire region with data release 1 currently reaching $\sim30 \mu$Jy beam$^{-1}$ rms. Here we present 466 new spectroscopic redshifts for radio sources in ATLAS as part of our optical follow-up program. Of the 466 radio sources with new spectroscopic redshifts, 142 have star-forming optical spectra, 282 show evidence for AGN in their optical spectra, 10 have stellar spectra and 32 have spectra revealing redshifts, but with insufficient features to classify. We compare our spectroscopic classifications with two mid-infrared diagnostics and find them to be in broad agreement. We also construct the radio luminosity function for star-forming galaxies to z $= 0.5$ and for AGN to z $= 0.8$. The radio luminosity function for star-forming galaxies appears to be in good agreement with previous studies. The radio luminosity function for AGN appears higher than previous studies of the local AGN radio luminosity function. We explore the possibility of evolution, cosmic variance and classification techniques affecting the AGN radio luminosity function. ATLAS is a pathfinder for the forthcoming EMU survey and the data presented in this paper will be used to guide EMU's survey design and early science papers.

We present a detailed study of the infrared spectral energy distribution of the high-redshift radio galaxy MRC 1138-26 at z = 2.156, also known as the Spiderweb Galaxy. By combining photometry from Spitzer, Herschel and LABOCA we fit the rest-frame 5-300 um emission using a two component, starburst and active galactic nucleus (AGN), model. The total infrared (8 - 1000 um) luminosity of this galaxy is (1.97+/-0.28)x10^13 Lsun with (1.17+/-0.27) and (0.79+/-0.09)x10^13 Lsun due to the AGN and starburst components respectively. The high derived AGN accretion rate of \sim20% Eddington, and the measured star formation rate (SFR) of 1390pm150 Msun/yr, suggest that this massive system is in a special phase of rapid central black hole and host galaxy growth, likely caused by a gas rich merger in a dense environment. The accretion rate is sufficient to power both the jets and the previously observed large outflow. The high SFR and strong outflow suggest this galaxy could potentially exhaust its fuel for stellar growth in a few tens of Myr, although the likely merger of the radio galaxy with nearby satellites suggest bursts of star formation may recur again on time scales of several hundreds of Myr. The age of the radio lobes implies the jet started after the current burst of star formation, and therefore we are possibly witnessing the transition from a merger-induced starburst phase to a radio-loud AGN phase. We also note tentative evidence for [CII]158um emission. This paper marks the first results from the Herschel Galaxy Evolution Project (Project HeRGE), a systematic study of the evolutionary state of 71 high redshift, 1 < z < 5.2, radio galaxies.

We report new, sensitive observations of two z ~ 3-3.5 FIR-luminous radio galaxies, 6C1909+72 and B3J2330+3927, in 12CO J=1-0 with the Karl Jansky VLA and at 100-500um using Herschel, alongside new and archival 12CO J=4-3 observations from IRAM PdBI. We introduce a new colour-colour diagnostic plot to constrain the redshifts of several distant, dusty galaxies in our target fields. A bright SMG near 6C1909+72 likely shares the same node or filament as the signpost AGN, but it is not detected in CO despite ~20,000 km/s of velocity coverage. Also in the 6C1909+72 field, a large, red dust feature spanning ~500 kpc is aligned with the radio jet. We suggest several processes by which metal-rich material may have been transported, favouring a collimated outflow reminiscent of the jet-oriented metal enrichment seen in local cluster environments. Our interferometric imaging reveals a gas-rich companion to B3J2330+3927; indeed, all bar one of the eight z >~ 2 radio galaxies (or companions) detected in CO provide some evidence that starburst activity in radio-loud AGN at high redshift is driven by the interaction of two or more gas-rich systems in which a significant mass of stars has already formed, rather than via steady accretion of cold gas from the cosmic web. We find that the CO Tb ratios in radio-loud AGN host galaxies are significantly higher than those seen in similarly intense starbursts where AGN activity is less pronounced. Our most extreme example, where L'(CO4-3)/L'(CO1-0) > 2.7, provides evidence that significant energy is being deposited rapidly into the molecular gas via X-rays and/or mechanical (`quasar-mode') feedback from the AGN, leading to a high degree of turbulence globally and a low optical depth in 12CO - feedback that may lead to the cessation of star formation on a timescale commensurate with that of the jet activity, <~10 Myr.

A deep Spitzer Infrared Spectrograph map of the PKS 1138-26 galaxy protocluster reveals ultraluminous PAH emission from obscured star formation in three protocluster galaxies, including Halpha-emitter (HAE) 229, HAE 131, and the central Spiderweb Galaxy. Star formation rates of 500-1100 Msun/yr are estimated from the 7.7 micron PAH feature. At such prodigious formation rates, the galaxy stellar masses will double in 0.6-1.1 Gyr. We are viewing the peak epoch of star formation for these protocluster galaxies. However, it appears that extinction of Halpha is much greater (factor of 40) in the two ULIRG HAEs compared to the Spiderweb. This may be attributed to different spatial distributions of star formation--nuclear star formation in the HAEs versus extended star formation in accreting satellite galaxies in the Spiderweb. We find extremely luminous mid-IR rotational line emission from warm molecular hydrogen in the Spiderweb Galaxy, with L(H2 0-0 S(3))= 1.4E44 erg/s (3.7E10 Lsun), 20 times more luminous than any previously known H2 emission galaxy (MOHEG). Depending on temperature, this corresponds to a very large mass of >9E6-2E9 Msun of T>300 K molecular gas, heated by the PKS 1138-26 radio jet, acting to quench nuclear star formation. There is >8 times more warm H2 at these temperatures in the Spiderweb than what has been seen in low-redshift (z<0.2) radio galaxies, indicating that the Spiderweb may have a larger reservoir of molecular gas than more evolved radio galaxies. This is the highest redshift galaxy yet in which warm molecular hydrogen has been directly detected.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

Taking advantage of the impressive sensitivity of Spitzer to detect massive galaxies at high redshift, we study the mid-infrared environments of powerful, high-redshift radio galaxies at 1.2<z<3. Galaxy cluster member candidates were isolated using a single Spitzer/IRAC mid-infrared color criterion, [3.6]-[4.5]>-0.1 (AB), in the fields of 48 radio galaxies at 1.2<z<3. This simple IRAC color selection is effective at identifying galaxies at z>1.2. Using a counts-in-cell analysis, we identify a field as overdense when 15 or more red IRAC sources are found within 1arcmin (i.e.,~0.5Mpc at 1.2<z<3) of the radio galaxy to the 5sigma flux density limits of our IRAC data (f3.6=11.0uJy, f4.5=13.4uJy). We find that radio galaxies lie preferentially in medium to dense regions, with 73% of the targeted fields denser than average. Our (shallow) 120s data permit the rediscovery of previously known clusters and protoclusters associated with radio galaxies as well as the discovery of new promising galaxy cluster candidates at z>1.2.

We present a statistical study of the environments of 63 high-redshift radio galaxies (HzRGs) between redshifts 1<z<5.2, using the 24um, waveband of the MIPS instrument aboard the Spitzer Space Telescope. Using a counts-in-cell analysis, a statistically significant source overdensity is found in 1.75arcmin radius circular cells centred on the HzRGs when compared to reference fields. We report an average overdensity of delta (= {N}_{targets} / {N}_{reference}) = 2.2 +/- 1.2 at a flux density cut of f24um=0.3mJy. This result implies that HzRGs are likely to lie in protoclusters of active and star-forming galaxies at high redshift. Over 95% of our targeted HzRGs lie in higher than average density fields. Further, 20 (32%) of our selected fields are found to be overdense to at least a 3sigma significance, of which 9 are newly identified protocluster candidates. We observe a weak correlation between redshift and 24um, source density, and discuss the populations being probed at different redshifts. In our uniformly selected sample, which was designed to cover two orders of magnitude in radio luminosity throughout z=1-4, we find that the 24um, source density does not depend on radio luminosity. We also compare this result with recent work describing IRAC source overdensities around the same HzRGs and find correlations between the results.

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 describe the generation of single-band point source catalogues from submillimetre Herschel-SPIRE observations taken as part of the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found by means of peak-finding and the fitting of a Gaussian point-response function. With highly-confused images, careful checks must be made on the completeness and flux density accuracy of the detected sources. This is done by injecting artificial sources into the images and analysing the resulting catalogues. Measured flux densities at which 50 per cent of injected sources result in good detections at (250, 350, 500) {\mu}m range from (11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depending on the depth of the observation (where a `good' detection is taken to be one with positional offset less than one full-width half-maximum of the point-response function, and with the measured flux density within a factor of 2 of the flux density of the injected source). This paper acts as a reference for the 2010 July HerMES public data release.

Ultra Steep Spectrum (USS) radio sources have been successfully used to select powerful radio sources at high redshifts (z>~2). Typically restricted to large-sky surveys and relatively bright radio flux densities, it has gradually become possible to extend the USS search to sub-mJy levels, thanks to the recent appearance of sensitive low-frequency radio facilities. Here a first detailed analysis of the nature of the faintest USS sources is presented. By using Giant Metrewave Radio Telescope and Very Large Array radio observations of the Lockman Hole at 610 MHz and 1.4 GHz, a sample of 58 USS sources, with 610 MHz integrated fluxes above 100 microJy, is assembled. Deep infrared data at 3.6 and 4.5 micron from the Spitzer Extragalactic Representative Volume Survey (SERVS) is used to reliably identify counterparts for 48 (83%) of these sources, showing an average total magnitude of [3.6](AB)=19.8 mag. Spectroscopic redshifts for 14 USS sources, together with photometric redshift estimates, improved by the use of the deep SERVS data, for a further 19 objects, show redshifts ranging from z=0.1 to z=2.8, peaking at z~0.6 and tailing off at high redshifts. The remaining 25 USS sources, with no redshift estimate, include the faintest [3.6] magnitudes, with 10 sources undetected at 3.6 and 4.5 micron (typically [3.6]>22-23 mag, from local measurements), which suggests the likely existence of higher redshifts among the sub-mJy USS population. The comparison with the Square Kilometre Array Design Studies Simulated Skies models indicate that Fanaroff-Riley type I radio sources and radio-quiet Active Galactic Nuclei may constitute the bulk of the faintest USS population, and raises the possibility that the high efficiency of the USS technique for the selection of high redshift sources remains even at the sub-mJy level.

For the first time, we investigate the X-ray/infrared (IR) correlation for star-forming galaxies at z~1, using SPIRE submm data from the recently-launched Herschel Space Observatory and deep X-ray data from the 2Ms Chandra deep field north (CDFN) survey. We examine the X-ray/IR correlation in the soft X-ray (SX, 0.5-2 keV) and hard X-ray (HX, 2-10 keV) bands by comparing our z~1 SPIRE-detected star-forming galaxies (SFGs) to equivalently IR-luminous (L_IR >10^10 L_sun) samples in the local/low redshift Universe. Our results suggest that the X-ray/IR properties of the SPIRE SFGs are on average similar to those of their local counterparts, as we find no evidence for evolution in the L_SX/L_IR and L_HX/L_IR ratios with redshift. We note however, that at all redshifts, both L_SX/L_IR and L_HX/L_IR are strongly dependent on IR luminosity, with luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs,L_IR >10^11 L_sun) having up to an order of magnitude lower values than normal infrared galaxies (L_IR <10^11 L_sun). We derive a L_SX-L_IR relation and confirm the applicability of an existing L_HX-L_IR relation for both local and distant LIRGs and ULIRGs, consistent with a scenario where X-ray luminosity is correlated with the star-formation rate (SFR).

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

Infrared-faint radio sources (IFRSs) are a rare class of object which are relatively bright at radio wavelengths but very faint at infrared and optical wavelengths. Here we present sensitive near-infrared observations of a sample of these sources taken as part of the Spitzer Extragalactic Representative Volume Survey (SERVS). Nearly all the IFRSs are undetected at a level of ~ 1 \mu$Jy in these new deep observations, and even the detections are consistent with confusion with unrelated galaxies. A stacked image implies that the median flux density is $S_{3.6\mu m} ~ 0.2$ \mu$Jy or less, giving extreme values of the radio-infrared flux density ratio. Comparison of these objects with known classes of object suggests that the majority are probably high-redshift radio-loud galaxies, possibly suffering from significant dust extinction.

As part of the Herschel Multi-tiered Extragalactic Survey we have investigated the rest-frame far-infrared (FIR) properties of a sample of more than 4800 Lyman Break Galaxies (LBGs) in the Great Observatories Origins Deep Survey North field. Most LBGs are not detected individually, but we do detect a sub-sample of 12 objects at 0.7 < z < 1.6 and one object at z ~ 2.0. The ones detected by Herschel SPIRE have redder observed NUV-U and U-R colors than the others, while the undetected ones have colors consistent with average LBGs at z > 2.5. The UV-to-FIR spectral energy distributions of the objects detected in the rest-frame FIR are investigated using the code CIGALE to estimate physical parameters. We find that LBGs detected by SPIRE are high mass, luminous infrared galaxies. It appears that LBGs are located in a triangle-shaped region in the A_FUV vs. Log L_FUV diagram limited by A_FUV=0 at the bottom and by a diagonal following the temporal evolution of the most massive galaxies from the bottom-right to the top-left of the diagram. This upper envelop can be used as upper limits for the UV dust attenuation as a function of L_FUV}. The limits of this region are well explained using a closed-box model, where the chemical evolution of galaxies produces metals, which in turn lead to higher dust attenuation when the galaxies age.

This paper presents a study of the environments of SDSS Quasi-Stellar-Objects (QSOs) in the Spitzer Extragalactic Representative Volume Survey (SERVS). We concentrate on the high-redshift QSOs as these have not been studied in large numbers with data of this depth before. We use the IRAC 3.6-4.5{\mu}m colour of objects and ancillary r-band data to filter out as much foreground contamination as possible. This technique allows us to find a significant (> 4-{\sigma}) over-density of galaxies around QSOs in a redshift bin centred on z ~ 2.0 and a (> 2-{\sigma}) over-density of galaxies around QSOs in a redshift bin centred on z ~ 3.3. We compare our findings to the predictions of a semi-analytic galaxy formation model, based on the {\Lambda}CDM millennium simulation, and find for both redshift bins that the model predictions match well the source-density we have measured from the SERVS data.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

We report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

We examine the rest-frame far-infrared emission from powerful radio sources with 1.4GHz luminosity densities of 25<=log(L_1.4/WHz^-1)<=26.5 in the extragalactic Spitzer First Look Survey field. We combine Herschel/SPIRE flux densities with Spitzer/IRAC and MIPS infrared data to obtain total (8-1000um) infrared luminosities for these radio sources. We separate our sources into a moderate, 0.4<z<0.9, and a high, 1.2<z<3.0, redshift sub-sample and we use Spitzer observations of a z<0.1 3CRR sample as a local comparison. By comparison to numbers from the SKA Simulated Skies we find that our moderate redshift sample is complete and our high redshift sample is 14per cent complete. We constrain the ranges of mean star formation rates (SFRs) to be 3.4-4.2, 18-41 and 80-581Msun/yr for the local, moderate and high redshift samples respectively. Hence, we observe an increase in the mean SFR with increasing redshift which we can parameterise as ~(1+z)^Q, where Q=4.2+/-0.8. However we observe no trends of mean SFR with radio luminosity within the moderate or high redshift bins. We estimate that radio-loud AGN in the high redshift sample contribute 0.1-0.5per cent to the total SFR density at that epoch. Hence, if all luminous starbursts host radio-loud AGN we infer a radio-loud phase duty cycle of 0.001-0.005.

We report the discovery of a significant excess of candidate Halpha emitters (HAEs) in the field of the radio galaxy 4C 23.56 at z=2.483. Using the MOIRCS near-infrared imager on the Subaru Telescope we found 11 candidate emission-line galaxies to a flux limit of ~7.5 10^-17 erg s-1 cm-2, which is about 5 times excess from the expected field counts with ~3-sigma significance. Three of these are spectroscopically confirmed as redshifted Halpha at z=2.49. The distribution of candidate emitters on the sky is tightly confined to a 1.2-Mpc-radius area at z=2.49, locating 4C 23.56 at the western edge of the distribution. Analysis of the deep Spitzer MIPS 24 mu m imaging shows that there is also an excess of faint MIPS sources. All but two of the 11 HAEs are also found in the MIPS data. The inferred star-formation rate (SFR) of the HAEs based on the extinction-corrected Halpha luminosity (median SFR >~100 M_solar yr-1) is similar to those of HAEs in random fields at z~2. On the other hand, the MIPS-based SFR for the HAEs is on average 3.6 times larger, suggesting the existence of the star-formation significanly obscured by dust. The comparison of the Halpha-based star-formation activities of the HAEs in the 4C 23.56 field to those in another proto-cluster around PKS 1138-262 at z=2.16 reveals that the latter tend to have fainter Halpha emission despite similar K-band magnitudes. This suggests that star-formation may be suppressed in the PKS 1138-262 protocluster relative to the 4C 23.56 protocluster. This difference among the HAEs in the two proto-clusters at z > 2 may imply that some massive cluster galaxies are just forming at these epochs with some variation among clusters.

We present intensive quasi-simultaneous X-ray and radio monitoring of the narrow line Seyfert 1 galaxy NGC 4051, over a 16 month period in 2000-2001. Observations were made with the Rossi Timing X-ray Explorer (RXTE) and the Very Large Array (VLA) at 8.4 and 4.8 GHz. In the X-ray band NGC 4051 behaves much like a Galactic black hole binary (GBH) system in a `soft-state'. In such systems, there has so far been no firm evidence for an active, radio-emitting jet like those found in `hard state' GBHs. VLBI observations of NGC 4051 show three co-linear compact components. This structure resembles the core and outer hot spots seen in powerful, jet-dominated, extragalactic radio sources and suggests the existence of a weak jet. Radio monitoring of the core of NGC 4051 is complicated by the presence of surrounding extended emission and by the changing array configurations of the VLA. Only in the A configuration is the core reasonably resolved. We have carefully removed the contaminations of the core by extended emission in the various arrays. The resulting lightcurve shows no sign of large amplitude variability (i.e. factor 50 %) over the 16 month period. Within the most sensitive configuration (A array) we see marginal evidence for radio core variability of ~25% (~0.12 mJy at 8.4GHz) on a 2-week timescale, correlated with X-ray variations. Even if the radio variations in NGC 4051 are real, the percentage variability is much less than in the X-ray band. Within the B configuration observations, where sensitivity is reduced, there is no sign of correlated X-ray/radio variability. The lack of radio variability in NGC 4051, which we commonly see in `hard state' GBHs, may be explained by orientation effects. Another possibility is that the radio emission arises from the X-ray corona, although the linear structure of the compact radio components here is hard to explain.

Infrared-faint radio sources (IFRS) are objects that have flux densities of several mJy at 1.4GHz, but that are invisible at 3.6um when using sensitive Spitzer observations with uJy sensitivities. Their nature is unclear and difficult to investigate since they are only visible in the radio. High-resolution radio images and comprehensive spectral coverage can yield constraints on the emission mechanisms of IFRS and can give hints to similarities with known objects. We imaged a sample of 17 IFRS at 4.8GHz and 8.6GHz with the Australia Telescope Compact Array to determine the structures on arcsecond scales. We added radio data from other observing projects and from the literature to obtain broad-band radio spectra. We find that the sources in our sample are either resolved out at the higher frequencies or are compact at resolutions of a few arcsec, which implies that they are smaller than a typical galaxy. The spectra of IFRS are remarkably steep, with a median spectral index of -1.4 and a prominent lack of spectral indices larger than -0.7. We also find that, given the IR non-detections, the ratio of 1.4GHz flux density to 3.6um flux density is very high, and this puts them into the same regime as high-redshift radio galaxies. The evidence that IFRS are predominantly high-redshift sources driven by active galactic nuclei (AGN) is strong, even though not all IFRS may be caused by the same phenomenon. Compared to the rare and painstakingly collected high-redshift radio galaxies, IFRS appear to be much more abundant, but less luminous, AGN-driven galaxies at similar cosmological distances.

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

We present results from a comprehensive imaging survey of 70 radio galaxies at redshifts 1<z<5.2 using all three cameras onboard the Spitzer Space Telescope. The resulting spectral energy distributions unambiguously show a stellar population in 46 sources and hot dust emission associated with the active nucleus in 59. Using a new restframe S_3um/S_1.6um versus S_um/S_3um criterion, we identify 42 sources where the restframe 1.6um emission from the stellar population can be measured. For these radio galaxies, the median stellar mass is high, 2x10^11 M_sun, and remarkably constant within the range 1<z<3. At z>3, there is tentative evidence for a factor of two decrease in stellar mass. This suggests that radio galaxies have assembled the bulk of their stellar mass by z~3, but confirmation by more detailed decomposition of stellar and AGN emission is needed. The restframe 500 MHz radio luminosities are only marginally correlated with stellar mass but are strongly correlated with the restframe 5um hot dust luminosity. This suggests that the radio galaxies have a large range of Eddington ratios. We also present new Very Large Array 4.86 and 8.46 GHz imaging of 14 radio galaxies and find that radio core dominance --- an indicator of jet orientation --- is strongly correlated with hot dust luminosity. While all of our targets were selected as narrow-lined, type 2 AGNs, this result can be understood in the context of orientation-dependent models if there is a continuous distribution of orientations from obscured type 2 to unobscured type 1 AGNs rather than a clear dichotomy. Finally, four radio galaxies have nearby (<6") companions whose mid-IR colors are suggestive of their being AGNs. This may indicate an association between radio galaxy activity and major mergers.

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

We examine the selection characteristics of infrared and sub-mm surveys with IRAS, Spitzer, BLAST, Herschel and SCUBA and identify the range of dust temperatures these surveys are sensitive to, for galaxies in the ULIRG luminosity range (12<log(LIR)<13), between z=0 and z=4. We find that the extent of the redshift range over which surveys are unbiased is a function of the wavelength of selection, flux density limit and ULIRG luminosity. Short wavelength (<200{\mu}m) surveys with IRAS, Spitzer/MIPS and Herschel/PACS are sensitive to all SED types in a large temperature interval (17-87K), over a substantial fraction of their accessible redshift range. On the other hand, long wavelength (>200{\mu}m) surveys with BLAST, Herschel/ SPIRE and SCUBA are significantly more sensitive to cold ULIRGs, disfavouring warmer SEDs even at low redshifts. We evaluate observations in the context of survey selection effects, finding that the lack of cold ULIRGs in the local (z<0.1) Universe is not a consequence of selection and that the range of ULIRG temperatures seen locally is only a subset of a much larger range which exists at high redshift. We demonstrate that the local luminosity-temperature (L-T) relation, which indicates that more luminous sources are also hotter, is not applicable in the distant Universe when extrapolated to the ULIRG regime, because the scatter in observed temperatures is too large. Finally, we show that the difference between the ULIRG temperature distributions locally and at high redshift is not the result of galaxies becoming colder due to an L-T relation which evolves as a function of redshift. Instead, they are consistent with a picture where the evolution of the infrared luminosity function is temperature dependent, i.e. cold galaxies evolve at a faster rate than their warm counterparts.

(abridged) We present the first study of the farIR properties of high redshift, radio-selected ULIRGs using deep observations obtained with SPIRE from the Herschel Multi-tiered Extragalactic Survey (HerMES). These galaxies span a large range of 850um fluxes from submillimetre-luminous ~10mJy "SCUBA galaxies -- SMGs" to ~1.5mJy from stacked SCUBA non-detections, thus likely representing a complete distribution of ULIRG spectral energy distributions. From Keck spectroscopic surveys in the Lockman-North field we identified a sample of 31 SMGs and 37 submillimetre-faint, optically-faint radio galaxies (OFRGs), all with radio-inferred IR luminosities >10^12 Lsun. These galaxies were cross-identified with SPIRE 250, 350 and 500um catalogs based on fluxes extracted at 24um positions in the SWIRE survey, yielding a sample of more than half of the galaxies well detected in at least two of the SPIRE bandpasses. By fitting greybody dust models to the SPIRE photometry together with SCUBA 850um measurements, we infer dust temperatures and far-infrared luminosities. The OFRGs detected by SPIRE have median <T_d>= 41+-5 K and the SMGs have <T_d>= 34+-5 K, both in reasonable agreement with previous (pre-Herschel) estimates, reaffirming that the local FIR/radio correlation holds (at least for this subset of high-z ULIRGs) at high redshift. Our observations firstly confirm that a substantial fraction of OFRGs exhibit large infrared luminosities corresponding to SFRs of ~400 Msun/yr. The SPIRE observations secondly confirm the higher dust temperatures for these OFRGs than similarly selected SMGs, consistent with early predictions of the submm-faint radio populations. Our observations also clearly confirm the large infrared luminosities of most SMGs selected with S850um>5 mJy and radio and strong 24um detections, corresponding to SFRs of ~700 Msun/yr.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

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.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

We study the environments of 6 radio galaxies at 2.2 < z < 2.6 using wide-field near-infrared images. We use colour cuts to identify galaxies in this redshift range, and find that three of the radio galaxies are surrounded by significant surface overdensities of such galaxies. The excess galaxies that comprise these overdensities are strongly clustered, suggesting they are physically associated. The colour distribution of the galaxies responsible for the overdensity are consistent with those of galaxies that lie within a narrow redshift range at z ~ 2.4. Thus the excess galaxies are consistent with being companions of the radio galaxies. The overdensities have estimated masses in excess of 10^14 solar masses, and are dense enough to collapse into virizalised structures by the present day: these structures may evolve into groups or clusters of galaxies. A flux-limited sample of protocluster galaxies with K < 20.6 mag is derived by statistically subtracting the fore- and background galaxies. The colour distribution of the protocluster galaxies is bimodal, consisting of a dominant blue sequence, comprising 77 +/- 10% of the galaxies, and a poorly populated red sequence. The blue protocluster galaxies have similar colours to local star-forming irregular galaxies (U -V ~ 0.6), suggesting most protocluster galaxies are still forming stars at the observed epoch. The blue colours and lack of a dominant protocluster red sequence implies that these cluster galaxies form the bulk of their stars at z < 3.

We present a study of protoclusters associated with high redshift radio galaxies. We imaged MRC1017-220 (z=1.77) and MRC0156-252 (z=2.02) using the near-infrared wide-field (7.5'x7.5') imager VLT/HAWK-I in the Y, H and Ks bands. We present the first deep Y-band galaxy number counts within a large area (200 arcmin2). We then develop a purely near-infrared colour selection technique to isolate galaxies at 1.6<z<3 that may be associated with the two targets, dividing them into (i) red passively evolving or dusty star-forming galaxies or (ii) blue/star-formation dominated galaxies with little or no dust. Both targeted fields show an excess of star-forming galaxies with respect to control fields. No clear overdensity of red galaxies is detected in the surroundings of MRC1017-220 although the spatial distribution of the red galaxies resembles a filament-like structure within which the radio galaxy is embedded. In contrast, a significant overdensity of red galaxies is detected in the field of MRC0156-252, ranging from a factor of 2-3 times the field density at large scales (2.5Mpc, angular distance) up to a factor of 3-4 times the field density within a 1Mpc radius of the radio galaxy. Half of these red galaxies have colours consistent with red sequence models at z~2, with a large fraction being bright (Ks<21.5, i.e. massive). In addition, we also find a small group of galaxies within 5" of MRC0156-252 suggesting that the radio galaxy has multiple companions within ~50 kpc. We conclude that the field of MRC0156-252 shows many remarkable similarities with the well-studied protocluster surrounding PKS1138-262 (z=2.16) suggesting that MRC0156-252 is associated with a galaxy protocluster at z~2.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

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 use the multi-epoch, mid-infrared Spitzer Deep, Wide-Field Survey to investigate the variability of 474,179 objects in 8.1 deg^2 of the NDWFS Bootes field. We perform a Difference Image Analysis of the four available epochs between 2004 and 2008, focusing on the deeper 3.6 and 4.5 micron bands. We find that 1.1% of the studied sample meet our standard selection criteria for being classed as a variable source. We require that the 3.6 and 4.5 micron light-curves are strongly correlated (r>0.8) and that their joint variance exceeds that for all sources with the same magnitude by 2 sigma. We then examine the mid-IR colors of the variable sources and match them with X-ray sources from the XBootes survey, radio catalogs, 24 micron-selected AGN candidates, and spectroscopically identified AGNs from the AGN and Galaxy Evolution Survey (AGES). Based on their mid-IR colors, most of the variable sources are AGNs (76%), with smaller contributions from stars (11%), galaxies (6%), and unclassified objects. Most of the stellar, galaxy and unclassified sources are false positives. For our standard selection criteria, 11-12% of the mid-IR counterparts to X-ray sources, 24 micron-selected AGN candidates and spectroscopically identified AGNs show variability. Mid-IR AGN variability can be well described by a single power-law structure function with a power-law index of 0.5 at both 3.6 and 4.5 microns, and an amplitude of 0.1 mag on rest-frame time scales of 2 years. The variability amplitude is higher for shorter rest-frame wavelengths and lower luminosities. (Abridged)

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

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 report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

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.

The current deepest radio surveys detect hundreds of sources per square degree below 0.1mJy. There is a growing consensus that a large fraction of these sources are dominated by star formation although the exact proportion has been debated in the literature. However, the low luminosity of these galaxies at most other wavelengths makes determining the nature of individual sources difficult. If future, deeper surveys performed with the next generation of radio instrumentation are to reap high scientific reward we need to develop reliable methods of distinguishing between radio emission powered by active galactic nuclei (AGN) and that powered by star formation. In particular, we believe that such discriminations should be based on purely radio, or relative to radio, diagnostics. These diagnostics include radio morphology, radio spectral index, polarisation, variability, radio luminosity and flux density ratios with non-radio wavelengths e.g. with different parts of the infrared (IR) regime. We discuss the advantages and limitations of these various diagnostics methods with current and future surveys. However, weeding AGN out of deep radio surveys can already provide several insights into the star formation at high redshift. As well as reproducing the well known rise with redshift in the comoving star formation rate density, we also see evidence for the continued dominance of LIRGs and ULIRGs to the total star forming budget across redshifts 1-3. Additionally, while we see that the IR-radio relation for star forming galaxies does hold to high redshifts (z>1) there is a mild deviation depending on the IR waveband used and the range of IR SEDs found. We will discuss the possible reasons behind this change in properties.

In this paper we examine the contribution of galaxies with different infrared (IR) spectral energy distributions (SEDs) to the comoving infrared luminosity density, a proxy for the comoving star formation rate (SFR) density. We characterise galaxies as having either a cold or hot IR SED depending upon whether the rest-frame wavelength of their peak IR energy output is above or below 90um. Our work is based on a far-IR selected sample both in the local Universe and at high redshift, the former consisting of IRAS 60um-selected galaxies at z<0.07 and the latter of Spitzer 70um selected galaxies across 0.1<z<1. We find that the total IR luminosity densities for each redshift/luminosity bin agree well with results derived from other deep mid/far-IR surveys. At z<0.07 we observe the previously known results: that moderate luminosity galaxies (L_IR<10^11 Lsun) dominate the total luminosity density and that the fraction of cold galaxies decreases with increasing luminosity, becoming negligible at the highest luminosities. Conversely, above z=0.1 we find that luminous IR galaxies (L_IR>10^11 Lsun), the majority of which are cold, dominate the IR luminosity density. We therefore infer that cold galaxies dominate the IR luminosity density across the whole 0<z<1 range, hence appear to be the main driver behind the increase in SFR density up to z~1 whereas local luminous galaxies are not, on the whole, representative of the high redshift population.

We present a spectroscopic campaign to follow-up red colour-selected candidate massive galaxies in two high redshift proto-clusters surrounding radio galaxies. We observed a total of 57 galaxies in the field of MRC0943-242 (z=2.93) and 33 in the field of PKS1138-262 (z=2.16) with a mix of optical and near-infrared multi-object spectroscopy. We confirm two red galaxies in the field of PKS1138-262 at the redshift of the radio galaxy. Based on an analysis of their spectral energy distributions, and their derived star formation rates from the H-alpha and 24um flux, one object belongs to the class of dust-obscured star-forming red galaxies, while the other is evolved with little ongoing star formation. This result represents the first red and mainly passively evolving galaxy to be confirmed as companion galaxies in a z>2 proto-cluster. Both red galaxies in PKS1138-262 are massive, of the order of 4-6x10^11 M_Sol. They lie along a Colour-Magnitude relation which implies that they formed the bulk of their stellar population around z=4. In the MRC0943-242 field we find no red galaxies at the redshift of the radio galaxy but we do confirm the effectiveness of our JHK_s selection of galaxies at 2.3<z<3.1, finding that 10 out of 18 (56%) of JHK_s-selected galaxies whose redshifts could be measured fall within this redshift range. We also serendipitously identify an interesting foreground structure of 6 galaxies at z=2.6 in the field of MRC0943-242. This may be a proto-cluster itself, but complicates any interpretation of the red sequence build-up in MRC0943-242 until more redshifts can be measured.

We explore the environments of two radio galaxies at z~1.5, 7C1751+6809 and 7C1756+6520, using deep optical and near-infrared imaging. Our data cover 15'x15' fields around the radio galaxies. We develop and apply BzK color criteria to select cluster member candidates around the radio galaxies and find no evidence of an overdensity of red galaxies within 2Mpc of 7C1751+6809. In contrast, 7C1756+6520 shows a significant overdensity of red galaxies within 2Mpc of the radio galaxy, by a factor of 3.1 +/- 0.8 relative to the four MUSYC fields. At small separation (r<6''), this radio galaxy also has one z~1.4 evolved galaxy candidate, one z~1.4 star-forming galaxy candidate, and an AGN candidate (at indeterminate redshift). This is suggestive of several close-by companions. Several concentrations of red galaxies are also noticed in the full 7C1756+6520 field, forming a possible large-scale structure of evolved galaxies with a NW-SE orientation. We construct the color-magnitude diagram of red galaxies found near (r<2Mpc) 7C1756+6520, and find a clear red sequence that is truncated at Ks~21.5 (AB). We also find an overdensity of mid-IR selected AGN in the surroundings of 7C1756+6520. These results are suggestive of a proto-cluster at high redshift.

In order to examine the far-IR/radio correlation at high redshift we have studied the Spitzer 70um/far-infrared (far-IR) properties of sub-mJy radio sources from the 13^H XMM-Newton/Chandra Deep Field by redshift and galaxy type: active galactic nucleus (AGN) or star forming galaxy (SFG). We directly detect 70um counterparts (at >3sigma significance) for 22.5% (92/408) of the radio sources, while for the rest we perform stacking analysis by redshift and galaxy type. For the sources detected at 70um we find that the median and scatter of the observed flux density ratio, q70, are similar to previous results in the literature, but with a slight decrease in q70 towards higher redshifts. Of the radio sources detected at 70um 8/92 were already classified as AGN, but two of which maybe SFGs. For the stacked sources we obtain a detection for the SFGs at every redshift bin which implies they have mean flux densities a factor ~5 below the original 70um detection limit. For the stacked AGN we obtain a detection only in our highest redshift bin (1<~z<~5) where we may be sampling hot dust associated with the AGN at rest-frame 12-35um. The combined observed mean value of q70 for the SFGs (detected and non-detected at 70um) decreases gradually with redshift, consistent with tracks derived from empirical spectral energy distributions (SEDs) of local SFGs. Upon closer inspection and when comparing with tracks of appropriate luminosity, the values of q70 broadly agree at low redshift. However, at z~1, the observed q70 (for ULIRGs) is 2sigma below the value seen for local ULIRGs tracks, implying a difference in the SED between local and z~1 ULIRGs. At higher redshifts, the convergence of the tracks and the higher uncertainties in q70 do not allow us to determine if such a difference persists.

The Spitzer Deep, Wide-Field Survey (SDWFS) is a four-epoch infrared survey of ten square degrees in the Bootes field of the NOAO Deep Wide-Field Survey using the IRAC instrument on the Spitzer Space Telescope. SDWFS, a Cycle four Spitzer Legacy project, occupies a unique position in the area-depth survey space defined by other Spitzer surveys. The four epochs that make up SDWFS permit -- for the first time -- the selection of infrared-variable and high proper motion objects over a wide field on timescales of years. Because of its large survey volume, SDWFS is sensitive to galaxies out to z~3 with relatively little impact from cosmic variance for all but the richest systems. The SDWFS datasets will thus be especially useful for characterizing galaxy evolution beyond z~1.5. This paper explains the SDWFS observing strategy and data processing, presents the SDWFS mosaics and source catalogs, and discusses some early scientific findings. The publicly-released, full-depth catalogs contain 6.78, 5.23, 1.20, and 0.96 x 10e5 distinct sources detected to the average 5-sigma, 4" diameter, aperture-corrected limits of 19.77, 18.83, 16.50, and 15.82 Vega mag at 3.6, 4.5, 5.8, and 8.0 micron, respectively. The SDWFS number counts and color-color distribution are consistent with other, earlier Spitzer surveys. At the 6 min integration time of the SDWFS IRAC imaging, more than 50% of isolated FIRST radio sources and more than 80% of on-axis XBootes sources are detected out to 8.0 micron. Finally, we present the four highest proper motion IRAC-selected sources identified from the multi-epoch imaging, two of which are likely field brown dwarfs of mid-T spectral class.

We show that the far-IR properties of distant Luminous and Ultraluminous InfraRed Galaxies (LIRGs and ULIRGs) are on average divergent from analogous sources in the local Universe. Our analysis is based on Spitzer MIPS and IRAC data of L_IR>10^10 L_solar, 70um-selected objects in the 0.1<z<2 redshift range and supported by a comparison with the IRAS Bright Galaxy Sample. The majority of the objects in our sample are described by Spectral Energy Distributions (SEDs) which peak at longer wavelengths than local sources of equivalent total infrared luminosity. This shift in SED peak wavelength implies a noticeable change in the dust and/or star-forming properties from z~0 to the early Universe, tending towards lower dust temperatures, indicative of strong evolution in the cold dust, `cirrus', component. We show that these objects are potentially the missing link between the well-studied local IR-luminous galaxies, Spitzer IR populations and SCUBA sources -- the z<1 counterparts of the cold z>1 SubMillimetre Galaxies (SMGs) discovered in blank-field submillimetre surveys. The Herschel Space Observatory is well placed to fully characterise the nature of these objects, as its coverage extends over a major part of the far-IR/submm SED for a wide redshift range.

We report the detection of luminous CO(3-2) line emission in the halo of the z=2.6 radio galaxy (HzRG) TXS0828+193, which has no detected counterpart at optical to mid-infrared wavelengths implying a stellar mass < few x10^9 M_sun and relatively low star-formation rates. With the IRAM PdBI we find two CO emission line components at the same position at ~80 kpc distance from the HzRG along the axis of the radio jet, with different blueshifts of few 100 km s^-1 relative to the HzRG and a total luminosity of ~2x10^10 K km s^-1 pc^2 detected at 8 sigma significance. HzRGs have significant galaxy overdensities and extended halos of metal-enriched gas often with embedded clouds or filaments of denser material, and likely trace very massive dark-matter halos. The CO emission may be associated with a gas-rich, low-mass satellite galaxy with little on-going star formation, in contrast to all previous CO detections of galaxies at similar redshifts. Alternatively, the CO may be related to a gas cloud or filament and perhaps jet-induced gas cooling in the outer halo, somewhat in analogy with extended CO emission found in low-redshift galaxy clusters.

Broad absorption line quasars (commonly termed BALQSOs) contain the most dramatic examples of AGN-driven winds. The high absorbing columns in these winds, ~10^24 cm^-2, ensure that BALQSOs are generally X-ray faint. This high X-ray absorption means that almost all BALQSOs have been discovered through optical surveys, and so what little we know about their X-ray properties is derived from very bright optically-selected sources. A small number of X-ray selected BALQSOs (XBALQSOs) have, however, recently been found in deep X-ray survey fields. In this paper we investigate the X-ray and rest-frame UV properties of five XBALQSOs for which we have obtained XMM-Newton EPIC X-ray spectra and deep optical imaging and spectroscopy. We find that, although the XBALQSOs have an alpha_ox steeper by ~0.5 than normal QSOs, their median alpha_ox is nevertheless flatter by 0.30 than that of a comparable sample of optically selected BALQSOs (OBALQSOs). We rule out the possibility that the higher X-ray to optical flux ratio is due to intrinsic optical extinction. We find that the amount of X-ray and UV absorption due to the wind in XBALQSOs is similar, or perhaps greater than, the corresponding wind absorption in OBALQSOs, so the flatter alpha_ox cannot be a result of weaker wind absorption. We conclude that these XBALQSOs have intrinsically higher X-ray to optical flux ratios than the OBALQSO sample with which we compare them.

We present spectroscopy and photometry of the He-rich supernova (SN) 2008ax. The early-time spectra show prominent P-Cygni H lines, which decrease with time and disappear completely about two months after the explosion. In the same period He I lines become the most prominent spectral features. SN 2008ax displays the ordinary spectral evolution of a type IIb supernova. A stringent pre-discovery limit constrains the time of the shock breakout of SN 2008ax to within only a few hours. Its light curve, which peaks in the B band about 20 days after the explosion, strongly resembles that of other He-rich core-collapse supernovae. The observed evolution of SN 2008ax is consistent with the explosion of a young Wolf-Rayet (of WNL type) star, which had retained a thin, low-mass shell of its original H envelope. The overall characteristics of SN 2008ax are reminiscent of those of SN 1993J, except for a likely smaller H mass. This may account for the findings that the progenitor of SN 2008ax was a WNL star and not a K supergiant as in the case of SN 1993J, that a prominent early-time peak is missing in the light curve of SN 2008ax, and that Halpha is observed at higher velocities in SN 2008ax than in SN 1993J.

We present the first mid-infrared Spitzer/Infrared Spectrograph (IRS) observations of powerful radio galaxies at z>2. These radio galaxies, 4C +23.56 (z=2.48) and 6C J1908+7220 (z=3.53), both show strong mid-infrared continua, but with 6C J1908+7220 also showing strong PAH emission at rest-frame 6.2 and 7.7um. In 4C+23.56 we see no obvious PAH features above the continuum. The PAH emission in 6C J1908+7220 is the amongst the most distant observed to date and implies that there is a large instantaneous star formation rate (SFR). This is consistent with the strong detection of 6C J1908+7220 at far-IR and sub-mm wavelengths, indicative of large amounts of cold dust, ~10^9Msun. Powerful radio galaxies at lower redshifts tend to have weak or undetectable PAH features and typically have lower far-IR luminosities. In addition, 4C 23.56 shows moderate silicate absorption as seen in less luminous radio galaxies, indicating tau_{9.7um}=0.3+/-0.05. This feature is shifted out of the observed wavelength range for 6C J1908+7220. The correlation of strong PAH features with large amounts of cold dust, despite the presence of a powerful AGN, is in agreement with other recent results and implies that star formation at high redshift is, in some cases at least, associated with powerful, obscured AGN.

We have used several different methods (radio morphology, radio spectral index, mid-IR to radio and near-IR to radio flux density ratios) to discriminate between AGN and SFGs in faint, sub-mJy radio surveys. We find that the latter two methods are the most powerful with current multi-wavelength data, but that future radio surveys with eMERLIN, LOFAR etc. (and ultimately the SKA) will greatly increase the power of the morphology and spectral index methods. As an example of the science possible we derive the IR luminosity density from the radio-selected SFGs using the radio/IR luminosity correlation. We also examine the contribution by luminosity to the total IR luminosity density and find evidence that fraction of LIRGs remains constant or decreases above z=1 while the relative fraction of ULIRGs continues to increase up to z=2.5.

Discerning the exact nature of the sub-mJy radio population has been historically difficult due to the low luminosity of these sources at most wavelengths. Using deep ground based optical follow-up and observations from the Spitzer Space Telescope we are able to disentangle the radio-selected Active Galactic Nuclei (AGN) and Star Forming Galaxy (SFG) populations for the first time in a deep multi-frequency VLA/MERLIN Survey of the 13^H XMM-Newton/Chandra Deep Field. The discrimination diagnostics include radio morphology, radio spectral index, radio/near-IR and mid-IR/radio flux density ratios. We are now able to calculate the extragalactic Euclidean normalised source counts separately for AGN and SFGs. We find that while SFGs dominate at the faintest flux densities and account for the majority of the up-turn in the counts, AGN still make up around one quarter of the counts at ~5 uJy (1.4 GHz). Using radio luminosity as an unobscured star formation rate (SFR) measure we are then able to examine the comoving SFR density of the Universe up to z=3 which agrees well with measures at other wavelengths. We find a rough correlation of SFR with stellar mass for both the sample presented here and a sample of local radio-selected SFGs from the 6df-NVSS survey. This work also confirms the existence of, and provides alternative evidence for, the evolution of distribution of star formation by galaxy mass: ``downsizing''. As both these samples are SFR-selected, this result suggests that there is a maximum SFR for a given galaxy that depends linearly on its stellar mass. The low ``characteristic times'' (inverse specific SFR) of the SFGs in our sample are similar to those of the 6dF-NVSS sample, implying that most of these sources are in a current phase of enhanced star formation.

We present Hubble Space Telescope images, and spectral energy distributions from optical to infrared wavelengths for a sample of six 0.3<z<0.8 type-2 quasars selected in the mid-infrared using data from the Spitzer Space Telescope. All the host galaxies show some signs of disturbance. Most seem to possess dusty, star-forming disks. The disk inclination, estimated from the axial ratio of the hosts, correlates with the depth of the silicate feature in the mid-infrared spectra, implying that at least some of the reddening towards the AGN arises in the host galaxy. The star formation rates in these objects, as inferred from the strengths of the PAH features and far-infrared continuum, range from 3-90 Msun/yr, but are mostly much larger than those inferred from the [OII]3727 emission line luminosity, due to obscuration. Taken together with studies of type-2 quasar hosts from samples selected in the optical and X-ray, this is consistent with previous suggestions that two types of extinction processes operate within the type-2 quasar population, namely a component due to the dusty torus in the immediate environment of the AGN, and a more extended component due to a dusty, star forming disk.

We present a detailed reduction of a mid-infrared 12um (LW10 filter) ISOCAM open time observation performed on the ESO-Sculptor Survey field (Arnouts et al. 1997). A complete catalogue of 142 sources (120 galaxies and 22 stars), detected with high significance (equivalent to 5sigma), is presented above an integrated flux density of 0.24mJy. Star/galaxy separation is performed by a detailed study of colour-colour diagrams. The catalogue is complete to 1mJy and below this flux density the incompleteness is corrected using two independent methods. The first method uses stars and the second uses optical counterparts of the ISOCAM galaxies; these methods yield consistent results. We also apply an empirical flux density calibration using stars in the field. For each star, the 12um flux density is derived by fitting optical colours from a multi-band chi^2 to stellar templates (BaSel-2.0) and using empirical optical-IR colour-colour relations. This article is a companion analysis to Rocca-Volmerange 2007 et al. where the 12um faint galaxy counts are presented and analysed by galaxy type with the evolutionary code PEGASE.3.

Multi-wavelength galaxy number counts provide clues on the nature of galaxy evolution. The interpretation per galaxy type with the new code PEGASE.3 of the 12mu (ISO) and 24mu (Spitzer) faint galaxy counts provide new constraints on the dust and stellar emission. It also reveals the nature of ULIRGs (L/Lsun > 10^12) and informs on the cosmic star formation history and the time-scales for mass build-up. We firstly present the faint galaxy counts at 12um derived from the catalogue of the ISO-ESO-Sculptor Survey (ISO-ESS) in a companion paper (Seymour et al. 2007) which go down to 0.24 mJy after corrections for incompleteness. We check that they are consistent with ISO number counts at 15um. Secondly we show that the ``normal'' scenarios which fit the deep UV-optical-near-IR counts with PEGASE.2 (Fioc et al. 1999), are unsuccessful in modelling the strong excess simultaneously observed in the cumulative and differential counts at 12um, 15um and 24um. Based on observed 12mu and 25mu IRAS luminosity functions and optical/mid-IR colors, we finally succeed in modelling cumulative and differential counts by only changing 9% of normal galaxies (1/3 of the ellipticals) into ultra-bright dusty elliptical galaxies, interpreted as ULIRGs. This ULIRG population has similarities with high-z radio-galaxy hosts. No number density evolution is included. The Herschel observatory will hopefully confirm these results.

We present the results of a deep 610 MHz survey of the 1^H XMM/Chandra survey area with the GMRT. The resulting maps have a resolution of ~7 arcsec and an rms noise limit of 60 microJy. To a 5 sigma detection limit of 300 microJy we detect 223 sources within a survey area of diameter 64 arcmin. We compute the 610 MHz source counts and compare them to those measured at other radio wavelengths. The well know flattening of the Euclidean-normalised 1.4 GHz source counts below ~2 mJy, usually explained by a population of starburst galaxies undergoing luminosity evolution, is seen at 610 MHz. The 610 MHz source counts can be modelled by the same populations that explain the 1.4 GHz source counts, assuming a spectral index of -0.7 for the starburst galaxies and the steep spectrum AGN population. We find a similar dependence of luminosity evolution on redshift for the starburst galaxies at 610 MHz as is found at 1.4 GHz (i.e. 'Q'= 2.45 (+0.3,-0.4)).

We present the results of a comprehensive Spitzer survey of 69 radio galaxies across 1<z<5.2. Using IRAC (3.6-8.0um), IRS (16um) and MIPS (24-160um) imaging, we decompose the rest-frame optical to infrared spectral energy distributions into stellar, AGN, and dust components and determine the contribution of host galaxy stellar emission at rest-frame H-band. Stellar masses derived from rest-frame near-IR data, where AGN and young star contributions are minimized, are significantly more reliable than those derived from rest-frame optical and UV data. We find that the fraction of emitted light at rest-frame H-band from stars is >60% for ~75% the high redshift radio galaxies. As expected from unified models of AGN, the stellar fraction of the rest-frame H-band luminosity has no correlation with redshift, radio luminosity, or rest-frame mid-IR (5um) luminosity. Additionally, while the stellar H-band luminosity does not vary with stellar fraction, the total H-band luminosity anti-correlates with the stellar fraction as would be expected if the underlying hosts of these radio galaxies comprise a homogeneous population. The resultant stellar luminosities imply stellar masses of 10^{11-11.5}Msun even at the highest redshifts. Powerful radio galaxies tend to lie in a similar region of mid-IR color-color space as unobscured AGN, despite the stellar contribution to their mid-IR SEDs at shorter-wavelengths. The mid-IR luminosities alone classify most HzRGs as LIRGs or ULIRGs with even higher total-IR luminosities. As expected, these exceptionally high mid-IR luminosities are consistent with an obscured, highly-accreting AGN. We find a weak correlation of stellar mass with radio luminosity.

We present the results of a radio polarimetric study of the high-redshift radio galaxy PKS B0529-549 (z=2.575), based on high-resolution 12 mm and 3 cm images obtained with the Australia Telescope Compact Array (ATCA). The source is found to have a rest-frame Faraday rotation measure of -9600 rad m^{-2}, the largest seen thus far in the environment of a z > 2 radio galaxy. In addition, the rest-frame Faraday dispersion in the screen responsible for the rotation is calculated to be 5800 rad m^{-2}, implying rotation measures as large as -15400 rad m^{-2}. Using supporting near-IR imaging from the Very Large Telescope (VLT), we suggest that the rotation measure originates in the Ly-alpha halo surrounding the host galaxy, and estimate the magnetic field strength to be ~10 microGauss. We also present a new optical spectrum of PKS B0529-549 obtained with the New Technology Telescope (NTT), and propose that the emission-line ratios are best described by a photoionization model. Furthermore, the host galaxy is found to exhibit both hot dust emission at 8.0 microns and significant internal visual extinction (~1.6 mag), as inferred from Spitzer Space Telescope near/mid-IR imaging.

Discerning the exact nature of the faint (sub-mJy) radio population has been historically difficult due to the low luminosity of these sources at most wavelengths. Using deep observations from Chandra/XMM-Newton/Spitzer and ground based follow up we are able to disentangle the AGN and star-forming populations for the first time in a deep multi-frequency GMRT/VLA/MERLIN Survey. The many diagnostics include radio luminosity, morphology, radio to mid-IR flux density ratios, radio to optical flux density ratios and radio spectral indices. Further diagnostics, e.g. optical spectra X-ray spectra/hardness ratios, IR colours indicate the presence of the AGN {\em independent} of whether the radio emission is powered by AGN or star-formation. We are able to examine the star-formation history of the universe up to z=2.5 in a unique way based on an unbiased star-formation rate indicator, radio luminosity. This work provides an alternative perspective on the distribution of star-formation by mass, ``downsizing'' and allows us to examine the prevalence of AGN in star-bursts.

We present the results of a comprehensive Spitzer survey of 70 radio galaxies across 1<z<5.2. Using IRAC, IRS and MIPS imaging we determine the rest-frame AGN contribution to the stellar emission peak at 1.6um. The stellar luminosities are found to be consistent with that of a giant elliptical with a stellar mass of 10^11-12Msun. The mean stellar mass remains constant at \~10^11.5Msun up to z=3 indicating that the upper end of the mass function is already in place by this redshift. The mid-IR luminosities imply bolometric IR luminosities that would classify all sources as ULIRGs. The mid-IR to radio luminosity generally correlate implying a common origin for these emissions. The ratio is higher than that found for lower redshift, ie z<1, radio galaxies.

We present the X-ray spectra of 86 optically-identified sources in the 13H XMM-Newton/Chandra deep field which have >70 X-ray counts. The sample consists of 50 broad line AGN, 25 narrow emission line galaxies, 6 absorption line galaxies, and 5 Galactic stars. The majority (42/50) of the broad line AGN have X-ray spectra which are consistent with a power law shape. They have a mean photon index of Gamma = 2.0 +- 0.1 and an intrinsic dispersion sigma = 0.4 +- 0.1. Five BLAGN show a deficit of soft X-rays, indicating absorption. Significant absorption is more common in the narrow emission line galaxies (13/25) and absorption line galaxies (2/6) than in the broad line AGN (5/50), but is not universal in any of these classes of object. The majority of the 20 absorbed sources have X-ray spectra consistent with a simple cold photoelectric absorber, but 6/20 require more complex models with either an additional component of soft X-ray emitting plasma, or an ionised absorber. Of the 16 galaxies which do not show evidence for X-ray absorption, only 2 objects are likely to be powered by star formation, and both have 2-10 keV X-ray luminosities of <= 10^40 cgs. The X-ray emission in the other 14 unabsorbed NELGs and galaxies is most likely powered by AGN, which are not detected in the optical because they are outshone by their luminous host galaxies. The Galactic stars show multi-temperature thermal spectra which peak between 0.5 and 1 keV. Star/AGN discrimination is possible for 4 of the 5 stars solely from their X-ray spectra.

We present: 1) a kinematic and morphological study of the giant Lya nebula associated with the radio galaxy MRC2104-242 (z = 2.49) based on integral field spectroscopic VIMOS data from VLT; 2) a photometric study of the host (proto?) galaxy based on Spitzer Space Telescope data. The galaxy appears to be embedded in a giant (>120 kpc) gas reservoir that surrounds it completely. The kinematic properties of the nebula suggest that it is a rotating structure, which would imply a lower limit to the dynamical mass of ~3 x 1e11 Msun. An alternate scenario is that the gas is infalling. Such a process would be able to initiate and sustain significant central starburst activity, although it is likely to contribute with less than 10% of the total stellar mass of the galaxy. The near- to mid-IR spectral energy distribution of the radio galaxy suggests the existence of a reddened, E(B-V)=0.4(+/-)0 .1$, evolved stellar population of age >~ 1.8 Gyr and mass 5(+/-)2 x 1e11 Msun. The implied formation redshift is z_f >~ 6. This stellar mass is similar to the stellar masses found for massive early-type galaxies at z~2 in deep, near-infrared surveys.

Despite the plethora of deep (sub-mJy) radio surveys there remains considerable doubt as to the exact nature of the galaxies contributing to the source counts. Current evidence suggests that starformation in moderately luminous 'normal' galaxies is responsible for the bulk of the emission below 1mJy. However given the sensitivities of these surveys we would expect a fraction of these sources to be distant radio galaxies. Using deep VLA and GMRT data we have found ~20 high-z candidate radio galaxies in two fields using the classical ultra-steep radio spectrum technique (De Breuck et al., 2000) and selecting galaxies with faint (i'>25) optical counterparts. Several of these sources have X-ray detections in our deep XMM/Chandra observations and have fluxes high enough to put them in the quasar regime if they lie above redshift 3. Recently performed Spitzer GTO observations and upcoming near-infrared observations will help reveal the nature of these sources.

We present the results of a deep X-ray survey conducted with XMM-Newton, centred on the UK ROSAT 13H deep field area. This region covers 0.18 deg^2 and is the first of two areas covered with XMM-Newton as part of an extensive multi-wavelength survey designed to study the nature and evolution of the faint X-ray source population. We have produced detailed Monte-Carlo simulations to obtain a quantitative characterisation of the source detection procedure and to assess the reliability of the resultant sourcelist. We use the simulations to establish a likelihood threshold above which we expect less than 7 (3%) of our sources to be spurious. We present the final catalogue of 225 sources. Within the central 9 arcmin, 68 per cent of source positions are accurate to 2 arcsec, making optical follow-up relatively straightforward. We construct the N(>S) relation in four energy bands: 0.2-0.5 keV, 0.5-2 keV, 2-5 keV and 5-10 keV. In all but our highest energy band we find that the source counts can be represented by a double powerlaw with a bright end slope consistent with the Euclidean case and a break around 10^-14 cgs. Below this flux the counts exhibit a flattening. Our source counts reach densities of 700, 1300, 900 and 300 deg^-2 at fluxes of 4.1x10^-16, 4.5x10^-16, 1.1x10^-15 and 5.3x10^-15 cgs in the 0.2-0.5, 0.5-2, 2-5 and 5-10 keV energy bands respectively. We have compared our source counts with those in the two Chandra deep fields and Lockman hole and find our source counts to be amongst the highest of these fields in all energy bands. We resolve >51% (>50%) of the X-ray background emission in the 1-2 keV (2-5 keV) energy bands.

The up-turn in Euclidean normalised source counts below 1mJy at 1.4GHz is well established in many deep radio surveys. There are strong reasons, observationally and theoretically, to believe that this up-turn is due to strong evolution of the starforming population up to z=2. However this hypothesis needs further confirmation spectroscopically and the examples in the literature are sparse. Theoretically the up-turn is well modelled by the evolution of the local radio starforming population and is consistent with the up-turn seen in recent mid-infrared source counts at 15um (ISOCAM) and 24um (Spitzer) and the tight correlation of the radio and MIR Luminosities of starforming galaxies.

In order to determine the relationship between the faint X-ray and faint radio source populations, and hence to help understand the X-ray and radio emission mechanisms in those faint source populations, we have made a deep 1.4GHz Very Large Array radio survey of the 13h +38deg XMM/ROSAT X-ray Survey Area. From a combined data set of 10hours B configuration data and 14hours A configuration data, maps with 3.35'' resolution and a noise limit of 7.5uJy were constructed. A complete sample of 449 sources was detected within a 30' diameter region above a 4sigma detection limit of 30uJy, at the map centre, making this one of the deepest radio surveys at this frequency. The differential source count shows a significant upturn at sub-milliJansky flux densities, similar to that seen in other deep surveys at 1.4GHz (eg the Phoenix survey), but larger than that seen in the HDF which may have been selected to be under-dense. This upturn is well modelled by the emergence of a population of medium redshift starforming galaxies which dominate at faint flux densities. The brighter source counts are well modelled by AGNs.

We show from the recent extrapolation of the radio-FIR correlation to the MIR that the 20 cm and 15 um differential source counts are likely to come from the same parent population.

As part of a multi-wavelength campaign to observe the 2000--2002 outburst of the prototypical symbiotic star Z Andromedae, we observed this object six times each with the Multi-Element Radio Linked Interferometer Network (MERLIN) and Very Large Array (VLA). The radio flux varied significantly during the course of the optical outburst at all three observation frequencies (1.4, 5, and 15 GHz). A jet-like extension was present in the 2001 September MERLIN image and appeared to be aligned perpendicularly to the plane of the binary orbit. Assuming that the ejection took place at the beginning of the optical outburst, the 0.06 arcsec separation between the peak of the extended emission and the central core implies that the ejected material was moving with a velocity of \~400 km/s. This extended emission faded on a timescale of ~months and was not detected at any other epoch. We consider the implications of jets being a component of a ``prototypical'' symbiotic system and compare properties of the observed jet of Z And with those of the jets in X-ray binary systems.

Our VLA observations of the XMM-Newton/Chandra 13hr deep survey field (see Page et al., this proceedings) result in one of the two deepest 1.4GHz radio maps ever made. Within the 15' radius field covered by the deep X-ray data (0.19 sq deg), a total of 556 radio sources are detected, down to a 4 sigma flux limit of 28uJy. Of the 214 Chandra sources, 55 have radio counterparts. The sub-arcsecond accuracy of the VLA and Chandra positions enable us to determine with high confidence the sources common to both surveys. Here we present the relationship between the X-ray and radio source populations at the faintest radio flux limits yet probed by such a study. We discuss how the X-ray/radio relationship differs as a function of optical morphology, ie between unresolved `stellar' objects and well resolved galaxies. We then discuss the origin of the X-ray and radio emission, ie AGN, starburst or a mixture of both, in these two classes of object.

Of the cosmological gamma-ray bursts, GRB 011121 has the lowest redshift, z=0.36. More importantly, the multi-color excess in the afterglow detected in the Hubble Space Telescope (HST) light curves is compelling observational evidence for an underlying supernova. Here we present near-infrared and radio observations of the afterglow. We undertake a comprehensive modeling of these observations and those reported in the literature and find good evidence favoring a wind-fed circumburst medium. In detail, we infer the progenitor had a mass loss rate of Mdot ~ 10^-7 / v_w3 Mo/yr where v_w3 is the speed of the wind from the progenitor in units of 10^3 km/s. This mass loss rate is similar to that inferred for the progenitor of SN 1998bw which has been associated with GRB 980425. Our data, taken in conjunction with the HST results of Bloom et al. (2002), provide a consistent picture: the long duration GRB 011121 had a massive star progenitor which exploded as a supernova at about the same time as the GRB event.

We describe the first results from a program of multi-wavelength monitoring of symbiotic binary stars in outburst. The prototypical symbiotic star Z Andromedae was observed extensively at radio, optical, far-UV, and X-ray wavelengths during its current (ongoing) outburst. These observations provide evidence for an initial obscuration or the source of ionizing radiation by an optically thick shell of material, outflow with velocities of 100s of km/s, plus collisional excitation and shock heating of the red-giant wind. In addition, the optical light curve reveals variability on time scales as short as days on both the rise to and decline from optical maximum.