Search SciRate

We use the Atacama Large sub/Millimetre Array (ALMA) to efficiently observe spectral lines across Bands 3, 4, 5, 6, 7, and 8 at high-resolution (0.5" - 0.1") for 16 bright southern Herschel sources at $1.5 < z < 4.2$. With only six and a half hours of observations, we reveal 66 spectral lines in 17 galaxies. These observations detect emission from CO (3-2) to CO(18-17), as well as atomic ([CI](1-0), (2-1), [OI] 145 $\mu$m and [NII] 205 $\mu$m) lines. Additional molecular lines are seen in emission (${\rm H_2O}$ and ${\rm H_2O^+}$) and absorption (OH$^+$ and CH$^+$). The morphologies based on dust continuum ranges from extended sources to strong lensed galaxies with magnifications between 2 and 30. CO line transitions indicate a diverse set of excitation conditions with a fraction of the sources ($\sim 35$%) showcasing dense, warm gas. The resolved gas to star-formation surface densities vary strongly per source, and suggest that the observed diversity of dusty star-forming galaxies could be a combination of lensed, compact dusty starbursts and extended, potentially-merging galaxies. The predicted gas depletion timescales are consistent with 100 Myr to 1 Gyr, but require efficient fueling from the extended gas reservoirs onto the more central starbursts, in line with the Doppler-shifted absorption lines that indicate inflowing gas for two out of six sources. This pilot paper explores a successful new method of observing spectral lines in large samples of galaxies, supports future studies of larger samples, and finds that the efficiency of this new observational method will be further improved with the planned ALMA Wideband Sensitivity Upgrade.

High-redshift protoclusters consisting of dusty starbursts are thought to play an important role in galaxy evolution. Their dusty nature makes them bright in the FIR/submm but difficult to find in optical/NIR surveys. Radio observations are an excellent way to study these dusty starbursts, as dust is transparent in the radio and there is a tight correlation between the FIR and radio emission of a galaxy. Here, we present MeerKAT 1.28 GHz radio imaging of 3 Herschel candidate protoclusters, with a synthesised beam size of ~$7.5'' \times 6.6''$ and a central thermal noise down to $4.35~\mu$Jy/beam. Our source counts are consistent with other radio counts with no evidence of overdensities. Around 95% of the Herschel sources have 1.28 GHz IDs. Using the Herschel 250 micron primary beam size as the searching radius, we find 54.2% Herschel sources have multiple 1.28 GHz IDs. Our average FIR-radio correlation coefficient $q_{250\mu m}$ is $2.33\pm 0.26$. Adding $q_{250\mu m}$ as a new constraint, the probability of finding chance-aligned sources is reduced by a factor of ~6, but with the risk of discarding true identifications of radio-loud/quiet sources. With accurate MeerKAT positions, we cross-match our Herschel sources to optical/NIR data followed by photometric redshift estimations. By removing $z<1$ sources, the density contrasts of two of the candidate protoclusters increase, suggestive of them being real protoclusters at $z>1$. There is also potentially a $0.9<z<1.2$ overdensity associated with one candidate protocluster. In summary, photometric redshifts from radio-optical cross-identifications have provided some tentative evidence of overdensities aligning with two of the candidate protoclusters.

The ALMACAL survey is based on a database of reprocessed ALMA calibration scans suitable for scientific analysis, observed as part of regular PI observations. We present all the data accumulated from the start of ALMA operations until May 2022 for 1047 calibrator fields across the southern sky spanning ALMA Bands 3 to 10 (~ 84 - 950 GHz), so-called ALMACAL-22. Encompassing over 1000 square arcmin and accumulating over 2000 hours of integration time, ALMACAL is not only one of the largest ALMA surveys to date, but it continues to grow with each new scientific observation. We outline the methods for processing and imaging a subset of the highest-quality data ('pruned sample'). Using deconvolution techniques within the visibility data (uv plane), we created data cubes as the final product for further scientific analysis. We describe the properties and shortcomings of ALMACAL and compare its area and sensitivity with other sub-millimetre surveys. Notably, ALMACAL overcomes limitations of previous sub-millimetre surveys, such as small sky coverage and the effects of cosmic variance. Moreover, we discuss the improvements introduced by the latest version of this dataset that will enhance our understanding of dusty star-forming galaxies, extragalactic absorption lines, active galactic nucleus physics, and ultimately the evolution of molecular gas.

In this Roadmap, we present a vision for the future of submillimetre and millimetre astronomy in the United Kingdom over the next decade and beyond. This Roadmap has been developed in response to the recommendation of the Astronomy Advisory Panel (AAP) of the STFC in the AAP Astronomy Roadmap 2022. In order to develop our stragetic priorities and recommendations, we surveyed the UK submillimetre and millimetre community to determine their key priorities for both the near-term and long-term future of the field. We further performed detailed reviews of UK leadership in submillimetre/millimetre science and instrumentation. Our key strategic priorities are as follows: 1. The UK must be a key partner in the forthcoming AtLAST telescope, for which it is essential that the UK remains a key partner in the JCMT in the intermediate term. 2. The UK must maintain, and if possible enhance, access to ALMA and aim to lead parts of instrument development for ALMA2040. Our strategic priorities complement one another: AtLAST (a 50m single-dish telescope) and an upgraded ALMA (a large configurable interferometric array) would be in synergy, not competition, with one another. Both have identified and are working towards the same overarching science goals, and both are required in order to fully address these goals.

Magnetic fields are prevalent on almost all astronomical scales, but their importance in different systems and over cosmic time is yet to be understood. Our current knowledge of the evolution of magnetic fields is limited by scarce observations in the distant Universe, where galaxies have recently been found to be more evolved than most model predictions. In this study, we conducted rest-frame 131 um full-polarisation observations of dust emission in a strongly lensed dusty star-forming galaxy, SPT0346-52, at z=5.6, when the Universe was only 1 Gyr old. Dust grains can become aligned with local magnetic fields, resulting in the emission of linearly polarised thermal infrared radiation. Our observations have revealed a median polarisation level of $0.9\pm0.2$ per cent with a variation of $\pm0.4$ per cent across the regions with polarisation detection, similar to that of local starburst galaxies. The polarised dust emission is patchy and mostly overlaps with the [C II] emission at a velocity of about -150 km/s, and extends over 3 kiloparsecs. The bimodal distribution of field orientations, their spatial distribution, and the connection with the cold gas kinematics further emphasise the complexity of the magnetic environment in this galaxy and the potential role of mergers in shaping its magnetic fields. Such early formation of ordered galactic magnetic fields also suggests that both small-scale and large-scale dynamos could be efficient in early galaxies. Continued observations of magnetic fields in early galaxies, as well as expanding surveys to a wider galaxy population, are essential for a comprehensive understanding of the prevalence and impact of magnetic fields in the evolving Universe.

Protoclusters at $z>2$ are gas-rich regions characterized by high star-formation activity. The same physical properties that enhance star formation in protoclusters are also thought to boost the growth of SMBHs. We aim to test this scenario by probing the AGN content of SPT2349-56, a massive, gas-rich, and highly star-forming protocluster core at $z=4.3$ discovered as an overdensity of DSFGs, via Chandra (200 ks) observations, and comparing the results with the field environment. We detected two protocluster members, corresponding to an AGN fraction among DSFGs of $\approx10\%$. This value is consistent with other protoclusters at $z=2-4$, but higher than the AGN incidence among DSFGs in the field environment. Both AGN are heavily obscured sources and hosted in star-forming galaxies with $\approx3\times10^{10}\,\mathrm{M_\odot}$ stellar masses. We estimate that the ISM in the host galaxies can contribute significantly to the nuclear obscuration. One of the two AGN is highly luminous ($L_X=2\times10^{45}\,\mathrm{erg\,s^{-1}}$) and Compton-thick ($N_H=2\times10^{24}\,\mathrm{cm^{-2}}$), and likely powered by a $M_{BH}>6\times10^8\,\mathrm{M_\odot}$ SMBH. Its high accretion rate suggests that it is in the phase of efficient growth required to explain the presence of extremely massive SMBHs in the centers of local galaxy clusters. Considering SPT2349-56 and DRC, a similar protocuster at $z=4$, we find that gas-rich protocluster cores at $z\approx4$ enhance the triggering of luminous (log$\frac{L_X}{\mathrm{erg\,s^{-1}}}=45-46$) AGN by 3-5 orders of magnitude with respect to the field environment. Our results indicate that gas-rich protoclusters at high redshift boost the growth of SMBHs, which will likely impact the subsequent evolution of the structures, and thus represent key science targets to obtain a complete understanding of the relation between environment and galaxy evolution.

Recent observations and theories have presented a strong challenge to the universality of the stellar initial mass function (IMF) in extreme environments. A notable example has been found for starburst conditions, where evidence favours a top-heavy IMF, i.e. there is a bias toward massive stars compared to the IMF that is responsible for the stellar mass function and elemental abundances observed in the Milky Way. Local starburst galaxies have star-formation rates similar to those in high-redshift main-sequence galaxies, which appear to dominate the stellar mass budget at early epochs. However, the IMF of high-redshift main-sequence galaxies is yet to be probed. Since $^{13}$CO and C$^{18}$O isotopologues are sensitive to the IMF, we have observed these lines towards four strongly-lensed high-redshift main-sequence galaxies using the Atacama Large Millimeter/sub-millimeter Array. Of our four targets, SDSS J0901+1814, at $z \approx 2.26$, is seen clearly in $^{13}$CO and C$^{18}$O, the first detection of CO isotopologues in the high-redshift main-sequence galaxy population. The observed $^{13}$C/$^{18}$O ratio, $2.4 \pm 0.8$, is significantly lower than that of local main-sequence galaxies. We estimate the isotope ratio, oxygen abundance and stellar mass using a series of chemical evolution models with varying star-formation histories and IMFs. All models favour an IMF that is more top-heavy than that of the Milky Way. Thus, as with starburst galaxies, main-sequence galaxies in the high-redshift Universe have a greater fraction of massive stars than a Milky-Way IMF would imply.

We report on deep SCUBA-2 observations at 850$\mu$m and NOEMA spectroscopic measurements at 2 mm of the environment surrounding the luminous, massive ($M_{*} \approx 2 \times 10^{11}$ M$_{\odot}$) Herschel-selected source HerBS-70. This source was revealed by previous NOEMA observations to be a binary system of dusty star-forming galaxies at $z= 2.3$, with the East component (HerBS-70E) hosting an Active Galactic Nucleus (AGN). The SCUBA-2 observations detected, in addition to the binary system, twenty-one sources at $> 3.5 \sigma$ over an area of $\sim 25$ square comoving Mpc with a sensitivity of $\sigma_{850} = 0.75$ mJy. The surface density of continuum sources around HerBS-70 is three times higher than for field galaxies. The NOEMA spectroscopic measurements confirm the protocluster membership of three of the nine brightest sources through their CO(4 - 3) line emission, yielding a volume density 36 times higher than for field galaxies. All five confirmed sub-mm galaxies in the HerBS-70 system have relatively short gas depletion times ($80 - 500$ Myr), indicating the onset of quenching for this protocluster core due to the depletion of gas. The dark matter halo mass of the HerBS-70 system is estimated around $5 \times{} 10^{13}$ M$_{\odot}$, with a projected current-day mass of $10^{15}$ M$_{\odot}$, similar to the local Virgo and Coma clusters. These observations support the claim that DSFGs, in particular the ones with observed multiplicity, can trace cosmic overdensities.

