Search

Massive, intensely star-forming galaxies at high redshift require a supply of molecular gas from their gas reservoirs, replenished by infall from the surrounding circumgalactic medium, to sustain their immense star-formation rates. However, our knowledge of the extent and morphology of their cold-gas reservoirs is still in its infancy. We present the results of stacking 80 hours of JVLA observations of CO(1--0) emission -- which traces the cold molecular gas -- in 19 $z=2.0-4.5$ dusty, star-forming galaxies from the AS2VLA survey. The visibility-plane stack reveals extended emission with a half-light radius of $3.8\pm0.5$~kpc, 2--3$\times$ more extended than the dust-obscured star formation and $1.4\pm0.2\times$ more extended than the stellar emission. Similarly, stacking the [CI](1--0) observations for a subsample of our galaxies yields sizes consistent with CO(1--0). The CO(1--0) size is comparable to the [CII] halos detected around high-redshift star-forming galaxies.The bulk (up to 80\%) of molecular gas resides outside the star-forming region; only a small part of their molecular gas reservoir directly contributes to their current star formation. Photon-dissociation region modelling indicates that the extended CO(1--0) emission arises from clumpy, dense clouds rather than smooth, diffuse gas., Comment: Submitted to A&A. 9 pages, 5 figures

In the present-day universe, the most massive galaxies are ellipticals located in the cores of galaxy clusters, harboring the heaviest super-massive black holes (SMBHs). However the mechanisms that drive the early growth phase and subsequent transformation of these morphology and kinematics of galaxies remain elusive. Here we report (sub)kiloparsec scale observations of stars, gas, and dust in ADF22.A1, a bright dusty starburst galaxy at z=3.1, hosting a heavily obscured active galactic nucleus and residing in a proto-cluster core. ADF22.A1 is a giant spiral galaxy with the kinematics of a rotating disk with rotation velocity Vrot=530+/-10km/s and diameter larger than 30 kpc. The high specific stellar angular momentum of this system, j*=3400+/-600 kpc km/s, requires a mechanism to effectively spin-up ADF22.A1, indicating the importance of accretion from the cosmic web to supply both gas and angular momentum to galaxies in their early gas-rich starburst phase. In its inner region, gas flows along dust lanes in a bar connected with the bright dusty core and the estimated mass ratio of a bulge to SMBH matches the local relation, suggesting that bars are a key mechanism to shape the early co-evolution of these components. Comparison with cosmological simulations shows that ADF22.A1 will likely evolve into a massive elliptical at the present day, experiencing a significant reduction in angular momentum associated with subsequent galaxy mergers., Comment: 13 pages, 12 figures. Submitted

