An ALMA survey of the S2CLS UDS field: optically invisible submillimetre galaxies

We analyse a robust sample of 30 near-infrared-faint (KAB > 25.3, 5σ) submillimetre galaxies (SMGs) selected from a 0.96 deg2 field to investigate their properties and the cause of their faintness in optical/near-infrared wavebands. Our analysis exploits precise identifications based on Atacama Large Millimeter Array (ALMA) 870-μm continuum imaging, combined with very deep near-infrared imaging from the UKIDSS Ultra Deep Survey. We estimate that SMGs with KAB > 25.3 mag represent 15 ± 2 per cent of the total population brighter than S870 = 3.6 mJy, with a potential surface density of ∼450 deg−2 above S870 ≥ 1 mJy. 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 SMGs represent the tail of the broader submillimetre population, with comparable dust and stellar masses to KAB ≤ 25.3 mag SMGs, 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 KAB ≤ 25.3 mag galaxies, as measured by ALMA, which suggests their high attenuation is related to their compact sizes. We identify a correlation of dust attenuation with star formation rate surface density (SFR), with the K-faint SMGs representing the higher 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 galaxies are completely missed in ultraviolet, optical, and even near-infrared surveys

An ALMA Survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS Field: The Far-infrared/Radio Correlation for High-redshift Dusty Star-forming Galaxies

We study the radio properties of 706 submillimeter galaxies (SMGs) selected at 870 μm with the Atacama Large Millimeter Array from the SCUBA-2 Cosmology Legacy Survey map of the Ultra Deep Survey field. We detect 273 SMGs at >4σ in deep Karl G. Jansky Very Large Array 1.4 GHz observations, of which a subset of 45 SMGs are additionally detected in 610 MHz Giant Metre-Wave Radio Telescope imaging. We quantify the far-infrared/radio correlation (FIRRC) through parameter q IR, defined as the logarithmic ratio of the far-infrared and radio luminosity, and include the radio-undetected SMGs through a stacking analysis. We determine a median q IR = 2.20 ± 0.03 for the full sample, independent of redshift, which places these z ~ 2.5 dusty star-forming galaxies 0.44 ± 0.04 dex below the local correlation for both normal star-forming galaxies and local ultra-luminous infrared galaxies (ULIRGs). Both the lack of redshift evolution and the offset from the local correlation are likely the result of the different physical conditions in high-redshift starburst galaxies, compared to local star-forming sources. We explain the offset through a combination of strong magnetic fields (B gsim 0.2 mG), high interstellar medium (ISM) densities and additional radio emission generated by secondary cosmic rays. While local ULIRGs are likely to have similar magnetic field strengths, we find that their compactness, in combination with a higher ISM density compared to SMGs, naturally explains why local and high-redshift dusty star-forming galaxies follow a different FIRRC. Overall, our findings paint SMGs as a homogeneous population of galaxies, as illustrated by their tight and nonevolving far-infrared/radio correlation

An ALMA Survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS Field: Number Counts of Submillimeter Galaxies

We report the first results of AS2UDS, an 870 μm continuum survey with the Atacama Large Millimeter/Submillimeter Array (ALMA) of a total area of ~50 arcmin2 comprising a complete sample of 716 submillimeter sources drawn from the SCUBA-2 Cosmology Legacy Survey (S2CLS) map of the UKIDSS/UDS field. The S2CLS parent sample covers a 0.96 degree2 field at σ 850 = 0.90 ± 0.05 mJy beam−1. Our deep, high-resolution ALMA observations with σ 870 ~ 0.25 mJy and a 0farcs15–0farcs30 FWHM synthesized beam, provide precise locations for 695 submillimeter galaxies (SMGs) responsible for the submillimeter emission corresponding to 606 sources in the low-resolution, single-dish map. We measure the number counts of SMGs brighter than S 870 ≥ 4 mJy, free from the effects of blending and show that the normalization of the counts falls by 28% ± 2% in comparison with the SCUBA-2 published counts, but that the shape remains unchanged. We determine that ${44}_{-14}^{+16}$% of the brighter single-dish sources with S 850 ≥ 9 mJy consist of a blend of two or more ALMA-detectable SMGs brighter than S 870 ~ 1 mJy (corresponding to a galaxy with a total-infrared luminosity of L IR gsim 1012 L ⊙), in comparison with 28% ± 2% for the single-dish sources at S 850 ≥ 5 mJy. Using the 46 single-dish submillimeter sources that contain two or more ALMA-detected SMGs with photometric redshifts, we show that there is a significant statistical excess of pairs of SMGs with similar redshifts (<1% probability of occurring by chance), suggesting that at least 30% of these blends arise from physically associated pairs of SMGs

An ALMA survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS field: number counts of submillimeter galaxies