The CO(1--0) and [\ionC1](1--0) emission lines are well-established tracers of cold molecular gas mass in local galaxies. At high redshift, where the interstellar medium (ISM) is likely to be denser, there have been limited direct comparisons of both ground state transitions. Here we present a study of CO(1--0) and [\ionC1](1--0) emission in a sample of 20 unlensed dusty, star-forming galaxies at $z=2-5$. The CO(1--0)/[\ionC1](1--0) ratio is constant up to at least $z=5$, supporting the use of [CI](1-0) as a gas mass tracer. PDR modelling of the available data indicates a median H$_2$ density of log$(n~[$cm$^{-3}])=4.7\pm0.2$, and UV radiation field log$(G_{\mathrm{UV}} [G$_0$])=3.2\pm0.2$. We use the CO(1--0), [\ionC1](1--0) and 3mm dust continuum measurements to cross--calibrate the respective gas mass conversion factors, finding no dependence of these factors on either redshift or infrared luminosity. Assuming a variable CO conversion factor then implies [\ionC1] and dust conversion factors that differ from canonically assumed values but are consistent with the solar/super-solar metallicities expected for our sources. Radiative transfer modelling shows that the warmer CMB at high redshift can significantly affect the [\ionC1] as well as CO emission, which can change the derived molecular gas masses by up to 70\% for the coldest kinetic gas temperatures expected. Nevertheless, we show that the magnitude of the effect on the ratio of the tracers is within the known scatter of the $L'_\mathrm{CO}-L'_\mathrm{[CI]}$ relation. Further determining the absolute decrease of individual line intensities will require well-sampled spectral line energy distributions (SLEDs) to model the gas excitation conditions in more detail.

Deep, unbiased surveys are essential to decipher the cosmic evolution of galaxies. The submillimetre (submm) and millimetre (mm) windows complement the UV/optical waveband and are key to revealing the cold and dusty Universe. Traditional ways of conducting deep surveys resort to either lensed fields or target small areas for ultra-long integrations. These surveys have greatly advanced our understanding of dusty star-forming galaxies (DSFGs), but are susceptible to lensing uncertainties and cosmic variance and will be expensive to expand. Here, we summarise our recent multi-wavelength survey of DSFGs in the vicinity of ALMA's calibrators: the ALMACAL survey. These fields have accumulated many hundreds of hours of on-source time, reaching depths and effective areas that are competitive with bespoke cosmological surveys. We summarise the multi-wavelength number counts from ALMACAL and the resolved fraction of the Cosmic Infrared Background (CIB) from submm to mm wavelengths. Meanwhile, combining all available ALMA observations in each field results in impressive frequency coverage, which often yields the redshifts of these DSFGs. The ALMACAL survey has demonstrated the scientific value of calibration scans for all submm/mm and radio telescopes, existing and planned.

During the most active period of star formation in galaxies, which occurs in the redshift range 1<z<3, strong bursts of star formation result in significant quantities of dust, which obscures new stars being formed as their UV/optical light is absorbed and then re-emitted in the infrared, which redshifts into the mm/sub-mm bands for these early times. To get a complete picture of the high-z galaxy population, we need to survey a large patch of the sky in the sub-mm with sufficient angular resolution to resolve all galaxies, but we also need the depth to fully sample their cosmic evolution, and therefore obtain their redshifts using direct mm spectroscopy with a very wide frequency coverage. This requires a large single-dish sub-mm telescope with fast mapping speeds at high sensitivity and angular resolution, a large bandwidth with good spectral resolution and multiplex spectroscopic capabilities. The proposed 50-m Atacama Large Aperture Submillimeter Telescope (AtLAST) will deliver these specifications. We discuss how AtLAST allows us to study the whole population of high-z galaxies, including the dusty star-forming ones which can only be detected and studied in the sub-mm, and obtain a wealth of information for each of these up to z~7: gas content, cooling budget, star formation rate, dust mass, and dust temperature. We present worked examples of surveys that AtLAST can perform, both deep and wide, and also focused on galaxies in proto-clusters. In addition we show how such surveys with AtLAST can measure the growth rate and the Hubble constant with high accuracy, and demonstrate the power of the line-intensity mapping method in the mm/sub-mm wavebands to constrain the cosmic expansion history at high redshifts, as good examples of what can uniquely be done by AtLAST in this research field.

Variations in the dust emissivity index, $\beta$, within and between galaxies, are evidence that the chemistry and physics of dust must vary on large scales, although the nature of the physical and/or chemical variations is still unknown. In this paper we estimate values of $\beta$ and dust temperature for a sample of 109 dusty star-forming galaxies (DSFGs) over the range, $2 < z < 6$. We compare the results obtained with both an optically-thin model and a general opacity model, finding that our estimates of $\beta$ are similar between the models but our estimates of dust temperature are not. We find no evidence of a change in $\beta$ with redshift, with a median value of $\beta = 1.96$ for the optically-thin model with a confidence interval (16 - 84%) of 1.67 to 2.35 for the population. Using simulations, we estimate the measurement errors from our procedure and show that the variation of $\beta$ in the population results from intrinsic variations in the properties of the dust in DSFGs. At a fixed far-infrared luminosity, we find no evidence for a change in dust temperature, $T_\textrm{dust}$, with redshift. After allowing for the effects of correlated measurement errors, we find an inverse correlation between $\beta$ and $T_\textrm{dust}$ in DSFGs, for which there is also evidence in low-redshift galaxies.

Magnetic fields are fundamental to the evolution of galaxies, playing a key role in the astrophysics of the interstellar medium and star formation. Large-scale ordered magnetic fields have been mapped in the Milky Way and nearby galaxies, but it is not known how early in the Universe such structures form. Here we report the detection of linearly polarized thermal emission from dust grains in a strongly lensed, intrinsically luminous galaxy that is forming stars at a rate more than a thousand times that of the Milky Way at redshift 2.6, within 2.5 Gyr of the Big Bang. The polarized emission arises from the alignment of dust grains with the local magnetic field. The median polarization fraction is of order one per cent, similar to nearby spiral galaxies. Our observations support the presence of a 5 kiloparsec-scale ordered magnetic field with a strength of around 500uG or lower, orientated parallel to the molecular gas disk. This confirms that such structures can be rapidly formed in galaxies, early in cosmic history.

Understanding the nature of high-$z$ dusty galaxies requires a comprehensive view of their ISM and molecular complexity. However, the molecular ISM at high-$z$ is commonly studied using only a few species beyond CO, limiting our understanding. In this paper, we present the results of deep 3 mm spectral line surveys using the NOEMA targeting two lensed dusty galaxies: APM 08279+5255 (APM), a quasar at redshift $z=3.911$, and NCv1.143 (NC), a $z=3.565$ starburst galaxy. The spectral line surveys cover rest-frame frequencies from about 330-550 GHz. We report the detection of 38 and 25 emission lines in APM and NC, respectively. The spectra reveal the chemical richness and the complexity of the physical properties of the ISM. By comparing the spectra of the two sources and combining the gas excitation analysis, we find that the physical properties and the chemical imprints of the ISM are different between them: the molecular gas is more excited in APM, exhibiting higher molecular-gas temperatures and densities compared to NC; the chemical abundances in APM are akin to the values of local AGN, showing boosted relative abundances of the dense gas tracers that might be related to high-temperature chemistry and/or XDRs, while NC more closely resembles local starburst galaxies. The most significant differences are found in H2O, where the 448GHz H2O line is significantly brighter in APM, likely linked to the intense far-infrared radiation from the dust powered by AGN. Our astrochemical model suggests that at such high column densities, FUV radiation is less important in regulating the ISM, while CRs (X-rays/shocks) are the key players in shaping the abundance of the molecules and the initial conditions of star formation. Such deep spectral line surveys open a new window to study the physical and chemical properties of the ISM and the radiation field of galaxies in the early Universe. (abridged)

(Abridged) Using the IRAM NOEMA interferometer, we measures the redshifts of 126 bright galaxies detected in the Herschel H-ATLAS, HeLMS, and HerS surveys. We report reliable spectroscopic redshifts for a total of 124 of the Herschel-selected galaxies. The redshifts are estimated from scans of the 3 and 2-mm bands (and, in one case, the 1-mm band) and are based on the detection of at least two emission lines. Together with the Pilot Programme (Neri et al. 2020), including spectroscopic redshifts of 11 sources, our survey has derived precise redshifts for 135 bright Herschel-selected galaxies, making it the largest sample of high-z galaxies with robust redshifts to date. Most emission lines detected are from 12CO (mainly from J=2-1 to 5-4), with some sources seen in [CI] and H2O emission lines. The spectroscopic redshifts are in the range 0.8<z<6.55 with a median value of z=2.56 +/- 0.10. The line widths of the sources are large, with a mean value for the full width at half maximum Delta(V) of 590 +/- 25 km/s and with 35% of the sources having widths of 700 km/s < Delta(V) < 1800 km/s. Most of the sources are unresolved or barely resolved on scales of 2 to 3 arcsec (or linear sizes of 15-25 kpc, unlensed). Some fields reveal double or multiple sources and, in some cases, sources at different redshifts. Taking these sources into account, there are, in total, 165 individual sources with robust spectroscopic redshifts, including lensed galaxies, binary systems, and over-densities. We present an overview of the z-GAL survey and provide the observed properties of the emission lines, the derived spectroscopic redshifts, and an atlas of the entire sample. The data presented here will serve as a foundation for the other z-GAL papers in this series reporting on the dust emission, the molecular and atomic gas properties, and a detailed analysis of the nature of the sources.

(Abridged) We present the dust properties of 125 bright Herschel galaxies selected from the z-GAL survey. The large instantaneous bandwidth of NOEMA provides an exquisite sampling of the underlying dust continuum emission at 2 and 3 mm in the observed frame, with flux densities in at least four side bands for each source. Together with the available Herschel 250, 350, and 500 micron and SCUBA-2 850 micron flux densities, the spectral energy distribution of each source can be analyzed from the far-infrared to the millimeter, with a fine sampling of the Rayleigh-Jeans tail. This wealth of data provides a solid basis to derive robust dust properties, in particular the dust emissivity index, beta, and the dust temperature, T(dust). In order to demonstrate our ability to constrain the dust properties, we used a flux-generated mock catalog and analyzed the results under the assumption of an optically thin and optically thick modified black body emission. For the z-GAL sources, we report a range of dust emissivities with beta ~ 1.5 - 3 estimated up to high precision with relative uncertainties that vary in the range 7% - 15%, and an average of 2.2 +/- 0.3. We find dust temperatures varying from 20 to 50 K with an average of T(dust) ~ 30 K for the optically thin case and ~38 K in the optically thick case. For all the sources, we estimate the dust masses and apparent infrared luminosities (based on the optically thin approach). An inverse correlation is found between T(dust) and beta, which is similar to what is seen in the local Universe. Finally, we report an increasing trend in the dust temperature as a function of redshift at a rate of 6.5 +/- 0.5 K/z for this 500 micron-selected sample. Based on this study, future prospects are outlined to further explore the evolution of dust temperature across cosmic time.

