Aims: We investigate the properties of the polarised radio population in the central 6.5 deg$^{2}$ of the XXL-South field observed at 2.1 GHz using the Australia Telescope Compact Array (ATCA) in 81 pointings with a synthesised beam of FWHM 5.2''. We also investigate the ATCA's susceptibility to polarisation leakage. Methods: We performed a survey of a 5.6 deg$^{2}$ subregion and calculated the number density of polarised sources. We derived the total and polarised spectral indices, in addition to comparing our source positions with those of X-ray-detected clusters. We measured the polarisation of sources in multiple pointings to examine leakage in the ATCA. Results: We find 39 polarised sources, involving 50 polarised source components, above a polarised flux density limit of 0.2 mJy at 1.332 GHz. The number density of polarised source components is comparable with recent surveys, although there is an indication of an excess at $\sim1$ mJy. We find that those sources coincident with X-ray clusters are consistent in their properties with regard to the general population. In terms of the ATCA leakage response, we find that ATCA mosaics with beam separation of $\lesssim 2/3$ of the primary beam FWHM have off-axis linear polarisation leakage $\lesssim 1.4$ % at 1.332 GHz.

We present the first systematic study of the radio continuum size evolution of star-forming galaxies (SFGs) over the redshift range $0.35<z<2.25$. We use the VLA COSMOS 3GHz map (noise $\rm rms=2.3\,\mu Jy \,beam^{-1}$, $\theta_{\rm beam}=0.75\,\rm arcsec$) to construct a mass-complete sample of 3184 radio-selected SFGs that reside on and above the main-sequence (MS) of SFGs. We find no clear dependence between the radio size and stellar mass, $M_{\star}$, of SFGs with $10.5\lesssim\log(M_\star/\rm M_\odot)\lesssim11.5$. Our analysis suggests that MS galaxies are preferentially extended, while SFGs above the MS are always compact. The median effective radius of SFGs on (above) the MS of $R_{\rm eff}=1.5\pm0.2$ ($1.0\pm0.2$) kpc remains nearly constant with cosmic time; a parametrization of the form $R_{\rm eff}\propto(1+z)^\alpha$ yields a shallow slope of only $\alpha=-0.26\pm0.08\,(0.12\pm0.14)$ for SFGs on (above) the MS. The size of the stellar component of galaxies is larger than the extent of the radio continuum emission by a factor $\sim$2 (1.3) at $z=0.5\,(2)$, indicating star formation is enhanced at small radii. The galactic-averaged star formation rate surface density $(\Sigma_{\rm SFR})$ scales with the distance to the MS, except for a fraction of MS galaxies ($\lesssim10\%$) that harbor starburst-like $\Sigma_{\rm SFR}$. These "hidden" starbursts might have experienced a compaction phase due to disk instability and/or merger-driven burst of star formation, which may or may not significantly offset a galaxy from the MS. We thus propose to jointly use $\Sigma_{\rm SFR}$ and distance to the MS to better identify the galaxy population undergoing a starbursting phase.

In this data paper we present and characterise the multi-component radio sources identified in the VLA-COSMOS Large Project at 3 GHz (0.75 arcsec resolution, 2.3 {\mu}Jy/beam rms), i.e. the radio sources which are composed of two or more radio blobs.The classification of objects into multi-components was done by visual inspection of 351 of the brightest and most extended blobs from a sample of 10,899 blobs identified by the automatic code blobcat. For that purpose we used multi-wavelength information of the field, such as the 1.4-GHz VLA-COSMOS data and the UltraVISTA stacked mosaic available for COSMOS. We have identified 67 multi-component radio sources at 3 GHz: 58 sources with AGN powered radio emission and 9 star-forming galaxies. We report 8 new detections that were not observed by the VLA-COSMOS Large Project at 1.4 GHz, due to the slightly larger area coverage at 3 GHz. The increased spatial resolution of 0.75 arcsec has allowed us to resolve (and isolate) multiple emission peaks of 28 extended radio sources not identified in the 1.4-GHz VLA-COSMOS map. We report the multi-frequency flux densities (324 MHz, 325 MHz, 1.4 GHz & 3 GHz), star-formation-rates, and stellar masses of these objects. Multi-component objects at 3-GHz VLA-COSMOS inhabit mainly massive galaxies (>10^10.5 Msun). The majority of the multi-component AGN lie below the main-sequence of star-forming galaxies (SFGs), in the green valley and the quiescent region. We provide detailed description of the objects: amongst the AGN there are 2 head-tail, 10 core-lobe, 9 wide-angle-tail (WAT), 8 double-double or Z-/X-shaped, 3 bent-tail radio sources, and 26 symmetric sources, while amongst the SFGs we find the only star-forming ring seen in radio emission in COSMOS. We report a large number (32/58) of disturbed/bent multi-component AGN, 18 of which do not lie within X-ray groups in COSMOS (0.08 < z < 1.53). [abridged]

