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

LineStacker is a new open access and open source tool for stacking of spectral lines in interferometric data. LineStacker is an ensemble of CASA tasks, and can stack both 3D cubes or already extracted spectra. The algorithm is tested on increasingly complex simulated data sets, mimicking Atacama Large Millimeter/submillimeter Array and Karl G. Jansky Very Large Array observations of [CII] and CO(3-2) emission lines, from $z\sim7$ and $z\sim4$ galaxies respectively. We find that the algorithm is very robust, successfully retrieving the input parameters of the stacked lines in all cases with an accuracy $\gtrsim90$\%. However, we distinguish some specific situations showcasing the intrinsic limitations of the method. Mainly that high uncertainties on the redshifts ($\Delta z > 0.01$) can lead to poor signal to noise ratio improvement, due to lines being stacked on shifted central frequencies. Additionally we give an extensive description of the embedded statistical tools included in LineStacker: mainly bootstrapping, rebinning and subsampling. Velocity rebinning {is applied on the data before stacking and} proves necessary when studying line profiles, in order to avoid artificial spectral features in the stack. Subsampling is useful to sort the stacked sources, allowing to find a subsample maximizing the searched parameters, while bootstrapping allows to detect inhomogeneities in the stacked sample. LineStacker is a useful tool for extracting the most from spectral observations of various types.

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.

The gas content of galaxies is a key factor for their growth, starting from star formation and black hole accretion to galaxy mergers. Thus, characterising its properties via observations of tracers like the CO emission line is of big importance in order to understand the bigger picture of galaxy evolution. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of dust continuum, CO(5-4) and CO(8-7) line emission in the quasar--star-forming companion system SMM J04135+10277 (z=2.84). Earlier low-$J$ CO studies of this system found a huge molecular gas reservoir associated to the companion galaxy, while the quasar appeared gas-poor. Our CO observations revealed that the host galaxy of the quasar is also gas-rich, with an estimated molecular gas mass of $\sim(0.7-2.3)\times10^{10}$ M$_{\odot}$. The CO line profiles of the companion galaxy are broad ($\sim1000$ km s$^{-1}$), and show signs of rotation of a compact, massive system. In contrast to previous far-infrared observations, we resolve the continuum emission and detect both sources, with the companion galaxy dominating the dust continuum and the quasar having a $\sim25\%$ contribution to the total dust emission. By fitting the infrared spectral energy distribution of the sources with \textsc{MR-MOOSE} and empirical templates, the infrared luminosities of the quasar and the companion are in the range of $L_{\rm IR, QSO}\sim(2.1-9.6)\times10^{12}$ L$_{\odot}$ and $L_{\rm IR, Comp.}\sim(2.4-24)\times10^{12}$ L$_{\odot}$, while the estimated star formation rates are $\sim210-960$ M$_{\odot}$ yr$^{-1}$ and $\sim240-2400$ M$_{\odot}$ yr$^{-1}$, respectively. Our results demonstrate that non-detection of low-$J$ CO transition lines in similar sources does not necessarily imply the absence of massive molecular gas reservoir but that the excitation conditions favour the excitation of high-$J$ transitions.

We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$\mu$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $\tau_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr).

We present ALMA observations and multiwavelength spectral energy distribution (SED) analysis in a WISE-selected, hyperluminous dust-obscured quasar W0533-3401 at $z=2.9$. We derive its physical properties of each component, such as molecular gas, stars, dust and the central supermassive black hole (SMBH). Both the dust continuum at 3 mm and the CO(3-2) line are detected. The derived molecular gas mass $M_{\rm gas}=8.4\times10^{10}\ M_\odot$ and its fraction $f_{\rm gas}=0.7$ suggest that W0533-3401 is gas-rich. The star formation rate (SFR) has been estimated to be $\sim3000-7000\ M_\odot$ yr$^{-1}$ by using different methods. The high values of SFR and specific SFR suggest that W0533-3401 is a maximum-starburst. The corresponding gas depletion timescales are very short ($t_{\rm depl}\sim12-28$ Myr). The CO(3-2) emission line is marginally resolved and has a velocity gradient, which is possibly due to a rotating gas disk, gas outflow or merger. Finally, we infer the black hole mass growth rate of W0533-3401 (${\dot{M}}_{\rm BH}$ = 49 $M_\odot$ yr$^{-1}$), which suggests a rapid growth of the central SMBH. The observed black hole to stellar mass ratio $M_{\rm BH}/M_\star$ of W0533-3401, which is dependent on the adopted Eddington ratio, is over one order of magnitude higher than the local value, and is evolving towards the evolutionary trend of unobscured quasars. Our results are consistent with the scenario that \wobs, with both a gas-rich maximum-starburst and a rapid black hole growth, is experiencing a short transition phase towards an unobscured quasar.

Outflows in quasars during the early epochs of galaxy evolution are an important part of the feedback mechanisms potentially affecting the evolution of the host galaxy. However, systematic observations of outflows are only now becoming possible with the advent of sensitive mm telescopes. In this study we use spectral stacking methods to search for faint high velocity outflow signal in a sample of [C II] detected, z ~ 6 quasars. We search for broad emission line signatures from high-velocity outflows for a sample of 26 z ~ 6 quasars observed with ALMA, with a detection of the [C II] line. The observed emission lines of the sources are dominated by the host galaxy, and outflow emission is not detected for the individual sources. We use a spectral line stacking analysis developed for interferometric data to search for outflow emission. We stack both extracted spectra and the full spectral cubes. We also investigate the possibility that only a sub-set of our sample contributes to the stacked outflow emission. We find only a tentative detection of a broad emission line component in the stacked spectra. When taking a region of about 2 arcsec around the source central position of the stacked cubes, the stacked line shows an excess emission due to a broad component of 1.1-1.5 sigma, but the significance drops to 0.4-0.7 sigma when stacking the extracted spectra from a smaller region. The broad component can be characterised by a line width of full width half max FWHM > 700 km/s. Furthermore, we find a sub-sample of 12 sources the stack of which maximises the broad component emission. The stack of this sub-sample shows an excess emission due to a broad component of 1.2-2.5 sigma. The stacked line of these sources has a broad component of FWHM > 775 km/s. Deeper ALMA observations are necessary to confirm the presence of a broad component in the individual spectra.

We present sub-kpc-scale mapping of the 870 $\mu$m ALMA continuum emission in six luminous ($L_{\rm IR}~\sim~5~\times10^{12}$ L$_{\odot}$) submillimeter galaxies (SMGs) from the ALESS survey of the Extended Chandra Deep Field South. Our high-fidelity 0.07$''$-resolution imaging ($\sim$500 pc) reveals robust evidence for structures with deconvolved sizes of $\lesssim$0.5-1 kpc embedded within (dominant) exponential dust disks. The large-scale morphologies of the structures within some of the galaxies show clear curvature and/or clump-like structures bracketing elongated nuclear emission, suggestive of bars, star-forming rings, and spiral arms. In this interpretation, the ratio of the `ring' and `bar' radii (1.9$\pm$0.3) agrees with that measured for such features in local galaxies. These potential spiral/ring/bar structures would be consistent with the idea of tidal disturbances, with their detailed properties implying flat inner rotation curves and Toomre-unstable disks (Q<1). The inferred one-dimensional velocity dispersions ($\sigma_{\rm r}\lesssim$ 70-160 km s$^{-1}$) are marginally consistent with the limits implied if the sizes of the largest structures are comparable to the Jeans length. We create maps of the star formation rate density ($\Sigma_{\rm SFR}$) on $\sim$500 pc scales and show that the SMGs are able to sustain a given (galaxy-averaged) $\Sigma_{\rm SFR}$ over much larger physical scales than local (ultra-)luminous infrared galaxies. However, on 500 pc scales, they do not exceed the Eddington limit set by radiation pressure on dust. If confirmed by kinematics, the potential presence of non-axisymmetric structures would provide a means for net angular momentum loss and efficient star formation, helping to explain the very high star formation rates measured in SMGs.

We present new ALMA observations and physical properties of a Lyman Break Galaxy at z=7.15. Our target, B14-65666, has a bright ultra-violet (UV) absolute magnitude, $M_{\rm UV}\approx-22.4$, and has been spectroscopically identified in Ly$\alpha$ with a small rest-frame equivalent width of $\approx4$ \AA. Previous HST image has shown that the target is comprised of two spatially separated clumps in the rest-frame UV. With ALMA, we have newly detected spatially resolved [OIII] 88 $\mu$m, [CII] 158 $\mu$m, and their underlying dust continuum emission. In the whole system of B14-65666, the [OIII] and [CII] lines have consistent redshifts of $7.1520\pm0.0003$, and the [OIII] luminosity, $(34.4\pm4.1)\times10^{8}L_{\rm \odot}$, is about three times higher than the [CII] luminosity, $(11.0\pm1.4)\times10^{8}L_{\rm \odot}$. With our two continuum flux densities, the dust temperature is constrained to be $T_{\rm d}\approx50-60$ K under the assumption of the dust emissivity index of $\beta_{\rm d}=2.0-1.5$, leading to a large total infrared luminosity of $L_{\rm TIR}\approx1\times10^{12}L_{\rm \odot}$. Owing to our high spatial resolution data, we show that the [OIII] and [CII] emission can be spatially decomposed into two clumps associated with the two rest-frame UV clumps whose spectra are kinematically separated by $\approx200$ km s$^{-1}$. We also find these two clumps have comparable UV, infrared, [OIII], and [CII] luminosities. Based on these results, we argue that B14-65666 is a starburst galaxy induced by a major-merger. The merger interpretation is also supported by the large specific star-formation rate (defined as the star-formation rate per unit stellar mass), sSFR$=260^{+119}_{-57}$ Gyr$^{-1}$, inferred from our SED fitting. Probably, a strong UV radiation field caused by intense star formation contributes to its high dust temperature and the [OIII]-to-[CII] luminosity ratio.