The z-GAL survey observed 137 bright Herschel-selected targets with the IRAM NOrthern Extended Millimeter Array, with the aim to measure their redshift and study their properties. Several of them have been resolved into multiple sources. Consequently, robust spectroscopic redshifts have been measured for 165 individual galaxies in the range 0.8<z<6.5. In this paper we analyse the millimetre spectra of the z-GAL sources, using both their continuum and line emission to derive their physical properties. At least two spectral lines are detected for each source, including transitions of 12CO, [CI], and H2O. The observed 12CO line ratios and spectral line energy distributions of individual sources resemble those of local starbursts. In seven sources the para-H2O(2_11-2_02) transition is detected and follows the IR versus H2O luminosity relation of sub-millimetre galaxies. The molecular gas mass of the z-GAL sources is derived from their 12CO, [CI], and sub-millimetre dust continuum emission. The three tracers lead to consistent results, with the dust continuum showing the largest scatter when compared to 12CO. The gas-to-dust mass ratio of these sources was computed by combining the information derived from 12CO and the dust continuum and has a median value of 107, similar to star-forming galaxies of near-solar metallicity. The same combined analysis leads to depletion timescales in the range between 0.1 and 1.0 Gyr, which place the z-GAL sources between the `main sequence' of star formation and the locus of starbursts. Finally, we derived a first estimate of stellar masses - modulo possible gravitational magnification - by inverting known gas scaling relations: the z-GAL sample is confirmed to be mostly composed by starbursts, whereas ~25% of its members lie on the main sequence of star-forming galaxies (within +/- 0.5 dex).

It may be unsurprising that the most common approach to finding proto-clusters is to search for over-densities of galaxies. Upgrades to submillimetre (submm) interferometers and the advent of the James Webb Space Telescope will soon offer the opportunity to find more distant candidate proto-clusters in deep sky surveys without any spectroscopic confirmation. In this letter, we report the serendipitous discovery of an extremely dense region centred on the blazar, J0217-0820, at z=0.6 in the ALMACAL sky survey. Its density is eight times higher than that predicted by blind submm surveys. Among the seven submm-bright galaxies, three are as bright as conventional single-dish submm galaxies, with S_870um > 3mJy. The over-density is thus comparable to the densest known and confirmed proto-cluster cores. However, their spectra betray a wide range of redshifts. We investigate the likelihood of line-of-sight projection effects using light cones from cosmological simulations, finding that the deeper we search, the higher the chance that we will suffer from such projection effects. The extreme over-density around J0217-0820 demonstrates the strong cosmic variance we may encounter in the deep submm surveys. Thus, we should also question the fidelity of galaxy proto-cluster candidates selected via over-densities of galaxies, where the negative K correction eases the detection of dusty galaxies along an extraordinarily extended line of sight.

We analyse the molecular and atomic emission lines of 71 bright Herschel-selected galaxies between redshifts 1.4 to 4.6 detected by the Atacama Large Millimetre/submillimetre Array. These lines include a total of 156 CO, [C I], and H2O emission lines. For 46 galaxies, we detect two transitions of CO lines, and for these galaxies we find gas properties similar to those of other dusty star-forming galaxy (DSFG) samples. A comparison to photo-dissociation models suggests that most of Herschel-selected galaxies have similar interstellar medium conditions as local infrared-luminous galaxies and high-redshift DSFGs, although with denser gas and more intense far-ultraviolet radiation fields than normal star-forming galaxies. The line luminosities agree with the luminosity scaling relations across five orders of magnitude, although the star-formation and gas surface density distributions (i.e., Schmidt-Kennicutt relation) suggest a different star-formation phase in our galaxies (and other DSFGs) compared to local and low-redshift gas-rich, normal star-forming systems. The gas-to-dust ratios of these galaxies are similar to Milky Way values, with no apparent redshift evolution. Four of 46 sources appear to have CO line ratios in excess of the expected maximum (thermalized) profile, suggesting a rare phase in the evolution of DSFGs. Finally, we create a deep stacked spectrum over a wide rest-frame frequency (220-890 GHz) that reveals faint transitions from HCN and CH, in line with previous stacking experiments.

We present a statistical study of 180 dust continuum sources identified in 33 massive cluster fields by the ALMA Lensing Cluster Survey (ALCS) over a total of 133 arcmin$^{2}$ area, homogeneously observed at 1.2 mm. ALCS enables us to detect extremely faint mm sources by lensing magnification, including near-infrared (NIR) dark objects showing no counterparts in existing \it Hubble Space Telescope and \it Spitzer images. The dust continuum sources belong to a blind sample ($N=141$) with S/N $\gtrsim$ 5.0 (a purity of $>$ 0.99) or a secondary sample ($N=39$) with S/N= $4.0-5.0$ screened by priors. With the blind sample, we securely derive 1.2-mm number counts down to $\sim7$ $\mu$Jy, and find that the total integrated 1.2mm flux is 20.7$^{+8.5}_{-6.5}$ Jy deg$^{-2}$, resolving $\simeq$ 80 % of the cosmic infrared background light. The resolved fraction varies by a factor of $0.6-1.1$ due to the completeness correction depending on the spatial size of the mm emission. We also derive infrared (IR) luminosity functions (LFs) at $z=0.6-7.5$ with the $1/V_{\rm max}$ method, finding the redshift evolution of IR LFs characterized by positive luminosity and negative density evolution. The total (=UV+IR) cosmic star-formation rate density (SFRD) at $z>4$ is estimated to be $161^{+25}_{-21}$ % of the established measurements, which were almost exclusively based on optical$-$NIR surveys. Although our general understanding of the cosmic SFRD is unlikely to change beyond a factor of 2, these results add to the weight of evidence for an additional ($\approx 60$ %) SFRD component contributed by the faint-mm population, including NIR dark objects.

We present initial results of an ongoing survey with the Karl G. Jansky Very Large Array targeting the CO($J$ = 1-0) transition in a sample of 30 submillimeter-selected, dusty star-forming galaxies at $z =$ 2-5 with existing mid--$J$ CO detections from ALMA and NOEMA, of which 17 have been fully observed. We detect CO(1-0) emission in 11 targets, along with three tentative ($\sim$1.5-2$\sigma$) detections; three galaxies are undetected. Our results yield total molecular gas masses of 6-23$\times$10$^{10}$ ($\alpha_\mathrm{CO}$/1) M$_\odot$, with gas mass fractions, $f_\mathrm{gas}$=$M_\mathrm{mol}$/($M_*$+$M_\mathrm{mol}$), of 0.1-0.8 and a median depletion time of (140$\pm$70) Myr. We find median CO excitation ratios of $r_{31}$ = 0.75$\pm$0.39 and $r_{41}$ = 0.63$\pm$0.44, with a significant scatter. We find no significant correlation between the excitation ratio and a number of key parameters such as redshift, CO(1-0) line width or $\Sigma_\mathrm{SFR}$. We only find a tentative positive correlation between $r_{41}$ and the star-forming efficiency, but we are limited by our small sample size. Finally, we compare our results to predictions from the SHARK semi-analytical model, finding a good agreement between the molecular gas masses, depletion times and gas fractions of our sources and their SHARK counterparts. Our results highlight the heterogeneous nature of the most massive star-forming galaxies at high-redshift, and the importance of CO(1--0) observations to robustly constrain their total molecular gas content and ISM properties.

We present the results of a survey of CO(1-0) emission in 14 infrared luminous dusty star forming galaxies (DSFGs) at 2 < z < 4 with the NSF's Karl G. Jansky Very Large Array. All sources are detected in CO(1-0), with an ~1arcsec angular resolution. Seven sources show extended and complex structure. We measure CO luminosities of $({\mu})L'_{CO(1-0)}=0.4-2.9x10^{11}$ K km s$^{-1}$ pc$^2$, and molecular gas masses of (${\mu}$)M$_{H2}$ = 1.3 - 8.6 x 10$^{11}$ Mo, where (\mu) is the magnification factor. The derived molecular gas depletion times of t$_{\rm dep}$ = 40 - 460 Myr, cover the expected range of both normal star forming galaxies and starbursts. Comparing to the higher-J CO transitions previously observed for the same sources, we find CO temperature brightness ratios of r$_{32/10}$ = 0.4 - 1.4, r$_{43/10}$ = 0.4 - 1.7, and r$_{54/10}$ = 0.3 - 1.3. We find a wide range of CO spectral line energy distributions (SLEDs), in agreement with other high-z DSFGs, with the exception of three sources that are most comparable to the Cloverleaf and APM08279+5255. Based on radiative transfer modelling of the CO SLEDs we determine densities of n$_{H2}$ = 0.3 - 8.5 x 10$^3$ cm$^{-3}$ and temperatures of T$_K$ = 100 - 200 K. Lastly, four sources are detected in the continuum, three have radio emission consistent with their infrared derived star formation rates, while HerBS-70E requires an additional synchrotron radiation component from an active galactic nucleus. Overall, we find that even though the sample is similarly luminous in the infrared, by tracing the CO(1-0) emission a diversity of galaxy and excitation properties are revealed, demonstrating the importance of CO(1-0) observations in combination to higher-J transitions.

We present a pilot, untargeted extragalactic carbon monoxide (CO) emission-line survey using ALMACAL, a project utilizing ALMA calibration data for scientific purposes. In 33 deep (Texp > 40 min) ALMACAL fields we report six CO emission-line detections above S/N > 4, one-third confirmed by MUSE observations. With this pilot survey, we probe a cosmologically significant volume of ~10^5 cMpc^3, widely distributed over many pointings in the southern sky, making the survey largely insusceptible to the effects of cosmic variance. We derive the redshift probability of the CO detections using probability functions from the Shark semi-analytical model of galaxy formation. By assuming typical CO excitations for the detections, we put constraints on the cosmic molecular gas mass density evolution over the redshift range 0 < z < 1.5. The results of our pilot survey are consistent with the findings of other untargeted emission-line surveys and the theoretical model predictions and currently cannot rule out a non-evolving molecular gas mass density. Our study demonstrates the potential of using ALMA calibrator fields as a multi-sightline untargeted CO emission line survey. Applying this approach to the full ALMACAL database will provide an accurate, free of cosmic variance, measurement of the molecular luminosity function as a function of redshift.

Wide, deep, blind continuum surveys at submillimetre/millimetre (submm/mm) wavelengths are required to provide a full inventory of the dusty, distant Universe. However, conducting such surveys to the necessary depth, with sub-arcsec angular resolution, is prohibitively time-consuming, even for the most advanced submm/mm telescopes. Here, we report the most recent results from the ALMACAL project, which exploits the 'free' calibration data from the Atacama Large Millimetre/submillimetre Array (ALMA) to map the lines of sight towards and beyond the ALMA calibrators. ALMACAL has now covered 1,001 calibrators, with a total sky coverage around 0.3 deg2, distributed across the sky accessible from the Atacama desert, and has accumulated more than 1,000h of integration. The depth reached by combining multiple visits to each field makes ALMACAL capable of searching for faint, dusty, star-forming galaxies (DSFGs), with detections at multiple frequencies to constrain the emission mechanism. Based on the most up-to-date ALMACAL database, we report the detection of 186 DSFGs with flux densities down to S870um ~ 0.2mJy, comparable with existing ALMA large surveys but less susceptible to cosmic variance. We report the number counts at five wavelengths between 870um and 3mm, in ALMA bands 3, 4, 5, 6 and 7, providing a benchmark for models of galaxy formation and evolution. By integrating the observed number counts and the best-fitting functions, we also present the resolved fraction of the cosmic infrared background (CIB) and the CIB spectral shape. Combining existing surveys, ALMA has currently resolved about half of the CIB in the submm/mm regime.

We report the detection of the ground state rotational emission of ammonia, ortho-NH$_3$ $(J_K=1_0\rightarrow0_0)$ in a gravitationally lensed, intrinsically hyperluminous, star-bursting galaxy at $z=2.6$. The integrated line profile is consistent with other molecular and atomic emission lines which have resolved kinematics well-modelled by a 5 kpc-diametre rotating disc. This implies that the gas responsible for NH$_3$ emission is broadly tracing the global molecular reservoir, but likely distributed in pockets of high density ($n\gtrsim5\times10^4$ cm$^{-3}$). With a luminosity of $2.8\times10^{6}$ $L_\odot$, the NH$_3$ emission represents $2.5\times10^{-7}$ of the total infrared luminosity of the galaxy, comparable to the ratio observed in the Kleinmann-Low nebula in Orion and consistent with sites of massive star formation in the Milky Way. If $L_{\rm NH_3}/L_{\rm IR}$ serves as a proxy for the 'mode' of star formation, this hints that the nature of star formation in extreme starbursts in the early Universe is similar to that of Galactic star-forming regions, with a large fraction of the cold interstellar medium in this state, plausibly driven by a storm of violent disc instabilities in the gas-dominated disc. This supports the 'full of Orions' picture of star formation in the most extreme galaxies seen close to the peak epoch of stellar mass assembly.

We present a self-consistent cross-calibration of the three main molecular gas mass tracers in galaxies, the $\rm ^{12}CO$(1-0), [CI]($^3P_1$-$^3P_0$) lines, and the submm dust continuum emission, using a sample of 407 galaxies, ranging from local disks to submillimetre-selected galaxies (SMGs) up to $z \approx 6$. A Bayesian method is used to produce galaxy-scale universal calibrations of these molecular gas indicators, that hold over 3-4 orders of magnitude in infrared luminosity, $L_{\rm IR}$. Regarding the dust continuum, we use a mass-weighted dust temperature, $T_{\rm mw}$, determined using new empirical relations between temperature and luminosity. We find the average $L/M_{\rm mol}$ gas mass conversion factors to be $\alpha_{850}= 6.9\times10^{12}\,\rm W\,Hz^{-1}\,M_{\odot}^{-1}$, $\alpha_{\rm CO} = \rm 4\,M_{\odot} (K\,km\,s^{-1}\,pc^2)^{-1}$ and $\alpha_{\rm CI} = \rm 17.0 \,M_{\odot} (K\,km\,s^{-1}\,pc^2)^{-1}$, based on the assumption that the mean dust properties of the sample ($\kappa_H$ = gas-to-dust ratio/dust emissivity) will be similar to those of local metal rich galaxies and the MW. The tracer with the least intrinsic scatter is [CI](1-0), while CO(1-0) has the highest. The conversion factors show a weak but significant correlation with $L_{\rm IR}$. Assuming dust properties typical of metal-rich galaxies, we infer a neutral carbon abundance $X_{\rm CI} = [C^0/\rm mol]=1.6\times 10^{-5}$, similar to that in the MW. We find no evidence for bimodality of $\alpha_{\rm CO}$ between main-sequence (MS) galaxies and those with extreme star-formation intensity, i.e. ULIRGs and SMGs. The means of the three conversion factors are found to be similar between MS galaxies and ULIRGs/SMGs, to within 10-20%. We show that for metal-rich galaxies, near-universal average values for $\alpha_{\rm CO}$, $X_{\rm CI}$ and $\kappa_H$ are adequate for global molecular gas estimates.

We examine the far-IR properties of a sample of 5391 optically selected QSOs in the 0.5<z<2.65 redshift range down to log[nuLnu,2500 (erg/s)]>44.7, using SPIRE data from Herschel-ATLAS. We split the sample in a grid of 74 luminosity-redshift bins and compute the average optical-infrared spectral energy distribution (SED) in each bin. By normalising an intrinsic AGN template to the AGN optical power (at 5100A) we decompose the total infrared emission (L_IR; 8-1000um) into an AGN (L_IR,AGN) and star-forming component (L_IR,SF). We find that the AGN contribution to L_IR increases as a function of AGN power which manifests as a reduction of the `far-IR bump' in the average QSO SEDs. We note that L_IR,SF does not correlate with AGN power; the mean star formation rates (SFRs) of AGN host galaxies are a function of redshift only and they range from ~6 Msun/yr at z~0 to a plateau of <200 Msun/yr at z~2.6. Our results indicate that the accuracy of far-IR emission as a proxy for SFR decreases with increasing AGN luminosity. We show that, at any given redshift, observed trends between infrared luminosity (whether monochromatic or total) and AGN power (in the optical or X-rays) can be explained by a simple model which is the sum of two components: (A) the infrared emission from star-formation, uncorrelated with AGN power and (B) the infrared emission from AGN, directly proportional to AGN power in the optical or X-rays.