We construct the average radio spectral energy distribution (SED) of highly star-forming galaxies (HSFGs) up to z~4. Infrared and radio luminosities are bound by a tight correlation that is defined by the so-called q parameter. This infrared-radio correlation provides the basis for the use of radio luminosity as a star-formation tracer. Recent stacking and survival analysis studies find q to be decreasing with increasing redshift. It was pointed out that a possible cause of the redshift trend could be the computation of rest-frame radio luminosity via a single power-law assumption of the star-forming galaxies' (SFGs) SED.To test this, we constrained the shape of the radio SED of a sample of HSFGs. To achieve a broad rest-frame frequency range, we combined previously published VLA observations of the COSMOS field at 1.4 GHz and 3 GHz with unpublished GMRT observations at 325 MHz and 610 MHz by employing survival analysis to account for non-detections in the GMRT maps. We selected a sample of HSFGs in a broad redshift range (0.3<z<4,SFR>100M0/yr) and constructed the average radio SED. By fitting a broken power-law, we find that the spectral index changes from $\alpha_1=0.42\pm0.06$ below a rest-frame frequency of 4.3 GHz to $\alpha_2=0.94\pm0.06$ above 4.3 GHz. Our results are in line with previous low-redshift studies of HSFGs (SFR>10M0/yr) that show the SED of HSFGs to differ from the SED found for normal SFGs (SFR<10M0/yr). The difference is mainly in a steeper spectrum around 10 GHz, which could indicate a smaller fraction of thermal free-free emission. Finally, we also discuss the impact of applying this broken power-law SED in place of a simple power-law in K-corrections of HSFGs and a typical radio SED for normal SFGs drawn from the literature. We find that the shape of the radio SED is unlikely to be the root cause of the q-z trend in SFGs.

We study a sample of 1,604 moderate-to-high radiative luminosity active galactic nuclei (HLAGN) selected at 3 GHz within the VLA-COSMOS 3 GHz Large Project. These were classified by combining multiple AGN diagnostics: X-ray data, mid-infrared data and broad-band spectral energy distribution fitting. We decompose the total radio 1.4 GHz luminosity ($\mathrm{L_{1.4\ GHz,TOT}}$) into the emission originating from star formation and AGN activity by measuring the excess in $\mathrm{L_{1.4\ GHz,TOT}}$ relative to the infrared-radio correlation of star-forming galaxies. To quantify the excess, for each source we calculate the AGN fraction ($\mathrm{f_{AGN}}$), the fractional contribution of AGN activity to $\mathrm{L_{1.4\ GHz,TOT}}$. The majority of the HLAGN, $(68.0\pm1.5)\%$, are dominated by star-forming processes ($f_{AGN}\leq0.5$), while $(32.0\pm1.5)\%$ are dominated by AGN-related radio emission ($0.5<f_{AGN}\leq1$). We use the AGN-related 1.4 GHz emission to derive the 1.4 GHz AGN luminosity functions of HLAGN. By assuming pure density and pure luminosity evolution models we constrain their cosmic evolution out to $z\sim6$, finding $\mathrm{\Phi^* (z) \propto (1+z)^{(2.64\pm0.10)+(-0.61\pm0.04) z}}$ and $\mathrm{L^* (z) \propto (1+z)^{(3.97\pm0.15) + (-0.92\pm0.06)z}}$. These evolutionary laws show that the number and luminosity density of HLAGN increased from higher redshifts ($z\sim6$) up to a maximum in the redshift range $ 1<z<2.5$, followed by a decline towards local values. By scaling the 1.4 GHz AGN luminosity to kinetic luminosity using the standard conversion, we estimate the kinetic luminosity density as a function of redshift. We compare our result to the semi-analytic models of radio mode feedback finding that this feedback could have played an important role in the context of AGN-host coevolution in HLAGN which show evidence of AGN-related radio emission ($f_{AGN}>0$).

In this paper we present the optical, near-infrared (NIR) and X-ray identifications of the 6287 radio sources detected in the 2.1 GHz deep radio survey down to a median rms of ~ 41microJy/beam obtained with the Australia Telescope Compact Array (ATCA) in the XXL-S field. The goal of this paper is to provide a multi wavelength catalogue of the counterparts of the radio sources to be used in further studies. For the optical and NIR identification of the radio sources, we used the likelihood ratio (LR) technique, slightly modified in order to take into account the presence of a large number of relatively bright counterparts close to the radio sources. This procedure led to the identification of optical/NIR counterparts for 4770 different radio sources (~77% of the whole radio sample), 414 of which also have an X-ray counterpart. This fraction of identification is in agreement with previous radio-optical association studies at a similar optical magnitude depth, but is relatively low in comparison to recent work conducted in other radio fields using deeper optical and NIR data. The analysis of optical and NIR properties of radio sources shows that, regardless of the radio flux limit of a radio survey, the nature of the identified sources is strongly dependent on the depth of the optical/NIR used in the identification process. Only with deep enough optical/NIR data will we be able to identify a significant fraction of radio sources with red (z_{DEC}-K) counterparts whose radio emission is dominated by nuclear activity rather than starburst activity.