We report the first results of AS2UDS, an 870 μm continuum survey with the Atacama Large Millimeter/Submillimeter Array (ALMA) of a total area of ~50 arcmin2 comprising a complete sample of 716 submillimeter sources drawn from the SCUBA-2 Cosmology Legacy Survey (S2CLS) map of the UKIDSS/UDS field. The S2CLS parent sample covers a 0.96 degree2 field at σ 850 = 0.90 ± 0.05 mJy beam−1. Our deep, high-resolution ALMA observations with σ 870 ~ 0.25 mJy and a 0farcs15–0farcs30 FWHM synthesized beam, provide precise locations for 695 submillimeter galaxies (SMGs) responsible for the submillimeter emission corresponding to 606 sources in the low-resolution, single-dish map. We measure the number counts of SMGs brighter than S 870 ≥ 4 mJy, free from the effects of blending and show that the normalization of the counts falls by 28% ± 2% in comparison with the SCUBA-2 published counts, but that the shape remains unchanged. We determine that ${44}_{-14}^{+16}$% of the brighter single-dish sources with S 850 ≥ 9 mJy consist of a blend of two or more ALMA-detectable SMGs brighter than S 870 ~ 1 mJy (corresponding to a galaxy with a total-infrared luminosity of L IR gsim 1012 L ⊙), in comparison with 28% ± 2% for the single-dish sources at S 850 ≥ 5 mJy. Using the 46 single-dish submillimeter sources that contain two or more ALMA-detected SMGs with photometric redshifts, we show that there is a significant statistical excess of pairs of SMGs with similar redshifts (<1% probability of occurring by chance), suggesting that at least 30% of these blends arise from physically associated pairs of SMGs

An ALMA survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS field: High-resolution dust continuum morphologies and the link between sub-millimetre galaxies and spheroid formation

We present an analysis of the morphology and profiles of the dust continuum emission in 153 bright sub-millimetre galaxies (SMGs) detected with ALMA at signal to noise ratios of >8 in high-resolution 0′′.18 (∼1 kpc) 870 μm maps. We measure sizes, shapes and light profiles for the rest-frame far-infrared emission from these luminous star-forming systems and derive a median effective radius (Re) of 0′′.10 ± 0′′.04 for our sample with a median flux of S870=5.6 ± 0.2 mJy. We find that the apparent axial ratio (b/a) distribution of the SMGs peaks at b/a ∼ 0.63 ± 0.02 and is best described by triaxial morphologies, while their emission profiles are best fit by a Sérsic model with n ≃ 1.0 ± 0.1, similar to exponential discs. This combination of triaxiality and n ∼ 1 Sérsic index are characteristic of bars and we suggest that the bulk of the 870 μm dust continuum emission in the central ∼2 kpc of these galaxies arises from bar-like structures. As such we caution against using the orientation of shape of the bright dust continuum emission at ≂ resolution to assess either the orientation of any disc on the sky or tits inclination. By stacking our 870 μm maps we recover faint extended dust continuum emission on ∼4 kpc scales which contributes 13 ± 1% of the total 870 μm emission. The scale of this extended emission is similar to that seen for the molecular gas and rest-frame optical light in these systems, suggesting that it represents an extended dust and gas disc at radii larger than the more active bar component. Including this component in our estimated size of the sources we derive a typical effective radius of ≃ 0′′.15 ± 0′′.05 or 1.2 ± 0.4 kpc. Our results suggest that kpc-scale bars are ubiquitous features of high star-formation rate systems at z ≫ 1, while these systems also contain fainter and more extended gas and stellar envelopes. We suggest that these features, seen some 10–12 Gyrs ago, represent the formation phase of the earliest galactic-scale components: stellar bulges

An ALMA Survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS Field: Identifying Candidate z ∼ 4.5 [C II] Emitters

We report the results of a search for serendipitous [C ii] 157.74 μm emitters at z sime 4.4–4.7 using the Atacama Large Millimeter/submillimeter Array (ALMA). The search exploits the AS2UDS continuum survey, which covers ~50 arcmin2 of the sky toward 695 luminous (S 870 gsim 1 mJy) submillimeter galaxies (SMGs), selected from the SCUBA-2 Cosmology Legacy Survey 0.96 deg2 Ultra Deep Survey (UDS) field. We detect 10 candidate line emitters, with an expected false detection rate of 10%. All of these line emitters correspond to 870 μm continuum-detected sources in AS2UDS. The emission lines in two emitters appear to be high-J CO, but the remainder have multi-wavelength properties consistent with [C ii] from z sime 4.5 galaxies. Using our sample, we place a lower limit of $\gt 5\times {10}^{-6}\,{\mathrm{Mpc}}^{-3}$ on the space density of luminous (L IR sime 1013 ${L}_{\odot }$) SMGs at z = 4.40–4.66, suggesting $\geqslant 7$% of SMGs with ${S}_{870\mu {\rm{m}}}\gtrsim 1$ mJy lie at 4 < z < 5. From stacking the high-resolution (~0farcs15 full-width half maximum) ALMA 870 μm imaging, we show that the [C ii] line emission is more extended than the continuum dust emission, with an average effective radius for the [C ii] of ${r}_{{\rm{e}}}={1.7}_{-0.2}^{+0.1}$ kpc, compared to r e = 1.0 ± 0.1 kpc for the continuum (rest-frame 160 μm). By fitting the far-infrared photometry for these galaxies from 100 to 870 μm, we show that SMGs at z ~ 4.5 have a median dust temperature of T d = 55 ± 4 K. This is systematically warmer than 870 μm selected SMGs at z sime 2, which typically have temperatures around 35 K. These z sime 4.5 SMGs display a steeper trend in the luminosity-temperature plane than z ≤ 2 SMGs. We discuss the implications of this result in terms of the selection biases of high-redshift starbursts in far-infrared/submillimeter surveys