While unobscured and radio-quiet active galactic nuclei are regularly being found at redshifts $z > 6$, their obscured and radio-loud counterparts remain elusive. We build upon our successful pilot study, presenting a new sample of low-frequency-selected candidate high-redshift radio galaxies (HzRGs) over a sky area twenty times larger. We have refined our selection technique, in which we select sources with curved radio spectra between 72-231 MHz from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. In combination with the requirements that our GLEAM-selected HzRG candidates have compact radio morphologies and be undetected in near-infrared $K_{\rm s}$-band imaging from the Visible and Infrared Survey Telescope for Astronomy Kilo-degree Infrared Galaxy (VIKING) survey, we find 51 new candidate HzRGs over a sky area of approximately 1200 deg$^2$. Our sample also includes two sources from the pilot study: the second-most distant radio galaxy currently known, at $z=5.55$, with another source potentially at $z \sim 8$. We present our refined selection technique and analyse the properties of the sample. We model the broadband radio spectra between 74 MHz and 9 GHz by supplementing the GLEAM data with both publicly available data and new observations from the Australia Telescope Compact Array at 5.5 and 9 GHz. In addition, deep $K_{\rm s}$-band imaging from the High-Acuity Widefield $K$-band Imager (HAWK-I) on the Very Large Telescope and from the Southern Herschel Astrophysical Terahertz Large Area Survey Regions $K_{\rm s}$-band Survey (SHARKS) is presented for five sources. We discuss the prospects of finding very distant radio galaxies in our sample, potentially within the epoch of reionisation at $z \gtrsim 6.5$.

In this paper we present the third data release (DR3) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We identify likely near-infrared counterparts to submillimetre sources in the South Galactic Pole (SGP) field using the VISTA VIKING survey. We search for the most probable counterparts within 15 arcsec of each Herschel source using a probability measure based on the ratio between the likelihood the true counterpart is found close to the submillimetre source and the likelihood that an unrelated object is found in the same location. For 110 374 (57.0$\%$) sources we find galaxies on the near-infrared images where the probability that the galaxy is associated to the source is greater than 0.8. We estimate the false identification rate to be 4.8$\%$, with a probability that the source has an associated counterpart on the VIKING images of 0.835$\pm$0.009. We investigate the effects of gravitational lensing and present 41 (0.14 deg$^{-2}$) candidate lensed systems with observed flux densities > 100 mJy at 500 $\mu$m. We include in the data release a probability that each source is gravitationally lensed and discover an additional 5 923 sources below 100 mJy that have a probability greater than 0.94 of being gravitationally lensed. We estimate that $\sim$ 400 - 1 000 sources have multiple true identifications in VIKING based on the similarity of redshift estimates for multiple counterparts close to a Herschel source. The data described in this paper can be found at the H-ATLAS website.

We present spectroscopic measurements for 71 galaxies associated with 62 of the brightest high-redshift submillimeter sources from the Southern fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), while targeting 85 sources which resolved into 142. We have obtained robust redshift measurements for all sources using the 12-m Array and an efficient tuning of ALMA to optimise its use as a redshift hunter, with 73 per cent of the sources having a robust redshift identification. Nine of these redshift identifications also rely on observations from the Atacama Compact Array. The spectroscopic redshifts span a range $1.41<z<4.53$ with a mean value of 2.75, and the CO emission line full-width at half-maxima range between $\rm 110\,km\,s^{-1} < FWHM < 1290\,km\,s^{-1}$ with a mean value of $\sim$ 500kms$^{-1}$, in line with other high-$z$ samples. The derived CO(1-0) luminosity is significantly elevated relative to line-width to CO(1-0) luminosity scaling relation, which is suggestive of lensing magnification across our sources. In fact, the distribution of magnification factors inferred from the CO equivalent widths is consistent with expectations from galaxy-galaxy lensing models, though there is a hint of an excess at large magnifications that may be attributable to the additional lensing optical depth from galaxy groups or clusters.

We have obtained 1 and 3 mm spectral scans of ASXDF1100.053.1 using NOEMA. ASXDF1100.053.1 is an unlensed optically dark millimetre-bright SMG with $K_{\rm AB}>25.7$ ($2\sigma$), which was expected to lie at $z=$5-7 based on its radio-submm photo-$z$. Our data detected line emission due to $^{12}$CO($J=$5-4) and ($J=$6-5), providing a $z_{\rm CO}= 5.2383\pm0.0005$. Energy-coupled SED modelling indicates properties of $L_{\rm IR}=8.3^{+1.5}_{-1.4}\times10^{12}$ L$_{\odot}$, SFR $=630^{+260}_{-380}$ M$_{\odot}$ yr$^{-1}$, $M_{\rm dust}=4.4^{+0.4}_{-0.3}\times10^{8}$ M$_{\odot}$, $M_{\rm stellar}=3.5^{+3.6}_{-1.4}\times10^{11}$ M$_{\odot}$, and $T_{\rm dust}=37.4^{+2.3}_{-1.8}$ K. The CO luminosity allows us to estimate a gas mass $M_{\rm gas}=3.1\pm0.3\times10^{10}$ M$_{\odot}$, suggesting a gas-to-dust mass ratio of around 70, fairly typical for $z\sim2$ SMGs. ASXDF1100.053.1 has $R_{\rm e, mm}=1.0^{+0.2}_{-0.1}$ kpc, so its surface $L_{\rm IR}$ density $\Sigma_{\rm IR}$ is $1.2^{+0.1}_{-0.2}\times10^{12}$ L$_{\odot}$ kpc$^{-2}$. These properties indicate that ASXDF1100.053.1 is a massive dusty star-forming (SF) galaxy with an unusually compact starburst. It lies close to the SF main sequence at $z\sim5$, with low $M_{\rm gas}$/$M_{\rm stellar}=0.09$, SFR/SFR$_{\rm MS} (R_{\rm SB})=0.6$, and a gas-depletion time $\tau_{\rm dep}$ of $\approx 50$ Myr, modulo assumptions about the stellar initial mass function in such objects. ASXDF1100.053.1 has extreme values of $M_{\rm gas}/M_{\rm stellar}$, $R_{\rm SB}$, and $\tau_{\rm dep}$ compared to SMGs at $z\sim$2-4, and those of ASXDF1100.053.1 are the smallest among SMGs at $z>5$. ASXDF1100.053.1 is likely a late-stage dusty starburst prior to passivisation. The number of $z=$5.1-5.3 unlensed SMGs now suggests a number density $dN/dz=30.4\pm19.0$ deg$^{-2}$, barely consistent with the latest cosmological simulations.

We introduce an ALMA band 3 spectroscopic survey, targeting the brightest submillimeter galaxies (SMGs) in the COSMOS field. Here we present the first results based on the 18 primary SMGs that have 870 $\mu$m flux densities of $S_{870}=12.4-19.3$ mJy and are drawn from a parent sample of 260 ALMA-detected SMGs from the AS2COSMOS survey. We detect emission lines in 17 and determine their redshifts to be in the range of $z=2-5$ with a median of ${3.3\pm0.3}$. We confirm that SMGs with brighter $S_{870}$ are located at higher redshifts. The data additionally cover five fainter companion SMGs, and we obtain line detection in one. Together with previous studies, our results indicate that for SMGs that satisfy our selection, their brightest companion SMGs are physically associated with their corresponding primary SMGs in $\ge40$% of the time, suggesting that mergers play a role in the triggering of star formation. By modeling the foreground gravitational fields, $<10$% of the primary SMGs can be strongly lensed with a magnification $\mu>2$. We determine that about 90\% of the primary SMGs have lines that are better described by double Gaussian profiles, and the median separation of the two Gaussian peaks is 430$\pm$40 km s$^{-1}$. This allows estimates of an average baryon mass, which together with the line dispersion measurements puts our primary SMGs on the similar mass-$\sigma$ correlation found on local early-type galaxies. Finally, the number density of our $z>4$ primary SMGs is found to be $1^{+0.9}_{-0.6}\times10^6$ cMpc$^{-3}$, suggesting that they can be the progenitors of $z\sim3-4$ massive quiescent galaxies.

We have carried out the first spatially-resolved investigation of the multi-phase interstellar medium (ISM) at high redshift, using the z=4.24 strongly-lensed sub-millimetre galaxy H-ATLASJ142413.9+022303 (ID141). We present high-resolution (down to ~350 pc) ALMA observations in dust continuum emission and in the CO(7-6), H_2O (2_1,1 - 2_0,2), CI(1-0) and CI(2-1) lines, the latter two allowing us to spatially resolve the cool phase of the ISM for the first time. Our modelling of the kinematics reveals that the system appears to be dominated by a rotationally-supported gas disk with evidence of a nearby perturber. We find that the CI(1-0) line has a very different distribution to the other lines, showing the existence of a reservoir of cool gas that might have been missed in studies of other galaxies. We have estimated the mass of the ISM using four different tracers, always obtaining an estimate in the range (3.2-3.8) x 10^11 M_sol, significantly higher than our dynamical mass estimate of (0.8-1.3) x 10^11 M_sol. We suggest that this conflict and other similar conflicts reported in the literature is because the gas-to-tracer ratios are ~4 times lower than the Galactic values used to calibrate the ISM in high-redshift galaxies. We demonstrate that this could result from a top-heavy initial mass function and strong chemical evolution. Using a variety of quantitative indicators, we show that, extreme though it is at z=4.24, ID141 will likely join the population of quiescent galaxies that appears in the Universe at z~3.

