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.

We present the detection of CO(5-4) with S/N> 7 - 13 and a lower CO transition with S/N > 3 (CO(4-3) for 4 galaxies, and CO(3-2) for one) with ALMA in band 3 and 4 in five main sequence star-forming galaxies with stellar masses 3-6x10^10 M/M_sun at 3 < z < 3.5. We find a good correlation between the total far-infrared luminosity LFIR and the luminosity of the CO(5-4) transition L'CO(5-4), where L'CO(5-4) increases with SFR, indicating that CO(5-4) is a good tracer of the obscured SFR in these galaxies. The two galaxies that lie closer to the star-forming main sequence have CO SLED slopes that are comparable to other star-forming populations, such as local SMGs and BzK star-forming galaxies; the three objects with higher specific star formation rates (sSFR) have far steeper CO SLEDs, which possibly indicates a more concentrated episode of star formation. By exploiting the CO SLED slopes to extrapolate the luminosity of the CO(1-0) transition, and using a classical conversion factor for main sequence galaxies of alpha_CO = 3.8 M_sun(K km s^-1 pc^-2)^-1, we find that these galaxies are very gas rich, with molecular gas fractions between 60 and 80%, and quite long depletion times, between 0.2 and 1 Gyr. Finally, we obtain dynamical masses that are comparable with the sum of stellar and gas mass (at least for four out of five galaxies), allowing us to put a first constraint on the alpha_CO parameter for main sequence galaxies at an unprecedented redshift.

We combine observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) to assess the redshift and to study the star formation conditions in AzTEC2: one of the brightest sub-millimeter galaxies (SMGs) in the COSMOS field ($S_{\rm 1.1mm}=10.5\pm1.4$mJy). Our high-resolution observations confirm that AzTEC2 splits into two components (namely AzTEC2-A and AzTEC2-B) for which we detect [C$\,$II] and $^{12}$CO(5$\to$4) line emission, implying a redshift of $4.626\pm0.001$ ($4.633\pm0.001$) for AzTEC2-A (AzTEC2-B) and ruling out previous associations with a galaxy at $z\sim1$. We use the $^{12}$CO(5$\to$4) line emission and adopt typical SMG-like gas excitation conditions to estimate the molecular gas mass, which is $M_{\rm gas}(\alpha_{\rm CO}/2.5)=2.1\pm0.4 \times10^{11}{\rm M}_\odot$ for AzTEC2-A, and a factor four lower for AzTEC2-B. With the infrared-derived star formation rate of AzTEC2-A ($1920\pm100 \,M_\odot{\rm \, yr}^{-1}$) and AzTEC2-B ($710\pm 35\,M_\odot{\rm \,yr}^{-1}$), they both will consume their current gas reservoir within $(30-200)\,$Myr. We find evidence of a rotation-dominated [C$\,$II] disk in AzTEC2-A, with a de-projected rotational velocity of $v_{\rm rot}(i=39^\circ)=660\pm130{\rm \,km\,s}^{-1}$, velocity dispersion $\lesssim100{\rm \,km\,s}^{-1}$, and dynamical mass of $M_{\rm dyn}(i=39^\circ)=2.6^{+1.2}_{-0.9}\times10^{11}\,M_\odot$. We propose that an elevated gas accretion rate from the cosmic web might be the main driver of the intense levels of star formation in AzTEC2-A, which might be further enhanced by gravitational torques induced by its minor companion (AzTEC2-B). These results strengthen the picture whereby the population of single-dish selected SMGs is rather heterogeneous, including a population of pairs of massive, highly-active galaxies in a pre-coalescence phase.

The rich information on (sub)millimeter dust continuum emission from distant galaxies in the public Atacama Large Millimeter/submillimeter Array (ALMA) archive is contained in thousands of inhomogeneous observations from individual PI-led programs. To increase the usability of these data for studies deepening our understanding of galaxy evolution, we have developed automated mining pipelines for the ALMA archive in the COSMOS field (A3COSMOS) that efficiently exploit the available information for large numbers of galaxies across cosmic time and keep the data products in sync with the increasing public ALMA archive: (a) a dedicated ALMA continuum imaging pipeline, (b) two complementary photometry pipelines for both blind source extraction and prior source fitting, (c) a counterpart association pipeline utilizing the multiwavelength data available (including quality assessment based on machine-learning techniques), (d) an assessment of potential (sub)millimeter line contribution to the measured ALMA continuum, and (e) extensive simulations to provide statistical corrections to biases and uncertainties in the ALMA continuum measurements. Application of these tools yields photometry catalogs with ~1000 (sub)millimeter detections (spurious fraction ~8%-12%) from over 1500 individual ALMA continuum images. Combined with ancillary photometric and redshift catalogs and the above quality assessments, we provide robust information on redshift, stellar mass, and star formation rate for ~700 galaxies at redshifts 0.5-6 in the COSMOS field (with undetermined selection function). The ALMA photometric measurements and galaxy properties are released publicly within our blind extraction, prior fitting, and galaxy property catalogs, plus the images. These products will be updated on a regular basis in the future.

We identify multi-wavelength counterparts to 1,147 submillimeter sources from the S2COSMOS SCUBA-2 survey of the COSMOS field by employing a recently developed radio$+$machine-learning method trained on a large sample of ALMA-identified submillimeter galaxies (SMGs), including 260 SMGs identified in the AS2COSMOS pilot survey. In total, we identify 1,222 optical/near-infrared(NIR)/radio counterparts to the 897 S2COSMOS submillimeter sources with S$_{850}$>1.6mJy, yielding an overall identification rate of ($78\pm9$)%. We find that ($22\pm5$)% of S2COSMOS sources have multiple identified counterparts. We estimate that roughly 27% of these multiple counterparts within the same SCUBA-2 error circles very likely arise from physically associated galaxies rather than line-of-sight projections by chance. The photometric redshift of our radio$+$machine-learning identified SMGs ranges from z=0.2 to 5.7 and peaks at $z=2.3\pm0.1$. The AGN fraction of our sample is ($19\pm4$)%, which is consistent with that of ALMA SMGs in the literature. Comparing with radio/NIR-detected field galaxy population in the COSMOS field, our radio+machine-learning identified counterparts of SMGs have the highest star-formation rates and stellar masses. These characteristics suggest that our identified counterparts of S2COSMOS sources are a representative sample of SMGs at z<3. We employ our machine-learning technique to the whole COSMOS field and identified 6,877 potential SMGs, most of which are expected to have submillimeter emission fainter than the confusion limit of our S2COSMOS surveys (S$_{850}$<1.5mJy). We study the clustering properties of SMGs based on this statistically large sample, finding that they reside in high-mass dark matter halos ($(1.2\pm0.3)\times10^{13}\,h^{-1}\,\rm M_{\odot}$), which suggests that SMGs may be the progenitors of massive ellipticals we see in the local Universe.

We present deep continuum observations at a wavelength of 2mm centered on the COSMOS field using the Goddard IRAM Superconducting Millimeter Observer (GISMO) at the IRAM 30m-telescope. These data constitute the widest deep 2mm survey to-date, reaching a uniform $\sigma\sim0.23$mJy beam$^{-1}$ sensitivity over $\sim250$arcmin$^2$ at $\sim\,24''$ resolution. We detect four sources at high significance (S/N$\,\geq\,$4.4) with an expected number of false detection of 0.09 sources, and five sources at $4.4>\,$S/N$\,\geq\,3.7$ with an expected number of false detection of 1.65 sources. Combined with deep GISMO observations in GOODS-N, we constrain the 2mm number counts over one decade in flux density. These measurements agree with most galaxy evolution models tested here, except those with large population of dusty star-forming galaxies at $z>7$. Five GISMO sources have counterparts in (sub-)millimeter catalogs available in COSMOS. Their redshifts suggest that all but one lie above $z\sim3$. These four high-redshift ($z>3$) galaxies have $\tilde{z}=3.9$, SFRs $\sim\,$400 - 1200M$_{\odot}\,$yr$^{-1}$ and $M_{\rm dust}\sim10^{9.5}\,$M$_{\odot}$. They provide a relatively complete selection ($\sim66\%$) of the most luminous ($L_{\rm IR}>10^{12.6}\,$L$_{\odot}$) and highest redshift ($z>3$) galaxies detected within our survey area by AzTEC at 1.1mm. We thus conclude that 2mm surveys favor the selection of massive, vigorously star-forming, high-redshift galaxies. This is corroborated by GISMO-C4, a source with a low false detection probability ($\sim\,6.2\%$), for which the absence of a (sub-)millimeter counterpart supports a high redshift origin ($z\gtrsim3$).

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\,\mu Jy \,beam^{-1}$, $\theta_{\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)^\alpha$ yields a shallow slope of only $\alpha=-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 $(\Sigma_{\rm SFR})$ scales with the distance to the MS, except for a fraction of MS galaxies ($\lesssim10\%$) that harbor starburst-like $\Sigma_{\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 $\Sigma_{\rm SFR}$ and distance to the MS to better identify the galaxy population undergoing a starbursting phase.

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 {\mu}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]

We present 0.15-arcsec (1 kpc) resolution ALMA observations of the [CII] 157.74 um line and rest-frame 160-um continuum emission in two z~3 dusty, star-forming galaxies - ALESS 49.1 and ALESS 57.1, combined with resolved CO(3-2) observations. In both sources, the [CII] surface brightness distribution is dominated by a compact core $\leq$1 kpc in radius, a factor of 2-3 smaller than the extent of the CO(3-2) emission. In ALESS 49.1, we find an additional extended (8-kpc radius), low surface-brightness [CII] component. Based on an analysis of mock ALMA observations, the [CII] and 160-um continuum surface brightness distributions are inconsistent with a single-Gaussian surface brightness distribution with the same size as the CO(3-2) emission. The [CII] rotation curves flatten at $\simeq$2 kpc radius, suggesting the kinematics of the central regions are dominated by a baryonic disc. Both galaxies exhibit a strong [CII]/FIR deficit on 1-kpc scales, with FIR-surface-brightness to [CII]/FIR slope steeper than in local star-forming galaxies. A comparison of the [CII]/CO(3-2) observations with PDR models suggests a strong FUV radiation field ($G_0\sim10^4$) and high gas density ($n\mathrm{(H)}\sim10^4-10^5$ cm$^{-3}$) in the central regions of ALESS 49.1 and 57.1. The most direct interpretation of the pronounced [CII]/FIR deficit is a thermal saturation of the C+ fine-structure levels at temperatures $\geq$500 K, driven by the strong FUV field.

Theoretical models of galaxy formation suggest that the presence of an active galactic nucleus (AGN) is required to regulate the growth of its host galaxy through feedback mechanisms, produced by e.g. AGN-driven outflows. Although such outflows are common both at low and high redshift, a comprehensive picture is still missing. The peak epoch of galaxy assembly (1<z<3) has been poorly explored so far, and current observations in this redshift range are mostly limited to targets with high chances to be in an outflowing phase. This paper introduces SUPER (a SINFONI Survey for Unveiling the Physics and Effect of Radiative feedback), an ongoing ESO's VLT/SINFONI Large Programme. SUPER will perform the first systematic investigation of ionized outflows in a sizeable and blindly-selected sample of 39 X-ray AGN at z~2, which reaches high spatial resolutions (~2 kpc) thanks to the adaptive optics-assisted IFU observations. The outflow morphology and star formation in the host galaxy will be mapped through the broad component of [OIII] and the narrow component of Ha emission lines. The main aim of our survey is to infer the impact of outflows on the on-going star formation and to link the outflow properties to a number of AGN and host galaxy properties. We describe here the survey characteristics and goals, as well as the selection of the target sample. Moreover, we present a full characterization of its multi-wavelength properties: we measure, via spectral energy distribution fitting of UV-to-FIR photometry, stellar masses (4x10^9-2x10^11 Msun), star formation rates (25-680 Msun yr^-1) and AGN bolometric luminosities (2x10^44-8x10^47 erg s^-1), along with obscuring column densities (<2x10^24 cm^-2) and 2-10 keV luminosities (2x10^43-6x10^45 erg s^-1) derived through X-ray spectral analysis. Finally, we classify our AGN as jetted or non-jetted according to their radio and FIR emission.

We use ALMA observations of four sub-millimetre galaxies (SMGs) at $z\sim2-3$ to investigate the spatially resolved properties of the inter-stellar medium (ISM) at scales of 1--5 kpc (0.1--0.6$''$). The velocity fields of our sources, traced by the $^{12}$CO($J$=3-2) emission, are consistent with disk rotation to first order, implying average dynamical masses of $\sim$3$\times10^{11}$M$_{\odot}$ within two half-light radii. Through a Bayesian approach we investigate the uncertainties inherent to dynamically constraining total gas masses. We explore the covariance between the stellar mass-to-light ratio and CO-to-H$_{2}$ conversion factor, $\alpha_{\rm CO}$, finding values of $\alpha_{\rm CO}=1.1^{+0.8}_{-0.7}$ for dark matter fractions of 15 \%. We show that the resolved spatial distribution of the gas and dust continuum can be uncorrelated to the stellar emission, challenging energy balance assumptions in global SED fitting. Through a stacking analysis of the resolved radial profiles of the CO(3-2), stellar and dust continuum emission in SMG samples, we find that the cool molecular gas emission in these sources (radii $\sim$5--14 kpc) is clearly more extended than the rest-frame $\sim$250 $\mu$m dust continuum by a factor $>2$. We propose that assuming a constant dust-to-gas ratio, this apparent difference in sizes can be explained by temperature and optical-depth gradients alone. Our results suggest that caution must be exercised when extrapolating morphological properties of dust continuum observations to conclusions about the molecular gas phase of the ISM.

We present deep high resolution (0.03", 200pc) ALMA Band 7 observations covering the dust continuum and [CII] $\lambda157.7\mu$m emission in four $z\sim4.4-4.8$ sub-millimeter galaxies (SMGs) selected from the ALESS and AS2UDS surveys. The data show that the rest-frame 160$\mu$m (observed 345 GHz) dust emission is consistent with smooth morphologies on kpc scales for three of the sources. One source, UDS47.0, displays apparent substructure but this is also consistent with a smooth morphology, as indicated by simulations showing that smooth exponential disks can appear clumpy when observed at high angular resolution (0.03") and depth of these observations ($\sigma_{345\text{GHz}} \sim27-47\mu$Jy beam$^{-1}$). The four SMGs are bright [CII] emitters, and we extract [CII] spectra from the high resolution data, and recover $\sim20-100$% of the [CII] flux and $\sim40-80$% of the dust continuum emission, compared to the previous lower resolution observations. When tapered to 0.2" resolution our maps recover $\sim80-100$% of the continuum emission, indicating that $\sim60$% of the emission is resolved out on $\sim200$pc scales. We find that the [CII] emission in high-redshift galaxies is more spatially extended than the rest-frame 160$\mu$m dust continuum by a factor of $1.6\pm0.4$. By considering the $L_{\text{[CII]}}$/$L_{\text{FIR}}$ ratio as a function of the star-formation rate surface density ($\Sigma_{\text{SFR}}$) we revisit the [CII] deficit, and suggest that the decline in the $L_{\text{[CII]}}$/$L_{\text{FIR}}$ ratio as a function of $\Sigma_{\text{SFR}}$ is consistent with local processes. We also explore the physical drivers that may be responsible for these trends and can give rise to the properties found in the densest regions of SMGs.

