Search

We present a multi-line survey of the interstellar medium (ISM) in two $z>6$ quasar (QSO) host galaxies, PJ231-20 ($z=6.59$) and PJ308-21 ($z=6.23$), and their two companion galaxies. Observations were carried out using the Atacama Large (sub-)Millimeter Array (ALMA). We targeted eleven transitions including atomic fine structure lines (FSLs) and molecular lines: [NII]$_{\rm 205\mu m}$, [CI]$_{\rm 369\mu m}$, CO ($J_{\rm up} = 7, 10, 15, 16$), H$_2$O $3_{12}-2_{21}$, $3_{21}-3_{12}$, $3_{03}-2_{12}$, and the OH$_{\rm 163\mu m}$ doublet. The underlying far-infrared (FIR) continuum samples the Rayleigh-Jeans tail of the respective dust emission. By combining this information with our earlier ALMA [CII]$_{\rm 158\mu m}$ observations, we explore the effects of star formation and black hole feedback on the galaxies' ISM using the CLOUDY radiative transfer models. We estimate dust masses, spectral indexes, IR luminosities, and star-formation rates from the FIR continuum. The analysis of the FSLs indicates that the [CII]$_{\rm 158\mu m}$ and [CI]$_{\rm 369\mu m}$ emission arises predominantly from the neutral medium in photodissociation regions (PDRs). We find that line deficits are in agreement with those of local luminous infrared galaxies. The CO spectral line energy distributions (SLEDs), reveal significant high-$J$ CO excitation in both quasar hosts. Our CO SLED modeling of the quasar PJ231-20 shows that PDRs dominate the molecular mass and CO luminosities for $J_{\rm up}\le 7$, while the $J_{\rm up}\ge10$ CO emission is likely driven by X-ray dissociation regions produced by the active galactic nucleus (AGN) at the very center of the quasar host [abridged]., Comment: 26 pages, 17 figures, accepted for publication in A&A

We explore the kinematics of 27 z~6 quasar host galaxies observed in [CII]-158 micron ([CII]) emission with the Atacama Large Millimeter/sub-millimeter Array at a resolution of ~0.25''. We find that nine of the galaxies show disturbed [CII] emission, either due to a close companion galaxy or recent merger. Ten galaxies have smooth velocity gradients consistent with the emission arising from a gaseous disk. The remaining eight quasar host galaxies show no velocity gradient, suggesting that the gas in these systems is dispersion-dominated. All galaxies show high velocity dispersions with a mean of 129+-10 km/s. To provide an estimate of the dynamical mass within twice the half-light radius of the quasar host galaxy, we model the kinematics of the [CII] emission line using our publicly available kinematic fitting code, qubefit. This results in a mean dynamical mass of 5.0+-0.8(+-3.5) x 10^10 Msun. Comparison between the dynamical mass and the mass of the supermassive black hole reveals that the sample falls above the locally derived bulge mass--black hole mass relation at 2.4sigma significance. This result is robust even if we account for the large systematic uncertainties. Using several different estimators for the molecular mass, we estimate a gas mass fraction of >10%, indicating gas makes up a large fraction of the baryonic mass of z~6 quasar host galaxies. Finally, we speculate that the large variety in [CII] kinematics is an indication that gas accretion onto z~6 super massive black holes is not caused by a single precipitating factor., 33 pages. Kinematic fitting code is available at https://github.com/mneeleman/qubefit

The interactions between radio jets and the interstellar medium play a defining role for the co-evolution of central supermassive black holes and their host galaxies, but observational constraints on these feedback processes are still very limited at redshifts $z > 2$. We investigate the radio-loud quasar PSO J352.4034-15.3373 at $z \sim 6$ at the edge of the Epoch of Reionization. This quasar is among the most powerful radio emitters and the first one with direct evidence of extended radio jets ($\sim$1.6 kpc) at these high redshifts. We analyze NOEMA and ALMA millimeter data targeting the CO (6-5) and [CII] far-infrared emission lines, respectively, and the underlying continuum. The broad $440\pm 80$ km s$^{-1}$ and marginally resolved [CII] emission line yields a systemic redshift of $z\!=\!5.832 \pm 0.001$. Additionally, we report a strong 215 MHz radio continuum detection, $88\pm 7$ mJy, using the GMRT. This measurement significantly improves the constraints at the low-frequency end of the spectral energy distribution of this quasar. In contrast to what is typically observed in high-redshift radio-quiet quasars, we show that cold dust emission alone cannot reproduce the millimeter continuum measurements. This is evidence that the strong synchrotron emission from the quasar contributes substantially to the emission even at millimeter (far-infrared in the rest-frame) wavelengths. This quasar is an ideal system to probe the effects of radio jets during the formation of a massive galaxy within the first Gyr of the Universe., 14 pages, 5 figures

Due to an error in the analysis of the ALMA data, the [C II] line fluxes Fline in 6 out of the 13 analyzed quasars were underestimated, and thus also based the inferred star formation rate (SFR) on the [C II] emission. We corrected these values and show the updated measurements for the six affected quasars in the corrected Table 2. The dust continuum estimates and their derived quantities were not affected by the error. Additionally, the observed peak frequency of the [C II] emission line nobs and thus the inferred systemic redshifts of the quasars were mildly offset, causing a shift of Dnobs ? 0.02 GHz and - Dzsub mm ? 0.0005. The derived proximity zone sizes were hence affected byDRp ? 0.04 pMpc. All updated values are shown in the corrected Tables 2 and 4. These changes did not affect any of the conclusions of the paper, nor did the figures in the paper change in any noticeable way. Furthermore, since the publication of our paper the offset between the peak of the Mg II emission line and the systemic redshift derived from submillimeter emission lines estimated by Schindler et al. (2020) changed slightly in their published version compared to the preliminary value we had used. Thus we use the now published velocity offset of Dv(Mg II-[C II]) = - - 391+ km s- 455 256 1 to estimate the systemic redshifts of the three quasars in our sample, for which no submillimeter data are available and therefore the redshift estimate is based on the Mg II emission. This resulted in minor changes of Dz ? 0.001, which are updated in Table 4.

