Taking advantage of HST CANDELS data, we analyze the lowest redshift (z<0.5) massive galaxies in order to disentangle their structural constituents and study possible faint non-axis-symmetric features. Due to the excellent HST spatial resolution for intermediate-z objects, they are hard to model by purely automatic parametric fitting algorithms. We performed careful single and double S\'ersic fits to their galaxy surface brightness profiles. We also compare the model color profiles with the observed ones and also derive multi-component global effective radii attempting to obtain a better interpretation of the mass-size relation. Additionally, we test the robustness of our measured structural parameters via simulations. We find that the S\'ersic index does not offer a good proxy for the visual morphological type for our sample of massive galaxies. Our derived multi-component effective radii give a better description of the size of our sample galaxies than those inferred from single S\'ersic models with GALFIT. Our galaxy population lays on the scatter of the local mass-size relation, indicating that these massive galaxies do not experience a significant growth in size since z~0.5. Interestingly the few outliers are late-type galaxies, indicating that spheroids must reach the local mass-size relation earlier. For most of our sample galaxies, both single and multi-component S\'ersic models with GALFIT show substantial systematic deviations from the observed SBPs in the outskirts. These residuals may be partly due to several factors, namely a non-optimal data reduction for low surface brightness features, the existence of prominent stellar haloes for massive galaxies and could also arise from conceptual shortcomings of parametric 2D image decomposition tools. They consequently propagate into galaxy color profiles.

We present a detailed study of a high-redshift iron low-ionization broad absorption line (FeLoBAL) quasar (SDSS1214 at $z = 1.046$), including new interferometric $^{12}$CO $J$=2-1 observations, optical through far-infrared photometry, and mid-infrared spectroscopy. The CO line is well-fit by a single Gaussian centered 40 kms$^{-1}$ away from the systemic velocity and implies a total molecular gas mass of $M_\textrm{gas} = 7.3 \times 10^{10} \textrm{M}_\odot$. The infrared SED requires three components: an active galactic nucleus (AGN) torus, an AGN polar dust component, and a starburst. The starburst dominates the infrared emission with a luminosity of log($L_\textrm{SB}[\textrm{L}_\odot]) = 12.91^{+0.02}_{-0.02}$, implying a star formation rate of about 2000 $\textrm{M}_{\odot}$yr$^{-1}$, the highest known among FeLoBAL quasars. The AGN torus and polar dust components are less luminous, at log($L_\textrm{AGN}[\textrm{L}_\odot]) = 12.36^{+0.14}_{-0.15}$ and log($L_\textrm{dust}[\textrm{L}_\odot]) = 11.75^{+0.26}_{-0.46}$, respectively. If all of the molecular gas is used to fuel the ongoing star formation, then the lower limit on the subsequent duration of the starburst is 40 Myr. We do not find conclusive evidence that the AGN is affecting the CO gas reservoir. The properties of SDSS1214 are consistent with it representing the endpoint of an obscured starburst transitioning through a LoBAL phase to that of a classical quasar.

We present the science case for mapping several thousand galaxy (proto)clusters at z=1-10 with a large aperture single dish sub-mm facility, producing a high-redshift counterpart to local large surveys of rich clusters like the well-studied Abell catalogue. Principal goals of a large survey of distant clusters are the evolution of galaxy clusters over cosmic time and the impact of environment on the evolution and formation of galaxies. To make a big leap forward in this emerging research field, the community would benefit from a large-format, wide-band, direct-detection spectrometer (e.g., based on MKID technology), covering a wide field of ~1 square degree and a frequency coverage from 70 to 700 GHz.

We present an exploration of the expected detection of the earliest Active Galactic Nuclei (AGN) in the Universe from state-of-art galaxy formation and evolution semi-analytic models and hydro-dynamical simulations. We estimate the number and radiative characteristics of Super Massive Black Holes (SMBHs) at $z\geq 6$, a redshift range that will be intensively explored by the next generation of telescopes, in particular in the radio through the Square Kilometre Array (SKA) and at high energies with ESA's Athena X-ray Observatory. We find that Athena will be able to observe over 5000 AGN/$\rm deg^2$ at the Epoch of Re-ionization (EoR), $6\leq z \leq 10$. Similarly, for the same redshift range the models/simulations suggest that SKA will detect at least 400 AGN/$\rm deg^2$. Additionally, we stress the importance of the volume of the simulation box as well as the initial physical conditions of the models/simulations on their effect on the luminosity functions (LFs) and the creation of the most massive SMBHs that we currently observe at the EoR. Furthermore, following the evolution of the accretion mode of the SMBHs in each model/simulation, we show that, while the quasar dominates over the radio mode at the EoR, detection at radio wavelengths still reaches significant numbers even at the highest redshifts. Finally, we present the effect that the radiative efficiency has on the LFs by comparing results produced with a constant value for the radiative efficiency and more complex calculations based on the spin of each SMBH.

We present a catalog of 26 faint submillimeter galaxies (SMGs) in the XMM-LSS field identified by cross-matching serendipitously detected sources in archival ALMA Band 6 and 7 data with multi-band near-infrared (NIR) and optical data from the Spitzer Extragalactic Representative Volume Survey, the VISTA Deep Extragalactic Survey, the Canada-France-Hawaii Telescope Legacy Large Survey, and the Hyper Suprime-Cam Subaru Strategic Program. Of the 26 SMGs in our sample, 15 are identified here for the first time. The majority of the sources in our sample (16/26) have faint submm fluxes ($0.1\,{\rm mJy} < S_{\rm 1\,mm} < 1\,{\rm mJy}$). In addition to the 26 SMGs with multi-band optical and NIR detections, there are 60 highly-reliable ($>5\sigma$) ALMA sources with no counterpart in any other band down to an IRAC [4.5] $AB$ magnitude of $\approx 23.7$. To further characterize the 26 galaxies with both ALMA and optical/NIR counterparts, we provide 13-band forced photometry for the entire catalog using the Tractor and calculate photometric redshifts and rest-frame colors. The median redshift of our sample is $\langle z \rangle = 2.66$. We find that our sample galaxies have bluer colors compared to bright SMGs, and the UVJ color plot indicates that their colors are consistent with main sequence star-forming galaxies. Our results provide new insights into the nature of the faint population of SMGs, and also highlight opportunities for galaxy evolution studies based on archival ALMA data.

A fraction of Narrow Line Seyfert 1 galaxies (NLS1) are hosted by galaxies that present a disturbed morphology, in some cases hinting for merger processes, that are putative sources of gas replenishment. We have been investigating the poorly studied population of radio loud NLS1 (RL-NLS1) showing a flat radio spectrum, assumed to be the manifestation of the presence of a radio jet. In some of the objects the infrared emission is well fitted by a combination of an AGN component and an "active" host galaxy component like M82, the estimate SFR being in the LIRG/ULIRG range (10-500 Msun/year). In order to better characterize that component, we have been investigating the sub-millimeter/millimeter emission of the sources using APEX. Here we present the results concerning a pilot sample of 2 representative objects.

There is significant scientific value to be gained from combining AKARI fluxes with data at other far-IR wavelengths from the IRAS and Herschel missions. To be able to do this we must ensure that there are no systematic differences between the datasets that need to be corrected before the fluxes are compatible with each other. One such systematic effect identified in the BSCv1 data is the issue of beam corrections. We determine these for the BSCv2 data by correlating ratios of appropriate IRAS and AKARI bands with the difference in 2MASS J band extended and point source magnitudes for sources cross matched between the IRAS FSC, AKARI BSCv2 and 2MASS catalogs. We find significant correlations (p<< 10^-13) indicating that beam corrections are necessary in the 65 and 90 micron bands. We then use these corrected fluxes to supplement existing data in spectral energy distribution (SED) fits for ultraluminous infrared galaxies in the HERUS survey. The addition of AKARI fluxes makes little difference to the results of simple (T,beta) fits to the SEDs of these sources, though there is a general decrease in reduced chi^2 values. The utility of the extra AKARI data, however, is in allowing physically more realistic SED models with more parameters to be fit to the data. We also extend our analysis of beam correction issues in the AKARI data by examining the Herschel Reference Sample, which have Herschel photometry from 100 to 500 microns and which are more spatially extended than the HERUS ULIRGs. Thirty-four of the HRS sources have good Herschel SEDs and matching data from AKARI. This investigation finds that our simple 2MASS-based beam correction scheme is inadequate for these larger and more complex sources. There are also indications that additional beam corrections at 140 and 160 microns are needed for these sources, extended on scales >1 arcminute.

When combined with the huge collecting area of the ELT, MOSAIC will be the most effective and flexible Multi-Object Spectrograph (MOS) facility in the world, having both a high multiplex and a multi-Integral Field Unit (Multi-IFU) capability. It will be the fastest way to spectroscopically follow-up the faintest sources, probing the reionisation epoch, as well as evaluating the evolution of the dwarf mass function over most of the age of the Universe. MOSAIC will be world-leading in generating an inventory of both the dark matter (from realistic rotation curves with MOAO fed NIR IFUs) and the cool to warm-hot gas phases in z=3.5 galactic haloes (with visible wavelenth IFUs). Galactic archaeology and the first massive black holes are additional targets for which MOSAIC will also be revolutionary. MOAO and accurate sky subtraction with fibres have now been demonstrated on sky, removing all low Technical Readiness Level (TRL) items from the instrument. A prompt implementation of MOSAIC is feasible, and indeed could increase the robustness and reduce risk on the ELT, since it does not require diffraction limited adaptive optics performance. Science programmes and survey strategies are currently being investigated by the Consortium, which is also hoping to welcome a few new partners in the next two years.