Dust-enshrouded, starbursting, submillimeter galaxies (SMGs) at $z \geq 3$ have been proposed as progenitors of $z \geq 2$ compact quiescent galaxies (cQGs). To test this connection, we present a detailed spatially resolved study of the stars, dust and stellar mass in a sample of six submillimeter-bright starburst galaxies at $z \sim 4.5$. The stellar UV emission probed by HST is extended, irregular and shows evidence of multiple components. Informed by HST, we deblend Spitzer/IRAC data at rest-frame optical finding that the systems are undergoing minor mergers, with a typical stellar mass ratio of 1:6.5. The FIR dust continuum emission traced by ALMA locates the bulk of star formation in extremely compact regions (median $r_{\rm{e}} = 0.70 \pm 0.29$ kpc) and it is in all cases associated with the most massive component of the mergers (median $\log (M_{*}/M_{\odot}) = 10.49 \pm 0.32$). We compare spatially resolved UV slope ($\beta$) maps with the FIR dust continuum to study the infrared excess ($\rm{IRX} = L_{\rm{IR}}/L_{\rm{UV}}$)-$\beta$ relation. The SMGs display systematically higher $\rm{IRX}$ values than expected from the nominal trend, demonstrating that the FIR and UV emissions are spatially disconnected. Finally, we show that the SMGs fall on the mass-size plane at smaller stellar masses and sizes than cQGs at $z = 2$. Taking into account the expected evolution in stellar mass and size between $z = 4.5$ and $z = 2$ due to the ongoing starburst and mergers with minor companions, this is in agreement with a direct evolutionary connection between the two populations.

We studied the molecular gas properties of AzTEC/C159, a star-forming disk galaxy at $z=4.567$. We secured $^{12}$CO molecular line detections for the $J=2\to1$ and $J=5\to4$ transitions using the Karl G. Jansky VLA and the NOEMA interferometer. The broad (FWHM$\sim750\,{\rm km\,s}^{-1}$) and tentative double-peaked profiles of both $^{12}$CO lines are consistent with an extended molecular gas reservoir, which is distributed in a rotating disk as previously revealed from [CII] 158$\mu$m line observations. Based on the $^{12}$CO(2$\to$1) emission line we derived $L'_{\rm{CO}}=(3.4\pm0.6)\times10^{10}{\rm \,K\,km\,s}^{-1}{\rm \,pc}^{2}$, that yields a molecular gas mass of $M_{\rm H_2 }(\alpha_{\rm CO}/4.3)=(1.5\pm0.3)\times 10^{11}{\rm M}_\odot$ and unveils a gas-rich system with $\mu_{\rm gas}(\alpha_{\rm CO}/4.3)\equiv M_{\rm H_2}/M_\star=3.3\pm0.7$. The extreme star formation efficiency (SFE) of AzTEC/C159, parametrized by the ratio $L_{\rm{IR}}/L'_{\rm{CO}}=(216\pm80)\, {\rm L}_{\odot}{\rm \,(K\,km\,s}^{-1}{\rm \,pc}^{2})^{-1}$, is comparable to merger-driven starbursts such as local ultra-luminous infrared galaxies (ULIRGs) and SMGs. Likewise, the $^{12}$CO(5$\to$4)/CO(2$\to$1) line brightness temperature ratio of $r_{52}= 0.55\pm 0.15$ is consistent with high excitation conditions, similar to that observed in SMGs. We constrained the value for the $L'_{\text{CO}}-{\rm H}_2$ mass conversion factor in AzTEC/C159, i.e. $\alpha_{\text{CO}}=3.9^{+2.7}_{-1.3}{\rm \,M}_{\odot}{\rm \,K}^{-1}{\rm \,km}^{-1}{\rm \,s\,pc}^{-2}$, that is consistent with a self-gravitating molecular gas distribution as observed in local star-forming disk galaxies. Cold gas streams from cosmological filaments might be fueling a gravitationally unstable gas-rich disk in AzTEC/C159, which breaks into giant clumps forming stars as efficiently as in merger-driven systems and generate high gas excitation.

We use a sample of 4178 Lyman break galaxies (LBGs) at z = 3, 4 and 5 in the UKIRT Infrared Deep Sky Survey (UKIDSS) Ultra Deep Survey (UDS) field to investigate the relationship between the observed slope of the stellar continuum emission in the ultraviolet, {\beta}, and the thermal dust emission, as quantified via the so-called 'infrared excess' (IRX = LIR/LUV). Through a stacking analysis we directly measure the 850-{\mu}m flux density of LBGs in our deep (0.9mJy) James Clerk Maxwell Telescope (JCMT) SCUBA-2 850-{\mu}m map, as well as deep public Herschel/SPIRE 250-, 350- and 500-{\mu}m imaging. We establish functional forms for the IRX-{\beta} relation to z ~ 5, confirming that there is no significant redshift evolution of the relation and that the resulting average IRX-{\beta} curve is consistent with a Calzetti-like attenuation law. We compare our results with recent work in the literature, finding that discrepancies in the slope of the IRX-{\beta} relation are driven by biases in the methodology used to determine the ultraviolet slopes. Consistent results are found when IRX-{\beta} is evaluated by stacking in bins of stellar mass, M, and we argue that the near-linear IRX-M relationship is a better proxy for correcting observed UV luminosities to total star formation rates, provided an accurate handle on M can be had, and also gives clues as to the physical driver of the role of dust-obscured star formation in high-redshift galaxies.

We report the detection of CO(1 - 0) line emission from seven Planck and Herschel selected hyper luminous (LIR(8-1000um) > 10^13Lsun) infrared galaxies with the Green Bank Telescope (GBT). CO(1 - 0) measurements are a vital tool to trace the bulk molecular gas mass across all redshifts. Our results place tight constraints on the total gas content of these most apparently luminous high-z star-forming galaxies (apparent IR luminosities of LIR > 10^(13-14) Lsun), while we confirm their predetermined redshifts measured using the Large Millimeter Telescope, LMT (zCO = 1.33 - 3.26). The CO(1 - 0) lines show similar profiles as compared to Jup = 2 -4 transitions previously observed with the LMT. We report enhanced infrared to CO line luminosity ratios of <LIR/L'CO(1-0)> = 110 (pm 22) Lsun(K km s^-1 pc^-2)^-1 compared to normal star-forming galaxies, yet similar to those of well-studied IR-luminous galaxies at high-z. We find average brightness temperature ratios of <r21> = 0.93 (2 sources), <r31> = 0.34 (5 sources), and <r41> = 0.18 (1 source). The r31 and r41 values are roughly half the average values for SMGs. We estimate the total gas mass content as uMH2 = (0.9 - 27.2) x 10^11(alphaCO/0.8)Msun, where u is the magnification factor and alphaCO is the CO line luminosity to molecular hydrogen gas mass conversion factor. The rapid gas depletion times are, on average, tau = 80 Myr, which reveal vigorous starburst activity, and contrast the Gyr depletion timescales observed in local, normal star-forming galaxies.

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.

Recent literature suggests that there are two modes through which galaxies grow their stellar mass - a normal mode characterized by quasi-steady star formation, and a highly efficient starburst mode possibly triggered by stochastic events such as galaxy mergers. While these differences are established for extreme cases, the population of galaxies in-between these two regimes is poorly studied and it is not clear where the transition between these two modes of star formation occurs. We utilize ALMA observations of the CO J=3-2 line luminosity in a sample of 20 infrared luminous galaxies that lie in the intermediate range between normal and starburst galaxies at z ~ 0.25-0.6 in the COSMOS field to examine the gas content and star formation efficiency of these galaxies. We compare these quantities to the galaxies' deviation from the well-studied "main sequence" correlation between star formation rate and stellar mass (MS) and find that at log($SFR/SFR_{MS}$) < 0.6, a galaxy's distance to the main sequence is mostly driven by increased gas content, and not a more efficient star formation process.

Until recently, determining the rotational properties of galaxies in the early universe (z>4, Universe age <1.5Gyr) was impractical, with the exception of a few strongly lensed systems. Combining the high resolution and sensitivity of ALMA at (sub-) millimeter wavelengths with the typically high strength of the [CII] 158um emission line from galaxies and long-developed dynamical modeling tools raises the possibility of characterizing the gas dynamics in both extreme starburst galaxies and normal star forming disk galaxies at z~4-7. Using a procedure centered around GIPSY's ROTCUR task, we have fit tilted ring models to some of the best available ALMA [CII] data of a small set of galaxies: the MS galaxies HZ9 & HZ10, the Damped Lyman-alpha Absorber (DLA) host galaxy ALMA J0817+1351, the submm galaxies AzTEC/C159 and COSMOS J1000+0234, and the quasar host galaxy ULAS J1319+0950. This procedure directly derives rotation curves and dynamical masses as functions of radius for each object. In one case, we present evidence for a dark matter halo of O(10^11) solar masses. We present an analysis of the possible velocity dispersions of AzTEC/C159 and ULAS J1319+0950 based on matching simulated observations to the integrated [CII] line profiles. Finally, we test the effects of observation resolution and sensitivity on our results. While the conclusions remain limited at the resolution and signal-to-noise ratios of these observations, the results demonstrate the viability of the modeling tools at high redshift, and the exciting potential for detailed dynamical analysis of the earliest galaxies, as ALMA achieves full observational capabilities.

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 $\mu$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.

We present detailed studies of a $z=2.12$ submillimeter galaxy, ALESS67.1, using sub-arcsecond resolution ALMA, AO-aided VLT/SINFONI, and HST/CANDELS data to investigate the kinematics and spatial distributions of dust emission (870 $\mu$m continuum), $^{12}$CO($J$=3-2), strong optical emission lines, and visible stars. Dynamical modelling of the optical emission lines suggests that ALESS67.1 is not a pure rotating disk but a merger, consistent with the apparent tidal features revealed in the HST imaging. Our sub-arcsecond resolution dataset allow us to measure half-light radii for all the tracers, and we find a factor of 4-6 smaller sizes in dust continuum compared to all the other tracers, including $^{12}$CO, and UV and H$\alpha$ emission is significantly offset from the dust continuum. The spatial mismatch between UV continuum and the cold dust and gas reservoir supports the explanation that geometrical effects are responsible for the offset of dusty galaxy on the IRX-$\beta$ diagram. Using a dynamical method we derive an $\alpha_{\rm CO}=1.8\pm1.0$, consistent with other SMGs that also have resolved CO and dust measurements. Assuming a single $\alpha_{\rm CO}$ value we also derive resolved gas and star-formation rate surface densities, and find that the core region of the galaxy ($\lesssim5$ kpc) follows the trend of mergers on the Schmidt-Kennicutt relationship, whereas the outskirts ($\gtrsim5$ kpc) lie on the locus of normal star-forming galaxies, suggesting different star-formation efficiencies within one galaxy. Our results caution against using single size or morphology for different tracers of the star-formation activity and gas content of galaxies, and therefore argue the need to use spatially-resolved, multi-wavelength observations to interpret the properties of SMGs, and perhaps even for $z>1$ galaxies in general.

