Showing 1–50 of 69 results for author: Novak, M

A super-linear "radio-AGN main sequence'' links mean radio-AGN power and galaxy stellar mass since z$\sim$3

Authors: I. Delvecchio, E. Daddi, M. T. Sargent, J. Aird, J. R. Mullaney, B. Magnelli, D. Elbaz, L. Bisigello, L. Ceraj, S. Jin, B. S. Kalita, D. Liu, M. Novak, I. Prandoni, J. F. Radcliffe, C. Spingola, G. Zamorani, V. Allevato, G. Rodighiero, V. Smolcic

Abstract: Mapping the average AGN luminosity across galaxy populations and over time encapsulates important clues on the interplay between supermassive black hole (SMBH) and galaxy growth. This paper presents the demography, mean power and cosmic evolution of radio AGN across star-forming galaxies (SFGs) of different stellar masses (${M_{*}}$). We exploit deep VLA-COSMOS 3 GHz data to build the rest-frame 1… ▽ More Mapping the average AGN luminosity across galaxy populations and over time encapsulates important clues on the interplay between supermassive black hole (SMBH) and galaxy growth. This paper presents the demography, mean power and cosmic evolution of radio AGN across star-forming galaxies (SFGs) of different stellar masses (${M_{*}}$). We exploit deep VLA-COSMOS 3 GHz data to build the rest-frame 1.4 GHz AGN luminosity functions at 0.1$\leq$$z$$\leq$4.5 hosted in SFGs. Splitting the AGN luminosity function into different ${M_{*}}$ bins reveals that, at all redshifts, radio AGN are both more frequent and more luminous in higher ${M_*}$ than in lower ${M_*}$ galaxies. The cumulative kinetic luminosity density exerted by radio AGN in SFGs peaks at $z$$\sim$2, and it is mostly driven by galaxies with 10.5$\leq$$\log$(${M_{*}}$/${M_{\odot}}$)$<$11. Averaging the cumulative radio AGN activity across all SFGs at each (${M_{*}}$,$z$) results in a "radio-AGN main sequence" that links the time-averaged radio-AGN power $\langle$$L_{1.4}^{AGN}$$\rangle$ and galaxy stellar mass, in the form: $\log$$\langle$[$L_{1.4}^{AGN}$/ W Hz$^{-1}]\rangle$ = (20.97$\pm$0.16) + (2.51$\pm$0.34)$\cdot$$\log$(1+$z$) + (1.41$\pm$0.09)$\cdot$($\log$[${M_{*}}$/${M_{\odot}}$] -10). The super-linear dependence on ${M_{*}}$, at fixed redshift, suggests enhanced radio-AGN activity in more massive SFGs, as compared to star formation. We ascribe this enhancement to both a higher radio AGN duty cycle and a brighter radio-AGN phase in more massive SFGs. A remarkably consistent ${M_{*}}$ dependence is seen for the evolving X-ray AGN population in SFGs. This similarity is interpreted as possibly driven by secular cold gas accretion fueling both radio and X-ray AGN activity in a similar fashion over the galaxy's lifetime. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: Accepted for publication in A&A. 20 pages + Appendices

Journal ref: A&A 668, A81 (2022)

The Decoupled Kinematics of high-z QSO Host Galaxies and their Lya halos

Authors: Alyssa Drake, Marcel Neeleman, Bram P. Venemans, Mladen Novak, Fabian Walter, Eduardo Banados, Roberto Decarli, Emanuele Paolo Farina, Chiara Mazzucchelli, Maxime Trebitsch

Abstract: 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 inf… ▽ More 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. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: Published in ApJ

Journal ref: 2022ApJ...929...86D

The evolution of the radio luminosity function of group galaxies in COSMOS

Authors: E. Vardoulaki, G. Gozaliasl, A. Finoguenov, M. Novak, H. G. Khosroshahi

Abstract: [Abridged] To understand the role of the galaxy group environment on galaxy evolution, we present a study of radio luminosity functions (RLFs) of group galaxies based on the VLA-COSMOS 3 GHz Large Project. The radio-selected sample of 7826 COSMOS galaxies with robust optical/near-infrared counterparts, excellent photometric coverage, and the COSMOS X-ray galaxy groups (M_200c > 10^13.5 M_sun) enab… ▽ More [Abridged] To understand the role of the galaxy group environment on galaxy evolution, we present a study of radio luminosity functions (RLFs) of group galaxies based on the VLA-COSMOS 3 GHz Large Project. The radio-selected sample of 7826 COSMOS galaxies with robust optical/near-infrared counterparts, excellent photometric coverage, and the COSMOS X-ray galaxy groups (M_200c > 10^13.5 M_sun) enables us to construct the RLF of group galaxies (GGs) and their contribution to the total RLF since z ~ 2.3. Using the Markov chain Monte Carlo algorithm, we fit a redshift-dependent pure luminosity evolution model and a linear and power-law model to the luminosity functions. We compare it with past RLF studies from VLA-COSMOS on individual populations of radio-selected star-forming galaxies (SFGs) and galaxies hosting active galactic nuclei (AGN). These populations are classified based on the presence or absence of a radio excess concerning the star-formation rates derived from the infrared emission. We find that the density of radio galaxies in groups is low compared to the field at z ~ 2 down to z ~ 1.25, followed by a sharp increase at z ~ 1 by a factor of 6, and then a smooth decline towards low redshifts. This trend is caused by both decrease in the volume abundance of massive groups at high-z and the changes in the halo occupation of radio AGN, which are found by other studies to reside at smaller halo mass groups. This indicates that the bulk of high-$z$ log10(M_200c/M_sun) > 13.5 groups must have been forming recently, and so the cooling has not been established as yet. The slope of the GG RLF is flatter compared to the field, with excess at high radio luminosities. The evolution in the GG RLF is driven mainly by satellite galaxies in groups. A drop in occurrence of AGN in groups at z > 1 by a factor of 6, manifests an important detail on the processes governing galaxy evolution. △ Less

Submitted 16 May, 2024; v1 submitted 5 April, 2022; originally announced April 2022.

Comments: A&A in review; 30 pages, 8 figures, 12 tables

Physical constraints on the extended interstellar medium of the z=6.42 quasar J1148+5251: [CII] 158um, [NII] 205um and [OI] 146um observations

Authors: Romain A. Meyer, Fabian Walter, Claudia Cicone, Pierre Cox, Roberto Decarli, Roberto Neri, Mladen Novak, Antonio Pensabene, Dominik Riechers, Axel Weiss

Abstract: 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 $θ= 2.51^{+0.46}_{-0.25}\ \rm{arcsec}$, corresponding to… ▽ More 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 $θ= 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σ$) 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. △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: Accepted for publication in The Astrophysical Journal

The Impact of Powerful Jets on the Far-infrared Emission of an Extreme Radio Quasar at z~6

Authors: Sofía Rojas-Ruiz, Eduardo Bañados, Marcel Neeleman, Thomas Connor, Anna-Christina Eilers, Bram P. Venemans, Yana Khusanova, Christopher L. Carilli, Chiara Mazzucchelli, Roberto Decarli, Emmanuel Momjian, Mladen Novak

Abstract: 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… ▽ More 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. △ Less

Submitted 9 August, 2021; originally announced August 2021.

Comments: 14 pages, 5 figures

The XXL Survey: XLIII. The quasar radio loudness dichotomy exposed via radio luminosity functions obtained by combining results from COSMOS and XXL-S X-ray selected quasars

Authors: Lana Ceraj, Vernesa Smolčić, Ivan Delvecchio, Andrew Butler, Krešimir Tisanić, Jacinta Delhaize, Cathy Horellou, Jeyhan Kartaltepe, Konstantinos Kolokythas, Sarah Leslie, Stefano Marchesi, Mladen Novak, Marguerite Pierre, Manolis Plionis, Eleni Vardoulaki, Giovanni Zamorani