Exploring how radio-emitting active galactic nuclei (AGN) behave and evolve with time is critical for understanding how AGN feedback impacts galaxy evolution. In this work, we investigate the relationship between 1.4 GHz radio continuum AGN luminosity ($L^{\rm AGN}_{\rm 1.4}$), specific black hole accretion rate (s-BHAR, defined as the accretion luminosity relative to the galaxy stellar mass) and redshift, for a luminosity-complete sample of radio-selected AGN in the VLA COSMOS 3 GHz Large Project. The sample was originally selected from radio-continuum observations at 3 GHz, and includes about 1800 radio AGN identified via ($>2\sigma$) radio-excess relative to the infrared-radio correlation of star-forming galaxies. We further select a subsample of over 1200 radio AGN that is complete in $L^{\rm AGN}_{\rm 1.4}$ over different redshift ranges, out to z~4, and use X-ray stacking to calculate the average s-BHAR in each $L^{\rm AGN}_{\rm 1.4}$-$z$ bin. We find that the average s-BHAR is independent of $L^{\rm AGN}_{\rm 1.4}$, at all redshifts. However, we see a strong increase of s-BHAR with redshift, at fixed $L^{\rm AGN}_{\rm 1.4}$. This trend resembles the strong increase in the fraction of star-forming host galaxies (based on the $(NUV-r)$ / $(r-J)$ colours) with redshift, at fixed $L^{\rm AGN}_{\rm 1.4}$. A possible explanation for this similarity might imply a link between average AGN radiative power and availability of cold gas supply within the host galaxy. This study corroborates the idea that radio-selected AGN become more radiatively efficient towards earlier epochs, independently of their radio power.

We show how the XXL multiwavelength survey can be used to shed light on radio galaxies and their environment. Two prominent radio galaxies were identified in a visual examination of the mosaic of XXL-North obtained with the Giant Metrewave Radio Telescope at 610MHz. Counterparts were searched for in other bands. Spectroscopic redshifts from the GAMA database were used to identify clusters and/or groups of galaxies, estimate their masses with the caustic method, and quantify anisotropies in the surrounding galaxy distribution via a Fourier analysis. Both radio galaxies are of FR I type and are hosted by early-type galaxies at a redshift of 0.138. The first radio source, named the Exemplar, has a physical extent of about 400 kpc; it is located in the cluster XLSSC112, which has a temperature of about 2 keV, a total mass of about $10^{14} M_\odot$, and resides in an XXL supercluster with eight known members. The second source, named the Double Irony, is a giant radio galaxy with a total length of about 1.1 Mpc. Its core coincides with a cataloged point-like X-ray source, but no extended X-ray emission from a surrounding galaxy cluster was detected. However, from the optical data we determined that the host is the brightest galaxy in a group that is younger, less virialized, and less massive than the Exemplar's cluster. A friends-of-friends analysis showed that the Double Irony's group is a member of the same supercluster as the Exemplar. There are indications that the jets and plumes of the Double Irony have been deflected by gas associated with the surrounding galaxy distribution. Another overdensity of galaxies (the tenth) containing a radio galaxy was found to be associated with the supercluster. Radio galaxies can be used to find galaxy clusters/groups that are below the current sensitivity of X-ray surveys.

We present a "super-deblended" far-infrared to (sub)millimeter photometric catalog in the Cosmic Evolution Survey (COSMOS), prepared with the method recently developed by Liu et al. 2018, with key adaptations. We obtain point spread function (PSF) fitting photometry at fixed prior positions including 88,008 galaxies detected in either VLA 1.4~GHz, 3~GHz and/or MIPS 24~$\mu$m images. By adding a specifically carved mass-selected sample (with an evolving stellar mass limit), a highly complete prior sample of 194,428 galaxies is achieved for deblending FIR/(sub)mm images. We performed ``active' removal of non relevant priors at FIR/(sub)mm bands using spectral energy distribution (SED) fitting and redshift information. In order to cope with the shallower COSMOS data we subtract from the maps the flux of faint non-fitted priors and explicitly account for the uncertainty of this step. The resulting photometry (including data from Spitzer, Herschel, SCUBA2, AzTEC, MAMBO and NSF's Karl G. Jansky Very Large Array at 3~GHz and 1.4~GHz) displays well behaved quasi-Gaussian uncertainties, calibrated from Monte Carlo simulations and tailored to observables (crowding, residual maps). Comparison to ALMA photometry for hundreds of sources provide a remarkable validation of the technique. We detect 11,220 galaxies over the 100--1200~$\mu$m range, extending to $z_{\rm phot}\sim7$. We conservatively selected a sample of 85 $z>4$ high redshift candidates, significantly detected in the FIR/(sub)mm, often with secure radio and/or Spitzer/IRAC counterparts. This provides a chance to investigate the first generation of vigorous starburst galaxies (SFRs$\sim1000\mathrm{M}_\odot$~yr$^{-1}$). The photometric and value added catalogs are publicly released.