We present modelling of ~0.1arcsec resolution Atacama Large Millimetre/sub-millimeter Array imaging of seven strong gravitationally lensed galaxies detected by the Herschel Space Observatory. Four of these systems are galaxy-galaxy scale strong lenses, with the remaining three being group-scale lenses. Through careful modelling of visibilities, we infer the mass profiles of the lensing galaxies and by determining the magnification factors, we investigate the intrinsic properties and morphologies of the lensed sub-millimetre sources. We find that these sub-millimetre sources all have ratios of star formation rate to dust mass that is consistent with or in excess of the mean ratio for high-redshift sub-millimetre galaxies and low redshift ultra-luminous infrared galaxies. The contribution to the infrared luminosity from possible AGN is not quantified and so could be biasing our star formation rates to higher values. The majority of our lens models have mass density slopes close to isothermal, but some systems show significant differences.

We present observations of a giant Lyman-alpha blob in the SSA22 proto-cluster at z=3.1, SSA22-LAB1, taken with the Atacama Large Millimeter/submillimeter Array (ALMA). Dust continuum, along with [C II]158um, and CO(4-3) line emission have been detected in LAB1, showing complex morphology and kinematics across a ~100 kpc central region. Seven galaxies at z=3.0987-3.1016 in the surroundings are identified in [C II] and dust continuum emission, with two of them potential companions or tidal structures associated with the most massive galaxies. Spatially resolved [C II] and infrared luminosity ratios for the widely distributed media (L[C II]/LIR~0.01-0.001) suggest that the observed extended interstellar media are likely to have originated from star-formation activity and the contribution from shocked gas is probably not dominant. LAB1 is found to harbour a total molecular gas mass Mmol=(8.7+/-2.0)e+10 Msun, concentrated in the core region of the Ly-alpha-emitting area. While (primarily obscured) star-formation activity in the LAB1 core is one of the most plausible power sources for the Ly-alpha emission, multiple major-mergers found in the core may also play a role in making LAB1 exceptionally bright and extended in Ly-alpha as a result of cooling radiation induced by gravitational interactions.

We present LMT/AzTEC 1.1mm observations of $\sim100$ luminous high-redshift dusty star-forming galaxy candidates from the $\sim600\,$sq.deg $Herschel$-ATLAS survey, selected on the basis of their SPIRE red far-infrared colours and with $S_{500\mu\rm m}=35-80$ mJy. With an effective $\theta_{\rm FWHM}\approx9.5\,$ arcsec angular resolution, our observations reveal that at least 9 per cent of the targets break into multiple systems with SNR $\geq 4$ members. The fraction of multiple systems increases to $\sim23\,$ per cent (or more) if some non-detected targets are considered multiples, as suggested by the data. Combining the new AzTEC and deblended $Herschel$ photometry we derive photometric redshifts, IR luminosities, and star formation rates. While the median redshifts of the multiple and single systems are similar $(z_{\rm med}\approx3.6)$, the redshift distribution of the latter is skewed towards higher redshifts. Of the AzTEC sources $\sim85\,$ per cent lie at $z_{\rm phot}>3$ while $\sim33\,$ per cent are at $z_{\rm phot}>4$. This corresponds to a lower limit on the space density of ultra-red sources at $4<z<6$ of $\sim3\times10^{-7}\, \textrm{Mpc}^{-3}$ with a contribution to the obscured star-formation of $\gtrsim 8\times10^{-4}\, \textrm{M}_\odot \textrm{yr}^{-1} \textrm{Mpc}^{-3}$. Some of the multiple systems have members with photometric redshifts consistent among them suggesting possible physical associations. Given their angular separations, these systems are most likely galaxy over-densities and/or early-stage pre-coalescence mergers. Finally, we present 3mm LMT/RSR spectroscopic redshifts of six red-$Herschel$ galaxies at $z_{\rm spec}=3.85-6.03$, two of them (at $z \sim 4.7$) representing new redshift confirmations. Here we release the AzTEC and deblended $Herschel$ photometry as well as catalogues of the most promising interacting systems and $z>4$ galaxies.

We present a comparative analysis of the $\rm CH^+$(1-0) and $\rm Ly \alpha$ lines, observed with the Atacama Large Millimeter Array (ALMA) and Keck telescope respectively, in the field of the submillimetre-selected galaxy (SMG) SMM\u2009J02399$-$0136 at $z\sim2.8$, which comprises a heavily obscured starburst galaxy and a broad absorption line quasar, immersed in a large $\rm Ly \alpha$ nebula. This comparison highlights the critical role played by turbulence in channeling the energy across gas phases and scales, splitting the energy trail between hot/thermal and cool/turbulent phases in the circum-galactic medium (CGM). The unique chemical and spectroscopic properties of $\rm CH^+$ are used to infer the existence of a massive ($\sim 3.5 \times 10^{10}$ ${\rm M}_\odot$), highly turbulent reservoir of diffuse molecular gas of radius $\sim 20\,$kpc coinciding with the core of the $\rm Ly \alpha$ nebula. The whole cool and cold CGM is shown to be inflowing towards the galaxies at a velocity $\sim$ 400 km$\,s^{-1}$. Several kpc-scale shocks are detected tentatively in $\rm CH^+$ emission. Their specific location in space and velocity with respect to the high-velocity $\rm Ly \alpha$ emission suggests that they lie at the interface of the inflowing CGM and the high-velocity $\rm Ly \alpha$ emission, and signpost the feeding of CGM turbulence by AGN- and stellar-driven outflows. The mass and energy budgets of the CGM require net mass accretion at a rate commensurate with the star formation rate (SFR). From this similarity, we infer that the merger-driven burst of star formation and black-hole growth are ultimately fuelled by large-scale gas accretion.

We present new empirical constraints on the evolution of $\rho_{\rm H_2}$, the cosmological mass density of molecular hydrogen, back to $z\approx2.5$. We employ a statistical approach measuring the average observed $850\mu{\rm m}$ flux density of near-infrared selected galaxies as a function of redshift. The redshift range considered corresponds to a span where the $850\mu{\rm m}$ band probes the Rayleigh-Jeans tail of thermal dust emission in the rest-frame, and can therefore be used as an estimate of the mass of the interstellar medium (ISM). Our sample comprises of ${\approx}150,000$ galaxies in the UKIDSS-UDS field with near-infrared magnitudes $K_{\rm AB}\leq25$ mag and photometric redshifts with corresponding probability distribution functions derived from deep 12-band photometry. With a sample approximately 2 orders of magnitude larger than in previous works we significantly reduce statistical uncertainties on $\rho_{\rm H_2}$ to $z\approx2.5$. Our measurements are in broad agreement with recent direct estimates from blank field molecular gas surveys, finding that the epoch of molecular gas coincides with the peak epoch of star formation with $\rho_{\rm H_2}\approx2\times10^7\,{\rm M_\odot}\,{\rm Mpc^{-3}}$ at $z\approx2$. We demonstrate that $\rho_{\rm H_2}$ can be broadly modelled by inverting the star-formation rate density with a fixed or weakly evolving star-formation efficiency. This 'constant efficiency' model shows a similar evolution to our statistically derived $\rho_{\rm H_2}$, indicating that the dominant factor driving the peak star formation history at $z\approx2$ is a larger supply of molecular gas in galaxies rather than a significant evolution of the star-formation rate efficiency within individual galaxies.

We present an analysis of the spatial clustering of a large sample of high-resolution, interferometically identified, submillimetre galaxies (SMGs). We measure the projected cross-correlation function of ~350 SMGs in the UKIDSS Ultra Deep-Survey Field across a redshift range of $z=1.5-3$ utilising a method that incorporates the uncertainties in the redshift measurements for both the SMGs and cross-correlated galaxies through sampling their full probability distribution functions. By measuring the absolute linear bias of the SMGs we derive halo masses of $\log_{10}(M_{\rm halo}[{h^{-1}\,\rm M_{\odot}}])\sim12.8$ with no evidence of evolution in the halo masses with redshift, contrary to some previous work. From considering models of halo mass growth rates we predict that the SMGs will reside in haloes of mass $\log_{10}(M_{\rm halo}[{h^{-1}\,\rm M_{\odot}}])\sim13.2$ at $z=0$, consistent with the expectation that the majority of $z=1.5-3$ SMGs will evolve into present-day spheroidal galaxies. Finally, comparing to models of stellar-to-halo mass ratios, we show that SMGs may correspond to systems that are maximally efficient at converting their gas reservoirs into stars. We compare them to a simple model for gas cooling in halos that suggests that the unique properties of the SMG population, including their high levels of star-formation and their redshift distribution, are a result of the SMGs being the most massive galaxies that are still able to accrete cool gas from their surrounding intragalactic medium.

(Abridged) Exploiting the sensitivity and broad band width of NOEMA, we have studied the molecular gas and dust in the galaxy HerBS-89a, at z=2.95. High angular resolution images reveal a partial 1.0" diameter Einstein ring in the dust continuum emission and the molecular emission lines of 12CO(9-8) and H2O(2_02-1_11). We report the detection of the three fundamental transitions of the molecular ion OH+, seen in absorption; the molecular ion CH+(1-0) seen in absorption (and tentatively in emission); two transitions of amidogen (NH2), seen in emission; and HCN(11-10) and/or NH(1_2-0_1) seen in absorption. The NOEMA data are complemented with VLA data tracing the 12CO(1-0) emission line, which provides a measurement of the total mass of molecular gas and an anchor for a CO excitation analysis. In addition, we present HST imaging that reveals the foreground lensing galaxy in the near-infrared. Together with data from the GTC, we derive a photometric redshift of z(phot)~0.9 for the foreground lensing galaxy. Modelling the lensing of HerBS-89a, we reconstruct the dust continuum and molecular emission lines (magnified by a factor ~4-5) in the source plane. The 12CO(9-8) and H2O emission lines have comparable spatial and kinematic distributions; the source-plane reconstructions do not clearly distinguish between a one-component and a two-component scenario, but the latter accounts for the observed broad line widths. HerBS-89a is a powerful star forming galaxy with a dust-to-gas ratio delta(GDR)~80, a SFR = 614 +/- 59 Msun/yr and a depletion timescale tau(depl) = (3.4 +/- 1.0) 1e8 years. The OH+ and CH+ absorption lines, all have their main velocity component red-shifted by ∆(V)~100 km/s relative to the global CO reservoir. We argue that these absorption lines trace a rare example of gas inflow towards the center of the galaxy.

We present Atacama Large Millimeter/Submillimeter Array (ALMA) continuum observations of a sample of nine star-forming galaxies at redshifts 1.47 and 2.23 selected from the High-$z$ Emission Line Survey (HiZELS). Four galaxies in our sample are detected at high significance by ALMA at a resolution of 0.25'' at rest-frame 355 $\mu$m. Together with the previously observed H$\alpha$ emission, from adaptive optics-assisted integral-field-unit spectroscopy (0.15'' resolution), and F606W and F140W imaging from the Hubble Space Telescope (0.2'' resolution), we study the star-formation activity, stellar and dust mass in these high-redshift galaxies at $\sim$kpc-scale resolution. We find that ALMA detection rates are higher for more massive galaxies ($M_*>10^{10.5}$ M$_\odot$) and higher [N \sc ii]/H$\alpha$ ratios ($>0.25$, a proxy for gas-phase metallicity). The dust extends out to a radius of 8 kpc, with a smooth structure, even for those galaxies presenting clumpy H$\alpha$ morphologies. The half-light radii ($R_{\rm dust}$) derived for the detected galaxies are of the order $\sim$4.5 kpc, more than twice the size of submillimetre-selected galaxies at a similar redshift. Our global star-formation rate estimates -- from far-IR and extinction-corrected H$\alpha$ luminosities -- are in good agreement. However, the different morphologies of the different phases of the interstellar medium suggest complex extinction properties of the high-redshift normal galaxies.