Abstract: 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 rad… ▽ More 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 <= z <= 3.7. The lower-1.4 GHz luminosity end shows a higher normalization than expected only from AGN contribution in all studied redshift bins. The found "bump" is mostly dominated by emission due to star-forming (SF) processes within the XQSO host galaxies. As expected, AGN-related radio emission dominates at the higher-luminosity end of RLF. The evolution of XQSO RLF was studied via combination of analytic forms from the literature to constrain the lower-luminosity "bump" and the higher-luminosity AGN part of the RLF. We defined two 1.4 GHz luminosity thresholds, L_th,SF and L_th,AGN, below and above which more than 80% of sources contributing to the RLF are dominated by SF and AGN-related activity, respectively. These thresholds evolve with redshift, most likely due to the strong evolution of SFRs of the XQSO host galaxies. △ Less

Submitted 16 June, 2021; originally announced June 2021.

Comments: 14 pages, 11 figures

Journal ref: A&A 642, A125 (2020)

An ALMA multi-line survey of the ISM in two quasar host-companion galaxy pairs at $z > 6$

Authors: A. Pensabene, R. Decarli, E. Bañados, B. Venemans, F. Walter, F. Bertoldi, X. Fan, E. P. Farina, J. Li, C. Mazzucchelli, M. Novak, D. Riechers, H. -W. Rix, M. A. Strauss, R. Wang, A. Weiß, J. Yang, Y. Yang

Abstract: 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μm}$, [CI]… ▽ More 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μm}$, [CI]$_{\rm 369μ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μ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μ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μm}$ and [CI]$_{\rm 369μ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]. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 26 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 652, A66 (2021)

The Kinematics of z ~ 6 Quasar Host Galaxies

Authors: Marcel Neeleman, Mladen Novak, Bram P. Venemans, Fabian Walter, Roberto Decarli, Melanie Kaasinen, Jan-Torge Schindler, Eduardo Banados, Chris L. Carilli, Alyssa B. Drake, Xiaohui Fan, Hans-Walter Rix

Abstract: 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 gase… ▽ More 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. △ Less

Submitted 10 February, 2021; originally announced February 2021.

Comments: 33 pages. Kinematic fitting code is available at https://github.com/mneeleman/qubefit

No evidence for [CII] halos or high-velocity outflows in z>6 quasar host galaxies

Authors: Mladen Novak, Bram P. Venemans, Fabian Walter, Marcel Neeleman, Melanie Kaasinen, Lichen Liang, Robert Feldmann, Eduardo Banados, Chris Carilli, Roberto Decarli, Alyssa B. Drake, Xiaohui Fan, Emanuele P. Farina, Chiara Mazzucchelli, Hans-Walter Rix, Ran Wang

Abstract: 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… ▽ More 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. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: 32 pages, 20 figures, 1 table, accepted for publication in ApJ

Kiloparsec-scale ALMA Imaging of [CII] and Dust Continuum Emission of 27 Quasar Host Galaxies at z~6

Authors: Bram Venemans, Fabian Walter, Marcel Neeleman, Mladen Novak, Justin Otter, Roberto Decarli, Eduardo Bañados, Alyssa Drake, Emanuele Farina, Melanie Kaasinen, Chiara Mazzucchelli, Chris Carilli, Xiaohui Fan, Hans-Walter Rix, Ran Wang

Abstract: 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 fi… ▽ More 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 <0.1"). This suggests that the central black holes are located at the bottom of the gravitational wells of the dark matter halos in which the z>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. △ Less

Submitted 28 November, 2020; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: 24 pages, 11 figures. Published in ApJ

Journal ref: The Astrophysical Journal, 2020, Volume 904, page 130 (20pp)

A Closer Look at Two of the Most Luminous Quasars in the Universe

Authors: Jan-Torge Schindler, Xiaohui Fan, Mladen Novak, Bram Venemans, Fabian Walter, Feige Wang, Jinyi Yang, Minghao Yue, Eduardo Banados, Yun-Hsin Huang

Abstract: Ultra-luminous quasars ($M_{1450} \leq -29$) provide us with a rare view into the nature of the most massive and most rapidly accreting supermassive black holes (SMBHs). Following the discovery of two of these extreme sources, J0341${+}$1720 ($M_{1450}=-29.56$, $z=3.71$) and J2125${-}$1719 ($M_{1450}=-29.39$, $z=3.90$), in the Extremely Luminous Quasar Survey (ELQS) and its extension to the Pan-ST… ▽ More Ultra-luminous quasars ($M_{1450} \leq -29$) provide us with a rare view into the nature of the most massive and most rapidly accreting supermassive black holes (SMBHs). Following the discovery of two of these extreme sources, J0341${+}$1720 ($M_{1450}=-29.56$, $z=3.71$) and J2125${-}$1719 ($M_{1450}=-29.39$, $z=3.90$), in the Extremely Luminous Quasar Survey (ELQS) and its extension to the Pan-STARRS\,1 footprint (PS-ELQS), we herein present an analysis of their rest-frame UV to optical spectroscopy. Both quasars harbor very massive SMBHs with $M_{\rm{BH}}=6.73_{-0.83}^{+0.75}\times10^{9}\,M_{\odot}$ and $M_{\rm{BH}}=5.45_{-0.55}^{+0.60}\times10^{9}\,M_{\odot}$, respectively, showing evidence of accretion above the Eddington limit ($L_{\rm{bol}}/L_{\rm{Edd}}=2.74_{-0.27}^{+0.39}$ and $L_{\rm{bol}}/L_{\rm{Edd}}=3.01_{-0.30}^{+0.34}$). NOEMA 3 millimeter observations of J0341${+}$1720 reveal a highly star-forming ($\rm{SFR}\approx1500\,M_{\odot}\,\rm{yr}^{-1}$), ultra-luminous infrared galaxy ($L_{\rm{TIR}}\approx1.0\times10^{13}\,L_{\odot}$) host, which, based on an estimate of its dynamical mass, is only ${\sim}30$ times more massive than the SMBH it harbors at its center. As examples of luminous super-Eddington accretion, these two quasars provide support for theories, which explain the existence of billion solar mass SMBHs ${\sim}700$ million years after the Big Bang by moderate super-Eddington growth from standard SMBH seeds. △ Less

Submitted 27 October, 2020; originally announced October 2020.

Comments: Accepted for publication in the Astrophysical Journal (ApJ)

The infrared-radio correlation of star-forming galaxies is strongly M$_{\star}$-dependent but nearly redshift-invariant since z$\sim$4

Authors: I. Delvecchio, E. Daddi, M. T. Sargent, M. J. Jarvis, D. Elbaz, S. Jin, D. Liu, I. H. Whittam, H. Algera, R. Carraro, C. D'Eugenio, J. Delhaize, B. S. Kalita, S. Leslie, D. Cs. Molnar, M. Novak, I. Prandoni, V. Smolcic, Y. Ao, M. Aravena, F. Bournaud, J. D. Collier, S. M. Randriamampandry, Z. Randriamanakoto, G. Rodighiero , et al. (3 additional authors not shown)