The classification of the host galaxies of the radio sources in the 25 deg$^2$ ultimate XMM extragalactic survey south field (XXL-S) is presented. XXL-S was surveyed at 2.1 GHz with the Australia Telescope Compact Array (ATCA) and is thus far the largest area radio survey conducted down to rms flux densities of $\sigma \sim 41$ $\mu$Jy beam$^{-1}$. Of the 6287 radio sources in XXL-S, 4758 (75.7%) were cross-matched to an optical counterpart using the likelihood ratio technique. There are 1110 spectroscopic redshifts and 3648 photometric redshifts available for the counterparts, of which 99.4% exist out to $z \sim 4$. A number of multiwavelength diagnostics, including X-ray luminosities, mid-infrared colours, spectral energy distribution fits, radio luminosities, and optical emission lines and colours, were used to classify the sources into three types: low-excitation radio galaxies (LERGs), high-excitation radio galaxies (HERGs), and star-forming galaxies (SFGs). The final sample contains 1729 LERGs (36.3%), 1159 radio-loud HERGs (24.4%), 296 radio-quiet HERGs (6.2%), 558 SFGs (11.7%), and 1016 unclassified sources (21.4%). The LERGs tend to exist in the most massive galaxies with low star formation rates and redder colours, whereas the HERGs and SFGs exist in galaxies of lower mass, higher star formation rates, and bluer colours. The fraction of blue host galaxies is higher for radio-quiet HERGs than for radio-loud HERGs. LERGs and radio-loud HERGs are found at all radio luminosities, but radio-loud HERGs tend to be more radio luminous than LERGs at a given redshift. These results are consistent with the emerging picture in which LERGs exist in the most massive quiescent galaxies typically found in clusters with hot X-ray halos and HERGs are associated with ongoing star formation in their host galaxies via the accretion of cold gas.

We present the 25 square-degree GMRT-XXL-N 610 MHz radio continuum survey, conducted at 50~cm wavelength with the Giant Metrewave Radio Telescope (GMRT) towards the XXL Northern field (XXL-N). We combined previously published observations of the XMM-Large Scale Structure (XMM-LSS) field, located in the central part of XXL-N, with newly conducted observations towards the remaining XXL-N area, and imaged the combined data-set using the Source Peeling and Atmospheric Modeling ({\sc SPAM}) pipeline. The final mosaic encompasses a total area of $30.4$ square degrees, with ${\rm rms}<150~\mu$Jy/beam over 60\% of the area. The ${\rm rms}$ achieved in the inner 9.6 square degree area, enclosing the XMM-LSS field, is about $200~\mu$Jy/beam, while that over the outer 12.66 square degree area (which excludes the noisy edges) is about $45~\mu$Jy/beam. The resolution of the final mosaic is 6.5 arcsec. We present a catalogue of 5\,434 sources detected at $\geq7\times {\rm rms}$. We verify, and correct the reliability of, the catalog in terms of astrometry, flux, and false detection rate. Making use of the (to date) deepest radio continuum survey over a relatively large (2 square degree) field, complete at the flux levels probed by the GMRT-XXL-N survey, we also assess the survey's incompleteness as a function of flux density. The radio continuum sensitivity reached over a large field with a wealth of multi-wavelength data available makes the GMRT-XXL-N 610 MHz survey an important asset for studying the physical properties, environments and cosmic evolution of radio sources, in particular radio-selected active galactic nuclei (AGN).

We present an investigation of radio luminosity functions (LFs) and number counts based on the Karl G. Jansky Very Large Array-COSMOS 3 GHz Large Project. The radio-selected sample of 7826 galaxies with robust optical/near-infrared counterparts with excellent photometric coverage allows us to construct the total radio LF since z~5.7. Using the Markov chain Monte Carlo algorithm, we fit the redshift dependent pure luminosity evolution model to the data and compare it with previously published VLA-COSMOS LFs obtained on individual populations of radio-selected star-forming galaxies and galaxies hosting active galactic nuclei classified on the basis of presence or absence of a radio excess with respect to the star-formation rates derived from the infrared emission. We find they are in excellent agreement, thus showing the reliability of the radio excess method in selecting these two galaxy populations at radio wavelengths. We study radio number counts down to submicrojansky levels drawn from different models of evolving LFs. We show that our evolving LFs are able to reproduce the observed radio sky brightness, even though we rely on extrapolations toward the faint end. Our results also imply that no new radio-emitting galaxy population is present below 1 ujy. Our work suggests that selecting galaxies with radio flux densities between 0.1 and 10 ujy will yield a star-forming galaxy in 90-95 % of the cases with a high percentage of these galaxies existing around a redshift of z~2, thus providing useful constraints for planned surveys with the Square Kilometer Array and its precursors.