The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure.

We present ALMA observations of the mid-J 12CO emission from six single-dish selected 870-micron sources in the Extended Chandra Deep Field-South (ECDFS) and UKIDSS Ultra-Deep Survey (UDS) fields. These six single-dish submillimetre sources were selected based on previous ALMA continuum observations, which showed that each comprised a blend of emission from two or more individual submillimetre galaxies (SMGs), separated on 5--10 arcsec scales. The six single-dish submillimetre sources targeted correspond to a total of 14 individual SMGs, of which seven have previously-measured robust optical/near-infrared spectroscopic redshifts, which were used to tune our ALMA observations. We detect CO(3-2) or CO(4-3) at z=2.3--3.7 in seven of the 14 SMGs, and in addition serendipitously detect line emission from three gas-rich companion galaxies, as well as identify four new 3.3-mm selected continuum sources in the six fields. Joint analysis of our CO spectroscopy and existing data suggests that 64 \pm 18% of the SMGs in blended submillimetre sources are unlikely to be physically associated. However, three of the SMG fields (50%) contain new, serendipitously-detected CO-emitting (but submillimetre-faint) sources at similar redshifts to the 870-micron selected SMGs we targeted. These data suggest that the SMGs inhabit overdense regions, but that these are not sufficiently overdense on ~100 kpc scales to influence the source blending given the short lifetimes of SMGs. We find that 21 \pm 12% of SMGs have spatially-distinct and kinematically-close companion galaxies (~8--150 kpc and <~300 km/s), which may have enhanced their star-formation via gravitational interactions.

We present the star formation rates (SFRs) of a sample of 109 galaxies with X-ray selected active galactic nuclei (AGN) with moderate to high X-ray luminosities (L(2-8keV)= 10^42-10^45 erg/s), at redshifts 1 < z < 4.7, that were selected to be faint or undetected in the Herschel bands. We combine our deep ALMA continuum observations with deblended 8-500{\mu}m photometry from Spitzer and Herschel, and use infrared (IR) SED fitting and AGN - star formation decomposition methods. The addition of the ALMA photometry results in an order of magnitude more X-ray AGN in our sample with a measured SFR (now 37%). The remaining 63% of the sources have SFR upper limits that are typically a factor of 2-10 times lower than the pre-ALMA constraints. With the improved constraints on the IR SEDs, we can now identify a mid-IR (MIR) AGN component in 50% of our sample, compared to only ~1% previously. We further explore the F870{\mu}m/F24{\mu}m-redshift plane as a tool for the identification of MIR emitting AGN, for three different samples representing AGN dominated, star formation dominated, and composite sources. We demonstrate that the F870{\mu}m/F24{\mu}m-redshift plane can successfully split between AGN and star formation dominated sources, and can be used as an AGN identification method.

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

We report the results of a pilot study of CO$(4-3)$ emission line of three {\it WISE}-selected hyper-luminous, dust-obscured quasars (QSOs) with sensitive ALMA Band 3 observations. These obscured QSOs with $L_{\rm bol}>10^{14}L_\odot$ are among the most luminous objects in the universe. All three QSO hosts are clearly detected both in continuum and in CO$(4-3)$ emission line. Based on CO$(4-3)$ emission line detection, we derive the molecular gas masses ($\sim 10^{10-11}$ M$_\odot$), suggesting that these QSOs are gas-rich systems. We find that three obscured QSOs in our sample follow the similar $L'_{\rm CO}- L_{\rm FIR}$ relation as unobscured QSOs at high redshifts. We also find the complex velocity structures of CO$(4-3)$ emission line, which provide the possible evidence for gas-rich merger in W0149+2350 and possible molecular outflow in W0220+0137 and W0410$-$0913. Massive molecular outflow can blow away the obscured interstellar medium (ISM) and make obscured QSOs evolve towards the UV/optical bright, unobscured phase. Our result is consistent with the popular AGN feedback scenario involving the co-evolution between the SMBH and host galaxy.

Hot dust-obscured galaxies (Hot DOGs) are a luminous, dust-obscured population recently discovered in the WISE All-Sky survey. Multiwavelength follow-up observations suggest that they are mainly powered by accreting supermassive black holes (SMBHs), lying in dense environments, and being in the transition phase between extreme starburst and UV-bright quasars. Therefore, they are good candidates for studying the interplay between SMBHs, star formation and environment. W2246$-$0526 (thereafter, W2246), a Hot DOG at $z\sim4.6$, has been taken as the most luminous galaxy known in the Universe. Revealed by the multiwavelength images, the previous Herschel SPIRE photometry of W2246 is contaminated by a foreground galaxy (W2246f), resulting in an overestimation of its total IR luminosity by a factor of about 2. We perform the rest-frame UV/optical-to-far-IR spectral energy distribution (SED) analysis with SED3FIT and re-estimate its physical properties. The derived stellar mass $M_\star = 4.3\times10^{11}~M_\odot$ makes it be among the most massive galaxies with spectroscopic redshift $z>4.5$. Its structure is extremely compact and requires an effective mechanism to puff-up. Most of ($>95\%$) its IR luminosity is from AGN torus emission, revealing the rapid growth of the central SMBH. We also predict that W2246 may have a significant molecular gas reservoir based on the dust mass estimation.

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

Hot dust-obscured galaxies (Hot DOGs) are a new population recently discovered in the \wise All-Sky survey. Multiwavelength follow-up observations suggest that they are luminous, dust-obscured quasars at high redshift. Here we present the JCMT SCUBA-2 850 $\mu m$ follow-up observations of 10 Hot DOGs. Four out of ten Hot DOGs have been detected at $>3\sigma$ level. Based on the IR SED decomposition approach, we derive the IR luminosities of AGN torus and cold dust components. Hot DOGs in our sample are extremely luminous with most of them having $L_{\rm IR}^{\rm tot}>10^{14} L_\odot$. The torus emissions dominate the total IR energy output. However, the cold dust contribution is still non-negligible, with the fraction of the cold dust contribution to the total IR luminosity $(\sim 8-24\%)$ being dependent on the choice of torus model. The derived cold dust temperatures in Hot DOGs are comparable to those in UV bright quasars with similar IR luminosity, but much higher than those in SMGs. Higher dust temperatures in Hot DOGs may be due to the more intense radiation field caused by intense starburst and obscured AGN activities. Fourteen and five submillimeter serendipitous sources in the 10 SCUBA-2 fields around Hot DOGs have been detected at $>3\sigma$ and $>3.5\sigma$ levels, respectively. By estimating their cumulative number counts, we confirm the previous argument that Hot DOGs lie in dense environments. Our results support the scenario in which Hot DOGs are luminous, dust-obscured quasars lying in dense environments, and being in the transition phase between extreme starburst and UV-bright quasars.

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

We investigate the mean star formation rates (SFRs) in the host galaxies of ~3000 optically selected QSOs from the SDSS survey within the Herschel-ATLAS fields, and a radio-luminous sub-sample, covering the redshift range of z = 0.2-2.5. Using WISE & Herschel photometry (12 - 500{\mu}m) we construct composite SEDs in bins of redshift and AGN luminosity. We perform SED fitting to measure the mean infrared luminosity due to star formation, removing the contamination from AGN emission. We find that the mean SFRs show a weak positive trend with increasing AGN luminosity. However, we demonstrate that the observed trend could be due to an increase in black hole (BH) mass (and a consequent increase of inferred stellar mass) with increasing AGN luminosity. We compare to a sample of X-ray selected AGN and find that the two populations have consistent mean SFRs when matched in AGN luminosity and redshift. On the basis of the available virial BH masses, and the evolving BH mass to stellar mass relationship, we find that the mean SFRs of our QSO sample are consistent with those of main sequence star-forming galaxies. Similarly, the radio-luminous QSOs have mean SFRs that are consistent with both the overall QSO sample and with star-forming galaxies on the main sequence. In conclusion, on average QSOs reside on the main sequence of star-forming galaxies, and the observed positive trend between the mean SFRs and AGN luminosity can be attributed to BH mass and redshift dependencies.

We present a new exploration of the cosmic star-formation history and dust obscuration in massive galaxies at redshifts $0.5< z<6$. We utilize the deepest 450 and 850$\mu$m imaging from SCUBA-2 CLS, covering 230arcmin$^2$ in the AEGIS, COSMOS and UDS fields, together with 100-250$\mu$m imaging from Herschel. We demonstrate the capability of the T-PHOT deconfusion code to reach below the confusion limit, using multi-wavelength prior catalogues from CANDELS/3D-HST. By combining IR and UV data, we measure the relationship between total star-formation rate (SFR) and stellar mass up to $z\sim5$, indicating that UV-derived dust corrections underestimate the SFR in massive galaxies. We investigate the relationship between obscuration and the UV slope (the IRX-$\beta$ relation) in our sample, which is similar to that of low-redshift starburst galaxies, although it deviates at high stellar masses. Our data provide new measurements of the total SFR density (SFRD) in $M_\ast>10^{10}M_\odot$ galaxies at $0.5<z<6$. This is dominated by obscured star formation by a factor of $>10$. One third of this is accounted for by 450$\mu$m-detected sources, while one fifth is attributed to UV-luminous sources (brighter than $L^\ast_{UV}$), although even these are largely obscured. By extrapolating our results to include all stellar masses, we estimate a total SFRD that is in good agreement with previous results from IR and UV data at $z\lesssim3$, and from UV-only data at $z\sim5$. The cosmic star-formation history undergoes a transition at $z\sim3-4$, as predominantly unobscured growth in the early Universe is overtaken by obscured star formation, driven by the build-up of the most massive galaxies during the peak of cosmic assembly.