We present images taken using the Gemini South Adaptive Optics Imager (GSAOI) with the Gemini Multiconjugate Adaptive Optics System (GeMS) in three 2 arcmin$^2$ fields in the Spitzer Extragalactic Representative Volume Survey. These GeMS/GSAOI observations are among the first $\approx 0.1^{''}$ resolution data in the near-infrared spanning extragalactic fields exceeding $1.5^{\prime}$ in size. We use these data to estimate galaxy sizes, obtaining results similar to those from studies with the Hubble Space Telescope, though we find a higher fraction of compact star forming galaxies at $z>2$. To disentangle the star-forming galaxies from active galactic nuclei (AGN), we use multiwavelength data from surveys in the optical and infrared, including far-infrared data from Herschel, as well as new radio continuum data from the Australia Telescope Compact Array and Very Large Array. We identify ultraluminous infrared galaxies (ULIRGs) at $z \sim 1-3$, which consist of a combination of pure starburst galaxies and Active Galactic Nuclei (AGN)/starburst composites. The ULIRGs show signs of recent merger activity, such as highly disturbed morphologies and include a rare candidate triple AGN. We find that AGN tend to reside in hosts with smaller scale sizes than purely star-forming galaxies of similar infrared luminosity. Our observations demonstrate the potential for MCAO to complement the deeper galaxy surveys to be made with the James Webb Space Telescope.

We present the VIMOS Spectroscopic Survey of a Supercluster in the COSMOS field (VIS3COS) at z~0.84. We use VIMOS high-resolution spectra (GG475 filter) to spectroscopically select 490 galaxies in and around the super-structure and an additional 481 galaxies in the line of sight. We present the redshift distribution, the catalog to be made public and the first results on the properties of individual galaxies and stacked spectra (3500{\AA}<\lambda<4200{\AA} rest-frame). We probe a wide range of densities and environments (from low-density field to clusters and rich groups). We find a decrease in the median star formation rate from low to high density environments in all bins of stellar mass and a sharp rise of the quenched fraction (from ~10% to ~40-60%) of intermediate stellar mass galaxies $(10<\log_{10}\left(M_\star/\mathrm{M_\odot}\right)<10.75)$ from filaments to clusters. The quenched fraction for massive galaxies shows little dependence on the environment being constant at ~30-40%. We find a break in the strength of the [OII] emission, with nearly constant line equivalent widths at lower densities (~-11{\AA}) and then a drop to ~-2.5{\AA} towards higher densities. The break in the [OII] line strength happens at similar densities ($\log_{10}(1+\delta)\sim0.0-0.5$) as the observed rise in the quenched fraction. Our results may provide further clues regarding the different environmental processes affecting galaxies with different stellar masses and highlight the advantages of a single dataset in the COSMOS field probing a wide range of stellar masses and environments. We hypothesize that quenching mechanisms are enhanced in high-density regions.

We report molecular gas observations of IRAS 20100-4156 and IRAS 03158+4227, two local ultraluminous infrared galaxies (ULIRGs) hosting some of the fastest and most massive molecular outflows known. Using ALMA and PdBI observations, we spatially resolve the CO(1-0) emission from the outflowing molecular gas in both and find maximum outflow velocities of $ v_{\rm max} \sim 1600$ and $\sim 1700$ km/s for IRAS 20100-4156 and IRAS 03158+4227, respectively. We find total gas mass outflow rates of $\dot M_{\rm OF} \sim 670$ and $\sim 350$ Msun/yr, respectively, corresponding to molecular gas depletion timescales $\tau^{\rm dep}_{\rm OF} \sim 11$ and $\sim 16$ Myr. This is nearly 3 times shorter than the depletion timescales implied by star formation, $\tau^{\rm dep}_{\rm SFR} \sim 33$ and $\sim 46$ Myr, respectively. To determine the outflow driving mechanism, we compare the starburst ($L_{*}$) and AGN ($L_{\rm AGN}$) luminosities to the outflowing energy and momentum fluxes, using mid-infrared spectral decomposition to discern $L_{\rm AGN}$. Comparison to other molecular outflows in ULIRGs reveals that outflow properties correlate similarly with $L_{*}$ and $L_{\rm IR}$ as with $L_{\rm AGN}$, indicating that AGN luminosity alone may not be a good tracer of feedback strength and that a combination of AGN and starburst activity may be driving the most powerful molecular outflows. We also detect the OH 1.667 GHz maser line from both sources and demonstrate its utility in detecting molecular outflows.

We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local Ultraluminous Infrared Galaxies (ULIRGs), selected at 60${\mu}$m by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350 and 500${\mu}$m. We describe these observations, present the results, and combine the new observations with data from IRAS to examine the far-IR spectral energy distributions (SEDs) of these sources. We fit the observed SEDs of HERUS objects with a simple parameterised modified black body model where temperature and emissivity $\beta$ are free parameters. We compare the fitted values to those of non-ULIRG local galaxies, and find, in agreement with earlier results, that HERUS ULIRGs have warmer dust (median temperature T = 37.9+/-4.7 K compared to 21.3+/-3.4 K) but a similar $\beta$ distribution (median $\beta$ = 1.7 compared to 1.8) to the Herschel reference sample (HRS, Cortese et al., 2014) galaxies. Dust masses are found to be in the range of 10^7.5 to 10^9 Msun significantly higher than that of Herschel Reference Sample (HRS) sources. We compare our results for local ULIRGs with higher redshift samples selected at 250 and 850${\mu}$m. These latter sources generally have cooler dust and/or redder 100-to-250 ${\mu}$m colours than our 60${\mu}$m-selected ULIRGs. We show that this difference may in part be the result of the sources being selected at different wavelengths rather than being a simple indication of rapid evolution in the properties of the population.