Abstract: Several works in the past decade have used the ratio between total (rest 8-1000$μ$m) infrared and radio (rest 1.4~GHz) luminosity in star-forming galaxies (q$_{IR}$), often referred to as the "infrared-radio correlation" (IRRC), to calibrate radio emission as a star formation rate (SFR) indicator. Previous studies constrained the evolution of q$_{IR}$ with redshift, finding a mild but significant… ▽ More Several works in the past decade have used the ratio between total (rest 8-1000$μ$m) infrared and radio (rest 1.4~GHz) luminosity in star-forming galaxies (q$_{IR}$), often referred to as the "infrared-radio correlation" (IRRC), to calibrate radio emission as a star formation rate (SFR) indicator. Previous studies constrained the evolution of q$_{IR}$ with redshift, finding a mild but significant decline, that is yet to be understood. For the first time, we calibrate q$_{IR}$ as a function of \textit{both} stellar mass (M$_{\star}$) and redshift, starting from an M$_{\star}$-selected sample of $>$400,000 star-forming galaxies in the COSMOS field, identified via (NUV-r)/(r-J) colours, at redshifts 0.1$<$z$<$4.5. Within each (M$_{\star}$,z) bin, we stack the deepest available infrared/sub-mm and radio images. We fit the stacked IR spectral energy distributions with typical star-forming galaxy and IR-AGN templates, and carefully remove radio AGN candidates via a recursive approach. We find that the IRRC evolves primarily with M$_{\star}$, with more massive galaxies displaying systematically lower q$_{IR}$. A secondary, weaker dependence on redshift is also observed. The best-fit analytical expression is the following: q$_{IR}$(M$_{\star}$,z)=(2.646$\pm$0.024)$\times$(1+z)$^{(-0.023\pm0.008)}$-(0.148$\pm$0.013)$\times$($\log~M_{\star}$/M$_{\odot}$-10). The lower IR/radio ratios seen in more massive galaxies are well described by their higher observed SFR surface densities. Our findings highlight that using radio-synchrotron emission as a proxy for SFR requires novel M$_{\star}$-dependent recipes, that will enable us to convert detections from future ultra deep radio surveys into accurate SFR measurements down to low-SFR, low-M$_{\star}$ galaxies. △ Less

Submitted 22 January, 2021; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: Accepted for publication in A&A. 22 pages + Appendices. 24 figures, 4 tables

Journal ref: A&A 647, A123 (2021)

A Multi-wavelength Analysis of the Faint Radio Sky (COSMOS-XS): the Nature of the Ultra-faint Radio Population

Authors: H. S. B. Algera, D. Van der Vlugt, J. A. Hodge, I. Smail, M. Novak, J. F. Radcliffe, D. A. Riechers, H. Röttgering, V. Smolčić, F. Walter

Abstract: Ultra-deep radio surveys are an invaluable probe of dust-obscured star formation, but require a clear understanding of the relative contribution from radio AGN to be used to their fullest potential. We study the composition of the $μ$Jy radio population detected in the Karl G. Jansky Very Large Array COSMOS-XS survey based on a sample of 1540 sources detected at 3 GHz over an area of… ▽ More Ultra-deep radio surveys are an invaluable probe of dust-obscured star formation, but require a clear understanding of the relative contribution from radio AGN to be used to their fullest potential. We study the composition of the $μ$Jy radio population detected in the Karl G. Jansky Very Large Array COSMOS-XS survey based on a sample of 1540 sources detected at 3 GHz over an area of $\sim350\text{arcmin}^2$. This ultra-deep survey consists of a single pointing in the well-studied COSMOS field at both 3 and 10 GHz and reaches RMS-sensitivities of $0.53$ and $0.41μ$Jy beam$^{-1}$, respectively. We find multi-wavelength counterparts for $97\%$ of radio sources, based on a combination of near-UV/optical to sub-mm data, and through a stacking analysis at optical/near-infrared wavelengths we further show that the sources lacking such counterparts are likely to be high-redshift in nature (typical $z\sim4-5$). Utilizing the multi-wavelength data over COSMOS, we identify AGN through a variety of diagnostics and find these to make up $23.2\pm1.3\%$ of our sample, with the remainder constituting uncontaminated star-forming galaxies. However, more than half of the AGN exhibit radio emission consistent with originating from star-formation, with only $8.8\pm0.8\%$ of radio sources showing a clear excess in radio luminosity. At flux densities of $\sim30μ$Jy at 3 GHz, the fraction of star-formation powered sources reaches $\sim90\%$, and this fraction is consistent with unity at even lower flux densities. Overall, our findings imply that ultra-deep radio surveys such as COSMOS-XS constitute a highly effective means of obtaining clean samples of star-formation powered radio sources. △ Less

Submitted 28 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJ; 23 pages + appendices; 15 figures in main text. This is a companion paper to Van der Vlugt et al., on arXiv today

An Ultra-deep Multi-band VLA Survey of the Faint Radio Sky (COSMOS-XS): Source Catalog and Number Counts

Authors: D. Van der Vlugt, H. S. B. Algera, J. A. Hodge, M. Novak, J. F. Radcliffe, D. A. Riechers, H. Röttgering, V. Smolčić, F. Walter

Abstract: We present ultra-deep, matched-resolution Karl G. Jansky Very Large Array (VLA) observations at 10 and $3$ GHz in the COSMOS field: the COSMOS-XS survey. The final 10 and $3$ GHz images cover $\sim16\rm{arcmin}^{2}$ and $\sim180\rm{arcmin}^{2}$ and reach median rms values of $0.41μ\rm{Jy\,beam}^{-1}$ and $0.53μ\rm{Jy\,beam}^{-1}$, respectively. Both images have an angular resolution of… ▽ More We present ultra-deep, matched-resolution Karl G. Jansky Very Large Array (VLA) observations at 10 and $3$ GHz in the COSMOS field: the COSMOS-XS survey. The final 10 and $3$ GHz images cover $\sim16\rm{arcmin}^{2}$ and $\sim180\rm{arcmin}^{2}$ and reach median rms values of $0.41μ\rm{Jy\,beam}^{-1}$ and $0.53μ\rm{Jy\,beam}^{-1}$, respectively. Both images have an angular resolution of $\sim 2.0''$. To fully account for the spectral shape and resolution variations across the broad bands, we image all data with a multi-scale, multi-frequency synthesis algorithm. We present source catalogs for the 10 and $3$ GHz image with 91 and 1498 sources, respectively, above a peak brightness threshold of $5σ$. We present source counts with completeness corrections included that are computed via Monte Carlo simulations. Our corrected radio counts at $3$ GHz with direct detections down to $\sim2.8μ$Jy are consistent within the uncertainties with other results at 3 and 1.4 GHz, but extend to fainter flux densities than previous direct detections. The ultra-faint $3$ GHz number counts are found to exceed the counts predicted by the semi-empirical radio sky simulations developed in the framework of the SKA Simulated Skies project, consistent with previous P(D) analyses. Our measured source counts suggest a steeper luminosity function evolution for these faint star-forming sources. The semi-empirical Tiered Radio Extragalactic Continuum Simulation (T-RECS) predicts this steeper evolution and is in better agreement with our results. The $10$ GHz radio number counts also agree with the counts predicted by the T-RECS simulation within the expected variations from cosmic variance. In summary, the multi-band, matched-resolution COSMOS-XS survey in the well-studied COSMOS field provides a high-resolution view of the ultra-faint radio sky that can help guide next generation radio facilities. △ Less

Submitted 28 September, 2020; originally announced September 2020.

Comments: Resubmitted to ApJ after incorporating suggestions from the referee; 20 pages, 17 figures. This is a companion paper to Algera et al., on arXiv today

The ALMA Spectroscopic Survey in the HUDF: Multi-band constraints on line luminosity functions and the cosmic density of molecular gas

Authors: Roberto Decarli, Manuel Aravena, Leindert Boogaard, Chris Carilli, Jorge González-López, Fabian Walter, Paulo C. Cortes, Pierre Cox, Elisabete da Cunha, Emanuele Daddi, Tanio Díaz-Santos, Jacqueline A. Hodge, Hanae Inami, Marcel Neeleman, Mladen Novak, Pascal Oesch, Gergö Popping, Dominik Riechers, Ian Smail, Bade Uzgil, Paul van der Werf, Jeff Wagg, Axel Weiss