The gravitationally-lensed galaxy A1689-zD1 is one of the most distant spectroscopically confirmed sources ($z=7.5$). It is the earliest known galaxy where the interstellar medium (ISM) has been detected; dust emission was detected with the Atacama Large Millimetre Array (ALMA). A1689-zD1 is also unusual among high-redshift dust emitters as it is a sub-L* galaxy and is therefore a good prospect for the detection of gaseous ISM in a more typical galaxy at this redshift. We observed A1689-zD1 with ALMA in bands 6 and 7 and with the Green Bank Telescope (GBT) in band $Q$. To study the structure of A1689-zD1, we map the mm thermal dust emission and find two spatial components with sizes about $0.4-1.7$\,kpc (lensing-corrected). The rough spatial morphology is similar to what is observed in the near-infrared with {\it HST} and points to a perturbed dynamical state, perhaps indicative of a major merger or a disc in early formation. The ALMA photometry is used to constrain the far-infrared spectral energy distribution, yielding a dust temperature ($T_{\rm dust} \sim 35$--$45$\,K for $\beta = 1.5-2$). We do not detect the CO(3-2) line in the GBT data with a 95\% upper limit of 0.3\,mJy observed. We find a slight excess emission in ALMA band~6 at 220.9\,GHz. If this excess is real, it is likely due to emission from the [CII] 158.8\,$\mu$m line at $z_{\rm [CII]} = 7.603$. The stringent upper limits on the [CII]/$L_{\rm FIR}$ luminosity ratio suggest a [CII] deficit similar to several bright quasars and massive starbursts.

(Abridged) In the last decade several massive molecular gas reservoirs were found <100 kpc distance from active galactic nuclei (AGNs), residing in gas-rich companion galaxies. The study of AGN-gas-rich companion systems opens the opportunity to determine whether the stellar mass of massive local galaxies was formed in their host after a merger event or outside of their host galaxy in a close starbursting companion and later incorporated via mergers. We study the quasar-companion galaxy system of SMM J04135+10277 (z=2.84) and investigate the expected frequency of quasar-starburst galaxy pairs at high redshift using a cosmological galaxy formation model. We use archive data and new APEX ArTeMiS data to construct and model the spectral energy distribution of SMM J04135. We also carry out a comprehensive analysis of the cosmological galaxy formation model GALFORM with the aim of characterising how typical the system of SMM J04135 is and whether quasar-star-forming galaxy pairs may constitute an important stage in galaxy evolution. The companion galaxy of SMM J04135 is a heavily dust-obscured starburst galaxy with a median star formation rate (SFR) of $700\,\rm{M_{\odot}\,yr^{-1}}$, median dust mass of $5.1\times 10^9\,\rm{M_{\odot}}$ and median dust luminosity of $\textrm 9.3 \times 10^{12}\,\rm{L_{\odot}}$. Our simulations, performed at z=2.8, suggest that SMM J04135 is not unique. In fact, at a distance of <100 kpc, 22% of our simulated quasar sample have at least one companion galaxy of a stellar mass $>10^8\, \rm{M_{\odot}}$, and 0.3% have at least one highly star-forming companion ($\rm{SFR}>100\,\rm{M_{\odot}\,yr^{-1}}$). Our results suggest that quasar-gas-rich companion galaxy systems are common phenomena in the early Universe and the high incidence of companions makes the study of such systems crucial to understand the growth and hierarchical build-up of galaxies and black holes.

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

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

We have observed the dust continuum of ten z=3.1 Lyman Break Galaxies with the Atacama Large Millimeter/Submillimeter Array at ~450 mas resolution in Band 7. We detect and resolve the 870um emission in one of the targets with an integrated flux density of S(870)=(192+/-57) uJy, and measure a stacked 3-sigma signal of S(870)=(67+/-23) uJy for the remaining nine. The total infrared luminosities estimated from full spectral energy distribution fits are L(8-1000um)=(8.4+/-2.3)x10^10 Lsun for the detection and L(8-1000um)=(2.9+/-0.9)x10^10 Lsun for the stack. With HST ACS I-band imaging we map the rest-frame UV emission on the same scale as the dust, effectively resolving the 'infrared excess' (IRX=L_FIR/L_UV) in a normal galaxy at z=3. Integrated over the galaxy we measure IRX=0.56+/-0.15, and the galaxy-averaged UV slope is beta=-1.25+/-0.03. This puts the galaxy a factor of ~10 below the IRX-beta relation for local starburst nuclei of Meurer et al. (1999). However, IRX varies by more than a factor of 3 across the galaxy, and we conclude that the complex relative morphology of the dust relative to UV emission is largely responsible for the scatter in the IRX-beta relation at high-z. A naive application of a Meurer-like dust correction based on the UV slope would dramatically over-estimate the total star formation rate, and our results support growing evidence that when integrated over the galaxy, the typical conditions in high-z star-forming galaxies are not analogous to those in the local starburst nuclei used to establish the Meurer relation.

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

The far-infrared fine-structure line [CII] at 1900.5\,GHz is known to be one of the brightest cooling lines in local galaxies, and therefore it has been suggested to be an efficient tracer for star-formation in very high-redshift galaxies. However, recent results for galaxies at $z>6$ have yielded numerous non-detections in star-forming galaxies, except for quasars and submillimeter galaxies. We report the results of ALMA observations of two lensed, star-forming galaxies at $z = 6.029$ and $z=6.703$. The galaxy A383-5.1 (star formation rate [SFR] of 3.2 M$_\odot$ yr$^{-1}$ and magnification of $\mu = 11.4\pm1.9$) shows a line detection with $L_{\rm [CII]} = 8.9\times10^{6}$ L$_\odot$, making it the lowest $L_{\rm [CII]}$ detection at $z>6$. For MS0451-H (SFR = 0.4 M$_\odot$ yr$^{-1}$ and $\mu = 100\pm20$) we provide an upper limit of $L_{\rm [CII]} < 3\times10^{5}$ L$_\odot$, which is 1\,dex below the local SFR-$L_{\rm [CII]}$ relations. The results are consistent with predictions for low-metallicity galaxies at $z>6$, however, other effects could also play a role in terms of decreasing $L_{\rm [CII]}$. The detection of A383-5.1 is encouraging and suggests that detections are possible, but much fainter than initially predicted.

We present CO(1-0) observations obtained at the Karl G. Jansky Very Large Array (VLA) for 14 z~2 galaxies with existing CO(3-2) measurements, including 11 galaxies which contain active galactic nuclei (AGN) and three submillimeter galaxies (SMGs). We combine this sample with an additional 15 z~2 galaxies from the literature that have both CO(1-0) and CO(3-2) measurements in order to evaluate differences in CO excitation between SMGs and AGN host galaxies, measure the effects of CO excitation on the derived molecular gas properties of these populations, and to look for correlations between the molecular gas excitation and other physical parameters. With our expanded sample of CO(3-2)/CO(1-0) line ratio measurements, we do not find a statistically significant difference in the mean line ratio between SMGs and AGN host galaxies as found in the literature, instead finding r_3,1=1.03+/-0.50 for AGN host galaxies and r_3,1=0.78+/-0.27 for SMGs (or r_3,1=0.90+/-0.40 for both populations combined). We also do not measure a statistically significant difference between the distributions of the line ratios for these populations at the p=0.05 level, although this result is less robust. We find no excitation dependence on the index or offset of the integrated Schmidt-Kennicutt relation for the two CO lines, and obtain indices consistent with N=1 for the various sub-populations. However, including low-z "normal" galaxies increases our best-fit Schmidt-Kennicutt index to N~1.2. While we do not reproduce correlations between the CO line width and luminosity, we do reproduce correlations between CO excitation and star formation efficiency.

We utilize a Bayesian approach to fit the observed mid-IR-to-submm/mm spectral energy distributions (SEDs) of 22 WISE-selected and submm-detected, hyperluminous hot dust-obscured galaxies. By adopting the Torus+GB model, we decompose the observed IR SEDs of Hot DOGs into torus and cold dust components. The main results are: 1) Hot DOGs in our submm-detected sample are hyperluminous, with torus emission dominating the IR energy output. However, cold dust emission is non-negligible, averagely contributing ~24% of total IR luminosity. 2) Compared to QSO and starburst SED templates, the median SED of Hot DOGs shows the highest luminosity ratio between mid-IR and submm at rest-frame, while it is very similar to that of QSOs at 10-50um suggesting that the heating sources of Hot DOGs should be buried AGNs. 3) Hot DOGs have both high dust temperatures ~73K and IR luminosity of cold dust. The T-L relation of Hot DOGs suggests that the increase in IR luminosity for Hot DOGs is mostly due to the increase of the dust temperature, rather than dust mass. Hot DOGs have lower dust masses than those of submillimeter galaxies (SMGs) and QSOs within the similar redshift range. Both high IR luminosity of cold dust and relatively low dust mass in Hot DOGs can be expected by their relatively high dust temperatures. 4) Hot DOGs have high dust covering factors, which deviate the previously proposed trend of the dust covering factor decreasing with increasing bolometric luminosity. Finally, we can reproduce the observed properties in Hot DOGs by employing a physical model of galaxy evolution. The result suggests that Hot DOGs may lie at or close to peaks of both star formation and black hole growth histories, and represent a transit phase during the evolution of massive galaxies, transforming from the dusty starburst dominated phase to the optically bright QSO phase. (abridged)

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

