J. Delhaize

The XXL Survey I. Scientific motivations - XMM-Newton observing plan - Follow-up observations and simulation programme

Context. The quest for the cosmological parameters that describe our universe continues to motivate the scientific community to undertake very large survey initiatives across the electromagnetic spectrum. Over the past two decades, the Chandra and XMM-Newton observatories have supported numerous studies of X-ray-selected clusters of galaxies, active galactic nuclei (AGNs), and the X-ray background. The present paper is the first in a series reporting results of the XXL-XMM survey; it comes at a time when the Planck mission results are being finalised. Aims. 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 deg(2) each at a point-source sensitivity of similar to 5 x 10(-15) erg s(-1) cm(-2) in the [0.5-2] keV band (completeness limit). The surveys 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. Methods. 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. Results. 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. Conclusions. 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 VLA-COSMOS 3 GHz Large Project: Multiwavelength counterparts and the composition of the faint radio population

We study the composition of the faint radio population selected from the Karl G. Jansky Very Large Array Cosmic Evolution Survey (VLA-COSMOS) 3 GHz Large Project, which is a radio continuum survey performed at 10 cm wavelength. The survey covers a 2.6 square degree area with a mean rms of ~ 2.3 μJy/beam, cataloging 10 830 sources above 5σ, and enclosing the full 2 square degree COSMOS field. By combining these radio data with optical, near-infrared (UltraVISTA), and mid-infrared (Spitzer/IRAC) data, as well as X-ray data (Chandra), we find counterparts to radio sources for ~93% of the total radio sample reaching out to z ≲ 6; these sources are found in the unmasked areas of the COSMOS field, i.e., those not affected by saturated or bright sources in the optical to near-infrared (NIR) bands. We further classify the sources as star-forming galaxies or AGN based on various criteria, such as X-ray luminosity; observed mid-infrared color; UV–far-infrared spectral energy distribution; rest-frame, near-UV optical color that is corrected for dust extinction; and radio excess relative to that expected from the star formation rate of the hosts. We separate the AGN into subsamples dominated by low-to-moderate and moderate-to-high radiative luminosity AGN, i.e., candidates for high-redshift analogs to local low- and high-excitation emission line AGN, respectively. We study the fractional contributions of these subpopulations down to radio flux levels of ~11 μJy at 3 GHz (or ~20 μJy at 1.4 GHz assuming a spectral index of –0.7). We find that the dominant fraction at 1.4 GHz flux densities above ~200 μJy is constituted of low-to-moderate radiative luminosity AGN. Below densities of ~100 μJy 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 Square Kilometer Array (SKA). Our estimates suggest that at the faint flux limits to be reached by the (Wide, Deep, and UltraDeep) 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.

The VLA-COSMOS 3 GHz Large Project: The infrared-radio correlation of star-forming galaxies and AGN to z ≲ 6

We examine the behaviour of the infrared-radio correlation (IRRC) over the range 0 <z ≲ 6 using new, highly sensitive 3 GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared data from the Herschel Space Observatory in the 2 deg^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.4 GHz radio luminosity ratio (q_(TIR)), decreases with increasing redshift: q_(TIR)(z) = (2.88 ± 0.03)(1 + z)^(−0.19 ± 0.01). This is consistent with several previous results from the literature. Moderate-to-high radiative luminosity AGN do not follow the same q_(TIR)(z) trend as star-forming galaxies, 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_(TIR)(z) trend of the star-forming galaxy population. We demonstrate that the choice of the average radio spectral index directly affects the normalisation, as well as the derived trend with redshift of the IRRC. An increasing fractional contribution to the observed 3 GHz 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 trend of decreasing q_(TIR) with increasing redshift. A more detailed understanding of the radio spectrum is therefore required for robust K corrections in the radio and to fully understand the cosmic evolution of the IRRC. Lastly, we present a redshift-dependent relation between rest-frame 1.4 GHz radio luminosity and star formation rate taking the derived redshift trend into account.

The VLA-COSMOS 3 GHz Large Project: Cosmic evolution of radio AGN and implications for radio- mode feedback since z 5

Based on a sample of over 1800 radio AGN at redshifts out to z 5, which have typical stellar masses within 3 × (1010 - 1011)M⊙, and 3 GHz radio data in the COSMOS field, we derived the 1.4 GHz radio luminosity functions for radio AGN (L1.4 GHz 1022 - 1027 W Hz-1) 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 Φ∗ ∝ (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.

SMBH accretion properties of radio-selected AGN out to z ∼ 4

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 (LAGN1.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σ) 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 LAGN1.4 over different redshift ranges, out to z~4, and use X-ray stacking to calculate the average s-BHAR in each LAGN1.4-z bin. We find that the average s-BHAR is independent of LAGN1.4, at all redshifts. However, we see a strong increase of s-BHAR with redshift, at fixed LAGN1.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 LAGN1.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.

Constraints on submicrojansky radio number counts based on evolving VLA-COSMOS luminosity functions

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.

The VLA-COSMOS 3 GHz large project: the cosmic star-formation history since z ~ 5

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 μJy beam^(-1). 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.4 GHz) ∝ (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 1000 M⊙ yr^(-1)) 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 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 VLA-COSMOS 3 GHz Large Project: Cosmic star formation history since z ~ 5

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 μJy beam^(-1). 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.4 GHz) ∝ (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 1000 M⊙ yr^(-1)) 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 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 XXL Survey - XI. ATCA 2.1 GHz continuum observations

We present 2.1 GHz imaging with the Australia Telescope Compact Array (ATCA) of a 6.5 deg2 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″ × 4:2″ in the final radio continuum map with a median rms noise level of 50 μJy/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 deg2 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 Donnees astronomiques de Strasbourg (CDS).

The XXL Survey. XXXI. Classification and host galaxy properties of 2.1 GHz ATCA XXL-S radio sources

The classification of the host galaxies of the radio sources in the 25 deg