Abstract: 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… ▽ More 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, $ρ_{\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. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: 23 pages, 11 figures, 5 tables. Accepted for publication in the Astrophysical Journal

The ALMA Spectroscopic Survey Large Program: The Infrared Excess of z=1.5-10 UV-selected Galaxies and the Implied High-Redshift Star Formation History

Authors: Rychard Bouwens, Jorge Gonzalez-Lopez, Manuel Aravena, Roberto Decarli, Mladen Novak, Mauro Stefanon, Fabian Walter, Leindert Boogaard, Chris Carilli, Ugne Dudzeviciute, Ian Smail, Emanuele Daddi, Elisabete da Cunha, Rob Ivison, Themiya Nanayakkara, Paulo Cortes, Pierre Cox, Hanae Inami, Pascal Oesch, Gergo Popping, Dominik Riechers, Paul van der Werf, Axel Weiss, Yoshi Fudamoto, Jeff Wagg

Abstract: 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 s… ▽ More 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.25} Msolar) galaxies over the ASPECS footprint, we find a mean continuum flux of -0.1+/-0.4 microJy/beam, implying a hard upper limit on the obscured SFR of <0.6 Msolar/yr (4 sigma) in a typical low-mass galaxy. The correlation between the infrared excess IRX of UV-selected galaxies (L(IR)/L(UV)) and the UV-continuum slope is also seen in our ASPECS data and shows consistency with a Calzetti-like relation at >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. △ Less

Submitted 19 October, 2020; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: 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

Probing the Nature of High Redshift Weak Emission Line Quasars: A Young Quasar with a Starburst Host Galaxy

Authors: Irham Taufik Andika, Knud Jahnke, Masafusa Onoue, Eduardo Bañados, Chiara Mazzucchelli, Mladen Novak, Anna-Christina Eilers, Bram P. Venemans, Jan-Torge Schindler, Fabian Walter, Marcel Neeleman, Robert A. Simcoe, Roberto Decarli, Emanuele Paolo Farina, Victor Marian, Antonio Pensabene, Thomas M. Cooper, Alejandra F. Rojas

Abstract: 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… ▽ More 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$α$+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. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 28 pages, 16 figures, 3 tables, accepted for publication in ApJ

Journal ref: The Astrophysical Journal, 903:34 (20pp), 2020 November 1

The Ionised- and Cool-Gas Content of The BR1202-0725 System as seen by MUSE and ALMA

Authors: A. B. Drake, F. Walter, M. Novak, E. P. Farina, M. Neeleman, D. Riechers, C. Carilli, R. Decarli, C. Mazzucchelli, M. Onoue

Abstract: 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 optica… ▽ More 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$α$ 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$α$ 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$α$ emission, rather, its profile is more consistent with being part of the QSO's extended Ly$α$ halo. The properties of LAE3 are typical of high-redshift LAEs; we measure F$_{\rm{Lyα}}$(LAE3) = $0.24\pm$0.03E-16 erg/s/cm$^2$, corresponding to SFR$_{\rm{Lyα}}\approx\ $5.0$\pm$0.5 M${_\odot}$/yr. The velocity width is $Δv$(LAE3) $\approx 400$ km/s, and equivalent width EW$_0$(Ly$α_{\,5σ}^{\,lim})\geq 34.05$ $\\A$, consistent with star formation being the primary driver of Ly$α$ 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. △ Less

Submitted 28 July, 2020; originally announced July 2020.

Comments: 23 pages, accepted for publication in ApJ

A Comparison of the Stellar, CO and Dust-Continuum Emission from Three, Star-Forming HUDF Galaxies at $z\sim 2$

Authors: Melanie Kaasinen, Fabian Walter, Mladen Novak, Marcel Neeleman, Ian Smail, Leindert Boogaard, Elisabete da Cunha, Axel Weiss, Daizhong Liu, Roberto Decarli, Gergö Popping, Tanio Diaz-Santos, Paulo Cortés, Manuel Aravena, Paul van der Werf, Dominik Riechers, Hanae Inami, Jacqueline A. Hodge, Hans-Walter Rix, Pierre Cox

Abstract: 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 sta… ▽ More 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. △ Less

Submitted 8 September, 2020; v1 submitted 7 July, 2020; originally announced July 2020.

The VLA-COSMOS 3 GHz Large Project: Evolution of specific star formation rates out to $z\sim5$

Authors: Sarah Leslie, Eva Schinnerer, Daizhong Liu, Benjamin Magnelli, Hiddo Algera, Alexander Karim, Iary Davidzon, Ghassem Gozaliasl, Eric F. Jiménez-Andrade, Philipp Lang, Mark Sargent, Mladen Novak, Brent Groves, Vernesa Smolčić, Giovanni Zamorani, Mattia Vaccari, Andrew Battisti, Eleni Vardoulaki, Yingjie Peng, Jeyhan Kartaltepe

Abstract: 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… ▽ More 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$)$<$10) and a flattening at high stellar mass that becomes more prominent at low redshift ($z<1.5$). We find that the SFR density peaks at $1.5<z<2$ and at each epoch there is a characteristic stellar mass ($M_* = 1 - 4 \times 10^{10}\mathrm{M}_\odot$) that contributes the most to the overall SFR density. This characteristic mass increases with redshift, at least to $z\sim2.5$. We find no significant evidence for variations in the MS relation for galaxies in different environments traced by the galaxy number density at $0.3<z<3$, nor for galaxies in X-ray groups at $z\sim0.75$. We confirm that massive bulge-dominated galaxies have lower SFRs than disk-dominated galaxies at a fixed stellar mass at $z<1.2$. As a consequence, the increase in bulge-dominated galaxies in the local star-forming population leads to a flattening of the MS at high stellar masses. This indicates that "mass-quenching" is linked with changes in the morphological composition of galaxies at a fixed stellar mass. △ Less

Submitted 24 June, 2020; originally announced June 2020.

Comments: 48 pages (main paper 25 pages, 15 figures). Accepted for publication in The Astrophysical Journal

The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: The nature of the faintest dusty star-forming galaxies

Authors: Manuel Aravena, Leindert Boogaard, Jorge Gónzalez-López, Roberto Decarli, Fabian Walter, Chris L. Carilli, Ian Smail, Axel Weiss, Roberto J. Assef, Franz Erik Bauer, Rychard J. Bouwens, Paulo C. Cortes, Pierre Cox, Elisabete da Cunha, Emanuele Daddi, Tanio Díaz-Santos, Hanae Inami, Rob Ivison, Mladen Novak, Gergö Popping, Dominik Riechers, Paul van der Werf, Jeff Wagg

Abstract: We present a characterization of the physical properties of a sample of 35 securely-detected, dusty galaxies in the deep ALMA 1.2-mm image obtained as part of the ALMA Spectroscopic Survey in the {\it Hubble} Ultra Deep Field (ASPECS) Large Program. This sample is complemented by 26 additional sources identified via an optical/infrared source positional prior. Using their well-characterized spectr… ▽ More We present a characterization of the physical properties of a sample of 35 securely-detected, dusty galaxies in the deep ALMA 1.2-mm image obtained as part of the ALMA Spectroscopic Survey in the {\it Hubble} Ultra Deep Field (ASPECS) Large Program. This sample is complemented by 26 additional sources identified via an optical/infrared source positional prior. Using their well-characterized spectral energy distributions, we derive median stellar masses and star formation rates (SFR) of $4.8\times10^{10}~M_\odot$ and 30 $M_\odot$ yr$^{-1}$, and interquartile ranges of $(2.4-11.7)\times10^{10}~M_\odot$ and $20-50~M_\odot$ yr$^{-1}$, respectively. We derive a median spectroscopic redshift of 1.8 with an interquartile range $1.1-2.6$, significantly lower than submillimeter galaxies detected in shallower, wide-field surveys. We find that 59\%$\pm$13\%, 6\%$\pm$4\%, and 34\%$\pm$9\% of our sources are within, above and below $\pm0.4$ dex from the SFR-stellar mass relation or main-sequence (MS), respectively. The ASPECS galaxies closely follow the SFR-molecular gas mass relation and other previously established scaling relations, confirming a factor of five increase of the gas-to-stellar mass ratio from $z=0.5$ to $z=2.5$ and a mild evolution of the gas depletion timescales with a typical value of 0.7 Gyr at $z=1-3$. ASPECS galaxies located significantly below the MS, a poorly exploited parameter space, have low gas-to-stellar-mass ratios of $\sim0.1-0.2$ and long depletion timescales $>1$ Gyr. Galaxies along the MS dominate the cosmic density of molecular gas at all redshifts. Systems above the main sequence have an increasing contribution to the total gas reservoirs from $z<1$ to $z=2.5$, while the opposite is found for galaxies below the MS. △ Less