We present new observations with the Atacama Large Millimeter/sub-millimeter Array of the 122um and 205um fine-structure line emission of singly-ionised nitrogen in a strongly lensed starburst galaxy at z=2.6. The 122/205um [NII] line ratio is sensitive to electron density, n_e, in the ionised interstellar medium, and we use this to measure n_e~300cm^-3 averaged across the galaxy. This is over an order of magnitude higher than the Milky Way average, but comparable to localised Galactic star-forming regions. Combined with observations of the atomic carbon (CI(1-0)) and carbon monoxide (CO(4-3)) in the same system, we reveal the conditions in this intensely star-forming system. The majority of the molecular interstellar medium has been driven to high density, and the resultant conflagration of star formation produces a correspondingly dense ionised phase, presumably co-located with myriad HII regions that litter the gas-rich disk.

We report new observations toward the hyper-luminous dusty starbursting major merger ADFS-27 (z=5.655), using ATCA and ALMA. We detect CO 2-1, 8-7, 9-8, 10-9 and H2O(321-221) emission, and a P-Cygni-shaped OH+(11-01) absorption/emission feature. We also tentatively detect H2O(321-312) and OH+(12-01) emission and CH+(1-0) absorption. We find a total cold molecular mass of M_gas = (2.1+/-0.2) x 10^11 (alpha_CO/1.0) Msun. We also find that the excitation of the star-forming gas is overall moderate for a z>5 dusty starburst, which is consistent with its moderate dust temperature. A high density, high kinetic temperature gas component embedded in the gas reservoir is required to fully explain the CO line ladder. This component is likely associated with the "maximum starburst" nuclei in the two merging galaxies, which are separated by only (140+/-13) km/s along the line of sight and 9.0 kpc in projection. The kinematic structure of both components is consistent with galaxy disks, but this interpretation remains limited by the spatial resolution of the current data. The OH+ features are only detected towards the northern component, which is also the one that is more enshrouded in dust and thus remains undetected up to 1.6 um even in our sensitive new HST/WFC3 imaging. The absorption component of the OH+ line is blueshifted and peaks near the CO and continuum emission peak while the emission is redshifted and peaks offset by 1.7 kpc from the CO and continuum emission peak, suggesting that the gas is associated with a massive molecular outflow from the intensely star-forming nucleus that supplies 125 Msun/yr of enriched gas to its halo.

We analyse a robust sample of 30 near-infrared-faint (K>25.3, 5 sigma) submillimetre galaxies selected across a 0.96 deg^2 field, to investigate their properties and the cause of their lack of detectable optical/near-infrared emission. Our analysis exploits precise identifications based on ALMA 870um continuum imaging, combined with the very deep near-infrared imaging from the UKIDSS-UDS survey. We estimate that K>25.3 submillimetre galaxies represent 15+/-2 per cent of the total population brighter than S870=3.6mJy, with an expected surface density of ~450/deg^2 above S870>1mJy. As such they pose a source of contamination in surveys for both high-redshift "quiescent" galaxies and very-high-redshift Lyman-break galaxies. We show that these K-faint submillimetre galaxies are simply the tail of the broader submillimetre population, with comparable dust and stellar masses to K<25.3 mag submillimetre galaxies, but lying at significantly higher redshifts (z=3.44+/-0.06 versus z=2.36+/-0.11) and having higher dust attenuation (Av=5.2+/-0.3 versus Av=2.9+/-0.1). We investigate the origin of the strong dust attenuation and find indications that these K-faint galaxies have smaller dust continuum sizes than the K<25.3 galaxies, as measured by ALMA, which suggests their high attenuation is related to their compact sizes. We find a correlation of dust attenuation with star-formation rate surface density (Sigma_SFR), with the K-faint submillimetre galaxies representing the higher-Sigma_SFR and highest-Av galaxies. The concentrated, intense star-formation activity in these systems is likely to be associated with the formation of spheroids in compact galaxies at high redshifts, but as a result of their high obscuration these are completely missed in UV, optical and even near-infrared surveys.

We make use of sensitive (9.3 microJy/beam RMS) 1.2mm-continuum observations from the ASPECS ALMA large program of the Hubble Ultra Deep Field (HUDF) to probe dust-enshrouded star formation from 1362 Lyman-break galaxies spanning the redshift range z=1.5-10 (to ~7-28 Msolar/yr at 4 sigma over the entire range). We find that the fraction of ALMA-detected galaxies in our z=1.5-10 samples increases steeply with stellar mass, with the detection fraction rising from 0% at 10^9 Msolar to 85(-18)(+9)% at >10^10 Msolar. Moreover, stacking all 1253 low-mass (<10^9.25 Msolar) galaxies over the ASPECS footprint, we find a mean continuum flux of -0.1+/-0.4 microJy/beam, implying a hard upper limit on the obscured SFR of <0.6 Msolar/yr (4 sigma) in a typical low-mass galaxy. The correlation between the infrared excess IRX of UV-selected galaxies (L(IR)/L(UV)) and the UV-continuum slope is also seen in our ASPECS data and shows consistency with a Calzetti-like relation at >10^9.5 M_solar and a SMC-like relation at lower masses. Using stellar-mass and beta measurements for z~2 galaxies over CANDELS, we derive a new empirical relation between beta and stellar mass and then use this correlation to show that our IRX-beta and IRX-stellar mass relations are consistent with each other. We then use these constraints to express the infrared excess as a bivariate function of beta and stellar mass. Finally, we present updated estimates of star-formation rate density determinations at z>3, leveraging current improvements in the measured infrared excess and recent probes of ultra-luminous far-IR galaxies at z>2.

We combine the recent determination of the evolution of the cosmic density of molecular gas (H_2) using deep, volumetric surveys, with previous estimates of the cosmic density of stellar mass, star formation rate and atomic gas (HI), to constrain the evolution of baryons associated with galaxies averaged over cosmic time and space. The cosmic HI and H_2 densities are roughly equal at z~1.5. The H_2 density then decreases by a factor 6^+3_-2 to today's value, whereas the HI density stays approximately constant. The stellar mass density is increasing continuously with time and surpasses that of the total gas density (HI and H_2) at redshift z~1.5. The growth in stellar mass cannot be accounted for by the decrease in cosmic H_2 density, necessitating significant accretion of additional gas onto galaxies. With the new H_2 constraints, we postulate and put observational constraints on a two step gas accretion process: (i) a net infall of ionized gas from the intergalactic/circumgalactic medium to refuel the extended HI reservoirs, and (ii) a net inflow of HI and subsequent conversion to H_2 in the galaxy centers. Both the infall and inflow rate densities have decreased by almost an order of magnitude since z~2. Assuming that the current trends continue, the cosmic molecular gas density will further decrease by about a factor of two over the next 5 Gyr, the stellar mass will increase by approximately 10%, and cosmic star formation activity will decline steadily toward zero, as the gas infall and accretion shut down.

We present a survey of the molecular gas in 61 submillimetre galaxies (SMGs) selected from 870$\mu$m continuum surveys of the COSMOS, UDS and ECDFS fields, using the Atacama Large Millimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA). 46 $^{12}$CO ($J=$2-5) emission lines are detected in 45 of the targets at $z=$1.2-4.8, with redshifts indicating that those which are submillimetre bright and undetected/faint in the optical/near-infrared typically lie at higher redshifts, with a gradient of $\Delta z/\Delta S_{870}=$0.11$\pm$0.04mJy$^{-1}$. We also supplement our data with literature sources to construct a statistical CO spectral line energy distribution and find the $^{12}$CO line luminosities in SMGs peak at $J_{\rm up}\sim$6, consistent with the Cosmic Eyelash, among similar studies. Our SMGs lie mostly on or just above the main sequence, displaying a decrease in their gas depletion timescales $t_{\rm dep} = M_{\rm gas}/{\rm SFR}$ with redshift in the range $z\sim$1-5 and a median of 200$\pm$50Myr at $z\sim$2.8. This coincides with an increase in molecular gas fraction $\mu_{\rm gas} = M_{\rm gas}/M_\ast$ across the same redshift range. Finally we demonstrate that the $M_{\rm baryon}$-$\sigma$ distribution of our SMGs is consistent with that followed by early-type galaxies in the Coma cluster, providing strong support to the suggestion that SMGs are progenitors of massive local spheroidal galaxies. On the basis of this we suggest that the SMG populations above and below an 870-$\mu$m flux limit of $S_{870}\sim$5mJy may correspond to the division between slow- and fast-rotators seen in local early-type galaxies.

We resolve the host galaxies of seven gravitationally lensed quasars at redshift 1.5 to 2.8 using observations with the Atacama Large (sub-)Millimetre Array. Using a visibility-plane lens modelling technique, we create pixellated reconstructions of the dust morphology, and CO line morphology and kinematics. We find that the quasar hosts in our sample can be distinguished into two types: 1) galaxies characterised by clumpy, extended dust distributions ($R_{\rm eff}\sim2$ kpc) and mean star formation rate surface densities comparable to sub-mm-selected dusty star-forming galaxies ($\Sigma_{\rm SFR}\sim3$ M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$); 2) galaxies that have sizes in dust emission similar to coeval passive galaxies and compact starbursts ($R_{\rm eff}\sim0.5$ kpc), with high mean star formation rate surface densities ($\Sigma_{\rm SFR}=$ 400$-$4500 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$) that may be Eddington-limited or super-Eddington. The small size of some quasar hosts suggests that we observe them at a stage in their transformation into compact spheroids, where a high density of dynamically unstable gas leads to efficient star formation and black hole accretion. For the one system where we probe the mass of the gas reservoir, we find a gas fraction of just $0.06 \pm 0.04$ and a depletion timescale of $50 \pm 40$ Myr, suggesting it is transitioning into quiescence. In general, we expect that the extreme level of star formation in the compact quasar host galaxies will rapidly exhaust their gas reservoirs and could quench with or without help from active galactic nuclei feedback.

Galaxy clusters in the local universe descend from high-redshift overdense regions known as protoclusters. The large gas reservoirs and high rate of galaxy interaction in protoclusters are expected to trigger star-formation activity and luminous SMBH accretion in the host galaxies. We investigated the AGN content of a gas-rich and starbursting protocluster at $z=4$, known as the Distant Red Core (DRC). We observed with Chandra (139 ks) the 13 identified members of the structure, and searched for luminous and possibly obscured AGN among them. We also tested whether a hidden AGN can power the $Ly\alpha$ blob (LAB) detected with VLT/MUSE in the DRC. We detected obscured X-ray emission from the two most gas-rich members of the DRC, named DRC-1 and DRC-2. Both of them are resolved into multiple interacting clumps in high-resolution ALMA and HST observations. In particular, DRC-2 is found to host a luminous ($L_{2-10\,\mathrm{keV}}\approx3\times10^{45}\,\mathrm{erg\,s^{-1}}$) Compton-thick ($N_H\gtrsim10^{24}\,\mathrm{cm^{-2}}$) QSO, comparable to the most luminous QSOs known at all cosmic times. The AGN fraction among DRC members is consistent with results found for lower redshift protoclusters. However, X-ray stacking analysis reveals that SMBH accretion is likely also taking place in other DRC galaxies that are not detected individually by Chandra. Our results point toward the presence of a strong link between large gas reservoirs, galaxy interactions, and luminous and obscured nuclear activity in protocluster members. The powerful and obscured QSO detected in DRC-2 is likely powering the nearby LAB detected with VLT/MUSE, possibly through photoionization; however, we propose that the diffuse $Ly\alpha$ emission may be due to gas shocked by a massive outflow launched by DRC-2 over a $\approx10$ kpc scale.