Candidates for the modest galaxies that formed most of the stars in the early universe, at redshifts $z > 7$, have been found in large numbers with extremely deep restframe-UV imaging. But it has proved difficult for existing spectrographs to characterise them in the UV. The detailed properties of these galaxies could be measured from dust and cool gas emission at far-infrared wavelengths if the galaxies have become sufficiently enriched in dust and metals. So far, however, the most distant UV-selected galaxy detected in dust emission is only at $z = 3.25$, and recent results have cast doubt on whether dust and molecules can be found in typical galaxies at this early epoch. Here we report thermal dust emission from an archetypal early universe star-forming galaxy, A1689-zD1. We detect its stellar continuum in spectroscopy and determine its redshift to be $z = 7.5\pm0.2$ from a spectroscopic detection of the Ly{\alpha} break. A1689-zD1 is representative of the star-forming population during reionisation, with a total star-formation rate of about 12M$_\odot$ yr$^{-1}$. The galaxy is highly evolved: it has a large stellar mass, and is heavily enriched in dust, with a dust-to-gas ratio close to that of the Milky Way. Dusty, evolved galaxies are thus present among the fainter star-forming population at $z > 7$, in spite of the very short time since they first appeared.

Stacking as a tool for studying objects that are not individually detected is becoming popular even for radio interferometric data, and will be widely used in the SKA era. Stacking is typically done using imaged data rather than directly using the visibilities (the uv-data). We have investigated and developed a novel algorithm to do stacking using the uv-data. We have performed exten- sive simulations comparing to image-stacking, and summarize the results of these simulations. Furthermore, we disuss the implications in light of the vast data volume produced by the SKA. Having access to the uv-stacked data provides a great advantage, as it allows the possibility to properly analyse the result with respect to calibration artifacts as well as source properties such as size. For SKA the main challenge lies in archiving the uv-data. For purposes of robust stacking analysis, it would be strongly desirable to either keep the calibrated uv-data at least in an aver- age form, or implement a stacking queue where stacking positions could be provided prior to the observations and the uv-stacking is done almost in real time.

We present a new algorithm for stacking radio interferometric data in the uv-domain. The performance of uv-stacking is compared to the stacking of fully imaged data using simulated Atacama Large Millimeter/sub-millimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA) deep extragalactic surveys. We find that image- and uv-stacking produce similar results, however, uv-stacking is typically the more robust method. An advantage of the uv-stacking algorithm is the availability of uv-data post stacking, which makes it possible to identify and remove problematic baselines. For deep VLA surveys uv-stacking yields a signal-to-noise ratio that is up to 20 per cent higher than image-stacking. Furthermore, we have investigated stacking of resolved sources with a simulated VLA data set where 1.5" (10-12 kpc at z ~ 1-4) sources are stacked. We find that uv-stacking, where a model is fitted directly to the visibilities, significantly improves the accuracy and robustness of the size estimates. While scientific motivation for this work is studying faint, high-z galaxies, the algorithm analysed here would also be applicable in other fields of astronomy. Stacking of radio interferometric data is also expected to play a big role for future surveys with telescopes such as LOFAR and Square Kilometre Array (SKA).

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

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

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

We analyze morphologies of the host galaxies of 35 X-ray selected active galactic nucleus (AGNs) at $z\sim2$ in the Cosmic Evolution Survey (COSMOS) field using Hubble Space Telescope/WFC3 imaging taken from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). We build a control sample of 350 galaxies in total, by selecting ten non-active galaxies drawn from the same field with the similar stellar mass and redshift for each AGN host. By performing two dimensional fitting with GALFIT on the surface brightness profile, we find that the distribution of S$\`e$rsic index (n) of AGN hosts does not show a statistical difference from that of the control sample. We measure the nonparametric morphological parameters (the asymmetry index A, the Gini coefficient G, the concentration index C and the M20 index) based on point source subtracted images. All the distributions of these morphological parameters of AGN hosts are consistent with those of the control sample. We finally investigate the fraction of distorted morphologies in both samples by visual classification. Only $\sim$15% of the AGN hosts have highly distorted morphologies, possibly due to a major merger or interaction. We find there is no significant difference in the distortion fractions between the AGN host sample and control sample. We conclude that the morphologies of X-ray selected AGN hosts are similar to those of nonactive galaxies and most AGN activity is not triggered by major merger.

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

In the course of our 870um APEX/LABOCA follow up of the Herschel Lensing Survey we have detected a source in AS1063 (RXC J2248.7-4431), that has no counterparts in any of the Herschel PACS/SPIRE bands, it is a Herschel 'drop-out' with S_870/S_500>0.5. The 870um emission is extended and centered on the brightest cluster galaxy suggesting either a multiply imaged background source or substructure in the Sunyaev-Zel'dovich (SZ) increment due to inhomogeneities in the hot cluster gas of this merging cluster. We discuss both interpretations with emphasis on the putative lensed source. Based on the observed properties and on our lens model we find that this source could be the first SMG with a moderate far infrared luminosity (L_FIR<10^12 L_sol) detected so far at z>4. In deep HST observations we identified a multiply imaged z~6 source and we measured its spectroscopic redshift z=6.107 with VLT/FORS. This source could be associated with the putative SMG but it is most likely offset spatially by 10-30kpc and they could be interacting galaxies. With a FIR luminosity in the range [5-15]x10^{11} L_sol corresponding to a star formation rate in the range [80-260]M_sol/yr, this SMG would be more representative than the extreme starbursts usually detected at z>4. With a total magnification of ~25 it would open a unique window to the 'normal' dusty galaxies at the end of the epoch of reionization.

We present ALMA observations of the [C II] 158 micron fine structure line and dust continuum emission from the host galaxies of five redshift 6 quasars. We also report complementary observations of 250 GHz dust continuum and CO (6-5) line emission from the z=6.00 quasar SDSS J231038.88+185519.7. The ALMA observations were carried out in the extended array at 0.7" resolution. We have detected the line and dust continuum in all five objects. The derived [C II] line luminosities are 1.6x10^{9} to 8.8x10^{9} Lsun and the [C II]-to-FIR luminosity ratios are 3.0-5.6x10^{-4}, which is comparable to the values found in other high-redshift quasar-starburst systems and local ultra-luminous infrared galaxies. The sources are marginally resolved and the intrinsic source sizes (major axis FWHM) are constrained to be 0.3" to 0.6" (i.e., 1.7 to 3.5 kpc) for the [C II] line emission and 0.2" to 0.4" (i.e., 1.2 to 2.3 kpc) for the continuum. These measurements indicate that there is vigorous star formation over the central few kpc in the quasar host galaxies. The ALMA observations also constrain the dynamical properties of the atomic gas in the starburst nuclei. The intensity-weighted velocity maps of three sources show clear velocity gradients. Such velocity gradients are consistent with a rotating, gravitationally bound gas component, although they are not uniquely interpreted as such. Under the simplifying assumption of rotation, the implied dynamical masses within the [C II]-emitting regions are of order 10^{10} to 10^{11} Msun. Given these estimates, the mass ratios between the SMBHs and the spheroidal bulge are an order of magnitude higher than the mean value found in local spheroidal galaxies, which is in agreement with results from previous CO observations of high redshift quasars.

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

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

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

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

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

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

[Abridged] Observations of molecular gas at all redshifts are critical for measuring the cosmic evolution in molecular gas density and understanding the star-formation history of the Universe. The 12CO molecule (J=1-0 transition = 115.27 GHz) is the best proxy for extragalactic H2, which is the gas reservoir from which star formation occurs, and has been detected out to z~6. Typically, redshifted high-J lines are observed at mm-wavelengths, the most commonly targeted systems exhibiting high SFRs (e.g. submm galaxies), and far-IR-bright QSOs. While the most luminous objects are the most readily observed, detections of more typical galaxies with modest SFRs are essential for completing the picture. ALMA will be revolutionary in terms of increasing the detection rate and pushing the sensitivity limit down to include such galaxies, however the limited FoV when observing at such high frequencies makes it difficult to use ALMA for studies of the large-scale structure traced out by molecular gas in galaxies. This article introduces a strategy for a systematic search for molecular gas during the EoR (z~7 and above), capitalizing on the fact that the J=1-0 transition of 12CO enters the upper bands of cm-wave instruments at high-z. The FoV advantage gained by observing at such frequencies, coupled with modern broadband correlators allows significant cosmological volumes to be probed on reasonable timescales. In this article we present an overview of our future observing programme which has been awarded 6,500 hours as one of the Large Survey Projects for MeerKAT, the forthcoming South African SKA pathfinder instrument. Its large FoV and correlator bandwidth, and high-sensitivity provide unprecedented survey speed for such work. An existing astrophysical simulation is coupled with instrumental considerations to demonstrate the feasibility of such observations and predict detection rates.

We present radio and infrared (3.6-24um) counterparts to submillimetre galaxies (SMGs) detected in the Extended Chandra Deep Field South with the LABOCA 870-um bolometer camera on the 12-m Atacama Pathfinder Experiment. Using the Very Large Array at 1.4GHz and Spitzer we have identified secure counterparts to 79 of the 126 SMGs (SNR>3.7, S_870>4.4mJy) in the field, 62 via their radio and/or 24-um emission, the remainder using a colour-flux cut on IRAC 3.6- and 5.8-um sources chosen to maximise the number of secure, coincident radio and 24-um counterparts. In constructing our radio catalogue, we have corrected for the effects of `flux boosting', then used the corrected flux densities to estimate the redshifts of the SMGs based on the radio/submm spectral indices. The effect of the boosting correction is to increase the median redshift by 0.2 resulting in a value of z=2.2 (+0.7-0.8) (1-sigma errors) for the secure radio counterparts, in agreement with other studies, both spectroscopic and photometric.

