Mladen Novak

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.

New insights from deep VLA data on the potentially recoiling black hole CID-42 in the COSMOS field

We present deep 3-GHz Karl G. Jansky Very Large Array (VLA) observations of the potentially recoiling black hole CID-42 in the Cosmic Evolution Survey (COSMOS) field. This galaxy shows two optical nuclei in the Hubble Space Telescope/Advanced Camera for Surveys (HST/ACS) image and a large velocity offset of ≈1300 km s^(−1) between the broad and narrow Hβ emission line although the spectrum is not spacially resolved (Civano et al. 2010). The new 3 GHz VLA data have a bandwidth of 2 GHz and to correctly interpret the flux densities imaging was done with two different methods: multiscale multifrequency (MSMF) synthesis and spectral windows (SPWs) stacking. The final resolutions and sensitivities of these maps are 0.7 arcsec with rms = 4.6 μJy beam^(−1) and 0.9 arcsec with rms = 4.8 μJy beam^(−1), respectively. With a 7σ detection, we find that the entire observed 3-GHz radio emission can be associated with the south-eastern component of CID-42, coincident with the detected X-ray emission. We use our 3 GHz data combined with other radio data from the literature ranging from 320 MHz to 9 GHz, which include the VLA, Very Long Baseline Array (VLBA) and Giant Metrewave Radio Telescope (GMRT) data, to construct a radio synchrotron spectrum of CID-42. The radio spectrum suggests a type I unobscured radio-quiet flat-spectrum active galactic nucleus (AGN) in the south-eastern component which may be surrounded by a more extended region of old synchrotron electron population or shocks generated by the outflow from the supermassive black hole (SMBH). Our data are consistent with the recoiling black hole picture but cannot rule out the presence of an obscured and radio-quiet SMBH in the north-western component.

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.

(Sub)millimetre interferometric imaging of a sample of COSMOS/AzTEC submillimetre galaxies. IV. Physical properties derived from spectral energy distributions

We characterise the physical nature of a 1.1 mm-selected, flux-limited, and interferometrically followed up sample of SMGs in COSMOS. We used the MAGPHYS code to fit the multiwavelength (UV-radio) SEDs of 16 of the target SMGs. We also constructed the pure radio SEDs of our SMGs using three different radio bands (325 MHz, 1.4 GHz, and 3 GHz). Moreover, since two SMGs in our sample, AzTEC1 and AzTEC3, benefit from previous CO line observations, we studied their properties in more detail. We found that 63% of our target SMGs lie above the galaxy main-sequence by more than a factor of 3, and hence are starbursts. The 3 GHz radio sizes we have previously measured for the target SMGs were compared with the present stellar mass estimates, and we found that the z>3 SMGs are fairly consistent with the mass-size relationship of z~2 compact, quiescent galaxies (cQGs). The median IR-radio correlation parameter is found to be q=2.27, which is lower than measured locally (median q=2.64). AzTEC1 is found to have a sub-Eddington SFR surface density (by a factor of 2.6), while AzTEC3 appears to be an Eddington-limited starburster. The gas reservoir in these two high-z SMGs would be exhausted in only ~86 and 19 Myr at the current SFR, respectively. A comparison of the MAGPHYS-based properties of our SMGs with those of equally bright ALESS SMGs suggests that the two populations share fairly similar physical characteristics, including the q parameter. A hint of negative correlation is found between the 3 GHz size and the level of starburstiness, and hence cosmic-ray electrons in more compact starbursts might be more susceptible to free-free absorption. Some of the derived low and high q values (compared to the local median) could be the result of a specific merger/post-starburst phase of galaxy evolution. Overall, our results support the scenario where z>3 SMGs evolve into todays giant ellipticals.

An ALMA survey of submillimetre galaxies in the COSMOS field: The extent of the radio-emitting region revealed by 3 GHz imaging with the Very Large Array