Combining the catalogue of galaxy morphologies in the COSMOS field and the sample of H$\alpha$ emitters at redshifts $z=0.4$ and $z=0.84$ of the HiZELS survey, we selected $\sim$ 220 star-forming bulgeless systems (S\'ersic index $n \leq 1.5$) at both epochs. We present their star formation properties and we investigate their contribution to the star formation rate function (SFRF) and global star formation rate density (SFRD) at $z < 1$. For comparison, we also analyse H$\alpha$ emitters with more structurally evolved morphologies that we split into two classes according to their S\'ersic index $n$: intermediate ($ 1.5 < n \leq 3 $) and bulge-dominated ($n > 3$). At both redshifts the SFRF is dominated by the contribution of bulgeless galaxies and we show that they account for more than 60% of the cosmic SFRD at $z < 1$. The decrease of the SFRD with redshift is common to the three morphological types but it is stronger for bulge-dominated systems. Star-forming bulgeless systems are mostly located in regions of low to intermediate galaxy densities ($\Sigma \sim 1 - 4$ Mpc$^{-2}$) typical of field-like and filament-like environments and their specific star formation rates (sSFRs) do not appear to vary strongly with local galaxy density. Only few bulgeless galaxies in our sample have high (sSFR $>$ 10$^{-9}$ yr$^{-1}$) and these are mainly low-mass systems. Above $M_* \sim 10^{10}$ M$_{\odot}$ bulgeless are evolving at a "normal" rate (10$^{-9}$ yr$^{-1} <$ sSFR $<$10$^{-10}$ yr$^{-1}$) and in the absence of an external trigger (i.e. mergers/strong interactions) they might not be able to develop a central classical bulge.

We discuss the ground-breaking science that will be possible with a wide area survey, using the MeerKAT telescope, known as MeerKLASS (MeerKAT Large Area Synoptic Survey). The current specifications of MeerKAT make it a great fit for science applications that require large survey speeds but not necessarily high angular resolutions. In particular, for cosmology, a large survey over $\sim 4,000 \, {\rm deg}^2$ for $\sim 4,000$ hours will potentially provide the first ever measurements of the baryon acoustic oscillations using the 21cm intensity mapping technique, with enough accuracy to impose constraints on the nature of dark energy. The combination with multi-wavelength data will give unique additional information, such as exquisite constraints on primordial non-Gaussianity using the multi-tracer technique, as well as a better handle on foregrounds and systematics. Such a wide survey with MeerKAT is also a great match for HI galaxy studies, providing unrivalled statistics in the pre-SKA era for galaxies resolved in the HI emission line beyond local structures at z > 0.01. It will also produce a large continuum galaxy sample down to a depth of about 5\,$\mu$Jy in L-band, which is quite unique over such large areas and will allow studies of the large-scale structure of the Universe out to high redshifts, complementing the galaxy HI survey to form a transformational multi-wavelength approach to study galaxy dynamics and evolution. Finally, the same survey will supply unique information for a range of other science applications, including a large statistical investigation of galaxy clusters as well as produce a rotation measure map across a huge swathe of the sky. The MeerKLASS survey will be a crucial step on the road to using SKA1-MID for cosmological applications and other commensal surveys, as described in the top priority SKA key science projects (abridged).

We present new results on the evolution of rest-frame blue/UV sizes and Sersic indices of H$\alpha$-selected star-forming galaxies over the last 11 Gyrs. We investigate how the perceived evolution can be affected by a range of biases and systematics such as cosmological dimming and resolution effects. We use GALFIT and an artificial redshifting technique, which includes the luminosity evolution of H$\alpha$-selected galaxies, to quantify the change on the measured structural parameters with redshift. We find typical sizes of 2 to 3 kpc and Sersic indices of n~1.2, close to pure exponential disks all the way from z=2.23 to z=0.4. At z=0 we find typical sizes of 4-5 kpc. Our results show that, when using GALFIT, cosmological dimming has a negligible impact on the derived effective radius for galaxies with <10 kpc, but we find a ~20% bias on the estimate of the median Sersic indices, rendering galaxies more disk-like. Star-forming galaxies have grown on average by a factor of 2-3 in the last 11 Gyrs with $r_e\propto(1+z)^{-0.75}$. By exploring the evolution of the stellar mass-size relation we find evidence for a stronger size evolution of the most massive star-forming galaxies since z~2, as they grow faster towards z~0 when compared to the lower stellar mass counterparts. As we are tracing the rest-frame blue/UV, we are likely witnessing the growth of disks where star formation is ongoing in galaxies while their profiles remain close to exponential disks, n<1.5, across the same period.

We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectra. Special treatment in the data reduction has been given to any observation suffering from artefacts in the data caused by anomalous instrumental effects to improve the final spectra. Complete spectra are shown covering $200 - 671\mu$m with photometry in the SPIRE bands at 250$\mu$m, 350$\mu$m and 500$\mu$m. The spectra include near complete CO ladders for over half of our sample, as well as fine structure lines from [CI] 370 $\mu$m, [CI] 609 $\mu$m, and [NII] 205 $\mu$m. We also detect H$_{2}$O lines in several objects. We construct CO Spectral Line Energy Distributions (SLEDs) for the sample, and compare their slopes with the far-infrared colours and luminosities. We show that the CO SLEDs of ULIRGs can be broadly grouped into three classes based on their excitation. We find that the mid-J (5$<$J$<$8) lines are better correlated with the total far-infrared luminosity, suggesting that the warm gas component is closely linked to recent star-formation. The higher J transitions do not linearly correlate with the far-infrared luminosity, consistent with them originating in hotter, denser gas unconnected to the current star-formation. {\bf We conclude that in most cases more than one temperature components are required to model the CO SLEDs.}

In recent years there has been a paradigm shift from centralised to geographically distributed resources. Individual entities are no longer able to host or afford the necessary expertise in-house, and, as a consequence, society increasingly relies on widespread collaborations. Although such collaborations are now the norm for scientific projects, more technical structures providing support to a distributed scientific community without direct financial or other material benefits are scarce. The network of European ALMA Regional Centre (ARC) nodes is an example of such an internationally distributed user support network. It is an organised effort to provide the European ALMA user community with uniform expert support to enable optimal usage and scientific output of the ALMA facility. The network model for the European ARC nodes is described in terms of its organisation, communication strategies and user support.

We combine a deep 0.5~deg$^2$, 1.4~GHz deep radio survey in the Lockman Hole with infrared and optical data in the same field, including the SERVS and UKIDSS near-infrared surveys, to make the largest study to date of the host galaxies of radio sources with typical radio flux densities $\sim 50 \;\mu$Jy. 87% (1274/1467) of radio sources have identifications in SERVS to $AB\approx 23.1$ at 3.6 or 4.5$\mu$m, and 9% are blended with bright objects (mostly stars), leaving only 4% (59 objects) which are too faint to confidently identify in the near-infrared. We are able to estimate photometric redshifts for 68% of the radio sources. We use mid-infrared diagnostics to show that the source population consists of a mixture of star forming galaxies, rapidly accreting (cold mode) AGN and low accretion rate, hot mode AGN, with neither AGN nor starforming galaxies clearly dominating. We see the breakdown in the $K-z$ relation in faint radio source samples, and show that it is due to radio source populations becoming dominated by sources with radio luminosities $\sim 10^{23}\;{\rm WHz^{-1}}$. At these luminosities, both the star forming galaxies and the cold mode AGN have hosts with stellar luminosities about a factor of two lower than those of hot mode AGN, which continue to reside in only the most massive hosts. We show that out to at least $z\sim 2$, galaxies with stellar masses $>10^{11.5}\, M_{\odot}$ have a radio-loud fraction up to $\sim 30$%. This is consistent with there being a sufficient number of radio sources that radio-mode feedback could play a role in galaxy evolution.

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

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.

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

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.

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.

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.

Understanding the evolution of accretion activity is fundamental to our understanding of how galaxies form and evolve over the history of the Universe. We analyse a complete sample of 27 radio galaxies which includes both high-excitation (HEGs) and low excitation galaxies (LEGs), spanning a narrow redshift range of 0.9 < z < 1.1 and covering a factor of ~1000 in radio luminosity. Using data from the Spitzer Space Telescope combined with ground-based optical and near-infrared imaging, we show that the host galaxies have masses in the range of 10.7 < log (M /M_sun) < 12.0 with HEGs and LEGs exhibiting no difference in their mass distributions. We also find that HEGs accrete at significantly higher rates than LEGs, with the HEG/LEG division lying at an Eddington ratio of ~0.04, which is in excellent agreement with theoretical predictions of where the accretion rate becomes radiatively inefficient, thus supporting the idea of HEGs and LEGs being powered by different modes of accretion. Our study also shows that at least up to L_151MHz ~3x10^27 W /Hz /sr, HEGs and LEGs are indistinguishable in terms of their radio properties. From this result we infer that, at least for the lower radio luminosity range, another factor besides accretion rate must play an important role in the process of triggering jet activity.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

Over the past 18 months we have revisited the science requirements for a multi-object spectrograph (MOS) for the European Extremely Large Telescope (E-ELT). These efforts span the full range of E-ELT science and include input from a broad cross-section of astronomers across the ESO partner countries. In this contribution we summarise the key cases relating to studies of high-redshift galaxies, galaxy evolution, and stellar populations, with a more expansive presentation of a new case relating to detection of exoplanets in stellar clusters. A general requirement is the need for two observational modes to best exploit the large (>40 sq. arcmin) patrol field of the E-ELT. The first mode ('high multiplex') requires integrated-light (or coarsely resolved) optical/near-IR spectroscopy of >100 objects simultaneously. The second ('high definition'), enabled by wide-field adaptive optics, requires spatially-resolved, near-IR of >10 objects/sub-fields. Within the context of the conceptual study for an ELT-MOS called MOSAIC, we summarise the top-level requirements from each case and introduce the next steps in the design process.

We present the results of a new search for bright star-forming galaxies at z ~ 7 within the UltraVISTA DR2 and UKIDSS UDS DR10 data, which together provide 1.65 sq deg of near-infrared imaging with overlapping optical and Spitzer data. Using a full photo-z analysis to identify high-z galaxies and reject contaminants, we have selected a sample of 34 luminous (-22.7 < M_UV < -21.2) galaxies with 6.5 < z < 7.5. Crucially, the deeper imaging provided by UltraVISTA DR2 confirms all of the robust objects previously uncovered by Bowler et al. (2012), validating our selection technique. Our sample includes the most massive galaxies known at z ~ 7, with M_* ~ 10^{10} M_sun, and the majority are resolved, consistent with larger sizes (r_{1/2} ~ 1 - 1.5 kpc) than displayed by less massive galaxies. From our final sample, we determine the form of the bright end of the rest-frame UV galaxy luminosity function (LF) at z ~ 7, providing strong evidence that the bright end of the z = 7 LF does not decline as steeply as predicted by the Schechter function fitted to fainter data. We consider carefully, and exclude the possibility that this is due to either gravitational lensing, or significant contamination of our galaxy sample by AGN. Rather, our results favour a double power-law form for the galaxy LF at high z, or, more interestingly, a LF which simply follows the form of the dark-matter halo mass function at bright magnitudes. This suggests that the physical mechanism which inhibits star-formation activity in massive galaxies (i.e. AGN feedback or some other form of `mass quenching') has yet to impact on the observable galaxy LF at z ~ 7, a conclusion supported by the estimated masses of our brightest galaxies which have only just reached a mass comparable to the critical `quenching mass' of M_* = 10 ^{10.2} M_sun derived from studies of the mass function of star-forming galaxies at lower z.

We present a catalogue of bulgeless galaxies, which includes 19225 objects selected in four of the deepest, largest multi-wavelength datasets available -- COSMOS, AEGIS, GEMS and GOODS -- at intermediate redshift ($0.4 \leq z \leq 1.0$). The morphological classification was provided by the Advanced Camera for Surveys General Catalogue (ACS-GC), which used publicly available data obtained with the ACS instrument on the Hubble Space Telescope. Rest-frame photometric quantities were derived using kcorrect. We analyse the properties of the sample and the evolution of pure-disc systems with redshift. Very massive [$\log (M_\star/M_{\odot}) > 10.5$] bulgeless galaxies contribute to ~30% of the total galaxy population number density at $z \geq 0.7$, but their number density drops substantially with decreasing redshift. We show that only a negligible fraction of pure discs appear to be quiescent systems, and red sequence bulgeless galaxies show indications of dust-obscured star formation. X-ray catalogues were used to search for X-ray emission within our sample. After visual inspection and detailed parametric morphological fitting we identify 30 AGN that reside in galaxies without a classical bulge. The finding of such peculiar objects at intermediate redshift shows that while AGN growth in merger-free systems is a rare event (0.2% AGN hosts in this sample of bulgeless galaxies), it can indeed happen relatively early in the Universe history.

A VLA Sky Survey of the extragalactic sky at S band (2-4 GHz) with polarization information can uniquely probe the magneto-ionic medium in a wide range of astrophysical environments over cosmic time. For a shallow all-sky survey, we expect to detect over 4 million sources in total intensity $>$ 0.45 mJy beam$^{-1}$ and over 2.2$\times$10$^5$ sources in polarized intensity. With these new observations, we expect to discover new classes of polarized radio sources in very turbulent astrophysical environments and those with extreme values of Faraday depth. Moreover, by determining reliable Faraday depths and by modeling depolarization effects, we can derive properties of the magneto-ionic medium associated with AGNs, absorption line systems and galaxies, addressing the following unresolved questions: (1) What is the covering fraction, the degree of turbulence and the origin of absorption line systems? (2) What is the thermal content in AGNs and radio galaxies? (3) How do AGNs and galaxies evolve over cosmic time? (4) What causes the increase in percentage polarization with decreasing flux densities at the low flux density end of the polarized source count? (5) What is the growth rate of large-scale magnetic fields in galaxies?

AIMS. This work addresses the AGN IR-selection dependency on intrinsic source luminosity and obscuration, in order to identify and charaterise biases which could affect conclusions in studies. METHODS. We study IR-selected AGN in the Chandra Deep Field South (CDFS) survey and in the Cosmological Survey (COSMOS). The AGN sample is divided into low and high X-ray luminosity classes, and into unobscured (type-1) and obscured (type-2) classes by means of X-ray and optical spectroscopy data. Specifically in the X-rays regime, we adopt the intrinsic luminosity taking into account the estimated column density (N_H). We also take the chance to highlight important differences resulting from adopting different methods to assess AGN obscuration. RESULTS. In agreement with previous studies, we also find that AGN IR-selection efficiency shows a decrease with decreasing source AGN X-ray luminosity. For the intermediate-luminosity AGN population (43.3<~log(L_X[erg/s])<~44$), the efficiency also worsens with increasing obscuration. The same sample also shows an evolution with cosmic time of the obscured fraction at the highest X-ray luminosities, independently of adopted type-1/type-2 classification method. CONCLUSIONS. We confirm that AGN IR-selection is genuinely biased toward unobscured AGN, but only at intermediate luminosities. At the highest luminosities, where AGN IR-selection is more efficient, there is no obscuration bias. We show that type-1 AGN are intrinsically more luminous than type-2 AGN only at z<~1.6, thus resulting in more type-1 AGN being selected the shallower the IR survey is. [ABRIDGED: Conclusion on survey strategy. See manuscript.] Finally, when the James Webb Space Telescope comes online, the broad-band filters 2.0um, 4.4um, 7.7um, and 18um will be key to disentangle AGN from non-AGN dominated SEDs at depths where spectroscopy becomes impractical.

We present a statistical method based on a maximum likelihood approach to constrain the number counts of extragalactic sources below the nominal flux-density limit of continuum imaging surveys. We extract flux densities from a radio map using positional information from an auxiliary catalogue and show that we can model the number counts of this undetected population down to flux density levels well below the detection threshold of the radio survey. We demonstrate the capabilities that our method will have with future generation wide-area radio surveys by performing simulations over various sky areas. We show that it is possible to accurately constrain the number counts of the simulated distribution down to one tenth of the flux noise rms with just a sky area of 100 deg$^2$. We then test the application of our method using data from the Faint Images of the Radio Sky at Twenty-Centimeters survey (FIRST). We extract flux densities from the FIRST map, sensitive to 150 $\mu$Jy/beam (1 $\sigma$), using the positional information from a catalogue in the same field, also acquired at the same frequency, sensitive to 12 $\mu$Jy/beam (1 $\sigma$). Implementing our method, with known source positions, we are able to recover the right differential number counts of the noise-dominated FIRST map fluxes down to a flux density level which is one-tenth the FIRST detection threshold.

We present Herschel photometry and spectroscopy, carried out as part of the Herschel ULIRG survey (HERUS), and a model for the infrared to submillimetre emission of the ultraluminous infrared galaxy IRAS 08572+3915. This source shows one of the deepest known silicate absorption features and no polycyclic aromatic hydrocarbon (PAH) emission. The model suggests that this object is powered by an active galactic nucleus (AGN) with a fairly smooth torus viewed almost edge-on and a very young starburst. According to our model the AGN contributes about 90% of the total luminosity of 1.1 x 10^13 Lo, which is about a factor of five higher than previous estimates. The large correction of the luminosity is due to the anisotropy of the emission of the best fit torus. Similar corrections may be necessary for other local and high-z analogs. This correction implies that IRAS 08572+3915 at a redshift of 0.05835 may be the nearest hyperluminous infrared galaxy and probably the most luminous infrared galaxy in the local (z < 0.2) Universe. IRAS 08572+3915 shows a low ratio of [CII] to IR luminosity (log L_[CII]/L_{IR} < -3.8) and a [OI]63um to [CII]158um line ratio of about 1 that supports the model presented in this letter.

[Abridged] We study the contribution of hot-dust to the luminosity density of galaxies and its evolution with cosmic time. Using the Spitzer-IRAC data in the COSMOS field, we estimate the contribution from hot-dust at rest-frame 4.2um (from ~0 < z < ~0.2 up to ~0.5 < z < ~0.9). This wavelength corresponds to black-body temperature of ~690K. The contribution due to stellar emission is estimated from the rest-frame 1.6um luminosity (assumed to result from stellar emission alone) and subtracted from the mid-infrared luminosity of galaxies to measure hot-dust emission. In order to attempt the study of the 3.3um-PAH feature, we use the rest-frame 4.2um to infer the hot-dust flux at 3.3um. This study is performed for different spectral types of galaxies: early-type, late-type, starburst, and IR-selected AGN. We find that: (a) the decrease of the hot-dust luminosity density since ~0.5 < z < ~1 is steeper (by at least ~0.5dex) compared to that of the cold-dust, giving support to the scenario where galaxy obscuration increases with redshift, as already proposed in the literature; (b) hot-dust and PAH emission evolution seems to be correlated with stellar mass, where rest-frame 1.6um luminous non-AGN galaxies (i.e., massive systems) show a stronger decrement (with decreasing redshift) in hot-dust and PAH emission than the less luminous (less massive) non-AGN galaxies; (c) despite comprising < ~3% of the total sample, AGN contribute as much as a third to the hot-dust luminosity density at z < 1 and clearly dominate the bright-end of the total hot-dust Luminosity Density Function at ~0.5 < z < ~0.9; (d) the average dust-to-total luminosity ratio increases with redshift, while PAH-to-total luminosity ratio remains fairly constant; (e) at M_1.6 > -25, the dust-to-total and PAH-to-total luminosity ratios increase with decreasing luminosity, but deeper data is required to confirm this result.

We present Herschel observations of six fine-structure lines in 25 Ultraluminous Infrared Galaxies at z<0.27. The lines, [O III]52, [N III]57, [O I]63, [N II]122, [O I]145, and [C II]158, are mostly single gaussians with widths <600 km s-1 and luminosities of 10^7 - 10^9 Solar. There are deficits in the [O I]63/L_IR, [N II]/L_IR, [O I]145/L_IR, and [C II]/L_IR ratios compared to lower luminosity systems. The majority of the line deficits are consistent with dustier H II regions, but part of the [C II] deficit may arise from an additional mechanism, plausibly charged dust grains. This is consistent with some of the [C II] originating from PDRs or the ISM. We derive relations between far-IR line luminosities and both IR luminosity and star formation rate. We find that [N II] and both [O I] lines are good tracers of IR luminosity and star formation rate. In contrast, [C II] is a poor tracer of IR luminosity and star formation rate, and does not improve as a tracer of either quantity if the [C II] deficit is accounted for. The continuum luminosity densities also correlate with IR luminosity and star formation rate. We derive ranges for the gas density and ultraviolet radiation intensity of 10^1 < n < 10^2.5 and 10^2.2 < G_0 < 10^3.6, respectively. These ranges depend on optical type, the importance of star formation, and merger stage. We do not find relationships between far-IR line properties and several other parameters; AGN activity, merger stage, mid-IR excitation, and SMBH mass. We conclude that these far-IR lines arise from gas heated by starlight, and that they are not strongly influenced by AGN activity.

We report on our observations of the 79 and 119um doublet transitions of OH for 24 local (z<0.262) ULIRGs observed with Herschel-PACS as part of the Herschel ULIRG Survey (HERUS). Some OH119 profiles display a clear P-Cygni shape and therefore imply outflowing OH gas, other profiles are predominantly in absorption or are completely in emission. We find that the relative strength of the OH emission component decreases as the silicate absorption increases. This locates the OH outflows inside the obscured nuclei. The maximum outflow velocities for our sources range from less than 100 to 2000 km/s, with 15/24 (10/24) sources showing OH absorption at velocities exceeding 700 km/s (1000 km/s). Three sources show maximum OH outflow velocities exceeding that of Mrk231. Since outflow velocities above 500-700 km/s are thought to require an active galactic nucleus (AGN) to drive them, about 2/3 of our ULIRG sample may host AGN-driven molecular outflows. This finding is supported by the correlation we find between the maximum OH outflow velocity and the IR-derived bolometric AGN luminosity. No such correlation is found with the IR-derived star formation rate. The highest outflow velocities are found among sources which are still deeply embedded. We speculate that the molecular outflows in these sources may be in an early phase of disrupting the nuclear dust veil before these sources evolve into less obscured AGN. Four of our sources show high-velocity wings in their [C II] fine-structure line profiles implying neutral gas outflow masses of at least 2-4.5 x 10^8 Msun.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Observations in the last decade have provided strong evidence that the growth of supermassive black holes at the centres of galaxies is among the most influential processes in galaxy evolution. Open questions that relate to our current understanding of black hole growth and its relation to the build-up of galaxies at redshifts z=1-4, when most black holes and stars we see in present-day galaxies were put in place, include: what is the nature of AGN feedback and whether it plays a significant role in the evolution of galaxies? what is the dominant population of accreting AGN at that critical epoch? is it dominated by obscured objects as required by many current observations and models? The Athena+ mission concept will provide the technological leap required for a breakthrough in our understanding of AGN and galaxy evolution at the heyday of the Universe. The high throughput of Athena+ will allow the systematic study of the incidence, nature and energetics of AGN feedback processes to z~4 via the identification and measurement of blue-shifted X-ray absorption lines. The excellent survey and spectral capabilities of the Athena+ Wide Field Imager will complete the census of black hole growth by yielding samples of up to 100 times larger than is currently possible of the most heavily obscured, including Compton thick, AGN to redshifts z~4. The demographics of this population relative to their hosts is fundamental for understanding how major black hole growth events relate to the build-up of galaxies.

A crucial challenge in astrophysics over the coming decades will be to understand the origins of supermassive black holes (SMBHs) that lie at the centres of most, if not all, galaxies. The processes responsible for the initial formation of these SMBHs and their early growth via accretion - when they are seen as Active Galactic Nuclei (AGN) - remain unknown. To address this challenge, we must identify low luminosity and obscured z>6 AGNs, which represent the bulk of early SMBH growth. Sensitive X-ray observations are a unique signpost of accretion activity, uncontaminated by star formation processes, which prevent reliable AGN identification at other wavelengths (e.g. optical, infrared). The Athena+ Wide Field Imager will enable X-ray surveys to be carried out two orders of magnitude faster than with Chandra or XMM-Newton, opening a new discovery space and identifying over 400 z>6 AGN, including obscured sources. Athena+ will also play a fundamental role to enhance the scientific return of future multiwavelength facilities that will probe the physical conditions within the host galaxies of early SMBHs, which is vital for understanding how SMBHs form, what fuels their subsequent growth, and to assess their impact on the early Universe. Follow-up of samples of z>6 galaxies with the Athena+ X-ray Integral Field Unit could also reveal the presence of highly obscured AGNs, thanks to the detection of strong iron lines. Thus, Athena+ will enable the first quantitative measurements of the extent and distribution of SMBH accretion in the early Universe.

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

MOONS is a new conceptual design for a Multi-Object Optical and Near-infrared Spectrograph for the Very Large Telescope (VLT), selected by ESO for a Phase A study. The baseline design consists of 1000 fibers deployable over a field of view of 500 square arcmin, the largest patrol field offered by the Nasmyth focus at the VLT. The total wavelength coverage is 0.8um-1.8um and two resolution modes: medium resolution and high resolution. In the medium resolution mode (R=4,000-6,000) the entire wavelength range 0.8um-1.8um is observed simultaneously, while the high resolution mode covers simultaneously three selected spectral regions: one around the CaII triplet (at R=8,000) to measure radial velocities, and two regions at R=20,000 one in the J-band and one in the H-band, for detailed measurements of chemical abundances. The grasp of the 8.2m Very Large Telescope (VLT) combined with the large multiplex and wavelength coverage of MOONS - extending into the near-IR - will provide the observational power necessary to study galaxy formation and evolution over the entire history of the Universe, from our Milky Way, through the redshift desert and up to the epoch of re-ionization at z>8-9. At the same time, the high spectral resolution mode will allow astronomers to study chemical abundances of stars in our Galaxy, in particular in the highly obscured regions of the Bulge, and provide the necessary follow-up of the Gaia mission. Such characteristics and versatility make MOONS the long-awaited workhorse near-IR MOS for the VLT, which will perfectly complement optical spectroscopy performed by FLAMES and VIMOS.

[Abridged] It is widely accepted that observations at mid-infrared (mid-IR) wavelengths enable the selection of galaxies with nuclear activity, which may not be revealed even in the deepest X-ray surveys. In this work new near- and mid-IR color diagnostics are explored, aiming for improved efficiency - better completeness and less contamination - in selecting AGN out to very high redshifts. We restrict our study to the James Webb Space Telescope wavelength range (0.6-27um). The criteria are created based on the predictions by state-of-the-art galaxy and AGN templates covering a wide variety of galaxy properties, and tested against control samples with deep multi-wavelength coverage (ranging from the X-rays to radio frequencies). We show that the colors Ks-[4.5], [4.5]-[8.0], and [8.0]-[24] are ideal as AGN/non-AGN diagnostics at, respectively, z<~1, 1<~z<~2.5, and z>~2.5-3. However, when the source redshift is unknown, these colors should be combined. We thus develop an improved IR criterion (using Ks and IRAC bands, KI) as a new alternative at z<~2.5. KI does not show improved completeness (50-60% overall) in comparison to commonly used IRAC-based AGN criteria, but is less affected by non-AGN contamination (revealing a >50-90% level of successful AGN selection). We also propose KIM (using Ks, IRAC, and MIPS-24um bands, KIM), which aims to select AGN hosts from local distances to as far back as the end of reionization (0<z<~7) with reduced non-AGN contamination. However, the necessary testing-constraints and the small control-sample sizes prevent the confirmation of its improved efficiency at z<~2.5. Overall, KIM shows a ~30-40% completeness and a >70-90% level of successful AGN selection. KI and KIM are built to be reliable against a ~10-20% error in flux, are based on existing filters, and are suitable for immediate use.

Multifrequency radio continuum observations (1.4-22 GHz) of a sample of reddened QSOs are presented. We find a high incidence (13/16) of radio spectral properties, such as low frequency turnovers, high frequency spectral breaks or steep power-law slopes, similar to those observed in powerful compact steep spectrum (CSS) and gigahertz-peaked spectrum (GPS) sources. The radio data are consistent with relatively young radio jets with synchotron ages <1e6-1e7yr. This calculation is limited by the lack of high resolution (milli-arcsec) radio observations. For the one source in the sample that such data are available a much younger radio age is determined, <2e3yr, similar to those of GPS/CSS sources. These findings are consistent with claims that reddened QSOs are young systems captured at the first stages of the growth of their supermassive black holes. It also suggests that expanding radio lobes may be an important feedback mode at the early stages of the evolution of AGN.

We use SDSS spectra and optical to far-infrared photometry for a sample of 31 FeLoBAL QSOs to study the relationship between the AGN-driven outflows, and obscured star formation in their host galaxies. We find that FeLoBAL QSOs invariably have IR luminosities exceeding 10^12 Solar luminosities. The AGN supplies 75% of the total IR emission, on average, but with a range from 20% to 100%. We find a clear anticorrelation between the strength of the AGN-driven outflows and the contribution from star formation to the total IR luminosity, with a much higher chance of seeing a starburst contribution in excess of 25% in systems with weak outflows than in systems with strong outflows. Moreover, we find no evidence that this effect is driven by the IR luminosity of the AGN. We conclude that radiatively driven outflows from AGN act to curtail obscured star formation in the host galaxies of reddened QSOs to less than ~25% of the total IR luminosity. This is the most direct evidence yet obtained for `quasar mode' AGN feedback.

We present optical to far-infrared photometry of 31 reddened QSOs that show evidence for radiatively driven outflows originating from AGN in their rest-frame UV spectra. We use these data to study the relationships between the AGN-driven outflows, and the AGN and starburst infrared luminosities. We find that FeLoBAL QSOs are invariably IR-luminous, with IR luminosities exceeding 10^{12} Solar luminosities in all cases. The AGN supplies 76% of the total IR emission, on average, but with a range from 20% to 100%. We find no evidence that the absolute luminosity of obscured star formation is affected by the AGN-driven outflows. Conversely, we find an anticorrelation between the strength of AGN-driven outflows, as measured from the range of outflow velocities over which absorption exceeds a minimal threshold, and the contribution from star formation to the total IR luminosity, with a much higher chance of seeing a starburst contribution in excess of 25% in systems with weak outflows than in systems with strong outflows. Moreover, we find no convincing evidence that this effect is driven by the IR luminosity of the AGN. We conclude that radiatively driven outflows from AGN can have a dramatic, negative impact on luminous star formation in their host galaxies. We find that such outflows act to curtail star formation such that star formation contributes less than ~25% of the total IR luminosity. We also propose that the degree to which termination of star formation takes place is not deducible from the IR luminosity of the AGN.

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.

At very low frequencies, the new pan-European radio telescope LOFAR is opening the last unexplored window of the electromagnetic spectrum for astrophysical studies. The revolutionary APERTIF phased arrays that are about to be installed on the Westerbork radio telescope (WSRT) will dramatically increase the survey speed for the WSRT. Combined surveys with these two facilities will deeply chart the northern sky over almost two decades in radio frequency from \sim 15 up to 1400 MHz. Here we briefly describe some of the capabilities of these new facilities and what radio surveys are planned to study fundamental issues related the formation and evolution of galaxies and clusters of galaxies. In the second part we briefly review some recent observational results directly showing that diffuse radio emission in clusters traces shocks due to cluster mergers. As these diffuse radio sources are relatively bright at low frequencies, LOFAR should be able to detect thousands of such sources up to the epoch of cluster formation. This will allow addressing many question about the origin and evolution of shocks and magnetic fields in clusters. At the end we briefly review some of the first and very preliminary LOFAR results on clusters.

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.

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.

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.

We present a multi-wavelength analysis of the properties of Extremely Red Galaxy (ERG) populations, selected in the GOODS-South/Chandra Deep Field South field. By using all the photometric and spectroscopic information available on large deep samples of EROs (645 sources), IEROs (294 sources), and DRGs (350 sources), we derive redshift distributions, identify AGN powered and Star-formation powered galaxies, and, using the radio observations of this field, estimate robust (AGN- and dust-unbiased) Star Formation Rate Densities (SFRD) for these populations. We also investigate the properties of "pure" (galaxies that conform to only one of the three ERG criteria considered) and "combined" (galaxies that verify all three criteria) sub-populations. Overall, a large number of AGN are identified (up to ~30%, based on X-rays, and mid-infrared criteria), the majority of which are type-2 (obscured) objects. Among ERGs with no evidence for AGN activity, we identify sub-populations covering a wide range of average star-formation rates, from below 10 Mo/yr to as high as 200 Mo/yr. Applying a redshift separation (1<z<2 and 2<z<3) we find significant evolution (an increase of a factor of 2 or higher) of SFRD for EROs and DRGs, while none is observed for IEROs. The former populations can contribute more than 20\% to the global SFRD at 2<z<3. The emission from AGN activity is typically not strong in the ERG population, with AGN increasing the average radio luminosity of ERG sub-populations by, nominally, less than 20%. AGN are common, however, and, if no discrimination is attempted, this could significantly increase the SFRD estimate (by over 100% in some cases). Thus, and while the contribution of star forming processes to the radio luminosity in galaxies with AGN remains uncertain, a comprehensive identification of AGN in these populations is necessary to obtain meaningful results.

We present mid-infrared spectra of six FeLoBAL QSOs at 1<z<1.8, taken with the Spitzer space telescope. The spectra span a range of shapes, from hot dust dominated AGN with silicate emission at 9.7 microns, to moderately obscured starbursts with strong Polycyclic Aromatic Hydrocarbon (PAH) emission. The spectrum of one object, SDSS 1214-0001, shows the most prominent PAHs yet seen in any QSO at any redshift, implying that the starburst dominates the mid-IR emission with an associated star formation rate of order 2700 solar masses per year. With the caveats that our sample is small and not robustly selected, we combine our mid-IR spectral diagnostics with previous observations to propose that FeLoBAL QSOs are at least largely comprised of systems in which (a) a merger driven starburst is ending, (b) a luminous AGN is in the last stages of burning through its surrounding dust, and (c) which we may be viewing over a restricted line of sight range.

We present the results of a study of a sample of 375 Extremely Red Galaxies (ERGs) in the Phoenix Deep Survey, 273 of which constitute a subsample which is 80% complete to K_s = 18.5 over an area of 1160 arcmin^2. The angular correlation function for ERGs is estimated, and the association of ERGs with faint radio sources explored. We find tentative evidence that ERGs and faint radio sources are associated at z > 0.5. A new overdensity-mapping algorithm has been used to characterize the ERG distribution, and identify a number of cluster candidates, including a likely cluster containing ERGs at 0.5 < z < 1. Our algorithm is also used in an attempt to probe the environments in which faint radio sources and ERGs are associated. We find limited evidence that the I - K_s > 4 criterion is more efficient than R - K_s > 5 at selecting dusty star-forming galaxies, rather than passively evolving ERGs.

We present mid/far-infrared photometry of nine FeLoBAL QSOs, taken using the Spitzer space telescope. All nine objects are extremely bright in the infrared, with rest-frame 1-1000 micron luminosities comparable to those of Ultraluminous Infrared Galaxies. Furthermore, a significant fraction of the infrared emission from many, and possibly all of the sample is likely to arise from star formation, with star formation rates of order several hundred solar masses per year. We combine these results with previous work to propose that FeLoBALs mark galaxies and QSOs in which an extremely luminous starburst is approaching its end, and in which a rapidly accreting supermassive black hole is in the last stages of casting off its dust cocoon. FeLoBAL signatures in high redshift QSOs and galaxies may thus be an efficient way of selecting sources at a critical point in their evolution.

We explore the nature of the radio emission of X-ray selected AGN by combining deep radio (1.4GHz; 60micro-Jy) and X-ray data with multiwavelength (optical, mid-infrared) observations in the Extended Chandra Deep Field South (E-CDFS). The fraction of radio detected X-ray sources increases from 9% in the E-CDFS to 14% in the central region of this field, which has deeper X-ray coverage from the 1Ms CDFS. We find evidence that the radio emission of up to 60% of the hard X-ray/radio matched AGN is likely associated with star-formation in the host galaxy. Firstly, the mid-IR (24micron) properties of these sources are consistent with the infrared/radio correlation of starbursts. Secondly, most of them are found in galaxies with blue rest-frame optical colours (U-V), suggesting a young stellar population. On the contrary, X-ray/radio matched AGN which are not detected in the mid-infrared have red U-V colours suggesting their radio emission is associated with AGN activity. We also find no evidence for a population of heavily obscured radio-selected AGN that are not detected in X-rays. Finally, we do no confirm previous claims for a correlation between radio emission and X-ray obscuration. Assuming that the radio continuum measures star-formation, this finding is against models where the dust and gas clouds associated with circumnuclear starbursts are spherically blocking our view to the central engine.

We present the first results from the Australia Telescope Large Area Survey (ATLAS), which consist of deep radio observations of a 3.7 square degree field surrounding the Chandra Deep Field South, largely coincident with the infrared Spitzer Wide-Area Extragalactic (SWIRE) Survey. We also list cross-identifications to infrared and optical photometry data from SWIRE, and ground-based optical spectroscopy. A total of 784 radio components are identified, corresponding to 726 distinct radio sources, nearly all of which are identified with SWIRE sources. Of the radio sources with measured redshifts, most lie in the redshift range 0.5-2, and include both star-forming galaxies and active galactic nuclei (AGN). We identify a rare population of infrared-faint radio sources which are bright at radio wavelengths but are not seen in the available optical, infrared, or X-ray data. Such rare classes of sources can only be discovered in wide, deep surveys such as this.

We estimate the star-formation rates and the stellar masses of the Extremely Red objects (EROs) detected in a 180arcmin2 Ks-band survey (Ks~20mag). This sample is complemented by sensitive 1.4GHz radio observations (12micro-Jy; 1sigma rms) and multiwaveband photometric data (UBVRIJ) as part of the Phoenix Deep Survey. For bright K<19.5mag EROs in this sample (I-K>4mag; total of 177) we use photometric methods to discriminate dust-enshrouded active systems from early-type galaxies and to constrain their redshifts. Radio stacking is then employed to estimate mean radio flux densities of 8.6 (3sigma) and 6.4micro-Jy (2.4sigma) for the dusty and early-type subsamples respectively. Assuming that dust enshrouded active EROs are powered by star-formation the above radio flux density at the median redshift of z=1 translates to a radio luminosity of 4.5e22W/Hz and a star-formation rate of SFR=25Mo/yr. Combining this result with photometric redshift estimates we find a lower limit to the star-formation rate density of ~0.02Mo/yr/Mpc^3 for the K<19.5mag dusty EROs in the range z=0.85-1.35. Comparison with the SFR density estimated from previous ERO samples (with similar selection criteria) using optical emission lines, suffering dust attenuation, suggests a mean dust reddening of at least E(B-V)~0.5 for this population. We further use the Ks-band luminosity as proxy to stellar mass and argue that the dust enshrouded EROs in our sample are massive systems, M>5e10Mo. We also find that EROs represent a sizable fraction (~50%) of the number density of galaxies more massive than M=5e10Mo at z~1, with almost equal contributions from dusty and early types. Similarly, we find that EROs contribute about half of the mass density of the Universe at z~1 after taking into account incompleteness because of the limit K=19.5mag.

We present optical and X-ray identifications for the sixty-four radio sources in the GOODS-S ACS field revealed in the ATCA 1.4 GHz survey of the Chandra Deep Field South. Optical identifications are made using the ACS images and catalogs, while the X-ray view is provided by the Chandra X-ray Observatory 1 Ms observations. Redshifts for the identified sources are drawn from publicly available catalogs of spectroscopic observations and multi-band photometric-based estimates. Using this multiwavelength information we provide a first characterization of the faint radio source population in this region. The sample contains a mixture of star-forming galaxies and active galactic nuclei, as identified by their X-ray properties and optical spectroscopy. A large number of morphologically disturbed galaxies is found, possibly related to the star-formation phenomena. In spite of the very deep optical data available in this field, seven of the sixty-four radio sources have no optical identification to z(850)~28 mag. Only one of these is identified in the X-rays.

We present deep submillimetre photometry for 14 galaxies at z=0.5 that are hosts of type 1a supernovae, with the aim of examining the evolution of dust mass and extinction in normal galaxies. We combine these results with our previous observations of 17 z~0.5 SN1a hosts to look for any evolution in the dust content of normal galaxies between z=0 and z=0.5. The average observed frame 850micron flux of SN1a hosts in the full sample, excluding 2 bright individually detected objects, is 0.44 +/- 0.22 mJy. This flux level is consistent with there being little or no evolution in the dust content, or optical extinction, of normal galaxies from z=0 to z=0.5. One galaxy, the host of SN1996cf, is detected individually, and we also present a deep HST STIS image for this object. It appears to be an edge on disk system, similar to the submm bright host of SN1997ey. We thus examine the dust properties of these and one other individually detected object. 450-to-850 micron flux ratios and limits suggest that the dust in the two brightest submm sources, SN1996cf and SN1997ey, is cold, T ~20K, implying that they contain a substantial mass of dust ~10^9 Msun. The presence of two bright (F850 >7mJy) submm sources at z~0.5 in a sample of ostensibly normal galaxies is surprising, and has important implications. It supports the idea that a substantial part of the Cosmic Infrared Background (CIB) may be produced at z<1, while also suggesting that 'foreground' objects such as these may be a significant 'contaminant' in submm surveys. Finally, we examine the overall submm luminosity distribution at z=0.5 implied by our results, and conclude that either there is substantial evolution in the submm luminosity function from z=0 to 0.5, or our submm detected sources are somehow not representative of the bulk of galaxies at this redshift.

The Phoenix Deep Survey is a multi-wavelength survey based on deep 1.4 GHz radio imaging, reaching well into the sub-100 microJy level. One of the aims of this survey is to characterize the sub-mJy radio population, exploring its nature and evolution. In this paper we present the catalog and results of the spectroscopic observations aimed at characterizing the optically ``bright'' (R<~ 21.5 mag) counterparts of faint radio sources. Out of 371 sources with redshift determination, 21% have absorption lines only, 11% show AGN signatures, 32% are star-forming galaxies, 34% show narrow emission lines that do not allow detailed spectral classification (due to poor signal-to-noise ratio and/or lack of diagnostic emission lines) and the remaining 2% are identified with stars. For the star-forming galaxies with a Balmer decrement measurement we find a median extinction of A(Ha)=1.9 mag, higher than that of optically selected samples. This is a result of the radio selection, which is not biased against dusty systems. Using the available spectroscopic information, we estimate the radio luminosity function of star-forming galaxies in two independent redshift bins at z~0.1 and 0.3 respectively. We find direct evidence for strong luminosity evolution of these systems consistent with L(1.4 GHz) ~ (1+z)^(2.7).

In this paper we explore the clustering properties and the environment of the Extremely Red Objects (EROs; I-K>4mag) detected in a ~180arcmin^2 deep (Ks~20mag) Ks-band survey of a region within the Phoenix Deep Survey, an on-going multiwavelength program aiming to investigate the nature and the evolution of faint radio sources. Using our complete sample of 289 EROs brighter than Ks=20mag we estimate a statistically significant (~3.7sigma) angular correlation function signal with amplitude Aw=8.7^{+2.1}_{-1.7}x10^{-3} consistent with earlier work based on smaller samples. This amplitude suggests a clustering length in the range ro=12-17h^{-1}Mpc, implying that EROs trace regions of enhanced density. Using a novel method we further explore the association of EROs with galaxy overdensities by smoothing the K-band galaxy distribution using the matched filter algorithm of Postman et al. (1996) and then cross-correlating the resulting density maps with the ERO positions. Our analysis provides direct evidence that EROs are associated with overdensities at redshifts z>1. We also exploit the deep radio 1.4GHz data (limiting flux 60microJy) available to explore the association of EROs and faint radio sources and whether the two populations trace similar large scale structures. Cross-correlation of the two samples (after excluding 17EROs with radio counterparts) gives a 2sigma signal only for the sub-sample of high-z radio sources (z>0.6). Although the statistics are poor this suggests that it is the high-z radio sub-sample that traces similar structures with EROs.

The Phoenix Deep Survey is a multi-wavelength galaxy survey based on deep 1.4 GHz radio imaging (Hopkins et al., 2003). The primary goal of this survey is to investigate the properties of star formation in galaxies and to trace the evolution in those properties to a redshift z=1, covering a significant fraction of the age of the Universe. By compiling a sample of star-forming galaxies based on selection at radio wavelengths we eliminate possible biases due to dust obscuration, a significant issue when selecting objects at optical and ultraviolet wavelengths. In this paper, we present the catalogs and results of deep optical (UBVRI) and near-infrared (Ks) imaging of the deepest region of the existing decimetric radio imaging. The observations and data-processing are summarised and the construction of the optical source catalogs described, together with the details of the identification of candidate optical counterparts to the radio catalogs. Based on our UBVRIKs imaging, photometric redshift estimates for the optical counterparts to the radio detections are explored.

The radio properties of hard (2-8keV) X-ray selected sources are explored by combining a single 50ks XMM-Newton pointing with the ultra-deep and homogeneous Phoenix radio (1.4GHz) survey (Hopkins et al. 2003). A total of 43 sources are detected above the X-ray flux limit f_X(2-8keV)=7.7e-15cgs with 14 of them exhibiting radio emission above ~40muJy (3sigma). The X-ray/radio matched population lies in the borderline between radio loud and quiet AGNs and comprises sources with both soft and hard X-ray spectral properties suggesting both obscured and unobscured systems. The spectroscopically identified sub-sample (total of 6 X-ray/radio matches) comprises narrow emission line AGNs (4) with hard X-ray spectral properties and broad line sources (2) with soft X-ray spectra. We find evidence that the fraction of X-ray/radio matches increases from ~20% for sources with rest-frame column density N_H<1e22cm^-2 to \~50% for more absorbed systems. Poor statistics however, limit the significance of the above result to the ~2sigma level. Also, the X-ray/radio matched sources have flatter coadded X-ray spectrum (Gamma=1.78^{+0.05}_{-0.03}) compared to sources without radio emission (Gamma=2.00^{+0.03}_{-0.04}). A possible explanation for the higher fraction of absorbed sources with radio emission at the muJy level is the presence of circum-nuclear starburst activity that both feeds and obscures the central engine. For a small sub-sample of z~0.4 radio emitting AGNs with N_H>1e22cm^-2 their combined spectrum exhibits a soft X-ray component that may be associated with star-formation activity, although other possibilities cannot be excluded. We also find that radio emitting AGNs make up about 13-20% of the hard-band X-ray background depending on the adopted normalisation.

Insensitive to dust obscuration, radio wavelengths are ideal to study star-forming galaxies free of dust induced biases. Using data from the Phoenix Deep Survey, we have identified a sample of star-forming extremely red objects (EROs). Stacking of the radio images of the radio-undetected star-forming EROs revealed a significant radio detection. Using the expected median redshift, we estimate an average star-formation rate of 61 M_sun/yr for these galaxies.

We present deep sub-millimetre observations of seventeen galaxies at z=0.5, selected through being hosts of a type 1a supernova. Two galaxies are detected directly, and the sample is detected statistically with a mean 850 micron flux of 1.01mJy +/- 0.33mJy, which is 25% - 135% higher than locally. We infer that the mean value of A_v in normal galaxies at z=0.5 is comparable to or greater than the mean A_v in local normal galaxies, in agreement with galaxy chemical evolution models and indirect observational evidence. Scaling from the local value given by Rowan-Robinson (2003) gives a mean extinction at z=0.5 of A_v = 0.56 +/- 0.17. The dust in the brightest sub-mm object in our sample is best interpreted as normal `cirrus' dust similar to that seen locally. The detection rate of our sample suggests that some sources found in blank-field sub-mm surveys may not be high redshift starbursts, but rather cirrus galaxies at moderate redshifts and with lower star formation rates. Finally, an increase in host dust extinction with redshift may impact the cosmological results from distant supernova searches. This emphasizes the need to carefully monitor dust extinction when using type Ia supernovae to measure the cosmological parameters.

The Phoenix Deep Survey (PDS) is a multiwavelength survey based on deep 1.4 GHz radio observations used to identify a large sample of star forming galaxies to z=1. Photometric redshifts are estimated for the optical counterparts to the radio-detected galaxies, and their uncertainties quantified by comparison with spectroscopic redshift measurements. The photometric redshift estimates and associated best-fitting spectral energy distributions are used in a stacking analysis exploring the mean radio properties of U-band selected galaxies. Average flux densities of a few microJy are measured.

The Phoenix Deep Survey (PDS) is a multiwavelength survey based on deep 1.4 GHz radio observations used to identify a large sample of star forming galaxies to z=1. Here we present an exploration of the evolutionary constraints on the star-forming population imposed by the 1.4 GHz source counts, followed by an analysis of the average properties of extremely red galaxies in the PDS, by using the "stacking" technique.

In this paper we use a 50ks XMM-Newton pointing overlapping with the Phoenix Deep Survey, a homogeneous radio survey reaching muJy sensitivities, to explore the X-ray properties and the evolution of star-forming galaxies. UV, optical and NIR photometry is available and is used to estimate photometric redshifts and spectral types for radio sources brighter than R=21.5mag (total of 82). Sources with R<21.5mag and spiral galaxy SEDs (34) are grouped into two redshift bins with a median of z=0.240 and 0.455 respectively. Stacking analysis for both the 0.5-2 and 2-8keV bands is performed on these subsamples. A high confidence level signal (>3.5sigma) is detected in the 0.5-2keV band corresponding to a mean flux of ~3e-16cgs for both subsamples. This flux translates to mean luminosities of ~5e40 and 1.5e41cgs for the z=0.240 and 0.455 subsamples respectively. Only a marginally significant signal (2.6sigma) is detected in the 2-8keV band for the z=0.455 subsample. We argue that the stacked signal above is dominated by star-formation. The mean L_X/L_B ratio and the mean L_X of the two subsamples are found to be higher than optically selected spirals and similar to starbursts. We also find that the mean L_X and L_1.4 of the faint radio sources studied here are consistent with the L_X-L_1.4 correlation of local star-forming galaxies. Moreover, the X-ray emissivity of sub-mJy sources to z~0.3 is found to be elevated compared to local HII galaxies. The observed increase is consistent with L_X evolution of the form (1+z)^3. Assuming that our sample is indeed dominated by starbursts this is direct evidence for evolution of such systems at X-ray wavelengths. Using an empirical L_X to SFR conversion we estimate a global SFR density at z~0.3 of \~0.029M_o/yr/Mpc in agreement with previous studies.

Many investigations of star formation rates (SFRs) in galaxies have explored details of dust obscuration, with a number of recent analyses suggesting that obscuration appears to increase in systems with high rates of star formation. To date these analyses have been primarily based on nearby (z < 0.03) or UV selected samples. Using 1.4 GHz imaging and optical spectroscopic data from the Phoenix Deep Survey, the SFR-dependent obscuration is explored. The use of a radio selected sample shows that previous studies exploring SFR-dependent obscurations have been biased against obscured galaxies. The observed relation between obscuration and SFR is found to be unsuitable to be used as an obscuration measure for individual galaxies. Nevertheless, it is shown to be successful as a first order correction for large samples of galaxies where no other measure of obscuration is available, out to intermediate redshifts (z ~ 0.8).

(Abridged) We examine the power source of 41 local Ultraluminous Infrared Galaxies using archival infrared and optical photometry. We fit the observed Spectral Energy Distributions (SEDs) with starburst and AGN components; each component being drawn from a family of templates. We find all of the sample require a starburst, whereas only half require an AGN. In 90% of the sample the starburst provides over half the IR emission, with a mean fractional luminosity of 82%. When combined with other galaxy samples we find that starburst and AGN luminosities correlate over 6 decades in IR luminosity suggesting that a common factor governs both luminosities, plausibly the gas masses in the nuclear regions. We find that the mid-IR 7.7 micron line-continuum ratio is no indication of the starburst luminosity, or the fractional AGN luminosity, and therefore that this ratio is not a reliable diagnostic of the power source in ULIRGs. We propose that the scatter in the radio-IR correlation in ULIRGs is due to a skewed starburst IMF and/or relic relativistic electrons from a previous starburst, rather than contamination from an obscured AGN. We show that most ULIRGs undergo multiple starbursts during their lifetime, and by inference that mergers between more than two galaxies may be common amongst ULIRGs. Our results support the evolutionary model for ULIRGs proposed by Farrah et al 2001, where they can follow many different evolutionary paths of starburst and AGN activity in transforming merging spiral galaxies into elliptical galaxies, but that most do not go through an optical QSO phase. The lower level of AGN activity in our local sample than in z~1 HLIRGs implies that the two samples are distinct populations. We postulate that different galaxy formation processes at high-z are responsible for this difference.

The initial Phoenix Deep Survey (PDS) observations with the Australia Telescope Compact Array have been supplemented by additional 1.4 GHz observations over the past few years. Here we present details of the construction of a new mosaic image covering an area of 4.56 square degrees, an investigation of the reliability of the source measurements, and the 1.4 GHz source counts for the compiled radio catalogue. The mosaic achieves a 1-sigma rms noise of 12 microJy at its most sensitive, and a homogeneous radio-selected catalogue of over 2000 sources reaching flux densities as faint as 60 microJy has been compiled. The source parameter measurements are found to be consistent with the expected uncertainties from the image noise levels and the Gaussian source fitting procedure. A radio-selected sample avoids the complications of obscuration associated with optically-selected samples, and by utilising complementary PDS observations including multicolour optical, near-infrared and spectroscopic data, this radio catalogue will be used in a detailed investigation of the evolution in star-formation spanning the redshift range 0 < z < 1. The homogeneity of the catalogue ensures a consistent picture of galaxy evolution can be developed over the full cosmologically significant redshift range of interest. The 1.4 GHz mosaic image and the source catalogue are available on the web at http://www.atnf.csiro.au/~ahopkins/phoenix/ or from the authors by request.

We report the discovery of the optical and near-infrared counterparts to GRB 001011. The GRB 001011 error box determined by Beppo-SAX was simultaneously imaged in the near-infrared by the 3.58-m New Technology Telescope and in the optical by the 1.54-m Danish Telescope ~8 hr after the gamma-ray event. Here we implement the colour-colour discrimination technique proposed by Rhoads (2001) and extend it using near-IR data as well. We present the results provided by an automatic colour-colour discrimination pipe-line developed to discern the different populations of objects present in the GRB 001011 error box. Our software revealed three candidates based on single-epoch images. Second-epoch observations carried out ~3.2 days after the burst revealed that the most likely candidate had faded, thus identifying it with the counterpart to the GRB. In deep R-band images obtained 7 months after the burst a faint (R=25.38+/-0.25) elongated object, presumably the host galaxy of GRB 001011, was detected at the position of the afterglow. The GRB 001011 afterglow is the first discovered with the assistance of colour-colour diagram techniques. We discuss the advantages of using this method and its application to error boxes determined by future missions.

In the course of the follow-up multiwavelength study of a deep radio survey we have discovered that the milli-Jansky radio source PDFJ011423 is a low-redshift (z = 0.65) extremely red galaxy (ERG) with K=15.3, R-K = 5.8 and J-K = 3.1. Optical, infrared and radio photometry, together with optical and near-infrared spectroscopy, reveal a heavily obscured galaxy (A_V=5-6, from the observed Balmer decrement) undergoing vigorous star formation and presenting an active galactic nucleus (AGN). PDFJ011423 is a representative member of the dusty ERG population, providing a local counterpart for studying more distant ERGs.

The Infrared Space Observatory (ISO) is used to carry out mid-IR (7 and 15 micron) and far-IR (90 micron) observations of a sample of star-forming sub-mJy radio sources. By selecting the sample at radio wavelengths, one avoids biases due to dust obscuration. It is found that the mid-IR luminosities, covering the PAH features, measure the star formation rate for galaxies with $P_{1.4 GHz} < 10^{23}$ W Hz$^{-1}$. This is further confirmed using the H$\alpha$ luminosities. The far-IR emission is also found to trace the SFR over the whole range of radio and H$\alpha$ luminosities. The implication of the mid-IR measurements in estimating the SFRs from the future infrared space missions (SIRTF and ASTRO-F) is discussed.

The cluster DC0107-46 (Abell 2877) lies within the Phoenix Deep Survey, made at 1.4 GHz with the Australia Telescope Compact Array. Of 89 known optical cluster members, 70 lie within the radio survey area. Of these 70 galaxies, 15 (21%) are detected, with luminosities as faint as 10^20 W/Hz. Spectroscopic observations are available for 14/15 of the radio-detected cluster galaxies. Six galaxies show only absorption features and are typical low-luminosity AGN radio sources. One galaxy hosts a Seyfert 2 nucleus, two are star-forming galaxies, and the remaining five may be star-forming galaxies, AGNs, or both.

The redshift dependence of the luminosity density in certain wavebands (e.g. UV and H-alpha) can be used to infer the history of star formation in the populations of galaxies producing this luminosity. This history is a useful datum in studies of galaxy evolution. It is therefore important to understand the errors that attend the inference of star formation rate densities from luminosity densities. This paper explores the self-consistency of star formation rate diagnostics by reproducing commonly used observational procedures in a model with known galaxy populations, evolutionary histories and spectral emission properties. The study reveals a number of potential sources of error in the diagnostic processes arising from the differential evolution of different galaxy types. We argue that multi-wavelength observations can help to reduce these errors.

The study of the nature of faint radio sources is of great importance since a significant fraction of these objects is thought to be composed of actively star-forming galaxies. Due to the increased sensitivity of radio telescopes, we are now not only able to catalogue large numbers of these sources in the sub-millijansky regime, but also to start the study of the nature of increasingly fainter microjansky sources. This paper presents a new very deep 1.4 GHz radio survey made as a part of the Phoenix Deep Survey, a project aimed to study the nature of the faintest radio sources. With a limiting sensitivity of 45 microJy, this new survey has allowed us to assemble a large number of sources with 1.4 GHz flux densities below 100 microJy. The resulting source counts and the analysis of the optical properties of the faintest radio sources are presented.

We present a deep 1.4 GHz survey made with the Australia Telescope Compact Array (ATCA), having a background RMS of 9 microJy near the image phase centre, up to 25 microJy at the edge of a 50' field of view. Over 770 radio sources brighter than 45 microJy have been catalogued in the field. The differential source counts in the deep field provide tentative support for the growing evidence that the microjansky radio population exhibits significantly higher clustering than found at higher flux density cutoffs. The optical identification rate on CCD images is approximately 50% to R=22.5, and the optical counterparts of the faintest radio sources appear to be mainly single galaxies close to this optical magnitude limit.