Submitted 22 September, 2020; v1 submitted 7 June, 2020; originally announced June 2020.

Comments: ApJ published on Sept. 24, 2020. Contains 36 pages, 17 figures. Final version matching the proof reviewed article

The XXL Survey. XLI. Radio AGN luminosity functions based on the GMRT $610 \ \mathrm{MHz}$ continuum observations

Authors: B. Slaus, V. Smolcic, M. Novak, S. Fotopoulou, P. Ciliegi, N. Jurlin, L. Ceraj, K. Tisanic, M. Birkinshaw, M. Bremer, L. Chiappetti, C. Horellou, M. Huynh, H. Intema, K. Kolokythas, M. Pierre, S. Raychaudhury, H. Rottgering

Abstract: We study the space density evolution of active galactic nuclei (AGN) using the $610 \ \mathrm{MHz}$ radio survey of the XXL-North field, performed with the Giant Metrewave Radio Telescope (GMRT). The survey covers an area of $30.4 \ \mathrm{deg}^2$, with a beamsize of $6.5 \ \mathrm{arcsec}$. The survey is divided into two parts, one covering an area of $11.9 \ \mathrm{deg}^2$ with $1 σ$ rms noise… ▽ More We study the space density evolution of active galactic nuclei (AGN) using the $610 \ \mathrm{MHz}$ radio survey of the XXL-North field, performed with the Giant Metrewave Radio Telescope (GMRT). The survey covers an area of $30.4 \ \mathrm{deg}^2$, with a beamsize of $6.5 \ \mathrm{arcsec}$. The survey is divided into two parts, one covering an area of $11.9 \ \mathrm{deg}^2$ with $1 σ$ rms noise of $200\ \mathrm{μJy \ beam^{-1}}$ and the other spanning $18.5 \ \mathrm{deg}^2$ with rms noise of $45\ \mathrm{μ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 (MCVCM). A flux density cut above $1\ \mathrm{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 \ \mathrm{GHz}$ radio luminosity functions of these sources using the maximum volume method. This survey allows us to probe luminosities of $ 23 \lesssim \log(L_{1.4 \ \mathrm{GHz}}[\mathrm{W/Hz}]) \lesssim 28$ up to redshifts of $z \approx 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. △ Less

Submitted 3 May, 2020; originally announced May 2020.

The ALMA Spectroscopic Survey in the HUDF: The Cosmic Dust and Gas Mass Densities in Galaxies up to $z\sim3$

Authors: Benjamin Magnelli, Leindert Boogaard, Roberto Decarli, Jorge Gónzalez-López, Mladen Novak, Gergö Popping, Ian Smail, Fabian Walter, Manuel Aravena, Roberto J. Assef, Franz Erik Bauer, Frank Bertoldi, Chris Carilli, Paulo C. Cortes, Elisabete da Cunha, Emanuele Daddi, Tanio Díaz-Santos, Hanae Inami, Robert J. Ivison, Olivier Le Fèvre, Pascal Oesch, Dominik Riechers, Hans-Walter Rix, Mark T. Sargent, Paul van der Werf , et al. (2 additional authors not shown)

Abstract: 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,… ▽ More 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. △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: 32 pages, 13 figures, Accepted for publication in ApJ

The ALMA Spectroscopic Survey in the HUDF: Deep 1.2 mm continuum number counts

Authors: Jorge González-López, Mladen Novak, Roberto Decarli, Fabian Walter, Manuel Aravena, Leindert Boogaard, Gergö Popping, Axel Weiss, Roberto J. Assef, Franz Erik Bauer, Rychard Bouwens, Paulo C. Cortes, Pierre Cox, Emanuele Daddi, Elisabete da Cunha, Tanio Díaz-Santos, Rob Ivison, Benjamin Magnelli, Dominik Riechers, Ian Smail, Paul van der Werf, Jeff Wagg

Abstract: We present the results from the 1.2 mm continuum image obtained as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin$^2$ within a primary beam response of 50% (10%) and a rms value of $9.3\thinspace{\rmμJy\thinspace beam^{-1}}$. We detect 35 sources at high significance (Fidelity $\geq0.5$), 32 of these have wel… ▽ More We present the results from the 1.2 mm continuum image obtained as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin$^2$ within a primary beam response of 50% (10%) and a rms value of $9.3\thinspace{\rmμJy\thinspace beam^{-1}}$. We detect 35 sources at high significance (Fidelity $\geq0.5$), 32 of these have well characterized near-infrared HST counterparts. We estimate the 1.2 mm number counts to flux levels of $<30\thinspace{\rmμJy}$ in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below $S_ν\sim 0.1$ mJy. In a second approach, we constrain the number counts by using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jack-knifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting $\sim$93% ($\sim$100% if we include the lower fidelity detections) of the total continuum dust emission associated to galaxies in the HUDF. The ancillary data allows us to study the dependence of the 1.2 mm number counts on redshift ($z=0-4$), galaxy dust mass (${\rm M}_{\rm dust}=10^{7}-10^{9}{\rm M}_{\odot}$), stellar mass (${\rm M}_{*}=10^{9}-10^{12}{\rm M}_{\odot}$), and star-formation rate (${\rm SFR}=1-1000\thinspace{\rm M}_{\odot}\thinspace{\rm yr^{-1}}$). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: Accepted for publication in ApJ

Detecting and Characterizing Young Quasars I: Systemic Redshifts and Proximity Zones Measurements

Authors: Anna-Christina Eilers, Joseph F. Hennawi, Roberto Decarli, Frederick B. Davies, Bram Venemans, Fabian Walter, Eduardo Bañados, Xiaohui Fan, Emanuele P. Farina, Chiara Mazzucchelli, Mladen Novak, Jan-Torge Schindler, Robert A. Simcoe, Feige Wang, Jinyi Yang

Abstract: 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… ▽ More 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}$. △ Less

Submitted 25 August, 2020; v1 submitted 5 February, 2020; originally announced February 2020.

Comments: Accepted for publication in ApJ

The z=7.54 Quasar ULAS J1342+0928 Is Hosted by a Galaxy Merger

Authors: Eduardo Banados, Mladen Novak, Marcel Neeleman, Fabian Walter, Roberto Decarli, Bram P. Venemans, Chiara Mazzucchelli, Chris Carilli, Feige Wang, Xiaohui Fan, Emanuele P. Farina, Hans-Walter Rix

Abstract: We present ALMA 0.28"x0.20" (1.4kpc x 1.0kpc) resolution observations of the [CII] 158um 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 [CII] and continuum are detected and spatially resolved. The [CII] line is ~1.5 times more extended than the continuum emission, showing an elongated and complex structure of… ▽ More We present ALMA 0.28"x0.20" (1.4kpc x 1.0kpc) resolution observations of the [CII] 158um 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 [CII] and continuum are detected and spatially resolved. The [CII] line is ~1.5 times more extended than the continuum emission, showing an elongated and complex structure of approximately 3.2kpc x 6.4kpc. Two separate peaks are clearly seen (>6sigma each) in three 100 km/s width [CII] channel maps. The velocity field of the [CII] gas does not show evidence of a coherent rotation field but rather chaotic motion, reminiscent of an on-going merger. △ Less

Submitted 30 August, 2019; originally announced September 2019.

Resolved [CII] Emission from z>6 Quasar Host -- Companion Galaxy Pairs

Authors: Marcel Neeleman, Eduardo Bañados, Fabian Walter, Roberto Decarli, Bram P. Venemans, Chris L. Carilli, Xiaohui Fan, Emanuele P. Farina, Chiara Mazzucchelli, Mladen Novak, Dominik A. Riechers, Hans-Walter Rix, Ran Wang

Abstract: We report on ~0.35"(~2 kpc) resolution observations of the [CII] and dust continuum emission from five z>6 quasar host-companion galaxy pairs obtained with the Atacama Large Millimeter/submillimeter Array. The [CII] emission is resolved in all galaxies, with physical extents of 3.2-5.4 kpc. The dust continuum is on-average 40% more compact, which results in larger [CII] deficits in the center of t… ▽ More We report on ~0.35"(~2 kpc) resolution observations of the [CII] and dust continuum emission from five z>6 quasar host-companion galaxy pairs obtained with the Atacama Large Millimeter/submillimeter Array. The [CII] emission is resolved in all galaxies, with physical extents of 3.2-5.4 kpc. The dust continuum is on-average 40% more compact, which results in larger [CII] deficits in the center of the galaxies. However, the measured [CII] deficits are fully consistent with those found at lower redshifts. Four of the galaxies show [CII] velocity fields that are consistent with ordered rotation, while the remaining six galaxies show no clear velocity gradient. All galaxies have high (~80-200 km/s) velocity dispersions, consistent with the interpretation that the interstellar medium (ISM) of these high redshift galaxies is turbulent. By fitting the galaxies with kinematic models, we estimate the dynamical mass of these systems, which range between (0.3 -> 5.4) x 1E10 Msun. For the three closest separation galaxy pairs, we observe dust and [CII] emission from gas in between and surrounding the galaxies, which is an indication that tidal interactions are disturbing the gas in these systems. Although gas exchange in these tidal interactions could power luminous quasars, the existence of quasars in host galaxies without nearby companions suggests that tidal interactions are not the only viable method for fueling their active centers. These observations corroborate the assertion that accreting supermassive black holes do not substantially contribute to the [CII] and dust continuum emission of the quasar host galaxies, and showcase the diverse ISM properties of galaxies when the universe was less than one billion years old. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Comments: 30 pages, 19 figures, Accepted for publication in ApJ

Far-Infrared Properties of the Bright, Gravitationally Lensed Quasar J0439+1634 at z=6.5

Authors: Jinyi Yang, Bram Venemans, Feige Wang, Xiaohui Fan, Mladen Novak, Roberto Decarli, Fabian Walter, Minghao Yue, Emmanuel Momjian, Charles R. Keeton, Ran Wang, Ann Zabludoff, Xue-Bing Wu, Fuyan Bian

Abstract: 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 $μ$m, [CI] 369 $μ$m, [OI] 146 $μ$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 yie… ▽ More 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 $μ$m, [CI] 369 $μ$m, [OI] 146 $μ$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}$ $μ^{-1}$ $L_{\odot}$, where $μ$ $\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}$ $μ^{-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} μ^{-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. △ Less

Submitted 30 June, 2019; originally announced July 2019.

Comments: 13 pages, 6 figures. Accepted by ApJ

An ALMA multi-line survey of the interstellar medium of the redshift 7.5 quasar host galaxy J1342+0928

Authors: Mladen Novak, Eduardo Banados, Roberto Decarli, Fabian Walter, Bram Venemans, Marcel Neeleman, Emanuele Paolo Farina, Chiara Mazzucchelli, Chris Carilli, Xiaohui Fan, Hans-Walter Rix, Feige Wang

Abstract: 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… ▽ More 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 <~5x10^4cm^-3. The data suggest that ~16% of hydrogen is in ionized form and that the HII regions have high electron densities of n_e>180 cm^-3. Our observations favor a low gas-to-dust ratio of <100 and a metallicity of the interstellar medium comparable to the Solar value. All the measurements presented here suggest that the host galaxy of J1342+0928 is highly enriched in metal and dust, despite being observed just 680 Myr after the Big Bang. △ Less

Submitted 20 June, 2019; originally announced June 2019.

Comments: 18 pages, 11 figures, 3 tables, accepted for publication in ApJ

ALMA and HST kiloparsec-scale imaging of a quasar-galaxy merger at $z\approx 6.2$

Authors: Roberto Decarli, Massimo Dotti, Eduardo Bañados, Emanuele Paolo Farina, Fabian Walter, Chris Carilli, Xiaohui Fan, Chiara Mazzucchelli, Marcel Neeleman, Mladen Novak, Dominik Riechers, Michael A. Strauss, Bram P. Venemans, Yujin Yang, Ran Wang

Abstract: 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 ext… ▽ More 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. △ Less

Submitted 12 June, 2019; originally announced June 2019.

Comments: 9 pages, 6 figures. Accepted for publication in ApJ

Radio continuum size evolution of star-forming galaxies over 0.35 < z < 2.25

Authors: E. F. Jiménez-Andrade, B. Magnelli, A. Karim, G. Zamorani, M. Bondi, E. Schinnerer, M. Sargent, E. Romano-Díaz, M. Novak, P. Lang, F. Bertoldi, E. Vardoulaki, S. Toft, V. Smolčić, K. Harrington, S. Leslie, J. Delhaize, D. Liu, C. Karoumpis, J. Kartaltepe, A. M. Koekemoer

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

Submitted 28 March, 2019; originally announced March 2019.

Comments: 21 pages, 19 figures, accepted for publication in A&A

Journal ref: A&A 625, A114 (2019)

400 pc Imaging of a Massive Quasar Host Galaxy at a Redshift of 6.6

Authors: Bram Venemans, Marcel Neeleman, Fabian Walter, Mladen Novak, Roberto Decarli, Joseph Hennawi, Hans-Walter Rix

Abstract: 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 t… ▽ More 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. △ Less

Submitted 16 April, 2019; v1 submitted 21 March, 2019; originally announced March 2019.

Comments: Published in ApJ Letters

Journal ref: The Astrophysical Journal Letters, Volume 874, Issue 2, article id. L30, 9 pp. (2019)

The ALMA Spectroscopic Survey in the HUDF: CO luminosity functions and the molecular gas content of galaxies through cosmic history

Authors: Roberto Decarli, Fabian Walter, Jorge Gónzalez-López, Manuel Aravena, Leindert Boogaard, Chris Carilli, Pierre Cox, Emanuele Daddi, Gergö Popping, Dominik Riechers, Bade Uzgil, Axel Weiss, Roberto J. Assef, Roland Bacon, Franz Erik Bauer, Frank Bertoldi, Rychard Bouwens, Thierry Contini, Paulo C. Cortes, Elisabete da Cunha, Tanio Díaz-Santos, David Elbaz, Hanae Inami, Jacqueline Hodge, Rob Ivison , et al. (12 additional authors not shown)

Abstract: 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 $ρ$(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… ▽ More 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 $ρ$(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 $ρ$(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 $ρ$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $ρ_{\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 $ρ$(H$_2$) and $ρ_{\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). △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: 22 pages, 11 figures. Paper re-submitted to ApJ after addressing the first round of comments by the referee

The evolution of the cosmic molecular gas density

Authors: Fabian Walter, Chris Carilli, Roberto Decarli, Dominik Riechers, Manuel Aravena, Franz Erik Bauer, Frank Bertoldi, Alberto Bolatto, Leindert Boogaard, Rychard Bouwens, Denis Burgarella, Caitlin Casey, Asantha Cooray, Paolo Cortes, Pierre Cox, Emanuele Daddi, Jeremy Darling, Bjorn Emonts, Jorge Gonzalez Lopez, Jacqueline Hodge, Hanae Inami, Rob Ivison, Ely Kovetz, Olivier Le Fevre, Benjamin Magnelli , et al. (16 additional authors not shown)

Abstract: One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it… ▽ More One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it is converted to stars. Over the last decade, important progress has been made in understanding the relative importance of these two factors along with the role of feedback, and the first measurements of the volume density of cold gas out to redshift 4, (the "cold gas history of the Universe") has been obtained. To match the precision of measurements of the star formation and black-hole accretion histories over the coming decades, a two orders of magnitude improvement in molecular line survey speeds is required compared to what is possible with current facilities. Possible pathways towards such large gains include significant upgrades to current facilities like ALMA by 2030 (and beyond), and eventually the construction of a new generation of radio-to-millimeter wavelength facilities, such as the next generation Very Large Array (ngVLA) concept. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: 7 pages, 2 figures, Science White paper submitted to Astro2020 Decadal Survey

A closer look at the deep radio sky: Multi-component radio sources at 3-GHz VLA-COSMOS

Authors: E. Vardoulaki, E. F. Jiménez Andrade, A. Karim, M. Novak, S. K. Leslie, K. Tisanić, V. Smolčić, E. Schinnerer, M. T. Sargent, M. Bondi, G. Zamorani, B. Magnelli, F. Bertoldi, N. Herrera Ruiz, K. P. Mooley, J. Delhaize, S. T. Myers, S. Marchesi, A. M. Koekemoer, G. Gozaliasl, A. Finoguenov, E. Middleberg, P. Ciliegi

Abstract: 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… ▽ More 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] △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: A&A accepted; 32 pages, 14 figures, 5 tables, 3 Appendices

Journal ref: A&A 627, A142 (2019)

The VLA-COSMOS 3 GHz Large Project: Average radio spectral energy distribution of highly star-forming galaxies

Authors: K. Tisanić, V. Smolčić, J. Delhaize, M. Novak, H. Intema, I. Delvecchio, E. Schinnerer, G. Zamorani, M. Bondi, E. Vardoulaki

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

Submitted 8 December, 2018; originally announced December 2018.

Comments: 21 pages, 15 figures

Journal ref: A&A 621, A139 (2019)

No evidence for enhanced [OIII] 88um emission in a z~6 quasar compared to its companion starbursting galaxy

Authors: Fabian Walter, Dominik Riechers, Mladen Novak, Roberto Decarli, Carl Ferkinhoff, Bram Venemans, Eduardo Banados, Frank Bertoldi, Chris Carilli, Xiaohui Fan, Emanuele Farina, Chiara Mazzucchelli, Marcel Neeleman, Hans-Walter Rix, Michael Strauss, Bade Uzgil, Ran Wang

Abstract: 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 com… ▽ More 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. △ Less

Submitted 30 November, 2018; originally announced November 2018.

Comments: ApJ Letters, accepted

The VLA-COSMOS 3 GHz Large Project: Star formation properties and radio luminosity functions of AGN with moderate-to-high radiative luminosities out to $z\sim6$

Authors: Lana Ceraj, Vernesa Smolčić, Ivan Delvecchio, Mladen Novak, Giovanni Zamorani, Jacinta Delhaize, Eva Schinnerer, Eleni Vardoulaki, Noelia Herrera Ruiz

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

Submitted 7 November, 2018; originally announced November 2018.

Comments: 20 pages, 14 figures

Journal ref: A&A 620, A192 (2018)

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

Authors: I. Delvecchio, V. Smolcic, G. Zamorani, D. J. Rosario, M. Bondi, S. Marchesi, T. Miyaji, M. Novak, M. T. Sargent, D. M. Alexander, J. Delhaize

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

Submitted 24 September, 2018; v1 submitted 17 September, 2018; originally announced September 2018.

Comments: 13 pages, 6 figures, 2 tables. Accepted for publication in MNRAS

"Super-deblended" Dust Emission in Galaxies: II. Far-IR to (sub)millimeter photometry and high redshift galaxy candidates in the full COSMOS field

Authors: Shuowen Jin, Emanuele Daddi, Daizhong Liu, Vernesa Smolcic, Eva Schinnerer, Antonello Calabrò, Qiusheng Gu, J. Delhaize, Ivan Delvecchio, Yu Gao, Mara Salvato, Annagrazia Puglisi, Mark Dickinson, Frank Bertoldi, Mark Sargent, M. Novak, G. E. Magdis, Itziar Aretxaga, Grant W. Wilson, Peter Capak

Abstract: 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 spec… ▽ More 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. △ Less

Submitted 12 July, 2018; originally announced July 2018.

Comments: 49 pages (+appendix), 48 figures, 3 tables, accepted for publication in ApJ

VLBA+GBT observations of the COSMOS field and radio source counts at 1.4 GHz

Authors: N. Herrera Ruiz, E. Middelberg, A. Deller, V. Smolčić, R. P. Norris, M. Novak, I. Delvecchio, P. N. Best, E. Schinnerer, E. Momjian, R. -J. Dettmar, W. Brisken, A. M. Koekemoer, N. Z. Scoville

Abstract: 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… ▽ More 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. △ Less

Submitted 27 April, 2018; originally announced April 2018.

Comments: 14 pages, 8 figures, 4 tables. Accepted for publication in A&A

Journal ref: A&A 616, A128 (2018)

The XXL Survey: XXIX. GMRT 610 MHz continuum observations

Authors: Vernesa Smolcic, Huib Intema, Bruno Slaus, Somak Raychaudhury, Mladen Novak, Cathy Horellou, Lucio Chiappetti, Jacinta Delhaize, Mark Birkinshaw, Marco Bondi, Malcolm Bremer, Paolo Ciliegi, Chiara Ferrari, Konstantinos Kolokythas, Chris Lidman, Sean L. McGee, Ray Norris, Marguerite Pierre, Huub Rottgering, Cyril Tasse, Wendy Williams

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

Submitted 16 March, 2018; originally announced March 2018.

Comments: 11 pages, 13 figures, 1 table; to appear in A&A (XXL special issue)

Journal ref: A&A 620, A14 (2018)

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

Authors: M. Novak, V. Smolcic, E. Schinnerer, G. Zamorani, I. Delvecchio, M. Bondi, J. Delhaize

Abstract: 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 redshif… ▽ More 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. △ Less

Submitted 5 March, 2018; originally announced March 2018.

Comments: Accepted for publication in A&A, 12 pages, 7 figures, 3 tables

Journal ref: A&A 614, A47 (2018)

Average radio spectral energy distribution of highly star-forming galaxies

Authors: Krešimir Tisanić, Vernesa Smolčić, Jacinta Delhaize, Mladen Novak, Huib Intema, Ivan Delvecchio, Eva Schinnerer, Giovanni Zamorani

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

Submitted 25 January, 2018; originally announced January 2018.

Comments: 4 pages, 1 figure, Proceedings of Science, IAU Symposium 333: "Peering towards Cosmic Dawn"

Cosmic evolution of AGN with moderate-to-high radiative luminosity in the COSMOS field

Authors: Lana Ceraj, Vernesa Smolčić, Ivan Delvecchio, Jacinta Delhaize, Mladen Novak

Abstract: 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$σ$ sensitivity of 2.3 $\mathrm{μJy}$/beam across the field. This provides the simultaneously largest and deepest radio continuum survey available to date with exquisite multi-wavelength coverage. The… ▽ More 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$σ$ sensitivity of 2.3 $\mathrm{μ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$. △ Less

Submitted 28 December, 2017; originally announced December 2017.

Comments: 4 pages, 1 figure, Proceedings of Science, IAU Symposium 333: "Peering towards Cosmic Dawn"

The clustering and bias of radio-selected AGN and star-forming galaxies in the COSMOS field

Authors: C. L. Hale, M. J. Jarvis, I. Delvecchio, P. W. Hatfield, M. Novak, V. Smolcic, G. Zamorani

Abstract: Dark matter haloes in which galaxies reside are likely to have a significant impact on their evolution. We investigate the link between dark matter haloes and their constituent galaxies by measuring the angular two-point correlation function of radio sources, using recently released 3 GHz imaging over $\sim 2 \ \mathrm{deg}^2$ of the COSMOS field. We split the radio source population into Star For… ▽ More Dark matter haloes in which galaxies reside are likely to have a significant impact on their evolution. We investigate the link between dark matter haloes and their constituent galaxies by measuring the angular two-point correlation function of radio sources, using recently released 3 GHz imaging over $\sim 2 \ \mathrm{deg}^2$ of the COSMOS field. We split the radio source population into Star Forming Galaxies (SFGs) and Active Galactic Nuclei (AGN), and further separate the AGN into radiatively efficient and inefficient accreters. Restricting our analysis to $z<1$, we find SFGs have a bias, $b = 1.5 ^{+0.1}_{-0.2}$, at a median redshift of $z=0.62$. On the other hand, AGN are significantly more strongly clustered with $b = 2.1\pm 0.2$ at a median redshift of 0.7. This supports the idea that AGN are hosted by more massive haloes than SFGs. We also find low-accretion rate AGN are more clustered ($b = 2.9 \pm 0.3$) than high-accretion rate AGN ($b = 1.8^{+0.4}_{-0.5}$) at the same redshift ($z \sim 0.7$), suggesting that low-accretion rate AGN reside in higher mass haloes. This supports previous evidence that the relatively hot gas that inhabits the most massive haloes is unable to be easily accreted by the central AGN, causing them to be inefficient. We also find evidence that low-accretion rate AGN appear to reside in halo masses of $M_{h} \sim 3-4 \times 10^{13}h^{-1}$M$_{\odot}$ at all redshifts. On the other hand, the efficient accreters reside in haloes of $M_{h} \sim 1-2 \times 10^{13}h^{-1}$M$_{\odot}$ at low redshift but can reside in relatively lower mass haloes at higher redshifts. This could be due to the increased prevalence of cold gas in lower mass haloes at $z \ge 1$ compared to $z<1$. △ Less

Submitted 14 November, 2017; originally announced November 2017.

Comments: 20 pages, 10 figures, 1 table, accepted by MNRAS

The infrared-radio correlation of spheroid- and disc-dominated star-forming galaxies to z $\sim$ 1.5 in the COSMOS field

Authors: Daniel Cs. Molnar, Mark T. Sargent, Jacinta Delhaize, Ivan Delvecchio, Vernesa Smolcic, Mladen Novak, Eva Schinnerer, Giovanni Zamorani, Marco Bondi, Noelia Herrera-Ruiz, Eric J. Murphy, Eleni Vardoulaki, Alexander Karim, Sarah Leslie, Benjamin Magnelli, C. Marcella Carollo, Enno Middelberg

Abstract: 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… ▽ More 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. △ Less

Submitted 20 October, 2017; originally announced October 2017.

Comments: Resubmitted to MNRAS after addressing reviewer feedback

Radio selection of the most distant galaxy clusters

Authors: E. Daddi, S. Jin, V. Strazzullo, M. T. Sargent, T. Wang, C. Ferrari, E. Schinnerer, V. Smolcic, A. Calabro, R. Coogan, J. Delhaize, I. Delvecchio, D. Elbaz, R. Gobat, Q. Gu, D. Liu, M. Novak, F. Valentino

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

Submitted 23 August, 2017; originally announced August 2017.

Comments: ApJ Letters in press

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

Authors: Drew Brisbin, Oskari Miettinen, Manuel Aravena, Vernesa Smolčić, Ivan Delvecchio, Chunyan Jiang, Benjamin Magnelli, Marcus Albrecht, Alejandra Muñoz Arancibia, Hervé Aussel, Nikola Baran, Frank Bertoldi, Matthieu Béthermin, Peter Capak, Caitlin M. Casey, Francesca Civano, Christopher C. Hayward, Olivier Ilbert, Alexander Karim, Olivier Le Fevre, Stefano Marchesi, Henry Joy McCracken, Felipe Navarrete, Mladen Novak, Dominik Riechers , et al. (6 additional authors not shown)

Abstract: 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 av… ▽ More 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.25 mJy, have a median redshift of $\tilde{z}$=2.18$\pm$0.09, while the brightest sources, 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. △ Less

Submitted 7 September, 2017; v1 submitted 18 August, 2017; originally announced August 2017.

Comments: Accepted for publication in Astronomy & Astrophysics for review 17 pages in body + 26 in appendices

Journal ref: A&A 608, A15 (2017)

An ALMA survey of submillimetre galaxies in the COSMOS field: Physical properties derived from energy balance spectral energy distribution modelling

Authors: Oskari Miettinen, Ivan Delvecchio, Vernesa Smolčić, Manuel Aravena, Drew Brisbin, Alexander Karim, Benjamin Magnelli, Mladen Novak, Eva Schinnerer, Marcus Albrecht, Hervé Aussel, Frank Bertoldi, Peter L. Capak, Caitlin M. Casey, Christopher C. Hayward, Olivier Ilbert, Huib T. Intema, Chunyan Jiang, Olivier Le Fèvre, Henry J. McCracken, Alejandra M. Muñoz Arancibia, Felipe Navarrete, Nelson D. Padilla, Dominik A. Riechers, Mara Salvato , et al. (3 additional authors not shown)

Abstract: We determine the physical properties of a sample of SMGs in the COSMOS field that were pre-selected at the observed wavelength of $λ_{\rm obs}=1.1$ mm, and followed up at $λ_{\rm obs}=1.3$ mm with ALMA. We used MAGPHYS to fit the panchromatic (ultraviolet to radio) SEDs of 124 of the target SMGs, 19.4% of which are spectroscopically confirmed. The SED analysis was complemented by estimating the ga… ▽ More We determine the physical properties of a sample of SMGs in the COSMOS field that were pre-selected at the observed wavelength of $λ_{\rm obs}=1.1$ mm, and followed up at $λ_{\rm obs}=1.3$ mm with ALMA. We used MAGPHYS to fit the panchromatic (ultraviolet to radio) SEDs of 124 of the target SMGs, 19.4% of which are spectroscopically confirmed. The SED analysis was complemented by estimating the gas masses of the SMGs by using the $λ_{\rm obs}=1.3$ mm emission as a tracer of the molecular gas. The sample median and 16th-84th percentile ranges of the stellar masses, SFRs, dust temperatures, and dust and gas masses were derived to be $\log(M_{\star}/{\rm M}_{\odot})=11.09^{+0.41}_{-0.53}$, ${\rm SFR}=402^{+661}_{-233}$ ${\rm M}_{\odot}~{\rm yr}^{-1}$, $T_{\rm dust}=39.7^{+9.7}_{-7.4}$ K, $\log(M_{\rm dust}/{\rm M}_{\odot})=9.01^{+0.20}_{-0.31}$, and $\log(M_{\rm gas}/{\rm M}_{\odot})=11.34^{+0.20}_{-0.23}$, respectively. The median gas-to-dust ratio and gas fraction were found to be $120^{+73}_{-30}$ and $0.62^{+0.27}_{-0.23}$, respectively. We found that 57.3% of our SMGs populate the main sequence (MS) of star-forming galaxies, while 41.9% of the sources lie above the MS by a factor of >3 (one source lies below the MS). The largest 3 GHz radio sizes are found among the MS sources. Those SMGs that appear irregular in the rest-frame UV are predominantly starbursts, while the MS SMGs are mostly disk-like. The larger radio-emitting sizes of the MS SMGs compared to starbursts is a likely indication of their more widespread, less intense star formation. The irregular UV morphologies of the starburst SMGs are likely to echo their merger nature. Our results suggest that the transition from high-$z$ SMGs to local ellipticals via compact, quiescent galaxies (cQGs) at $z \sim 2$ might not be universal, and the latter population might also descend from the so-called blue nuggets. △ Less

Submitted 3 July, 2017; originally announced July 2017.

Comments: 43 pages (incl. 4 appendices), 20 figures, 5 tables, accepted for publication in A&A, abstract abridged for arXiv

Journal ref: A&A 606, A17 (2017)