We present Atacama Large Millimeter/sub-millimeter Array (ALMA) observations towards 27 low-redshift ($0.02< z<0.2$) star-forming galaxies taken from the Valparaíso ALMA/APEX Line Emission Survey (VALES). We perform stacking analyses of the $^{12}$CO($1-0$), $^{13}$CO($1-0$) and C$^{18}$O($1-0$) emission lines to explore the $L'$ ($^{12}$CO($1-0$))/$L'$($^{13}$CO($1-0$))) (hereafter $L'$($^{12}$CO)/$L'$($^{13}$CO)) and $L'$($^{13}$CO($1-0$))/$L'$(C$^{18}$O($1-0$)) (hereafter $L'$($^{13}$CO)/$L'$(C$^{18}$O) line luminosity ratio dependence as a function of different global galaxy parameters related to the star formation activity. The sample has far-IR luminosities $10^{10.1-11.9}$L$_{\odot}$ and stellar masses of $10^{9.8-10.9}$M$_{\odot}$ corresponding to typical star-forming and starburst galaxies at these redshifts. On average we find a $L'$($^{12}$CO)/$L'$($^{13}$CO) line luminosity ratio value of 16.1$\pm$2.5. Galaxies with evidences of possible merging activity tend to show higher $L'$($^{12}$CO)/$L'$($^{13}$CO) ratios by a factor of two, while variations of this order are also found in galaxy samples with higher star formation rates or star formation efficiencies. We also find an average $L'$($^{13}$CO)/$L'$(C$^{18}$O) line luminosity ratio of 2.5$\pm$0.6, which is in good agreement with those previously reported for starburst galaxies. We find that galaxy samples with high $L_{\text{IR}}$, SFR and SFE show low $L'$($^{13}$CO)/$L'$(C$^{18}$O) line luminosity ratios with high $L'$($^{12}$CO)/$L'$($^{13}$CO) line luminosity ratios, suggesting that these trends are produced by selective enrichment of massive stars in young starbursts.

We report observations with the Atacama Large Millimetre Array (ALMA) of six submillimetre galaxies (SMGs) within 3 arcmin of the Distant Red Core (DRC) at $z=4.0$, a site of intense cluster-scale star formation, first reported by Oteo et al. (2018). We find new members of DRC in three SMG fields; in two fields, the SMGs are shown to lie along the line of sight towards DRC; one SMG is spurious. Although at first sight this rate of association is consistent with earlier predictions, associations with the bright SMGs are rarer than expected, which suggests caution when interpreting continuum over-densities. We consider the implications of all 14 confirmed DRC components passing simultaneously through an active phase of star formation. In the simplest explanation, we see only the tip of the iceberg in terms of star formation and gas available for future star formation, consistent with our remarkable finding that the majority of newly confirmed DRC galaxies are not the brightest continuum emitters in their immediate vicinity. Thus while ALMA continuum follow-up of SMGs identifies the brightest continuum emitters in each field, it does not necessarily reveal all the gas-rich galaxies. To hunt effectively for protocluster members requires wide and deep spectral-line imaging to uncover any relatively continuum-faint galaxies that are rich in atomic or molecular gas. Searching with short-baseline arrays or single-dish facilities, the true scale of the underlying gas reservoirs may be revealed.

Using the EMIR instrument on the IRAM 30m telescope, we conducted a spectroscopic redshift search of seven z$_{\rm phot}$ $\sim$ 4 sub-millimetre bright galaxies selected from the Herschel Bright Sources (HerBS) sample with fluxes at 500 $\mu$m greater than 80 mJy. For four sources, we obtained spectroscopic redshifts between 3.4 < z < 4.1 through the detection of multiple CO-spectral lines with J $\leq$ 3. Later, we detected low-J transitions for two of these sources with the GBT including the CO(1-0) transition. For the remaining three sources, more data are needed to determine the spectroscopic redshift unambiguously. The measured CO luminosities and line widths suggest that all these sources are gravitationally lensed. These observations demonstrate that the 2 mm window is indispensable to confirm robust spectroscopic redshifts for z < 4 sources. Finally, we present an efficient graphical method to correctly identify spectroscopic redshifts.

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

We make use of the ASAGAO, deep 1.2 mm continuum observations of a 26 arcmin$^2$ region in the GOODS-South field obtained with ALMA, to probe dust-enshrouded star formation in $K$-band selected (i.e., stellar mass selected) galaxies, which are drawn from the ZFOURGE catalog. Based on the ASAGAO combined map, which was created by combining ASAGAO and ALMA archival data in the GOODS-South field, we find that 24 ZFOURGE sources have 1.2 mm counterparts with a signal-to-noise ratio $>$ 4.5 (1$\sigma\simeq$ 30 - 70 $\mu$Jy beam$^{-1}$ at 1.2 mm). Their median redshift is estimated to be $z_\mathrm{median}=$ 2.38 $\pm$ 0.14. They generally follow the tight relationship of the stellar mass versus star formation rate (i.e., the main sequence of star-forming galaxies). ALMA-detected ZFOURGE sources exhibit systematically larger infrared (IR) excess (IRX $\equiv L_\mathrm{IR}/L_\mathrm{UV}$) compared to ZFOURGE galaxies without ALMA detections even though they have similar redshifts, stellar masses, and star formation rates. This implies the consensus stellar-mass versus IRX relation, which is known to be tight among rest-frame-UV-selected galaxies, can not fully predict the ALMA detectability of stellar-mass-selected galaxies. We find that ALMA-detected ZFOURGE sources are the main contributors to the cosmic IR star formation rate density at $z$ = 2 - 3.

Using the NSF's Karl G. Jansky Very Large Array (VLA), we report six detections of CO(J=1-0) emission and one upper limit in z=2-3 galaxies originally detected in higher-J CO emission in the Atacama Large submillimeter/Millimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). From the CO(J=1-0) line strengths, we measure total cold molecular gas masses of M_gas = 2.4-11.6 x 10^10 (alpha_CO/3.6) Msun. We also measure a median CO(J=3-2) to CO(J=1-0) line brightness temperature ratio of r_31 = 0.84 +/- 0.26, and a CO(J=7-6) to CO(J=1-0) ratio range of r_71 <0.05 to 0.17. These results suggest that CO(J=3-2) selected galaxies may have a higher CO line excitation on average than CO(J=1-0) selected galaxies, based on the limited, currently available samples from the ASPECS and VLA CO Luminosity Density at High Redshift (COLDz) surveys. This implies that previous estimates of the cosmic density of cold gas in galaxies based on CO(J=3-2) measurements should be revised down by a factor of ~=2 on average based on assumptions regarding CO excitation alone. This correction further improves the agreement between the best currently existing constraints on the cold gas density evolution across cosmic history from line scan surveys, and the implied characteristic gas depletion times.

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5\mu$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

Measurements of the cosmic far-infrared background (CIB) indicate that emission from many extragalactic phenomena, including star formation and black hole accretion, in the Universe can be obscured by dust. Resolving the CIB to study the population of galaxies in which this activity takes place is a major goal of submillimetre astronomy. Here, we present interferometric 650$\mu$m submillimetre number counts. Using the Band 8 data from the ALMACAL survey, we have analysed 81 ALMA calibrator fields together covering a total area of 5.5~arcmin$^2$. The typical central rms in these fields is $\sim 100 \mu$Jy~beam$^{-1}$ with the deepest maps reaching $\sigma = 47 \mu$Jy~beam$^{-1}$ at sub-arcsec resolution. Multi-wavelength coverage from ALMACAL allows us to exclude contamination from jets associated with the calibrators. However, residual contamination by jets and lensing remain a possibility. Using a signal-to-noise threshold of $4.5\sigma$, we find 21 dusty, star-forming galaxies with 650$\mu$m flux densities of $\geq 0.7 $mJy. At the detection limit we resolve $\simeq 100$ per cent of the CIB at 650$\mu$m, a significant improvement compared to low resolution studies at similar wavelength. We have therefore identified all the sources contributing to the EBL at 650 microns and predict that the contribution from objects with flux 0.7<mJy will be small.

Historically, GHz radio emission has been used extensively to characterize the star-formation activity in galaxies. In this work, we look for empirical relations amongst the radio luminosity, the infrared luminosity, and the CO-based molecular gas mass. We assemble a sample of 278 nearby galaxies with measurements of radio continuum and total infrared emission, and the $^{12}$CO (J = 1-0) emission line. We find a correlation between the radio continuum and the CO emission line (with a scatter of 0.36 dex), in a large sample of different kind of galaxies. Making use of this correlation, we explore the evolution of the molecular gas mass function and the cosmological molecular gas mass density in six redshift bins up to $z = 1.5$. These results agree with previous semi-analytic predictions and direct measurements: the cosmic molecular gas density increases up to $z=1.5$. In addition, we find a single plane across five orders of magnitude for the explored luminosities, with a scatter of 0.27 dex. These correlations are sufficiently robust to be used for samples where no CO measurements exist.

We report the detection of CO($J$=2$\to$1) emission from three massive dusty starburst galaxies at $z$$>$5 through molecular line scans in the NSF's Karl G. Jansky Very Large Array (VLA) CO Luminosity Density at High Redshift (COLDz) survey. Redshifts for two of the sources, HDF 850.1 ($z$=5.183) and AzTEC-3 ($z$=5.298), were previously known. We revise a previous redshift estimate for the third source GN10 ($z$=5.303), which we have independently confirmed through detections of CO $J$=1$\to$0, 5$\to$4, 6$\to$5, and [CII] 158 $\mu$m emission with the VLA and the NOrthern Extended Milllimeter Array (NOEMA). We find that two currently independently confirmed CO sources in COLDz are "optically dark", and that three of them are dust-obscured galaxies at $z$$>$5. Given our survey area of $\sim$60 arcmin$^2$, our results appear to imply a $\sim$6-55 times higher space density of such distant dusty systems within the first billion years after the Big Bang than previously thought. At least two of these $z$$>$5 galaxies show star-formation rate surface densities consistent with so-called "maximum" starbursts, but we find significant differences in CO excitation between them. This result may suggest that different fractions of the massive gas reservoirs are located in the dense, star-forming nuclear regions - consistent with the more extended sizes of the [CII] emission compared to the dust continuum and higher [CII]-to-far-infrared luminosity ratios in those galaxies with lower gas excitation. We thus find substantial variations in the conditions for star formation between $z$$>$5 dusty starbursts, which typically have dust temperatures $\sim$57%$\pm$25% warmer than starbursts at $z$=2-3 due to their enhanced star formation activity.

Recent simulations and observations of massive galaxy cluster evolution predict that the majority of stellar mass build up happens within cluster members by $z=2$, before cluster virialization. Protoclusters rich with dusty, star-forming galaxies (DSFGs) at $z>3$ are the favored candidate progenitors for these massive galaxy clusters at $z\sim0$. We present here the first study analyzing stellar emission along with cold dust and gas continuum emission in a spectroscopically confirmed $z=4.002$ protocluster core rich with DSFGs, the Distant Red Core (DRC). We combine new HST and Spitzer data with existing Gemini, Herschel, and ALMA observations to derive individual galaxy-level properties, and compare them to coeval field and other protocluster galaxies. All of the protocluster members are massive ($>10^{10}$ M$_\odot$), but not significantly more so than their coeval field counterparts. Within uncertainty, all are nearly indistinguishable from galaxies on the star-forming vs. stellar mass main-sequence relationship, and on the star formation efficiency plane. Assuming no future major influx of fresh gas, we estimate that these gaseous DSFGs will deplete their gas reservoirs in $\sim300$ Myr, becoming the massive quiescent ellipticals dominating cluster cores by $z\sim3$. Using various methodologies, we derive a total $z=4$ halo mass of $\sim10^{14}$ M$_\odot$, and estimate that the DRC will evolve to become an ultra-massive cluster core of mass $\gtrsim10^{15}$ M$_\odot$ by $z=0$.

With the spatial resolution of the Atacama Large Millimetre Array (ALMA), dusty galaxies in the distant Universe typically appear as single, compact blobs of dust emission, with a median half-light radius, $\approx$ 1 kpc. Occasionally, strong gravitational lensing by foreground galaxies or galaxy clusters has probed spatial scales 1-2 orders of magnitude smaller, often revealing late-stage mergers, sometimes with tantalising hints of sub-structure. One lensed galaxy in particular, the Cosmic Eyelash at $z=$ 2.3, has been cited extensively as an example of where the interstellar medium exhibits obvious, pronounced clumps, on a spatial scale of $\approx$ 100 pc. Seven orders of magnitude more luminous than giant molecular clouds in the local Universe, these features are presented as circumstantial evidence that the blue clumps observed in many $z\sim$ 2-3 galaxies are important sites of ongoing star formation, with significant masses of gas and stars. Here, we present data from ALMA which reveal that the dust continuum of the Cosmic Eyelash is in fact smooth and can be reproduced using two Sérsic profiles with effective radii, 1.2 and 4.4 kpc, with no evidence of significant star-forming clumps down to a spatial scale of $\approx$ 80 pc and a star-formation rate of $<$ 3 M$_\odot$ yr$^{-1}$.

We present an ALMA study of the ~180 brightest sources in the SCUBA-2 map of the COSMOS field from the S2COSMOS survey, as a pilot study for AS2COSMOS - a full survey of the ~1,000 sources in this field. In this pilot we have obtained 870-um continuum maps of an essentially complete sample of the brightest 182 sub-millimetre sources (S_850um=6.2mJy) in COSMOS. Our ALMA maps detect 260 sub-millimetre galaxies (SMGs) spanning a range in flux density of S_870um=0.7-19.2mJy. We detect more than one SMG counterpart in 34+/-2 per cent of sub-millimetre sources, increasing to 53+/-8 per cent for SCUBA-2 sources brighter than S_850um>12mJy. We estimate that approximately one-third of these SMG-SMG pairs are physically associated (with a higher rate for the brighter secondary SMGs, S_870um>3mJy), and illustrate this with the serendipitous detection of bright [CII] 157.74um line emission in two SMGs, AS2COS0001.1 & 0001.2 at z=4.63, associated with the highest significance single-dish source. Using our source catalogue we construct the interferometric 870um number counts at S_870um>6.2mJy. We use the extensive archival data of this field to construct the multiwavelength spectral energy distribution of each AS2COSMOS SMG, and subsequently model this emission with MAGPHYS to estimate their photometric redshifts. We find a median photometric redshift for the S_870um>6.2mJy AS2COSMOS sample of z=2.87+/-0.08, and clear evidence for an increase in the median redshift with 870-um flux density suggesting strong evolution in the bright-end of the 870um luminosity function.

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

We present a new semi-empirical model for the dust continuum number counts of galaxies at 1.1 millimeter and 850 \micron. Our approach couples an observationally motivated model for the stellar mass and SFR distribution of galaxies with empirical scaling relations to predict the dust continuum flux density of these galaxies. Without a need to tweak the IMF, the model reproduces the currently available observations of the 1.1 millimeter and 850 \micron number counts, including the observed flattening in the 1.1 millimeter number counts below 0.3 mJy \citepGonzalez2019numbercounts and the number counts in discrete bins of different galaxy properties. Predictions of our work include : (1) the galaxies that dominate the number counts at flux densities below 1 mJy (3 mJy) at 1.1 millimeter (850 $\mu$m) have redshifts between $z=1$ and $z=2$, stellar masses of $\sim 5\times10^{10}~\rm{M}_\odot$, and dust masses of $\sim 10^{8}~\rm{M}_\odot$; (2) the flattening in the observed 1.1 millimeter number counts corresponds to the knee of the 1.1 millimeter luminosity function. A similar flattening is predicted for the number counts at 850 $\mu$m; (3) the model reproduces the redshift distribution of current 1.1 millimeter detections; (4) to efficiently detect large numbers of galaxies through their dust continuum, future surveys should scan large areas once reaching a 1.1 millimeter flux density of 0.1 mJy rather than integrating to fainter fluxes. Our modeling framework also suggests that the amount of information on galaxy physics that can be extracted from the 1.1 millimeter and 850 $\mu$m number counts is almost exhausted.

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

We analyse 870um Atacama Large Millimetre Array (ALMA) dust continuum detections of 41 canonically-selected z~3 Lyman-break galaxies (LBGs), as well as 209 ALMA-undetected LBGs, in follow-up of SCUBA-2 mapping of the UKIDSS Ultra Deep Survey (UDS) field. We find that our ALMA-bright LBGs lie significantly off the locally calibrated IRX-beta relation and tend to have relatively bluer rest-frame UV slopes (as parametrised by beta), given their high values of the 'infrared excess' (IRX=L_IR/L_UV), relative to the average 'local' IRX-beta relation. We attribute this finding in part to the young ages of the underlying stellar populations but we find that the main reason behind the unusually blue UV slopes are the relatively shallow slopes of the corresponding dust attenuation curves. We show that, when stellar masses are being established via SED fitting, it is absolutely crucial to allow the attenuation curves to vary (rather than fixing it on Calzetti-like law), where we find that the inappropriate curves may underestimate the resulting stellar masses by a factor of ~2-3x on average. In addition, we find these LBGs to have relatively high specific star-formation rates (sSFRs), dominated by the dust component, as quantified via the fraction of obscured star formation ( f_obs = SFR_IR/SFR_(UV+IR)). We conclude that the ALMA-bright LBGs are, by selection, massive galaxies undergoing a burst of a star formation (large sSFRs, driven, for example, by secular or merger processes), with a likely geometrical disconnection of the dust and stars, responsible for producing shallow dust attenuation curves.

Using the IRAM NOrthern Extended Millimeter Array (NOEMA), we conducted a program to measure redshifts for 13 bright galaxies detected in the Herschel Astrophysical Large Area Survey (H-ATLAS) with $S_{500{\mu}\rm m}\ge$80 mJy. We report reliable spectroscopic redshifts for 12 individual sources, which are derived from scans of the 3 and 2 mm bands, covering up to 31 GHz in each band, and are based on the detection of at least two emission lines. The spectroscopic redshifts are in the range $2.08<z<4.05$ with a median value of $z=2.9\pm$0.6. The sources are unresolved or barely resolved on scales of 10 kpc. In one field, two galaxies with different redshifts were detected. In two cases the sources are found to be binary galaxies with projected distances of ~140 kpc. The linewidths of the sources are large, with a mean value for the full width at half maximum of 700$\pm$300 km/s and a median of 800 km/s. We analyse the nature of the sources with currently available ancillary data to determine if they are lensed or hyper-luminous $L_{\rm FIR} > 10^{13}\,L_\odot$ galaxies. We also present a reanalysis of the spectral energy distributions including the continuum flux densities measured at 3 and 2 mm to derive the overall properties of the sources. Future prospects based on these efficient measurements of redshifts of high-z galaxies using NOEMA are outlined, including a comprehensive survey of all the brightest Herschel galaxies.

We present sensitive 850$\mu$m imaging of the COSMOS field using 640hr of new and archival observations taken with SCUBA-2 at the East Asian Observatory's James Clerk Maxwell Telescope. The SCUBA-2 COSMOS survey (S2COSMOS) achieves a median noise level of $\sigma_{850\mu{\mathrm{m}}}$=1.2mJy/beam over an area of 1.6 sq. degree (MAIN; HST/ACS footprint), and $\sigma_{850\mu{\mathrm{m}}}$=1.7mJy/beam over an additional 1 sq. degree of supplementary (SUPP) coverage. We present a catalogue of 1020 and 127 sources detected at a significance level of >4$\sigma$ and >4.3$\sigma$ in the MAIN and SUPP regions, respectively, corresponding to a uniform 2% false-detection rate. We construct the single-dish 850$\mu$m number counts at $S_{850}$>2mJy and show that these S2COSMOS counts are in agreement with previous single-dish surveys, demonstrating that degree-scale fields are sufficient to overcome the effects of cosmic variance in the $S_{850}$=2-10mJy population. To investigate the properties of the galaxies identified by S2COSMOS sources we measure the surface density of near-infrared-selected galaxies around their positions and identify an average excess of 2.0$\pm$0.2 galaxies within a 13$''$ radius (~100kpc at $z$~2). The bulk of these galaxies represent near-infrared-selected SMGs and/or spatially-correlated sources and lie at a median photometric redshift of $z$=2.0$\pm$0.1. Finally, we perform a stacking analysis at sub-millimeter and far-infrared wavelengths of stellar-mass-selected galaxies ($M_{\star}$=10$^{10}$-10$^{12}{\rm M_{\odot}}$) from $z$=0-4, obtaining high-significance detections at 850um in all subsets (SNR=4-30), and investigate the relation between far-infrared luminosity, stellar mass, and the peak wavelength of the dust SED. The publication of this survey adds a new deep, uniform sub-millimeter layer to the wavelength coverage of this well-studied COSMOS field.

The cosmic infrared background (CIB) provides a fundamental observational constraint on the star-formation history of galaxies over cosmic history. We estimate the contribution to the CIB from catalogued galaxies in the COSMOS field by using a novel map fitting technique on the \textitHerschel SPIRE maps. Prior galaxy positions are obtained using detections over a large range in wavelengths in the $K_{\rm s}$--3\u2009GHz range. Our method simultaneously fits the galaxies, the system foreground, and the leakage of flux from galaxies located in masked areas and corrects for an "over-fitting" effect not previously accounted for in stacking methods. We explore the contribution to the CIB as a function of galaxy survey wavelength and depth. We find high contributions to the CIB with the deep $r$ ($m_{\rm AB} \le 26.5$), $K_{\rm s}$ ($m_{\rm AB} \le 24.0$) and 3.6\,$\mu$m ($m_{\rm AB} \le 25.5$) catalogues. We combine these three deep catalogues and find a total CIB contributions of 10.5 $\pm$ 1.6, 6.7 $\pm$ 1.5 and 3.1 $\pm$ 0.7\u2009nWm$^{-2}$sr$^{-1}$ at 250, 350 and 500\,$\mu$m, respectively. Our CIB estimates are consistent with recent phenomenological models, prior based SPIRE number counts and with (though more precise than) the diffuse total measured by FIRAS. Our results raise the interesting prospect that the CIB contribution at $\lambda \le 500\,\mu$m from known galaxies has converged. Future large-area surveys like those with the Large Synoptic Survey Telescope are therefore likely to resolve a substantial fraction of the population responsible for the CIB at 250\,$\mu$m $\leq \lambda \leq$ 500\,$\mu$m.

We analyse the physical properties of a large, homogeneously selected sample of ALMA-located sub-mm galaxies (SMGs). This survey, AS2UDS, identified 707 SMGs across the ~1 sq.deg. field, including ~17 per cent, which are undetected at $K$>~25.7 mag. We interpret their ultraviolet-to-radio data using MAGPHYS and determine a median redshift of z=2.61+-0.08 (1$\sigma$ range of z=1.8-3.4) with just ~6 per cent at z>4. Our survey provides a sample of massive dusty galaxies at z>~1, with median dust and stellar masses of $M_d$=(6.8+-0.3)x10$^{8}$M$_\odot$ (thus, gas masses of ~10$^{11}$M$_\odot$) and $M_\ast=$(1.26+-0.05)x10$^{11}$M$_\odot$. We find no evolution in dust temperature at a constant far-infrared luminosity across z~1.5-4. The gas mass function of our sample increases to z~2-3 and then declines at z>3. The space density and masses of SMGs suggest that almost all galaxies with $M_\ast$>~3x10$^{11}$M$_\odot$ have passed through an SMG-like phase. The redshift distribution is well fit by a model combining evolution of the gas fraction in halos with the growth of halo mass past a threshold of $M_h$~6x10$^{12}$M$_\odot$, thus SMGs may represent the highly efficient collapse of gas-rich massive halos. We show that SMGs are broadly consistent with simple homologous systems in the far-infrared, consistent with a centrally illuminated starburst. Our study provides strong support for an evolutionary link between the active, gas-rich SMG population at z>1 and the formation of massive, bulge-dominated galaxies across the history of the Universe.