Context. The observed spatial scale of the radio continuum emission from star-forming galaxies can be used to investigate the spatial extent of active star formation, constrain the importance of cosmic-ray transport, and examine the effects of galaxy interactions. Aims. We determine the radio size distribution of a large sample of 152 submillimetre galaxies (SMGs) in the COSMOS field that were pre-selected at 1.1 mm, and later detected with the Atacama Large Millimetre/submillimetre Array (ALMA) in the observed-frame 1.3 mm dust continuum emission at a signal-to-noise ratio (S/N) of ≥5. Methods. We used the deep, subarcsecond-resolution (1σ = 2.3μJy beam^(-1); . 75) centimetre radio continuum observations taken by the Karl G. Jansky Very Large Array (VLA)-COSMOS 3 GHz Large Project. Results. One hundred and fifteen of the 152 target SMGs (76% ± 7%) were found to have a 3 GHz counterpart (≥ 4.2σ), which renders the radio detection rate notably high. The median value of the deconvolved major axis full width at half maximum (FWHM) size at 3 GHz is derived to be 0.59 ± 0.05 , or 4.6 ± 0.4 kpc in physical units, where the median redshift of the sources is z = 2.23 ± 0.13 (23% are spectroscopic and 77% are photometric values). The radio sizes are roughly log-normally distributed, and they show no evolutionary trend with redshift, or difference between different galaxy morphologies. We also derived the spectral indices between 1.4 and 3 GHz, and 3 GHz brightness temperatures for the sources, and the median values were found to be α_(1.4 GHz)^(3 GHz) = -0.67 (S_ν ∝ ν^α) and T_B = 12.6 ± 2 K. Three of the target SMGs, which are also detected with the Very Long Baseline Array (VLBA) at 1.4 GHz (AzTEC/C24b, 61, and 77a), show clearly higher brightness temperatures than the typical values, reaching T_B(3 GHz) > 10^(4.03) K for AzTEC/C61. Conclusions. The derived median radio spectral index agrees with a value expected for optically thin non-thermal synchrotron radiation, and the low median 3 GHz brightness temperature shows that the observed radio emission is predominantly powered by star formation and supernova activity. However, our results provide a strong indication of the presence of an active galactic nucleus in the VLBA and X-ray-detected SMG AzTEC/C61 (high TB and an inverted radio spectrum). The median radio-emitting size we have derived is ~ 1.5–3 times larger than the typical far-infrared dust-emitting sizes of SMGs, but similar to that of the SMGs’ molecular gas component traced through mid-J line emission of carbon monoxide. The physical conditions of SMGs probably render the diffusion of cosmic-ray electrons inefficient, and hence an unlikely process to lead to the observed extended radio sizes. Instead, our results point towards a scenario where SMGs are driven by galaxy interactions and mergers. Besides triggering vigorous starbursts, galaxy collisions can also pull out the magnetised fluids from the interacting disks, and give rise to a taffy-like synchrotron-emitting bridge. This provides an explanation for the spatially extended radio emission of SMGs, and can also cause a deviation from the well-known infrared-radio correlation owing to an excess radio emission. Nevertheless, further high-resolution observations are required to examine the other potential reasons for the very compact dust-emitting sizes of SMGs, such as the radial dust temperature and metallicity gradients.

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.

An ALMA survey of submillimeter galaxies in the COSMOS field: Multiwavelength counterparts and redshift distribution

We carried out targeted ALMA observations of 129 fields in the COSMOS region at 1.25 mm, detecting 152 galaxies at S/N ≥ 5with an average continuum RMS of 150 μJy. These fields represent a S/N-limited sample of AzTEC/ASTE sources with 1.1 mm S/N≥ 4 over an area of 0.72 square degrees. Given ALMA’s fine resolution and the exceptional spectroscopic and multiwavelength photometric data available in COSMOS, this survey allows us unprecedented power in identifying submillimeter galaxy counterparts and determining their redshifts through spectroscopic or photometric means. In addition to 30 sources with prior spectroscopic redshifts, we identified redshifts for 113 galaxies through photometric methods and an additional nine sources with lower limits, which allowed a statistically robust determination of the redshift distribution. We have resolved 33 AzTEC sources into multi-component systems and our redshifts suggest that nine are likely to be physically associated. Our overall redshift distribution peaks at z ~ 2.0 with a high-redshift tail skewing the median redshift to z = 2.48 ± 0.05. We find that brighter millimeter sources are preferentially found at higher redshifts. Our faintestsources, with S_(1.25 mm) 1.8 mJy, have a median redshift of z = 3.08 ± 0.17. After accounting for spectral energy distribution shape and selection effects, these results are consistent with several previous submillimeter galaxy surveys, and moreover, support the conclusion that the submillimeter galaxy redshift distribution is sensitive to survey depth.

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