We present the first results of the JWST Emission Line Survey (JELS). Utilising the first NIRCam narrow-band imaging at 4.7$\mu$m, over 63 arcmin$^{2}$ in the PRIMER/COSMOS field, we identified 609 emission line galaxy candidates. From these, we robustly selected 35 H$\alpha$ star-forming galaxies at $z \sim 6.1$, with H$\alpha$ star-formation rates ($\rm{SFR_{H\alpha}}$) $\sim0.9-15\ \rm{M_{\odot} \ yr^{-1}}$. Combining our unique H$\alpha$ sample with the exquisite panchromatic data in the field, we explored their physical properties and star-formation histories, and compared these to a broad-band selected sample at $z\sim 6$ which offered vital new insights into the nature of high-redshift galaxies. UV-continuum slopes ($\beta$) were considerably redder for our H$\alpha$ sample ($\langle\beta\rangle\sim-1.92$) compared to the broad-band sample ($\langle\beta\rangle\sim-2.35$). This was not due to dust attenuation as our H$\alpha$ sample was relatively dust-poor (median $A_V=0.23$); instead, we argued the reddened slopes could be due to nebular continuum. We compared $\rm{SFR_{H\alpha}}$ and the UV-continuum-derived $\rm{SFR_{UV}}$ to SED-fitted measurements averaged over canonical timescales of 10 and 100 Myr ($\rm{SFR_{10}}$ and $\rm{SFR_{100}}$). We found an increase in recent SFR for our sample of H$\alpha$ emitters, particularly at lower stellar masses ($<10^9 \ \rm{M_{\odot}}$). We also found $\rm{SFR_{H\alpha}}$ strongly traced SFR averaged over 10 Myr timescales, whereas the UV-continuum over-predicts SFR on 100 Myr timescales at low stellar masses. These results point to our H$\alpha$ sample undergoing `bursty' star formation. Our F356W $z \sim 6$ sample showed a larger scatter in $\rm{SFR_{10}/SFR_{100}}$ across all stellar masses, which highlighted how narrow-band photometric selections of H$\alpha$ emitters are key to quantifying the burstiness of star-formation activity., Comment: Submitted to MRNAS; 27 pages, 19 figures. Accompanying data will be publicly released upon acceptance of this manuscript

We present the JWST Emission Line Survey (JELS), a Cycle 1 JWST imaging programme exploiting the wavelength coverage and sensitivity of NIRCam to extend narrow-band rest-optical emission line selection into the epoch of reionization (EoR) for the first time, and to enable unique studies of the resolved ionised gas morphology in individual galaxies across cosmic history. The primary JELS observations comprise $\sim4.7\mu$m narrow-band imaging over $\sim63$ arcmin$^{2}$ designed to enable selection of H$\alpha$ emitters at $z\sim6.1$, as well as the selection of a host of novel emission-line samples including [OIII] at $z\sim8.3$ and Pa $\alpha/\beta$ at $z\sim1.5/2.8$. For the prime F466N and F470N narrow-band observations, the emission-line sensitivities achieved are up to $\sim2\times$ more sensitive than current slitless spectroscopy surveys (5$\sigma$ limits of 1.1-1.6$\times10^{-18}\text{erg s}^{-1}\text{cm}^{-2}$), corresponding to unobscured H$\alpha$ star-formation rates (SFRs) of 1-1.6 $\text{M}_{\odot}\,\text{yr}^{-1}$ at $z\sim6.1$ and extending emission-line selections in the EoR to fainter populations. Simultaneously, JELS also obtained F200W broadband and F212N narrow-band imaging (H$\alpha$ at $z\sim2.23$) that probes SFRs $\gtrsim5\times$ fainter than previous ground-based narrow-band studies ($\sim0.2 \text{M}_{\odot}\text{yr}^{-1}$), offering an unprecedented resolved view of star formation at cosmic noon. In this paper we describe the detailed JELS survey design, key data processing steps specific to the survey observations, and demonstrate the exceptional data quality and imaging sensitivity achieved. We then summarise the key scientific goals of JELS and present some early science results, including examples of spectroscopically confirmed H$\alpha$ and [OIII] emitters discovered by JELS that illustrate the novel parameter space probed., Comment: Submitted to MRNAS; 13 pages, 8 figures. Accompanying data will be publicly released upon acceptance of this and companion paper (Pirie et al., arXiv:2410.11808). V2 updated to reflect revised author list and cross-citation to Pirie et al

We characterize, for the first time, the average extended emission in multiple lines ([OII], [OIII], and Hbeta) around a statistical sample of 560 galaxies at z~0.25-0.85. By stacking the Multi Unit Spectroscopic Explorer (MUSE) 3D data from two large surveys, the MUSE Analysis of Gas around Galaxies (MAGG) and the MUSE Ultra Deep Field (MUDF), we detect significant [OII] emission out to ~40 kpc, while [OIII] and Hbeta emission is detected out to ~30 kpc. Via comparisons with the nearby average stellar continuum emission, we find that the line emission at 20-30 kpc likely arises from the disk-halo interface. Combining our results with that of our previous study at z~1, we find that the average [OII] surface brightness increases independently with redshift over z~0.4-1.3 and with stellar mass over M* ~10^{6-12} Msun, which is likely driven by the star formation rate as well as the physical conditions of the gas. By comparing the observed line fluxes with photoionization models, we find that the ionization parameter declines with distance, going from log q (cm/s) ~7.7 at <=5 kpc to ~7.3 at 20-30 kpc, which reflects a weaker radiation field in the outer regions of galaxies. The gas-phase metallicity shows no significant variation over 30 kpc, with a metallicity gradient of ~0.003 dex/kpc, which indicates an efficient mixing of metals on these scales. Alternatively, there could be a significant contribution from shocks and diffuse ionized gas to the line emission in the outer regions., Comment: 9 pages, 6 figures, A&A in press

We present the first detailed radio study of the galaxy cluster XMMXCS J2215.9-1738 at $z$ = 1.46 using MeerKAT $L$-band (1.3 GHz) observations. We combine our radio observation with archival optical and infrared data to investigate the star formation and AGN population within $R_{200}$ ($R = $ 0.8 Mpc) of the cluster centre. Using three selection criteria; the radio luminosity, the far-infrared radio ratio ($q_{\rm{IR}}$) and the mid-infrared colour, we distinguish galaxies with radio emission predominantly powered by star formation from that powered by AGNs. We selected 24 cluster members within $R_{\rm{200}}$ in the MeerKAT image based on either their photometric or spectroscopic redshift. We classified 12/24 ($50\%$) as galaxies whose radio emission is dominated by star-formation activity, 6/24 ($25\%$) as intermediate star-forming galaxies and 6/24 ($25\%$) as AGN-dominated galaxies. Using the radio continuum luminosities of the star-forming cluster galaxies, we estimated an integrated star formation rate (SFR) value of 1700 $\pm$ 330 M$_{\odot}$yr$^{-1}$ within $R_{200}$. We derived a mass-normalized integrated SFR value of $(570 \pm 110) \times 10^{-14}$ yr$^{-1}$. This supports previous observational and theoretical studies that indicated a rapid increase in star formation activity within the core of high-redshift clusters. We also show that the high AGN fraction within the cluster core is consistent with previous cluster observations at $z >$ 1.5.

We present an analysis of the JWST NIRCam and MIRI morphological properties of 80 massive ($\log_{10}(M_\ast[M_{\odot}])$=11.2$\pm$0.1) dusty star-forming galaxies at $z$$=$2.7$^{+1.2}_{-0.7}$, identified as sub-millimetre galaxies (SMGs) by ALMA, that have been observed as part of the JWST PRIMER project. To compare the structure of these massive, active galaxies to more typical less actively star-forming galaxies, we define two comparison samples. The first of 850 field galaxies matched in specific star-formation rate and redshift and the second of 80 field galaxies matched in stellar mass. We identify 20$\pm$5% of the SMGs as candidate late-stage major mergers, a further 40$\pm$10% as potential minor mergers and 40$\pm$10% which have comparatively undisturbed disk-like morphologies, with no obvious massive neighbours. These rates are comparable to those for the field samples and indicate that the majority of the sub-millimetre-detected galaxies are not late-stage major mergers, but have interaction rates similar to the less-active population at $z$$\sim$2-3. We establish that SMGs have comparable near-infrared sizes to the less active populations, but exhibit lower S\'ersic indices, consistent with bulge-less disks and have more structured morphologies at 2$\mu$m relative to 4$\mu$m. We find evidence for dust reddening as the origin of the morphological differences between the populations, identifying a strong correlation between the F200W$-$F444W pixel colour and the 870$\mu$m surface brightness. We conclude that SMGs and less active galaxies at the same epochs share a common disk-like structure, but the weaker bulge components of the SMGs results in a lower dynamical stability. Consequently, instabilities triggered either secularly or by minor external perturbations result in higher levels of activity (and dust content) in SMGs compared to typical star-forming galaxies. [Abridged], Comment: 23 pages, 13 figures. Accepted for publication in A&A

The CO(1--0) and [\ion{C}{1}](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 [\ion{C}{1}](1--0) emission in a sample of 20 unlensed dusty, star-forming galaxies at $z=2-5$. The CO(1--0)/[\ion{C}{1}](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), [\ion{C}{1}](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 [\ion{C}{1}] 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 [\ion{C}{1}] 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.

We present a new method for modelling the kinematics of galaxies from interferometric observations by performing the optimization of the kinematic model parameters directly in visibility-space instead of the conventional approach of fitting velocity fields produced with the CLEAN algorithm in real-space. We demonstrate our method on ALMA observations of $^{12}$CO (2$-$1), (3$-$2) or (4$-$3) emission lines from an initial sample of 30 massive 850$\mu$m-selected dusty star-forming galaxies with far-infrared luminosities $\gtrsim$$\,10^{12}\,$L$_{\odot}$ in the redshift range $z \sim\,$1.2$-$4.7. Using the results from our modelling analysis for the 12 sources with the highest signal-to-noise emission lines and disk-like kinematics, we conclude the following: (i) Our sample prefers a CO-to-$H_2$ conversion factor, of $\alpha_{\rm CO} = 0.92 \pm 0.36$; (ii) These far-infrared luminous galaxies follow a similar Tully$-$Fisher relation between the circularized velocity, $V_{\rm circ}$, and baryonic mass, $M_{\rm b}$, as more typical star-forming samples at high redshift, but extend this relation to much higher masses $-$ showing that these are some of the most massive disk-like galaxies in the Universe; (iii) Finally, we demonstrate support for an evolutionary link between massive high-redshift dusty star-forming galaxies and the formation of local early-type galaxies using the both the distributions of the baryonic and kinematic masses of these two populations on the $M_{\rm b}\,-\,\sigma$ plane and their relative space densities.

We present a multi-wavelength analysis using SMA, JCMT, NOEMA, JWST, HST, and SST of two dusty strongly star-forming galaxies, 850.1 and 850.2, seen through the massive cluster lens A1489. These SMA-located sources both lie at z=4.26 and have bright dust continuum emission, but 850.2 is a UV-detected Lyman-break galaxy, while 850.1 is undetected at <2um, even with deep JWST/NIRCam observations. We investigate their stellar, ISM, and dynamical properties, including a pixel-level SED analysis to derive sub-kpc-resolution stellar-mass and Av maps. We find that 850.1 is one of the most massive and highly obscured, Av~5, galaxies known at z>4 with M*~10^11.8 Mo (likely forming at z>6), and 850.2 is one of the least massive and least obscured, Av~1, members of the z>4 dusty star-forming population. The diversity of these two dust-mass-selected galaxies illustrates the incompleteness of galaxy surveys at z>3-4 based on imaging at <2um, the longest wavelengths feasible from HST or the ground. The resolved mass map of 850.1 shows a compact stellar mass distribution, Re(mass)~1kpc, but its expected evolution to z~1.5 and then z~0 matches both the properties of massive, quiescent galaxies at z~1.5 and ultra-massive early-type galaxies at z~0. We suggest that 850.1 is the central galaxy of a group in which 850.2 is a satellite that will likely merge in the near future. The stellar morphology of 850.1 shows arms and a linear bar feature which we link to the active dynamical environment it resides within., Comment: Submitted to ApJ, comments welcome!

We present spatially resolved kinematics of 27 ALMA-identified dust-obscured star-forming galaxies (DSFGs) at $z\sim$1.3-2.6, as traced by H$\alpha$ emission using VLT/KMOS near-infrared integral field spectroscopy from the "KMOS-ALMA Observations of Submillimetre Sources" (KAOSS) Large Programme. We derive H$\alpha$ rotation curves and velocity dispersion profiles for the DSFGs, and find that among the 27 sources with bright, spatially extended H$\alpha$ emission, 24 display evidence for disc-like kinematics. We measure a median inclination-corrected velocity at 2.2$R_{\rm d}$ of $v_{\rm rot}=$190$\pm$40kms$^{-1}$ and intrinsic velocity dispersion of $\sigma_0=$87$\pm$6kms$^{-1}$ for these disc-like sources. The kinematics yield median circular velocities of $v_{\rm circ}=$230$\pm$20kms$^{-1}$ and dynamical masses within 2$R_{\rm e}$ ($\sim$7kpc radius) of $M_{\rm dyn}=$(1.1$\pm$0.2)$\times$10$^{11}$M$_\odot$. Compared to less actively star-forming galaxies, KAOSS DSFGs are both faster rotating with higher intrinsic velocity dispersions, but have similar $v_{\rm rot}/\sigma_0$ ratios, median $v/\sigma_0=$2.5$\pm$0.5. We suggest that the kinematics of the DSFGs are primarily rotation supported but with a non-negligible contribution from pressure support, which may be driven by star formation or mergers/interactions. We estimate the normalisation of the stellar mass Tully-Fisher relation (sTFR) for the disc-like DSFGs and compare it with local studies, finding no evolution at fixed slope between $z\sim$2 and $z\sim$0. Finally, we show that the kinematic properties of the DSFG population are consistent with them evolving into massive early-type galaxies, the dominant $z\sim$0 population at these masses., Comment: 24 pages, 15 figures, accepted for publication in MNRAS, updated author list

The James Webb Space Telescope Time-Domain Field (JWST-TDF) is an $\sim$14$'$ diameter field near the North Ecliptic Pole that will be targeted by one of the JWST Guaranteed Time Observations programs. Here, we describe our James Clerk Maxwell Telescope SCUBA-2 850 $\mu$m imaging of the JWST-TDF and present the submillimeter source catalog and properties. We also present a catalog of radio sources from Karl J. Jansky Very Large Array 3 GHz observations of the field. These observations were obtained to aid JWST's study of the dust-obscured galaxies that contribute significantly to the cosmic star formation at high redshifts. Our deep 850 $\mu$m map covers the JWST TDF at a noise level of $\sigma_{850}$ = 1.0 mJy beam$^{-1}$, detecting 83/31 sources in the main/supplementary signal-to-noise ratio (S/N $>$ 4 / S/N = 3.5 - 4) sample respectively. The 3 GHz observations cover a 24$'$ diameter field with a 1 $\sigma$ noise of 1$\mu$Jy beam$^{-1}$ at a 0$.\!\!^{\prime\prime}$7 FWHM. We identified eighty-five 3 GHz counterparts to sixty-six 850 $\mu$m sources and then matched these with multiwavelength data from the optical to the mid-infrared wave bands. We performed spectral energy distribution fitting for 61 submillimeter galaxies (SMGs) matched with optical/near-infrared data, and found that SMGs at S/N $>$ 4 have a median value of $z_{phot} = $2.22 $\pm$ 0.12, star formation rates of 300 $\pm$ 40 M$_{\odot}\,{\rm yr^{-1}}$ (Chabrier initial mass function), and typical cold dust masses of 5.9 $\pm$ 0.7 $ \times$ 10$^{8} $M$_{\odot}$, in line with bright SMGs from other surveys. The large cold dust masses indicate correspondingly large cool gas masses, which we suggest are a key factor necessary to drive the high star formation rates seen in this population, Comment: 22 pages, 18 figures (1 extended figure and the full catalog are included), Published in ApJS

We exploit the unprecedented depth of integral field data from the KMOS Ultra-deep Rotational Velocity Survey (KURVS) to analyse the strong (H$\alpha$) and forbidden ([NII], [SII]) emission line ratios in 22 main-sequence galaxies at $z\approx1.5$. Using the [NII]/H$\alpha$ emission-line ratio we confirm the presence of the stellar mass $-$ gas-phase metallicity relation at this epoch, with galaxies exhibiting on average 0.13$\pm$0.04 dex lower gas-phase metallicity (12+log(O/H)$_{\rm M13}$=8.40$\pm$0.03) for a given stellar mass ($\log_{10}$($M_{\rm *}$[$M_{\odot}$]=10.1$\pm$0.1) than local main-sequence galaxies. We determine the galaxy-integrated [SII] doublet ratio, with a median value of [SII]$\lambda$6716/$\lambda$6731=1.26$\pm$0.14 equivalent to an electron density of log$_{10}$($n_{\rm e}$[cm$^{-3}$])=1.95$\pm$0.12. Utilising CANDELS $HST$ multi-band imaging we define the pixel surface-mass and star-formation rate density in each galaxy and spatially resolve the fundamental metallicity relation at $z\approx1.5$, finding an evolution of 0.05$\pm$0.01 dex compared to the local relation. We quantify the intrinsic gas-phase metallicity gradient within the galaxies using the [NII]/H$\alpha$ calibration, finding a median annuli-based gradient of $\Delta$Z/$\Delta$R=$-$0.015$\pm$0.005 dex kpc$^{-1}$. Finally we examine the azimuthal variations in gas-phase metallicity, which show a negative correlation with the galaxy integrated star-formation rate surface density ($r_{\rm s}$=$-$0.40, $p_{\rm s}$=0.07) but no connection to the galaxies kinematic or morphological properties nor radial variations in stellar mass surface density or star formation rate surface density. This suggests both the radial and azimuthal variations in interstellar medium properties are connected to the galaxy integrated density of recent star formation., Comment: Submitted for publication in MNRAS. Comments welcome. 18 pages

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., Comment: 30 pages, 16 figures, ApJ in press

We present detailed stellar specific angular momentum ($j_*$) measurements of ten star-forming galaxies at $z\sim1.5-2$ using both high and low spatial resolution integral field spectroscopic data. We developed a code that simultaneously models the adaptive optics (AO) assisted observations from OSIRIS/SINFONI along with their natural seeing (NS) counterparts from KMOS at spatial resolutions of [$0.1-0.4$] arcsec and [$0.6-1.0$] arcsec respectively. The AO data reveals 2/10 systems to be mergers and for the remaining eight the mean uncertainties $\bar \Delta j_*$ decrease from 49% (NS), and 26.5% (AO), to 16% in the combined analysis. These $j_*$ measurements agree within 20% with simple estimates ($\tilde{j_*}$) calculated from the Hubble Space Telescope photometry and NS kinematics, however higher resolution kinematics are required to first identify these disks. We find that the choice of surface mass density model and the measurement of effective radius from photometry are the key sources of systematic effects in the measurement of $j_*$ between different analyses. Fitting the $j_*$ vs $M_*$ relations (Fall, 1983) with a fixed power-law slope of $\beta=2/3$, we find a zero-point consistent with prior NS results at $z\geq1$ within $\sim 0.3$ dex. Finally, we find a $\sim 0.38$ dex scatter about that relation that remains high despite the AO data so we conclude it is intrinsic to galaxies at $z>1$. This compares to a scatter of $\leq 0.2$ dex for disks at $z\simeq0$ pointing to a settling of the Fall relation with cosmic time., Comment: 29 pages (21 for the paper and 8 for a supplementary section), 31 figures. Published in Monthly Notices of the Royal Astronomical Society

Giant star-forming regions (clumps) are widespread features of galaxies at $z \approx 1-4$. Theory predicts that they can play a crucial role in galaxy evolution if they survive to stellar feedback for > 50 Myr. Numerical simulations show that clumps' survival depends on the stellar feedback recipes that are adopted. Up to date, observational constraints on both clumps' outflows strength and gas removal timescale are still uncertain. In this context, we study a line-emitting galaxy at redshift $z \simeq 3.4$ lensed by the foreground galaxy cluster Abell 2895. Four compact clumps with sizes $\lesssim$ 280 pc and representative of the low-mass end of clumps' mass distribution (stellar masses $\lesssim 2\times10^8\ {\rm M}_\odot$) dominate the galaxy morphology. The clumps are likely forming stars in a starbursting mode and have a young stellar population ($\sim$ 10 Myr). The properties of the Lyman-$\alpha$ (Ly$\alpha$) emission and nebular far-ultraviolet absorption lines indicate the presence of ejected material with global outflowing velocities of $\sim$ 200-300 km/s. Assuming that the detected outflows are the consequence of star formation feedback, we infer an average mass loading factor ($\eta$) for the clumps of $\sim$ 1.8 - 2.4 consistent with results obtained from hydro-dynamical simulations of clumpy galaxies that assume relatively strong stellar feedback. Assuming no gas inflows (semi-closed box model), the estimates of $\eta$ suggest that the timescale over which the outflows expel the molecular gas reservoir ($\simeq 7\times 10^8\ \text{M}_\odot$) of the four detected low-mass clumps is $\lesssim$ 50 Myr., Comment: 19 pages, 11 figures + Appendix. Published in MNRAS

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., Comment: 23 pages, 11 figures, accepted for publication in ApJ

We present Atacama Large Millimetre Array (ALMA) 2mm continuum observations of a complete and unbiased sample of 99 870micron-selected sub-millimeter galaxies (SMGs) in the Extended Chandra Deep Field South (ALESS). Our observations of each SMG reach average sensitivities of 53 microJy/beam. We measure the flux densities for 70 sources, for which we obtain a typical 870micron-to-2mm flux ratio of 14 +/- 5. We do not find a redshift dependence of this flux ratio, which would be expected if the dust emission properties of our SMGs were the same at all redshifts. By combining our ALMA measurements with existing Herschel/SPIRE observations, we construct a (biased) subset of 27 galaxies for which the cool dust emission is sufficiently well sampled to obtain precise constraints on their dust properties using simple isothermal models. Thanks to our new 2mm observations, the dust emissivity index is well-constrained and robust against different dust opacity assumptions. The median dust emissivity index of our SMGs is $\beta\simeq1.9\pm0.4$, consistent with the emissivity index of dust in the Milky Way and other local and high-redshift galaxies, as well as classical dust grain model predictions. We also find a negative correlation between the dust temperature and $\beta$, similar to low-redshift observational and theoretical studies. Our results indicate that $\beta\simeq2$ in high-redshift dusty star-forming galaxies, implying little evolution in dust grain properties between our SMGs and local dusty galaxy samples, and suggesting these high-mass and high-metallicity galaxies have dust reservoirs driven by grain growth in their ISM., Comment: 31 pages, 20 figures, accepted for publication in the Astrophysical Journal

We present the KMOS Galaxy Evolution Survey (KGES), a $K$-band Multi-Object Spectrograph (KMOS) study of the H$\alpha$ and [NII] emission from 288 $K$ band-selected galaxies at $1.2 \lesssim z \lesssim 1.8$, with stellar masses in the range $\log_{10}(M_{*}/\rm{M}_{\odot})\approx$9-11.5. In this paper, we describe the survey design, present the sample, and discuss the key properties of the KGES galaxies. We combine KGES with appropriately matched samples at lower redshifts from the KMOS Redshift One Spectroscopic Survey (KROSS) and the SAMI Galaxy Survey. Accounting for the effects of sample selection, data quality, and analysis techniques between surveys, we examine the kinematic characteristics and angular momentum content of star-forming galaxies at $z\approx1.5$, $\approx1$ and $\approx0$. We find that stellar mass, rather than redshift, most strongly correlates with the disc fraction amongst star-forming galaxies at $z \lesssim 1.5$, observing only a modest increase in the prevalence of discs between $z\approx1.5$ and $z\approx0.04$ at fixed stellar mass. Furthermore, typical star-forming galaxies follow the same median relation between specific angular momentum and stellar mass, regardless of their redshift, with the normalisation of the relation depending more strongly on how disc-like a galaxy's kinematics are. This suggests that massive star-forming discs form in a very similar manner across the $\approx$ 10 Gyr encompassed by our study and that the inferred link between the angular momentum of galaxies and their haloes does not change significantly across the stellar mass and redshift ranges probed in this work., Comment: 24 pages, 9 figures, 3 tables. MNRAS in press

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., Comment: 13 pages, 7 figures, submitted to MNRAS

We present ~0.15'' spatial resolution imaging of SHiZELS-14, a massive (M*~10^11 M_sol), dusty, star-forming galaxy at z=2.24. Our rest-frame ~1kpc-scale, matched-resolution data comprise four different widely used tracers of star formation: the H-alpha emission line (from SINFONI/VLT), rest-frame UV continuum (from HST F606W imaging), the rest-frame far-infrared (from ALMA), and the radio continuum (from JVLA). Although originally identified by its modest H-alpha emission line flux, SHiZELS-14 appears to be a vigorously star-forming (SFR~1000 M_sol/yr) example of a submillimeter galaxy, probably undergoing a merger. SHiZELS-14 displays a compact, dusty central starburst, as well as extended emission in $\rm{H}\alpha$ and the rest-frame optical and FIR. The UV emission is spatially offset from the peak of the dust continuum emission, and appears to trace holes in the dust distribution. We find that the dust attenuation varies across the spatial extent of the galaxy, reaching a peak of at least A_H-alpha~5 in the most dusty regions, although the extinction in the central starburst is likely to be much higher. Global star-formation rates inferred using standard calibrations for the different tracers vary from ~10-1000 M_sol/yr, and are particularly discrepant in the galaxy's dusty centre. This galaxy highlights the biased view of the evolution of star-forming galaxies provided by shorter wavelength data., Comment: 18 pages, 11 figures. Accepted for publication in MNRAS

We present an analysis of the chemical abundance properties of $\approx$650 star-forming galaxies at $z \approx0.6-1.8$. Using integral-field observations from the $K$-band Multi-Object Spectrograph (KMOS), we quantify the [NII]/H$\alpha$ emission-line ratio, a proxy for the gas-phase Oxygen abundance within the interstellar medium. We define the stellar mass-metallicity relation at $z \approx0.6-1.0$ and $z \approx1.2-1.8$ and analyse the correlation between the scatter in the relation and fundamental galaxy properties (e.g. H$\alpha$ star-formation rate, H$\alpha$ specific star-formation rate, rotation dominance, stellar continuum half-light radius and Hubble-type morphology). We find that for a given stellar mass, more highly star-forming, larger and irregular galaxies have lower gas-phase metallicities, which may be attributable to their lower surface mass densities and the higher gas fractions of irregular systems. We measure the radial dependence of gas-phase metallicity in the galaxies, establishing a median, beam smearing-corrected, metallicity gradient of $ \Delta Z/ \Delta R=0.002 \pm0.004$ dex kpc$^{-1}$, indicating on average there is no significant dependence on radius. The metallicity gradient of a galaxy is independent of its rest-frame optical morphology, whilst correlating with its stellar mass and specific star-formation rate, in agreement with an inside-out model of galaxy evolution, as well as its rotation dominance. We quantify the evolution of metallicity gradients, comparing the distribution of $\Delta Z/ \Delta R$ in our sample with numerical simulations and observations at $z \approx0-3$. Galaxies in our sample exhibit flatter metallicity gradients than local star-forming galaxies, in agreement with numerical models in which stellar feedback plays a crucial role redistributing metals., Comment: Accepted for publication in MNRAS. Comments welcome. 20 pages, 12 figures. Revised version with updated introduction and FMR evolution following comments

We analyse the physical properties of 121 SNR $\geq$ 5 sub-millimetre galaxies (SMGs) from the STUDIES 450-$\mu$m survey. We model their UV-to-radio spectral energy distributions using MAGPHYS+photo-$z$ and compare the results to similar modelling of 850-$\mu$m-selected SMG sample from AS2UDS, to understand the fundamental physical differences between the two populations at the observed depths. The redshift distribution of the 450-$\mu$m sample has a median of $z$ = 1.85 $\pm$ 0.12 and can be described by strong evolution of the far-infrared luminosity function. The fainter 450-$\mu$m sample has $\sim$14 times higher space density than the brighter 850-$\mu$m sample at $z$ $\lesssim$2, and a comparable space density at $z$ = 2-3, before rapidly declining, suggesting LIRGs are the main obscured population at $z$ $\sim$ 1-2, while ULIRGs dominate at higher redshifts. We construct rest-frame $\sim$ 180-$\mu$m-selected and dust-mass-matched samples at $z$ = 1-2 and $z$ = 3-4 from the 450-$\mu$m and 850-$\mu$m samples, respectively, to probe the evolution of a uniform sample of galaxies spanning the cosmic noon era. Using far-infrared luminosity, dust masses and an optically-thick dust model, we suggest that higher-redshift sources have higher dust densities due to inferred dust continuum sizes which are roughly half of those for the lower-redshift population at a given dust mass, leading to higher dust attenuation. We track the evolution in the cosmic dust mass density and suggest that the dust content of galaxies is governed by a combination of both the variation of gas content and dust destruction timescale., Comment: 20 pages, 7 figures, MNRAS accepted

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., Comment: 9 pages, 4 figures. Submitted to MNRAS, comments welcome

We study the radio properties of 706 sub-millimeter galaxies (SMGs) selected at 870$\mu$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\sigma$ 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 through parameter $q_\text{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_\text{IR} = 2.20\pm0.03$ for the full sample, independent of redshift, which places these $z\sim2.5$ dusty star-forming galaxies $0.44\pm0.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\gtrsim0.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 far-infrared/radio correlation. Overall, our findings paint SMGs as a homogeneous population of galaxies, as illustrated by their tight and non-evolving far-infrared/radio correlation., Comment: Accepted for publication in ApJ; 18 pages + appendices; 7 figures in main text

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., Comment: Submitted to MNRAS

Deep surveys with Atacama Large Millimeter Array (ALMA) have uncovered a population of dusty star-forming galaxies which are faint or even undetected at optical to near-infrared wavelengths. Their faintness at short wavelengths makes detailed characterization of the population challenging. Here we present a spectroscopic redshift identification and characterization of one of such near-infrared-dark galaxy discovered by an ALMA deep survey. Detection of [CI](1-0) and CO(4-3) emission lines determines the precise redshift of the galaxy, ADF22.A2, to be z=3.9913+/-0.0008. On the basis of multi-wavelength analysis, ADF22.A2 is found to be a massive, star-forming galaxy with stellar mass Mstar = $1.1_{-0.6}^{+1.3}$ x 10^{11} Msun and SFR = $430_{-150}^{+230}$ Msun/yr. The molecular gas mass is derived to be M ($H_2$) = 5.9 +/- 1.5x10^{10} Msun, indicating a gas fraction of ~35%, and the ratios of $L_{\rm [CI](1-0)}/L_{\rm IR}$ and $L_{\rm [CI](1-0)}/L_{\rm CO(4-3)}$ suggests that the nature of the interstellar medium in ADF22.A2 is in accordance with those of other bright submillimeter galaxies. The properties of ADF22.A2, including redshift, star-formation rate, stellar mass, and depletion time scale (tau~0.1-0.2 Gyr), also suggest that ADF22.A2 has the characteristics expected for the progenitors of quiescent galaxies at z>3. Our results demonstrate the power of ALMA contiguous mapping and line scan to obtain an unbiased view of galaxy formation in the early Universe., Comment: Accepted for publication in A&A Letters, 7 pages, 7 figures, 1 table

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., Comment: 21 pages, 12 colour figures, 3 tables -- accepted for publication in MNRAS

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

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., Comment: Submitted to MNRAS, comments welcome. Figures 1, 4, A1 degraded to comply with size limits

Using data from ALMA and near-infrared (NIR) integral field spectrographs including both SINFONI and KMOS on the VLT, we investigate the two-dimensional distributions of H$\alpha$ and rest-frame far-infrared (FIR) continuum in six submillimeter galaxies at $z\sim2$. At a similar spatial resolution ($\sim$0.5" FWHM; $\sim$4.5 kpc at $z=2$), we find that the half-light radius of H$\alpha$ is significantly larger than that of the FIR continuum in half of the sample, and on average H$\alpha$ is a median factor of $2.0\pm0.4$ larger. Having explored various ways to correct for the attenuation, we find that the attenuation-corrected H$\alpha$-based SFRs are systematically lower than the IR-based SFRs by at least a median factor of $3\pm1$, which cannot be explained by the difference in half-light radius alone. In addition, we find that in 40% of cases the total $V$-band attenuation ($A_V$) derived from energy balance modeling of the full ultraviolet(UV)-to-FIR spectral energy distributions (SEDs) is significantly higher than that derived from SED modeling using only the UV-to-NIR part of the SEDs, and the discrepancy appears to increase with increasing total infrared luminosity. Finally, considering all our findings along with the studies in the literature, we postulate that the dust distributions in SMGs, and possibly also in less IR luminous $z\sim2$ massive star-forming galaxies, can be decomposed into three main components; the diffuse dust heated by older stellar populations, the more obscured and extended young star-forming HII regions, and the heavily obscured central regions that have a low filling factor but dominate the infrared luminosity in which the majority of attenuation cannot be probed via UV-to-NIR emissions., Comment: 15 pages, 8 figures, A&A in press

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., Comment: 17 pages, 12 figures. Resubmitted to MNRAS after referee report

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., Comment: Published in ApJ July 2019

We present an analysis of the gas dynamics of star-forming galaxies at z~1.5 using data from the KMOS Galaxy Evolution Survey (KGES). We quantify the morphology of the galaxies using $HST$ CANDELS imaging parametrically and non-parametrically. We combine the H$\alpha$ dynamics from KMOS with the high-resolution imaging to derive the relation between stellar mass (M$_{*}$) and stellar specific angular momentum (j$_{*}$). We show that high-redshift star-forming galaxies at z~1.5 follow a power-law trend in specific stellar angular momentum with stellar mass similar to that of local late-type galaxies of the form j$_*$$\propto$M$_*^{0.53 \pm 0.10}$. The highest specific angular momentum galaxies are mostly disc-like, although generally, both peculiar morphologies and disc-like systems are found across the sequence of specific angular momentum at a fixed stellar mass. We explore the scatter within the j$_{*}$-M$_{*}$ plane and its correlation with both the integrated dynamical properties of a galaxy (e.g. velocity dispersion, Toomre Q$_{\rm g}$, H$\alpha$ star formation rate surface density $\Sigma_{\rm SFR}$) and its parameterised rest-frame UV/optical morphology (e.g. S\'ersic index, bulge to total ratio, Clumpiness, Asymmetry and Concentration). We establish that the position in the j$_{*}$-M$_{*}$ plane is strongly correlated with the star-formation surface density and the Clumpiness of the stellar light distribution. Galaxies with peculiar rest-frame UV/optical morphologies have comparable specific angular momentum to disc-dominated galaxies of the same stellar mass, but are clumpier and have higher star-formation rate surface densities. We propose that the peculiar morphologies in high--redshift systems are driven by higher star formation rate surface densities and higher gas fractions leading to a more clumpy inter-stellar medium., Comment: Accepted for publication in MNRAS. Comments welcome. 25 pages, 14 figures. The full pdf, which includes an additional 17 pages of supplementary figures is available on this http://astro.dur.ac.uk/KROSS/publications.html

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., Comment: 33 pages, 16 figures, MNRAS accepted version

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

Cosmological simulations predict the Universe contains a network of intergalactic gas filaments, within which galaxies form and evolve. However, the faintness of any emission from these filaments has limited tests of this prediction. We report the detection of rest-frame ultraviolet Lyman-alpha radiation from multiple filaments extending more than one megaparsec between galaxies within the SSA 22 proto-cluster at a redshift of 3.1. Intense star formation and supermassive black-hole activity is occurring within the galaxies embedded in these structures, which are the likely sources of the elevated ionizing radiation powering the observed Lyman-alpha emission. Our observations map the gas in filamentary structures of the type thought to fuel the growth of galaxies and black holes in massive proto-clusters., Comment: Published in Science on 2019 October 3 (Main: 3 figures, Supp: 10 figures, 2 tables)

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 S/N ratios of $>8$ in high-resolution $0.18''$ ($\sim1$kpc) 870$\mu$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 ($R_e$) of $0.10''\pm0.04''$ for our sample with a median flux of $S_{870}=5.6\pm0.2$mJy. We find that the apparent axial ratio ($b/a$) distribution of the SMGs peaks at $b/a\sim0.63\pm0.24$ and is best described by triaxial morphologies, while their emission profiles are best fit by a Sersic model with $n\simeq1.0\pm0.1$, similar to exponential discs. This combination of triaxiality and $n\sim1$ Sersic index are characteristic of bars and we suggest that the bulk of the 870$\mu$m dust continuum emission in the central $\sim2$kpc of these galaxies arises from bar-like structures. By stacking our 870$\mu$m maps we recover faint extended dust continuum emission on $\sim4$kpc scales which contributes $13\pm1$% of the total 870$\mu$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 $\simeq0.15''\pm0.05''$ or $1.2\pm0.4$kpc. Our results suggest that kpc-scale bars are ubiquitous features of high star-formation rate systems at $z\gg1$, 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., Comment: 14 pages, 10 figures, accepted for publication in MNRAS

We present Atacama Large Millimeter/submillimeter Array (ALMA) CO($J=2-1$) observations of two main-sequence star-forming galaxies at $z\sim1.47$ taken from the High-Z Emission Line Survey (HiZELS). These two systems have been previously reported to be molecular gas rich $f_{\rm H_2} \equiv M_{\rm H_2}/(M_{\rm H_2} + M_\star) \sim 0.8$. We carried out a follow-up study to resolve, at $\sim$kpc-scales, the CO emission. These new observations are combined with our earlier ALMA observations (sensitive to diffuse CO emission) and compared with our previous H$\alpha$-based study at matched spatial resolution. One target is marginally resolved in CO(2-1), showing complex dynamics with respect to the ionised gas traced by H$\alpha$. While the other source is spatially resolved, enabling a detailed exploration of its internal dynamical properties. In this system, both gaseous phases show similar spatial extension, rotation velocities and velocity dispersions ($V_{\rm rot} \sim \sigma_v \sim 100$\,km\,s$^{-1}$) suggesting a rotational velocity to velocity dispersion ratio consistent with unity. By comparing the ionized and molecular gas tracers through the use of a two-dimensional kinematic model, we estimate a median depletion time $\tau_{\rm dep}=2.3 \pm 1.2$\,Gyr for the galaxy as a whole. This value is in agreement with the average $\tau_{\rm dep}$ value observed in local star-forming galaxies at similar spatial scales. Using a thick-disk dynamical modelling, we derive a dynamical mass $M_{\rm dyn} = (1.59\pm0.19) \times 10^{11}$\,$M_\odot$ within $\approx 6$\,kpc. This suggests a dark matter fraction ($f_{\rm DM} \equiv M_{\rm DM}/M_{\rm dyn}$) of $0.59\pm0.10$, in agreement with the average $f_{\rm DM}$ value derived from stacked rotation curve analysis of galaxies at similar redshift range., Comment: 15 Pages, 9 Figures and 4 Tables. Accepted for publication in MNRAS

We combine high-resolution ALMA and HST/CANDELS observations of 20 submillimeter galaxies (SMGs) predominantly from the AS2UDS survey at z~2 with bright rest-frame optical counterparts (Ks < 22.9) to investigate the resolved structural properties of their dust and stellar components. We derive two-dimensional stellar-mass distributions that are inferred from spatial mass-to-light ratio (M/L) corrections based on rest-frame optical colors. Due to the high central column densities of dust in our SMGs, our mass distributions likely represent a lower limit to the true central mass density. The centroid positions between the inferred stellar-mass and the dust distributions agree within 1.1 kpc, indicating an overall good spatial agreement between the two components. The majority of our sources exhibit compact dust configurations relative to the stellar component (with a median ratio of effective radii Re,dust/Re,Mstar = 0.6). This ratio does not change with specific star-formation rate (sSFR) over the factor of 30 spanned by our targets, sampling the locus of "normal" main sequence galaxies up to the starburst regime, log(sSFR/sSFRMS) > 0.5. Our results imply that massive SMGs are experiencing centrally enhanced star formation unlike typical spiral galaxies in the local Universe. The sizes and stellar densities of our SMGs are in agreement with those of the passive population at z=1.5, consistent with these systems being the descendants of z~2 SMGs., Comment: 26 pages, 11 figures, accepted for publication in ApJ

We analyse the multi-frequency radio spectral properties of $41$ 6GHz-detected ALMA-identified, submillimetre galaxies (SMGs), observed at 610MHz, 1.4GHz, 6GHz with GMRT and the VLA. Combining high-resolution ($\sim0.5''$) 6GHz radio and ALMA $870\,\mu$m imaging (tracing rest-frame $\sim20$GHz, and $\sim250\,\mu$m dust continuum), we study the far-infrared/radio correlation via the logarithmic flux ratio $q_{\rm IR}$, measuring $\langle q_{\rm IR}\rangle=2.20\pm 0.06$ for our sample. We show that the high-frequency radio sizes of SMGs are $\sim1.9\pm 0.4\times$ ($\sim2$-$3$kpc) larger than those of the cool dust emission, and find evidence for a subset of our sources being extended on $\sim 10$kpc scales at 1.4GHz. 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 show evidence of a spectral break, with a median $\langle\alpha^{1.4\,{\rm GHz}}_{610\,{\rm GHz}}\rangle=-0.60\pm 0.06$, but $\langle\alpha^{6\,{\rm GHz}}_{1.4\,{\rm 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\sim35\,\mu$G) and young ages ($t_{\rm SB}\sim40$--$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 half-way through their current episode of star formation, this implies starburst durations of $\lesssim100$Myr, in reasonable agreement with estimates derived via gas depletion timescales., Comment: 25 pages, 10 colour figures, 2 tables; accepted for publication in ApJ

We analyse 850um continuum observations of eight massive X-ray detected galaxy clusters at z~0.8-1.6 taken with SCUBA-2 on the James Clerk Maxwell Telescope. We find an average overdensity of 850um-selected sources of a factor of 4+/-2 per cluster within the central 1Mpc compared to the field. We investigate the multiwavelength properties of these sources and identify 34 infrared counterparts to 26 SCUBA-2 sources. Their colours suggest that the majority of these counterparts are probable cluster members. We use the multi-wavelength far-infrared photometry to measure the total luminosities and total cluster star-formation rates demonstrating that they are roughly three orders of magnitude higher than local clusters. We predict the H-band luminosities of the descendants of our cluster submillimetre galaxies and find that their stellar luminosity distribution is consistent with that of passive elliptical galaxies in z~0 clusters. Together, the faded descendants of the passive cluster population already in place at z~1 and the cluster submillimetre galaxies are able to account for the total luminosity function of early-type cluster galaxies at z~0. This suggests that the majority of the luminous passive population in z~0 clusters are likely to have formed at z>>1 through an extreme, dust-obscured starburst event., Comment: Accepted for publication in MNRAS. 14 pages, 9 figures, 3 tables

We analyze the interstellar medium properties of a sample of sixteen bright CO line emitting galaxies identified in the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Large Program. This CO$-$selected galaxy sample is complemented by a couple of additional CO line emitters in the UDF that are identified based on their MUSE optical spectroscopic redshifts. The ASPECS CO$-$selected galaxies cover a larger range of star-formation rates and stellar masses compared to literature CO emitting galaxies at $z>1$ for which scaling relations have been established previously. Most of ASPECS CO-selected galaxies follow these established relations in terms of gas depletion timescales and gas fractions as a function of redshift, as well as the star-formation rate-stellar mass relation (`galaxy main sequence'). However, we find that $\sim30\%$ of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with $\sim12\%$ (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low star-formation rates, and yet show substantial molecular gas reservoirs, yielding long gas depletion timescales. Capitalizing on the well-defined cosmic volume probed by our observations, we measure the contribution of galaxies above, below, and on the galaxy main sequence to the total cosmic molecular gas density at different lookback times. We conclude that main sequence galaxies are the largest contributor to the molecular gas density at any redshift probed by our observations (z$\sim$1$-$3). The respective contribution by starburst galaxies above the main sequence decreases from z$\sim$2.5 to z$\sim$1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequence., Comment: Accepted for publication in ApJ. Only minor differences respect to previous version