We have detected the CO(2-1) transition from the submillimetre galaxy (SMG) LESSJ033229.4-275619 at z=4.755 using the new Compact Array Broadband Backend system on the Australian Telescope Compact Array. These data have identified a massive gas reservoir available for star formation for the first time in an SMG at z~5. We use the luminosity and velocity width (FWHM of 160 km/s) of the CO(2--1) line emission to constrain the gas and dynamical mass of Mgas~1.6x10^10 Msun and Mdyn(<2kpc)~5x10^10 (0.25/sin^2(i)) Msun, respectively, similar to that observed for SMGs at lower redshifts of z~2-4, although we note that our observed CO FWHM is a factor of ~3 narrower than typically seen in SMGs. Together with the stellar mass we estimate a total baryonic mass of Mbary~1x10^11 Msun, consistent with the dynamical mass for this young galaxy within the uncertainties. Dynamical and baryonic mass limits of high-redshift galaxies are useful tests of galaxy formation models: using the known z~4-5 SMGs as examples of massive baryonic systems, we find that their space density is consistent with that predicted by current galaxy formation models. In addition, these observations have helped to confirm that z~4-5 SMGs possess the baryonic masses and gas consumption timescales necessary to be the progenitors of the luminous old red galaxies seen at z~3. Our results provide a preview of the science that ALMA will enable on the formation and evolution of the earliest massive galaxies in the Universe.

[abridged] We derive photometric redshifts from 17-band optical to mid-IR photometry of 74 robust counterparts to 68 of the 126 submillimetre galaxies (SMGs) selected at 870um by LABOCA observations in the ECDFS. The median photometric redshift of identified SMGs is z=2.2\pm0.1, the interquartile range is z=1.8-2.7 and we identify 10 (~15%) high-redshift (z>3) SMGs. We derive a simple redshift estimator for SMGs based on the 3.6 and 8um fluxes, which is accurate to Delta_z~0.4 for SMGs at z<4. A statistical analysis of sources around unidentified SMGs identifies a population of likely counterparts with a redshift distribution peaking at z=2.5\pm0.3, which likely comprises ~60% of the unidentified SMGs. This confirms that the bulk of the undetected SMGs are coeval with those detected in the radio/mid-IR. We conclude that ~15% of all the SMGs are below the flux limits of our survey and lie at z>3 and hence ~30% of all SMGs have z>3. We estimate that the full S_870um>4mJy SMG population has a median redshift of 2.5\pm0.6. In contrast to previous suggestions we find no significant correlation between S_870um and redshift. The median stellar mass of the SMGs derived from SED fitting is (9.2\pm0.9)x10^10Msun and the interquartile range is (4.7-14)x10^10Msun, although we caution that uncertainty in the star-formation histories results in a factor of ~5 uncertainty in these stellar masses. The median characteristic dust temperature of SMGs is 35.9\pm1.4K and the interquartile range is 28.5-43.3K. The infrared luminosity function shows that SMGs at z=2-3 typically have higher far-IR luminosities and luminosity density than those at z=1-2. This is mirrored in the evolution of the star-formation rate density (SFRD) for SMGs which peaks at z~2. The maximum contribution of bright SMGs to the global SFRD (~5% for SMGs with S_870um>4mJy; ~50% for SMGs with S_870um>1mJy) also occurs at z~2.

We present sensitive phase-referenced VLBI results on the radio continuum emission from the z=1.87 luminous submillimeter galaxy (SMG) GOODS 850-3. The observations were carried out at 1.4 GHz using the High Sensitivity Array (HSA). Our sensitive tapered VLBI image of GOODS 850-3 at 0.47 x 0.34 arcsec (3.9 x 2.9 kpc) resolution shows a marginally resolved continuum structure with a peak flux density of 148 \pm 38 uJy/beam, and a total flux density of 168 \pm 73 uJy, consistent with previous VLA and MERLIN measurements. The derived intrinsic brightness temperature is > 5 \pm 2 x 10^3 K. The radio continuum position of this galaxy coincides with a bright and extended near-infrared source that nearly disappears in the deep HST optical image, indicating a dusty source of nearly 9 kpc in diameter. No continuum emission is detected at the full VLBI resolution (13.2 x 7.2 mas, 111 x 61 pc), with a 4-sigma point source upper limit of 26 uJy/beam, or an upper limit to the intrinsic brightness temperature of 4.7 x 10^5 K. The extent of the observed continuum source at 1.4 GHz and the derived brightness temperature limits are consistent with the radio emission (and thus presumably the far-infrared emission) being powered by a major starburst in GOODS 850-3, with a star formation rate of ~2500 M_sun/yr. Moreover, the absence of any continuum emission at the full resolution of the VLBI observations indicates the lack of a compact radio AGN source in this z=1.87 SMG.

We present the identification of a bright submm source, SMMJ163555.5+661300, detected in the lensing cluster Abell2218, for which we have accurately determined the position using observations from the Submillimeter Array (SMA). The identified optical counterpart has a spectroscopic redshift of z=4.044+-0.001 if we attribute the single emission line detected at lambda=6140AA to Lyman-alpha. This redshift identification is in good agreement with the optical/near-infrared photometric redshift as well as the submm flux ratio S_450/S_850~1.6, the radio-submm flux ratio S_1.4/S_850 < 0.004, and the 24um to 850um flux ratio S_24/S_850 < 0.005. Correcting for the gravitational lensing amplification of ~5.5, we find that the source has a far-infrared luminosity of 1.3x10^12 Lsun, which implies a star formation rate of 230 Msun/yr. This makes it the lowest-luminosity SMG known at z>4 to date. Previous CO(4-3) emission line obserations yielded a non-detection, for which we derived an upper limit of the CO line luminosity of L'_CO = 0.3x10^10 K km/s/pc^2, which is not inconsistent with the L'_CO - L_FIR relation for starburst galaxies. The best fit model to the optical and near-infrared photometry give a stellar population with an age of 1.4 Gyr and a stellar mass of 1.6x10^10 Msun. The optical morphology is compact and in the source plane the galaxy has an extent of ~6kpc x 3kpc with individual star forming knots of <500 pc in size. J163556 is not resolved in the SMA data and we place a strict upper limit on the size of the starburst region of 8kpc x 3kpc, which implies a lower limit on the star formation rate surface density of 12 Msun/yr/kpc^2. The redshift of J163556 extends the redshift distribution of faint, lensed SMGs, and we find no evidence that these have a different redshift distribution than bright SMGs.

We present a sensitive 870 micron survey of the Extended Chandra Deep Field South (ECDFS) using LABOCA on the APEX telescope. The LABOCA ECDFS Submillimetre Survey (LESS) covers the full 30' x 30' field size of the ECDFS and has a uniform noise level of 1.2 mJy/beam. LESS is thus the largest contiguous deep submillimetre survey undertaken to date. The noise properties of our map show clear evidence that we are beginning to be affected by confusion noise. We present a catalog of 126 SMGs detected with a significance level above 3.7 sigma. The ECDFS exhibits a deficit of bright SMGs relative to previously studied blank fields but not of normal star-forming galaxies that dominate the extragalactic background light (EBL). This is in line with the underdensities observed for optically defined high redshift source populations in the ECDFS (BzKs, DRGs,optically bright AGN and massive K-band selected galaxies). The differential source counts in the full field are well described by a power law with a slope of alpha=-3.2, comparable to the results from other fields. We show that the shape of the source counts is not uniform across the field. The integrated 870 micron flux densities of our source-count models account for >65% of the estimated EBL from COBE measurements. We have investigated the clustering of SMGs in the ECDFS by means of a two-point correlation function and find evidence for strong clustering on angular scales <1'. Assuming a power law dependence for the correlation function and a typical redshift distribution for the SMGs we derive a spatial correlation length of r_0=13+/-6 h^-1 Mpc.

Using the 330hr ESO-MPG 870-micron survey of the Extended Chandra Deep Field South (ECDF-S) obtained with the Large Apex BOlometer CAmera (LABOCA) on the Atacama Pathfinder EXperiment (APEX), we have carried out a stacking analysis at submillimeter (submm) wavelengths of a sample of 8266 near-infra-red (near-IR) selected (K_vega <= 20) galaxies, including 893 BzK galaxies, 1253 extremely red objects (EROs) and 737 distant red galaxies (DRGs), selected from the Multi-wavelength Survey by Yale-Chile (MUSYC). We measure average 870-micron fluxes of 0.20+-0.01mJy (20.0sigma), 0.45+-0.04mJy (11.3sigma), 0.42+-0.03mJy (14.0sigma), and 0.41+-0.04mJy (10.3sigma) for the K_vega <= 20, BzK, ERO and DRG samples, respectively. For the BzK, ERO and DRG subsamples, which overlap to some degree and are like to be at z ~ 1-2, this implies an average far-IR luminosity of ~2-6x10^{11} Lsolar and star formation rate of ~40-100Msolar. Splitting the BzK galaxies up into star-forming (sBzK) and passive (pBzK) galaxies, the former is significantly detected (0.48+-0.04mJy, 12.0sigma) while the latter is only marginally detected (0.27+-0.10mJy, 2.7sigma), thus confirming that the sBzK/pBzK criteria do isolate obscured, star forming and truly passive galaxies. The K_vega <= 20 galaxies are found to contribute with 6.6+-0.3Jy deg^{-2} (~15%) to the 870-micron extragalactic background light (EBL). sBzK and pBzK galaxies contribute 1.7+-0.2Jy deg^{-2} (~4%) and 0.2+-0.1 Jy deg^{-2} (< 0.5%) to the EBL. [Abridged]