Bright Ly-$\alpha$ blobs (LABs) --- extended nebulae with sizes of $\sim$100kpc and Ly-$\alpha$ luminosities of $\sim$10$^{44}$erg s$^{-1}$ --- often reside in overdensities of compact Ly-$\alpha$ emitters (LAEs) that may be galaxy protoclusters. The number density, variance, and internal kinematics of LABs suggest that they themselves trace group-like halos. Here we test this hierarchical picture, presenting deep, wide-field Ly-$\alpha$ narrowband imaging of a 1$^\circ$ $\times$ 0.5$^\circ$ region around a LAB pair at $z$ = 2.3 discovered previously by a blind survey. We find 183 Ly-$\alpha$ emitters, including the original LAB pair and three new LABs with Ly-$\alpha$ luminosities of (0.9--1.3)$\times$10$^{43}$erg s$^{-1}$ and isophotal areas of 16--24 arcsec$^2$. Using the LAEs as tracers and a new kernel density estimation method, we discover a large-scale overdensity (Bo{\"o}tes J1430+3522) with a surface density contrast of $\delta_{\Sigma}$ = 2.7, a volume density contrast of $\delta$ $\sim$ 10.4, and a projected diameter of $\approx$ 20 comoving Mpc. Comparing with cosmological simulations, we conclude that this LAE overdensity will evolve into a present-day Coma-like cluster with $\log{(M/M_\odot)}$ $\sim$ $15.1\pm0.2$. In this and three other wide-field LAE surveys re-analyzed here, the extents and peak amplitudes of the largest LAE overdensities are similar, not increasing with survey size, implying that they were indeed the largest structures then and do evolve into rich clusters today. Intriguingly, LABs favor the outskirts of the densest LAE concentrations, i.e., intermediate LAE overdensities of $\delta_\Sigma = 1 - 2$. We speculate that these LABs mark infalling proto-groups being accreted by the more massive protocluster.

We determine the physical properties of a sample of SMGs in the COSMOS field that were pre-selected at the observed wavelength of $\lambda_{\rm obs}=1.1$ mm, and followed up at $\lambda_{\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 $\lambda_{\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.

The star formation rate per unit area correlates well with the gas surface density for different types of galaxies. However, this Kennicutt-Schmidt law has not yet been examined for a large, homogeneously selected sample of SMGs, which could provide useful SF implementation information for models of massive galaxy formation and evolution. We aim at determining the K-S law parameters for the first time for a well-selected, statistical sample of SMGs. We used ALMA to conduct a 0.2" resolution, 870 $\mu$m imaging survey of 40 SMGs in COSMOS, which were initially selected at 1.1 mm. We analysed a sample of 32/40 target SMGs, for which our new ALMA 870 $\mu$m data provide information about the spatial extent of dust emission, and all of which have dust-obscured SFR and dust-based gas mass estimates available from our previous study. We divided our sample into equally large subsamples of main-sequence objects and starbursts, and found their K-S relations to be of the form $\Sigma_{\rm SFR} \propto \Sigma_{\rm gas}^{0.81\pm0.01}$ and $\Sigma_{\rm SFR} \propto \Sigma_{\rm gas}^{0.84\pm0.39}$, respectively. The slightly sub-linear K-S slopes we derived suggest that the SF efficiency is nearly constant across the $\Sigma_{\rm gas}$ range probed. Under the assumption of a Galactic CO-to-H$_2$ conversion factor for the whole sample, the MS SMGs obey a constant global SFE of about 21% per 100 Myr, while that of starburst SMGs is about 27%. The corresponding gas depletion times are $\sim480$ Myr and 370 Myr. On average, our SMGs have $\Sigma_{\rm gas}\gtrsim10^{3.9}$ M$_{\odot}$ pc$^{-2}$, which suggests that they are Eddington-limited. This is consistent with the theoretical expectation of a linear K-S relation for such systems. However, size measurements of the CO-emitting regions of SMGs, and the $\alpha_{\rm CO}$ values of SMGs are needed to further constrain their $\Sigma_{\rm gas}$ values.

We present results on the dust attenuation of galaxies at redshift $\sim3-6$ by studying the relationship between the UV spectral slope ($\beta_{\rm UV}$) and the infrared excess (IRX; $L_{\rm IR}$/$L_{\rm UV}$) using ALMA far-infrared continuum observations. Our study is based on a sample of 67 massive, star-forming galaxies with a median mass of $M_{\ast}\sim 10^{10.7}\,M_{\rm \odot}$ spanning a redshift range $z=2.6-3.7$ (median $z=3.2$) that were observed with ALMA at $\lambda_{rest}=300\,{\rm \mu m}$. Both the individual ALMA detections (41 sources) and stacks including all galaxies show the IRX-$\beta_{\rm UV}$ relationship at $z\sim3$ is mostly consistent with that of local starburst galaxies on average. However, we find evidence for a large dispersion around the mean relationship by up to $\pm0.5$ dex. Nevertheless, the locally calibrated dust correction factors based on the IRX-$\beta_{\rm UV}$ relation are on average applicable to main-sequence $z\sim3$ galaxies. This does not appear to be the case at even higher redshifts, however. Using public ALMA observations of $z\sim4-6$ galaxies we find evidence for a significant evolution in the IRX-$\beta_{\rm UV}$ and the IRX-$M_{\ast}$ relations beyond $z\sim3$ toward lower IRX values. We discuss several caveats that could affect these results, including the assumed dust temperature. ALMA observations of larger $z>3$ galaxy samples will be required to confirm this intriguing redshift evolution.

We present spectroscopic redshifts of S(870)>2mJy submillimetre galaxies (SMGs) which have been identified from the ALMA follow-up observations of 870um detected sources in the Extended Chandra Deep Field South (the ALMA-LESS survey). We derive spectroscopic redshifts for 52 SMGs, with a median of z=2.4+/-0.1. However, the distribution features a high redshift tail, with ~25% of the SMGs at z>3. Spectral diagnostics suggest that the SMGs are young starbursts, and the velocity offsets between the nebular emission and UV ISM absorption lines suggest that many are driving winds, with velocity offsets up to 2000km/s. Using the spectroscopic redshifts and the extensive UV-to-radio photometry in this field, we produce optimised spectral energy distributions (SEDs) using Magphys, and use the SEDs to infer a median stellar mass of M*=(6+/-1)x10^{10}Msol for our SMGs with spectroscopic redshifts. By combining these stellar masses with the star-formation rates (measured from the far-infrared SEDs), we show that SMGs (on average) lie a factor ~5 above the main-sequence at z~2. We provide this library of 52 template fits with robust and well-sampled SEDs available as a resource for future studies of SMGs, and also release the spectroscopic catalog of ~2000 (mostly infrared-selected) galaxies targeted as part of the spectroscopic campaign.

We examine the behaviour of the infrared-radio correlation (IRRC) over the range $0<z<6$ using new, highly sensitive 3GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared data from the Herschel Space Observatory in the 2deg$^{2}$ COSMOS field. We distinguish between objects where emission is believed to arise solely from star-formation, and those where an active galactic nucleus (AGN) is thought to be present. We account for non-detections in the radio or in the infrared using a doubly-censored survival analysis. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4GHz radio luminosity ratio ($q_{\rm TIR}$), decreases with increasing redshift: $q_{\rm TIR}(z)=(2.88\pm0.03)(1+z)^{-0.19\pm0.01}$. Moderate-to-high radiative luminosity AGN do not follow the same $q_{\rm TIR}$$(z)$ trend, having a lower normalisation and steeper decrease with redshift. We cannot rule out the possibility that unidentified AGN contributions only to the radio regime may be steepening the observed $q_{\rm TIR}(z)$ trend of the star-forming population. An increasing fractional contribution to the observed 3GHz flux by free-free emission of star-forming galaxies may also affect the derived evolution. However, we find that the standard (M82-based) assumption of the typical radio spectral energy distribution (SED) for star-forming galaxies is inconsistent with our results. This suggests a more complex shape of the typical radio SED for star-forming galaxies, and that imperfect $K$ corrections in the radio may govern the derived redshift trend of $q_{\rm TIR}$. Lastly, we present a redshift-dependent relation between rest-frame 1.4GHz radio luminosity and star formation rate taking the derived redshift trend into account.

(abridged) We study the composition of the faint radio population selected from the VLA-COSMOS 3GHz Large Project. The survey covers a 2.6sq.deg. area with a mean rms of ~2.3uJy/b, cataloging 10830 sources (>5sigma). Combining these radio data with optical, near-infrared (UltraVISTA), mid-infrared (Spitzer/IRAC) data, and X-ray data (Chandra), we find counterparts to radio sources for ~93% of the radio sample (in the areas of the COSMOS field not affected by saturated or bright sources in the optical to NIR bands), reaching out to z<6. We further classify the sources as star forming galaxies or AGN based on various criteria, such as X-ray luminosity, observed MIR color, UV-FIR spectral-energy distribution, rest-frame NUV-optical color corrected for dust extinction, and radio-excess relative to that expected from the hosts' star-formation rate. We separate the AGN into sub-samples dominated by low-to-moderate and moderate-to-high radiative luminosity AGN, candidates for high-redshift analogues to local low- and high-excitation emission line AGN, respectively. We study the fractional contributions of these sub-populations down to radio flux levels of ~11uJy at 3GHz (or ~20uJy at 1.4GHz assuming a spectral index of -0.7). We find that the dominant fraction at 1.4GHz flux densities above ~200uJy is constituted of low-to-moderate radiative luminosity AGN. Below densities of ~100uJy the fraction of star-forming galaxies increases to ~60%, followed by the moderate-to-high radiative luminosity AGN (~20%), and low-to-moderate radiative luminosity AGN (~20%). Based on this observational evidence, we extrapolate the fractions down to sensitivities of the SKA. Our estimates suggest that at the faint flux limits to be reached by the SKA1 surveys, a selection based only on radio flux limits can provide a simple tool to efficiently identify samples highly (>75%) dominated by star-forming galaxies.

We make use of the deep Karl G. Jansky Very Large Array (VLA) COSMOS radio observations at 3 GHz to infer radio luminosity functions of star-forming galaxies up to redshifts of z~5 based on approximately 6000 detections with reliable optical counterparts. This is currently the largest radio-selected sample available out to z~5 across an area of 2 square degrees with a sensitivity of rms=2.3 ujy/beam. By fixing the faint and bright end shape of the radio luminosity function to the local values, we find a strong redshift trend that can be fitted with a pure luminosity evolution L~(1+z)^{(3.16 +- 0.2)-(0.32 +- 0.07) z}. We estimate star formation rates (SFRs) from our radio luminosities using an infrared (IR)-radio correlation that is redshift dependent. By integrating the parametric fits of the evolved luminosity function we calculate the cosmic SFR density (SFRD) history since z~5. Our data suggest that the SFRD history peaks between 2<z<3 and that the ultraluminous infrared galaxies (ULIRGs; 100 Msol/yr<SFR<1000 Msol/yr) contribute up to ~25% to the total SFRD in the same redshift range. Hyperluminous infrared galaxies (HyLIRGs; SFR>1000 Msol/yr) contribute an additional <2% in the entire observed redshift range. We find evidence of a potential underestimation of SFRD based on ultraviolet (UV) rest-frame observations of Lyman break galaxies (LBGs) at high redshifts (z>4) on the order of 15-20%, owing to appreciable star formation in highly dust-obscured galaxies, which might remain undetected in such UV observations.

We explore the multiwavelength properties of AGN host galaxies for different classes of radio-selected AGN out to z$\lesssim$6 via a multiwavelength analysis of about 7700 radio sources in the COSMOS field. The sources were selected with the Very Large Array (VLA) at 3 GHz (10 cm) within the VLA-COSMOS 3 GHz Large Project, and cross-matched with multiwavelength ancillary data. This is the largest sample of high-redshift (z$\lesssim$6) radio sources with exquisite photometric coverage and redshift measurements available. We constructed a sample of moderate-to-high radiative luminosity AGN (HLAGN) via spectral energy distribution (SED) decomposition combined with standard X-ray and mid-infrared diagnostics. Within the remainder of the sample we further identified low-to-moderate radiative luminosity AGN (MLAGN) via excess in radio emission relative to the star formation rates in their host galaxies. We show that AGN power in HLAGN occurs predominantly in radiative form, while MLAGN display a substantial mechanical AGN luminosity component. We found significant differences in the host properties of the two AGN classes, as a function of redshift. At z$<$1.5, MLAGN appear to reside in significantly more massive and less star-forming galaxies compared to HLAGN. At z$>$1.5, we observed a reversal in the behaviour of the stellar mass distributions with the HLAGN populating the higher stellar mass tail. We interpret this finding as a possible hint of the downsizing of galaxies hosting HLAGN, with the most massive galaxies triggering AGN activity earlier than less massive galaxies, and then fading to MLAGN at lower redshifts. Our conclusion is that HLAGN and MLAGN samples trace two distinct galaxy and AGN populations in a wide range of redshifts, possibly resembling the radio AGN types often referred to as radiative- and jet-mode (or high- and low-excitation), respectively.

We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys.

We determine the radio size distribution of a large sample of 152 SMGs in COSMOS that were detected with ALMA at 1.3 mm. For this purpose, we used the observations taken by the VLA-COSMOS 3 GHz Large Project. One hundred and fifteen of the 152 target SMGs were found to have a 3 GHz counterpart. The median value of the major axis FWHM at 3 GHz is derived to be $4.6\pm0.4$ kpc. The radio sizes show no evolutionary trend with redshift, or difference between different galaxy morphologies. We also derived the spectral indices between 1.4 and 3 GHz, and 3 GHz brightness temperatures for the sources, and the median values were found to be $\alpha=-0.67$ and $T_{\rm B}=12.6\pm2$ K. Three of the target SMGs, which are also detected with the VLBA, show clearly higher brightness temperatures than the typical values. Although the observed radio emission appears to be predominantly powered by star formation and supernova activity, our results provide a strong indication of the presence of an AGN in the VLBA and X-ray-detected SMG AzTEC/C61. The median radio-emitting size we have derived is 1.5-3 times larger than the typical FIR dust-emitting sizes of SMGs, but similar to that of the SMGs' molecular gas component traced through mid-$J$ line emission of CO. The physical conditions of SMGs probably render the diffusion of cosmic-ray electrons inefficient, and hence an unlikely process to lead to the observed extended radio sizes. Instead, our results point towards a scenario where SMGs are driven by galaxy interactions and mergers. Besides triggering vigorous starbursts, galaxy collisions can also pull out the magnetised fluids from the interacting disks, and give rise to a taffy-like synchrotron-emitting bridge. This provides an explanation for the spatially extended radio emission of SMGs, and can also cause a deviation from the well-known IR-radio correlation.

The observed evolution of the gas fraction and its associated depletion time in main sequence (MS) galaxies provides insights on how star formation proceeds over cosmic time. We report ALMA detections of the rest-frame $\sim$300$\mu$m continuum observed at 240 GHz for 45 massive ($\rm \langle log(M_{\star}(M_{\odot}))\rangle=10.7$), normal star forming ($\rm \langle log(sSFR(yr^{-1}))\rangle=-8.6$), i.e. MS, galaxies at $\rm z\approx3.2$ in the COSMOS field. From an empirical calibration between cold neutral, i.e. molecular and atomic, gas mass $\rm M_{gas}$ and monochromatic (rest-frame) infrared luminosity, the gas mass for this sample is derived. Combined with stellar mass $\rm M_{\star}$ and star formation rate (SFR) estimates (from {\sc MagPhys} fits) we obtain a median gas fraction of $\rm \mu_{gas}=M_{gas}/M_{\star}=1.65_{-0.19}^{+0.18}$ and a median gas depletion time $\rm t_{depl.}(Gyr)=M_{gas}/SFR=0.68_{-0.08}^{+0.07}$; correction for the location on the MS will only slightly change the values. The reported uncertainties are the $\rm 1\sigma$ error on the median. Our results are fully consistent with the expected flattening of the redshift evolution from the 2-SFM (2 star formation mode) framework that empirically prescribes the evolution assuming a universal, log-linear relation between SFR and gas mass coupled to the redshift evolution of the specific star formation rate (sSFR) of main sequence galaxies. While $\rm t_{dep.}$ shows only a mild dependence on location within the MS, a clear trend of increasing $\rm \mu_{gas}$ across the MS is observed (as known from previous studies). Further we comment on trends within the MS and (in)consistencies with other studies.

We present high-resolution (0.16$"$) 870um Atacama Large Millimeter/submillimeter Array (ALMA) imaging of 16 luminous (L_IR ~ 4 x 10^12 L_sun) submillimeter galaxies (SMGs) from the ALESS survey of the Extended Chandra Deep Field South. This dust imaging traces the dust-obscured star formation in these z~2.5 galaxies on ~1.3 kpc scales. The emission has a median effective radius of $R_e=0.24" \pm 0.02"$, corresponding to a typical physical size of $R_{e}=1.8\pm$0.2 kpc. We derive a median S\'ersic index of $n=0.9\pm0.2$, implying that the dust emission is remarkably disk-like at the current resolution and sensitivity. We use different weighting schemes with the visibilities to search for clumps on 0.12$"$ (~1.0 kpc) scales, but we find no significant evidence for clumping in the majority of cases. Indeed, we demonstrate using simulations that the observed morphologies are generally consistent with smooth exponential disks, suggesting that caution should be exercised when identifying candidate clumps in even moderate S/N interferometric data. We compare our maps to comparable-resolution HST H$_{160}$-band images, finding that the stellar morphologies appear significantly more extended and disturbed, and suggesting that major mergers may be responsible for driving the formation of the compact dust disks we observe. The stark contrast between the obscured and unobscured morphologies may also have implications for SED fitting routines that assume the dust is co-located with the optical/near-IR continuum emission. Finally, we discuss the potential of the current bursts of star formation to transform the observed galaxy sizes and light profiles, showing that the $z\sim0$ descendants of these SMGs are expected to have stellar masses, effective radii, and gas surface densities consistent with the most compact massive (M_* ~ 1-2 x 10^11 M_sun) early-type galaxies observed locally.

We present the measurement of the projected and redshift space 2-point correlation function (2pcf) of the new catalog of Chandra COSMOS-Legacy AGN at 2.9$\leq$z$\leq$5.5 ($\langle L_{bol} \rangle \sim$10$^{46}$ erg/s) using the generalized clustering estimator based on phot-z probability distribution functions (Pdfs) in addition to any available spec-z. We model the projected 2pcf estimated using $\pi_{max}$ = 200 h$^{-1}$ Mpc with the 2-halo term and we derive a bias at z$\sim$3.4 equal to b = 6.6$^{+0.60}_{-0.55}$, which corresponds to a typical mass of the hosting halos of log M$_h$ = 12.83$^{+0.12}_{-0.11}$ h$^{-1}$ M$_{\odot}$. A similar bias is derived using the redshift-space 2pcf, modelled including the typical phot-z error $\sigma_z$ = 0.052 of our sample at z$\geq$2.9. Once we integrate the projected 2pcf up to $\pi_{max}$ = 200 h$^{-1}$ Mpc, the bias of XMM and \textit{Chandra} COSMOS at z=2.8 used in Allevato et al. (2014) is consistent with our results at higher redshift. The results suggest only a slight increase of the bias factor of COSMOS AGN at z$\gtrsim$3 with the typical hosting halo mass of moderate luminosity AGN almost constant with redshift and equal to logM$_h$ = 12.92$^{+0.13}_{-0.18}$ at z=2.8 and log M$_h$ = 12.83$^{+0.12}_{-0.11}$ at z$\sim$3.4, respectively. The observed redshift evolution of the bias of COSMOS AGN implies that moderate luminosity AGN still inhabit group-sized halos at z$\gtrsim$3, but slightly less massive than observed in different independent studies using X-ray AGN at z$\leq2$.

We report interferometric measurements of [NII] 205 um fine-structure line emission from a representative sample of three galaxies at z=5-6 using the Atacama Large (sub)Millimeter Array (ALMA). These galaxies were previously detected in [CII] and far-infrared continuum emission and span almost two orders of magnitude in star formation rate (SFR). Our results show at least two different regimes of ionized inter-stellar medium properties for galaxies in the first billion years of cosmic time, separated by their L_[CII]/L_[NII] ratio. We find extremely low [NII] emission compared to [CII] (L_ [CII]/L_[NII]=68 [+200/-28]) from a "typical" L*_UV star-forming galaxy, likely directly or indirectly (by its effect on the radiation field) related to low dust abundance and low metallicity. The infrared-luminous modestly star-forming Lyman Break Galaxy (LBG) in our sample is characterized by an ionized-gas fraction (L_[CII]/L_[NII]<=20) typical of local star-forming galaxies and shows evidence for spatial variations in its ionized-gas fraction across an extended gas reservoir. The extreme SFR, warm and compact dusty starburst AzTEC-3 shows an ionized fraction higher than expected given its star-formation rate surface density (L_[CII]/L_[NII]=22+/-8) suggesting that [NII] dominantly traces a diffuse ionized medium rather than star-forming HII regions in this type of galaxy. This highest redshift sample of [NII] detections provides some of the first constraints on ionized and neutral gas modeling attempts and on the structure of the inter-stellar medium at z=5-6 in "normal" galaxies and starbursts.

We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

We present an analysis of a deep (1$\sigma$=13 $\mu$Jy) cosmological 1.2-mm continuum map based on ASPECS, the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field. In the 1 arcmin$^2$ covered by ASPECS we detect nine sources at $>3.5\sigma$ significance at 1.2-mm. Our ALMA--selected sample has a median redshift of $z=1.6\pm0.4$, with only one galaxy detected at z$>$2 within the survey area. This value is significantly lower than that found in millimeter samples selected at a higher flux density cut-off and similar frequencies. Most galaxies have specific star formation rates similar to that of main sequence galaxies at the same epoch, and we find median values of stellar mass and star formation rates of $4.0\times10^{10}\ M_\odot$ and $\sim40~M_\odot$ yr$^{-1}$, respectively. Using the dust emission as a tracer for the ISM mass, we derive depletion times that are typically longer than 300 Myr, and we find molecular gas fractions ranging from $\sim$0.1 to 1.0. As noted by previous studies, these values are lower than using CO--based ISM estimates by a factor $\sim$2. The 1\,mm number counts (corrected for fidelity and completeness) are in agreement with previous studies that were typically restricted to brighter sources. With our individual detections only, we recover $55\pm4\%$ of the extragalactic background light (EBL) at 1.2 mm measured by the Planck satellite, and we recover $80\pm7\%$ of this EBL if we include the bright end of the number counts and additional detections from stacking. The stacked contribution is dominated by galaxies at $z\sim1-2$, with stellar masses of (1-3)$\times$10$^{10}$ M$_\odot$. For the first time, we are able to characterize the population of galaxies that dominate the EBL at 1.2 mm.

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

We measure the sizes of redshift ~2 star-forming galaxies by stacking data from the Atacama Large Millimeter/submillimeter Array (ALMA). We use a uv-stacking algorithm in combination with model fitting in the uv-domain and show that this allows for robust measures of the sizes of marginally resolved sources. The analysis is primarily based on the 344 GHz ALMA continuum observations centred on 88 sub-millimeter galaxies in the LABOCA ECDFS Submillimeter Survey (ALESS). We study several samples of galaxies at z~2 with $M_* \sim{} 5\times{}10^{10} $ M$_\odot$ , selected using near-infrared photometry (distant red galaxies, extremely red objects, sBzK-galaxies, and galaxies selected on photometric redshift). We find that the typical sizes of these galaxies are ~0.6 arcsec which corresponds to ~5 kpc at z~2, this agrees well with the median sizes measured in the near-infrared z-band (~0.6 arcsec). We find errors on our size estimates of ~0.1-0.2 arcsec, which agree well with the expected errors for model fitting at the given signal-to-noise ratio. With the uv-coverage of our observations (18-160 m), the size and flux density measurements are sensitive to scales out to 2 . We compare this to a simulated ALMA Cycle 3 dataset with intermediate length baseline coverage, and we find that, using only these baselines, the measured stacked flux density would be an order of magnitude fainter. This highlights the importance of short baselines to recover the full flux density of high-redshift galaxies.

We present a search for [CII] line and dust continuum emission from optical dropout galaxies at $z>6$ using ASPECS, our ALMA Spectroscopic Survey in the Hubble Ultra-Deep Field (UDF). Our observations, which cover the frequency range $212-272$ GHz, encompass approximately the range $6<z<8$ for [CII] line emission and reach a limiting luminosity of L$_{\rm [CII]}\sim$(1.6-2.5)$\times$10$^{8}$ L$_{\odot}$. We identify fourteen [CII] line emitting candidates in this redshift range with significances $>$4.5 $\sigma$, two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60\% of our candidates are expected to be spurious. For one of our blindly selected [CII] line candidates, we tentatively detect the CO(6-5) line in our parallel 3-mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all [CII] candidates results in a tentative detection with $S_{1.2mm}=14\pm5\mu$Jy. This implies a dust-obscured star formation rate (SFR) of $(3\pm1)$ M$_\odot$ yr$^{-1}$. We find that the two highest--SFR objects have candidate [CII] lines with luminosities that are consistent with the low-redshift $L_{\rm [CII]}$ vs. SFR relation. The other candidates have significantly higher [CII] luminosities than expected from their UV--based SFR. At the current sensitivity it is unclear whether the majority of these sources are intrinsically bright [CII] emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for [CII] emitters at $6<z<8$ that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.

We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%.

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $\mu$Jy beam$^{-1}$ and 12.7 $\mu$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

We present observations of the 85.69 GHz continuum emission and H42alpha line emission from the central 30 arcsec within NGC 4945. Both sources of emission originate from nearly identical structures that can be modelled as exponential discs with a scale length of ~2.1 arcsec (or ~40 pc). An analysis of the spectral energy distribution based on combining these data with archival data imply that 84% +/- 10% of the 85.69 GHz continuum emission originates from free-free emission. The electron temperature is 5400 +/- 600 K, which is comparable to what has been measured near the centre of the Milky Way Galaxy. The star formation rate (SFR) based on the H42alpha and 85.69 GHz free-free emission (and using a distance of 3.8 Mpc) is 4.35 +/- 0.25 M/yr. This is consistent with the SFR from the total infrared flux and with previous measurements based on recombination line emission, and it is within a factor of ~2 of SFRs derived from radio data. The Spitzer Space Telescope 24 micron data and Wide-field Infrared Survey Explorer 22 micron data yield SFRs ~10x lower than the ALMA measurements, most likely because the mid-infrared data are strongly affected by dust attenuation equivalent to A_V=150. These results indicate that SFRs based on mid-infrared emission may be highly inaccurate for dusty, compact circumnuclear starbursts.

We exploit ALMA 870um observations to measure the star-formation rates (SFRs) of eight X-ray detected Active Galactic Nuclei (AGNs) in a z~3.1 protocluster, four of which reside in extended Ly-alpha haloes (often termed Ly-alpha blobs: LABs). Three of the AGNs are detected by ALMA and have implied SFRs of ~220-410~M_sun/yr; the non detection of the other five AGNs places SFR upper limits of <210 M_sun/yr. The mean SFR of the protocluster AGNs (~110-210 M_sun/yr) is consistent (within a factor of ~0.7-2.3) with that found for co-eval AGNs in the field, implying that galaxy growth is not significantly accelerated in these systems. However, when also considering ALMA data from the literature, we find evidence for elevated mean SFRs (up-to a factor of ~5.9 over the field) for AGNs at the protocluster core, indicating that galaxy growth is significantly accelerated in the central regions of the protocluster. We also show that all of the four protocluster LABs are associated with an ALMA counterpart within the extent of their Ly-alpha emission. The SFRs of the ALMA sources within the LABs (~150-410 M_sun/yr) are consistent with those expected for co-eval massive star-forming galaxies in the field. Furthermore, the two giant LABs (with physical extents of >100 kpc) do not host more luminous star formation than the smaller LABs, despite being an order of magnitude brighter in Ly-alpha emission. We use these results to discuss star formation as the power source of LABs.

We study the environment of 23 submillimetre galaxies (SMGs) drawn from the JCMT/AzTEC 1.1mm S/N-limited sample in the COSMOS field, as well as 4 COSMOS SMGs at z_spec>4.5, and 1 at z_spec=2.49, yielding a sample of 28 SMGs. We search for overdensities using the COSMOS photometric redshifts based on over 30 UV-NIR photometric bands, reaching an accuracy of sigma(Delta z/(1+z))=0.0067 (0.0155) at z<3.5 (>3.5). To identify overdensities we apply the Voronoi tessellation analysis, and estimate the overdensity estimator delta_g as a function of distance from the SMG and/or overdensity center. We test and validate our approach via simulations, X-ray detected groups, and spectroscopic verifications using VUDS and zCOSMOS catalogues showing that even with photometric redshifts in COSMOS we can efficiently retrieve overdensities out to z~5. Our results yield that 11/23 (48%) JCMT/AzTEC 1.1mm SMGs occupy overdense environments. Considering the entire JCMT/AzTEC 1.1mm S/N>4 sample, and accounting for the expected fraction of spurious detections, yields that 35-61% of the SMGs in the S/N-limited sample occupy overdense environments. We perform an X-ray stacking analysis in the 0.5-2keV band using a 32" aperture and our SMG positions, and find statistically significant detections. For our z<2 [z>2] subsample we find an average flux of (4.0+/-0.8)x10^{-16} [(1.3+/-0.5)x10^{-16}] erg/s/cm^2 and a corresponding total mass of M_200 = 2.8x10^{13} [2x10^{13}] MSol. Our results suggest a higher occurrence of SMGs occupying overdense environments at z>3, than at z<3. This may be understood if highly star forming galaxies can only be formed in the highest peaks of the density field tracing the most massive dark matter haloes at early cosmic epochs, while at later times cosmic structure may have matured sufficiently that more modest overdensities correspond to sufficiently massive haloes to form SMGs.

The COSMOS-Legacy survey is a 4.6 Ms Chandra program that has imaged 2.2 deg$^2$ of the COSMOS field with an effective exposure of $\simeq$160 ks over the central 1.5 deg$^2$ and of $\simeq$80 ks in the remaining area. The survey is the combination of 56 new observations, obtained as an X-ray Visionary Project, with the previous C-COSMOS survey. We describe the reduction and analysis of the new observations and the properties of 2273 point sources detected above a spurious probability of 2$\times 10^{-5}$. We also present the updated properties of the C-COSMOS sources detected in the new data. The whole survey includes 4016 point sources (3814, 2920 and 2440 in the full, soft and hard band). The limiting depths are 2.2 $\times$ 10$^{-16}$, 1.5 $\times$ 10$^{-15}$ and 8.9$\times$ 10$^{-16}$ ${\rm erg~cm}^{-2}~{\rm s}^{-1}$ in the 0.5-2, 2-10 and 0.5-10 keV bands, respectively. The observed fraction of obscured AGN with column density $> 10^{22}$ cm$^{-2}$ from the hardness ratio (HR) is $\sim$50$^{+17}_{-16}$%. Given the large sample, we compute source number counts in the hard and soft bands, significantly reducing the uncertainties of 5-10%. For the first time, we compute number counts for obscured (HR$>$-0.2) and unobscured (HR$<$-0.2) sources and find significant differences between the two populations in the soft band. Due to the un-precedent large exposure, COSMOS-Legacy area is 3 times larger than surveys at similar depth and its depth is 3 times fainter than surveys covering similar area. The area-flux region occupied by COSMOS-Legacy is likely to remain unsurpassed for years to come.

We present the catalog of optical and infrared counterparts of the Chandra COSMOS-Legacy Survey, a 4.6 Ms Chandra program on the 2.2 square degrees of the COSMOS field, combination of 56 new overlapping observations obtained in Cycle 14 with the previous C-COSMOS survey. In this Paper we report the i, K, and 3.6 micron identifications of the 2273 X-ray point sources detected in the new Cycle 14 observations. We use the likelihood ratio technique to derive the association of optical/infrared (IR) counterparts for 97% of the X-ray sources. We also update the information for the 1743 sources detected in C-COSMOS, using new K and 3.6 micron information not available when the C-COSMOS analysis was performed. The final catalog contains 4016 X-ray sources, 97% of which have an optical/IR counterpart and a photometric redshift, while 54% of the sources have a spectroscopic redshift. The full catalog, including spectroscopic and photometric redshifts and optical and X-ray properties described here in detail, is available online. We study several X-ray to optical (X/O) properties: with our large statistics we put better constraints on the X/O flux ratio locus, finding a shift towards faint optical magnitudes in both soft and hard X-ray band. We confirm the existence of a correlation between X/O and the the 2-10 keV luminosity for Type 2 sources. We extend to low luminosities the analysis of the correlation between the fraction of obscured AGN and the hard band luminosity, finding a different behavior between the optically and X-ray classified obscured fraction.

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

We consider the high radio frequency (15 GHz - 353 GHz) properties and variability of 35 Brightest Cluster Galaxies (BCGs). These are the most core-dominated sources drawn from a parent sample of more than 700 X-ray selected clusters, thus allowing us to relate our results to the general population. We find that >6.0% of our parent sample (>15.1% if only cool-core clusters are considered) contain a radio-source at 150 GHz of at least 3mJy (~1x10^23 W/Hz at our median redshift of z~0.13). Furthermore, >3.4% of the BCGs in our parent sample contain a peaked component (Gigahertz Peaked Spectrum, GPS) in their spectra that peaks above 2 GHz, increasing to >8.5% if only cool-core clusters are considered. We see little evidence for strong variability at 15 GHz on short (week-month) time-scales although we see variations greater than 20% at 150 GHz over 6-month times-frames for 4 of the 23 sources with multi-epoch observations. Much more prevalent is long-term (year-decade time-scale) variability, with average annual amplitude variations greater than 1% at 15 GHz being commonplace. There is a weak trend towards higher variability as the peak of the GPS-like component occurs at higher frequency. We demonstrate the complexity that is seen in the radio spectra of BCGs and discuss the potentially significant implications of these high-peaking components for Sunyaev-Zel'dovich cluster searches.

We present high-resolution 870-um ALMA continuum maps of 30 bright sub-millimeter sources in the UKIDSS UDS field. These sources are selected from deep, 1-square degrees 850-um maps from the SCUBA--2 Cosmology Legacy Survey, and are representative of the brightest sources in the field (median SCUBA2 flux S_850=8.7+/-0.4 mJy). We detect 52 sub-millimeter galaxies (SMGs) at >4-sigma significance in our 30 ALMA maps. In 61+/-17% of the ALMA maps the single-dish source comprises a blend of >=2 SMGs, where the secondary SMGs are Ultra--Luminous Infrared Galaxies (ULIRGs) with L_IR>10^12 Lo. The brightest SMG contributes on average 80+/-4% of the single-dish flux density, and in the ALMA maps containing >=2 SMGs the secondary SMG contributes 25+/-3% of the integrated ALMA flux. We construct source counts and show that multiplicity boosts the apparent single-dish cumulative counts by 20% at S_870>7.5mJy, and by 60% at S_870>12mJy. We combine our sample with previous ALMA studies of fainter SMGs and show that the counts are well-described by a double power-law with a break at 8.5+/-0.6mJy. The break corresponds to a luminosity of ~6x10^12Lsol or a star-formation rate of ~1000Mo/yr. For the typical sizes of these SMGs, which are resolved in our ALMA data with r=1.2+/-0.1kpc, this yields a limiting SFR density of ~100Msol/yr/kpc2. Finally, the number density of S_870>2mJy SMGs is 80+/-30 times higher than that derived from blank-field counts. An over-abundance of faint SMGs is inconsistent with line-of-sight projections dominating multiplicity in the brightest SMGs, and indicates that a significant proportion of these high-redshift ULIRGs must be physically associated.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

[abridged] The ALESS survey has followed-up a sample of 122 sub-millimeter sources in the Extended Chandra Deep Field South at 870um with ALMA, allowing to pinpoint the positions of sub-millimeter galaxies (SMGs) to 0.3'' and to find their precise counterparts at different wavelengths. This enabled the first compilation of the multi-wavelength spectral energy distributions (SEDs) of a statistically reliable survey of SMGs. In this paper, we present a new calibration of the MAGPHYS modelling code that is optimized to fit these UV-to-radio SEDs of z>1 star-forming galaxies using an energy balance technique to connect the emission from stellar populations, dust attenuation and dust emission in a physically consistent way. We derive statistically and physically robust estimates of the photometric redshifts and physical parameters for the ALESS SMGs. We find that they have a median stellar mass $M_\ast=(8.9\pm0.1)\times10^{10} M_\odot$, SFR$=280\pm70 M_\odot$/yr, overall V-band dust attenuation $A_V=1.9\pm0.2$ mag, dust mass $M_\rm{dust}=(5.6\pm1.0)\times10^8 M_\odot$, and average dust temperature Tdust~40 K. The average intrinsic SED of the ALESS SMGs resembles that of local ULIRGs in the IR range, but the stellar emission of our average SMG is brighter and bluer, indicating lower dust attenuation, possibly because they are more extended. We explore how the average SEDs vary with different parameters, and we provide a new set of SMG templates. To put the ALESS SMGs into context, we compare their stellar masses and SFRs with those of less actively star-forming galaxies at the same redshifts. At z~2, about half of the SMGs lie above the star-forming main sequence, while half are at the high-mass end of the sequence. At higher redshifts (z~3.5), the SMGs tend to have higher SFR and Mstar, but the fraction of SMGs that lie significantly above the main sequence decreases to less than a third.

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

We present the first spectroscopic confirmation of an ultra-massive galaxy at redshift z>3 using data from Keck-NIRSPEC, VLT-Xshooter, and GTC-Osiris. We detect strong [OIII] and Ly$\alpha$ emission, and weak [OII], CIV, and HeII, placing C1-23152 at a spectroscopic redshift of $z_{spec}$=3.351. The modeling of the emission-line corrected spectral energy distribution results in a best-fit stellar mass of $M_{*}=3.1^{+0.6}_{-0.7}\times10^{11} M_{\odot}$, a star-formation rate of <7 $M_{\odot} yr^{-1}$, and negligible dust extinction. The stars appear to have formed in a short intense burst ~300-500 Myr prior to the observation epoch, setting the formation redshift of this galaxy at z~4.1. From the analysis of the line ratios and widths, and the observed flux at 24$\mu$m, we confirm the presence of a luminous hidden active galactic nucleus (AGN), with bolometric luminosity of ~$10^{46}erg s^{-1}$. Potential contamination to the observed SED from the AGN continuum is constrained, placing a lower limit on the stellar mass of $2\times10^{11} M_{\odot}$. HST/WFC3 $H_{160}$ and ACS $I_{814}$ images are modeled, resulting in an effective radius of $r_{e}$~1 kpc in the $H_{160}$ band and a Sersic index n~4.4. This object may be a prototype of the progenitors of local most massive elliptical galaxies in the first 2 Gyr of cosmic history, having formed most of its stars at z>4 in a highly dissipative, intense, and short burst of star formation. C1-23152 is completing its transition to a post-starburst phase while hosting a powerful AGN, potentially responsible for the quenching of the star formation activity.

The relation between the stellar mass and the star formation rate characterizes how the instantaneous star formation is determined by the galaxy past star formation history and by the growth of the dark matter structures. We deconstruct the M-SFR plane by measuring the specific SFR functions in several stellar mass bins from z=0.2 out to z=1.4. Our analysis is primary based on a MIPS 24$\mu m$ selected catalogue combining the COSMOS and GOODS surveys. We estimate the SFR by combining mid- and far-infrared data for 20500 galaxies. The sSFR functions are derived in four stellar mass bins within the range 9.5<log(M/Msun)<11.5. First, we demonstrate the importance of taking into account selection effects when studying the M-SFR relation. Secondly, we find a mass-dependent evolution of the median sSFR with redshift varying as $sSFR \propto (1+z)^{b}$, with $b$ increasing from $b=2.88$ to $b=3.78$ between $M=10^{9.75}Msun$ and $M=10^{11.1}Msun$, respectively. At low masses, this evolution is consistent with the cosmological accretion rate and predictions from semi-analytical models (SAM). This agreement breaks down for more massive galaxies showing the need for a more comprehensive description of the star-formation history in massive galaxies. Third, we obtain that the shape of the sSFR function is invariant with time at z<1.4 but depends on the mass. We observe a broadening of the sSFR function ranging from 0.28 dex at $M=10^{9.75}Msun$ to 0.46 dex at $M=10^{11.1}Msun$. Such increase in the scatter of the M-SFR relation suggests an increasing diversity of SFHs as the stellar mass increases. Finally, we find a gradual decline of the sSFR with mass as $log(sSFR) \propto -0.17M$. We discuss the numerous physical processes, as gas exhaustion in hot gas halos or secular evolution, which can gradually reduce the sSFR and increase the SFH diversity.

Evolution in the measured rest frame ultraviolet spectral slope and ultraviolet to optical flux ratios indicate a rapid evolution in the dust obscuration of galaxies during the first 3 billion years of cosmic time (z>4). This evolution implies a change in the average interstellar medium properties, but the measurements are systematically uncertain due to untested assumptions, and the inability to measure heavily obscured regions of the galaxies. Previous attempts to directly measure the interstellar medium in normal galaxies at these redshifts have failed for a number of reasons with one notable exception. Here we report measurements of the [CII] gas and dust emission in 9 typical (~1-4L*) star-forming galaxies ~1 billon years after the big bang (z~5-6). We find these galaxies have >12x less thermal emission compared with similar systems ~2 billion years later, and enhanced [CII] emission relative to the far-infrared continuum, confirming a strong evolution in the interstellar medium properties in the early universe. The gas is distributed over scales of 1-8 kpc, and shows diverse dynamics within the sample. These results are consistent with early galaxies having significantly less dust than typical galaxies seen at z<3 and being comparable to local low-metallicity systems.

We used the PdBI to map a sample of 15 SMGs in the COSMOS field at the wavelength of 1.3 mm. The target SMGs were originally discovered in the JCMT/AzTEC 1.1 mm continuum survey at S/N=4-4.5. This paper presents, for the first time, interferometric millimetre-wavelength observations of these sources. The angular resolution of our observations, 1.8", allowed us to accurately determine the positions of the target SMGs. Using a detection threshold of S/N>4.5 regardless of multiwavelength counterpart association, and 4<S/N<=4.5 if a multiwavelength counterpart within 1.5" is also present, the total number of detections in our survey is 22. Three of our detected SMGs (AzTEC21, 27, and 28; which corresponds to 20%) are marginally resolved at our angular resolution, and these sources are found to have elongated or clumpy morphologies and/or multiple components. Using optical to NIR photometric redshifts, available spectroscopic redshifts, and redshifts estimated from the radio-to-submm spectral index we infer a median redshift of $\tilde{z}=3.20\pm0.25$ for our sample. To study the overall multiplicity and redshift distribution of flux-limited samples of SMGs we combined these sources with the 15 brightest AzTEC SMGs detected at 1.1 mm, AzTEC1-15, and studied previously. This constitutes a complete, flux- and S/N-limited 1.1-mm selected sample. We find that the median redshift for the 15 brightest AzTEC SMGs ($\tilde{z}=3.05\pm0.44$) is consistent with that for AzTEC16-30. This conforms to recent observational findings that SMGs do not exhibit any significant trend between the redshift and (sub)mm flux density. For the combined AzTEC1-30 sample we derive a median redshift of $\tilde{z}=3.17\pm0.27$, consistent with previous results based on mm-selected samples. We further infer that within the combined AzTEC1-30 sample $\sim25\pm9\%$ of sources separate into multiple components.

In this Chapter we present the motivation for undertaking both a wide and deep survey with the SKA in the context of studying AGN activity across cosmic time. With an rms down to 1 $\mu$Jy/beam at 1 GHz over 1,000 - 5,000 deg$^2$ in 1 year (wide tier band 1/2) and an rms down to 200 nJy/beam over 10 - 30 deg$^2$ in 2000 hours (deep tier band 1/2), these surveys will directly detect faint radio-loud and radio-quiet AGN (down to a 1 GHz radio luminosity of about $2\times10^{23}$ W/Hz at $z=6$). For the first time, this will enable us to conduct detailed studies of the cosmic evolution of radio AGN activity to the cosmic dawn ($z\gtrsim6$), covering all environmental densities.

We study the physical properties of a sample of 6 SMGs in the COSMOS field, spectroscopically confirmed to lie at z>4. We use new GMRT 325 MHz and 3 GHz JVLA data to probe the rest-frame 1.4 GHz emission at z=4, and to estimate the sizes of the star-forming (SF) regions of these sources, resp. Combining our size estimates with those available in the literature for AzTEC1 and AzTEC3 we infer a median radio-emitting size for our z>4 SMGs of (0.63"+/-0.12")x(0.35"+/-0.05") or 4.1x2.3 kpc^2 (major times minor axis; assuming z=4.5) or lower if we take the two marginally resolved SMGs as unresolved. This is consistent with the sizes of SF regions in lower-redshift SMGs, and local normal galaxies, yet higher than the sizes of SF regions of local ULIRGs. Our SMG sample consists of a fair mix of compact and more clumpy systems with multiple, perhaps merging, components. With an average formation time of ~280 Myr, derived through modeling of the UV-IR SEDs, the studied SMGs are young systems. The average stellar mass, dust temperature, and IR luminosity we derive are M*~1.4x10^11 M_sun, T_dust~43 K, and L_IR~1.3x10^13L_sun, resp. The average L_IR is up to an order of magnitude higher than for SMGs at lower redshifts. Our SMGs follow the correlation between dust temperature and IR luminosity as derived for Herschel-selected 0.1<z<2 galaxies. We study the IR-radio correlation for our sources and find a deviation from that derived for z<3 ULIRGs (<q_IR>=1.95+/-0.26 for our sample, compared to q~2.6 for IR luminous galaxies at z<2). In summary, we find that the physical properties derived for our z>4 SMGs put them at the high end of the L_IR-T_dust distribution of SMGs, and that our SMGs form a morphologically heterogeneous sample. Thus, further in-depth analyses of large, statistical samples of high-redshift SMGs are needed to fully understand their role in galaxy formation and evolution.

In this chapter, we highlight a number of science investigations that are enabled by the inclusion of Band~5 ($4.6-13.8$ GHz) for SKA1-MID science operations, while focusing on the astrophysics of star formation over cosmic time. For studying the detailed astrophysics of star formation at high-redshift, surveys at frequencies $\gtrsim$10 GHz have the distinct advantage over traditional $\sim$1.4 GHz surveys as they are able to yield higher angular resolution imaging while probing higher rest frame frequencies of galaxies with increasing redshift, where emission of star-forming galaxies becomes dominated by thermal (free-free) radiation. In doing so, surveys carried out at $\gtrsim$10 GHz provide a robust, dust-unbiased measurement of the massive star formation rate by being highly sensitive to the number of ionizing photons that are produced. To access this powerful star formation rate diagnostic requires that Band~5 be available for SKA1-MID. We additionally present a detailed science case for frequency coverage extending up to 30 GHz during full SKA2 operations, as this allows for highly diverse science while additionally providing contiguous frequency coverage between the SKA and ALMA, which will likely be the two most powerful interferometers for the coming decades. To enable this synergy, it is crucial that the dish design of the SKA be flexible enough to include the possibility of being fit with receivers operating up to 30 GHz.

Astronomy at or below the 'survey threshold' has expanded significantly since the publication of the original 'Science with the Square Kilometer Array' in 1999 and its update in 2004. The techniques in this regime may be broadly (but far from exclusively) defined as 'confusion' or 'P(D)' analyses (analyses of one-point statistics), and 'stacking', accounting for the flux-density distribution of noise-limited images co-added at the positions of objects detected/isolated in a different waveband. Here we discuss the relevant issues, present some examples of recent analyses, and consider some of the consequences for the design and use of surveys with the SKA and its pathfinders.

We analyse HST WFC3/$H_{160}$-band observations of a sample of 48 ALMA-detected submillimeter galaxies (SMGs) in the Extended Chandra Deep Field South field, to study their stellar morphologies and sizes. We detect 79$\pm$17% of the SMGs in the $H_{160}$-band imaging with a median sensitivity of 27.8 mag, and most (80%) of the non-detections are SMGs with 870$\mu$m fluxes of $S_{870} < $3 mJy. With a surface brightness limit of $\mu_H \sim $26 mag arcsec$^{-2}$, we find that 82$\pm$9% of the $H_{160}$-band detected SMGs at $z =$ 1-3 appear to have disturbed morphologies, meaning they are visually classified as either irregulars or interacting systems, or both. By determining a S\'ersic fit to the $H_{160}$ surface-brightness profiles we derive a median S\'ersic index of $n = $1.2$\pm$0.3 and a median half-light radius of $r_e = $4.4$^{+1.1}_{-0.5}$ kpc for our SMGs at $z = $1-3. We also find significant displacements between the positions of the $H_{160}$-component and 870$\mu$m emission in these systems, suggesting that the dusty star-burst regions and less-obscured stellar distribution are not co-located. We find significant differences in the sizes and the S\'ersic index between our $z = $2-3 SMGs and $z \sim $2 quiescent galaxies, suggesting a major transformation of the stellar light profile is needed in the quenching processes if SMGs are progenitors of the red-and-dead $z\sim$2 galaxies. Given the short-lived nature of SMGs, we postulate that the majority of the $z = $2-3 SMGs with $S_{870} \gtrsim $2 mJy are early/mid-stage major mergers.

We aim to measure the average dust and molecular gas content of massive star-forming galaxies ($\rm > 3 \times 10^{10}\,M_\odot$) up to z=4 in the COSMOS field to determine if the intense star formation observed at high redshift is induced by major mergers or caused by large gas reservoirs. Firstly, we measured the evolution of the average spectral energy distributions as a function of redshift using a stacking analysis of Spitzer, Herschel, LABOCA, and AzTEC data for two samples of galaxies: normal star-forming objects and strong starbursts, as defined by their distance to the main sequence. We found that the mean intensity of the radiation field $< U >$ heating the dust (strongly correlated with dust temperature) increases with increasing redshift up to z$\sim$4 in main-sequence galaxies. We can reproduce this evolution with simple models that account for the decrease of the gas metallicity with redshift. No evolution of $< U >$ with redshift is found in strong starbursts. We then deduced the evolution of the molecular gas fraction (defined here as $\rm M_{\rm mol}/(M_{\rm mol}+M_\star)$) with redshift and found a similar, steeply increasing trend for both samples. At z$\sim$4, this fraction reaches $\sim$60%. The average position of the main-sequence galaxies is on the locus of the local, normal star-forming disks in the integrated Schmidt-Kennicutt diagram (star formation rate versus mass of molecular gas), suggesting that the bulk of the star formation up to z=4 is dominated by secular processes.

We present high-resolution (0.3'') ALMA 870um imaging of 52 sub-millimeter galaxies (SMGs) in the Ultra Deep Survey (UDS) field and investigate the size and morphology of the sub-millimeter (sub-mm) emission on 2-10kpc scales. We derive a median intrinsic angular size of FWHM=0.30$\pm$0.04'' for the 23 SMGs in the sample detected at a signal-to-noise ratio (SNR) >10. Using the photometric redshifts of the SMGs we show that this corresponds to a median physical half-light diameter of 2.4$\pm$0.2kpc. A stacking analysis of the SMGs detected at an SNR <10 shows they have sizes consistent with the 870um-bright SMGs in the sample. We compare our results to the sizes of SMGs derived from other multi-wavelength studies, and show that the rest-frame ~250um sizes of SMGs are consistent with studies of resolved 12CO (J=3-2 to 7-6) emission lines, but that sizes derived from 1.4GHz imaging appear to be approximately two times larger on average, which we attribute to cosmic ray diffusion. The rest-frame optical sizes of SMGs are around four times larger than the sub-millimeter sizes, indicating that the star formation in these galaxies is compact relative to the pre-existing stellar distribution. The size of the starburst region in SMGs is consistent with the majority of the star formation occurring in a central region, a few kpc in extent, with a median star formation rate surface density of 90$\pm$30Msol/yr/kpc$^2$, which may suggest that we are witnessing an intense period of bulge growth in these galaxies.

We present stringent constraints on the average mid-, far-infrared and radio emissions of $\sim$14200 quiescent galaxies (QGs), identified out to $z=3$ in the COSMOS field via their rest-frame NUV$-$r and r$-$J colors, and with stellar masses $M_{\star}=10^{9.8-12.2} \,M_{\odot} $. Stacking in deep Spitzer (MIPS $24\,\mu$m), Herschel (PACS and SPIRE), and VLA (1.4 GHz) maps reveals extremely low dust-obscured star formation rates for QGs (SFR $<0.1-3\,M_{\odot}$yr$^{-1}$ at $z \leqslant 2$ and $<6-18\,M_{\odot}$yr$^{-1}$ at $z > 2$), consistent with the low unobscured SFRs ($<0.01-1.2\,M_{\odot}$yr$^{-1}$) inferred from modeling their ultraviolet-to-near-infrared photometry. The average SFRs of QGs are $>10\times$ below those of star-forming galaxies (SFGs) within the $M_{\star}$- and $z$-ranges considered. The stacked 1.4 GHz signals (S/N $> 5$) are, if attributed solely to star formation, in excess of the total (obscured plus unobscured) SFR limits, suggestive of a widespread presence of low-luminosity active galactic nuclei (AGN) among QGs. Our results reaffirm the existence of a significant population QGs out to $z = 3$, thus corroborating the need for powerful quenching mechanism(s) to terminate star formation in galaxies at earlier epochs.

[Abridged] We study the evolution of the radio spectral index and far-infrared/radio correlation (FRC) across the star-formation rate-stellar masse (i.e. SFR-M*) plane up to z 2. We start from a M*-selected sample of galaxies with reliable SFR and redshift estimates. We then grid the SFR-M* plane in several redshift ranges and measure the infrared luminosity, radio luminosity, radio spectral index, and ultimately the FRC index (i.e. qFIR) of each SFR-M*-z bin. The infrared luminosities of our SFR-M*-z bins are estimated using their stacked far-infrared flux densities inferred from observations obtained with Herschel. Their radio luminosities and radio spectral indices (i.e. alpha, where Snu nu^-alpha) are estimated using their stacked 1.4GHz and 610MHz flux densities from the VLA and GMRT, respectively. Our far-infrared and radio observations include the most widely studied blank extragalactic fields -GOODS-N/S, ECDFS, and COSMOS- covering a sky area of 2deg^2. Using this methodology, we constrain the radio spectral index and FRC index of star-forming galaxies with M*>10^10Msun and 0<z<2.3. We find that alpha^1.4GHz_610MHz does not evolve significantly with redshift or with the distance of a galaxy with respect to the main sequence (MS) of the SFR-M* plane (i.e. Delta_log(SSFR)_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). Instead, star-forming galaxies have a radio spectral index consistent with a canonical value of 0.8, which suggests that their radio spectra are dominated by non-thermal optically thin synchrotron emission. We find that qFIR displays a moderate but statistically significant redshift evolution as qFIR(z)=(2.35+/-0.08)*(1+z)^(-0.12+/-0.04), consistent with some previous literature. Finally, we find no significant correlation between qFIR and Delta_log(SSFR)_MS, though a weak positive trend, as observed in one of our redshift bins, cannot be firmly ruled out using our dataset.

We report interferometric imaging of [CII] and OH emission toward the center of the galaxy protocluster associated with the z=5.3 submillimeter galaxy (SMG) AzTEC-3, using the Atacama Large (sub)Millimeter Array (ALMA). We detect strong [CII], OH, and rest-frame 157.7 um continuum emission toward the SMG. The [CII] emission is distributed over a scale of 3.9 kpc, implying a dynamical mass of 9.7 x 10^10 Msun, and a star formation rate (SFR) surface density of Sigma_SFR = 530 Msun/yr/kpc2. This suggests that AzTEC-3 forms stars at Sigma_SFR approaching the Eddington limit for radiation pressure supported disks. We find that the OH emission is slightly blueshifted relative to the [CII] line, which may indicate a molecular outflow associated with the peak phase of the starburst. We also detect and dynamically resolve [CII] emission over a scale of 7.5 kpc toward a triplet of Lyman-break galaxies with moderate UV-based SFRs in the protocluster at ~95kpc projected distance from the SMG. These galaxies are not detected in the continuum, suggesting far-infrared SFRs of <18-54 Msun/yr, consistent with a UV-based estimate of 22 Msun/yr. The spectral energy distribution of these galaxies is inconsistent with nearby spiral and starburst galaxies, but resembles those of dwarf galaxies. This is consistent with expectations for young starbursts without significant older stellar populations. This suggests that these galaxies are significantly metal-enriched, but not heavily dust-obscured, "normal" star-forming galaxies at z>5, showing that ALMA can detect the interstellar medium in "typical" galaxies in the very early universe.

We present statistically significant detections at 850um of the Lyman Break Galaxy (LBG) population at z=3, 4, and 5 using data from the Submillimetre Common User Bolometer Array 2 (SCUBA-2) Cosmology Legacy Survey (S2CLS) in the United Kingdom Infrared Deep Sky Survey Ultra Deep Survey (UKIDSS-UDS) field. We employ a stacking technique to probe beneath the survey limit to measure the average 850um flux density of LBGs at z=3, 4, and 5 with typical ultraviolet luminosities of L(1700A)~10^29 erg/s/Hz. We measure 850um flux densities of (0.25 +/- 0.03, (0.41 +/- 0.06), and (0.88 +/- 0.23) mJy respectively, and find that they contribute at most 20 per cent to the cosmic far-infrared background at 850um. Fitting an appropriate range of spectral energy distributions to the z=3, 4, and 5 LBG stacked 24-850um fluxes, we derive infrared (IR) luminosities of L(8-1000um)~3.2, 5.5, and 11.0x10^11 Lsun (corresponding to star formation rates of ~50-200 Msun/yr) respectively. We find that the evolution in the IR luminosity density of LBGs is broadly consistent with model predictions for the expected contribution of luminous IR galaxy (LIRG) to ultraluminous IR galaxy (ULIRG) type systems at these epochs. We also see a strong positive correlation between stellar mass and IR luminosity. Our data are consistent with the main sequence of star formation showing little or no evolution from z=3 to 5. We have also confirmed that, for a fixed mass, the reddest LBGs (UV slope Beta -> 0) are indeed redder due to dust extinction, with SFR(IR)/SFR(UV) increasing by approximately an order of magnitude over -2<Beta<0 such that SFR(IR)/SFR(UV)~20 for the reddest LBGs. Furthermore, the most massive LBGs also tend to have higher obscured-to-unobscured ratio, hinting at a variation in the obscuration properties across the mass range.

We present 90 cm VLA imaging of the COSMOS field, comprising a circular area of 3.14 square degrees at 8.0"x6.0" angular resolution with an average rms of 0.5 mJy/beam. The extracted catalog contains 182 sources (down to 5.5sigma), 30 of which are multi-component sources. Using Monte Carlo artificial source simulations we derive the completeness of the catalog, and we show that our 90 cm source counts agree very well with those from previous studies. Using X-ray, NUV-NIR and radio COSMOS data to investigate the population mix of our 90 cm radio sample, we find that our sample is dominated by active galactic nuclei (AGN). The average 90-20 cm spectral index (S_nu~nu**alpha, where S_nu is the flux density at frequency nu, and alpha the spectral index) of our 90 cm selected sources is -0.70, with an interquartile range of -0.90 to -0.53. Only a few ultra-steep-spectrum sources are present in our sample, consistent with results in the literature for similar fields. Our data do not show clear steepening of the spectral index with redshift. Nevertheless, our sample suggests that sources with spectral indices steeper than -1 all lie at z>1, in agreement with the idea that ultra-steep-spectrum radio sources may trace intermediate-redshift galaxies (z>1).

We present deep continuum observations using the GISMO camera at a wavelength of 2 mm centered on the Hubble Deep Field (HDF) in the GOODS-N field. These are the first deep field observations ever obtained at this wavelength. The 1 sigma sensitivity in the innermost approx. 4 arcminutes of the 7 utes map is approx. 135 uJy/beam, a factor of three higher in flux/beam sensitivity than the deepest available SCUBA 850 um observations, and almost a factor of four higher in flux / beam sensitivity than the combined MAMBO/AzTEC 1.2 mm observations of this region. Our source extraction algorithm identifies 12 sources directly, and another 3 through correlation with known sources at 1.2 mm and 850 um. Five of the directly detected GISMO sources have counterparts in the MAMBO/AzTEC catalog, and four of those also have SCUBA counterparts. HDF850.1, one of the first blank-field detected submillimeter galaxies, is now detected at 2 mm. The median redshift of all sources with counterparts of known redshifts is med(z) = 2.91 +/- 0.94. Statistically, the detections are most likely real for 5 of the seven 2 mm sources without shorter wavelength counterparts, while the probability for none of them being real is negligible.

We present a study of the radio properties of 870$\mu$m-selected submillimetre galaxies (SMGs), observed at high resolution with ALMA in the Extended Chandra Deep Field South. From our initial sample of 76 ALMA SMGs, we detect 52 SMGs at $>3\sigma$ significance in VLA 1400MHz imaging, of which 35 are also detected at $>3\sigma$ in new 610MHz GMRT imaging. Within this sample of radio-detected SMGs, we measure a median radio spectral index $\alpha_{610}^{1400} = -0.79 \pm 0.06$, (with inter-quartile range $\alpha=[-1.16,-0.56]$) and investigate the far-infrared/radio correlation via the parameter $q_{\rm IR}$, the logarithmic ratio of the rest-frame 8-1000$\mu$m flux and monochromatic radio flux. Our median $q_{\rm IR} = 2.56 \pm 0.05$ (inter-quartile range $q_{\rm IR}=[2.42,2.78]$) is higher than that typically seen in single-dish 870$\mu$m-selected sources ($q_{\rm IR} \sim 2.4$), which may reflect the fact that our ALMA-based study is not biased to radio-bright counterparts, as previous samples were. Finally, we search for evidence that $q_{\rm IR}$ and $\alpha$ evolve with age in a co-dependent manner, as predicted by starburst models: the data populate the predicted region of parameter space, with the stellar mass tending to increase along tracks of $q_{\rm IR}$ versus $\alpha$ in the direction expected, providing the first observational evidence in support of these models.

We present the first photometric redshift distribution for a large unbiased sample of 870um selected submillimeter galaxies (SMGs) with robust identifications based on observations with the Atacama Large Millimeter Array (ALMA). In our analysis we consider 96 SMGs in the Extended Chandra Deep Field South, 77 of which have 4-19 band, optical-near-infrared, photometry. We model the Spectral Energy Distributions (SEDs) for these 77 SMGs, deriving a median photometric redshift of z=2.3+/-0.1. The remaining 19 SMGs have insufficient optical or near-infrared photometry to derive photometric redshifts, but a stacking analysis of IRAC and Herschel observations confirms they are not spurious. Assuming these sources have an absolute H-band magnitude distribution comparable to that of a complete sample of z~1-2 SMGs, we demonstrate that the undetected SMGs lie at higher redshifts, raising the median redshift for SMGs to z=2.5+/-0.2. More critically we show that the proportion of galaxies undergoing an SMG phase at z>3 is 35+/-5% of the total population. We derive a median stellar mass for SMGs of Mstar=(8+/-1)x10^10Mo, but caution that there are significant systematic uncertainties in our stellar mass estimate, up to x5 for individual sources. We compare our sample of SMGs to a volume-limited, morphologically classified sample of ellipticals in the local Universe. Assuming the star formation activity in SMGs has a timescale of ~100Myr we show that their descendants at z~0 would have a space density and M_H distribution which are in good agreement with those of local ellipticals. In addition the inferred mass-weighted ages of the local ellipticals broadly agree with the look-back times of the SMG events. Taken together, these results are consistent with a simple model that identifies SMGs as events that form most of the stars seen in the majority of luminous elliptical galaxies at the present day.

We exploit ALMA 870um (345GHz) observations of submillimetre sources in the Extended Chandra Deep Field South to investigate the far-infrared properties of high-redshift submillimetre galaxies (SMGs). Using the precisely located 870um ALMA positions of 99 SMGs, together with 24um and radio imaging of this field, we deblend the Herschel/SPIRE imaging of this region to extract their far-infrared fluxes and colours. The median photometric redshifts for ALMA LESS (ALESS) SMGs which are detected in at least two SPIRE bands increases with wavelength of the peak in their SEDs, with z=2.3+/-0.2, 2.5+/-0.3 and 3.5+/-0.5 for the 250, 350 and 500-um peakers respectively. We find that 34 ALESS SMGs do not have a >3-sigma counterpart at 250, 350 or 500-um. These galaxies have a median photometric redshift of z=3.3+/-0.5, which is higher than the full ALESS SMG sample; z=2.5+/-0.2. Using the photometric redshifts together with the 250-870um photometry, we estimate the far-infrared luminosities and characteristic dust temperature of each SMG. The median infrared luminosity of the S_870um>2mJy SMGs is L_IR=(3.0+/-0.3)x10^{12}Lo(SFR=300+/-30Mo/yr). At a fixed luminosity, the characteristic dust temperature of these high-redshift SMGs is 2-3K lower than comparably luminous galaxies at z=0, reflecting the more extended star formation occurring in these systems. By extrapolating the 870um number counts to S_ 870um=1mJy, we show that the contribution of S_870um>1mJy SMGs to the cosmic star formation budget is 20% of the total over the redshift range z~1-4. We derive a median dust mass for these SMGs of M_d=(3.6+/-0.3)x10^8Mo and by adopting an appropriate gas-to-dust ratio, we estimate an average molecular mass of M_H2=(4.2+/-0.4)x10^{10}Mo. Finally, we use our estimates of the H2 masses to show that SMGs with S_870um>1mJy contain ~10% of the z~2 volume-averaged H2 mass density at this epoch.

Three billion years after the big bang (at redshift z=2), half of the most massive galaxies were already old, quiescent systems with little to no residual star formation and extremely compact with stellar mass densities at least an order of magnitude larger than in low redshift ellipticals, their descendants. Little is known about how they formed, but their evolved, dense stellar populations suggest formation within intense, compact starbursts 1-2 Gyr earlier (at 3<z<6). Simulations show that gas-rich major mergers can give rise to such starbursts which produce dense remnants. Sub-millimeter selected galaxies (SMGs) are prime examples of intense, gas-rich, starbursts. With a new, representative spectroscopic sample of compact quiescent galaxies at z=2 and a statistically well-understood sample of SMGs, we show that z=3-6 SMGs are consistent with being the progenitors of z=2 quiescent galaxies, matching their formation redshifts and their distributions of sizes, stellar masses and internal velocities. Assuming an evolutionary connection, their space densities also match if the mean duty cycle of SMG starbursts is 42 (+40/-29) Myr (consistent with independent estimates), which indicates that the bulk of stars in these massive galaxies were formed in a major, early surge of star-formation. These results suggests a coherent picture of the formation history of the most massive galaxies in the universe, from their initial burst of violent star-formation through their appearance as high stellar-density galaxy cores and to their ultimate fate as giant ellipticals.

We analyse new SCUBA-2 submillimeter and archival SPIRE far-infrared imaging of a z=1.62 cluster, Cl0218.3-0510, which lies in the UKIDSS/UDS field of the SCUBA-2 Cosmology Legacy Survey. Combining these tracers of obscured star formation activity with the extensive photometric and spectroscopic information available for this field, we identify 31 far-infrared/submillimeter-detected probable cluster members with bolometric luminosities >1e12 Lo and show that by virtue of their dust content and activity, these represent some of the reddest and brightest galaxies in this structure. We exploit Cycle-1 ALMA submillimeter continuum imaging which covers one of these sources to confirm the identification of a SCUBA-2-detected ultraluminous star-forming galaxy in this structure. Integrating the total star-formation activity in the central region of the structure, we estimate that it is an order of magnitude higher (in a mass-normalised sense) than clusters at z~0.5-1. However, we also find that the most active cluster members do not reside in the densest regions of the structure, which instead host a population of passive and massive, red galaxies. We suggest that while the passive and active populations have comparable near-infrared luminosities at z=1.6, M(H)~-23, the subsequent stronger fading of the more active galaxies means that they will evolve into passive systems at the present-day which are less luminous than the descendants of those galaxies which were already passive at z~1.6 (M(H)~-20.5 and M(H)~-21.5 respectively at z~0). We conclude that the massive galaxy population in the dense cores of present-day clusters were already in place at z=1.6 and that in Cl0218.3-0510 we are seeing continuing infall of less extreme, but still ultraluminous, star-forming galaxies onto a pre-existing structure.

We study the sub-mm properties of color-selected galaxies via a stacking analysis applied for the first time to interferometric data at sub-mm wavelengths. We base our study on 344 GHz ALMA continuum observations of ~20''-wide fields centered on 86 sub-mm sources detected in the LABOCA Extended Chandra Deep Field South Sub-mm Survey (LESS). We select various classes of galaxies (K-selected, star-forming sBzK galaxies, extremely red objects and distant red galaxies) according to their optical/NIR fluxes. We find clear, >10-sigma detections in the stacked images of all these galaxy classes. We include in our stacking analysis Herschel/SPIRE data to constrain the dust SED of these galaxies. We find that their dust emission is well described by a modified black body with Tdust~30 K and beta=1.6 and IR luminosities of (5-11)x10^{11} Lsun, or implied star formation rates of 75-140 Msun/yr. We compare our results with those of previous studies based on single-dish observations at 870 micron and find that our flux densities are a factor 2-3 higher than previous estimates. The discrepancy is observed also after removing sources individually detected in ALESS maps. We report a similar discrepancy by repeating our analysis on 1.4\,GHz observations of the whole ECDFS. Hence we find tentative evidence that galaxies that are associated in projected and redshift space with sub-mm bright sources are brighter than the average population. Finally, we put our findings in the context of the cosmic star formation rate density as a function of redshift.

The large gas and dust reservoirs of submm galaxies (SMGs) could potentially provide ample fuel to trigger an Active Galactic Nucleus (AGN), but previous studies of the AGN fraction in SMGs have been controversial largely due to the inhomogeneity and limited angular resolution of the available submillimeter surveys. Here we set improved constraints on the AGN fraction and X-ray properties of the SMGs with ALMA and Chandra observations in the Extended Chandra Deep Field-South (E-CDF-S). This study is the first among similar works to have unambiguously identified the X-ray counterparts of SMGs; this is accomplished using the fully submm-identified, statistically reliable SMG catalog with 99 SMGs from the ALMA LABOCA E-CDF-S Submillimeter Survey (ALESS). We found 10 X-ray sources associated with SMGs (median redshift z = 2.3), of which 8 were identified as AGNs using several techniques that enable cross-checking. The other 2 X-ray detected SMGs have levels of X-ray emission that can be plausibly explained by their star-formation activity. 6 of the 8 SMG-AGNs are moderately/highly absorbed, with N_H > 10e23 cm-2. An analysis of the AGN fraction, taking into account the spatial variation of X-ray sensitivity, yields an AGN fraction of 17+16-6% for AGNs with rest-frame 0.5-8 keV absorption-corrected luminosity >7.8x10e42 erg s-1; we provide estimated AGN fractions as a function of X-ray flux and luminosity. ALMA's high angular resolution also enables direct X-ray stacking at the precise positions of SMGs for the first time, and we found 4 potential SMG-AGNs in our stacking sample.

We investigate the multi-wavelength properties of a sample of 450-\mu m selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4\sigma\ in deep SCUBA-2 450-\mu m observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of \sigma 450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450\mu m-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <\beta_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between \beta_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 \mu m sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

We estimate the galaxy stellar mass function and stellar mass density for star-forming and quiescent galaxies with 0.2<z<4. We construct a deep K<24 sample of 220000 galaxies selected using the UltraVISTA DR1 data release. Our analysis is based on precise 30-band photometric redshifts. By comparing these photometric redshifts with 10800 spectroscopic redshifts from the zCOSMOS bright and faint surveys, we find a precision of sigma(dz/(1+z))=0.008 at i<22.5 and sigma(dz/(1+zs))=0.03 at 1.5<z<4. We derive the stellar mass function and correct for the Eddington bias. We find a mass-dependent evolution of the global and star-forming populations. This mass-dependent evolution is a direct consequence of the star formation being quenched in galaxies more massive than M>10^10.7Msun. For the mass function of the quiescent galaxies, we do not find any significant evolution of the high-mass end at z<1; however we observe a clear flattening of the faint-end slope. From z~3 to z~1, the density of quiescent galaxies increases over the entire mass range. Their comoving stellar mass density increases by 1.6 dex between z~3 and z~1 and by less than 0.2dex at z<1. We infer the star formation history from the mass density evolution and we find an excellent agreement with instantaneous star formation rate measurements at z<1.5, while we find differences of 0.2dex at z>1.5 consistent with the expected uncertainties. We also develop a new method to infer the specific star formation rate from the mass function of star-forming galaxies. We find that the specific star formation rate of 10^10Msun galaxies increases continuously in the redshift range 1<z<4. Finally, we compare our results with a semi-analytical model and find that these models overestimate the density of low mass quiescent galaxies by an order of magnitude, while the density of low-mass star-forming galaxies is successfully reproduced.

We present an Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 survey of 126 submillimeter sources from the LABOCA ECDFS Submillimeter Survey (LESS). Our 870 micron survey with ALMA (ALESS) has produced maps ~3X deeper and with a beam area ~200X smaller than the original LESS observations, doubling the current number of interferometrically-observed submillimeter sources. The high resolution of these maps allows us to resolve sources that were previously blended and accurately identify the origin of the submillimeter emission. We discuss the creation of the ALESS submillimeter galaxy (SMG) catalog, including the main sample of 99 SMGs and a supplementary sample of 32 SMGs. We find that at least 35% (possibly up to 50%) of the detected LABOCA sources have been resolved into multiple SMGs, and that the average number of SMGs per LESS source increases with LESS flux density. Using the (now precisely known) SMG positions, we empirically test the theoretical expectation for the uncertainty in the single-dish source positions. We also compare our catalog to the previously predicted radio/mid-infrared counterparts, finding that 45% of the ALESS SMGs were missed by this method. Our ~1.6" resolution allows us to measure a size of ~9 kpc X 5 kpc for the rest-frame ~300 um emission region in one resolved SMG, implying a star formation rate surface density of 80 M_sol yr^-1 kpc^-2, and we constrain the emission regions in the remaining SMGs to be <10 kpc. As the first statistically reliable survey of SMGs, this will provide the basis for an unbiased multiwavelength study of SMG properties.

We present CO(2-1) observations of the submillimetre galaxy ALESS65.1 performed with the Australia Telescope Compact Array at 42.3 GHz. A previous ALMA study of submillimetre galaxies in the Extended Chandra Deep Field South detected [CII] 157.74 micron emission from this galaxy at a redshift of z = 4.44. No CO(2-1) emission was detected but we derive a firm upper limit to the cold gas mass in ALESS65.1 of M_gas < 1.7 x 10^10 M_odot. The estimated gas depletion timescale is <50 Myr, which is similar to other high redshift SMGs, and consistent with z > 4 SMGs being the likely progenitors of massive red-and-dead galaxies at z > 2. The ratio of the [CII], CO and far-infrared luminosities implies a strong far-ultraviolet field of G_0 > 10^3, as seen in Galactic star forming regions or local ULIRGs. The observed L_[CII]/L_FIR = 2.3 x 10^{-3} is high compared to local ULIRGs and, combined with L_[CII]/L_CO > 2700, it is consistent with ALESS65.1 either having an extended (several kpc) [CII] emitting region or lower than solar metallicity.

We report the first counts of faint submillimetre galaxies (SMG) in the 870-um band derived from arcsecond resolution observations with the Atacama Large Millimeter Array (ALMA). We have used ALMA to map a sample of 122 870-um-selected submillimetre sources drawn from the (0.5x0.5)deg^2 LABOCA Extended Chandra Deep Field South Submillimetre Survey (LESS). These ALMA maps have an average depth of sigma(870um)~0.4mJy, some ~3x deeper than the original LABOCA survey and critically the angular resolution is more than an order of magnitude higher, FWHM of ~1.5" compared to ~19" for the LABOCA discovery map. This combination of sensitivity and resolution allows us to precisely pin-point the SMGs contributing to the submillimetre sources from the LABOCA map, free from the effects of confusion. We show that our ALMA-derived SMG counts broadly agree with the submillimetre source counts from previous, lower-resolution single-dish surveys, demonstrating that the bulk of the submillimetre sources are not caused by blending of unresolved SMGs. The difficulty which well-constrained theoretical models have in reproducing the high-surface densities of SMGs, thus remains. However, our observations do show that all of the very brightest sources in the LESS sample, S(870um)>12mJy, comprise emission from multiple, fainter SMGs, each with 870-um fluxes of <9mJy. This implies a natural limit to the star-formation rate in SMGs of <10^3 M_Sun/yr, which in turn suggests that the space densities of z>1 galaxies with gas masses in excess of ~5x10^10 M_Sun is <10^-5 Mpc^-3. We also discuss the influence of this blending on the identification and characterisation of the SMG counterparts to these bright submillimetre sources and suggest that it may be responsible for previous claims that they lie at higher redshifts than fainter SMGs.

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

We present ALMA 870-um (345GHz) observations of two sub-millimetre galaxies (SMGs) drawn from an ALMA study of the 126 sub-millimeter sources from the LABOCA Extended Chandra Deep Field South Survey (LESS). The ALMA data identify the counterparts to these previously unidentified sub-millimeter sources and serendipitously detect bright emission lines in their spectra which we show are most likely to be [C II]157.74um emission yielding redshifts of z=4.42 and z=4.44. This blind detection rate within the 7.5-GHz bandpass of ALMA is consistent with the previously derived photometric redshift distribution of SMGs and suggests a modest, but not dominant (<25%), tail of 870-um selected SMGs at z>4. We find that the ratio of L_CII/L_FIR in these SMGs is much higher than seen for similarly far-infrared-luminous galaxies at z~0, which is attributed to the more extended gas reservoirs in these high-redshift ULIRGs. Indeed, in one system we show that the [C II] emission shows hints of extended emission on >3kpc scales. Finally, we use the volume probed by our ALMA survey to show that the bright end of the [C II] luminosity function evolves strongly between z=0 and z~4.4, reflecting the increased ISM cooling in galaxies as a result of their higher star-formation rates. These observations demonstrate that even with short integrations, ALMA is able to detect the dominant fine structure cooling lines from high-redshift ULIRGs, measure their energetics and trace their evolution with redshift.

We present Herschel-PACS spectroscopy of the [OI]63um far-infrared cooling line from a sample of six unlensed and spectroscopically-confirmed 870um-selected submillimetre (submm) galaxies (SMGs) at 1.1<z<1.6 from the LABOCA Extended Chandra Deep Field South (ECDFS) Submm Survey (LESS). This is the first survey of [OI]63um, one of the main photodissociation region (PDR) cooling lines, in SMGs. New high-resolution ALMA interferometric 870um continuum imaging confirms that these six Herschel-targeted SMG counterparts are bona fide sources of submm emission. We detect [OI]63um in two SMGs with a SNR >3, tentatively detect [OI]63um in one SMG, and constrain the line flux for the non-detections. We also exploit the combination of submm continuum photometry from 250-870um and our new PACS continuum measurements to constrain the far-infrared (FIR) luminosity, L_FIR, in these SMGs to < 30%. We find that SMGs do not show a deficit in their [OI]63um-to-far-infrared continuum luminosity ratios (with ratios ranging from ~0.5-1.5%), similar to what was seen previously for the [CII]158um-to-FIR ratios in SMGs. These observed ratios are about an order of magnitude higher than what is seen typically for local ultra luminous infrared galaxies (ULIRGs), which adds to the growing body of evidence that SMGs are not simply `scaled up' versions of local ULIRGs. Rather, the PDR line-to-L_FIR ratios suggest that the star formation modes of SMGs are likely more akin to that of local normal (lower-luminosity) star-forming galaxies, with the bulk of the star formation occurring in extended regions, galaxy-scale (~kpc) in size. These observations represent the first step towards a census of the major PDR cooling lines in typical SMGs that will be attainable with ALMA, enabling detailed modelling to probe the global properties of the star formation and the evolutionary status of SMGs.

We present new IRAM PdBI 1.3mm continuum observations at ~1.5" resolution of 28 SMGs previously discovered with the 870um bolometer LABOCA at APEX within the central 0.7deg2 of the COSMOS field. 19 out of the 28 LABOCA sources were detected with the PdBI at a >~3sigma level of ~1.4mJy/b. A combined analysis of this new sample with existing interferometrically identified SMGs in the COSMOS field yields the following results: 1) >~15%, and possibly up to ~40% of single-dish detected SMGs consist of multiple sources, 2) statistical identifications of multi-wavelength counterparts to the single-dish SMGs yield that only ~50% of these single-dish SMGs have real radio or IR counterparts, 3) ~18% of interferometric SMGs have only radio or even no multi-wavelength counterpart at all, and 4) ~50-70% of z>~3 SMGs have no radio counterparts down to an rms of 7-12uJy at 1.4GHz. Using the exact interferometric positions to identify proper multi-wavelength counterparts allows us to determine accurate photometric redshifts for these sources. The redshift distributions of the combined and the individual 1.1mm and 870um selected samples have a higher mean and broader width than the redshift distributions derived in previous studies. Our sample supports the previous tentative trend that on average brighter and/or mm-selected SMGs are located at higher redshifts. There is a tentative offset between the mean redshift for the 1.1mm (<z>=3.1+/-0.4) and 870um (<z>=2.6+/-0.4) selected samples, with the 1.1mm sources lying on average at higher redshifts. Based on our nearly complete sample of AzTEC 1.1mm SMGs within a uniform 0.15deg2 area we infer a higher surface density of z>~4 SMGs than predicted by current cosmological models. In summary, our findings imply that (sub-)millimeter interferometric identifications are crucial to build statistically complete and unbiased samples of SMGs.

We report on interferometric observations at 1.3 mm at 2"-3" resolution using the Combined Array for Research in Millimeter-wave Astronomy (CARMA). We identify multi-wavelength counterparts of three submillimeter galaxies (SMGs; F(1mm)>5.5 mJy) in the COSMOS field, initially detected with MAMBO and AzTEC bolometers at low, ~10"-30", resolution. All three sources -- AzTEC/C1, Cosbo-3 and Cosbo-8 -- are identified to coincide with positions of 20 cm radio sources. Cosbo-3, however, is not associated with the most likely radio counterpart, closest to the MAMBO source position, but that further away from it. This illustrates the need for intermediate-resolution (~2") mm-observations to identify the correct counterparts of single-dish detected SMGs. All of our three sources become prominent only at NIR wavelengths, and their mm-to-radio flux based redshifts suggest that they lie at redshifts z>~2. As a proof of concept, we show that photometric redshifts can be well determined for SMGs, and we find photometric-redshifts of 5.6+/-1.2, 1.9+0.9(-0.5), and ~4 for AzTEC/C1, Cosbo-3, and Cosbo-8, respectively. Using these we infer that these galaxies have radio-based star formation rates of >~1000 Msol/yr, and IR luminosities of ~10^13 Lsol consistent with properties of high-redshift SMGs. In summary, our sources reflect a variety of SMG properties in terms of redshift and clustering, consistent with the framework that SMGs are progenitors of z~2 and today's passive galaxies.

Analyses of high-redshift ultraluminous infrared (IR) galaxies traditionally use the observed optical to submillimeter spectral energy distribution (SED) and estimates of the dynamical mass as observational constraints to derive the star formation rate (SFR), the stellar mass, and age of these objects. An important observational constraint neglected in the analysis is the mass of dust giving rise to the IR emission. In this paper we add this constraint to the analysis of AzTEC-3. Adopting an upper limit to the mass of stars and a bolometric luminosity for this object, we construct stellar and chemical evolutionary scenarios, constrained to produce the inferred dust mass and observed luminosity before the associated stellar mass exceeds the observational limit. We find that the model with a Top Heavy IMF provided the most plausible scenario consistent with the observational constraints. In this scenario the dust formed over a period of ~200 Myr, with a SFR of ~500 Msun/yr. These values for the age and SFR in AzTEC-3 are significantly higher and lower, respectively, from those derived without the dust mass constraint. However, this scenario is not unique, and others cannot be completely ruled out because of the prevailing uncertainties in the age of the galaxy, its bolometric luminosity, and its stellar and dust masses. A robust result of our models is that all scenarios require most of the radiating dust mass to have been accreted in molecular clouds. Our new procedure highlights the importance of a multiwavelength approach, and of the use of dust evolution models in constraining the age and the star formation activity and history in galaxies.

We present results of our 1.4 GHz image stacking analysis of mass-selected galaxies in the COSMOS field. From the resulting median radio continuum flux density we have determined the evolution of the average star formation rate (SFR) of galaxies as a function of stellar mass, unbiased from effects of dust but also source confusion due to the 1.5" angular resolution achieved by the VLA. We find a power-law relation between specific SFR (SSFR) and stellar mass for star forming galaxies out to z=3. While higher mass systems exhibit lower SSFRs at any epoch, no differential, more rapid evolution of high mass galaxies is evident at least out to z~1.5 where our conclusions are most robust. Utilizing measured mass functions of star forming systems, the characteristic stellar mass for galaxies contributing most to the comoving SFR density appears not to evolve. These findings hence challenge 'downsizing' scenarios in which star formation has gradually shifted towards lower mass systems with cosmic time. Our analysis constitutes the to-date best determination of the cosmic star formation history (CSFH) since z=3 and yields indirect evidence for a rapid decline of the global mass density of molecular gas with time.

Massive clusters of galaxies have been found as early as 3.9 Billion years (z=1.62) after the Big Bang containing stars that formed at even earlier epochs. Cosmological simulations using the current cold dark matter paradigm predict these systems should descend from "proto-clusters" - early over-densities of massive galaxies that merge hierarchically to form a cluster. These proto-cluster regions themselves are built-up hierarchically and so are expected to contain extremely massive galaxies which can be observed as luminous quasars and starbursts. However, observational evidence for this scenario is sparse due to the fact that high-redshift proto-clusters are rare and difficult to observe. Here we report a proto-cluster region 1 billion years (z=5.3) after the Big Bang. This cluster of massive galaxies extends over >13 Mega-parsecs, contains a luminous quasar as well as a system rich in molecular gas. These massive galaxies place a lower limit of >4x10^11 solar masses of dark and luminous matter in this region consistent with that expected from cosmological simulations for the earliest galaxy clusters.

We explore the evolution of the specific star formation rate (SSFR) for 3.6um-selected galaxies of different M_* in the COSMOS field. The average SFR for sub-sets of these galaxies is estimated with stacked 1.4GHz radio continuum emission. We separately consider the total sample and a subset of galaxies (SF) that shows evidence for substantive recent star formation in the rest-frame optical SED. At 0.2<z<3 both populations show a strong and M_*-independent decrease in their SSFR towards z=0.2, best described by a power- law (1+z)^n, where n~4.3 for all galaxies and n~3.5 for SF sources. The decrease appears to have started at z>2, at least above 4x10^10M_Sun where our conclusions are most robust. We find a tight correlation with power-law dependence, SSFR (M_*)^beta, between SSFR and M_* at all z. It tends to flatten below ~10^10M_Sun if quiescent galaxies are included; if they are excluded a shallow index beta_SFG -0.4 fits the correlation. On average, higher M_* objects always have lower SSFRs, also among SF galaxies. At z>1.5 there is tentative evidence for an upper SSFR-limit that an average galaxy cannot exceed. It is suggested by a flattening of the SSFR-M_* relation (also for SF sources), but affects massive (>10^10M_Sun) galaxies only at the highest z. Below z=1.5 there thus is no direct evidence that galaxies of higher M_* experience a more rapid waning of their SSFR than lower M_* SF systems. In this sense, the data rule out any strong 'downsizing'. We combine our results with recent measurements of the galaxy (stellar) mass function in order to determine the characteristic mass of a SF galaxy (M_*=10^(10.6\pm0.4)M_Sun). In this sense, too, there is no 'downsizing'. Our analysis constitutes the most extensive SFR density determination with a single technique to z=3. Recent Herschel results are consistent with our results, but rely on far smaller samples.

We present observations at 1.2 mm with MAMBO-II of a sample of z>~2 radio-intermediate obscured quasars, as well as CO observations of two sources with the Plateau de Bure Interferometer. Five out of 21 sources (24%) are detected at a significance of >=3sigma. Stacking all sources leads to a statistical detection of <S_1.2mm>= 0.96+-0.11 mJy and stacking only the non-detections also yields a statistical detection, with <S_1.2mm>= 0.51+-0.13 mJy. This corresponds to a typical far-infrared luminosity L_FIR~4x10^12 Lsol. If the far-infrared luminosity is powered entirely by star-formation, and not by AGN-heated dust, then the characteristic inferred star-formation rate is ~700 Msol yr-1. This far-infrared luminosity implies a dust mass of M_dust~3x10^8 Msol. We estimate that such large dust masses on kpc scales can plausibly cause the obscuration of the quasars. We present dust SEDs for our sample and derive a mean SED for our sample. This mean SED is not well fitted by clumpy torus models, unless additional extinction and far-infrared re-emission due to cool dust are included. There is a hint that the host galaxies of obscured quasars must have higher far-infrared luminosities and cool-dust masses and are therefore often found at an earlier evolutionary phase than those of unobscured quasars. For one source at z=2.767, we detect the CO(3-2) transition, with S_CO Delta nu=630+-50 mJy km s-1, corresponding to L_CO(3-2)= 3.2x10^7 Lsol, or L'_CO(3-2)=2.4x10^10 K km s-1 pc2. For another source at z=4.17, the lack of detection of the CO(4-3) line yields a limit of L'_CO(4-3)<1x10^10 K km s-1 pc2. Molecular gas masses, gas depletion timescales and gas-to-dust ratios are estimated (Abridged).