The ALMA SPECtroscopic Survey in the {\it Hubble} Ultra Deep Field is an ALMA large program that obtained a frequency scan in the 3\,mm band to detect emission lines from the molecular gas in distant galaxies. We here present our search strategy for emission lines and continuum sources in the HUDF. We compare several line search algorithms used in the literature, and critically account for the line-widths of the emission line candidates when assessing significance. We identify sixteen emission lines at high fidelity in our search. Comparing these sources to multi-wavelength data we find that all sources have optical/infrared counterparts. Our search also recovers candidates that have lower significance that can be used statistically to derive, e.g. the CO luminosity function. We apply the same detection algorithm to obtain a sample of six 3 mm continuum sources. All of these are also detected in the 1.2 mm continuum with optical/near-infrared counterparts. We use the continuum sources to compute 3 mm number counts in the sub-mJy regime, and find them to be higher by an order of magnitude than expected for synchrotron-dominated sources. However, the number counts are consistent with those derived at shorter wavelengths (0.85--1.3\,mm) once extrapolating to 3\,mm with a dust emissivity index of $\beta=1.5$, dust temperature of 35\,K and an average redshift of $z=2.5$. These results represent the best constraints to date on the faint end of the 3 mm number counts., Comment: Accepted for publication in ApJ. 28 Pages and 13 Figures