This work is part of a series of studies focusing on the environment and the properties of the X-ray selected active galactic nuclei (AGN) population from the XXL survey. The present survey, given its large area, continuity, extensive multiwavelength coverage, and large-scale structure information, is ideal for this kind of study. Here, we focus on the XXL-South (XXL-S) field. Our main aim is to study the environment of the various types of X-ray selected AGN and investigate its possible role in AGN triggering and evolution. We studied the large-scale (>1 Mpc) environment up to redshift z=1 using the nearest neighbour distance method to compare various pairs of AGN types. We also investigated the small-scale environment (<0.4 Mpc) by calculating the local overdensities of optical galaxies. In addition, we built a catalogue of AGN concentrations with two or more members using the hierarchical clustering method and we correlated them with the X-ray galaxy clusters detected in the XXL survey. It is found that radio detected X-ray sources are more obscured than non-radio ones, though not all radio sources are obscured AGN. We did not find any significant differences in the large-scale clustering between luminous and faint X-ray AGN, or between obscured and unobscured ones, or between radio and non-radio sources. At local scales (<0.4 Mpc), AGN typically reside in overdense regions, compared to non-AGN; however, no differences were found between the various types of AGN. A majority of AGN concentrations with two or more members are found in the neighbourhood of X-ray galaxy clusters within <25-45 Mpc. Our results suggest that X-ray AGN are typically located in supercluster filaments, but they are also found in over- and underdense regions.

The infrared-radio correlation (IRRC) offers a way to assess star formation from radio emission. Multiple studies found the IRRC to decrease with increasing redshift. This may in part be due to the lack of knowledge about the possible radio spectral energy distributions (SEDs) of star-forming galaxies. We constrain the radio SED of a complete sample of highly star-forming galaxies ($SFR>100\,\mathrm{M_{\odot}}/\,\mathrm{yr}$) based on the VLA-COSMOS $1.4\,\mathrm{GHz}$ Joint and $3\,\mathrm{GHz}$ Large Project catalogs. We reduce archival GMRT $325\,\mathrm{MHz}$ and $610\,\mathrm{MHz}$ observations, broadening the rest-frame frequency range to $0.3-15\,\mathrm{GHz}$. Employing survival analysis and fitting a double power law SED, we find that the slope steepens from a spectral index of $\alpha_1=0.51\pm 0.04$ below $4.5\,\mathrm{GHz}$ to $\alpha_2=0.98\pm0.07$ above $4.5\,\mathrm{GHz}.$ Our results suggest that the use of a K-correction assuming a single power-law radio SED for star forming galaxies is likely not the root cause of the IRRC trend.

We study the moderate-to-high radiative luminosity active galactic nuclei (HLAGN) within the VLA-COSMOS 3 GHz Large Project. The survey covers 2.6 square degrees centered on the COSMOS field with a 1$\sigma$ sensitivity of 2.3 $\mathrm{\mu Jy}$/beam across the field. This provides the simultaneously largest and deepest radio continuum survey available to date with exquisite multi-wavelength coverage. The survey yields 10,830 radio sources with signal-to-noise ratios $\geq$5. A subsample of 1,604 HLAGN is analyzed here. These were selected via a combination of X-ray luminosity and mid-infrared colors. We derive luminosity functions for these AGN and constrain their cosmic evolution out to a redshift of $z\sim6$, for the first time decomposing the star formation and AGN contributions to the radio continuum emission in the AGN. We study the evolution of number density and luminosity density finding a peak at $z\sim1.5$ followed by a decrease out to a redshift $z\sim6$.

