Search

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance., Comment: Accepted for publication in the A&A special issue`Euclid on Sky'

Dusty, distant, massive ($M_*\gtrsim 10^{11}\,\rm M_\odot$) galaxies are usually found to show a remarkable star-formation activity, contributing on the order of $25\%$ of the cosmic star-formation rate density at $z\approx3$--$5$, and up to $30\%$ at $z\sim7$ from ALMA observations. Nonetheless, they are elusive in classical optical surveys, and current near-infrared surveys are able to detect them only in very small sky areas. Since these objects have low space densities, deep and wide surveys are necessary to obtain statistically relevant results about them. Euclid will be potentially capable of delivering the required information, but, given the lack of spectroscopic features at these distances within its bands, it is still unclear if it will be possible to identify and characterize these objects. The goal of this work is to assess the capability of Euclid, together with ancillary optical and near-infrared data, to identify these distant, dusty and massive galaxies, based on broadband photometry. We used a gradient-boosting algorithm to predict both the redshift and spectral type of objects at high $z$. To perform such an analysis we make use of simulated photometric observations derived using the SPRITZ software. The gradient-boosting algorithm was found to be accurate in predicting both the redshift and spectral type of objects within the Euclid Deep Survey simulated catalog at $z>2$. In particular, we study the analog of HIEROs (i.e. sources with $H-[4.5]>2.25$), combining Euclid and Spitzer data at the depth of the Deep Fields. We found that the dusty population at $3\lesssim z\lesssim 7$ is well identified, with a redshift RMS and OLF of only $0.55$ and $8.5\%$ ($H_E\leq26$), respectively. Our findings suggest that with Euclid we will obtain meaningful insights into the role of massive and dusty galaxies in the cosmic star-formation rate over time., Comment: 18 pages, 13 figures, accepted in A&A

Star formation histories (SFH) of early (6$9$, 87% of massive galaxies, ($\log(M_\ast/M_\odot)\gtrsim$9), have SFR gradients consistent with a stochastic star-formation activity during the last 100 Myr, while this fraction drops to 15% at $z<7$. On the other hand, we see an increasing fraction of galaxies with a star-formation activity following a common stream on the SFR-$M_\ast$ plane with cosmic time, indicating that a secular mode of star-formation is emerging. We place our results in the context of the observed excess of UV emission as probed by the UV luminosity function at $z\gtrsim10$, by estimating $\sigma_{UV}$, the dispersion of the UV absolute magnitude distribution, to be of the order of 1.2mag and compare it with predictions from the literature. In conclusion, we find a transition of star-formation mode happening around $z\sim9$: Galaxies with stochastic SFHs dominates at $z\gtrsim9$, although this level of stochasticity is too low to reach those invoked by recent models to reproduce the observed UV luminosity function., Comment: Submitted to A&A

Many efforts have been done in recent years to probe the gas fraction evolution of massive quiescent galaxies (QGs); however, a clear picture has not yet been established. Recent spectroscopic confirmations at z>3 offer the chance to measure the residual gas reservoirs of massive galaxies a few hundreds of Myr after their death and to study how fast quenching proceeds in a highly star-forming Universe. Even so, stringent constraints at z$>$2 remain hardly accessible with ALMA when adopting molecular gas tracers commonly used for the quenched population. In this letter, we propose overcoming this impasse by using the carbon [CII] 158 $\mu$m emission line to systematically probe the gaseous budget of unlensed QGs at z>2.8, when these galaxies could still host non-negligible star formation on an absolute scale and when the line becomes best observable with ALMA (Bands 8 and 7). So far predominantly used for star-forming galaxies, this emission line is the best choice to probe the gas budget of spectroscopically confirmed QGs at $z>3$, reaching 2-4 and 13-30 times deeper than dust continuum (ALMA band 7) and CO(2-1)/(1-0) (VLA K-K$\alpha$ bands), respectively, at fixed integration time. Exploiting archival ALMA observations, we place conservative 3$\sigma$ upper limits on the molecular gas fraction (f$_{\rm{mol}}=M_{\rm{H_2}}/M_{\star}$) of ADF22-QG1 (f$_{\rm{mol}}$<21%), ZF-COS-20115 (f$_{\rm{mol}}$<3.2%), two of the best-studied high-z QGs in the literature, and GS-9209 (f$_{\rm{mol}}$<72%), the most distant massive QG discovered to date. The deep upper limit found for ZF-COS-20115 is 3 times lower than previously anticipated for high-z QGs suggesting, at best, the existence of a large scatter in the f$_{\rm{mol}}$ distribution of the first QGs. Lastly, we discuss the current limitations of the method and propose ways to mitigate some of them by exploiting ALMA bands 9 and 10., Comment: 7 pages, 2 figures. A&A Letters in press

We present a 1.1mm stacking analysis of moderately massive (log($M_{*}$/$M_{\odot}$) = 10.7 $\pm$ 0.2) quiescent galaxies (QGs) at $\langle z\rangle \sim1.5$, searching for cold dust continuum emission, an excellent tracer of dust and gas mass. Using both the recent GOODS-ALMA survey as well as the full suite of ALMA Band-6 ancillary data in the GOODS-S field, we report the tentative detection of dust continuum equivalent of dust mass log($M_{dust}$/$M_{\odot}$) = 7.47 $\pm$ 0.13 and gas mass log($M_{gas}$/$M_{\odot}$) = 9.42 $\pm$ 0.14. The emerging gas fraction is $f_{gas}$ = 5.3 $\pm$ 1.8%, consistent with the results of previous stacking analyses based on lower resolution sub(mm) observations. Our results support the scenario where high-z QGs have an order of magnitude larger $f_{gas}$ compared to their local counterparts and have experienced quenching with a non negligible gas reservoir in their interstellar medium - i.e. with gas retention. Subsequent analysis yields an anti-correlation between the $f_{gas}$ and the stellar mass of QGs, especially in the high mass end where galaxies reside in the most massive haloes. The $f_{gas}$ - $M_{*}$ anti-correlation promotes the selection bias as a possible solution to the tension between the stacking results pointing towards gas retention in high-z QGs of moderate $M_{*}$ and the studies of individual targets that favour a fully depleted ISM in massive (log($M_{*}$/$M_{\odot}$) high-z QGs., Comment: Accepted to A&A on the 16/03/2023

Recently, a population of compact main sequence (MS) galaxies exhibiting starburst-like properties have been identified in the GOODS-ALMA blind survey at 1.1mm. Several evolution scenarios were proposed to explain their particular physical properties (e.g., compact size, low gas content, short depletion time). In this work, we aim at studying the star formation history (SFH) of the GOODS-ALMA galaxies to understand if the so-called ``starburst (SB) in the MS'' galaxies exhibit a different star formation activity over the last Gyr compared to MS galaxies that could explain their specificity. We use the CIGALE SED modelling code to which we add non-parametric SFHs. To compare quantitatively the recent SFH of the galaxies, we define a parameter, the star formation rate (SFR) gradient that provides the angle showing the direction that a galaxy has followed in the SFR vs stellar mass plane over a given period. We show that ``SB in the MS'' have positive or weak negative gradients over the last 100, 300, and 1000 Myr, at odds with a scenario where these galaxies would be transitioning from the SB region at the end of a strong starburst phase. Normal GOODS-ALMA galaxies and ``SB in the MS'' have the same SFR gradients distributions meaning that they have similar recent SFH, despite their different properties (compactness, low depletion time). The ``SBs in the MS'' manage to maintain a star-formation activity allowing them to stay within the MS. This points toward a diversity of galaxies within a complex MS., Comment: Submitted to A&A

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., Comment: Accepted for publication in A&A. 20 pages + Appendices

We analyse measurements of the evolving stellar mass (M0) at which the bending of the star-forming main sequence (MS) occurs over 01.4 (for which we propose a revision in agreement with latest simulations). We hence argue that the MS bending is primarily due to the lessening of cold-accretion causing a reduction in available cold gas in galaxies and supports predictions of gas feeding theory. In particular, the rapidly rising M0 with redshift at z>1 is confirming evidence for the cold-streams scenario. In this picture, a progressive fueling reduction rather than its sudden suppression in halos more massive than Mshock/Mstream produces a nearly constant star-formation rate in galaxies with stellar masses larger than M0, and not their quenching, for which other physical processes are thus required. Compared to the knee M* in the stellar mass function of galaxies, M0 is significantly lower at z<1.5, and higher at z>2, suggesting that the imprint of gas deprivation on the distribution of galaxy masses happened at early times (z>1.5-2). The typical mass at which galaxies inside the MS become bulge-dominated evolves differently from M0, consistent with the idea that bulge-formation is a distinct process from the phasing-out of cold-accretion., Comment: A&A Letters, in press

We use deep Herschel observations of the complete 2Jy sample of powerful radio AGNs in the local universe (0.05 < z < 0.7) to probe their cool interstellar medium (ISM) contents and star-forming properties, comparing them against other samples of nearby luminous AGNs and quiescent galaxies. This allows us to investigate triggering and feedback mechanisms. We find that the dust masses of the strong-line radio galaxies (SLRGs) in our sample are similar to those of radio-quiet quasars, and that their median dust mass (Mdust = 2 x 10^7 Msun) is enhanced by a factor ~200 compared to that of non-AGN ellipticals, but lower by a factor ~16 relative to that of local ultra-luminous infrared galaxies (UILRGs). Along with compelling evidence for merger signatures in optical images, the SLRGs in our sample also show relatively high star-formation efficiencies, despite the fact that many of them fall below the main sequence for star forming galaxies. Together, these results suggest that most of our SLRGs have been re-triggered by late-time mergers that are relatively minor in terms of their gas contents. In comparison with the SLRGs, the radio AGNs with weak optical emission lines (WLRGs) and edge-darkened radio jets (FRIs) have both lower cool ISM masses and star-formation rates (by a factor of >30), consistent with being fuelled by a different mechanism (e.g. the direct accretion of hot gas)., Comment: 20 pages, 8 figures, accepted for publication in MNRAS

We present Keck Cosmic Web Imager (KCWI) observations of giant Lya halos surrounding 9 galaxy groups and clusters at 2

Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived physical properties, such as star formation rates, gas fractions, depletion timescales, and dust temperatures for the galaxy sample built from the survey. There exists a subset of galaxies that exhibit starburst-like short depletion timescales, but they are located within the scatter of the so-called main sequence of SFGs. These are dubbed starbursts in the main sequence and display the most compact star formation and they are characterized by the shortest depletion timescales, lowest gas fractions, and highest dust temperatures of the galaxy sample, compared to typical SFGs at the same stellar mass and redshift. They are also very massive, accounting for $\sim 60\%$ of the most massive galaxies in the sample ($\log (M_{\rm{*}}/M_{\odot}) > 11.0$). We find trends between the areas of the ongoing star formation regions and the derived physical properties for the sample, unveiling the role of compact star formation as a physical driver of these properties. Starbursts in the main sequence appear to be the extreme cases of these trends. We discuss possible scenarios of galaxy evolution to explain the results drawn from our galaxy sample. Our findings suggest that the star formation rate is sustained in SFGs by gas and star formation compression, keeping them within the main sequence even when their gas fractions are low and they are presumably on the way to quiescence., Comment: Accepted by A&A. 26 pages, 13 figures

We study the radio spectral properties of 2,094 star-forming galaxies (SFGs) by combining our early science data from the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey with VLA, GMRT radio data, and rich ancillary data in the COSMOS field. These SFGs are selected at VLA 3GHz, and their flux densities from MeerKAT 1.3GHz and GMRT 325MHz imaging data are extracted using the "super-deblending" technique. The median radio spectral index is $\alpha_{\rm 1.3GHz}^{\rm 3GHz}=-0.80\pm0.01$ without significant variation across the rest-frame frequencies ~1.3-10GHz, indicating radio spectra dominated by synchrotron radiation. On average, the radio spectrum at observer-frame 1.3-3GHz slightly steepens with increasing stellar mass with a linear fitted slope of $\beta=-0.08\pm0.01$, which could be explained by age-related synchrotron losses. Due to the sensitivity of GMRT 325MHz data, we apply a further flux density cut at 3GHz ($S_{\rm 3GHz}\ge50\,\mu$Jy) and obtain a sample of 166 SFGs with measured flux densities at 325MHz, 1.3GHz, and 3GHz. On average, the radio spectrum of SFGs flattens at low frequency with the median spectral indices of $\alpha^{\rm 1.3GHz}_{\rm 325MHz}=-0.59^{+0.02}_{-0.03}$ and $\alpha^{\rm 3.0GHz}_{\rm 1.3GHz}=-0.74^{+0.01}_{-0.02}$. At low frequency, our stacking analyses show that the radio spectrum also slightly steepens with increasing stellar mass. By comparing the far-infrared-radio correlations of SFGs based on different radio spectral indices, we find that adopting $\alpha_{\rm 1.3GHz}^{\rm 3GHz}$ for $k$-corrections will significantly underestimate the infrared-to-radio luminosity ratio ($q_{\rm IR}$) for >17% of the SFGs with measured flux density at the three radio frequencies in our sample, because their radio spectra are significantly flatter at low frequency (0.33-1.3GHz)., Comment: 15 pages, 6 figures, 1 table; accepted for publication in MNRAS on 2021 Aug 4

Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas with inhomogeneous sensitivities, and suffer from modest $uv$ coverage coming from single array configurations. GOODS-ALMA is a 1.1mm galaxy survey over a continuous area of 72.42arcmin$^2$ at a homogeneous sensitivity. In this version 2.0, we present a new low resolution dataset and its combination with the previous high resolution dataset from the survey, improving the $uv$ coverage and sensitivity reaching an average of $\sigma = 68.4\mu$Jy beam$^{-1}$. A total of 88 galaxies are detected in a blind search (compared to 35 in the high resolution dataset alone), 50% at $S/N_{peak} \geq 5$ and 50% at $3.5 \leq S/N_{peak} \leq 5$ aided by priors. Among them, 13 out of the 88 are optically dark or faint sources ($H$- or $K$-band dropouts). The sample dust continuum sizes at 1.1mm are generally compact, with a median effective radius of $R_{e} = 0"10 \pm 0"05$ (a physical size of $R_{e} = 0.73 \pm 0.29$kpc at the redshift of each source). Dust continuum sizes evolve with redshift and stellar mass resembling the trends of the stellar sizes measured at optical wavelengths, albeit a lower normalization compared to those of late-type galaxies. We conclude that for sources with flux densities $S_{1.1mm} > 1$mJy, compact dust continuum emission at 1.1mm prevails, and sizes as extended as typical star-forming stellar disks are rare. The $S_{1.1mm} < 1$mJy sources appear slightly more extended at 1.1mm, although they are still generally compact below the sizes of typical star-forming stellar disks., Comment: Accepted by A&A. 30 pages, 18 figures

The mean size ( effective radius $R_e$) of Massive Galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$) is observed to increase steadily with cosmic time. It is still unclear whether this trend originates from the size growth of individual galaxies (via, e.g., mergers and/or AGN feedback) or from the inclusion of larger galaxies entering the selection at later epochs (progenitor bias). We here build a data-driven, flexible theoretical framework to probe the structural evolution of MGs. We assign galaxies to dark matter haloes via stellar mass-halo mass (SMHM) relations with varying high-mass slopes and scatters $\sigma_{\rm SMHM}$ in stellar mass at fixed halo mass, and assign sizes to galaxies using an empirically-motivated, constant and linear relationship between $R_e$ and the host dark matter halo radius $R_h$. We find that: 1) the fast mean size growth of MGs is well reproduced independently of the shape of the input SMHM relation; 2) the numbers of compact MGs grow steadily until $z\gtrsim2$ and fall off at lower redshifts, suggesting a lesser role of progenitor bias at later epochs; 3) a time-independent scatter $\sigma_{\rm SMHM}$ is consistent with a scenario in which compact starforming MGs transition into quiescent MGs in a few $10^8$yr with a negligible structural evolution during the compact phase, while a scatter increasing at high redshift implies significant size growth during the starforming phase. A robust measurement of the size function of MGs at high redshift can set strong constraints on the scatter of the SMHM relation and, by extension, on models of galaxy evolution., Comment: Accepted for publication in MNRAS. Comments still welcome

Studying the resolved stellar populations of the different structural components which build massive galaxies directly unveils their assembly history. We aim at characterizing the stellar population properties of a representative sample of bulges and pure spheroids in massive galaxies ($M_{\star}>10^{10}$ M$_{\odot}$) in the GOODS-N field. We take advantage of the spectral and spatial information provided by SHARDS and HST data to perform the multi-image spectro-photometrical decoupling of the galaxy light. We derive the spectral energy distribution separately for bulges and disks in the redshift range $0.14

In this work we study a system of two galaxies, Astarte and Adonis, at z $\sim $2 when the Universe was undergoing its peak of star formation activity. Astarte is a dusty star-forming galaxy at the massive-end of the main sequence (MS) and Adonis is a less-massive, bright in ultraviolet (UV), companion galaxy with an optical spectroscopic redshift. We analyse the physical properties of this system, and probe the gas mass of Astarte with its ALMA CO emission, to investigate whether this ultra-massive galaxy is quenching or not. We use CIGALE - a spectral energy distribution modeling code - to derive the key physical properties of Astarte and Adonis, mainly their star formation rates (SFRs), stellar masses, and dust luminosities. We inspect the variation of the physical parameters depending on the assumed dust attenuation law. We also estimate the molecular gas mass of Astarte from its CO emission, using different $\alpha_{CO}$ and transition ratios ($r_{31}$) and discuss the implication of the various assumptions on the gas mass derivation. We find that Astarte exhibits a MS-like star formation activity, while Adonis is undergoing a strong starburst (SB) phase. The molecular gas mass of Astarte is far below the gas fraction of typical star-forming galaxies at z=2. This low gas content and high SFR, result in a depletion time of $0.22\pm0.07$ Gyrs, slightly shorter than what is expected for a MS galaxy at this redshift. The CO luminosity versus the total IR luminosity suggests a MS-like activity assuming a galactic conversion factor and a low transition ratio. The SFR of Astarte is of the same order using different attenuation laws, unlike its stellar mass that increases using shallow attenuation laws. We discuss these properties and suggest that Astarte might be experiencing a recent decrease of star formation activity and is quenching through the MS following a SB epoch., Comment: 10 pages, 7 figures, accepted for publication in A&A. Abstract abridged for arXiv submission

We characterise the basic far-IR (FIR) properties and the gas mass fraction of massive ( ~ 11.0) quiescent galaxies (QGs) and explore how these evolve from z = 2.0 to the present day. We use robust, multi-wavelength (mid- to far-IR and sub-millimetre to radio) stacking ensembles of homogeneously selected and mass complete samples of log(M*/Msun) > 10.8 QGs. We find that the dust to stellar mass ratio (Md/M*) rises steeply as a function of redshift up to z~1.0 and then remains flat at least out to z = 2.0. Using Md as a proxy of gas mass (Mgas), we find a similar trend for the evolution of the gas mass fraction (fgas) with z > 1.0 QGs having fgas ~ 7.0% (for solar metallicity). This fgas is 3 - 10 times lower than that of normal star forming galaxies (SFGs) at their corresponding redshift but ~3 and ~10 times larger compared to that of z = 0.5 and local QGs. Furthermore, the inferred gas depletion time scales are comparable to that of local SFGs and systematically longer than that of main sequence galaxies at their corresponding redshifts. Our analysis also reveals that the average dust temperature (Td) of massive QGs remains roughly constant (< Td > = 21.0 \pm 2.0K) at least out to z ~ 2.0 and is substantially colder (~ 10K) compared to that of z > 0 SFGs. This motivated us to construct and release a redshift-invariant template IR SED, that we use to make predictions for ALMA observations and to explore systematic effects in the Mgas estimates of massive, high-z QGs. Finally, we discuss how a simple model that considers progenitor-bias can effectively reproduce the observed evolution of Md/M* and fgas. Our results indicate universal initial interstellar medium conditions for quenched galaxies and a large degree of uniformity in their internal processes across cosmic time., Comment: Accepted for publication in A&A

The high-frequency radio sky has historically remained largely unexplored due to the typical faintness of sources in this regime, and the modest survey speed compared to observations at lower frequencies. However, high-frequency radio surveys present an invaluable tracer of high-redshift star-formation, as they directly target the faint radio free-free emission. We present deep continuum observations at 34 GHz in the COSMOS and GOODS-North fields from the Karl G. Jansky Very Large Array (VLA), as part of the COLDz survey. The deep COSMOS mosaic spans $\sim10\text{arcmin}^2$ down to $\sigma=1.3\mu\text{Jy beam}^{-1}$, while the wider GOODS-N observations cover $\sim50\text{arcmin}^2$ to $\sigma=5.3\mu\text{Jy beam}^{-1}$. We present the deepest 34 GHz radio number counts to date, with five and thirteen continuum detections in COSMOS and GOODS-N, respectively. Nine galaxies show 34 GHz continuum emission that is originating from star-formation, although for two sources this is likely due to thermal emission from dust. Utilizing deep ancillary radio data at 1.4, 3, 5 and 10 GHz, we decompose the spectra of the remaining seven star-forming galaxies into their synchrotron and thermal free-free components, finding typical thermal fractions and synchrotron spectral indices comparable to those observed in local star-forming galaxies. Using calibrations from the literature, we determine free-free star-formation rates (SFRs), and show that these are in agreement with SFRs from spectral energy distribution fitting and the far-infrared/radio correlation. Our observations place strong direct constraints on the high-frequency radio emission in typical galaxies at high-redshift, and provide some of the first insight in what is set to become a key area of study with future radio facilities as the Square Kilometer Array Phase 1 and next-generation VLA., Comment: Submitted to ApJ; 22 pages + appendices; 13 figures in main text

We report the discovery of two new giant radio galaxies (GRGs) using the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey. Both GRGs were found within a 1 deg^2 region inside the COSMOS field. They have redshifts of z=0.1656 and z=0.3363 and physical sizes of 2.4Mpc and 2.0Mpc, respectively. Only the cores of these GRGs were clearly visible in previous high resolution VLA observations, since the diffuse emission of the lobes was resolved out. However, the excellent sensitivity and uv coverage of the new MeerKAT telescope allowed this diffuse emission to be detected. The GRGs occupy a unpopulated region of radio power - size parameter space. Based on a recent estimate of the GRG number density, the probability of finding two or more GRGs with such large sizes at z<0.4 in a ~1deg^2 field is only 2.7x10^-6, assuming Poisson statistics. This supports the hypothesis that the prevalence of GRGs has been significantly underestimated in the past due to limited sensitivity to low surface brightness emission. The two GRGs presented here may be the first of a new population to be revealed through surveys like MIGHTEE which provide exquisite sensitivity to diffuse, extended emission., Comment: 15 pages; 8 figures, accepted for publication in MNRAS on 2020 Dec 9

Tracing the emergence of the massive quiescent galaxy (QG) population requires the build-up of reliable quenched samples. We present Hubble Space Telescope WFC3/G141 grism spectra of 10 quiescent galaxy candidates selected at $2.5

Several works in the past decade have used the ratio between total (rest 8-1000$\mu$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., Comment: Accepted for publication in A&A. 22 pages + Appendices. 24 figures, 4 tables

We present Herschel-PACS spectroscopy of four main-sequence star-forming galaxies at z~1.5. We detect [OI]63micron line emission in BzK-21000 at z=1.5213, and measure a line luminosity, L([OI]63micron) = (3.9+/-0.7)x1.E+9 Lsun. Our PDR modelling of the interstellar medium in BzK-21000 suggests a UV radiation field strength, G~320 G0, and gas density, n~1800 cm-3, consistent with previous LVG modelling of the molecular CO line excitation. The other three targets in our sample are individually undetected in these data, and we perform a spectral stacking analysis which yields a detection of their average emission and an [OI]63micron line luminosity, L([OI]63micron) =(1.1+/-0.2)x1E+9 Lsun. We find that the implied luminosity ratio, L([OI]63micron)/L(IR), of the undetected BzK-selected star-forming galaxies broadly agrees with that of low-redshift star-forming galaxies, while BzK-21000 has a similar ratio to that of a dusty star-forming galaxy at z~6. The high [OI]63micron line luminosities observed in BzK-21000 and the $z \sim 1 -3$ dusty and submm luminous star-forming galaxies may be associated with extended reservoirs of low density, cool neutral gas., Comment: accepted for publication in MNRAS

The MOONS instrument possesses an exceptional combination of large multiplexing, high sensitivity, broad simultaneous spectral coverage (from optical to near-infrared bands), large patrol area and high fibre density. These properties provide the unprecedented potential of enabling, for the very first time, SDSS-like surveys around Cosmic Noon (z~1-2.5), when the star formation rate in the Universe peaked. The high-quality spectra delivered by MOONS will sample the same nebular and stellar diagnostics observed in extensive surveys of local galaxies, providing an accurate and consistent description of the evolution of various physical properties of galaxies, and hence a solid test of different scenarios of galaxy formation and transformation. Most importantly, by spectroscopically identifying hundreds of thousands of galaxies at high redshift, the MOONS surveys will be capable of determining the environments in which primeval galaxies lived and will reveal how such environments affected galaxy evolution. In this article, we specifically focus on the main Guaranteed Time Observation (GTO) MOONS extragalactic survey, MOONRISE, by providing an overview of its scientific goals and observing strategy., Comment: 6 pages, 6 figures, Published on the ESO Messenger No. 180

In this paper we study the properties of the six optically dark galaxies detected in the 69 arcmin^2 GOODS-ALMA 1.1mm continuum survey. While none of them are listed in the deepest H-band based CANDELS catalog in the GOODS-South field down to H=28.16AB, we were able to de-blend two of them from their bright neighbor and measure an $H$-band flux for them. We note that AGS4 and AGS15 have H=25.23, 27.11AB respectively. Their extreme proximity (0.50", 0.27") to a bright optical source and their extreme faintness prevented them from being included in the H-band catalog. We present the spectroscopic scan follow-up of five of the six sources with ALMA band 4. All are detected in the 2mm continuum with signal-to-noise ratios higher than eight. One emission line is detected in AGS4 (\nu_{obs} =151.44GHz with a S/N=8.58) and AGS17 (\nu_{obs} =154.78GHz with a S/N=10.23), which we interpret in both cases as being due to the CO(6-5) line at z^{AGS4}_{spec}=3.556 and z^{AGS4}_{spec}=3.467, respectively. These redshifts match both the probability distribution of the photometric redshifts derived from the UV to near-infrared spectral energy distributions (SEDs) and the far-infrared SEDs for typical dust temperatures of galaxies at these redshifts. We present evidence that nearly 70% (4/6 of galaxies) of the optically dark galaxies belong to the same overdensity of galaxies at z~3.5. overdensity The most massive one, AGS24 (M_{\star} = 10^{11.32^{+0.02}_{-0.19}} M_{\odot}), is the most massive galaxy without an active galactic nucleus (AGN) at $z$\,>\,3 in the GOODS-ALMA field. It falls in the very center of the peak of the galaxy surface density, which suggests that the surrounding overdensity is a proto-cluster in the process of virialization and that AGS24 is the candidate progenitor of the future brightest cluster galaxy (BCG)., Comment: 17 pages, 13 figures, 3 tables, accepted by A&A

We study the co-evolution between the black hole accretion rate (BHAR) and the star formation rate (SFR) in different galaxy life phases: main sequence star-forming galaxies, quiescent and starburst galaxies at different cosmic epochs. We take advantage of the X-ray data from the Chandra COSMOS-Legacy survey and of the extensive multiwavelength ancillary observations in the COSMOS field presented in the COSMOS2015 catalog. We perform an X-ray stacking analysis and combine it with detected sources, in a broad redshift interval ($0.1

We have discovered a 300kpc-wide giant Lya nebula centered on the massive galaxy group RO-1001 at z=2.91 in the COSMOS field. Keck Cosmic Web Imager observations reveal three cold gas filaments converging into the center of the potential well of its ~4x10^13Msun dark matter halo, hosting 1200Msun/yr of star formation as probed by ALMA and NOEMA observations. The nebula morphological and kinematics properties and the prevalence of blueshifted components in the Lya spectra are consistent with a scenario of gas accretion. The upper limits on AGN activity and overall energetics favor gravity as the primary Lya powering source and infall as the main source of gas flows to the system. Although interpretational difficulties remain, with outflows and likely also photoionization with ensuing recombination still playing a role, this finding provides arguably an ideal environment to quantitatively test models of cold gas accretion and galaxy feeding inside an actively star-forming massive halo at high redshift., Comment: 20 pages, 12 figures, 3 tables. Accepted for publication in A&A

We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA are among the most massive at $z$ = 2-4 (M$_{\star,med}$ = 8.5$ \times$ 10$^{10}$ M$_\odot$) and are either starburst or located in the upper part of the galaxy star-forming main sequence. A significant portion of our galaxy population ($\sim$ 40%), located at $z\sim$ 2.5-3, exhibits abnormally low gas fractions. The sizes of these galaxies, measured with ALMA, are compatible with the trend between $H$-band size and stellar mass observed for $z\sim2$ elliptical galaxies suggesting that they are building compact bulges. We show that there is a strong link between star formation surface density (at 1.1 mm) and gas depletion time: the more compact a galaxy's star-forming region is, the shorter its lifetime will be (without gas replenishment). The identified compact sources associated with relatively short depletion timescales ($\sim$100 Myr), are the ideal candidates to be the progenitors of compact elliptical galaxies at $z$ $\sim$ 2., Comment: 18 pages, 8 figures, Accepted for publication in A&A

In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $\sigma$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $\mu$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$\sigma$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$\sigma$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues., Comment: 18 pages, 13 figures. Accepted for publication in A&A

We have obtained spectroscopic confirmation with Hubble Space Telescope WFC3/G141 of a first sizeable sample of nine quiescent galaxies at 2.4

We present a new modeling of the X-ray luminosity function (XLF) of Active Galactic Nuclei (AGN) out to z$\sim$3, dissecting the contribution of main-sequence (MS) and starburst (SB) galaxies. For each galaxy population, we convolved the observed galaxy stellar mass (M$_{\star}$) function with a grid of M$_{\star}$-independent Eddington ratio ($\lambda_{\rm EDD}$) distributions, normalised via empirical black hole accretion rate (BHAR) to star formation rate (SFR) relations. Our simple approach yields an excellent agreement with the observed XLF since z$\sim$3. We find that the redshift evolution of the observed XLF can only be reproduced through an intrinsic flattening of the $\lambda_{\rm EDD}$ distribution, and with a positive shift of the break $\lambda^{*}$, consistent with an anti-hierarchical behavior. The AGN accretion history is predominantly made by massive (10$^{10}<$M$_{\star}<$10$^{11}$ M$_{\odot}$) MS galaxies, while SB-driven BH accretion, possibly associated with galaxy mergers, becomes dominant only in bright quasars, at $\log$(L$_{\rm X}$/erg s$^{-1}$)$>$44.36 + 1.28$\cdot$(1+z). We infer that the probability of finding highly-accreting ($\lambda_{\rm EDD}>$ 10%) AGN significantly increases with redshift, from 0.4% (3.0%) at z=0.5 to 6.5% (15.3%) at z=3 for MS (SB) galaxies, implying a longer AGN duty cycle in the early Universe. Our results strongly favor a M$_{\star}$-dependent ratio between BHAR and SFR, as BHAR/SFR $\propto$ M$_{\star}^{0.73[+0.22,-0.29]}$, supporting a non-linear BH buildup relative to the host. Finally, this framework opens potential questions on super-Eddington BH accretion and different $\lambda_{\rm EDD}$ prescriptions for understanding the cosmic BH mass assembly., Comment: Accepted for publication in ApJ. 25 pages, 15 figures, 2 tables. The best X-ray luminosity functions are available at https://github.com/idelvecchio/XLF_Delvecchio2020.git

[Abridged] We discovered in the Herschel Reference Survey an extremely bright IR source with $S_{500}$~120mJy (Red Virgo 4 - RV4). Based on IRAM/EMIR and IRAM/NOEMA detections of the CO(5-4), CO(4-3), and [CI] lines, RV4 is located at z=4.724, yielding a total observed L$_{IR}$ of 1.1+/-0.6x0$^{14}$L$_{\odot}$. At the position of the Herschel emission, three blobs are detected with the VLA at 10cm. The CO(5-4) line detection of each blob confirms that they are at the same redshift with the same line width, indicating that they are multiple images of the same source. In Spitzer and deep optical observations, two sources, High-z Lens 1 (HL1) West and HL1 East, are detected at the center of the three VLA/NOEMA blobs. These two sources are placed at z=1.48 with XSHOOTER spectra, suggesting that they could be merging and gravitationally lensing the emission of RV4. HL1 is the second most distant lens known to date in strong lensing systems. The Einstein radius of the lensing system is 2.2"+/-0.2 (20kpc). The high redshift of HL1 and the large Einstein radius are highly unusual for a strong lensing system. We present the ISM properties of the background source RV4. Different estimates of the gas depletion time yield low values suggesting that RV4 is a SB galaxy. Among all high-z SMGs, this source exhibits one of the lowest L$_{[CI]}$ to L$_{IR}$ ratios, 3.2+/-0.9x10$^{-6}$, suggesting an extremely short gas tdepl of only 14+/-5Myr. It also shows a relatively high L$_{[CI]}$ to L$_{CO(4-3)}$ ratio (0.7+/-0.2) and low L$_{CO(5-4)}$ to L$_{IR}$ ratio (only ~50% of the value expected for normal galaxies) hinting a low density of gas. Finally, we discuss that the short tdepl of RV4 can be explained by either a very high SFE, which is difficult to reconcile with major mergers simulations of high-z galaxies, or a rapid decrease of SF, which would bias the estimate of tdepl toward low value., Comment: 14 pages, 15 figures, accepted for publication in A&A

We present new measurements of the cosmic cold molecular gas evolution out to redshift 6 based on systematic mining of the ALMA public archive in the COSMOS deep field (A3COSMOS). Our A3COSMOS dataset contains ~700 galaxies (0.3 < z < 6) with high-confidence ALMA detections in the (sub-)millimeter continuum and multi-wavelength spectral energy distributions (SEDs). Multiple gas mass calibration methods are compared and biases in band conversions (from observed ALMA wavelength to rest-frame Rayleigh-Jeans(RJ)-tail continuum) have been tested. Combining our A3COSMOS sample with ~1,000 CO-observed galaxies at 0 < z < 4 (75% at z < 0.1), we parameterize galaxies' molecular gas depletion time and molecular gas to stellar mass ratio (gas fraction) each as a function of the stellar mass, offset from the star-forming main sequence (Delta MS) and cosmic age (or redshift). Our proposed functional form provides a statistically better fit to current data (than functional forms in the literature), and implies a "downsizing" effect (i.e., more-massive galaxies evolve earlier than less-massive ones) and "mass-quenching" (gas consumption slows down with cosmic time for massive galaxies but speeds up for low-mass ones). Adopting galaxy stellar mass functions and applying our function for gas mass calculation, we for the first time infer the cosmic cold molecular gas density evolution out to redshift 6 and find agreement with CO blind surveys as well as semi-analytic modeling. These together provide a coherent picture of cold molecular gas, SFR and stellar mass evolution in galaxies across cosmic time., Comment: 41 pages, 23 figures, 3 tables, resubmitted to ApJ after addressing referee's comments. Meanwhile we provide a Python package to calculate our gas evolution function ( https://ascl.net/code/v/2377 )

Supermassive black holes (SMBHs) are tightly correlated with their hosts but the origin of such connection remains elusive. To explore the cosmic build-up of this scaling relation, we present an empirically-motivated model that tracks galaxy and SMBH growth down to z=0. Starting from a random mass seed distribution at z=10, we assume that each galaxy evolves on the star-forming "main sequence" (MS) and each BH follows the recently-derived stellar mass (M$_{\star}$) dependent ratio between BH accretion rate and star formation rate, going as BHAR/SFR$\propto$M$_{\star}^{0.73[+0.22,-0.29]}$. Our simple recipe naturally describes the BH-galaxy build-up in two stages. At first, the SMBH lags behind the host that evolves along the MS. Later, as the galaxy grows in M$_{\star}$, our M$_{\star}$-dependent BHAR/SFR induces a super-linear BH growth, as M$_{\rm BH}$$\propto$M$_{\star}^{1.7}$. According to this formalism, smaller BH seeds increase their relative mass faster and earlier than bigger BH seeds, at fixed M$_{\star}$, thus setting along a gradually tighter M$_{\rm BH}$-M$_{\star}$ locus towards higher M$_{\star}$. Assuming reasonable values of the radiative efficiency $\epsilon \sim$0.1, our empirical trend agrees with both high-redshift model predictions and intrinsic M$_{\rm BH}$-M$_{\star}$ relations of local BHs. We speculate that the observed non-linear BH-galaxy build-up is reflected in a twofold behavior with dark matter halo mass (M$_{\rm DM}$), displaying a clear turnover at M$_{\rm DM}\sim$2$\times$10$^{12}$M$_{\odot}$. While Supernovae-driven feedback suppresses BH growth in smaller halos (BHAR/SFR$\propto$M$_{\rm DM}^{1.6}$), above the M$_{\rm DM}$ threshold cold gas inflows possibly fuel both BH accretion and star formation in a similar fashion (BHAR/SFR$\propto$M$_{\rm DM}^{0.3}$)., Comment: Accepted for publication in ApJ Letter. 8 pages, 5 figures

We present an enhanced version of the multiwavelength spectral modeling code MAGPHYS that allows the estimation of galaxy photometric redshift and physical properties (e.g., stellar mass, star formation rate, dust attenuation) simultaneously, together with robust characterization of their uncertainties. The self-consistent modeling over ultraviolet to radio wavelengths in MAGPHYS+photo-z is unique compared to standard photometric redshift codes. The broader wavelength consideration is particularly useful for breaking certain degeneracies in color vs. redshift for dusty galaxies with limited observer-frame ultraviolet and optical data (or upper limits). We demonstrate the success of the code in estimating redshifts and physical properties for over 4,000 infrared-detected galaxies at 0.4

We investigate the contribution of clumps and satellites to the galaxy mass assembly. We analyzed spatially-resolved Hubble Space Telescope observations (imaging and slitless spectroscopy) of 53 star-forming galaxies at z ~ 1 - 3. We created continuum and emission line maps and pinpointed residual "blobs" detected after subtracting the galaxy disk. Those were separated into compact (unresolved) and extended (resolved) components. Extended components have sizes ~ 2 kpc and comparable stellar mass and age as the galaxy disks, whereas the compact components are 1.5 dex less massive and 0.4 dex younger than the disks. Furthermore the extended blobs are typically found at larger distances from the galaxy barycenter than the compact ones. Prompted by these observations and by the comparison with simulations, we suggest that compact blobs are in-situ formed clumps, whereas the extended ones are accreting satellites. Clumps and satellites enclose respectively ~ 20% and ~< 80% of the galaxy stellar mass, ~ 30% and ~ 20% of its star formation rate. Considering the compact blobs, we statistically estimated that massive clumps (Mstar >~ 10^9 Msun) have lifetimes of ~ 650 Myr, and the less massive ones (10^8 < Mstar < 10^9 Msun) of ~ 145 Myr. This supports simulations predicting long-lived clumps (lifetime > 100 Myr). Finally, ~< 30% (13%) of our sample galaxies are undergoing single (multiple) merger(s), they have a projected separation ~< 10 kpc, and the typical mass ratio of our satellites is 1:5 (but ranges between 1:10 and 1:1), in agreement with literature results for close pair galaxies., Comment: 20 pages (plus appendix), 8 figures, 2 tables. Accepted for publication in MNRAS

We present Atacama Large Millimetre Array and Atacama Compact Array observations of the Sunyaev-Zel'dovich effect in the z = 2 galaxy cluster Cl J1449+0856, an X-ray-detected progenitor of typical massive clusters in the present day Universe. While in a cleaned but otherwise untouched 92 GHz map of this cluster, little to no negative signal is visible, careful subtraction of known sub-millimetre emitters in the uv plane reveals a decrement at 5$\sigma$ significance. The total signal is -190$\pm$36 $\mu$Jy, with a peak offset by 5"-9" ($\sim$50 kpc) from both the X-ray centroid and the still-forming brightest cluster galaxy. A comparison of the recovered uv-amplitude profile of the decrement with different pressure models allows us to derive total mass constraints consistent with the $\sim$6$\times$10$^{13}$ M$_{\odot}$ estimated from X-ray data. Moreover, we find no strong evidence for a deviation of the pressure profile with respect to local galaxy clusters, although a slight tension at small-to-intermediate spatial scales suggests a flattened central profile, opposite to what seen in a cool core and possibly an AGN-related effect. This analysis of the lowest mass single SZ detection so far illustrates the importance of interferometers when observing the SZ effect in high-redshift clusters, the cores of which cannot be considered quiescent, such that careful subtraction of galaxy emission is necessary., Comment: 11 pages, 6 figures (+appendices); accepted for publication in Astronomy & Astrophysics

We present ALMA Band~7 observations at $850\mu$m of 20 luminous ($\log\, L_{\rm bol}>46.9$ [erg s$^{-1}$]) unobscured quasars at $z\sim2$. We detect continuum emission for 19/20 quasars. After subtracting an AGN contribution, we measure the total far-IR luminosity for 18 quasars, assuming a modified blackbody model, and attribute the emission as indicative of the star formation rate (SFR). Our sample can be characterized with a log-normal SFR distribution having a mean of 140 $M_\odot$ yr$^{-1}$ and a dispersion of 0.5 dex. Based on an inference of their stellar masses, the SFRs are similar, in both the mean and dispersion, with star-forming main-sequence galaxies at the equivalent epoch. Thus, there is no evidence for a systematic enhancement or suppression (i.e., regulation or quenching) of star formation in the hosts of the most luminous quasars at $z\sim2$. These results are consistent with the Magneticum cosmological simulation, while in disagreement with a widely recognized phenomenological model that predicts higher SFRs than observed here based on the high bolometric luminosities of this sample. Furthermore, there is only a weak relation between SFR and accretion rate onto their supermassive black holes both for average and individual measurements. We interpret these results as indicative of star formation and quasar accretion being fed from the available gas reservoir(s) in their host with a disconnect due to their different physical sizes, temporal scales, and means of gas processing., Comment: 20 pages, 12 figures, re-submitted to MNRAS

We report Atacama Large Millimetre Array (ALMA) observations of four high-redshift dusty star-forming galaxy candidates selected from far-Infrared (FIR)/submm observations in the COSMOS field. We securely detect all galaxies in the continuum and spectroscopically confirm them at z=3.62--5.85 using ALMA 3mm line scans, detecting multiple CO and/or [CI] transitions. This includes the most distant dusty galaxy currently known in the COSMOS field, ID85001929 at z=5.847. These redshifts are lower than we had expected as these galaxies have substantially colder dust temperatures (i.e., their SEDs peak at longer rest frame wavelengths) than most literature sources at z>4. The observed cold dust temperatures are best understood as evidence for optically thick dust continuum in the FIR, rather than the result of low star formation efficiency with rapid metal enrichment. We provide direct evidence that, given their cold spectral energy distributions, CMB plays a significant role biasing their observed Rayleigh-Jeans (RJ) slopes to unlikely steep values and, possibly, reducing their CO fluxes by a factor of two. We recover standard RJ slopes when the CMB contribution is taken into account. High resolution ALMA imaging shows compact morphology and evidence for mergers. This work reveals a population of cold dusty star-forming galaxies that were under-represented in current surveys, and are even colder than typical Main Sequence galaxies at the same redshift. High FIR dust optical depth might be a widespread feature of compact starbursts at any redshift., Comment: 13 pages, 8 figures, 4 tables, accepted for publication in ApJ

We investigate a sample of 152 dusty sources at 1:5 < z < 2:5 to understand the connection of enhanced Star-Formation-Rate (SFR) and Black-Hole-Accretion-Rate (BHAR). The sources are Herschel-selected, having stellar masses M*>10^10Msun and SFR (100-1000Msun/yr) elevated(>4?) above the star-forming "main sequence", classifying them as Starbursts (SB). Through a multiwavelength fitting approach (including a dusty torus component), we divided the sample into active SBs (dominated by an AGN emission, SBs-AGN, ? 23% of the sample) and purely star-forming SBs (SBs-SFR). We visually inspected their HST/UV-restframe maps: SBs-SFR are generally irregular and composite systems; ? 50% of SBs-AGN are instead dominated by regular compact morphologies. We then found archival ALMA continuum counterparts for 33 galaxies (12 SBs-AGN and 21 SBs-SFR). For these sources we computed dust masses, and, with standard assumptions, we also guessed total molecular gas-masses. SBs turn to be gas rich systems (fgas = Mgas=Mgas/(Mgas+M*) \sim 20%-70%), and the gas fractions of the two SB classes are very similar (fgas = 43 +/-4% and fgas = 42+/- 2%). Our results show that SBs are consistent with a mixture of: 1) highly star-forming merging systems (dominating the SBs-SFR), and 2) primordial galaxies, rapidly growing their M* together with their Black Hole (mainly the more compact SBs-AGN). Anyway, feedback effects have not reduced their fgas yet. Indeed, SBs at z = 2, with relatively low bolometric AGN luminosities in the range 10^44 < Lbol(AGN) < 10^46 erg/s (compared to bright optical and X-ray quasars), are still relatively far from the epoch when the AGN feedback will quench the SFR in the host and will substantially depress the gas fractions., Comment: Accepted for publication in ApJ Letters

ALMA measurements for 93 $Herschel$-selected galaxies at $1.1 \leqslant z \leqslant 1.7$ in COSMOS reveal a sizable ($>29$\%) population with compact star formation (SF) sizes, lying on average $> \times 3.6$ below the optical stellar mass ($M_{\star}$)-size relation of disks. This sample widely spans the star-forming Main Sequence (MS), having $10^{8} \leqslant M_{\star} \leqslant 10^{11.5} \ M_{\odot}$ and $20 \leqslant SFR \leqslant 680 \ M_{\odot} \rm yr^{-1}$. The 32 size measurements and 61 upper limits are measured on ALMA images that combine observations of CO(5-4), CO(4-3), CO(2-1) and $\lambda_{\rm obs} \sim 1.1-1.3 \ \rm mm$ continuum, all tracing the star-forming molecular gas. These compact galaxies have instead normally extended $K_{band}$ sizes, suggesting strong specific $SFR$ gradients. Compact galaxies comprise the $50\pm18 \%$ of MS galaxies at $M_{\star} > 10^{11} M_{\odot}$. This is not expected in standard bi-modal scenarios where MS galaxies are mostly steadily-growing extended disks. We suggest that compact MS objects are early post-starburst galaxies in which the merger-driven boost of SF has subsided. They retain their compact SF size until either further gas accretion restores pre-merger galaxy-wide SF, or until becoming quenched. The fraction of merger-affected SF inside the MS seems thus larger than anticipated and might reach $\sim 50$\% at the highest $M_{\star}$. The presence of large galaxies above the MS demonstrates an overall poor correlation between galaxy SF size and specific $SFR$., Comment: 7 pages, 5 figures. Accepted for publication in ApJ Letters

We infer and compare the specific X-ray luminosity distributions for a sample of massive (i.e. $\log_{10} (M*/M\odot) > 10.5$) galaxies split according to their far-infrared-derived star-forming properties (i.e., starburst and non-starburst) and redshift. We model each distribution as a power-law with an upper and lower turnover, and adopt a maximum likelihood method to include information from non-detections in the form of upper limits. When we use our inferred distributions to calculate the ratios of high to low sLx AGN (corresponding to above and below $0.1\lambda_{\text{Edd}}$, respectively) we find that starbursts have significantly higher proportions of high sLx AGN compared to their non-starburst counterparts. These findings help explain the increase in average X-ray luminosity in bins of increasing SFR reported by previous studies., Comment: 14 pages, 10 figures, accepted for publication in MNRAS

We present high-fidelity, 30 milliarcsecond (200-pc) resolution ALMA rest-frame 240 $\mu$m observations of cold dust emission in three typical main-sequence star-forming galaxies (SFGs) at $z \sim 3$ in the Hubble Ultra-Deep Field (HUDF). The cold dust is distributed within the smooth disk-like central regions of star formation $1 - 3$ kpc in diameter, despite their complex and disturbed rest-frame UV and optical morphologies. No dust substructures or clumps are seen down to $\simeq 1- 3$ $M_\odot$yr$^{-1}$ (1$\sigma$) per 200-pc beam. No dust emission is observed at the locations of UV-emitting clumps, which lie $\simeq 2-10$ kpc from the bulk of star formation. Clumpy substructures can contribute no more than $1-7$% of the total star formation in these galaxies (3$\sigma$ upper limits). The lack of star-forming substructures in our HUDF galaxies is to be contrasted with the multiple substructures characteristic of submillimeter-selected galaxies (SMGs) at the same cosmic epoch, particularly the far-IR-bright SMGs with similarly high-fidelity ALMA observations of Hodge et al. (2019). Individual star-forming substructures in these SMGs contain $\sim10-30$% of their total star formation. A substructure in these SMGs is often comparably bright in the far-infrared as (or in some cases brighter than) our typical SFGs, suggesting that these SMGs originate from a class of disruptive event involving multiple objects at the scale of our HUDF galaxies. The scale of the disruptive event found in our main-sequence SFGs, characterized by the lack of star-forming substructures at our resolution and sensitivity, could be less violent, e.g., gas-rich disk instability or minor mergers., Comment: 22 pages, 7 figures; ApJ, accepted

We study a population of significantly sub-solar enrichment galaxies at z=1.99, to investigate how molecular gas, dust and star-formation relate in low-metallicity galaxies at the peak epoch of star-formation. We target our sample with several deep ALMA and VLA datasets, and find no individual detections of CO[4-3], CO[1-0] or dust, in stark contrast to the $>$60% detection rate expected for solar-enrichment galaxies with these MS Halpha SFRs. We find that both low and high density molecular gas (traced by CO[1-0] and CO[4-3] respectively) are affected by the low enrichment, showing sample average (stacked) luminosity deficits $>$0.5-0.7 dex below expectations. This is particularly pertinent for the use of high-J CO emission as a proxy of instantaneous star-formation rate. Our individual galaxy data and stacked constraints point to a strong inverse dependence $\propto Z^{\gamma}$ of gas-to-dust ratios (G/D) and CO-to-H$_{2}$ conversion factors (aco) on metallicity at z$\sim$2, with $\gamma_{\rm G/D}<$-2.2 and $\gamma_{\alpha_{\rm CO}}<$-0.8, respectively. We quantify the importance of comparing G/D and aco vs. metallicity trends from the literature on a common, suitably normalised metallicity scale. When accounting for systematic offsets between different metallicity scales, our z$\sim$2 constraints on these scaling relations are consistent with the corresponding relations for local galaxies. However, among those local relations, we favour those with a steep/double power-law dependence of G/D on metallicity. Finally, we discuss the implications of these findings for (a) gas mass measurements for sub-M$^{*}$ galaxies, and (b) efforts to identify the characteristic galaxy mass scale contributing most to the comoving molecular gas density at z=2., Comment: 15 pages, 5 figures, 3 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We explore the distribution of RMS=SFR/SFR_MS (where SFR_MS is the star formation rate of "Main Sequence" star-forming galaxies) for AGN hosts at z=1. We split our sample into two bins of X-ray luminosity divided at Lx=2x10^43erg s-1 to investigate whether the RMS distribution changes as a function of AGN power. Our main results suggest that, when the RMS distribution of AGN hosts is modelled as a log-normal distribution (i.e. the same shape as that of MS galaxies), galaxies hosting more powerful X-ray AGNs (i.e. Lx>2x10^43erg s-1) display a narrower RMS distribution that is shifted to higher values compared to their lower Lx counterparts. In addition, we find that more powerful X-ray AGNs have SFRs that are more consistent with that of MS galaxies compared to lower Lx AGNs. Despite this, the mean SFRs (as opposed to RMS) measured from these distributions are consistent with the previously observed flat relationship between SFR and Lx. Our results suggest that the typical star-forming properties of AGN hosts change with Lx , and that more powerful AGNs typically reside in more MS-like star-forming galaxies compared to lower Lx AGNs., Comment: 6 pages, 4 figures, Accepted for publication in MNRAS Letter

We present the first results of an ALMA survey of the lower fine structure line of atomic carbon [C I]$(^3P_1\,-\,^{3}P_0)$ in far infrared-selected galaxies on the main sequence at $z\sim1.2$ in the COSMOS field. We compare our sample with a comprehensive compilation of data available in the literature for local and high-redshift starbursting systems and quasars. We show that the [C I]($^3P_1$$\rightarrow$$^3P_0$) luminosity correlates on global scales with the infrared luminosity $L_{\rm IR}$ similarly to low-$J$ CO transitions. We report a systematic variation of $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$/$L_{\rm IR}$ as a function of the galaxy type, with the ratio being larger for main-sequence galaxies than for starbursts and sub-millimeter galaxies at fixed $L_{\rm IR}$. The $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$/$L'_{\rm CO(2-1)}$ and $M_{\rm{[C I]}}$/$M_{\rm dust}$ mass ratios are similar for main-sequence galaxies and for local and high-redshift starbursts within a 0.2 dex intrinsic scatter, suggesting that [C I] is a good tracer of molecular gas mass as CO and dust. We derive a fraction of $f_{\rm{[C\,I]}} = M_{\rm{[C\,I]}} / M_{\rm{C}}\sim3-13$% of the total carbon mass in the atomic neutral phase. Moreover, we estimate the neutral atomic carbon abundance, the fundamental ingredient to calibrate [C I] as a gas tracer, by comparing $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$ and available gas masses from CO lines and dust emission. We find lower [C I] abundances in main-sequence galaxies than in starbursting systems and sub-millimeter galaxies, as a consequence of the canonical $\alpha_{\rm CO}$ and gas-to-dust conversion factors. This argues against the application to different galaxy populations of a universal standard [C I] abundance derived from highly biased samples., Comment: 14 pages + Appendix. Accepted for publication in ApJ. All the data tables in Appendix will be also released in electronic format

We present near-IR spectroscopy in J- and H-band for a large sample of 243 X-ray selected moderate-luminosity type-1 AGN in the COSMOS, SXDS and E-CDF-S survey fields using the multi-object spectrograph Subaru/FMOS. Our sample covers the redshift range 0.50.4 dex). We extend a recently proposed correction based on the CIV blueshift to lower luminosities and black hole masses. While our sample shows an improvement in their CIV black hole mass estimates, the deficit of high blueshift sources reduces its overall importance for moderate-luminosity AGN, compared to the most luminous quasars. In addition, we revisit luminosity correlations between $L_{\rm{bol}}$, $L_X$, $L_{\rm{[OIII]}}$, $L_{5100}$ and $L_{\rm{H}\alpha}$ and find them to be consistent with a simple empirical model, based on a small number of well-established scaling relations. We finally highlight our highest redshift AGN, CID 781, at z=4.6 which shows the lowest black hole mass ($\sim10^8$M$_\odot$) among current near-IR samples at this redshift, and is in a state of fast growth., Comment: 25 pages, 15 figures, accepted for publication in ApJS