We present adaptive optics assisted integral field spectroscopy of 34 star-forming galaxies at $z$ = 0.8-3.3 selected from the HiZELS narrow-band survey. We measure the kinematics of the ionised interstellar medium on $\sim$1 kpc scales, and show that the galaxies are turbulent, with a median ratio of rotational to dispersion support of $v$/$\sigma$=0.82$\pm$0.13. We combine the dynamics with high-resolution rest-frame optical imaging and extract emission line rotation curves. We show that high-redshift star-forming galaxies follow a similar power-law trend in specific angular momentum with stellar mass as that of local late type galaxies. We exploit the high resolution of our data and examine the radial distribution of angular momentum within each galaxy by constructing total angular momentum profiles. Although the stellar mass of a typical star-forming galaxy is expected to grow by a factor $\sim$8 in the $\sim$5 Gyrs between $z$$\sim$3.3 and $z$$\sim$0.8, we show that the internal distribution of angular momentum becomes less centrally concentrated in this period i.e the angular momentum grows outwards. To interpret our observations, we exploit the EAGLE simulation and trace the angular momentum evolution of star forming galaxies from $z$$\sim$3 to $z$$\sim$0, identifying a similar trend of decreasing angular momentum concentration. This change is attributed to a combination of gas accretion in the outer disk, and feedback that preferentially arises from the central regions of the galaxy. We discuss how the combination of the growing bulge and angular momentum stabilises the disk and gives rise to the Hubble sequence., Comment: 31 pages, 18 figures including appendices. Accepted for publication in MNRAS 2019 March 12. Updated text and figures