We studied a sample of 274 radio and X-ray selected quasars (XQSOs) detected in the COSMOS and XXL-S radio surveys at 3 GHz and 2.1 GHz, respectively. This sample was identified by adopting a conservative threshold in X-ray luminosity, Lx [2-10\ keV] >= 10^44 erg/s, selecting only the most powerful quasars. Using available multiwavelength data, we examined various criteria for the selection of radio-loud (RL) and radio-quiet (RQ) XQSOs, finding that the number of RL/RQ XQSOs changes significantly depending on the chosen criterion. This discrepancy arises due to the different criteria tracing different physical processes and due to our sample being selected from flux-limited radio and X-ray surveys. Another approach to study the origin of radio emission in XQSOs is via their radio luminosity function (RLF). We constructed the XQSO 1.4 GHz radio luminosity functions (RLFs) in six redshift bins at 0.5, 14 pages, 11 figures

We make use of sensitive (9.3 microJy/beam RMS) 1.2mm-continuum observations from the ASPECS ALMA large program of the Hubble Ultra Deep Field (HUDF) to probe dust-enshrouded star formation from 1362 Lyman-break galaxies spanning the redshift range z=1.5-10 (to ~7-28 Msolar/yr at 4 sigma over the entire range). We find that the fraction of ALMA-detected galaxies in our z=1.5-10 samples increases steeply with stellar mass, with the detection fraction rising from 0% at 10^9 Msolar to 85(-18)(+9)% at >10^{10} Msolar. Moreover, stacking all 1253 low-mass (10^{9.5} M_{solar} and a SMC-like relation at lower masses. Using stellar-mass and beta measurements for z~2 galaxies over CANDELS, we derive a new empirical relation between beta and stellar mass and then use this correlation to show that our IRX-beta and IRX-stellar mass relations are consistent with each other. We then use these constraints to express the infrared excess as a bivariate function of beta and stellar mass. Finally, we present updated estimates of star-formation rate density determinations at z>3, leveraging current improvements in the measured infrared excess and recent probes of ultra-luminous far-IR galaxies at z>2., 32 pages, 22 figures, 11 pages, accepted for publication in the Astrophysical Journal, key results shown in Figures 7, 8, 11, and 20, updated to match version in press

We present the discovery of PSO J083.8371+11.8482, a weak emission line quasar with extreme star formation rate at $z=6.3401$. This quasar was selected from Pan-STARRS1, UHS, and unWISE photometric data. Gemini/GNIRS spectroscopy follow-up indicates a MgII-based black hole mass of $M_\mathrm{BH}=\left(2.0^{+0.7}_{-0.4}\right)\times10^9~M_\odot$ and an Eddington ratio of $L_\mathrm{bol}/L_\mathrm{Edd}=0.5^{+0.1}_{-0.2}$, in line with actively accreting supermassive black hole (SMBH) at $z\gtrsim6$. HST imaging sets strong constraint on lens-boosting, showing no relevant effect on the apparent emission. The quasar is also observed as a pure point-source with no additional emission component. The broad line region (BLR) emission is intrinsically weak and not likely caused by an intervening absorber. We found rest-frame equivalent widths of EW(Ly$\alpha$+NV) $=5.7\pm0.7$ Angstrom, EW(CIV) $\leq5.8$ Angstrom (3-sigma upper limit), and EW(MgII) $=8.7\pm0.7$ Angstrom. A small proximity zone size ($R_\mathrm{p}=1.2\pm0.4$ pMpc) indicates a lifetime of only $t_\mathrm{Q}=10^{3.4\pm0.7}$ years from the last quasar phase ignition. ALMA shows extended [CII] emission with a mild velocity gradient. The inferred far-infrared luminosity ($L_\mathrm{FIR}=(1.2\pm0.1)\times10^{13}\,L_\odot$) is one of the highest among all known quasar hosts at $z\gtrsim6$. Dust and [CII] emissions put a constraint on the star formation rate of SFR $=900-4900~M_\odot\,\mathrm{yr^{-1}}$, similar to that of hyper-luminous infrared galaxy. Considering the observed quasar lifetime and BLR formation timescale, the weak-line profile in the quasar spectrum is most likely caused by a BLR which is not yet fully formed rather than continuum boosting by gravitational lensing or a soft continuum due to super-Eddington accretion., Comment: 28 pages, 16 figures, 3 tables, accepted for publication in ApJ

In a multi-wavelength survey of $13$ quasars at $5.8\lesssim z\lesssim6.5$, that were pre-selected to be potentially young, we find five objects with extremely small proximity zone sizes that may imply UV-luminous quasar lifetimes of $\lesssim 100,000$ years. Proximity zones are regions of enhanced transmitted flux in the vicinity of the quasars that are sensitive to the quasars' lifetimes because the intergalactic gas has a finite response time to their radiation. We combine sub-mm observations from the Atacama Large Millimetre Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA), as well as deep optical and near-infrared spectra from medium-resolution spectrograph on the Very Large Telescope (VLT) and on the Keck telescopes, in order to identify and characterize these new young quasars, which provide valuable clues about the accretion behavior of supermassive black holes (SMBHs) in the early universe, and pose challenges on current black hole formation models to explain the rapid formation of billion solar mass black holes. We measure the quasars' systemic redshifts, black hole masses, Eddington ratios, emission line luminosities, and star formation rates of their host galaxies. Combined with previous results we estimate the fraction of young objects within the high-redshift quasar population at large to be $5\%\lesssim f_{\rm young}\lesssim 10\%$. One of the young objects, PSO J158-14, shows a very bright dust continuum flux ($F_{\rm cont}=3.46\pm 0.02\,\rm mJy$), indicating a highly star-bursting host galaxy with a star formation rate of approximately $1420\,M_{\odot}\,\rm yr^{-1}$., Accepted for publication in ApJ

We provide a coherent, uniform measurement of the evolution of the logarithmic star formation rate (SFR) - stellar mass ($M_*$) relation, called the main sequence of star-forming galaxies (MS), for galaxies out to $z\sim5$. We measure the MS using mean stacks of 3 GHz radio continuum images to derive average SFRs for $\sim$200,000 mass-selected galaxies at $z>0.3$ in the COSMOS field. We describe the MS relation adopting a new model that incorporates a linear relation at low stellar mass (log($M_*$/M$_\odot$)$, Comment: 48 pages (main paper 25 pages, 15 figures). Accepted for publication in The Astrophysical Journal

We study the space density evolution of active galactic nuclei (AGN) using the 610 MHz radio survey of the XXL-North field, performed with the Giant Metrewave Radio Telescope. The survey covers an area of 30.4 deg2, with a beamsize of 6.5 arcsec. The survey is divided into two parts, one covering an area of 11.9 deg2 with 1σ rms noise of 200 μJy beam−1 and the other spanning 18.5 deg2 with rms noise of 45 μJy beam−1. We extracted the catalog of radio components above 7σ. The catalog was cross-matched with a multi-wavelength catalog of the XXL-North field (covering about 80% of the radio XXL-North field) using a likelihood ratio method, which determines the counterparts based on their positions and their optical properties. The multi-component sources were matched visually with the aid of a computer code: Multi-Catalog Visual Cross-Matching. A flux density cut above 1 mJy selects AGN hosts with a high purity in terms of star formation contamination based on the available source counts. After cross-matching and elimination of observational biases arising from survey incompletenesses, the number of remaining sources was 1150. We constructed the rest-frame 1.4 GHz radio luminosity functions of these sources using the maximum volume method. This survey allows us to probe luminosities of 23 ≲ log(L1.4 GHz[W Hz−1]) ≲ 28 up to redshifts of z ≈ 2.1. Our results are consistent with the results from the literature in which AGN are comprised of two differently evolving populations, where the high luminosity end of the luminosity functions evolves more strongly than the low-luminosity end.

Using the deepest 1.2 mm continuum map to date in the Hubble Ultra Deep Field obtained as part of the ALMA Spectroscopic Survey (ASPECS) large program, we measure the cosmic density of dust and implied gas (H$_{2}+$H I) mass in galaxies as a function of look-back time. We do so by stacking the contribution from all $H$-band selected galaxies above a given stellar mass in distinct redshift bins, $��_{\rm dust}(M_\ast>M,z)$ and $��_{\rm gas}(M_\ast>M,z)$. At all redshifts, $��_{\rm dust}(M_\ast>M,z)$ and $��_{\rm gas}(M_\ast>M,z)$ grow rapidly as $M$ decreases down to $10^{10}\,M_\odot$, but this growth slows down towards lower stellar masses. This flattening implies that at our stellar mass-completeness limits ($10^8\,M_\odot$ and $10^{8.9}\,M_\odot$ at $z\sim0.4$ and $z\sim3$), both quantities converge towards the total cosmic dust and gas mass densities in galaxies. The cosmic dust and gas mass densities increase at early cosmic time, peak around $z\sim2$, and decrease by a factor $\sim4$ and 7, compared to the density of dust and molecular gas in the local universe, respectively. The contribution of quiescent galaxies -- i.e., with little on-going star-formation -- to the cosmic dust and gas mass densities is minor ($\lesssim10\%$). The redshift evolution of the cosmic gas mass density resembles that of the star-formation rate density, as previously found by CO-based measurements. This confirms that galaxies have relatively constant star-formation efficiencies (within a factor $\sim2$) across cosmic time. Our results also imply that by $z\sim0$, a large fraction ($\sim90\%$) of dust formed in galaxies across cosmic time has been destroyed or ejected to the intergalactic medium., 32 pages, 13 figures, Accepted for publication in ApJ

We present a comparison of the interstellar medium traced by [CII] (ALMA), and ionised halo gas traced by Lya (MUSE), in and around QSO host galaxies at z~6. To date, 18 QSOs at this redshift have been studied with both MUSE and high-resolution ALMA imaging; of these, 8 objects display a Lya halo. Using datacubes matched in velocity resolution, we compare and contrast the spatial and kinematic information of the Lya halos and the host galaxies' [CII] (and dust-continuum) emission. We find that the Lya halos extend typically 3-30 times beyond the interstellar medium of the host galaxies. The majority of the Lya halos do not show ordered motion in their velocity fields, whereas most of the [CII] velocity fields do. In those cases where a velocity gradient can be measured in Lya, the kinematics do not align with those derived from the [CII] emission. This implies that the Lya emission is not tracing the outskirts of a large rotating disk that is a simple extension of the central galaxy seen in [CII] emission. It rather suggests that the kinematics of the halo gas are decoupled from those of the central galaxy. Given the scattering nature of Lya, these results need to be confirmed with JWST IFU observations that can constrain the halo kinematics further using the non-resonant Ha line., Published in ApJ

We report new Northern Extended Millimeter Array (NOEMA) observations of the [CII], [NII] and [OI] atomic fine structure lines and dust continuum emission of J1148+5251, a z=6.42 quasar, that probe the physical properties of its interstellar medium (ISM). The radially-averaged [CII] and dust continuum emission have similar extensions (up to $\theta = 2.51^{+0.46}_{-0.25}\ \rm{arcsec}$, corresponding to $r= 9.8^{+3.3}_{-2.1}\ \rm{kpc}$ accounting for beam-convolution), confirming that J1148+5251 is the quasar with the largest [CII]-emitting has reservoir known at these epochs.Moreover, if the [CII] emission is examined only along its NE-SW axis, a significant excess ($>5.8\sigma$) of [CII] emission (with respect to the dust) is detected. The new wide--bandwidth observations enable us to accurately constrain the continuum emission, and do not statistically require the presence of broad [CII] line wings that were reported in previous studies. We also report the first detection of the [OI] and (tentatively) [NII] emission lines in J1148+5251. Using Fine Structure Lines (FSL) ratios of the [CII], [NII], [OI] and previously measured [CI] emission lines, we show that J1148+5251 has similar ISM conditions compared to lower--redshift (ultra)-luminous infrared galaxies. CLOUDY modelling of the FSL ratios exclude X--ray dominated regions (XDR) and favours photodissociation regions (PDR) as the origin of the FSL emission. We find that a high radiation field ($10^{3.5-4.5}\,G_0$), high gas density ($n \simeq 10^{3.5-4.5}\, \rm{cm}^{-3}$) and HI column density of $10^{23} \,\rm{cm^{-2}}$ reproduce the observed FSL ratios well., Comment: Accepted for publication in The Astrophysical Journal

23 pages, 11 figures, 5 tables. Accepted for publication in the Astrophysical Journal; International audience; We present a CO and atomic fine-structure line luminosity function analysis using the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). ASPECS consists of two spatially-overlapping mosaics that cover the entire ALMA 3mm and 1.2mm bands. We combine the results of a line candidate search of the 1.2mm data cube with those previously obtained from the 3mm cube. Our analysis shows that $\sim$80% of the line flux observed at 3mm arises from CO(2-1) or CO(3-2) emitters at $z$=1-3 (`cosmic noon'). At 1.2mm, more than half of the line flux arises from intermediate-J CO transitions ($J_{\rm up}$=3-6); $\sim12$% from neutral carbon lines; and $< 1$% from singly-ionized carbon, [CII]. This implies that future [CII] intensity mapping surveys in the epoch of reionization will need to account for a highly significant CO foreground. The CO luminosity functions probed at 1.2mm show a decrease in the number density at a given line luminosity (in units of $L'$) at increasing $J_{\rm up}$ and redshift. Comparisons between the CO luminosity functions for different CO transitions at a fixed redshift reveal sub-thermal conditions on average in galaxies up to $z\sim 4$. In addition, the comparison of the CO luminosity functions for the same transition at different redshifts reveals that the evolution is not driven by excitation. The cosmic density of molecular gas in galaxies, $\rho_{\rm H2}$, shows a redshift evolution with an increase from high redshift up to $z\sim1.5$ followed by a factor $\sim 6$ drop down to the present day. This is in qualitative agreement with the evolution of the cosmic star-formation rate density, suggesting that the molecular gas depletion time is approximately constant with redshift, after averaging over the star-forming galaxy population.

We present a study of the [CII] 158micron line and underlying far-infrared (FIR) continuum emission of 27 quasar host galaxies at z~6, traced by the Atacama Large Millimeter/submillimeter Array at a spatial resolution of ~1 physical kpc. The [CII] emission in the bright, central regions of the quasars have sizes of 1.0-4.8kpc. The dust continuum emission is typically more compact than [CII]. We find that 13/27 quasars (approximately one-half) have companion galaxies in the field, at projected separations of 3-90kpc. The position of dust emission and the Gaia-corrected positions of the central accreting black holes are cospatial (typical offsets 6 quasar hosts reside. Some outliers with offsets of ~500pc can be linked to disturbed morphologies, most likely due to ongoing or recent mergers. We find no correlation between the central brightness of the FIR emission and the bolometric luminosity of the accreting black hole. The FIR-derived star-formation rate densities (SFRDs) in the host galaxies peak at the galaxies' centers, at typical values between 100 and 1000 M_sun/yr/kpc^2. These values are below the Eddington limit for star formation, but similar to those found in local ultraluminous infrared galaxies. The SFRDs drop toward larger radii by an order of magnitude. Likewise, the [CII]/FIR luminosity ratios of the quasar hosts are lowest in their centers (few x10^-4) and increase by a factor of a few toward the galaxies' outskirts, consistent with resolved studies of lower-redshift sources., 24 pages, 11 figures. Published in ApJ

We present MUSE observations of the gas-rich major-merger BR1202-0725 at z~4.7, which constitutes one of the most overdense fields known in the early Universe. We utilise these data in conjunction with existing ALMA observations to compare and contrast the spatially resolved ionised- and cool-gas content of this system which hosts a quasar (QSO), a sub-millimeter galaxy (SMG), the two known optical companions ("LAE1", "LAE2"), and an additional companion discovered in this work "LAE3" just 5 arcsec to the North of the QSO. We find that QSO BR1202-0725 exhibits a large Ly$\alpha$ halo, covering $\approx55$ pkpc on-sky at surface brightness levels of SB$\geq$1E-17 erg/s/cm$^2$/arcsec$^2$. In contrast, the SMG, of similar far-infrared luminosity and star formation rate (SFR), does not exhibit such a Ly$\alpha$ halo. The QSO's halo exhibits high velocity widths ($\sim1000$ km/s) but the gas motion is to some extent kinematically coupled with the previously observed [CII] bridge between the QSO and the SMG. We note that the object known in the literature as LAE2 shows no local peak of Ly$\alpha$ emission, rather, its profile is more consistent with being part of the QSO's extended Ly$\alpha$ halo. The properties of LAE3 are typical of high-redshift LAEs; we measure F$_{\rm{Ly\alpha}}$(LAE3) = $0.24\pm$0.03E-16 erg/s/cm$^2$, corresponding to SFR$_{\rm{Ly\alpha}}\approx\ $5.0$\pm$0.5 M${_\odot}$/yr. The velocity width is $\Delta v$(LAE3) $\approx 400$ km/s, and equivalent width EW$_0$(Ly$\alpha_{\,5\sigma}^{\,lim})\geq 34.05$ $\\A$, consistent with star formation being the primary driver of Ly$\alpha$ emission. We also note a coherent absorption feature at $\sim -400$km/s in spectra from at least three objects; the QSO, LAE1 and "LAE2" which could imply the presence of an expanding neutral gas shell with an extent of at least $24$ pkpc., Comment: 23 pages, accepted for publication in ApJ

International audience; We compare the extent of the dust, molecular gas and stars in three star-forming galaxies, at $z= 1.4, 1.6$ and $2.7$, selected from the Hubble Ultra Deep Field based on their bright CO and dust-continuum emission as well as their large rest-frame optical sizes. The galaxies have high stellar masses, $\mathrm{M}_*>10^{11}\mathrm{M}_\odot$, and reside on, or slightly below, the main sequence of star-forming galaxies at their respective redshifts. We probe the dust and molecular gas using subarcsecond Atacama Large Millimeter/submillimeter Array observations of the 1.3 mm continuum and CO line emission, respectively, and probe the stellar distribution using \emph{Hubble Space Telescope} observations at 1.6 \textmu m. We find that for all three galaxies the CO emission appears $\gtrsim 30\%$ more compact than the stellar emission. For the $z= 1.4$ and $2.7$ galaxies, the dust emission is also more compact, by $\gtrsim 50\%$, than the stellar emission, whereas for the $z=1.6$ galaxy, the dust and stellar emission have similar spatial extents. This similar spatial extent is consistent with observations of local disk galaxies. However, most high redshift observations show more compact dust emission, likely due to the ubiquity of central starbursts at high redshift and the limited sensitivity of many of these observations. Using the CO emission line, we also investigate the kinematics of the cold interstellar medium in the galaxies, and find that all three have kinematics consistent with a rotation-dominated disk.

We present kpc-scale ALMA and HST imaging of the quasar PJ308-21 at $z$=$6.2342$, tracing dust, gas (via the [CII] 158 $��$m line) and young stars. At a resolution of $\sim0.3''$ ($\approx1.7$ kpc), the system is resolved over $>4''$ ($>$20 kpc). In particular, it features a main component, identified to be the quasar host galaxy, centered on the accreting supermassive black hole; and two other extended components on the West and East side, one redshifted and the other blueshifted relative to the quasar. The [CII] emission of the entire system stretches over $>$1500 km/s along the line of sight. All the components of the system are observed in dust, [CII], and rest-frame UV emission. The inferred [CII] luminosities [(0.9-4.6)$\times 10^9$ L$_\odot$], dust luminosities [(0.15-2.6)$\times10^{12}$ L$_\odot$], and rest-frame UV luminosities [(6.6-15)$\times10^{10}$ L$_\odot$], their ratios, and the implied gas/dust masses and star formation rates [11-290 M$_\odot$ yr$^{-1}$] are typical of high-redshift star-forming galaxies. A toy model of a single satellite galaxy that is tidally stripped by the interaction with the quasar host galaxy can account for the observed velocity and spatial extent of the two extended components. An outflow interpretation of the unique features in PJ308-21 is not supported by the data. PJ308-21 is thus one of the earliest galaxy mergers imaged at cosmic dawn., 9 pages, 6 figures. Accepted for publication in ApJ

We study the interstellar medium in a sample of 27 high-redshift quasar host galaxies at z>6, using the [CII] 158um emission line and the underlying dust continuum observed at ~1kpc resolution with ALMA. By performing uv-plane spectral stacking of both the high and low spatial resolution data, we investigate the spatial and velocity extent of gas, and the size of the dust-emitting regions. We find that the average surface brightness profile of both the [CII] and the dust continuum emission can be described by a steep component within a radius of 2kpc, and a shallower component with a scale length of 2kpc, detected up to ~10kpc. The surface brightness of the extended emission drops below ~1% of the peak at radius of ~5kpc, beyond which it constitutes 10-20% of the total measured flux density. Although the central component of the dust continuum emission is more compact than that of the [CII] emission, the extended components have equivalent profiles. The observed extended components are consistent with those predicted by hydrodynamical simulations of galaxies with similar infrared luminosities, where the dust emission is powered by star formation. The [CII] spectrum measured in the mean uv-plane stacked data can be described by a single Gaussian, with no observable [CII] broad-line emission (velocities in excess of >500km/s), that would be indicative of outflows. Our findings suggest that we are probing the interstellar medium and associated star formation in the quasar host galaxies up to radii of 10kpc, whereas we find no evidence for halos or outflows., Comment: 32 pages, 20 figures, 1 table, accepted for publication in ApJ

We present ALMA 0.″28 × 0.″20 (1.4 kpc × 1.0 kpc) resolution observations of the [C II] 158 μm line and underlying dust continuum in the host galaxy of the most distant quasar currently known, ULAS J1342+0928 at z = 7.5413. Both the [C II] and continuum are detected and spatially resolved. The [C II] line is ∼1.5 times more extended than the continuum emission, showing an elongated and complex structure of approximately 3.2 kpc × 6.4 kpc. Two separate peaks are clearly seen (≳6σ each) in three 100 km s-1 width [C II] channel maps. The velocity field of the [C II] gas does not show evidence of a coherent rotation field but rather chaotic motion reminiscent of an ongoing merger.

We present IRAM/NOEMA, JCMT/SCUBA-2 and VLA observations of the most distant known gravitationally lensed quasar J0439+1634 at z = 6.5. We detect strong dust emission, [CII] 158 $\mu$m, [CI] 369 $\mu$m, [OI] 146 $\mu$m, CO(6-5), CO(7-6), CO(9-8), CO(10-9), H$_{\rm 2}$O $3_{\rm 1,2}-2_{\rm 2,1}$, and H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ lines as well as a weak radio continuum. The strong [CII] line yields a systemic redshift of the host galaxy to be z=6.5188$\pm$0.0001. The magnification makes J0439+1634 the far-infrared (FIR) brightest quasar at z > 6 known, with the brightest [CII] line yet detected at this redshift. The FIR luminosity is (3.4$\pm$0.2)$\times$10$^{13}$ $\mu^{-1}$ $L_{\odot}$, where $\mu$ $\sim$ 2.6 - 6.6 is the magnification of the host galaxy, estimated based on the lensing configuration from HST imaging. We estimate the dust mass to be (2.2$\pm$0.1)$\times$10$^{9}$ $\mu^{-1}$ $M_{\odot}$. The CO Spectral Line Energy Distribution using four CO lines is best fit by a two-component model of the molecular gas excitation. The estimates of molecular gas mass derived from CO lines and atomic carbon mass are consistent, in the range of 3.9 - 8.9 $\times$10$^{10} \mu^{-1}$ $M_{\odot}$. The [CII]/[CI], [CII]/CO, and [OI]/[CII] line luminosity ratios suggest a photodissociation region model with more than one component. The ratio of H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ line luminosity to $L_{\rm TIR}$ is consistent with values in local and high redshift ultra-/hyper-luminous infrared galaxies. The VLA observations reveal an unresolved radio continuum source, and indicate that J0439+1634 is a radio quiet quasar with R = 0.05 - 0.17., Comment: 13 pages, 6 figures. Accepted by ApJ

We report high spatial resolution (~0.076", 410pc) Atacama Large Millimeter/submillimeter Array imaging of the dust continuum and the ionized carbon line [CII] in a luminous quasar host galaxy at z=6.6, 800 million years after the big bang. Based on previous studies, this galaxy hosts a ~1x10^9 M_sun black hole and has a star-formation rate of ~1500 M_sun/yr. The unprecedented high resolution of the observations reveals a complex morphology of gas within 3kpc of the accreting central black hole. The gas has a high velocity dispersion with little ordered motion along the line of sight, as would be expected from gas accretion that has yet to settle in a disk. In addition, we find the presence of [CII] cavities in the gas distribution (with diameters of ~0.5kpc), offset from the central black hole. This unique distribution and kinematics cannot be explained by a simple model. Plausible scenarios are that the gas is located in a truncated or warped disk, or the holes are created by interactions with nearby galaxies or due to energy injection into the gas. In the latter case, the energy required to form the cavities must originate from the central active galactic nucleus, as the required energy far exceeds the energy output expected from supernovae. This energy input into the gas, however, does not inhibit the high rate of star-formation. Both star-formation and black hole activity could have been triggered by interactions with satellite galaxies; our data reveal three additional companions detected in [CII] emission around the quasar., Comment: Published in ApJ Letters

We use ALMA observations of the host galaxy of the quasar ULAS-J1342+0928 at z=7.54 to study the dust continuum and far infrared lines emitted from its interstellar medium. The Rayleigh-Jeans tail of the dust continuum is well sampled with eight different spectral setups, and from a modified black body fit we obtain an emissivity coefficient of beta=1.85+-0.3. Assuming a standard dust temperature of 47 K we derive a dust mass of Mdust=0.35x10^8 Msol and a star formation rate of 150+-30 Msol/yr. We have >4sigma detections of the [CII]_158, [OIII]_88 and [NII]_122 atomic fine structure lines and limits on the [CI]_369, [OI]_146 and [NII]_205 emission. We also report multiple limits of CO rotational lines with Jup >= 7, as well as a tentative 3.3sigma detection of the stack of four CO lines (Jup=11, 10, 8, 7). We find line deficits that are in agreement with local ultra luminous infrared galaxies. Comparison of the [NII]_122 and [CII]_158 lines indicates that the [CII]_158 emission arises predominantly from the neutral medium, and we estimate that the photo-disassociation regions in J1342+0928 have densities 180 cm^-3. Our observations favor a low gas-to-dust ratio of, Comment: 18 pages, 11 figures, 3 tables, accepted for publication in ApJ

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 ��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], A&A accepted; 32 pages, 14 figures, 5 tables, 3 Appendices

We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift evolution of $\rho$(H$_2$). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at $z>1$. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at $z$=1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm SFR}$. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history)., Comment: 22 pages, 11 figures. Paper re-submitted to ApJ after addressing the first round of comments by the referee

We present very long baseline interferometry (VLBI) observations of 179 radio sources in the COSMOS field with extremely high sensitivity using the Green Bank Telescope (GBT) together with the Very Long Baseline Array (VLBA) (VLBA+GBT) at 1.4 GHz, to explore the faint radio population in the flux density regime of tens of $��$Jy. Here, the identification of active galactic nuclei (AGN) is based on the VLBI detection of the source, i.e., it is independent of X-ray or infrared properties. The milli-arcsecond resolution provided by the VLBI technique implies that the detected sources must be compact and have large brightness temperatures, and therefore they are most likely AGN (when the host galaxy is located at z$\geq$0.1). On the other hand, this technique allows us to only positively identify when a radio-active AGN is present, i.e., we cannot affirm that there is no AGN when the source is not detected. For this reason, the number of identified AGN using VLBI should be always treated as a lower limit. We present a catalogue containing the 35 radio sources detected with the VLBA+GBT, 10 of which were not previously detected using only the VLBA. We have constructed the radio source counts at 1.4 GHz using the samples of the VLBA and VLBA+GBT detected sources of the COSMOS field to determine a lower limit for the AGN contribution to the faint radio source population. We found an AGN contribution of >40-75% at flux density levels between 150 $��$Jy and 1 mJy. This flux density range is characterised by the upturn of the Euclidean-normalised radio source counts, which implies a contribution of a new population. This result supports the idea that the sub-mJy radio population is composed of a significant fraction of radio-emitting AGN, rather than solely by star-forming galaxies, in agreement with previous studies., 14 pages, 8 figures, 4 tables. Accepted for publication in A&A

Using spectroscopic observations taken for the VIMOS Ultra-Deep Survey (VUDS) we report here on the discovery of PCl J1001+0220, a massive proto-cluster located at $z_{spec}\sim4.57$ in the COSMOS field. The proto-cluster was initially detected as a $\sim12\sigma$ overdensity of typical star-forming galaxies in the blind spectroscopic survey of the early universe ($2, Comment: 27 pages, 14 figures, 1 table, accepted to A&A, minor changes from v1

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., Comment: Resubmitted to MNRAS after addressing reviewer feedback

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., Comment: Accepted for publication in A&A, 12 pages, 7 figures, 3 tables

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.50$)., Comment: 20 pages, 14 figures

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}, Comment: 11 pages, 13 figures, 1 table; to appear in A&A (XXL special issue)

We present ALMA band 8 observations of the [OIII] 88um line and the underlying thermal infrared continuum emission in the z=6.08 quasar CFHQS J2100-1715 and its dust-obscured starburst companion galaxy (projected distance: ~60 kpc). Each galaxy hosts dust-obscured star formation at rates > 100 M_sun/yr, but only the quasar shows evidence for an accreting 10^9 M_sun black hole. Therefore we can compare the properties of the interstellar medium in distinct galactic environments in two physically associated objects, ~1 Gyr after the Big Bang. Bright [OIII] 88um emission from ionized gas is detected in both systems; the positions and line-widths are consistent with earlier [CII] measurements, indicating that both lines trace the same gravitational potential on galactic scales. The [OIII] 88um/FIR luminosity ratios in both sources fall in the upper range observed in local luminous infrared galaxies of similar dust temperature, although the ratio of the quasar is smaller than in the companion. This suggests that gas ionization by the quasar (expected to lead to strong optical [OIII] 5008A emission) does not dominantly determine the quasar's FIR [OIII] 88um luminosity. Both the inferred number of photons needed for the creation of O++ and the typical line ratios can be accounted for without invoking extreme (top-heavy) stellar initial mass functions in the starbursts of both sources., ApJ Letters, accepted

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~$��$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~$��$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., 49 pages (+appendix), 48 figures, 3 tables, accepted for publication in ApJ

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$., Comment: 4 pages, 1 figure, Proceedings of Science, IAU Symposium 333: "Peering towards Cosmic Dawn"

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 z2.5. Samples of hundreds such high-redshift clusters could potentially constrain cosmological parameters and test cluster and galaxy formation models., ApJ Letters in press

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., 19 pages, 18 figures, 5 tables, accepted for publication in A&A

(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 z75%) dominated by star-forming galaxies., 21 pages, 22 figures, accepted for publication in A&A

We carried out targeted ALMA observations of 129 fields in the COSMOS region at 1.25 mm, detecting 152 galaxies at S/N$\geq$5 with 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$\geq$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\sim$2.0 with a high redshift tail skewing the median redshift to $\tilde{z}$=2.48$\pm$0.05. We find that brighter millimeter sources are preferentially found at higher redshifts. Our faintest sources, with S$_{1.25 \rm mm}$1.8 mJy, have a median redshift of $\tilde{z}$=3.08$\pm$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., Accepted for publication in Astronomy & Astrophysics for review 17 pages in body + 26 in appendices

We determine the radio size distribution of a large sample of 152 SMGs in COSMOS that were detected with ALMA at 1.3 mm. For this purpose, we used the observations taken by the VLA-COSMOS 3 GHz Large Project. One hundred and fifteen of the 152 target SMGs were found to have a 3 GHz counterpart. The median value of the major axis FWHM at 3 GHz is derived to be $4.6\pm0.4$ kpc. The radio sizes 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 $\alpha=-0.67$ and $T_{\rm B}=12.6\pm2$ K. Three of the target SMGs, which are also detected with the VLBA, show clearly higher brightness temperatures than the typical values. Although the observed radio emission appears to be predominantly powered by star formation and supernova activity, our results provide a strong indication of the presence of an AGN in the VLBA and X-ray-detected SMG AzTEC/C61. The median radio-emitting size we have derived is 1.5-3 times larger than the typical FIR dust-emitting sizes of SMGs, but similar to that of the SMGs' molecular gas component traced through mid-$J$ line emission of CO. 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 IR-radio correlation., Comment: 32 pages (incl. 5 appendices), 17 figures, 7 tables; accepted for publication in A&A; abstract abridged for arXiv

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 today's giant ellipticals., 26 pages (incl. 4 appendices), 13 figures, 8 tables; accepted for publication in A&A; abstract abridged for arXiv

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., Comment: 13 pages, 9 figures, 2 tables, to appear in A&A

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., Comment: 17 pages, accepted in A&A

We present observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward a sub-field of the XXL-North 25 deg2 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 × 0.̊1 on the sky. The aim of this paper is twofold. First, we present the 3 GHz VLA radio continuum observations, the final radio mosaic and radio source catalogue, and, second, we perform a detailed analysis of the supercluster in the optical and radio regimes using photometric redshifts from the CFHTLS survey and our new VLA-XXL data. Our final 3 GHz radio mosaic has a resolution of 3.̋2 × 1.̋9, and encompasses an area of 41' × 41' with rms noise level lower than ~ 20 μJy beam-1. The noise in the central 15' × 15' region is ≈ 11 μJy beam-1. From the mosaic we extract a catalogue of 155 radio sources with signal-to-noise ratio (S/N) ≥ 6, eight of which are large, multicomponent sources, and 123 (79%) of which 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 using photometric redshifts from the CFHTLS survey we identify a total of seventeen overdensities at zphot = 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 the brightest cluster galaxies of seven X-ray detected clusters. The full VTA-identified structure extends over ~ 0.̊6 × 0.̊2 on the sky, which corresponds to a physical size of ~ 12 × 4 Mpc2 at z = 0.43. No large radio galaxies are present within the overdensities, and we associate eight (S/N> 7) radio sources with potential group/cluster member galaxies. The spatial distribution of the red and blue VTA-identified potential group member galaxies, selected by their observed g - r colours, suggests that the clusters are not virialised yet, but are dynamically young, as expected for hierarchical structure growth in a ΛCDM universe. Further spectroscopic data are required to analyse the dynamical state of the groups. The full catalogue is available as a queryable database table XXL_VLA_15 via the XXL Master Catalogue browser http://cosmosdb.iasf-milano.inaf.it/XXL. A copy of the catalogue and the mosaic are also available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (http://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?IX/49A FITS file of the reduced image is only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (http://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/592/A8

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.3100M0/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, Comment: 21 pages, 15 figures

Radio emission at cm wavelengths from highly star-forming galaxies, such as SMGs, is dominated by synchrotron radiation arising from supernova activity. Using deep, high-resolution ($1\sigma=2.3$ $\mu$Jy beam$^{-1}$; $0.75^{"}$) cm radio-continuum observations taken by the VLA-COSMOS 3 GHz Large Project, we studied the radio-emitting sizes of a flux-limited sample of SMGs in the COSMOS field. Of the 39 SMGs studied here, 3 GHz emission was detected towards 18 of them ($\sim46\pm11\%$) with S/N ratios in the range of ${\rm S/N=4.2-37.4}$. Using 2D elliptical Gaussian fits, we derived a median deconvolved major axis FWHM size of $0.54^{"}\pm 0.11^{"}$ for our 18 SMGs detected at 3 GHz. For the 15 SMGs with known redshift we derived a median linear major axis FWHM of $4.2\pm0.9$ kpc. No clear correlation was found between the radio-emitting size and the 3 GHz or submm flux density, or the redshift of the SMG. However, there is a hint of larger radio sizes at $z\sim2.5-5$ compared to lower redshifts. The sizes we derived are consistent with previous SMG sizes measured at 1.4 GHz and in mid-$J$ CO emission, but significantly larger than those seen in the (sub)mm continuum emission. One possible scenario is that SMGs have i) an extended gas component with a low dust temperature, and which can be traced by low- to mid-$J$ CO line emission and radio continuum emission, and ii) a warmer, compact starburst region giving rise to the high-$J$ line emission of CO, which could dominate the dust continuum size measurements. Because of the rapid cooling of CR electrons in dense starburst galaxies ($\sim10^4-10^5$ yr), the more extended synchrotron radio-emitting size being a result of CR diffusion seems unlikely. Instead, if SMGs are driven by galaxy mergers the radio synchrotron emission might arise from more extended magnetised ISM around the starburst region., Comment: 18 pages (incl. 5 appendices), 7 figures, and 5 tables; accepted for publication in A&A; abstract abridged for arXiv

We used the PdBI to map a sample of 15 SMGs in the COSMOS field at the wavelength of 1.3 mm. The target SMGs were originally discovered in the JCMT/AzTEC 1.1 mm continuum survey at S/N=4-4.5. This paper presents, for the first time, interferometric millimetre-wavelength observations of these sources. The angular resolution of our observations, 1.8", allowed us to accurately determine the positions of the target SMGs. Using a detection threshold of S/N>4.5 regardless of multiwavelength counterpart association, and 4, Comment: 31 pages (incl. 5 appendices with 3 figures and 1 table), 6 figures and 4 tables in the main text. Accepted for publication in A&A

We present deep 3 GHz VLA observations of the potentially recoiling black hole CID-42 in the COSMOS field. This galaxy shows two optical nuclei in the HST/ACS image and a large velocity offset of ~ 1300 km/s between the broad and narrow H beta emission line although the spectrum is not spacially resolved (Civano et al. 2010). The new 3 GHz VLA data has a bandwidth of 2 GHz and to correctly interpret the flux densities imaging was done with two different methods: multi-scale multi-frequency synthesis and spectral windows stacking. The final resolutions and sensitivities of these maps are 0.7" with rms = 4.6 muJy/beam and 0.9" with rms = 4.8 muJy/beam respectively. With a 7 sigma 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, VLBA and GMRT data, to construct a radio synchrotron spectrum of CID-42. The radio spectrum suggests a type I unobscured radio-quiet flat-spectrum 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. 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., Accepted for publication in MNRAS; 7 pages, 4 figures, 2 tables