We report on the identification of the highest redshift submm-selected source currently known: LESSJ033229.4-275619. This source was detected in the Large Apex BOlometer CAmera (LABOCA) Extended Chandra Deep Field South (ECDFS) Submillimetre Survey (LESS), a sensitive 870-um survey (~1.2-mJy rms) of the full 30'x30' ECDFS with the LABOCA camera on the Atacama Pathfinder EXperiment (APEX) telescope. The submm emission is identified with a radio counterpart for which optical spectroscopy provides a redshift of z=4.76. We show that the bolometric emission is dominated by a starburst with a star formation rate of ~1000 Msun/yr, although we also identify a moderate luminosity Active Galactic Nucleus (AGN) in this galaxy. Thus it has characteristics similar to those of z~2 submm galaxies (SMGs), with a mix of starburst and obscured AGN signatures. This demonstrates that ultraluminous starburst activity is not just restricted to the hosts of the most luminous (and hence rare) QSOs at z~5, but was also occurring in less extreme galaxies at a time when the Universe was less than 10% of its current age. Assuming that we are seeing the major phase of star formation in this galaxy, then we demonstrate that it would be identified as a luminous distant red galaxy at z~3 and that the current estimate of the space density of z>4 SMGs is only sufficient to produce ~10% of the luminous red galaxy population at these early times. However, this leaves open the possibility that some of these galaxies formed through less intense, but more extended star formation events. If the progenitors of all of the luminous red galaxies at z~3 go through an ultraluminous starburst at z>4 then the required volume density of z>4 SMGs will exceed that predicted by current galaxy formation models by more than an order of magnitude.

Context: Submillimeter galaxies (SMGs) are distant, dusty galaxies undergoing star formation at prodigious rates. Recently there has been major progress in understanding the nature of the bright SMGs (i.e. S(850um)>5mJy). The samples for the fainter SMGs are small and are currently in a phase of being built up through identification studies. Aims: We study the molecular gas content in two SMGs, SMMJ163555 and SMMJ163541, at z=1.034 and z=3.187 with unlensed submm fluxes of 0.4mJy and 6.0mJy. Both SMGs are gravitationally lensed by the foreground cluster A2218. Methods: IRAM Plateau de Bure Interferometry observations at 3mm were obtained for the lines CO(2-1) for SMMJ163555 and CO(3-2) for SMMJ163541. Additionally we obtained CO(4-3) for the candidate z=4.048 SMMJ163556 with an unlensed submm flux of 2.7mJy. Results: CO(2-1) was detected for SMMJ163555 at z=1.0313 with an integrated line intensity of 1.2+-0.2Jy km/s and a line width of 410+-120 km/s. From this a gas mass of 1.6x10^9 Msun is derived and a star formation efficiency of 440Lsun/Msun is estimated. CO(3-2) was detected for SMMJ163541 at z=3.1824, possibly with a second component at z=3.1883, with an integrated line intensity of 1.0+-0.1 Jy km/s and a line width of 280+-50 km/s. From this a gas mass of 2.2x10^10 Msun is derived and a star formation efficiency of 1000 Lsun/Msun is estimated. For SMMJ163556 the CO(4-3) is undetected within the redshift range 4.035-4.082 down to a sensitivity of 0.15 Jy km/s. Conclusions: Our CO line observations confirm the optical redshifts for SMMJ163555 and SMMJ163541. The CO line luminosity L'_CO for both galaxies is consistent with the L_FIR-L'_CO relation. SMMJ163555 has the lowest FIR luminosity of all SMGs with a known redshift and is one of the few high redshift LIRGs whose properties can be estimated prior to ALMA.