We present the catalogue and properties of sources in AS2UDS, an 870-$\mu$m continuum survey with the Atacama Large Millimetre/sub-millimetre Array (ALMA) of 716 single-dish sub-millimetre sources detected in the UKIDSS/UDS field by the SCUBA-2 Cosmology Legacy Survey. In our sensitive ALMA follow-up observations we detect 708 sub-millimetre galaxies (SMGs) at $>$\,4.3$\sigma$ significance across the $\sim$\,1-degree diameter field. We combine our precise ALMA positions with the extensive multi-wavelength coverage in the UDS field to fit the spectral energy distributions of our SMGs to derive a median redshift of $z_{\rm phot}=$\,2.61$\pm$0.09. This large sample reveals a statistically significant trend of increasing sub-millimetre flux with redshift suggestive of galaxy downsizing. 101 ALMA maps do not show a $>$\,4.3$\sigma$ SMG, but we demonstrate from stacking {\it Herschel} SPIRE observations at these positions, that the vast majority of these blank maps correspond to real single-dish sub-millimetre sources. We further show that these blank maps contain an excess of galaxies at $z_{\rm phot}=$\,1.5--4 compared to random fields, similar to the redshift range of the ALMA-detected SMGs. In addition, we combine X-ray and mid-infrared active galaxy nuclei activity (AGN) indicators to yield a likely range for the AGN fraction of 8--28\,\% in our sample. Finally, we compare the redshifts of this population of high-redshift, strongly star-forming galaxies with the inferred formation redshifts of massive, passive galaxies being found out to $z\sim$\,2, finding reasonable agreement -- in support of an evolutionary connection between these two classes of massive galaxy., Comment: Submitted to MNRAS; Comments welcome. Full catalogue will be made publicly available on acceptance of paper