Using infrared data from the Herschel Space Observatory and Karl G. Jansky Very Large Array (VLA) 3 GHz observations in the COSMOS field, we investigate the redshift evolution of the infrared-radio correlation (IRRC) for star-forming galaxies (SFGs) we classify as either spheroid- or disc-dominated based on their morphology. The sample predominantly consists of disc galaxies with stellar mass ${\gtrsim}10^{10}\,M_{\odot}$, and residing on the star-forming main sequence (MS). After the removal of AGN using standard approaches, we observe a significant difference between the redshift-evolution of the median IR/radio ratio $\overline{q}_{\mathrm{TIR}}$ of (i) a sample of ellipticals, plus discs with a substantial bulge component (`spheroid-dominated' SFGs) and, (ii) virtually pure discs and irregular systems (`disc-dominated' SFGs). The spheroid-dominated population follows a declining $\overline{q}_{\mathrm{TIR}}$ vs. $z$ trend similar to that measured in recent evolutionary studies of the IRRC. However, for disc-dominated galaxies, where radio and IR emission should be linked to star formation in the most straightforward way, we measure very little change in $\overline{q}_{\mathrm{TIR}}$. This suggests that low-redshift calibrations of radio emission as an SFR-tracer may remain valid out to at least $z\,{\simeq}\,1\,{-}\,1.5$ for pure star-forming systems. We find that the different redshift-evolution of $q_{\rm TIR}$ for the spheroid- and disc-dominated sample is mainly due to an increasing radio excess for spheroid-dominated galaxies at $z\,{\gtrsim}\,$0.8, hinting at some residual AGN activity in these systems. This finding demonstrates that in the absence of AGN the IRRC is independent of redshift, and that radio observations can therefore be used to estimate SFRs at all redshifts for genuinely star-forming galaxies.

We show that the most distant X-ray detected cluster known to date, ClJ1001 at z=2.506, hosts a strong overdensity of radio sources. Six of them are individually detected (within 10") in deep 0.75" resolution VLA 3GHz imaging, with S(3GHz)>8uJy. Of the six, AGN likely affects the radio emission in two galaxies while star formation is the dominant source powering the remaining four. We searched for cluster candidates over the full COSMOS 2-square degree field using radio-detected 3GHz sources and looking for peaks in Sigma5 density maps. ClJ1001 is the strongest overdensity by far with >10sigma, with a simple z_phot>1.5 preselection. A cruder photometric rejection of z<1 radio foregrounds leaves ClJ1001 as the second strongest overdensity, while even using all radio sources ClJ1001 remains among the four strongest projected overdensities. We conclude that there are great prospects for future, deep and wide-area radio surveys to discover large samples of the first generation of forming galaxy clusters. In these remarkable structures widespread star formation and AGN activity of massive galaxy cluster members, residing within the inner cluster core, will ultimately lead to radio continuum as one of the most effective means for their identification, with detection rates expected in the ballpark of 0.1-1 per square degree at z>2.5. Samples of hundreds such high-redshift clusters could potentially constrain cosmological parameters and test cluster and galaxy formation models.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

Based on a sample of over 1,800 radio AGN at redshifts out to z~5, which have typical stellar masses within ~3x(10^{10}-10^{11}) Msol, and 3 GHz radio data in the COSMOS field, we derived the 1.4 GHz radio luminosity functions for radio AGN (L_1.4GHz ~ 10^{22}-10^{27} W/Hz) out to z~5. We constrained the evolution of this population via continuous models of pure density and pure luminosity evolutions, and we found best-fit parametrizations of Phi*~(1+z)^{(2.00+/-0.18)-(0.60+/-0.14)z}, and L*~(1+z)^{(2.88+/-0.82)-(0.84+/-0.34)z}, respectively, with a turnover in number and luminosity densities of the population at z~1.5. We converted 1.4 GHz luminosity to kinetic luminosity taking uncertainties of the scaling relation used into account. We thereby derived the cosmic evolution of the kinetic luminosity density provided by the AGN and compared this luminosity density to the radio-mode AGN feedback assumed in the Semi-Analytic Galaxy Evolution (SAGE) model, i.e., to the redshift evolution of the central supermassive black hole accretion luminosity taken in the model as the source of heating that offsets the energy losses of the cooling, hot halo gas, and thereby limits further stellar mass growth of massive galaxies. We find that the kinetic luminosity exerted by our radio AGN may be high enough to balance the radiative cooling of the hot gas at each cosmic epoch since z~5. However, although our findings support the idea of radio-mode AGN feedback as a cosmologically relevant process in massive galaxy formation, many simplifications in both the observational and semi-analytic approaches still remain and need to be resolved before robust conclusions can be reached.

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

(abridged) We study the composition of the faint radio population selected from the VLA-COSMOS 3GHz Large Project. The survey covers a 2.6sq.deg. area with a mean rms of ~2.3uJy/b, cataloging 10830 sources (>5sigma). Combining these radio data with optical, near-infrared (UltraVISTA), mid-infrared (Spitzer/IRAC) data, and X-ray data (Chandra), we find counterparts to radio sources for ~93% of the radio sample (in the areas of the COSMOS field not affected by saturated or bright sources in the optical to NIR bands), reaching out to z<6. We further classify the sources as star forming galaxies or AGN based on various criteria, such as X-ray luminosity, observed MIR color, UV-FIR spectral-energy distribution, rest-frame NUV-optical color corrected for dust extinction, and radio-excess relative to that expected from the hosts' star-formation rate. We separate the AGN into sub-samples dominated by low-to-moderate and moderate-to-high radiative luminosity AGN, candidates for high-redshift analogues to local low- and high-excitation emission line AGN, respectively. We study the fractional contributions of these sub-populations down to radio flux levels of ~11uJy at 3GHz (or ~20uJy at 1.4GHz assuming a spectral index of -0.7). We find that the dominant fraction at 1.4GHz flux densities above ~200uJy is constituted of low-to-moderate radiative luminosity AGN. Below densities of ~100uJy the fraction of star-forming galaxies increases to ~60%, followed by the moderate-to-high radiative luminosity AGN (~20%), and low-to-moderate radiative luminosity AGN (~20%). Based on this observational evidence, we extrapolate the fractions down to sensitivities of the SKA. Our estimates suggest that at the faint flux limits to be reached by the SKA1 surveys, a selection based only on radio flux limits can provide a simple tool to efficiently identify samples highly (>75%) dominated by star-forming galaxies.

We explore the multiwavelength properties of AGN host galaxies for different classes of radio-selected AGN out to z$\lesssim$6 via a multiwavelength analysis of about 7700 radio sources in the COSMOS field. The sources were selected with the Very Large Array (VLA) at 3 GHz (10 cm) within the VLA-COSMOS 3 GHz Large Project, and cross-matched with multiwavelength ancillary data. This is the largest sample of high-redshift (z$\lesssim$6) radio sources with exquisite photometric coverage and redshift measurements available. We constructed a sample of moderate-to-high radiative luminosity AGN (HLAGN) via spectral energy distribution (SED) decomposition combined with standard X-ray and mid-infrared diagnostics. Within the remainder of the sample we further identified low-to-moderate radiative luminosity AGN (MLAGN) via excess in radio emission relative to the star formation rates in their host galaxies. We show that AGN power in HLAGN occurs predominantly in radiative form, while MLAGN display a substantial mechanical AGN luminosity component. We found significant differences in the host properties of the two AGN classes, as a function of redshift. At z$<$1.5, MLAGN appear to reside in significantly more massive and less star-forming galaxies compared to HLAGN. At z$>$1.5, we observed a reversal in the behaviour of the stellar mass distributions with the HLAGN populating the higher stellar mass tail. We interpret this finding as a possible hint of the downsizing of galaxies hosting HLAGN, with the most massive galaxies triggering AGN activity earlier than less massive galaxies, and then fading to MLAGN at lower redshifts. Our conclusion is that HLAGN and MLAGN samples trace two distinct galaxy and AGN populations in a wide range of redshifts, possibly resembling the radio AGN types often referred to as radiative- and jet-mode (or high- and low-excitation), respectively.

We make use of the deep Karl G. Jansky Very Large Array (VLA) COSMOS radio observations at 3 GHz to infer radio luminosity functions of star-forming galaxies up to redshifts of z~5 based on approximately 6000 detections with reliable optical counterparts. This is currently the largest radio-selected sample available out to z~5 across an area of 2 square degrees with a sensitivity of rms=2.3 ujy/beam. By fixing the faint and bright end shape of the radio luminosity function to the local values, we find a strong redshift trend that can be fitted with a pure luminosity evolution L~(1+z)^{(3.16 +- 0.2)-(0.32 +- 0.07) z}. We estimate star formation rates (SFRs) from our radio luminosities using an infrared (IR)-radio correlation that is redshift dependent. By integrating the parametric fits of the evolved luminosity function we calculate the cosmic SFR density (SFRD) history since z~5. Our data suggest that the SFRD history peaks between 2<z<3 and that the ultraluminous infrared galaxies (ULIRGs; 100 Msol/yr<SFR<1000 Msol/yr) contribute up to ~25% to the total SFRD in the same redshift range. Hyperluminous infrared galaxies (HyLIRGs; SFR>1000 Msol/yr) contribute an additional <2% in the entire observed redshift range. We find evidence of a potential underestimation of SFRD based on ultraviolet (UV) rest-frame observations of Lyman break galaxies (LBGs) at high redshifts (z>4) on the order of 15-20%, owing to appreciable star formation in highly dust-obscured galaxies, which might remain undetected in such UV observations.

The 2.1 GHz radio source catalogue of the 25 deg$^2$ ultimate XMM extragalactic survey south (XXL-S) field, observed with the Australia Telescope Compact Array (ATCA), is presented. The final radio mosaic achieved a resolution of $\sim$$4.8"$ and a median rms noise of $\sigma \approx41$ $\mu$Jy/beam. To date, this is the largest area radio survey to reach this flux density level. A total of 6350 radio components above 5$\sigma$ are included in the component catalogue, 26.4% of which are resolved. Of these components, 111 were merged together to create 48 multiple-component radio sources, resulting in a total of 6287 radio sources in the source catalogue, 25.9% of which were resolved. A survival analysis revealed that the median spectral index of the Sydney University Molonglo Sky Survey (SUMSS) 843 MHz sources in the field is $\alpha$ = $-$0.75, consistent with the values of $-0.7$ to $-0.8$ commonly used to characterise radio spectral energy distributions of active galactic nuclei (AGN). The 2.1 GHz and 1.4 GHz differential radio source counts are presented and compared to other 1.4 GHz radio surveys. The XXL-S source counts show good agreement with the other surveys.

We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys.

We present observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz toward a sub-field of the XXL-North 25 deg$^{2}$ field targeting the first supercluster discovered in the XXL Survey. The structure has been found at a spectroscopic redshift of 0.43 and extending over $0.35\times0.1$ deg$^{2}$ on the sky. We present the 3 GHz VLA radio continuum observations, the radio mosaic and radio source catalogue, and, perform an analysis of the supercluster in the optical and radio regimes using photometric redshifts from the CFHTLS survey and our new VLA data. Our 3 GHz radio mosaic has a resolution of $3.2\times1.9$ arcsec$^2$, and covers an area of $41\times41$ arcmin$^2$ with rms noise level below $\sim20 \mu$Jy beam$^{-1}$ ($\sim11\mu$Jy beam$^{-1}$ in central $15\times15$ arcmin$^2$). We extract a catalogue of 155 radio sources $\geq6\sigma$, 8 of which are multicomponent sources, and 123 ($79\%$) can be associated with optical sources in the CFHTLS W1 catalogue. Applying Voronoi tessellation analysis (VTA) in the area around the X-ray identified supercluster we identify a total of 17 overdensities at $z=0.35-0.50$, 7 of which are associated with clusters detected in the $XMM-Newton$ XXL data. We find a mean photometric redshift of 0.43 for our overdensities, consistent with the spectroscopic redshifts of 7 X-ray detected clusters. The full VTA-identified structure extends over $\sim0.6\times0.2$ deg$^2$ on the sky, ($\sim12\times4$ Mpc$^{2}$ at $z=0.43$). We associate 8 radio sources with potential group/cluster member galaxies none of which are large radio galaxies. The spatial distribution of the red and blue potential group member galaxies suggests that the clusters are dynamically young and not virialised yet, as expected for hierarchical structure growth in a $\Lambda$CDM universe. Further spectroscopic data are required to analyze the dynamical state of the groups.

We study the composition of the faint radio population selected from the VLA-COSMOS 3GHz Large Project, a radio continuum survey performed at 10 cm wavelength. The survey covers the full 2 square degree COSMOS field with mean $rms\sim2.3$ $\mu$Jy/beam, cataloging 10,899 source components above $5\times rms$. By combining these radio data with UltraVISTA, optical, near-infrared, and Spitzer/IRAC mid-infrared data, as well as X-ray data from the Chandra Legacy, and Chandra COSMOS surveys, we gain insight into the emission mechanisms within our radio sources out to redshifts of $z\sim5$. From these emission characteristics we classify our souces as star forming galaxies or AGN. Using their multi-wavelength properties we further separate the AGN into sub-samples dominated by radiatively efficient and inefficient AGN, often referred to as high- and low-excitation emission line AGN. We compare our method with other results based on fitting of the sources' spectral energy distributions using both galaxy and AGN spectral models, and those based on the infrared-radio correlation. We study the fractional contributions of these sub-populations down to radio flux levels of $\sim$10 $\mu$Jy. We find that at 3 GHz flux densities above $\sim$400 $\mu$Jy quiescent, red galaxies, consistent with the low-excitation radio AGN class constitute the dominant fraction. Below densities of $\sim$200 $\mu$Jy star-forming galaxies begin to constitute the largest fraction, followed by the low-excitation, and X-ray- and IR-identified high-excitation radio AGN.

We present 2.1 GHz imaging with the Australia Telescope Compact Array (ATCA) of a 6.5 deg^2 region within the XXM-Newton XXL South field using a band of 1.1-3.1 GHz. We achieve an angular resolution of 4.7" x 4.2" in the final radio continuum map with a median rms noise level of 50 uJy/beam. We identify 1389 radio sources in the field with peak S/N >=5 and present the catalogue of observed parameters. We find that 305 sources are resolved, of which 77 consist of multiple radio components. These number counts are in agreement with those found for the COSMOS-VLA 1.4 GHz survey. We derive spectral indices by a comparison with the Sydney University Molongolo Sky Survey (SUMSS) 843MHz data. We find an average spectral index of -0.78 and a scatter of 0.28, in line with expectations. This pilot survey was conducted in preparation for a larger ATCA program to observe the full 25 deg^2 southern XXL field. When complete, the survey will provide a unique resource of sensitive, wide-field radio continuum imaging with complementary X-ray data in the field. This will facilitate studies of the physical mechanisms of radio-loud and radio-quiet AGNs and galaxy clusters, and the role they play in galaxy evolution. The source catalogue is publicly available online via the XXL Master Catalogue browser and the Centre de Donn\'ees astronomiques de Strasbourg (CDS).

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.

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

Using the Parkes radio telescope, we study the 21cm neutral hydrogen (HI) properties of a sample of galaxies with redshifts z<0.13 extracted from the optical 2dF Galaxy Redshift Survey (2dFGRS). Galaxies at 0.04<z<0.13 are studied using new Parkes observations of a 42deg^2 field near the South Galactic Pole (SGP). A spectral stacking analysis of the 3,277 2dFGRS objects within this field results in a convincing 12sigma detection. For the low-redshift sample at 0<z<0.04, we use the 15,093 2dFGRS galaxies observed by the HI Parkes All-Sky Survey (HIPASS) and find a 31sigma stacked detection. We measure average HI masses of (6.93 +/- 0.17)*10^9 h^{-2} Msun and (1.48 +/- 0.03)*10^9 h^{-2} Msun for the SGP and HIPASS samples, respectively. Accounting for source confusion and sample bias, we find a cosmic HI mass density of Omega_HI=(3.19_{-0.59} ^{+0.43})*10^{-4} h^{-1} for the SGP sample and (2.82_{-0.59}^{+0.30})*10^{-4} h^{-1} for the HIPASS sample. This suggests no (12 +/- 23%) evolution in the cosmic HI density over the last ~1h^{-1} Gyr. Due to the very large effective volumes, cosmic variance in our determination of Omega_HI is considerably lower than previous estimates. Our stacking analysis reproduces and quantifies the expected trends in the HI mass and mass-to-light ratio of galaxies with redshift, luminosity and colour.

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.