We have used the IRAM Plateau de Bure mm interferometer and the UKIRT 1-5 um Imager Spectrometer to test the connection between the major phases of spheroid growth and nuclear accretion by mapping CO emission in nine submm-detected QSOs at z=1.7-2.6 with black hole (BH) masses derived from near-infrared spectroscopy. With a previously published QSO, we present sensitive CO(3-2) observations of 10 submm-detected QSOs selected at the epoch of peak activity in both QSOs and submm galaxies (SMGs). CO is detected in 5/6 very optically luminous (M_B~-28) submm-detected QSOs with BH masses M_BH~10^9-10^10 Msun, confirming the presence of large gas reservoirs of M_gas~3.4x10^10 Msun. However, we find that their BH masses are ~30 times too large and their surface density is ~300 times too small to be related to typical SMGs in an evolutionary sequence. Conversely, we measure weaker CO emission in four fainter (M_B~-25) submm-detected QSOs with properties, BH masses (M_BH~5x10^8 Msun), and surface densities similar to SMGs. These QSOs appear to lie near the local M_BH/M_sph relation, making them plausible `transition objects' in the proposed evolutionary sequence linking QSOs to the formation of massive young galaxies and BHs at high-redshift. We show that SMGs have a higher incidence of bimodal CO line profiles than seen in our QSO sample, which we interpret as an effect of their relative inclinations, with the QSOs seen more face-on. Finally, we find that the gas masses of the four fainter submm-detected QSOs imply that their star formation episodes could be sustained for ~10 Myr, and are consistent with representing a phase in the formation of massive galaxies which overlaps a preceding SMG starburst phase, before subsequently evolving into a population of present-day massive ellipticals. [abridged]

We have conducted a submillimetre mapping survey of faint, gravitationally lensed sources, where we have targeted twelve galaxy clusters and additionally the NTT Deep Field. The total area surveyed is 71.5 arcmin^2 in the image plane; correcting for gravitational lensing, the total area surveyed is 40 arcmin^2 in the source plane for a typical source redshift z=2.5. In the deepest maps, an image plane depth of 1sigma r.m.s. ~0.8 mJy is reached. This survey is the largest survey to date to reach such depths. In total 59 sources were detected, including three multiply-imaged sources. The gravitational lensing makes it possible to detect sources with flux density below the blank field confusion limit. The lensing corrected fluxes ranges from 0.11 mJy to 19 mJy. After correcting for multiplicity there are 10 sources with fluxes <2 mJy of which 7 have sub-mJy fluxes, doubling the number of such sources known. Number counts are determined below the confusion limit. At 1 mJy the integrated number count is ~10^4 deg^-2, and at 0.5 mJy it is ~2x10^4 deg^-2. Based on the number counts, at a source plan flux limit of 0.1 mJy, essentially all of the 850um background emission has been resolved. The dominant contribution (> 50 per cent) to the integrated background arises from sources with fluxes S(850) between 0.4 and 2.5 mJy, while the bright sources S(850) > 6mJy contribute only 10 per cent.

We present a reconstruction of the mass distribution of galaxy cluster Abell 1689 at z = 0.18 using detected strong lensing features from deep HST/ACS observations and extensive ground based spectroscopy. Earlier analyses have reported up to 32 multiply imaged systems in this cluster, of which only 3 were spectroscopically confirmed. In this work, we present a parametric strong lensing mass reconstruction using 24 multiply imaged systems with newly determined spectroscopic redshifts, which is a major step forward in building a robust mass model. In turn, the new spectroscopic data allows a more secure identification of multiply imaged systems. The resultant mass model enables us to reliably predict the redshifts of additional multiply imaged systems for which no spectra are currently available, and to use the location of these systems to further constrain the mass model. In particular, we have detected 5 strong galaxy-galaxy lensing systems just outside the Einstein ring region, further constraining the mass profile. Our strong lensing mass model is consistent with that inferred from our large scale weak lensing analysis derived using CFH12k wide field images. Thanks to a new method for reliably selecting a well defined background lensed galaxy population, we resolve the discrepancy found between the strong and weak lensing mass models reported in earlier work. [ABRIDGED]

We present Spitzer/IRS mid-infrared spectra for 15 gravitationally lensed, 24 micron--selected galaxies, and combine the results with 4 additional very faint galaxies with IRS spectra in the literature. The median intrinsic 24 micron flux density of the sample is 130 microJy, enabling a systematic survey of the spectral properties of the very faint 24 micron sources that dominate the number counts of Spitzer cosmological surveys. Six of the 19 galaxy spectra (32%) show the strong mid-IR continuua expected of AGN; X-ray detections confirm the presence of AGN in three of these cases, and reveal AGNs in two other galaxies. These results suggest that nuclear accretion may contribute more flux to faint 24 micron--selected samples than previously assumed. Almost all the spectra show some aromatic (PAH) emission features; the measured aromatic flux ratios do not show evolution from z=0. In particular, the high S/N mid-IR spectrum of SMM J163554.2+661225 agrees remarkably well with low--redshift, lower--luminosity templates. We compare the rest-frame 8 micron and total infrared luminosities of star--forming galaxies, and find that the behavior of this ratio with total IR luminosity has evolved modestly from z=2 to z=0. Since the high aromatic--to--continuum flux ratios in these galaxies rule out a dominant contribution by AGN, this finding implies systematic evolution in the structure and/or metallicity of infrared sources with redshift. It also has implications for the estimates of star forming rates inferred from 24 micron measurements, in the sense that at z ~2, a given observed frame 24 micron luminosity corresponds to a lower bolometric luminosity than would be inferred from low-redshift templates of similar luminosity at the corresponding rest wavelength.

We have imaged the central ~1kpc of the circumnuclear starburst disk in the galaxy NGC253 in the HCN(1-0), HCO+(1-0), and CO(1-0) transitions at 60pc resolution using the Owens Valley Radio Observatory Millimeter-Wavelength Array (OVRO). We have also obtained Atacama Pathfinder Experiment (APEX) observations of the HCN(4-3) and the HCO+(4-3) lines of the starburst disk. We find that the emission from the HCN(1-0) and HCO+(1-0) transitions, both indicators of dense molecular gas, trace regions which are non-distinguishable within the uncertainties of our observations. Even though the continuum flux varies by more than a factor 10 across the starburst disk, the HCN/HCO+ ratio is constant throughout the disk, and we derive an average ratio of 1.1+/-0.2. From an excitation analysis we find that all lines from both molecules are subthermally excited and that they are optically thick. This subthermal excitation implies that the observed HCN/HCO+ line ratio is sensitive to the underlying chemistry. The constant line ratio thus implies that there are no strong abundance gradients across the starburst disk of NGC253. This finding may also explain the variations in L'(HCN)/L'(HCO+) between different star forming galaxies both nearby and at high redshifts.

We present CO(1-0) observations of the high-redshift quasi-stellar objects (QSOs) BR 1202-0725 (z=4.69), PSS J2322+1944 (z=4.12), and APM 08279+5255 (z=3.91) using the NRAO Green Bank Telescope (GBT) and the MPIfR Effelsberg 100m telescope. We detect, for the first time, the CO ground-level transition in BR 1202-0725. For PSS J2322+1944 and APM 08279+5255, our observations result in line fluxes that are consistent with previous NRAO Very Large Array (VLA) observations, but they reveal the full line profiles. We report a typical lensing-corrected velocity-integrated intrinsic CO(1-0) line luminosity of L'(CO) = 5 x 10^10 K km/s pc^2 and a typical total H_2 mass of M(H2) = 4 x 10^10 M_sun for the sources in our sample. The CO/FIR luminosity ratios of these high-z sources follow the same trend as seen for low-z galaxies, leading to a combined solution of log(L_FIR) = (1.39 +/- 0.05) x log(L(CO))-1.76. It has previously been suggested that the molecular gas reservoirs in some quasar host galaxies may exhibit luminous, extended CO(1-0) components that are not observed in the higher-J CO transitions. Utilizing the line profiles and the total intensities of our observations and large velocity gradient (LVG) models based on previous results for higher-J CO transitions, we derive that emission from all CO transitions is described well by a single gas component where all molecular gas is concentrated in a compact nuclear region. Thus, our observations and models show no indication of a luminous extended, low surface brightness molecular gas component in any of the high-redshift QSOs in our sample. If such extended components exist, their contribution to the overall luminosity is limited to at most 30%.

We present a search for 183 GHz H_2O(3_13-2_20) emission in the infrared-luminous quasar MG 0751+2716 with the NRAO Very Large Array (VLA). At z=3.200+/-0.001, this water emission feature is redshifted to 43.6 GHz. As opposed to the faint rotational transitions of HCN (the standard high-density tracer at high-z), H_2O(3_13-2_20) is observed with high maser amplification factors in Galactic star-forming regions. It therefore holds the potential to trace high-density star-forming regions in the distant universe. If indeed all star-forming regions in massively star-forming galaxies at z>3 have similar physical properties as e.g. the Orion or W49N molecular cloud cores, the flux ratio between the maser-amplified H_2O(3_13-2_20) and the thermally excited CO(1-0) transitions may be as high as factor of 20 (but has to be corrected by their relative filling factor). MG 0751+2716 is a strong CO(4-3) emitter, and therefore one of the most suitable targets to search for H_2O(3_13-2_20) at cosmological redshifts. Our search resulted in an upper limit in line luminosity of L'(H_2O) < 0.6 x 10^9 K km/s pc^2. Assuming a brightness temperature of T_b(H_2O) ~= 500 K for the maser emission and CO properties from the literature, this translates to a H_2O(3_13-2_20)/CO(4-3) area filling factor of less than 1%. However, this limit is not valid if the H_2O(3_13-2_20) maser emission is quenched, i.e. if the line is only thermally excited. We conclude that, if our results were to hold for other high-z sources, H_2O does not appear to be a more luminous alternative to HCN to detect high-density gas in star-forming environments at high redshift.

We report the detection of HCO+(1-0) emission towards the Cloverleaf quasar (z=2.56) through observations with the Very Large Array. This is the first detection of ionized molecular gas emission at high redshift (z>2). HCO+ emission is a star formation indicator similar to HCN, tracing dense molecular hydrogen gas (n(H_2) ~= 10^5 cm^{-3}) within star-forming molecular clouds. We derive a lensing-corrected HCO+ line luminosity of L'(HCO+) = 3.5 x 10^9 K km/s pc^2. Combining our new results with CO and HCN measurements from the literature, we find a HCO+/CO luminosity ratio of 0.08 and a HCO+/HCN luminosity ratio of 0.8. These ratios fall within the scatter of the same relationships found for low-z star-forming galaxies. However, a HCO+/HCN luminosity ratio close to unity would not be expected for the Cloverleaf if the recently suggested relation between this ratio and the far-infrared luminosity were to hold. We conclude that a ratio between HCO+ and HCN luminosity close to 1 is likely due to the fact that the emission from both lines is optically thick and thermalized and emerges from dense regions of similar volumes. The CO, HCN and HCO+ luminosities suggest that the Cloverleaf is a composite AGN--starburst system, in agreement with the previous finding that about 20% of the total infrared luminosity in this system results from dust heated by star formation rather than heating by the AGN. We conclude that HCO+ is potentially a good tracer for dense molecular gas at high redshift.

We present the identification of a bright submillimeter galaxy (SMG) in the field of Abell 2218. The galaxy has a spectroscopic redshift of ~4, and is currently the highest redshift SMG known. It is detected at all wavelengths from optical to submm, including the Spitzer IRAC bands. We discuss the properties of this galaxy, which is undergoing intense star formation at a rate~600 Msun/yr. We also compare the properties to those of radio-preselected submm-bright galaxies. The z~4 result extends the redshift distribution of SMGs.

Extracting sources with low signal-to-noise from maps with structured background is a non-trivial task which has become important in studying the faint end of the submillimetre number counts. In this article we study source extraction from submillimetre jiggle-maps from the Submillimetre Common-User Bolometer Array (SCUBA) using the Mexican Hat Wavelet (MHW), an isotropic wavelet technique. As a case study we use a large (11.8 arcmin^2) jiggle-map of the galaxy cluster Abell 2218, with a 850um 1sigma r.m.s. sensitivity of 0.6-1mJy. We show via simulations that MHW is a powerful tool for reliable extraction of low signal-to-noise sources from SCUBA jiggle-maps and nine sources are detected in the A2218 850um image. Three of these sources are identified as images of a single background source with an unlensed flux of 0.8mJy. Further, two single-imaged sources also have unlensed fluxes <2mJy, below the blank-field confusion limit. In this ultradeep map, the individual sources detected resolve nearly all of the extragalactic background light at 850um, and the deep data allow to put an upper limit of 44 sources per arcmin^2 to 0.2mJy at 850um.

We present a study of the submillimeter (submm) emission of Distant Red Galaxies (DRGs). The DRGs are selected by the criterion (J-K)>2.3, and are generally massive galaxies at redshifts higher than 2, with red rest-frame optical colors. Using a deep SCUBA submm image of a field centred on the cluster MS1054-03, we obtain a statistical detection of the DRGs at redshift z=2-3.5, with an average 850um flux density of 1.11+-0.28mJy. The detection implies an average star formation rate (SFR) of 127+-34 M_sun/yr (lensing corrected), assuming that the far-infrared (FIR) spectral energy distribution (SED) is well-described by a modified blackbody. The SFR derived from the submm agrees well with SFRs derived from SED fitting of optical-near-infrared data and average X-ray emission. Constant Star Formation models imply ages of 2Gyr, extinction A_V=2.4 mag, which is consistent with the FIR to rest-frame optical luminosity ratio of ~15. DRGs are older and have lower SFRs relative to optical luminosity than (ultra-)luminous infrared galaxies, although their FIR luminosities are similar. The DRGs at 2<z<3.5 and the Extremely Red Objects ((I-K)>4) at 1<z<2, which were also investigated, contribute 5.7 and 5.9 Jy deg^-2 respectively to the submm background. Simple estimates suggest that these populations contribute ~50% of the flux from sources with 0.5 < f_850 < 5 mJy, which is where the peak of energy is produced. We have therefore uncovered one of the most important populations of galaxies contributing to the sub-mm background.

We present the results of observations from the IRAM array of the submm galaxy SMMJ16359+6612 lying at z=2.516 behind the massive cluster A2218. The foreground gravitational lens produces 3 images with a total mag. of 45 of this faint submm galaxy, which has an intrinsic submm flux of f_850mic=0.8mJy placing it below the confusion limit of blank-field surveys. The substantial magnification provides a rare opportunity to probe the nature of a distant sub-mJy submm-selected galaxy, part of the population which produces the bulk of the submm cosmic far-infrared background. Our observations detect the CO(3-2) line in all 3 images, as well as the CO(7-6) line and the dust continuum at 1.3mm for the brightest image. The CO(3-2) velocity profile displays a double-peak profile which is well fit by two Gaussians with FWHM of 220km/s and separated by 280km/s. We estimate the dynamical mass of the system to be ~1.5 10^10 M_sun and an H2 gas mass of 2.6 10^9 M_sun. We identify a spatial offset of ~1'' between the two CO(3-2) velocity components, modeling of which indicates that the offset corresponds to just ~3kpc in projection at z=2.5. The spatial and velocity properties of these two components are closely related to features detected in previously published Halpha spectroscopy. We conclude that this source is likely to be a compact merger of 2 fairly typical Ly-break galaxies with a maximal separation between the two nuclei of ~3kpc. This system is much less luminous and massive than other high-z submm galaxies studied to date, but it bears a close similarity to similarly luminous, dusty starburst resulting from lower-mass mergers in the local Universe.

We report the detection of discrete, lensed radio emission from the multiply imaged, z=2.516 submillimetre selected galaxy, SMM J16359+6612. All three images are detected in deep WSRT 1.4 GHz and VLA 8.2 GHz observations, and the radio positions are coincident with previous sub-mm SCUBA observations of this system. This is the widest separation lens system to be detected in the radio so far, and the first time that multiply imaged lensed radio emission has been detected from a star forming galaxy -- all previous multiply-lensed radio systems being associated with radio-loud AGN. Taking into account the total magnification of ~45, the WSRT 1.4 GHz observations suggest a star formation rate of 500 Solar mass/yr. The source has a steep radio spectrum (alpha -0.7) and an intrinsic flux density of just 3 microJy at 8.2 GHz. Three other SCUBA sources in the field are also detected by the WSRT, including SMMJ16359+66118, a singly imaged (and magnified) arclet at z=1.034. Higher resolution radio observations of SMMJ16359+6612 (and other highly magnified star forming galaxies) provide a unique opportunity to study the general properties and radio morphology of intrinsically faint, distant and obscured star forming galaxies. They can also help to constrain the technical specification of next generation radio telescopes, such as the Square Kilometre Array.

We report the detection of CO ($J$=3$\to$2) line emission from all three multiple images (A,B and C) of the intrinsically faint ($\simeq$ 0.8 mJy) submillimeter-selected galaxy SMM J16359+6612. The brightest source of the submm continuum emission (B) also corresponds to the brightest CO emission, which is centered at $z$=2.5168, consistent with the pre-existing redshift derived from \Ha. The observed CO flux in the A, B and C images is 1.2, 3.5 and 1.6 Jy \kms respectively, with a linewidth of $500\pm 100$ \kms. After correcting for the lensing amplification, the CO flux corresponds to a molecular gas mass of $\sim 2\times 10 ^{10} h_{71}^{-2}$ \Msun, while the extent of the CO emission indicates that the dynamical mass of the system $\sim9\times10^{10}$ \Msun. Two velocity components are seen in the CO spectra; these could arise from either a rotating compact ring or disk of gas, or merging substructure. The star formation rate in this galaxy was previously derived to be $\sim$100--500 \Msun \yr. If all the CO emission arises from the inner few kpc of the galaxy and the galactic CO-to-H$_2$ conversion factor holds, then the gas consumption timescale is a relatively short 40 Myr, and so the submm emission from SMM J16359+6612 may be produced by a powerful, but short-lived circumnuclear starburst event in an otherwise normal and representative high-redshift galaxy.

We recently identified a substantial population of galaxies at z>2 with red rest-frame optical colors. These distant red galaxies (DRGs) are efficiently selected by the simple observed color criterion J-K>2.3. In this paper we present NIR spectroscopy with Keck/NIRSPEC of six DRGs at 2.4<z<3.2. We detect continuum emission and emission lines of all observed galaxies. Equivalent widths of H alpha are 20-30 Ang, smaller than measured for LBGs and nearby LIRGs, and comparable to normal nearby galaxies. The modest equivalent widths imply that the galaxies either have a decreasing star formation rate, or that they are very dusty. Fitting both the photometry and the H alpha lines, we find continuum extinction A_V=1-2 mag, ages 1-2.5 Gyr, star formation rates 200-400 solar masses/yr, and stellar masses 1-5x10^11 solar masses for models with constant star formation rates. From [NII]/H alpha ratios we infer that the metallicities are high, 1-1.5 x Solar. For four galaxies we can determine line widths from the optical emission lines. The widths are high, ranging from 130-240 km/s, and by combining data for LBGs and DRGs we find significant correlations between linewidth and restframe U-V color, and between linewidth and stellar mass. The latter correlation has a similar slope and offset as the ``baryonic Tully-Fisher relation'' for nearby galaxies. The median dynamical mass is ~2x10^11 solar masses, supporting the high stellar masses inferred from the photometry. We find that the median M/L_V ~ 0.8, a factor of ~5 higher than measured for LBGs. We infer from our small sample that DRGs are dustier, more metal rich, more massive, and have higher ages than z=3 LBGs of the same rest-frame V-band luminosity. Their high M/L ratios imply that they contribute significantly to the stellar mass density at z~2.5. [ABRIDGED]

We present observations of a remarkable submillimetre-selected galaxy, SMMJ16359+6612. This distant galaxy lies behind the core of a massive cluster of galaxies, A2218, and is gravitationally lensed by the foreground cluster into three discrete images which were identified in deep submillimetre maps of the cluster core at both 450 and 850micron. Subsequent follow-up using deep optical and NIR images identify a faint counterpart to each of the 3 images, with similar red optical--NIR colours and HST morphologies. By exploiting a detailed mass model for the cluster lens we estimate that the combined images of this galaxy are magnified by a factor of ~45, implying that this galaxy would have un-lensed magnitudes K_s=22.9 and I=26.1, and an un-lensed 850micron flux density of only 0.8mJy. Moreover, the highly constrained lens model predicted the redshift of SMMJ16359+6612 to be z=2.6+/-0.4. We confirm this estimate using deep optical and NIR Keck spectroscopy, measuring a redshift of z=2.516. SMMJ16359+6612 is the faintest submm-selected galaxy so far identified with a precise redshift. Thanks to the large gravitational magnification of this source, we identify 3 sub-components in this submm galaxy, which are also seen in the NIRSPEC data, arguing for either a strong dust (lane) absorption or a merger. Interestingly, there are 2 other highly-amplified galaxies at almost identical redshifts in this field (although neither is a strong submm emitter). The 3 galaxies lie within a ~100kpc region on the background sky, suggesting this submm galaxy is located in a dense high-redshift group.

We have conducted a new deep SCUBA survey, which has targetted 12 lensing galaxy clusters and one blank field. In this survey we have detected several sub-mJy sources after correcting for the gravitational lensing by the intervening clusters. We here present the preliminary results and point out two highlights.

We present the discovery of a dusty quasar, SMMJ04135+10277, located behind the galaxy cluster Abell 478. The quasar was discovered as the brightest source in a submillimetre survey of high redshift galaxies lensed by foreground rich clusters of galaxies in a project aimed at studying the cosmic star formation history of dusty galaxies. With submm fluxes of S_850 = 25+-2.8 mJy and S_450 = 55+-17 mJy this object is one of the brightest submillimetre sources known. Optical observations revealed a quasar with redshift z=2.837. The quasar was also detected at shorter infrared wavelengths with the Infrared Space Observatory. This object is the first quasar discovered by its submm emission. Given the general lack of overlap between deep submm and X-ray samples, usually interpreted as a low incidence of active nuclei in submm samples, this is an unusual object. Analysis of number counts of quasars and of submm galaxies bears out this suggestion. We compare the properties of SMMJ04135+10277 to those of optically selected quasars with submm emission, and argue that the optical faintness results from a large viewing angle with the direction of relativistic beaming, and not from abnormally high extinction. We also find indications that the bulk of the submm flux density is not powered by the quasar nucleus. This conclusion is supported by analysis of the infrared spectral energy distribution. These results are consistent with previous observations that quasars at higher redshift tend to have a more prominent cold dust component, most likely powered by extended star formation in the host galaxy. The quasar is found to have a total infrared luminosity of (2.9+-0.5)x10^13 L_sun, dominated by the emission from cool dust. (Abridged)

Can the submillimetre counts and background be produced by applying a locally derived extinction correction to the population of Lyman break galaxies? We investigate the submillimetre emission of two strongly lensed Lyman break galaxies (MS1512+36-cB58 and MS1358+62-G1) and find that the procedure that is used to predict the submillimetre emission of the Lyman break galaxy population overpredicts the observed 850 micron fluxes by a factor of up to 14. This result calls for caution in applying local correlations to distant galaxies. It also shows that large extinction corrections on Lyman break galaxies should be viewed with skepticism. It is concluded that the Lyman break galaxies may contribute to the submillimetre background at the 25 to 50% level. The brighter submillimetre galaxies making up the rest of the background are either not detected in optical surveys, or if they are detected, their submillimetre emission cannot be reliably estimated from their rest-frame ultraviolet properties.

Can the submillimeter counts and background be produced by applying a locally derived extinction correction to the population of Lyman break galaxies? We investigate the submillimeter emission of two strongly lensed Lyman break galaxies (MS1512+36-cB58 and MS1358+62-G1) and find that the procedure that is used to predict the submillimeter emission of the Lyman break galaxy population overpredicts the observed 850micron fluxes by up to a factor of 14. This result calls for caution in applying local correlations to distant galaxies. It also shows that large extinction corrections on Lyman break galaxies should be viewed with skepticism. It is concluded that the Lyman break galaxies may contribute to the submillimeter background at the 25 to 50% level. The brighter submillimeter galaxies making up the rest of the background are either not detected in optical surveys, or if they are detected, their submillimeter emission cannot be reliably estimated from their rest-frame ultraviolet properties.

We report the discovery of a z=2.83 quasar in the field of the cooling flow galaxy cluster Abell 478. This quasar was first detected in a submm survey of star forming galaxies at high redshifts, as the brightest source. We discuss the optical spectrum and far-IR spectral energy distribution (SED) of this object.