We present first results from Multi Unit Spectroscopic Explorer (MUSE) observations at the Very Large Telescope in the MUSE Ultra Deep Field (MUDF), a $\approx 1.2\times 1.4$ arcmin$^2$ region for which we are collecting $\approx$200 hours of integral field spectroscopy. The $\approx 40$-hour observation completed to date reveals the presence of a group of three Ly$\alpha$ nebulae associated with a bright quasar pair at $z\simeq3.23$ with projected separation of $\approx 500\rm~kpc$. Two of the nebulae are physically associated with the quasars which are likely powering the Ly$\alpha$ emission, and extend for $\gtrsim 100~\rm kpc$ at a surface brightness level of $\approx 6\times 10^{-19}~\rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$. A third smaller ($\approx$35 kpc) nebula lies at a velocity offset of $\approx 1550$ km s$^{-1}$. Despite their clustered nature, the two large nebulae have properties similar to those observed in isolated quasars and exhibit no sharp decline in flux at the current depth, suggesting an even more extended distribution of gas around the quasars. We interpret the shape and the alignment of the two brighter nebulae as suggestive of the presence of an extended structure connecting the two quasar host galaxies, as seen for massive galaxies forming within gas-rich filaments in cosmological simulations., Comment: MNRAS Letters in press, 3 figures, revised version: minor changes following proofs