Radio Spectra and Sizes of Atacama Large Millimeter/submillimeter Array-identified Submillimeter Galaxies: Evidence of Age-related Spectral Curvature and Cosmic-Ray Diffusion?

We analyze the multifrequency radio spectral properties of 41 6 GHz-detected Atacama Large Millimeter/submillimeter Array (ALMA)-identified, submillimeter galaxies (SMGs), observed at 610 MHz, 1.4 GHz, and 6 GHz with the Giant Metrewave Radio Telescope and the Very Large Array. Combining high-resolution (~0farcs5) 6 GHz radio and ALMA 870 μm imaging (tracing rest frame ~20 GHz, and ~250 μm dust continuum), we study the far-infrared/radio correlation via the logarithmic flux ratio q IR, measuring $\langle {q}_{\mathrm{IR}}\rangle =2.20\pm 0.06$ for our sample. We show that the high-frequency radio sizes of SMGs are ~1.9 ± 0.4× (~2–3 kpc) larger than those of the cool dust emission, and find evidence for a subset of our sources being extended on ~10 kpc scales at 1.4 GHz. By combining radio flux densities measured at three frequencies, we can move beyond simple linear fits to the radio spectra of high-redshift star-forming galaxies, and search for spectral curvature, which has been observed in local starburst galaxies. At least a quarter (10/41) of our sample shows evidence of a spectral break, with a median $\langle {\alpha }_{610\,\mathrm{GHz}}^{1.4\,\mathrm{GHz}}\rangle =-0.60\pm 0.06$, but $\langle {\alpha }_{1.4\,\mathrm{GHz}}^{6\,\mathrm{GHz}}\rangle =-1.06\pm 0.04$—a high-frequency flux deficit relative to simple extrapolations from the low-frequency data. We explore this result within this subset of sources in the context of age-related synchrotron losses, showing that a combination of weak magnetic fields (B ~ 35 μG) and young ages (t SB ~ 40–80 Myr) for the central starburst can reproduce the observed spectral break. Assuming these represent evolved (but ongoing) starbursts, and we are observing these systems roughly halfway through their current episode of star formation, this implies starburst durations of lesssim100 Myr, in reasonable agreement with estimates derived via gas depletion timescales

A Machine-learning Method for Identifying Multiwavelength Counterparts of Submillimeter Galaxies: Training and Testing Using AS2UDS and ALESS

We describe the application of supervised machine-learning algorithms to identify the likely multiwavelength counterparts to submillimeter sources detected in panoramic, single-dish submillimeter surveys. As a training set, we employ a sample of 695 (S 870μm gsim 1 mJy) submillimeter galaxies (SMGs) with precise identifications from the ALMA follow-up of the SCUBA-2 Cosmology Legacy Survey's UKIDSS-UDS field (AS2UDS). We show that radio emission, near-/mid-infrared colors, photometric redshift, and absolute H-band magnitude are effective predictors that can distinguish SMGs from submillimeter-faint field galaxies. Our combined radio + machine-learning method is able to successfully recover ~85% of ALMA-identified SMGs that are detected in at least three bands from the ultraviolet to radio. We confirm the robustness of our method by dividing our training set into independent subsets and using these for training and testing, respectively, as well as applying our method to an independent sample of ~100 ALMA-identified SMGs from the ALMA/LABOCA ECDF-South Survey (ALESS). To further test our methodology, we stack the 870 μm ALMA maps at the positions of those K-band galaxies that are classified as SMG counterparts by the machine learning but do not have a >4.3σ ALMA detection. The median peak flux density of these galaxies is S 870μm = (0.61 ± 0.03) mJy, demonstrating that our method can recover faint and/or diffuse SMGs even when they are below the detection threshold of our ALMA observations. In future, we will apply this method to samples drawn from panoramic single-dish submillimeter surveys that currently lack interferometric follow-up observations to address science questions that can only be tackled with large statistical samples of SMGs