Using ALMA, we report high angular-resolution observations of the redshift z=3.63 galaxy, G09v1.97, one of the most luminous strongly lensed galaxies discovered by the H-ATLAS survey. We present 0"2-0"4 resolution images of the rest-frame 188 and 419$\mu$m dust continuum and the CO(6-5), H2O(211-202) and J=2 H2O+ line emission. We also report the detection of H$_2^{18}$O in this source. The dust continuum and molecular gas emission are resolved into a nearly complete ~1"5 diameter Einstein ring plus a weaker image in the center, which is caused by a special dual deflector lensing configuration. The observed line profiles of the CO, H2O and H2O+ lines are strikingly similar. In the source plane, we reconstruct the dust continuum images and the spectral cubes of the line emission at sub-kpc scales. The reconstructed dust emission in the source plane is dominated by a compact disk with an effective radius of 0.7kpc plus an overlapping extended disk with a radius twice as large. While the average magnification for the dust continuum is $\mu$~10-11, the magnification of the line emission varies 5 to 22 across different velocity components. The emission lines have similar spatial and kinematic distributions. The molecular gas and dust content reveal that G09v1.97 is a gas-rich major merger in its pre-coalescence phase. Both of the merging companions are intrinsically ULIRGs with LIR reaching $\gtrsim 4\times10^{12}L_\odot$, and the total LIR of G09v1.97 is $1.4\times10^{13}L_\odot$. The approaching southern galaxy shows no obvious kinematic structure with a semi-major half-light radius a_s=0.4kpc, while the receding galaxy resembles an a_s=1.2kpc rotating disk. The two galaxies are separated by a projected distance of 1.3kpc, bridged by weak line emission that is co-spatially located with the cold-dust-emission peak, suggesting a large amount of cold ISM in the interacting region. (abridged), Comment: 19 pages, 12 figures, and 4-page appendix, A&A in press. Minor changes to match language editing

Active galactic nuclei play a crucial role in the accretion and ejection of gas in galaxies. Although their outflows are well studied, finding direct evidence of accretion has proved very difficult and has so far been done for very few sources. A promising way to study the significance of cold accretion is by observing the absorption of an active galactic nucleus's extremely bright radio emission by the cold gas lying along the line-of-sight. As such, we present ALMA CO(1-0) and CO(2-1) observations of the Hydra-A brightest cluster galaxy (z=0.054) which reveal the existence of cold, molecular gas clouds along the line-of-sight to the galaxy's extremely bright and compact mm-continuum source. They have apparent motions relative to the central supermassive black hole of between -43 and -4 km s$^{-1}$ and are most likely moving along stable, low ellipticity orbits. The identified clouds form part of a $\sim$$10^{9}$ $\text{M}_{\odot}$, approximately edge-on disc of cold molecular gas. With peak CO(2-1) optical depths of $\tau$=0.88 $^{+0.06}_{-0.06}$, they include the narrowest and by far the deepest absorption of this type which has been observed to date in a brightest cluster galaxy. By comparing the relative strengths of the lines for the most strongly absorbing region, we are able to estimate a gas temperature of $42^{+25}_{-11}$ K and line-of-sight column densities of $N_{CO}=2^{+3}_{-1}\times 10 ^{17} cm^{-2}$ and $N_{ H_{2} }=7^{+10}_{-4}\times 10 ^{20} cm^{-2}$., Comment: 11 pages, 8 figures, accepted to MNRAS

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