The South Pole Telescope (SPT) has systematically identified 81 high-redshift, strongly gravitationally lensed, dusty star-forming galaxies (DSFGs) in a 2500 square degree cosmological mm-wave survey. We present the final spectroscopic redshift survey of this flux-limited ($S_{870 \, \mathrm{\mu m}} > 25 \, \mathrm{mJy}$) sample. The redshift survey was conducted with the Atacama Large Millimeter/submillimeter Array across the 3 mm spectral window, targeting carbon monoxide line emission. The SPT sample is now spectroscopically complete, with redshifts spanning the range $1.9$$<$$z$$<$$6.9$, with a median of $z=3.9 \pm 0.2$. We present the mm through far-infrared photometry and spectral energy density fits for all sources, along with their inferred intrinsic properties. Comparing the properties of the SPT sources to the unlensed DSFG population, we demonstrate that the SPT-selected DSFGs represent the most extreme infrared-luminous galaxies, even after accounting for strong gravitational lensing. The SPT sources have a median star formation rate of $2.3(2)\times 10^3 \, \mathrm{M_\odot yr^{-1}}$ and a median dust mass of $1.4(1)\times10^9\, \mathrm{M_\odot}$. However, the inferred gas depletion timescales of the SPT sources are comparable to those of unlensed DSFGs, once redshift is taken into account. This SPT sample contains roughly half of the known spectroscopically confirmed DSFGs at $z$$>$$5$, making this the largest sample of high-redshift DSFGs to-date, and enables us to measure the "high-redshift tail" of the distribution of luminous DSFGs. Though galaxy formation models struggle to account for the SPT redshift distribution, the larger sample statistics from this complete and well-defined survey will help inform future theoretical efforts.

We present an extensive ALMA spectroscopic follow-up programme of the $z\,{=}\,4.3$ structure SPT2349$-$56, one of the most actively star-forming proto-cluster cores known, to identify additional members using their [C{\sc ii}] 158\,$\mu$m and \mbox{CO(4--3)} lines. In addition to robustly detecting the 14 previously published galaxies in this structure, we identify a further 15 associated galaxies at $z\,{=}\,4.3$, resolving 55$\,{\pm}\,$5\,per cent of the 870-$\mu$m flux density at 0.5\,arcsec resolution compared to 21\,arcsec single-dish data. These galaxies are distributed into a central core containing 23 galaxies extending out to 300\,kpc in diameter, and a northern extension, offset from the core by 400\,kpc, containing three galaxies. We discovered three additional galaxies in a red {\it Herschel\/}-SPIRE source 1.5\,Mpc from the main structure, suggesting the existence of many other sources at the same redshift as SPT2349$-$56 that are not yet detected in the limited coverage of our data. An analysis of the velocity distribution of the central galaxies indicates that this region may be virialized with a mass of (9$\pm$5)$\,{\times}\,$10$^{12}$\,M$_{\odot}$, while the two offset galaxy groups are about 30 and 60\,per cent less massive and show significant velocity offsets from the central group. We calculate the [C{\sc ii}] and far-infrared number counts, and find evidence for a break in the [C{\sc ii}] luminosity function. We estimate the average SFR density within the region of SPT2349$-$56 containing single-dish emission (a proper diametre of 720\,kpc), assuming spherical symmetry, to be roughly 4$\,{\times}\,10^4$\,M$_{\odot}$\,yr$^{-1}$\,Mpc$^{-3}$; this may be an order of magnitude greater than the most extreme examples seen in simulations.

We analyze an extremely deep 450-$\mu$m image ($1\sigma=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$\mu$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$\mu$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$\mu$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$\mu$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$\mu$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$\mu$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies.

We present the first predictions for the $L_{\rm [CII]}$ - SFR relation and [CII] luminosity function (LF) in the Epoch of Reionization (EoR) based on cosmological hydrodynamics simulations using the SIMBA suite plus radiative transfer calculations via S\'IGAME. The sample consists of 11,137 galaxies covering halo mass $\log M_{\rm halo}\in$[9, 12.4] $M_\odot$, star formation rate SFR$\in$[0.01, 330] $M_\odot$ yr$^{-1}$, and metallicity $<Z_{\rm gas}>_{\rm SFR}\in$[0.1, 0.9] $Z_\odot$. The simulated $L_{\rm [CII]}$-SFR relation is consistent with the range observed, but with a spread of $\simeq$0.3 dex at the high end of SFR ($>$100 $M_\odot$ yr$^{-1}$) and $\simeq$0.6 dex at the lower end, and there is tension between our predictions and the values of $L_{\rm [CII]}$ above 10$^{8.5}$ $L_\odot$ observed in some galaxies reported in the literature. The scatter in the $L_{\rm [CII]}$-SFR relation is mostly driven by galaxy properties, such that at a given SFR, galaxies with higher molecular gas mass and metallicity have higher $L_{\rm [CII]}$. The [CII] LF predicted by SIMBA is consistent with the upper limits placed by the only existing untargeted flux-limited [CII] survey at the EoR (ASPECS) and those predicted by semi-analytic models. We compare our results with existing models and discuss differences responsible for the discrepant slopes in the $L_{\rm [CII]}$-SFR relatiion.

We present Atacama Compact Array and Atacama Pathfinder Experiment observations of the [N II] 205 $\mu$m fine-structure line in 40 sub-millimetre galaxies lying at redshifts z = 3 to 6, drawn from the 2500 deg$^2$ South Pole Telescope survey. This represents the largest uniformly selected sample of high-redshift [N II] 205 $\mu$m measurements to date. 29 sources also have [C II] 158 $\mu$m line observations allowing a characterization of the distribution of the [C II] to [N II] luminosity ratio for the first time at high-redshift. The sample exhibits a median L$_{[C II]}$ /L$_{[N II]}$ $\approx$ 11 and interquartile range of 5.0 to 24.7. These ratios are similar to those observed in local (U)LIRGs, possibly indicating similarities in their interstellar medium. At the extremes, we find individual sub-millimetre galaxies with L$_{[C II]}$ /L$_{[N II]}$ low enough to suggest a smaller contribution from neutral gas than ionized gas to the [C II] flux and high enough to suggest strongly photon or X-ray region dominated flux. These results highlight a large range in this line luminosity ratio for sub-millimetre galaxies, which may be caused by variations in gas density, the relative abundances of carbon and nitrogen, ionization parameter, metallicity, and a variation in the fractional abundance of ionized and neutral interstellar medium.

To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission features tend to be centrally concentrated, which is in contrast to the shallower distribution of CO 1-0 and the stellar component. The dense-gas fraction ($f_\text{dense}$, traced by the velocity-integrated-intensity ratios of HCN/CO and HCO$^+$/CO) and the ratio $R_\text{31}$ (CO 3-2/1-0) decline towards larger galactocentric distances, but increase with higher SFR surface density. The radial variation and the large scatter of $f_\text{dense}$ and $R_\text{31}$ imply distinct physical conditions in different regions of the galactic disc. The relationships of $f_\text{dense}$ versus $\Sigma_\text{stellar}$, and SFE$_\text{dense}$ versus $\Sigma_\text{stellar}$ are explored. SFE$_\text{dense}$ increases with higher $\Sigma_\text{stellar}$ in this galaxy, which is inconsistent with previous work that used HCN 1-0 data. This implies that existing stellar components might have different effects on the high-$J$ HCN and HCO$^+$ than their low-$J$ emission. We also find that SFE$_\text{dense}$ seems to be decreasing with higher $f_\text{dense}$, which is consistent with previous works, and it suggests that the ability of the dense gas to form stars diminishes when the average density of the gas increases. This is expected in a scenario where only the regions with high-density contrast collapse and form stars.

We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K ($\beta=1.9$) derived under the common assumption of an optically thin far-infrared dust emission, but fully consistent with $T_{\rm d}=52\pm5$ K of a general opacity model where the optical depth ($\tau$) reaches unity at a wavelength of $\lambda_0=170\pm23$ $\mu$m. Moreover, the general opacity solution returns a factor of $\sim 2\times$ lower dust mass and, hence, a lower molecular gas mass for a fixed gas-to-dust ratio, than with the optically thin dust model. The derived properties of GN20 thus provide an appealing solution to the puzzling discovery of starbursts appearing colder than main-sequence galaxies above $z>2.5$, in addition to a lower dust-to-stellar mass ratio that approaches the physical value predicted for starburst galaxies.

We study the r31=L'CO(3-2)/L'CO(1-0) luminosity line ratio in a sample of nearby (z < 0.05) galaxies: 25 star-forming galaxies (SFGs) from the xCOLD GASS survey, 36 hard X-ray selected AGN host galaxies from BASS and 37 infrared luminous galaxies from SLUGS. We find a trend for r31 to increase with star-formation efficiency (SFE). We model r31 using the UCL-PDR code and find that the gas density is the main parameter responsible for variation of r31, while the interstellar radiation field and cosmic ray ionization rate play only a minor role. We interpret these results to indicate a relation between SFE and gas density. We do not find a difference in the r31 value of SFGs and AGN host galaxies, when the galaxies are matched in SSFR (<r31>= 0.52 +/- 0.04 for SFGs and <r31> = 0.53 +/- 0.06 for AGN hosts). According to the results of UCL-PDR models, the X-rays can contribute to the enhancement of the CO line ratio, but only for strong X-ray fluxes and for high gas density (nH > 10$^4$ cm-3). We find a mild tightening of the Kennicutt-Schmidt relation when we use the molecular gas mass surface density traced by CO(3-2) (Pearson correlation coefficient R=0.83), instead of the molecular gas mass surface density traced by CO(1-0) (R=0.78), but the increase in correlation is not statistically significant (p-value=0.06). This suggests that the CO(3-2) line can be reliably used to study the relation between SFR and molecular gas for normal SFGs at high redshift, and to compare it with studies of low-redshift galaxies, as is common practice.

Water ($\rm H_{2}O$), one of the most ubiquitous molecules in the universe, has bright millimeter-wave emission lines easily observed at high-redshift with the current generation of instruments. The low excitation transition of $\rm H_{2}O$, p$-$$\rm H_{2}O$(202 $-$ 111) ($\nu_{rest}$ = 987.927 GHz) is known to trace the far-infrared (FIR) radiation field independent of the presence of active galactic nuclei (AGN) over many orders-of-magnitude in FIR luminosity (L$_{\rm FIR}$). This indicates that this transition arises mainly due to star formation. In this paper, we present spatially ($\sim$0.5 arcsec corresponding to $\sim$1 kiloparsec) and spectrally resolved ($\sim$100 kms$^{-1}$) observations of p$-$$\rm H_{2}O$(202 $-$ 111) in a sample of four strong gravitationally lensed high-redshift galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA). In addition to increasing the sample of luminous ($ > $ $10^{12}$L$_{\odot}$) galaxies observed with $\rm H_{2}O$, this paper examines the L$_{\rm H_{2}O}$/L$_{\rm FIR}$ relation on resolved scales for the first time at high-redshift. We find that L$_{\rm H_{2}O}$ is correlated with L$_{\rm FIR}$ on both global and resolved kiloparsec scales within the galaxy in starbursts and AGN with average L$_{\rm H_{2}O}$/L$_{\rm FIR}$ =$2.76^{+2.15}_{-1.21}\times10^{-5}$. We find that the scatter in the observed L$_{\rm H_{2}O}$/L$_{\rm FIR}$ relation does not obviously correlate with the effective temperature of the dust spectral energy distribution (SED) or the molecular gas surface density. This is a first step in developing p$-$$\rm H_{2}O$(202 $-$ 111) as a resolved star formation rate (SFR) calibrator.

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.

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.

Understanding the nuclear growth and feedback processes in galaxies requires investigating their often obscured central regions. One way to do this is to use (sub)millimeter line emission from vibrationally excited HCN (HCN-vib). It has been suggested that the most intense HCN-vib emission from a galaxy is connected to a phase of nuclear growth that occurs before the nuclear feedback processes have been fully developed. We aim to investigate if there is a connection between the presence of strong HCN-vib emission and the development of feedback in (U)LIRGs. We collected literature and archival data to compare the luminosities of rotational lines of HCN-vib, normalized to the total infrared luminosity, to the median velocities of 119 {\mu}m OH absorption lines, potentially indicating outflows, in a total of 17 (U)LIRGs. The most HCN-vib luminous systems all lack signatures of significant molecular outflows in the far-infrared OH absorption lines. However, at least some of the systems with bright HCN-vib emission do have fast and collimated outflows that can be seen in spectral lines at longer wavelengths. We conclude that the galaxy nuclei with the highest L(HCN-vib)/L(IR) do not drive wide-angle outflows detectable using the median velocities of far-infrared OH absorption lines. It is possible that this is due to an orientation effect where sources which are oriented in such a way that their outflows are not along our line of sight also radiate a smaller proportion of their infrared luminosity in our direction. It could also be that massive wide-angle outflows destroy the deeply embedded regions responsible for bright HCN-vib emission, so that the two phenomena cannot coexist. This would strengthen the idea that vibrationally excited HCN traces a heavily obscured stage of evolution before nuclear feedback mechanisms are fully developed.

SPT0346-52 is one of the most most luminous and intensely star-forming galaxies in the universe, with L_FIR > 10^13 L_sol and Sigma_SFR ~ 4200 M_sol yr^-1 kpc^-2. In this paper, we present ~0.15'' ALMA observations of the [CII]158micron emission line in this z=5.7 dusty star-forming galaxy. We use a pixellated lensing reconstruction code to spatially and kinematically resolve the source-plane [CII] and rest-frame 158 micron dust continuum structure at ~700 pc (~0.12'') resolution. We discuss the [CII] deficit with a pixellated study of the L_[CII]/L_FIR ratio in the source plane. We find that individual pixels within the galaxy follow the same trend found using unresolved observations of other galaxies, indicating that the deficit arises on scales <700 pc. The lensing reconstruction reveals two spatially and kinematically separated components (~1 kpc and ~500 km s^-1 apart) connected by a bridge of gas. Both components are found to be globally unstable, with Toomre Q instability parameters << 1 everywhere. We argue that SPT0346-52 is undergoing a major merger, which is likely driving the intense and compact star formation.

The origin of the high SFR observed in high-z dusty star-forming galaxies is still unknown. Large fractions of dense molecular gas might provide part of the explanation, but there are few observational constraints on the amount of dense gas in high-z systems dominated by star formation. We present the results of our ALMA program targeting dense-gas tracers (HCN(5-4), HCO+(5-4), and HNC(5-4)) in 5 strongly lensed galaxies from the SPT SMG sample. We detected two of these lines (SNR>5) in SPT-125-47 at z=2.51 and tentatively detected all three (SNR~3) in SPT0551-50 at z=3.16. Since a significant fraction of our target lines is not detected, we developed a statistical method to derive unbiased mean properties taking into account both detections and non-detections. On average, the HCN(5-4) and HCO+(5-4) luminosities of our sources are a factor of ~1.7 fainter than expected, based on the local L'HCN(5-4)-LIR relation, but this offset corresponds to only ~2 sigma. We find that both the HCO+/HCN and HNC/HCN flux ratios are compatible with unity. The first ratio is expected for PDRs while the second is consistent with PDRs or XDRs and/or mid-IR pumping of HNC. Our sources are at the high end of the local relation between the star formation efficiency, determined using the LIR/[CI] and LIR/CO ratios, and the dense gas fraction, estimated using the HCN/[CI] and HCN/CO ratios. In SPT0125-47, we found that the velocity profiles of the lines tracing dense (HCN, HCO+) and lower-density (CO, [CI]) gas are similar. In addition to these lines, we obtained one robust and one tentative detection of 13CO(4-3) and found an average I12CO(4-3)/I13CO(4-3) flux ratio of 26.1$_{-3.5}^{+4.5}$, indicating a young but not pristine interstellar medium. We argue that the combination of large and slightly enriched gas reservoirs and high dense-gas fractions could explain the prodigious star formation in these systems.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where the dense gas is traced by the HCN J=4-3 and the HCO^+ J=4-3 emission, measured in our sample of spatially-resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, L_IR/L'_dense, shows systematic variations with L_IR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between L_IR/L'_dense ratio and the warm-dust temperature gauged by the 70 \mu m/100 \mu m flux ratio. We find the luminosity ratios of IR/HCN(4-3) and IR/HCO^+(4-3), which can be taken as a proxy for the efficiency of star formation in the dense molecular gas (SFE_dense), appears to be nearly independent of the dense-gas fraction (f_dense) for our sample of galaxies. The SFE of the total molecular gas (SFE_mol) is found to increase substantially with f_dense when combining our data with that on local (ultra)luminous infrared galaxies and high-z quasars. The mean L'_HCN(4-3)/L'_HCO^+(4-3) line ratio measured for the six targeted galaxies is 0.9+/-0.6. No significant correlation is found for the L'_HCN(4-3)/L'_HCO^+(4-3) ratio with the SFR as traced by L_IR, nor with the warm-dust temperature, for the different populations of galaxies.

Massive galaxy clusters are now found as early as 3 billion years after the Big Bang, containing stars that formed at even earlier epochs. The high-redshift progenitors of these galaxy clusters, termed 'protoclusters', are identified in cosmological simulations with the highest dark matter overdensities. While their observational signatures are less well defined compared to virialized clusters with a substantial hot intra-cluster medium (ICM), protoclusters are expected to contain extremely massive galaxies that can be observed as luminous starbursts. Recent claimed detections of protoclusters hosting such starbursts do not support the kind of rapid cluster core formation expected in simulations because these structures contain only a handful of starbursting galaxies spread throughout a broad structure, with poor evidence for eventual collapse into a protocluster. Here we report that the source SPT2349-56 consists of at least 14 gas-rich galaxies all lying at z = 4.31 based on sensitive observations of carbon monoxide and ionized carbon. We demonstrate that each of these galaxies is forming stars between 50 and 1000 times faster than our own Milky Way, and all are located within a projected region only $\sim$ 130 kiloparsecs in diameter. This galaxy surface density is more than 10 times the average blank field value (integrated over all redshifts) and $>$1000 times the average field volume density. The velocity dispersion ($\sim$ 410 km s$^{-1}$) of these galaxies and enormous gas and star formation densities suggest that this system represents a galaxy cluster core at an advanced stage of formation when the Universe was only 1.4 billion years old. A comparison with other known protoclusters at high redshifts shows that SPT2349-56 is a uniquely massive and dense system that could be building one of the most massive structures in the Universe today.

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

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

According to the current understanding of cosmic structure formation, the precursors of the most massive structures in the Universe began to form shortly after the Big Bang, in regions corresponding to the largest fluctuations in the cosmic density field. Observing these structures during their period of active growth and assembly - the first few hundred million years of the Universe - is challenging because it requires surveys that are sensitive enough to detect the distant galaxies that act as signposts for these structures and wide enough to capture the rarest objects. As a result, very few such objects have been detected so far. Here we report observations of a far-infrared-luminous object at redshift 6.900 (less than 800 Myr after the Big Bang) that was discovered in a wide-field survey. High-resolution imaging reveals this source to be a pair of extremely massive star-forming galaxies. The larger of these galaxies is forming stars at a rate of 2900 solar masses per year, contains 270 billion solar masses of gas and 2.5 billion solar masses of dust, and is more massive than any other known object at a redshift of more than 6. Its rapid star formation is probably triggered by its companion galaxy at a projected separation of just 8 kiloparsecs. This merging companion hosts 35 billion solar masses of stars and has a star-formation rate of 540 solar masses per year, but has an order of magnitude less gas and dust than its neighbor and physical conditions akin to those observed in lower-metallicity galaxies in the nearby Universe. These objects suggest the presence of a dark-matter halo with a mass of more than 400 billion solar masses, making it among the rarest dark-matter haloes that should exist in the Universe at this epoch.

Of the almost 40 star forming galaxies at z>~5 (not counting QSOs) observed in [CII] to date, nearly half are either very faint in [CII], or not detected at all, and fall well below expectations based on locally derived relations between star formation rate (SFR) and [CII] luminosity. Combining cosmological zoom simulations of galaxies with SIGAME (SImulator of GAlaxy Millimeter/submillimeter Emission) we have modeled the multi-phased interstellar medium (ISM) and its emission in [CII], [OI] and [OIII], from 30 main sequence galaxies at z~6 with star formation rates ~3-23Msun/yr, stellar masses ~(0.7-8)x10^9Msun, and metallicities ~(0.1-0.4)xZsun. The simulations are able to reproduce the aforementioned [CII]-faintness at z>5, match two of the three existing z>~5 detections of [OIII], and are furthermore roughly consistent with the [OI] and [OIII] luminosity relations with SFR observed for local starburst galaxies. We find that the [CII] emission is dominated by the diffuse ionized gas phase and molecular clouds, which on average contribute ~66% and ~27%, respectively. The molecular gas, which constitutes only ~10% of the total gas mass is thus a more efficient emitter of [CII] than the ionized gas making up ~85% of the total gas mass. A principal component analysis shows that the [CII] luminosity correlates with the star formation activity as well as average metallicity. The low metallicities of our simulations together with their low molecular gas mass fractions can account for their [CII]-faintness, and we suggest these factors may also be responsible for the [CII]-faint normal galaxies observed at these early epochs.

We report the discovery and constrain the physical conditions of the interstellar medium of the highest-redshift millimeter-selected dusty star-forming galaxy (DSFG) to date, SPT-S J031132-5823.4 (hereafter SPT0311-58), at $z=6.900 +/- 0.002$. SPT0311-58 was discovered via its 1.4mm thermal dust continuum emission in the South Pole Telescope (SPT)-SZ survey. The spectroscopic redshift was determined through an ALMA 3mm frequency scan that detected CO(6-5), CO(7-6) and [CI](2-1), and subsequently confirmed by detections of CO(3-2) with ATCA and [CII] with APEX. We constrain the properties of the ISM in SPT0311-58 with a radiative transfer analysis of the dust continuum photometry and the CO and [CI] line emission. This allows us to determine the gas content without ad hoc assumptions about gas mass scaling factors. SPT0311-58 is extremely massive, with an intrinsic gas mass of $M_{\rm gas} = 3.3 \pm 1.9 \times10^{11}\,M_{\odot}$. Its large mass and intense star formation is very rare for a source well into the Epoch of Reionization.

We present ALMA [CI]($1-0$) (rest frequency 492 GHz) observations for a sample of 13 strongly-lensed dusty star-forming galaxies originally discovered at 1.4mm in a blank-field survey by the South Pole Telescope. We compare these new data with available [CI] observations from the literature, allowing a study of the ISM properties of $\sim 30$ extreme dusty star-forming galaxies spanning a redshift range $2 < z < 5$. Using the [CI] line as a tracer of the molecular ISM, we find a mean molecular gas mass for SPT-DSFGs of $6.6 \times 10^{10}$ M$_{\odot}$. This is in tension with gas masses derived via low-$J$ $^{12}$CO and dust masses; bringing the estimates into accordance requires either (a) an elevated CO-to-H$_2$ conversion factor for our sample of $\alpha_{\rm CO} \sim 2.5$ and a gas-to-dust ratio $\sim200$, or (b) an high carbon abundance $X_{\rm CI} \sim 7\times10^{-5}$. Using observations of a range of additional atomic and molecular lines (including [CI], [CII], and multiple transitions of CO), we use a modern Photodissociation Region code (3D-PDR) to assess the physical conditions (including the density, UV radiation field strength, and gas temperature) within the ISM of the DSFGs in our sample. We find that the ISM within our DSFGs is characterised by dense gas permeated by strong UV fields. We note that previous efforts to characterise PDR regions in DSFGs may have significantly underestimated the density of the ISM. Combined, our analysis suggests that the ISM of extreme dusty starbursts at high redshift consists of dense, carbon-rich gas not directly comparable to the ISM of starbursts in the local Universe.

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.

The South Pole Telescope has discovered one hundred gravitationally lensed, high-redshift, dusty, star-forming galaxies (DSFGs). We present 0.5" resolution 870um Atacama Large Millimeter/submillimeter Array imaging of a sample of 47 DSFGs spanning z=1.9-5.7, and construct gravitational lens models of these sources. Our visibility-based lens modeling incorporates several sources of residual interferometric calibration uncertainty, allowing us to properly account for noise in the observations. At least 70% of the sources are strongly lensed by foreground galaxies (mu_870um > 2), with a median magnification mu_870um = 6.3, extending to mu_870um > 30. We compare the intrinsic size distribution of the strongly lensed sources to a similar number of unlensed DSFGs and find no significant differences in spite of a bias between the magnification and intrinsic source size. This may indicate that the true size distribution of DSFGs is relatively narrow. We use the source sizes to constrain the wavelength at which the dust optical depth is unity and find this wavelength to be correlated with the dust temperature. This correlation leads to discrepancies in dust mass estimates of a factor of 2 compared to estimates using a single value for this wavelength. We investigate the relationship between the [CII] line and the far-infrared luminosity and find that the same correlation between the [CII]L_FIR ratio and Sigma_FIR found for low-redshift star-forming galaxies applies to high-redshift galaxies and extends at least two orders of magnitude higher in Sigma_FIR. This lends further credence to the claim that the compactness of the IR-emitting region is the controlling parameter in establishing the "[CII] deficit."

We use the Atacama Large Millimeter/submillimeter Array (ALMA) in Cycle 1 to determine spectroscopic redshifts of high-redshift dusty star-forming galaxies (DSFGs) selected by their 1.4mm continuum emission in the South Pole Telescope (SPT) survey. We present ALMA 3mm spectral scans between 84-114GHz for 15 galaxies and targeted ALMA 1mm observations for an additional eight sources. Our observations yield 30 new line detections from CO, [CI] , [NII] , H_2O and NH_3. We further present APEX [CII] and CO mid-J observations for seven sources for which only a single line was detected in spectral-scan data from ALMA Cycle 0 or Cycle 1. We combine the new observations with previously published and new mm/submm line and photometric data of the SPT-selected DSFGs to study their redshift distribution. The combined data yield 39 spectroscopic redshifts from molecular lines, a success rate of >85%. Our sample represents the largest data set of its kind today and has the highest spectroscopic completeness among all redshift surveys of high-z DSFGs. The median of the redshift distribution is z=3.9+/-0.4, and the highest-redshift source in our sample is at z=5.8. We discuss how the selection of our sources affects the redshift distribution, focusing on source brightness, selection wavelength, and strong gravitational lensing. We correct for the effect of gravitational lensing and find the redshift distribution for 1.4mm-selected sources with a median redshift of z=3.1+/-0.3. Comparing to redshift distributions selected at shorter wavelengths from the literature, we show that selection wavelength affects the shape of the redshift distribution.

We present SIGAME (SImulator of GAlaxy Millimetre/submillimetre Emission), a new numerical code designed to simulate the 12CO rotational line emission spectrum of galaxies. Using sub-grid physics recipes to post-process the outputs of smoothed particle hydrodynamics (SPH) simulations, a molecular gas phase is condensed out of the hot and partly ionized SPH gas. The gas is subjected to far-UV radiation fields and cosmic ray ionization rates which are set to scale with the local star formation rate volume density. Level populations and radiative transport of the CO lines are solved with the 3-D radiative transfer code LIME. We have applied SIGAME to cosmological SPH simulations of three disc galaxies at z=2 with stellar masses in the range ~(0.5-2)x10^11 Msun and star formation rates ~40-140 Msun/yr. Global CO luminosities and line ratios are in agreement with observations of disc galaxies at z~2 up to and including J=3-2 but falling short of the few existing J=5-4 observations. The central 5 kpc regions of our galaxies have CO 3-2/1-0 and 7-6/1-0 brightness temperature ratios of ~0.55-0.65 and ~0.02-0.08, respectively, while further out in the disc the ratios drop to more quiescent values of ~0.5 and <0.01. Global CO-to-H2 conversion (alpha_CO) factors are ~=1.5 Msun*pc^2/(K km s/1), i.e. ~2-3 times below typically adopted values for disc galaxies, and alpha_CO increases with radius, in agreement with observations of nearby galaxies. Adopting a top-heavy Giant Molecular Cloud (GMC) mass spectrum does not significantly change the results. Steepening the GMC density profile leads to higher global line ratios for J_up>=3 and CO-to-H2 conversion factors [~=3.6 Msun*pc^2/(K km/s)].

We explore the recoverability of gas physical conditions with the Large Velocity Gradient (LVG) model, using the public code RADEX and the molecules HCN and CO. Examining a wide parameter range with a series of models of increasing complexity we use both grid and Monte Carlo Markov Chain (MCMC) methods to recover the input conditions, and quantify the inherent and noise induced uncertainties in the model results. We find that even with the benefit of generous assumptions the LVG models struggle to recover any parameter better than to within half a dex, although we find no evidence of systemic offsets. Examining isotopologue lines we demonstrate that it is always preferable to model the isotopologue abundance ratio as a free parameter, due to large biases introduced in all other parameters when an incorrect ratio is assumed. Finally, we explore the effects of the background radiation temperature on CO and HCN line ratios, with an emphasis on the effect of the CMB at $z>4$, and show that while the effect on the line ratios is minor, the effect on the SLED peak is significant and that the CO$(1-0)$ line luminosity to H$_2$ mass conversion factor ($\alpha_{\rm CO}$) needs to be altered to account for the loss of contrast against the hotter CMB as redshift increases.

Using the Australia Telescope Compact Array (ATCA), we conducted a survey of CO J=1-0 and J=2-1 line emission towards strongly lensed high-redshift dusty star forming galaxies (DSFGs) previously discovered with the South Pole Telescope (SPT). Our sample comprises 17 sources that had CO-based spectroscopic redshifts obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Atacama Pathfinder Experiment (APEX). We detect all sources with known redshifts in either CO J=1-0 or J=2-1. Twelve sources are detected in the 7-mm continuum. The derived CO luminosities imply gas masses in the range (0.5-11)x10^{10} M_sun and gas depletion timescales <200 Myr, using a CO to gas mass conversion factor alpha_CO=0.8 M_sun (K km/s pc^2)^{-1}. Combining the CO luminosities and dust masses, along with a fixed gas-to-dust ratio, we derive alpha_CO factors in the range 0.4-1.8, similar to what is found in other starbursting systems. We find small scatter in alpha_CO values within the sample, even though inherent variations in the spatial distribution of dust and gas in individual cases could bias the dust-based alpha_CO estimates. We find that lensing magnification factors based on the CO linewidth to luminosity relation (mu_CO) are highly unreliable, but particularly when mu<5. Finally, comparison of the gas and dynamical masses suggest that the average molecular gas fraction stays relatively constant at z=2-5 in the SPT DSFG sample.

We present ALMA detections of the [NII] 205$\mu$m and CO(12-11) emission lines, and the tentative detection of [CI] $^3$P$_1$ - $^3$P$_0$ for the strongly lensed ($\mu$=5.7$\pm$0.5) dusty, star-forming galaxy SPT-S J213242-5802.9 (hereafter SPT2132-58) at z=4.77. The [NII] and CO(12-11) lines are detected at 11.5 and 8.5 $\sigma$ level, respectively, by our band-6 observations. The [CI] line is detected at 3.2 $\sigma$ after a re-analysis of existing band-3 data. The [CI] luminosity implies a gas mass of 3.8$\pm$1.2$\times$10$^{10}$ M$_\odot$, and consequently a very short depletion timescale of 34$\pm$13 Myr and a CO-luminosity-to-gas-mass conversion factor $\alpha_{\rm CO}$ of 1.0$\pm$0.3 M$_\odot$ (K km s$^{-1}$ pc$^{2}$)$^{-1}$. SPT2132-58 is an extreme starburst with an intrinsic star formation rate of 1100$\pm$200 M$_\odot$/yr. We find a [CII]/[NII] ratio of 26$\pm$7, which is the highest reported at z$>$4. This suggests that SPT2132-58 hosts an evolved interstellar medium (0.5 Z$_\odot$<Z<1.5 Z$_\odot$), which may be dominated by photodissociation regions. The CO(2-1) and CO(5-4) transitions have lower CO-to-far-infrared ratios than local and high-redshift samples, while CO(12-11) is similar to these samples, suggesting the presence of an additional very excited component or an AGN.

We present S\'IGAME simulations of the [CII]157.7$\mu$ fine structure line emission from cosmological smoothed particle hydrodynamics (SPH) simulations of seven main sequence galaxies at z=2. Using sub-grid physics prescriptions the gas in our simulations is modeled as a multi-phased interstellar medium (ISM) comprised of molecular gas residing in giant molecular clouds, an atomic gas phase associated with photo-dissociation regions (PDRs) at the cloud surfaces, and a diffuse, ionized gas phase. Adopting logotropic cloud density profiles and accounting for heating by the local FUV radiation field and cosmic rays by scaling both with local star formation rate (SFR) volume density, we calculate the [CII] emission using a photon escape probability formalism. The [CII] emission peaks in the central $\lesssim$1 kpc of our galaxies as do the SFR radial profiles, with most [CII] ($\gtrsim$70%) originating in the molecular gas phase, whereas further out ($\gtrsim$2 kpc), the atomic/PDR gas dominates ($\gtrsim$90%) the [CII] emission, no longer tracing on-going star formation. Throughout, the ionized gas contribution is negligible ($\lesssim$3%). The [CII] luminosity vs. SFR ([CII]-SFR) relationship, integrated as well as spatially resolved (on scales of 1 kpc), delineated by our simulated galaxies is in good agreement with the corresponding relations observed locally and at high redshifts. In our simulations, the molecular gas dominates the [CII] budget at SFR$\gtrsim$20 M$_{\odot}$yr$^{-1}$ (SFR density $\gtrsim$0.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$), while atomic/PDR gas takes over at lower SFRs, suggesting a picture in which [CII] predominantly traces the molecular gas in high-density/pressure regions where star formation is on-going, and otherwise reveals the atomic/PDR gas phase.

We present the first observations of H$^{13}$CN$(1-0)$, H$^{13}$CO$^+(1-0)$ and SiO$(2-1)$ in NGC\,6240, obtained with the IRAM PdBI. Combining a Markov Chain Monte Carlo (MCMC) code with Large Velocity Gradient (LVG) modelling, and with additional data from the literature, we simultaneously fit three gas phases and six molecular species to constrain the physical condition of the molecular gas, including mass$-$luminosity conversion factors. We find $\sim10^{10}M_\odot$ of dense molecular gas in cold, dense clouds ($T_{\rm k}\sim10$\,K, $n_{{\rm H}_2}\sim10^6$\,cm$^{-3}$) with a volume filling factor $<0.002$, embedded in a shock heated molecular medium ($T_{\rm k}\sim2000$\,K, $n_{{\rm H}_2}\sim10^{3.6}$\,cm$^{-3}$), both surrounded by an extended diffuse phase ($T_{\rm k}\sim200$\,K, $n_{{\rm H}_2}\sim10^{2.5}$\,cm$^{-3}$). We derive a global $\alpha_{\rm CO}=1.5^{7.1}_{1.1}$ with gas masses $\log_{10}\left(M / [M_\odot]\right)=10.1_{10.0}^{10.8}$, dominated by the dense gas. We also find $\alpha_{\rm HCN} = 32^{89}_{13}$, which traces the cold, dense gas. The [$^{12}$C]/[$^{13}$C] ratio is only slightly elevated ($98^{230}_{65}$), contrary to the very high [CO]/[$^{13}$CO] ratio (300-500) reported in the literature. However, we find very high [HCN]/[H$^{13}$CN] and [HCO$^+$]/[H$^{13}$CO$^+$] abundance ratios $(300^{500}_{200})$ which we attribute to isotope fractionation in the cold, dense clouds.

We probe star formation in the environments of massive $\sim10^{13}\,M_{\odot}$ dark matter halos at redshifts of $z$$\sim$$1$. This star formation is linked to a sub-millimetre clustering signal which we detect in maps of the Planck High Frequency Instrument that are stacked at the positions of a sample of high-redshift ($z$$>$$2$) strongly-lensed dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope (SPT) 2500 deg$^2$ survey. The clustering signal has sub-millimetre colours which are consistent with the mean redshift of the foreground lensing halos ($z$$\sim$$1$). We report a mean excess of star formation rate (SFR) compared to the field, of $(2700\pm700)\,M_{\odot}\,{yr}^{-1}$ from all galaxies contributing to this clustering signal within a radius of 3.5' from the SPT DSFGs. The magnitude of the Planck excess is in broad agreement with predictions of a current model of the cosmic infrared background. The model predicts that 80$\%$ of the excess emission measured by Planck originates from galaxies lying in the neighbouring halos of the lensing halo. Using Herschel maps of the same fields, we find a clear excess, relative to the field, of individual sources which contribute to the Planck excess. The mean excess SFR compared to the field is measured to be ($370\pm40)$$\,M_{\odot}\,{yr}^{-1}$ per resolved, clustered source. Our findings suggest that the environments around these massive $z$$\sim$$1$ lensing halos host intense star formation out to about $2\,$Mpc. The flux enhancement due to clustering should also be considered when measuring flux densities of galaxies in Planck data.

To understand cosmic mass assembly in the Universe at early epochs, we primarily rely on measurements of stellar mass and star formation rate of distant galaxies. In this paper, we present stellar masses and star formation rates of six high-redshift ($2.8\leq z \leq 5.7$) dusty, star-forming galaxies (DSFGs) that are strongly gravitationally lensed by foreground galaxies. These sources were first discovered by the South Pole Telescope (SPT) at millimeter wavelengths and all have spectroscopic redshifts and robust lens models derived from ALMA observations. We have conducted follow-up observations, obtaining multi-wavelength imaging data, using {\it HST}, {\it Spitzer}, {\it Herschel} and the Atacama Pathfinder EXperiment (APEX). We use the high-resolution {\it HST}/WFC3 images to disentangle the background source from the foreground lens in {\it Spitzer}/IRAC data. The detections and upper limits provide important constraints on the spectral energy distributions (SEDs) for these DSFGs, yielding stellar masses, IR luminosities, and star formation rates (SFRs). The SED fits of six SPT sources show that the intrinsic stellar masses span a range more than one order of magnitude with a median value $\sim$ 5 $\times 10^{10}M_{\Sun}$. The intrinsic IR luminosities range from 4$\times 10^{12}L_{\Sun}$ to 4$\times 10^{13}L_{\Sun}$. They all have prodigious intrinsic star formation rates of 510 to 4800 $M_{\Sun} {\rm yr}^{-1}$. Compared to the star-forming main sequence (MS), these six DSFGs have specific SFRs that all lie above the MS, including two galaxies that are a factor of 10 higher than the MS. Our results suggest that we are witnessing the ongoing strong starburst events which may be driven by major mergers.

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($\nu_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $\nu_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($\nu$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.

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

We study the global SF law - the relation between gas and SFRs in a sample of 181 local galaxies with L_IR spanning almost five orders of magnitude, which includes 115 normal galaxies and 66 (U)LIRGs. We derive their atomic, molecular gas and dense molecular gas masses using newly available HI, CO and HCN data from the literature, and SFRs are determined both from total IR and 1.4 GHz radio continuum (RC) luminosities. In order to derive the disk-averaged surface densities of gas and SFRs, we have used high-resolution RC observations to measure the radio sizes for all galaxies. We find that dense molecular gas (as traced by HCN) has the tightest correlation with that of SFRs, and is linear in (N=1.01 +/- 0.02) across the full galaxy sample. The correlation between densities of molecular gas (traced by CO) and SFRs is sensitive to the adopted value of the alpha_CO used to infer molecular gas masses from CO luminosities. For a fixed value of alpha_CO, a slope of 1.14+/-0.02 is found. If instead we adopt values of 4.6 and 0.8 for disk galaxies and (U)LIRGs, respectively, we find the two distinct relations. If applying a continuously varying alpha_CO to our sample, we recover a single relation with slope of 1.60+/-0.03. The SFRs is a steeper function of total gas than that of molecular gas, and is tighter among low-luminosity galaxies. We find no correlation between SFRs and atomic gas.

We present [CII] observations of 20 strongly lensed dusty star forming galaxies at 2.1 < z < 5.7 using APEX and Herschel. The sources were selected on their 1.4 mm flux (S_1.4mm > 20 mJy) from the South Pole Telescope survey, with far-infrared (FIR) luminosities determined from extensive photometric data. The [CII] line is robustly detected in 17 sources, all but one being spectrally resolved. Eleven out of 20 sources observed in [CII] also have low-J CO detections from ATCA. A comparison with mid- and high-J CO lines from ALMA reveals consistent [CII] and CO velocity profiles, suggesting that there is little differential lensing between these species. The [CII], low-J CO and FIR data allow us to constrain the properties of the interstellar medium. We find [CII] to CO(1-0) luminosity ratios in the SPT sample of 5200 +- 1800, with significantly less scatter than in other samples. This line ratio can be best described by a medium of [CII] and CO emitting gas with a higher [CII] than CO excitation temperature, high CO optical depth tau_CO >> 1, and low to moderate [CII] optical depth tau_CII ~< 1. The geometric structure of photodissociation regions allows for such conditions.

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $\mu$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $\mu$m, [OI] 63 $\mu$m, and [OI] 145 $\mu$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $\mu$m, and [CI] 609 $\mu$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $\mu$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $\mu$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.

We present the results of interferometric spectral line observations of Arp 220 at 3.5mm and 1.2mm from the Plateau de Bure Interferometer (PdBI), imaging the two nuclear disks in H$^{13}$CN$(1 - 0)$ and $(3 - 2)$, H$^{13}$CO$^+(1 - 0)$ and $(3 - 2)$, and HN$^{13}$C$(3 - 2)$ as well as SiO$(2 - 1)$ and $(6 - 5)$, HC$^{15}$N$(3 - 2)$, and SO$(6_6 - 5_5)$. The gas traced by SiO$(6 - 5)$ has a complex and extended kinematic signature including a prominent P Cygni profile, almost identical to previous observations of HCO$^+(3 - 2)$. Spatial offsets $0.1''$ north and south of the continuum centre in the emission and absorption of the SiO$(6 - 5)$ P Cygni profile in the western nucleus (WN) imply a bipolar outflow, delineating the northern and southern edges of its disk and suggesting a disk radius of $\sim40$ pc, consistent with that found by ALMA observations of Arp 220. We address the blending of SiO$(6 - 5)$ and H$^{13}$CO$^+(3 - 2)$ by considering two limiting cases with regards to the H$^{13}$CO$^+$ emission throughout our analysis. Large velocity gradient (LVG) modelling is used to constrain the physical conditions of the gas and to infer abundance ratios in the two nuclei. Our most conservative lower limit on the [H$^{13}$CN]/[H$^{13}$CO$^+$] abundance ratio is 11 in the WN, cf. 0.10 in the eastern nucleus (EN). Comparing these ratios to the literature we argue on chemical grounds for an energetically significant AGN in the WN driving either X-ray or shock chemistry, and a dominant starburst in the EN.

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

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

We present FIR-CO luminosity relations ($\log L_{\rm FIR} = \alpha \log L'_{\rm CO} + \beta$) for the full CO rotational ladder from J=1-0 to J=13-12 for 62 local (z < 0.1) (Ultra) Luminous Infrared Galaxies (LIRGs) using data from Herschel SPIRE-FTS and ground-based telescopes. We extend our sample to high redshifts (z > 1) by including 35 (sub)-millimeter selected dusty star forming galaxies from the literature with robust CO observations. The addition of luminous starbursts at high redshifts enlarge the range of the FIR-CO luminosity relations towards the high-IR-luminosity end while also significantly increasing the small amount of mid-/high-J CO line data available prior to Herschel. This new data-set (both in terms of IR luminosity and J-ladder) reveals linear FIR-CO luminosity relations ($\alpha \sim 1$) for J=1-0 up to J=5-4, with a nearly constant normalisation ($\beta \sim 2$). This is expected from the (also) linear FIR-(molecular line) relations found for the dense gas tracer lines (HCN and CS), as long as the dense gas mass fraction does not vary strongly within our (merger/starburst)-dominated sample. However from J=6-5 and up to J=13-12 we find an increasingly sub-linear slope and higher normalization constant with increasing J. We argue that these are caused by a warm (~100K) and dense ($>10^4{\rm cm^{-3}}$) gas component whose thermal state is unlikely to be maintained by star formation powered far-UV radiation fields (and thus is no longer directly tied to the star formation rate). We suggest that mechanical heating (e.g., supernova driven turbulence and shocks), and not cosmic rays, is the more likely source of energy for this component. The global CO spectral line energy distributions (SLEDs), which remain highly excited from J=6-5 up to J=13-12, are found to be a generic feature of the (U)LIRGs in our sample, and further support the presence of this gas component.

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

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

We used the SPIRE/FTS instrument aboard the Herschel Space Observatory (HSO) to obtain the Spectral Line Energy Distributions (SLEDs) of CO from J=4-3 to J=13-12 of Arp 193 and NGC 6240, two classical merger/starbursts selected from our molecular line survey of local Luminous Infrared Galaxies (LIRGs: L_{IR}>=10^{11} L_{sol}). The high-J CO SLEDs are then combined with ground-based low-J CO, {13}CO, HCN, HCO+, CS line data and used to probe the thermal and dynamical states of their large molecular gas reservoirs. We find the two CO SLEDs strongly diverging from J=4-3 onwards, with NGC6240 having a much higher CO line excitation than Arp193, despite their similar low-J CO SLEDs and L_{FIR}/L_{CO,1-0}, L_{HCN}/L_{CO} (J=1-0) ratios (proxies of star formation efficiency and dense gas mass fraction). In Arp193, one of the three most extreme starbursts in the local Universe, the molecular SLEDs indicate a small amount ~(5-15)% of dense gas (n>=10^{4}cm^{-3}) unlike NGC6240 where most of the molecular gas (~(60-70)%) is dense n~(10^4-10^5)cm^{-3}. Strong star-formation feedback can drive this disparity in their dense gas mass fractions, and also induce extreme thermal and dynamical states for the molecular gas.In NGC6240, and to a lesser degree in Arp193, we find large molecular gas masses whose thermal states cannot be maintained by FUV photons from Photon Dominated Regions (PDRs). We argue that this may happen often in metal-rich merger/starbursts, strongly altering the initial conditions of star formation. ALMA can now directly probe these conditions across cosmic epoch, and even probe their deeply dust-enshrouded outcome, the stellar IMF averaged over galactic evolution.

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 present the average rest-frame spectrum of high-redshift dusty, star-forming galaxies from 250-770GHz. This spectrum was constructed by stacking ALMA 3mm spectra of 22 such sources discovered by the South Pole Telescope and spanning z=2.0-5.7. In addition to multiple bright spectral features of 12CO, [CI], and H2O, we also detect several faint transitions of 13CO, HCN, HNC, HCO+, and CN, and use the observed line strengths to characterize the typical properties of the interstellar medium of these high-redshift starburst galaxies. We find that the 13CO brightness in these objects is comparable to that of the only other z>2 star-forming galaxy in which 13CO has been observed. We show that the emission from the high-critical density molecules HCN, HNC, HCO+, and CN is consistent with a warm, dense medium with T_kin ~ 55K and n_H2 >~ 10^5.5 cm^-3. High molecular hydrogen densities are required to reproduce the observed line ratios, and we demonstrate that alternatives to purely collisional excitation are unlikely to be significant for the bulk of these systems. We quantify the average emission from several species with no individually detected transitions, and find emission from the hydride CH and the linear molecule CCH for the first time at high redshift, indicating that these molecules may be powerful probes of interstellar chemistry in high-redshift systems. These observations represent the first constraints on many molecular species with rest-frame transitions from 0.4-1.2mm in star-forming systems at high redshift, and will be invaluable in making effective use of ALMA in full science operations.

We present a CO(1-0) survey for cold molecular gas in a representative sample of 13 high-z radio galaxies (HzRGs) at 1.4<z<2.8, using the Australia Telescope Compact Array. We detect CO(1-0) emission associated with five sources: MRC 0114-211, MRC 0152-209, MRC 0156-252, MRC 1138-262 and MRC 2048-272. The CO(1-0) luminosities are in the range $L'_{\rm CO} \sim (5 - 9) \times 10^{10}$ K km/s pc$^{2}$. For MRC 0152-209 and MRC 1138-262 part of the CO(1-0) emission coincides with the radio galaxy, while part is spread on scales of tens of kpc and likely associated with galaxy mergers. The molecular gas mass derived for these two systems is M$_{\rm H2} \sim 6 \times 10^{10}\, {\rm M}_{\odot}$ (M$_{\rm H2}$/$L'_{\rm CO}$=0.8). For the remaining three CO-detected sources, the CO(1-0) emission is located in the halo (~50-kpc) environment. These three HzRGs are among the fainter far-IR emitters in our sample, suggesting that similar reservoirs of cold molecular halo gas may have been missed in earlier studies due to pre-selection of IR-bright sources. In all three cases the CO(1-0) is aligned along the radio axis and found beyond the brightest radio hot-spot, in a region devoid of 4.5$\mu$m emission in Spitzer imaging. The CO(1-0) profiles are broad, with velocity widths of ~ 1000 - 3600 km/s. We discuss several possible scenarios to explain these halo reservoirs of CO(1-0). Following these results, we complement our CO(1-0) study with detections of extended CO from the literature and find at marginal statistical significance (95% level) that CO in HzRGs is preferentially aligned towards the radio jet axis. For the eight sources in which we do not detect CO(1-0), we set realistic upper limits of $L'_{\rm CO} \sim 3-4 \times 10^{10}$ K km/s pc$^{2}$. Our survey reveals a CO(1-0) detection rate of 38%, allowing us to compare the CO(1-0) content of HzRGs with that of other types of high-z galaxies.

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.

We present observations of SPT-S J053816-5030.8, a gravitationally-lensed dusty star forming galaxy (DSFG) at z = 2.7817, first discovered at millimeter wavelengths by the South Pole Telescope. SPT 0538-50 is typical of the brightest sources found by wide-field millimeter-wavelength surveys, being lensed by an intervening galaxy at moderate redshift (in this instance, at z = 0.441). We present a wide array of multi-wavelength spectroscopic and photometric data on SPT 0538-50, including data from ALMA, Herschel PACS and SPIRE, Hubble, Spitzer, VLT, ATCA, APEX, and the SMA. We use high resolution imaging from HST to de-blend SPT 0538-50, separating DSFG emission from that of the foreground lens. Combined with a source model derived from ALMA imaging (which suggests a magnification factor of 21 +/- 4), we derive the intrinsic properties of SPT 0538-50, including the stellar mass, far-IR luminosity, star formation rate, molecular gas mass, and - using molecular line fluxes - the excitation conditions within the ISM. The derived physical properties argue that we are witnessing compact, merger-driven star formation in SPT 0538-50, similar to local starburst galaxies, and unlike that seen in some other DSFGs at this epoch.

We report the detection of CO(1-0) line emission in the bright, lensed star-forming galaxies SPT-S 233227-5358.5 (z=2.73) and SPT-S 053816-5030.8 (z=2.78), using the Australia Telescope Compact Array (ATCA). Both galaxies were discovered in a large-area millimeter survey with the South Pole Telescope (SPT) and found to be gravitationally lensed by intervening structures. The measured CO intensities imply galaxies with molecular gas masses of (3.2 \pm 0.5)x10^10 (mu/15)^{-1}(X_CO/0.8) M_sun and (1.7 \pm 0.3)x10^10 (mu/20)^{-1}(X_CO/0.8) M_sun, and gas depletion timescales of 4.9x10^7 (X_CO/0.8) yr and 2.6x10^7 (X_CO/0.8) yr, respectively, where mu corresponds to the lens magnification and X_CO is the CO luminosity to gas mass conversion factor. In the case of SPT-S 053816-5030.8, we also obtained significant detections of the rest-frame 115.7 and 132.4 GHz radio continuum. Based on the radio to infrared spectral energy distribution and an assumed synchrotron spectral index, we find that 42 \pm 10 % and 55 \pm 13 % of the flux at rest-frame 115.7 and 132.4 GHz arise from free-free emission. We find a radio-derived intrinsic star formation rate (SFR) of 470 \pm 170 M_sun/yr, consistent within the uncertainties with the infrared estimate. Based on the morphology of this object in the source plane, the derived gas mass and the possible flattening of the radio spectral index towards low frequencies, we argue that SPT-S 053816-5030.8 exhibits properties compatible with a scaled-up local ultra-luminous infrared galaxy.

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.

In the past decade, our understanding of galaxy evolution has been revolutionized by the discovery that luminous, dusty, starburst galaxies were 1,000 times more abundant in the early Universe than at present. It has, however, been difficult to measure the complete redshift 2 distribution of these objects, especially at the highest redshifts (z > 4). Here we report a redshift survey at a wavelength of three millimeters, targeting carbon monoxide line emission from the star-forming molecular gas in the direction of extraordinarily bright millimetrewave-selected sources. High-resolution imaging demonstrates that these sources are strongly gravitationally lensed by foreground galaxies. We detect spectral lines in 23 out of 26 sources and multiple lines in 12 of those 23 sources, from which we obtain robust, unambiguous redshifts. At least 10 of the sources are found to lie at z > 4, indicating that the fraction of dusty starburst galaxies at high redshifts is greater than previously thought. Models of lens geometries in the sample indicate that the background objects are ultra-luminous infrared galaxies, powered by extreme bursts of star formation.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 860 micrometer imaging of four high-redshift (z=2.8-5.7) dusty sources that were detected using the South Pole Telescope (SPT) at 1.4 mm and are not seen in existing radio to far-infrared catalogs. At 1.5 arcsec resolution, the ALMA data reveal multiple images of each submillimeter source, separated by 1-3 arcsec, consistent with strong lensing by intervening galaxies visible in near-IR imaging of these sources. We describe a gravitational lens modeling procedure that operates on the measured visibilities and incorporates self-calibration-like antenna phase corrections as part of the model optimization, which we use to interpret the source structure. Lens models indicate that SPT0346-52, located at z=5.7, is one of the most luminous and intensely star-forming sources in the universe with a lensing corrected FIR luminosity of 3.7 X 10^13 L_sun and star formation surface density of 4200 M_sun yr^-1 kpc^-2. We find magnification factors of 5 to 22, with lens Einstein radii of 1.1-2.0 arcsec and Einstein enclosed masses of 1.6-7.2x10^11 M_sun. These observations confirm the lensing origin of these objects, allow us to measure the their intrinsic sizes and luminosities, and demonstrate the important role that ALMA will play in the interpretation of lensed submillimeter sources.

Using the Atacama Large Millimeter/submillimeter Array (ALMA), we have conducted a blind redshift survey in the 3 mm atmospheric transmission window for 26 strongly lensd dusty star-forming galaxies (DSFGs) selected with the South Pole Telescope (SPT). The sources were selected to have S_1.4mm>20 mJy and a dust-like spectrum and, to remove low-z sources, not have bright radio (S_843MHz<6mJy) or far-infrared counterparts (S_100um<1 Jy, S_60um<200mJy). We robustly detect 44 line features in our survey, which we identify as redshifted emission lines of 12CO, 13CO, [CI], H2O, and H2O+. We find one or more spectral features in 23 sources yielding a ~90% detection rate for this survey; in 12 of these sources we detect multiple lines, while in 11 sources we detect only a single line. For the sources with only one detected line, we break the redshift degeneracy with additional spectroscopic observations if available, or infer the most likely line identification based on photometric data. This yields secure redshifts for ~70% of the sample. The three sources with no lines detected are tentatively placed in the redshift desert between 1.7<z<2.0. The resulting mean redshift of our sample is <z>=3.5. This finding is in contrast to the redshift distribution of radio-identified DSFGs, which have a significantly lower mean redshift of <z>=2.3 and for which only 10-15% of the population is expected to be at z>3. We discuss the effect of gravitational lensing on the redshift distribution and compare our measured redshift distribution to that of models in the literature.

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.

We present Herschel observations at 70, 160, 250, 350 and 500 micron of the environment of the radio galaxy 4C+41.17 at z = 3.792. About 65% of the extracted sources are securely identified with mid-IR sources observed with the Spitzer Space Telescope at 3.6, 4.5, 5.8, 8 and 24 micron. We derive simple photometric redshifts, also including existing 850 micron and 1200 micron data, using templates of AGN, starburst-dominated systems and evolved stellar populations. We find that most of the Herschel sources are foreground to the radio galaxy and therefore do not belong to a structure associated with 4C+41.17. We do, however, find that the SED of the closest (~ 25" offset) source to the radio galaxy is fully consistent with being at the same redshift as 4C+41.17. We show that finding such a bright source that close to the radio galaxy at the same redshift is a very unlikely event, making the environment of 4C+41.17 a special case. We demonstrate that multi-wavelength data, in particular on the Rayleigh-Jeans side of the spectral energy distribution, allow us to confirm or rule out the presence of protocluster candidates that were previously selected by single wavelength data sets.

Results from a large, multi-J CO, {13}CO, and HCN line survey of Luminous Infrared Galaxies (L_{IR}>=10^{10} L_{\odot}) in the local Universe (z<=0.1), complemented by CO J=4--3 up to J=13--12 observations from the Herschel Space Observatory (HSO), paints a new picture for the average conditions of the molecular gas of the most luminous of these galaxies with turbulence and/or large cosmic ray (CR) energy densities U_{CR} rather than far-UV/optical photons from star-forming sites as the dominant heating sources. Especially in ULIRGs (L_{IR}>10^{12} L_{\odot}) the Photon Dominated Regions (PDRs) can encompass at most \sim few% of their molecular gas mass while the large U_{CR} and the strong turbulence in these merger/starbursts, can volumetrically heat much of their molecular gas to T_{kin}\sim(100-200)K, unhindered by the high dust extinctions. Moreover the strong supersonic turbulence in ULIRGs relocates much of their molecular gas at much higher average densities than in isolated spirals. This renders low-J CO lines incapable of constraining the properties of the bulk of the molecular gas in ULIRGs, with substantial and systematic underestimates of its mass possible when only such lines are used. A comparative study of multi-J HCN lines and CO SLEDs from J=1--0 up to J=13--12 of NGC 6240 and Arp 193 offers a clear example of two merger/starbursts whose similar low-J CO SLEDs, and L_{IR}/L_{CO,1-0}, L_{HCN, 1-0}/L_{CO,1-0} ratios, yield no indications about their strongly diverging CO SLEDs beyond J=4--3, and ultimately the different physical conditions in their molecular ISM. The much larger sensitivity of ALMA and its excellent site in the Atacama desert now allows the observations necessary to ....

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

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

We present a detailed study of the infrared spectral energy distribution of the high-redshift radio galaxy MRC 1138-26 at z = 2.156, also known as the Spiderweb Galaxy. By combining photometry from Spitzer, Herschel and LABOCA we fit the rest-frame 5-300 um emission using a two component, starburst and active galactic nucleus (AGN), model. The total infrared (8 - 1000 um) luminosity of this galaxy is (1.97+/-0.28)x10^13 Lsun with (1.17+/-0.27) and (0.79+/-0.09)x10^13 Lsun due to the AGN and starburst components respectively. The high derived AGN accretion rate of \sim20% Eddington, and the measured star formation rate (SFR) of 1390pm150 Msun/yr, suggest that this massive system is in a special phase of rapid central black hole and host galaxy growth, likely caused by a gas rich merger in a dense environment. The accretion rate is sufficient to power both the jets and the previously observed large outflow. The high SFR and strong outflow suggest this galaxy could potentially exhaust its fuel for stellar growth in a few tens of Myr, although the likely merger of the radio galaxy with nearby satellites suggest bursts of star formation may recur again on time scales of several hundreds of Myr. The age of the radio lobes implies the jet started after the current burst of star formation, and therefore we are possibly witnessing the transition from a merger-induced starburst phase to a radio-loud AGN phase. We also note tentative evidence for [CII]158um emission. This paper marks the first results from the Herschel Galaxy Evolution Project (Project HeRGE), a systematic study of the evolutionary state of 71 high redshift, 1 < z < 5.2, radio galaxies.

We present APEX SABOCA 350micron and LABOCA 870micron observations of 11 representative examples of the rare, extremely bright (S_1.4mm > 15mJy), dust-dominated millimeter-selected galaxies recently discovered by the South Pole Telescope (SPT). All 11 sources are robustly detected with LABOCA with 40 < S_870micron < 130mJy, approximately an order of magnitude higher than the canonical submillimeter galaxy (SMG) population. Six of the sources are also detected by SABOCA at >3sigma, with the detections or upper limits providing a key constraint on the shape of the spectral energy distribution (SED) near its peak. We model the SEDs of these galaxies using a simple modified blackbody and perform the same analysis on samples of SMGs of known redshift from the literature. These calibration samples inform the distribution of dust temperature for similar SMG populations, and this dust temperature prior allows us to derive photometric redshift estimates and far infrared luminosities for the sources. We find a median redshift of <z> = 3.0, higher than the <z> = 2.2 inferred for the normal SMG population. We also derive the apparent size of the sources from the temperature and apparent luminosity, finding them to appear larger than our unlensed calibration sample, which supports the idea that these sources are gravitationally magnified by massive structures along the line of sight.

We report new, sensitive observations of two z ~ 3-3.5 FIR-luminous radio galaxies, 6C1909+72 and B3J2330+3927, in 12CO J=1-0 with the Karl Jansky VLA and at 100-500um using Herschel, alongside new and archival 12CO J=4-3 observations from IRAM PdBI. We introduce a new colour-colour diagnostic plot to constrain the redshifts of several distant, dusty galaxies in our target fields. A bright SMG near 6C1909+72 likely shares the same node or filament as the signpost AGN, but it is not detected in CO despite ~20,000 km/s of velocity coverage. Also in the 6C1909+72 field, a large, red dust feature spanning ~500 kpc is aligned with the radio jet. We suggest several processes by which metal-rich material may have been transported, favouring a collimated outflow reminiscent of the jet-oriented metal enrichment seen in local cluster environments. Our interferometric imaging reveals a gas-rich companion to B3J2330+3927; indeed, all bar one of the eight z >~ 2 radio galaxies (or companions) detected in CO provide some evidence that starburst activity in radio-loud AGN at high redshift is driven by the interaction of two or more gas-rich systems in which a significant mass of stars has already formed, rather than via steady accretion of cold gas from the cosmic web. We find that the CO Tb ratios in radio-loud AGN host galaxies are significantly higher than those seen in similarly intense starbursts where AGN activity is less pronounced. Our most extreme example, where L'(CO4-3)/L'(CO1-0) > 2.7, provides evidence that significant energy is being deposited rapidly into the molecular gas via X-rays and/or mechanical (`quasar-mode') feedback from the AGN, leading to a high degree of turbulence globally and a low optical depth in 12CO - feedback that may lead to the cessation of star formation on a timescale commensurate with that of the jet activity, <~10 Myr.

We present the results from a survey for 12CO emission in 40 luminous sub-millimetre galaxies (SMGs), with 850um fluxes of S850 = 4 - 20 mJy, conducted with the Plateau de Bure Interferometer. We detect 12CO emission in 32 SMGs at z~1.2 - 4.1, including 16 SMGs not previously published. Using multiple 12CO line (J_up =2 - 7) observations, we derive a median spectral line energy distribution for luminous SMGs and use this to estimate a mean gas mass of (5.3 +/- 1.0) \times 10^10 Msun. We report the discovery of a fundamental relationship between 12CO FWHM and 12CO line luminosity in high-redshift starbursts, which we interpret as a natural consequence of the baryon-dominated dynamics within the regions probed by our observations. We use far-infrared luminosities to assess the star-formation efficiency in our SMGs, finding a steepening of the L'CO-LFIR relation as a function of increasing 12CO J_up transition. We derive dynamical masses and molecular gas masses, and use these to determine the redshift evolution of the gas content of SMGs, finding that they do not appear to be significantly more gas rich than less vigorously star-forming galaxies at high redshifts. Finally, we collate X-ray observations, and study the interdependence of gas and dynamical properties of SMGs with their AGN activity and supermassive black hole masses (MBH), finding that SMGs lie significantly below the local M_BH-sigma relation. We conclude that SMGs represent a class of massive, gas-rich ultraluminous galaxies with somewhat heterogeneous properties, ranging from starbursting disc-like systems with L~10^12 L_sun, to the most highly star-forming mergers in the Universe.

We determine star formation rates (SFRs) in a sample of color-selected, star-forming (sBzK) galaxies (K(AB)<21.8) in the Extended Chandra Deep Field - South. To identify and avoid active galactic nuclei, we use X-ray, IRAC color, and IR/radio flux ratio selection methods. Photometric redshift-binned, average flux densities are measured with stacking analyses in Spitzer-MIPS IR, BLAST and APEX/LABOCA submillimeter, VLA and GMRT radio and Chandra X-ray data. We include averages of aperture fluxes in MUSYC UBVRIz'JHK images to determine UV-through-radio spectral energy distributions. We determine the total IR luminosities and compare SFR calibrations from FIR, 24 micron, UV, radio and X-ray wavebands. We find consistency with our best estimator, SFR(IR+UV), to within errors for the preferred radio SFR calibration. Our results imply that 24 micron only and X-ray SFR estimates should be applied to high redshift galaxies with caution. Average IR luminosities are consistent with luminous infrared galaxies. We find SFR(IR+UV) for our stacked sBzKs at median redshifts 1.4, 1.8, and 2.2 to be 55+/-6 (random error), 74+/-8 and 154+/-17 Msun yr^-1 respectively, with additional systematic uncertainty of a factor of ~2.

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

[abridged] We present interferometric CO observations of twelve z~2 submillimetre-faint, star-forming radio galaxies (SFRGs) which are thought to be ultraluminous infrared galaxies (ULIRGs) possibly dominated by warmer dust (T_dust ~> 40 K) than submillimetre galaxies (SMGs) of similar luminosities. Four other CO-observed SFRGs are included from the literature, and all observations are taken at the Plateau de Bure Interferometer (PdBI) in the compact configuration. Ten of the sixteen SFRGs observed in CO (63%) are detected at >4sigma with a mean inferred molecular gas mass of ~2*10^10 M_sun. SFRGs trend slightly above the local ULIRG L_FIR-L'_CO relation. Since SFRGs are about two times fainter in radio luminosity but exhibit similar CO luminosities to SMGs, this suggests SFRGs are slightly more efficient star formers than SMGs at the same redshifts. SFRGs also have a narrow mean CO line width, 320+-80km/s. SFRGs bridge the gap between properties of very luminous >5*10^12 L_sun SMGs and those of local ULIRGs and are consistent with intermediate stage major mergers. We suspect that more moderate-luminosity SMGs, not yet surveyed in CO, would show similar molecular gas properties to SFRGs. The AGN fraction of SFRGs is consistent with SMGs and is estimated to be 0.3+-0.1, suggesting that SFRGs are observed near the peak phase of star formation activity and not in a later, post-SMG enhanced AGN phase. This CO survey of SFRGs serves as a pilot project for the much more extensive survey of Herschel and SCUBA-2 selected sources which only partially overlap with SMGs. Better constraints on CO properties of a diverse high-z ULIRG population are needed from ALMA to determine the evolutionary origin of extreme starbursts, and what role ULIRGs serve in catalyzing the formation of massive stellar systems in the early Universe.

We identify and investigate the nature of the 20 brightest 250-micron sources detected by the BLAST within the central 150 sq. arcmin of the GOODS-South field. Aided by the available deep VLA radio imaging, reaching S_1.4 = 30 micro-Jy, we have identified radio counterparts for 17/20 of the 250-micron sources. The resulting enhanced positional accuracy of ~1 arcsec has then allowed us to exploit the deep multi-frequency imaging of GOODS-South to establish secure galaxy counterparts for the 17 radio-identified sources, and plausible galaxy candidates for the 3 radio-unidentified sources. Confusion is a serious issue for this deep BLAST 250-micron survey, due to the large size of the beam. Nevertheless, we argue that our chosen counterparts are significant, and often dominant contributors to the measured BLAST flux densities. For all of these 20 galaxies we have been able to determine spectroscopic (8) or photometric (12) redshifts. The result is the first `complete' redshift distribution for a deep 250-micron selected galaxy sample. This reveals that 250-micron surveys reaching detection limits of ~30 mJy contain not only low-redshift spirals/LIRGs, but also the extreme z~2 dust-enshrouded starburst galaxies previously discovered at sub-millimetre wavelengths. Based on their IRAC colours, we find that virtually all of the BLAST galaxy identifications appear better described as analogues of the M82 starburst galaxy, or Sc star-forming discs rather than highly obscured ULIRGs. Inspection of the LABOCA 870-micron imaging of the GOODS-South field yields detections of 7/11 of the z>1 BLAST sources, and reveals 250/870 flux-density ratios consistent with a standard 40K modified black-body fit with a dust emissivity index beta=1.5.

We present the first aperture synthesis unbiased spectral line survey toward an extragalactic object. The survey covered the 40 GHz frequency range between 202 and 242 GHz of the 1.3 mm atmospheric window. We find that 80% of the observed band shows molecular emission, with 73 features identified from 15 molecular species and 6 isotopologues. The 13C isotopic substitutions of HC3N and transitions from H2(18)O, 29SiO, and CH2CO are detected for the first time outside the Galaxy. Within the broad observed band, we estimate that 28% of the total measured flux is due to the molecular line contribution, with CO only contributing 9% to the overall flux. We present maps of the CO emission at a resolution of 2.9"x1.9" which, though not enough to resolve the two nuclei, recover all the single-dish flux. The 40 GHz spectral scan has been modelled assuming LTE conditions and abundances are derived for all identified species. The chemical composition of Arp 220 shows no clear evidence of an AGN impact on the molecular emission but seems indicative of a purely starburst-heated ISM. The overabundance of H2S and the low isotopic ratios observed suggest a chemically enriched environment by consecutive bursts of star formation, with an ongoing burst at an early evolutionary stage. The large abundance of water (~10^-5), derived from the isotopologue H2(18)O, as well as the vibrationally excited emission from HC3N and CH3CN are claimed to be evidence of massive star forming regions within Arp 220. Moreover, the observations put strong constraints on the compactness of the starburst event in Arp 220. We estimate that such emission would require ~2-8x10^6 hot cores, similar to those found in the Sgr B2 region in the Galactic center, concentrated within the central 700 pc of Arp 220.

We report the results of a pilot study with the EVLA of 12CO J=1-0 emission from four SMGs at z=2.2-2.5, each with an existing detection of CO J=3-2. Using the EVLA's most compact configuration we detect strong, broad J=1-0 line emission from all of our targets. The median line width ratio, sigma(1-0)/sigma(3-2) = 1.15 +/- 0.06, suggests that the J=1-0 is more spatially extended than the J=3-2 emission, a situation confirmed by our maps which reveal velocity structure in several cases and typical sizes of ~16 kpc FWHM. The median Tb ratio is r(3-2/1-0) = 0.55 +/- 0.05, noting that our value may be biased high because of the J=3-2-based sample selection. Naively, this suggests gas masses ~2x higher than estimates made using higher-J transitions of CO, with the discrepency due to the difference in assumed Tb ratio. We also estimate masses using the 12CO J=1-0 line and the observed global Tb ratios, assuming standard underlying Tb ratios as well as a limiting SFE, i.e. without calling upon X(CO). Using this new method, we find a median molecular gas mass of (2.5 +/- 0.8) x 10^10 Msun, with a plausible range stretching 3x higher. Even larger masses cannot be ruled out, but are not favoured by dynamical constraints: the median dynamical mass for our sample is (2.3 +/- 1.4) x 10^11 Msun. We examine the Schmidt-Kennicutt relation for all the distant galaxy populations for which CO J=1-0 or J=2-1 data are available, finding small systematic differences. These have previously been interpreted as evidence for different modes of star formation, but we argue that these differences are to be expected, given the still considerable uncertainties. Finally, we discuss the morass of degeneracies surrounding molecular gas mass estimates, the possibilities for breaking them, and the future prospects for imaging and studying cold, quiescent molecular gas at high redshifts [abridged].

We present a joint analysis of the overlapping BLAST 250, 350, 500um, and LABOCA 870um observations (from the LESS survey) of the ECDF-S. Out to z~3, the BLAST filters sample near the peak wavelength of far-infrared (FIR) emission from galaxies (rest-frame wavelengths ~60--200um), primarily produced by dust heated through absorption in star-forming clouds. However, identifying counterparts to individual BLAST peaks is very challenging, given the large beams (FWHM 36--60"). In contrast, the 870um observations have a significantly smaller 19" FWHM beam, and are sensitive to higher redshifts (z~1--5, and potentially beyond) due to the more favourable negative K-correction. We use the LESS data, as well as deep Spitzer and VLA imaging, to identify 118 individual sources that produce significant emission in the BLAST bands. We characterize the temperatures and FIR luminosities for a subset of 69 sources which have well-measured submm SEDs and redshift measurements out to z~3. For flux-limited sub-samples in each BLAST band, and a dust emissivity index \beta=2.0, we find a median temperature T=30K (all bands) as well as median redshifts: z=1.1 (interquartile range 0.2--1.9) for S250 > 40mJy; z=1.3 (interquartile range 0.6--2.1) for S350 > 30mJy; and z=1.6 (interquartile range 1.3--2.3) for S500 > 20mJy. Taking into account selection effects for our survey (a bias toward detecting lower-temperature galaxies), we find no evidence for evolution in the local FIR-temperature correlation out to z~2.5. Comparing with star-forming galaxy SED templates, about 8% of our sample appears to exhibit significant excesses in the radio and/or mid-IR, consistent with those sources harbouring an AGN. We describe the following techniques in two appendices: our `matched filter' for identifying sources in the presence of point-source confusion; and our approach for identifying counterparts using likelihood ratios.

We present a study of the cosmic infrared background, which is a measure of the dust obscured activity in all galaxies in the Universe. We venture to isolate the galaxies responsible for the background at 1mm; with spectroscopic and photometric redshifts we constrain the redshift distribution of these galaxies. We create a deep 1.16mm map (sigma ~ 0.5mJy) by combining the AzTEC 1.1mm and MAMBO 1.2mm datasets in GOODS-N. This combined map contains 41 secure detections, 13 of which are new. By averaging the 1.16mm flux densities of individually undetected galaxies with 24um flux densities > 25uJy, we resolve 31--45 per cent of the 1.16mm background. Repeating our analysis on the SCUBA 850um map, we resolve a higher percentage (40--64 per cent) of the 850um background. A majority of the background resolved (attributed to individual galaxies) at both wavelengths comes from galaxies at z > 1.3. If the ratio of the resolved submillimeter to millimeter background is applied to a reasonable scenario for the origins of the unresolved submillimeter background, 60--88 per cent of the total 1.16mm background comes from galaxies at z > 1.3.

We analyse subarcsecond resolution interferometric CO line data for twelve sub-millimetre-luminous (S850um > 5mJy) galaxies with redshifts between 1 and 3, presenting new data for four of them. Morphologically and kinematically most of the twelve systems appear to be major mergers. Five of them are well-resolved binary systems, and seven are compact or poorly resolved. Of the four binary systems for which mass measurements for both separate components can be made, all have mass ratios of 1:3 or closer. Furthermore, comparison of the ratio of compact to binary systems with that observed in local ULIRGs indicates that at least a significant fraction of the compact SMGs must also be late-stage mergers. In addition, the dynamical and gas masses we derive are most consistent with the lower end of the range of stellar masses published for these systems, favouring cosmological models in which SMGs result from mergers. These results all point to the same conclusion, that likely most of the bright SMGs with L_IR > 5x10e12L_sun are major mergers.

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.

We present SHARC-2 350 micron data on 20 luminous z ~ 2 starbursts with S(1.2mm) > 2 mJy from the Spitzer-selected samples of Lonsdale et al. and Fiolet et al. All the sources were detected, with S(350um) > 25 mJy for 18 of them. With the data, we determine precise dust temperatures and luminosities for these galaxies using both single-temperature fits and models with power-law mass--temperature distributions. We derive appropriate formulae to use when optical depths are non-negligible. Our models provide an excellent fit to the 6um--2mm measurements of local starbursts. We find characteristic single-component temperatures T1 ~ 35.5+-2.2 K and integrated infrared (IR) luminosities around 10^(12.9+-0.1) Lsun for the SWIRE-selected sources. Molecular gas masses are estimated at 4 x 10^(10) Msun, assuming kappa(850um)=0.15 m^2/kg and a submillimeter-selected galaxy (SMG)-like gas-to-dust mass ratio. The best-fit models imply >~2 kpc emission scales. We also note a tight correlation between rest-frame 1.4 GHz radio and IR luminosities confirming star formation as the predominant power source. The far-IR properties of our sample are indistinguishable from the purely submillimeter-selected populations from current surveys. We therefore conclude that our original selection criteria, based on mid-IR colors and 24 um flux densities, provides an effective means for the study of SMGs at z ~ 1.5--2.5.

The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

We present submillimetre and mid-infrared imaging observations of five fields centred on quasi-stellar objects (QSOs) at 1.7<z<2.8. All 5 QSOs were detected previously at submillimetre wavelengths. At 850 (450) um we detect 17 (11) submillimetre galaxies (SMGs) in addition to the QSOs. The total area mapped at 850 um is ~28 arcmin^2 down to RMS noise levels of 1-2 mJy/beam, depending on the field. Integral number counts are computed from the 850 um data using the same analytical techniques adopted by `blank-field' submillimetre surveys. We find that the `QSO-field' counts show a clear excess over the blank-field counts at deboosted flux densities of 2-4 mJy; at higher flux densities the counts are consistent with the blank-field counts. Robust mid-infrared counterparts are identified for all four submillimetre detected QSOs and ~60 per cent of the SMGs. The mid-infrared colours of the QSOs are similar to those of the local ULIRG/AGN Mrk 231 if placed at 1<z<3 whilst most of the SMGs have colours very similar to those of the local ULIRG Arp 220 at 1<z<3. Mid-infrared diagnostics therefore find no strong evidence that the SMGs host buried AGN although we cannot rule out such a possibility. Taken together our results suggest that the QSOs sit in regions of the early universe which are undergoing an enhanced level of major star-formation activity, and should evolve to become similarly dense regions containing massive galaxies at the present epoch. Finally, we find evidence that the level of star-formation activity in individual galaxies appears to be lower around the QSOs than it is around more powerful radio-loud AGN at higher redshifts.

We study the co-existence of star formation and AGN activity in X-ray selected AGN by analyzing stacked 870um submm emission from a deep and wide map of the Extended Chandra Deep Field South, obtained with LABOCA at the APEX telescope. The total X-ray sample of 895 sources with median redshift z~1 is detected at a mean submm flux of 0.49+-0.04mJy, corresponding to a typical star formation rate around 30Msun/yr for a T=35K, beta=1.5 greybody far-infrared SED. The good S/N permits stacking analyses for subgroups. We observe a trend of star formation rate increasing with redshift. An increase of star formation rate with AGN luminosity is indicated at the highest L_2-10>~1E44erg/s luminosities only. Increasing trends with X-ray obscuration as expected in some AGN evolutionary scenarios are not observed for the bulk of the X-ray AGN sample but may be present for the highest intrinsic luminosity objects. This suggests a transition between two modes in the coexistence of AGN activity and star formation. For the bulk of the sample, the X-ray luminosity and obscuration of the AGN are not intimately linked to the global star formation rate of their hosts. The hosts are likely massive and forming stars secularly, at rates similar to the pervasive star formation seen in massive galaxies without an AGN at similar redshifts. The change indicated towards more intense star formation, and a more pronounced increase in star formation rates between unobscured and obscured AGN at highest luminosities suggests that luminous AGN follow an evolutionary path on which obscured AGN activity and intense star formation are linked, possibly via merging. Comparison to local hard X-ray selected AGN supports this interpretation. [Abridged]

We present high spatial resolution (0.4", ~3.5 kpc) PdBI interferometric data on three ultra-luminous infrared galaxies (ULIRGs) at z~2: two submillimetre galaxies and one submillimetre faint star forming radio galaxy. The three galaxies have been ro- bustly detected in CO rotational transitions, either 12CO(J=4-3) or 12CO(J=3-2), allowing their sizes and gas masses to be accurately constrained. These are the highest spatial resolution observations observed to date (by a factor of ~2) for intermediate-excitation CO emission in z~2 ULIRGs. The galaxies appear extended over several resolution elements, having a mean radius of 3.7 kpc. High-resolution (0.3") combined MERLIN-VLA observations of their radio continua allow an analysis of the star formation behaviour of these galaxies, on comparable spatial scales to that of the CO observations. This 'matched beam' approach sheds light on the spatial distribution of both molecular gas and star formation, and we can therefore calculate accurate star formation rates and gas surface densities: this allows us to place the three systems in the context of a Kennicutt-Schmidt (KS)-style star formation law. We find a difference in size between the CO and radio emission regions, and as such we suggest that using the spatial extent of the CO emission region to estimate the surface density of star formation may lead to error. This size difference also causes the star formation efficiencies within systems to vary by up to a factor of 5. We also find, with our new accurate sizes, that SMGs lie significantly above the KS relation, indicating that stars are formed more efficiently in these extreme systems than in other high-z star forming galaxies.

Over the last decade observations at submillimetre (submm) and millimetre (mm) wavelengths, with their unique ability to trace molecular gas and dust, have attained a central role in our exploration of galaxies at all redshifts. Due to the limited sensitivities and angular resolutions of current submm/mm telescopes, however, only the most luminous objects have been uncovered at high redshifts, with interferometric follow-up observations succeeding in resolving the dust and gas reservoirs in only a handful of cases. The coming years will witness a drastic improvement in the current situation, thanks to the arrival of a new suite of powerful submm observatories (single-dish and interferometers) with an order of magnitude improvement in sensitivity and resolution. In this overview I outline a few of what I expect to be the major advances in the field of galaxy formation and evolution that these new ground-breaking facilities will facilitate.

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 IRAM Plateau-de-Bure Interferometer (PdBI) we have searched for the upper fine structure line of neutral carbon (CI(2->1), nu_rest = 809 GHz) and 12CO(J=7->6) (nu_rest=806GHz) towards the submillimetre galaxies (SMGs) GN20 (SMM J123711.9+622212, z = 4.055) and GN20.2 (SMM J123708.8+622202, z = 4.051). The far-infrared (FIR) continuum is detected at 8sigma significance in GN20, with a flux density of S_1.8mm = 1.9+-0.2mJy, while no continuum is detected in GN20.2. Both sources are statistically undetected in both CI and CO(7-6) lines; we derive line luminosity limits for both CI and CO of L'<~2x10^10K km s^-1 pc^2. Assuming carbon excitation temperatures of T_ex = 30K (the galaxies' measured dust temperatures), we infer CI mass limits of M_CI < 5.4x10^6 Msun (GN20) and M_CI < 6.8x10^6 Msun (GN20.2). The derived CI abundance limits are <1.8x10^-5 for GN20 and <3.8x10^-5 for GN20.2 implying that the systems have Milky Way level neutral carbon enrichment (X[CI]/X[H2]) or lower, similar to high-redshift carbon-detected systems (at 5x10^-5) but about 50 times less than the neutral carbon enrichment of local starburst galaxies. Observations of GN20 and GN20.2 in high-resolution MERLIN+VLA radio maps of GOODS-N are used to further constrain the sizes and locations of active regions. We conclude that the physical gas properties of young rapidly evolving systems like GN20 and GN20.2 are likely significantly different than starburst/ULIRG environments in the local Universe yet similar to z~2 SMGs. Unless gravitationally amplified examples can be found, observations of galaxies like GN20 will require the order of magnitude increase in sensitivity of the Atacama Large Millimetre Array (ALMA) to constrain their CI and high-J CO content, despite the fact that they are the brightest systems at z~4.

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.

The interstellar medium (ISM) is a key ingredient in galaxy formation and evolution as it provides the molecular gas reservoir which fuels star formation and supermassive black hole accretion. Yet the ISM is one of the least studied aspects of distant galaxies. Molecular and atomic transitions at (sub)millimetre wavelengths hold great promise in measuring macroscopic properties (e.g. masses, morphologies, star formation laws), as well as microscopic properties (e.g. gas densities, temperatures, cooling) of high-z galaxies. In this overview I summarize the growing number of high-z molecular line detections, highlighting some of the most intriguing results along the way. I end by discussing a few areas where future facilities (e.g. ALMA, EVLA, CCAT, LMT) will drastically improve on the current state of affairs.

We report single-dish multi-transition measurements of the 12^CO, HCN, and HCO^+ molecular line emission as well as HNC J=1-0 and HNCO in the two ultraluminous infra-red galaxies Arp220 and NGC6240. Using this new molecular line inventory, in conjunction with existing data in the literature, we compiled the most extensive molecular line data sets to date for such galaxies. The many rotational transitions, with their different excitation requirements, allow the study of the molecular gas over a wide range of different densities and temperatures with significant redundancy, and thus allow good constraints on the properties of the dense gas in these two systems. The mass (~(1-2) x 10^10 Msun) of dense gas (>10^5-6 cm^-3) found accounts for the bulk of their molecular gas mass, and is consistent with most of their IR luminosities powered by intense star bursts while self-regulated by O,B star cluster radiative pressure onto the star-forming dense molecular gas. The highly excited HCN transitions trace a gas phase ~(10-100)x denser than that of the sub-thermally excited HCO^+ lines (for both galaxies). These two phases are consistent with an underlying density-size power law found for Galactic GMCs (but with a steeper exponent), with HCN lines tracing denser and more compact regions than HCO^+. Whether this is true in IR-luminous, star forming galaxies in general remains to be seen, and underlines the need for observations of molecular transitions with high critical densities for a sample of bright (U)LIRGs in the local Universe -- a task for which the HI-FI instrument on board Herschel is ideally suited to do.

We use deep far-IR, submm, radio and X-ray imaging and mid-IR spectroscopy to explore the nature of a sample of Spitzer-selected dust-obscured galaxies (DOGs) in GOODS-N. A sample of 79 galaxies satisfy the criteria R-[24]>14 (Vega) down to S24>100 microJy. Twelve of these galaxies have IRS spectra available which we use to measure redshifts and classify these objects as being dominated by star formation or active galactic nuclei (AGN) activity in the mid-IR. The IRS spectra and Spitzer photometric redshifts confirm that the DOGs lie in a tight redshift distribution around z~2. Based on mid-IR colors, 80% of DOGs are likely dominated by star formation; the stacked X-ray emission from this sub-sample of DOGs is also consistent with star formation. Since only a small number of DOGs are individually detected at far-IR and submm wavelengths, we use a stacking analysis to determine the average flux from these objects and plot a composite IR (8-1000 microns) spectral energy distribution (SED). The average luminosity of these star forming DOGs is LIR~1e12 Lsun. We compare the average star forming DOG to the average bright (S850>5 mJy) submillimeter galaxy (SMG); the S24>100 microJy DOGs are 3 times more numerous but 8 times less luminous in the IR. The far-IR SED shape of DOGs is similar to that of SMGs (average dust temperature of around 30 K) but DOGs have a higher mid-IR to far-IR flux ratio. The average star formation-dominated DOG has a star formation rate of 200 Msun/yr which, given their space density, amounts to a contribution of 0.01 Msun/yr/Mpc3 (or 5-10%) to the star formation rate density at z~2. We use the composite SED to predict the average flux of DOGs in future Herschel/PACS 100 micron and SCUBA-2 450 micron surveys and show that the majority of them will be detected.

High-resolution submm imaging of the HzRG, 4C60.07, at z=3.8, has revealed two dusty components. Spitzer imaging shows that one of these components (B) is coincident with an extremely red AGN, offset by ~4" (~30 kpc) from the HzRG core. The other submm component (A) - resolved by our beam and devoid of emission at 3.6-8.0um - lies between B and the HzRG core. Since the HzRG was discovered via its young, steep-spectrum lobes and their creation was likely triggered by the interaction, we argue that we are witnessing an early-stage merger, prior to its eventual equilibrium state. The interaction is between the host galaxy of an actively-fueled BH, and a gas-rich starburst/AGN (B) marked by the compact submm component and coincident with broad CO emission. `A' is a plume of cold, dusty gas, associated with a narrow (~150 km/s) CO feature, and may represent a short-lived tidal structure. It has been claimed that HzRGs and SMGs differ only in the activity of their AGNs, but such complex submm morphologies are seen only rarely amongst SMGs. Our study has important implications: where a galaxy's gas is not aligned with its central BH, CO may be an unreliable probe of dynamical mass, affecting work on the co-assembly of BHs and spheroids. Our data support the picture wherein close binary AGN are induced by mergers. They also raise the possibility that some supposedly jet-induced starbursts may have formed co-evally with (yet independently of) the radio jets, both triggered by the same interaction. We note that the HzRG host would have gone unnoticed without its jets/companion, so there may be many other unseen BHs at high redshift, lost in the sea of ~5 x 10^8 similarly bright IRAC sources - sufficiently massive to drive a >10^27-W radio source, yet practically invisible unless actively fueled (abridged).

High-redshift, dust-obscured galaxies -- selected to be luminous in the radio but relatively faint at 850um -- appear to represent a different population from the ultra-luminous submillimeter- (submm-) bright population. They may be star-forming galaxies with hotter dust temperatures or they may have lower far-infrared luminosities and larger contributions from obscured active galactic nuclei (AGN). Here we present observations of three z~2 examples of this population, which we term submm-faint radio galaxies (SFRGs) in CO(3-2) using the IRAM Plateau de Bure Interferometer to study their gas and dynamical properties. We estimate the molecular gas mass in each of the three SFRGs (8.3x10^{9} M_odot, <5.6x10^{9} M_odot and 15.4x10^{9} M_odot, respectively) and, in the case of RG163655, a dynamical mass by measurement of the width of the CO(3-2) line (8x10^{10} csc^2i M_odot). While these gas masses are substantial, on average they are 4x lower than submm-selected galaxies (SMGs). Radio-inferred star formation rates (<SFR_radio>=970 M_odot\yr) suggest much higher star-formation efficiencies than are found for SMGs, and shorter gas depletion time scales (~11 Myr), much shorter than the time required to form their current stellar masses (~160 Myr; ~10^{11} M_odot). By contrast, SFRs may be overestimated by factors of a few, bringing the efficiencies in line with those typically measured for other ultraluminous star-forming galaxies and suggesting SFRGs are more like ultraviolet- (UV-)selected star-forming galaxies with enhanced radio emission. A tentative detection of \rga at 350um suggests hotter dust temperatures -- and thus similar gas-to-dust mass fractions -- as the SMGs. We conclude that SFRGs' radio luminosities are larger than would naturally scale from local ULIRGs given their gas masses or gas fractions.

We present a 1200-micron image of the Great Observatories Origin Deep Survey North (GOODS-N) field, obtained with the Max Planck Millimeter Bolometer array (MAMBO) on the IRAM 30-m telescope. The survey covers a contiguous area of 287 square arcmin to a near-uniform noise level of ~0.7mJy/beam. After Bayesian flux deboosting, a total of 30 sources are recovered (>=3.5sigma). An optimal combination of our 1200-micron data and an existing 850-micron image from the Submillimetre Common-User Bolometer Array (SCUBA) yielded 33 sources (>=4sigma). We combine our GOODS-N sample with those obtained in the Lockman Hole and ELAIS-N2 fields (Scott et al. 2002; Greve et al. 2004) in order to explore the degree of overlap between 1200-micron- and 850-micron-selected galaxies (hereafter SMGs), finding no significant difference between their 850-micron to 1200-micron flux density distributions. However, a noise-weighted stacking analysis yields a significant detection of the 1200-micron-blank SCUBA sources, whereas no significant 850-micron signal is found for the 850-micron-blank MAMBO sources. The hypothesis that the 850/1200-micron flux density distribution of SCUBA sources is also representative of the MAMBO population is rejected at the ~4sigma level, via Monte Carlo simulations. Therefore, although the populations overlap, galaxies selected at 850 and 1200micron are different, and there is compelling evidence for a significant 1200-micron-detected population which is not recovered at 850micron. These are submm drop-outs (SDOs), with S_850/S_1200 = 0.7-1.7, requiring very cold dust or unusual spectral energy distributions (T_d ~ 10K; beta ~ 1), unless SDOs reside beyond the redshift range observed for radio-identified SMGs, i.e. at z > 4.

We present high-resolution interferometric imaging of LH850.02, the brightest 850- and 1200-micron submillimetre (submm) galaxy in the Lockman Hole. Our observations were made at 890 micron with the Submillimetre Array (SMA). Our high-resolution submm imaging detects LH850.02 at >6-sigma as a single compact (size < 1 arcsec or < 8 kpc) point source and yields its absolute position to ~0.2-arcsec accuracy. LH850.02 has two alternative radio counterparts within the SCUBA beam (LH850.02N & S), both of which are statistically very unlikely to be so close to the SCUBA source position by chance. However, the precise astrometry from the SMA shows that the submm emission arises entirely from LH850.02N, and is not associated with LH850.02S (by far the brighter of the two alternative identifications at 24-micron). Fits to the optical-infrared multi-colour photometry of LH850.02N & S indicate that both lie at z~3.3, and are therefore likely to be physically associated. At these redshifts, the 24 micron--to--submm flux density ratios suggest that LH850.02N has an Arp220-type starburst-dominated far-IR SED, while LH850.02S is more similar to Mrk231, with less dust-enshrouded star-formation activity, but a significant contribution at 24-micron (rest-frame ~5-6 micron) from an active nucleus. This complex mix of star-formation and AGN activity in multi-component sources may be common in the high redshift ultraluminous galaxy population, and highlights the need for precise astrometry from high resolution interferometric imaging for a more complete understanding.

[Abridged] The detection of the rotational lines of CO in proto-galaxies in the early Universe provides one of the most promising ways of probing the fundamental physical properties of a galaxy, such as its size, dynamical mass, gas density, and temperature. While such observations are currently limited to the most luminous galaxies, the advent of ALMA will change the situation dramatically, resulting in the detection of numerous normal galaxies at high redshifts. Maps and spectra of rotational CO line emission were calculated from a cosmological N-body/hydrodynamical TreeSPH simulation of a z ~ 3 "Lyman break galaxy'' of UV luminosity about one order of magnitude below L*. To simulate a typical observation of our system with ALMA, we imposed characteristic noise, angular, and spectral resolution constraints. The CO line properties predicted by our simulation are in good agreement with the two Lyman break systems detected in CO to date. We find that while supernovae explosions from the ongoing star formation carve out large cavities in the molecular ISM, they do not generate large enough gas outflows to make a substantial imprint on the CO line profile. This implies that for most proto-galaxies - except possibly the most extreme cases - stellar feedback effects do not affect CO as a dynamical mass tracer. Detecting CO in sub-L* galaxies at z ~3 will push ALMA to the limits of its cababilities, and whether a source is detected or not may depend critically on its inclination angle. Both these effects (sensitivity and inclination) will severely impair the ability of ALMA to infer the gas kinematics and dynamical masses using line observations.

We present the optical-to-submm spectral energy distributions for 33 radio & mid-IR identified submillimetre galaxies discovered via the SHADES 850 micron SCUBA imaging in the Subaru-XMM Deep Field (SXDF). Optical data for the sources comes from the Subaru-XMM Deep Field (SXDF) and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6 and 4.5 micron fluxes. We then fit spectral energy distribution (SED) templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The infrared template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submm sources and on the SHADES SXDF field. Our template fitting finds that AGN, present in about 10% of sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fit by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submm galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.

We report sub-arcsecond resolution IRAM PdBI millimeter CO interferometry of four z~2 submillimeter galaxies (SMGs), and sensitive CO (3-2) flux limits toward three z~2 UV-/optically selected star forming galaxies. The new data reveal for the first time spatially resolved CO gas kinematics in the observed SMGs. Two of the SMGs show double or multiple morphologies, with complex, disturbed gas motions. The other two SMGs exhibit CO velocity gradients of ~500 km/s across 0.2 arcsec (1.6 kpc) diameter regions, suggesting that the star forming gas is in compact, rotating disks. Our data provide compelling evidence that these SMGs represent extreme, short-lived 'maximum' star forming events in highly dissipative mergers of gas rich galaxies. The resulting high mass surface and volume densities of SMGs are similar to those of compact quiescent galaxies in the same redshift range, and much higher than those in local spheroids. From the ratio of the comoving volume densities of SMGs and quiescent galaxies in the same mass and redshift ranges, and from the comparison of gas exhaustion time scales and stellar ages, we estimate that the SMG phase duration is about 100 Myrs. Our analysis of SMGs and optically/UV selected high redshift star forming galaxies supports a 'universal' Chabrier IMF as being valid over the star forming history of these galaxies. We find that the 12CO luminosity to total gas mass conversion factors at z~2-3 are probably similar to those assumed at z~0. The implied gas fractions in our sample galaxies range from 20 to 50%.

We present deep 350- and 1200-micron imaging of the region around 4C41.17 -- one of the most distant (z = 3.792) and luminous known radio galaxies -- obtained with the Submillimeter High Angular Resolution Camera (SHARC-II) and the Max Planck Millimeter Bolometer Array (MAMBO). The radio galaxy is robustly detected at 350- and 1200-micron, as are two nearby 850-micron-selected galaxies; a third 850-micron source is detected at 350-micron and coincides with a ~ 2-sigma feature in the 1200-micron map. Further away from the radio galaxy an additional nine sources are detected at 1200-micron, bringing the total number of detected (sub)millimeter selected galaxies (SMGs) in this field to 14. Using radio images from the Very Large Array (VLA) and Spitzer mid-infrared (mid-IR) data, we find statistically robust radio and/or 24-micron counterparts to eight of the 14 SMGs in the field around 4C41.17. Follow-up spectroscopy with Keck/LRIS has yielded redshifts for three of the eight robustly identified SMGs, placing them in the redshift range 0.5 < z < 2.7, i.e. well below that of 4C41.17. We infer photometric redshifts for a further four sources using their 1.6-micron (rest-frame) stellar feature as probed by the IRAC bands; only one of them is likely to be at the same redshift as 4C41.17. Thus at least four, and as many as seven, of the SMGs within the 4C41.17 field are physically unrelated to the radio galaxy. With the redshift information at hand we are able to constrain the observed over-densities of SMGs within radial bins stretching to R=50 and 100" (~ 0.4 and ~ 0.8Mpc at z ~ 3.8) from the radio galaxy to ~ 5x and ~ 2x that of the field, dropping off to the background value at R=150". [Abridged]

Determining an accurate position for a submm galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample - source counts and 2-D clustering - to an assessment of their detailed, multi-wavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or IR counterparts, and hence accurate positions, for over two thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-um and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parameterisation appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour-colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2-6", or 15-50/(sin i) kpc at z ~ 2, consistent with early bursts seen in merger simulations.

We present 350micron observations of 36 ultraluminous infrared galaxies (ULIRGs) at intermediate redshifts (0.089 <= z <= 0.926) using the Submillimeter High Angular Resolution Camera II (SHARC-II) on the Caltech Submillimeter Observatory (CSO). In total, 28 sources are detected at S/N >= 3, providing the first flux measurements longward of 100micron for a statistically significant sample of ULIRGs in the redshift range of 0.1 < z < 1.0. Combining our 350micron flux measurements with the existing IRAS 60 and 100micron data, we fit a single-temperature model to the spectral energy distribution (SED), and thereby estimate dust temperatures and far-IR luminosities. Assuming an emissivity index of beta = 1.5, we find a median dust temperature and far-IR luminosity of Td = 42.8+-7.1K and log(Lfir/Lsolar) = 12.2+-0.5, respectively. The far-IR/radio correlation observed in local star-forming galaxies is found to hold for ULIRGs in the redshift range 0.1 < z < 0.5, suggesting that the dust in these sources is predominantly heated by starbursts. We compare the far-IR luminosities and dust temperatures derived for dusty galaxy samples at low and high redshifts with our sample of ULIRGs at intermediate redshift. A general Lfir-Td relation is observed, albeit with significant scatter, due to differing selection effects and variations in dust mass and grain properties. The relatively high dust temperatures observed for our sample compared to that of high-z submillimeter-selected starbursts with similar far-IR luminosities suggest that the dominant star formation in ULIRGs at moderate redshifts takes place on smaller spatial scales than at higher redshifts.

A large CO, HCN multi-transition survey of 30 Luminous Infrared Galaxies ($\rm L_{IR}>10^{11} L_{\odot}$) is nearing completion with the James Clerk Maxwell Telescope (JCMT) on Mauna Kea (Hawaii), and the IRAM 30-meter telescope at Pico Veleta (Spain). The CO J=1--0, 2--1, 3--2, 4--3,6--5, $ ^{13}$CO J=2--1, HCN J=1--0, 3--2, 4--3 observations, resulting from $\sim 250$ hours of JCMT, $\sim 100$ hours of 30-m observing time and data from the literature constitute {\it the largest extragalactic molecular line survey to date}, and can be used to address a wide range of issues and eventually allow the construction of reliable Spectral Line Energy Distributions (SLEDs) for the molecular gas in local starbursts. First results suggest that: a) HCN and HCO$^+$ J=1--0 line luminosities can be poor mass estimators of dense molecular gas ($\rm n\geq 10^4 cm^{-3}$) unless their excitation is accounted for, b) CO cooling of such gas in ULIRGs may be comparable to that of the CII line at $\rm 158 \mu m$, and c) low excitation of the {\it global} molecular gas reservoir remains possible in such systems. In such cases the expected low CO $\rm J+1\to J$ line luminosities for $\rm J+1\geq 4$ can lead to a strong bias against their detection from ULIRGs at high redshifts.

Ultra Luminous Infra-Red Galaxies (ULIRGs) -- gas-rich mergers harboring the most extreme star-forming conditions encountered in the local Universe -- are thought to re-enact the galaxy formation processes we are only barely able to glimpse in the distant Universe. Combining new single-dish molecular line observations of 12CO, 13CO, HCO+, HCN, and CS towards the two ULIRGs Arp220 and NGC6240 with existing data in the literature, we have compiled the most extensive molecular line data-sets to date of such galaxies. The data allow us to put strong constraints on the properties of the dense star forming gas in these two systems, and compare the relative usefulness of CS, HCN and HCO+ as tracers of dense gas. In addition, we have build molecular line templates based on our observations, and demonstrate that Herschel/HI-FI will be able to detect the high-J transitions of most of the above molecules in a large sample of ULIRGs out to z<=0.5, assuming Arp220 and NGC6240 are representative of the ULIRG population at these redshifts.

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

We report deep K-band (18-27GHz) observations with the 100-m Green Bank Telescope of HCN(1-0) line emission towards the two submillimeter-selected galaxies (SMGs) SMMJ02399-0136 (z=2.81) and SMMJ16359+6612 (z=2.52). For both sources we have obtained spectra with channel-to-channel rms noise of <=0.5mJy, resulting in velocity-integrated line fluxes better than < 0.1 Jy km/s, although we do not detect either source. Such sensitive observations -- aided by gravitational lensing of the sources -- permit us to put upper limits of L_HCN(1-0) < 2x10^10 K km/s pc^2 on the intrinsic HCN(1-0) line luminosities of the two SMGs. The far-infrared (FIR) luminosities for all three SMGs with sensitive HCN(1-0) observations to date are found to be consistent with the tight FIR-HCN luminosity correlation observed in Galactic molecular clouds, quiescent spirals and (ultra) luminous infrared galaxies in the local Universe. Thus, the observed HCN luminosities remain in accordance with what is expected from the universal star formation efficiency per dense molecular gas mass implied by the aforementioned correlation, and more sensitive observations with today's large aperture radio telescopes hold the promise of detecting HCN(1-0) emission in similar objects in the distant Universe.

We report the first detection of CO(1-0) emission from a submillimeter-selected galaxy, using the Green Bank Telescope. We identify the line in the spectrum of SMM J13120+4242 as a broad emission feature at z=3.408, with Delta(V_FWHM)=1040 +/- 190 km/s. If the observed CO(1-0) line profile arises from a single object and not several merging objects, then the CO(4-3)/CO(1-0) brightness temperature ratio of ~0.26 suggests n(H_2) > 3-10 x 10^2 cm^-3 and the presence of sub-thermally excited gas. The integrated line flux implies a cold molecular gas mass M(H_2)=1.6 x 10^11 M_sun, comparable to the dynamical mass estimate and 4 times larger than the H_2 mass predicted from the CO(4-3) line assuming a brightness temperature ratio of 1.0. While our observations confirm that this submillimeter galaxy is massive and gas-rich, they also suggest that extrapolating gas masses from J_upper >= 3 transitions of CO leads to considerable uncertainties. We also report an upper limit to the mass of cold molecular gas in a second submillimeter galaxy, SMM J09431+4700, of M(H_2)< 4 x 10^10 M_sun.

We present sub-arcsecond resolution IRAM PdBI interferometry of eight submillimeter galaxies at redshifts from 2 to 3.4, where we detect continuum at 1mm and/or CO lines at 3 and 1 mm. The CO 3-2/4-3 line profiles in five of the sources are double-peaked, indicative of orbital motion either in a single rotating disk or of a merger of two galaxies. The millimeter line and continuum emission is compact; we marginally resolve the sources or obtain tight upper limits to their intrinsic sizes in all cases. The median FWHM diameter for these sources and the previously resolved sources, SMMJ023952-0136 and SMMJ140104+0252 is less than or equal to 0.5" (4 kpc). The compactness of the sources does not support a scenario where the far-IR/submm emission comes from a cold, very extended dust distribution. These measurements clearly show that the submillimeter galaxies we have observed resemble scaled-up and more gas rich versions of the local Universe, ultra-luminous galaxy (ULIRG) population. Their central densities and potential well depths are much greater than in other redshift 2-3 galaxy samples studied so far. They are comparable to those of elliptical galaxies or massive bulges. The SMG properties fulfill the criteria of 'maximal' starbursts, in which most of the available initial gas reservoir of 10^10-10^11 solar masses is converted to stars on a few dynamical timescales.

The modest significance of most sources detected in current submm surveys can potentially compromise some analyses due to the inclusion of spurious sources in catalogues typically selected at >3.0-3.5sigma. Here, we develop and apply a dual-survey extraction technique to SCUBA and MAMBO images of the Lockman Hole. Cut above 5sigma, our catalogue of SMGs is more robust than previous samples, with a reduced likelihood of real, but faint SMGs (beneath and around the confusion limit) entering via superposition with noise. The effective flux limit of the survey is well matched to our deep 1.4-GHz image. The former is sensitive to luminous, dusty galaxies at extreme redshifts whilst the latter probes the z<3 regime. A high fraction of our robust SMGs (~80%) have radio counterparts which, given the ~10% contamination by spurious sources, suggests that very distant SMGs (z>>3) are unlikely to make up more than ~10% of the bright SMG population. This implies that almost all of the bright SMG population is amenable to study via the deepest current radio imaging. We use these radio counterparts to provide an empirical calibration of the positional uncertainty in SMG catalogues. We then go on to outline the acquisition of redshifts for radio-identified SMGs, from sample selection in the submm, to counterpart selection in the radio and optical/IR, to slit placement on spectrograph masks. We determine a median of z=2.05 +/- 0.41 from a sample of six secure redshifts for unambigious radio-identified submm sources and z=2.14 +/- 0.27 when we include submm sources with multiple radio counterparts and/or less reliable redshifts. These figures are consistent with previous estimates, suggesting that our knowledge of the median redshift of bright SMGs population has not been biased by the low significance of the source catalogues employed.

Using the filled bolometer array SHARC-II on the Caltech Submillimeter Observatory (CSO), we have obtained deep (rms ~ 15mJy/beam), high-quality 350 micron maps of five of the most luminous high-z radio galaxies known. In all cases the central radio galaxy is detected at the >3sigma level, and in some cases the high resolution of SHARC-II (FWHM ~ 9") allows us to confirm the spatially extended nature of the dust emission. In PKS1138-262 (z=2.156), 8C1909+722 (z=3.538) and 4C41.17 (z=3.792), additional sources -- first discovered by SCUBA at 850 micron and believed to be dusty, merging systems associated with the central radio galaxy -- are detected at 350 micron. Furthermore, in PKS1138 and 4C41.17 additional SHARC-II sources are seen which were not detected at 850 micron, although the reality of these sources will have to be confirmed by independent submm observations. Thus, our observations seem to support the notion of extended star formation taking place in radio galaxies at high redshifts, and that these systems sit at the centers of overdense regions harbouring a complexity of dusty and vigorously star forming systems. At the redshift of the radio galaxies, the 350 micron observations sample very close to the rest-frame dust peak (typically at ~100 micron), and they therefore contribute a particularly important point to the spectral energy distributions of these sources, which we use this in conjunction with existing (sub)millimeter data to derive FIR luminosities, dust temperatures and spectral emissivities of the central radio galaxies.

We report sensitive imaging observations of the CO J=1-0 line emission in the powerful high-redshift radio galaxy 4C 41.17 (z=3.8) with the NRAO Very Large Array (VLA), conducted in order to detect the large concomitant H_2 gas reservoir recently unveiled in this system by De Breuck et al (2005) via the emission of the high excitation J=4-3 line. Our observations fail to detect the J=1-0 line but yield sensitive lower limits on the R_43=(4-3)/(1-0) brightness temperature ratio of R_43 ~ 0.55 - >1.0 for the bulk of the H_2 gas mass. Such high ratios are typical of the high-excitation molecular gas phase ``fueling'' the star formation in local starbursts, but quite unlike these objects, much of the molecular gas in 4C 41.17 seems to be in such a state, and thus participating in the observed starburst episode. The widely observed and unique association of highly excited molecular gas with star forming sites allows CO line emission with large (high-J)/(low-J) intensity ratios to serve as an excellent ``marker'' of the spatial distribution of star formation in distant dust-obscured starbursts, unaffected by extinction.

In this paper we present results from an IRAM Plateau de Bure millimetre-wave Interferometer (PdBI) survey for CO emission towards radio-detected submillimetre galaxies (SMGs) with known optical and near-infrared spectroscopic redshifts. Five sources in the redshift range z~1-3.5 were detected, nearly doubling the number of SMGs detected in CO. We summarise the properties of all 12 CO-detected SMGs, as well as 6 sources not detected in CO by our survey, and use this sample to explore the bulk physical properties of the SMG population as a whole. The median CO line luminosity of the SMGs is <L_CO> = (3.8 +- 2.0) x 10^10 K km/s pc^2. Using a CO-to-H_2 conversion factor appropriate for starburst galaxies, this corresponds to a molecular gas mass <M(H_2)> = (3.0 +- 1.6) x 10^10 Msun within a ~2kpc radius, about four times greater than the most luminous local ultraluminous infrared galaxies (ULIRGs) but comparable to that of the most extreme high-redshift radio galaxies and QSOs. The median CO fwhm linewidth is broad, <fwhm> = 780 +- 320 km/s, and the SMGs often have double peaked line profiles, indicative of either a merger or a disk. From their median gas reservoirs (~3 x 10^10 Msun) and star-formation rates (>700 Msun/yr) we estimate a lower limit on the typical gas-depletion time scale of >40Myr in SMGs. This is marginally below the typical age expected for the starbursts in SMGs, and suggests that negative feedback processes may play an important role in prolonging the gas consumption time scale. We find a statistically-significant correlation between the far-infrared and CO luminosities of the SMGs which extends the observed correlation for local ULIRGs to higher luminosities and higher redshifts. [ABRIDGED]

We present new sensitive wide-band measurements of the fine structure line 3^P_1 -> 3^P_0 (J=1-0, 492GHz) of neutral atomic carbon (CI) in the two typical Ultra Luminous Infrared Galaxies NGC6240 and Arp220. We then use them along with several other CI measurements in similar objects found in the literature to estimate their global molecular gas content under the assumption of a full CI-H_2 concomitance. We find excellent agreement between the H_2 gas mass estimated with this method and the standard methods using 12^CO. This may provide a new way to measure H_2 gas mass in galaxies, and one which may be very valuable in ULIRGs since in such systems the bright 12^CO emission is known to systematically overestimate the gas mass while their 13^CO emission is usually very weak. At redshifts z>=1 the CI J=1-0 line shifts to much more favorable atmospheric windows and can become a viable alternative tracer of the H_2 gas fueling starburst events in the distant Universe.

We present Spitzer observations in five wavebands between 3.6 and 24um of an unbiased sample of 9 luminous, dusty galaxies selected at 1200um by the MAMBO camera on the IRAM 30-m telescope, a population akin to the well-known submm or `SCUBA' galaxies (hereafter SMGs). Owing to the coarse resolution of submm/mm instrumentation, SMGs have traditionally been difficult to identify at other wavelengths. We compare our multi-wavelength catalogs to show that the overlap between 24 and 1200um must be close to complete at these flux levels. We find that all (4/4) of the most secure >=4sigma SMGs have robust >=4sigma counterparts at 1.4GHz, while the fraction drops to 7/9 using all >=3sigma SMGs. We show that combining mid-IR and marginal (>=3sigma) radio detections provides plausible identifications in the remaining cases, enabling us to identify the complete sample. Accretion onto an obscured central engine is betrayed by the shape of the mid-IR continuum emission for several sources, confirming Spitzer's potential to weed out active galaxies. We demonstrate the power of a S(24um)/S(8um) vs S(8um)/S(4.5um) color-color plot as a diagnostic for this purpose. However, we conclude that the majority (~75%) of SMGs have rest-frame mid-/far-IR SEDs commensurate with obscured starbursts. Sensitive 24-um observations are clearly a useful route to identify and characterize reliable counterparts to high-redshift far-IR-bright galaxies, complementing what is possible via deep radio imaging.

We present sub-millimeter statistical detections of galaxies discovered in the 5'x5' Spitzer Early Release Observations (to 4-15 microJy 5 sigma at 3.6-8 microns, 170 microJy at 24 microns) through a stacking analysis of our reanalysed SCUBA 8mJy survey maps, and a Spitzer identification of a new sub-millimeter point source in the 8mJy survey region. For sources detected at 5.8 or 8 microns (154 and 111 sources respectively), we detect positive skews in the sub-millimeter flux distributions at 99.2-99.8% confidence using Kolmogorov-Smirnov tests, at both 850 microns and 450 microns. We also marginally detect the Spitzer 24 micron galaxies at 850 microns at 97% confidence, and place limits on the mean sub-millimeter fluxes of the 3.6 and 4.5 micron sources. Integrating the sub-millimeter fluxes of the Spitzer populations, we find the 5.8 micron galaxies contribute 0.12 +/- 0.05 nW/m^2/sr to the 850 micron background, and 2.4 +/- 0.7 nW/m^2/sr to the 450 micron background; similar contributions are made by the 8 micron-selected sample. We infer that the populations dominating the 5.8 and 8 micron extragalactic background light also contribute around a quarter of the 850 micron background and the majority of the 450 micron background.

We present a deep, new 1200 micron survey of the ELAIS N2 and Lockman Hole fields using the Max Planck Millimeter Bolometer array (MAMBO). The areas surveyed are 160 arcmin^2 in ELAIS N2 and 197 arcmin^2 in the Lockman Hole, covering the entire SCUBA `8 mJy Survey'. In total, 27 (44) sources have been detected at a significance >=4.0sigma (>=3.5sigma). We present the 1200 micron number counts and find evidence of a fall at bright flux levels. Employing parametric models for the evolution of the local 60 micron IRAS luminosity function (LF), we are able to account simultaneously for the 1200 and 850 micron counts, suggesting that the MAMBO and SCUBA sources trace the same underlying population of high-redshift, dust-enshrouded galaxies. From a nearest-neighbour clustering analysis we find tentative evidence that the most significant MAMBO sources come in pairs, typically separated by ~23". Our MAMBO observations unambiguously confirm around half of the SCUBA sources. In a robust sub-sample of 13 SMGs detected by both MAMBO and SCUBA at a significance >=3.5sigma, only one has no radio counterpart. Furthermore, the distribution of 850/1200 micron flux density ratios for this sub-sample is consistent with the spectroscopic redshift distribution of radio-detected SMGs (Chapman et al. 2003). Finally, we have searched for evidence of a high-redshift tail of SMGs amongst the 18 MAMBO sources which are not detected by SCUBA. While we cannot rule out that some of them are SCUBA drop-outs at z >> 3, their overall 850-to-1200 micron flux distribution is statistically indistinguishable from that of the 13 SMGS which were robustly identified by both MAMBO and SCUBA.

Luminous submillimeter galaxies (SMGs) have been known for some time to produce a substantial fraction of the far-IR/submillimeter background. New observations show that these systems are not only very luminous, but also very massive. Here we present results from (a) a survey of SMGs' molecular line emission, which has confirmed Keck optical redshifts and revealed their high gas and dynamical masses; and (b) a detailed case study of the rest-frame optical emission from one SMG, which has revealed its high metallicity and old stellar population. We infer a high z~2.7 comoving number density of galaxies with baryonic masses greater than 10^11 solar masses; this is not accounted for in current semi-analytic models of galaxy evolution.

We report on the detection of the lowest CO J = 1 - 0 transition in the powerful high-redshift radio galaxy 4C60.07 at z = 3.79. The CO emission is distributed in two spatially and kinematically distinct components as was previously known from the observations of the higher excitation CO J = 4 - 3 line. The total molecular gas mass in 4C60.07 inferred from the CO J = 1 - 0 emission is M(H_2) = 1.3 x 10^11 M_sun, sufficient to fuel the inferred star-formation rate of ~1600 M_sun yr^-1 for 10^8 yrs. From our high-resolution CO J = 1 - 0 VLA maps we find the dynamical mass of 4C60.07 to be comparable to that of a giant elliptical at the present time. A significant fraction of the mass is in the form of molecular gas suggesting that 4C60.07 is in an early state of its evolution. The merging nature of 4C60.07 along with its large dynamical mass imply that this system is a giant elliptical caught in its formative stages.

Using the MPIfR Max Planck Millimeter Bolometer array (MAMBO) on the IRAM 30m Telescope we have mapped the ELAIS N2 and Lockman Hole East Fields at 1200 micron to a rms noise level of 0.8-1.0 mJy per 11" beam. The areas surveyed are 326 arcmin^2 in the ELAIS N2 field and 212 arcmin^2 in the Lockman Hole,and cover the 260 arcmin^2 previously observed by SCUBA [5].

We report IRAM Plateau de Bure, millimeter interferometry of three z=~2.4 to 3.4, SCUBA deep field galaxies. Our CO line observations confirm the rest-frame UV/optical redshifts, thus more than doubling the number of confirmed, published redshifts of the faint submillimeter population and proving their high-z nature. In all three sources our measurements of the intrinsic gas and dynamical mass are large (1e10 to 1e11 Msun). In at least two cases the data show that the submillimeter sources are part of an interacting system. Together with recent information gathered in the X-ray, optical and radio bands our observations support the interpretation that the submm-population consists of gas rich (gas to dynamical mass ratio ~0.5) and massive, composite starburst/AGN systems, which are undergoing a major burst of star formation and are evolving into m*-galaxies.

We report the first detection of the lowest CO transition in a sub-millimetre bright galaxy and extremely red object (ERO) at z = 1.44 using the Very Large Array. The total J = 1 - 0 line luminosity of ERO J164502+4626.4 is (7+-1) x 10^{10} K km s^{-1} pc^2, which yields a total molecular gas mass of ~6 x 10^{10} Msun. We also present a map of the 850-um continuum emission obtained using SCUBA, from which we infer a far-IR luminosity and dust mass of L_FIR ~ 9 x 10^{12} Lsun and M_d ~ 9 x 10^{8} Msun. We find tentative evidence that the CO and sub-mm dust emission is extended over several tens of kpc. If confirmed by high-resolution imaging, this implies that ERO J164502+4626.4 is not simply a high redshift counterpart of a typical Ultra Luminous Infrared Galaxy (ULIRG).

The SCUBA 8-mJy survey is the largest submm extragalactic mapping survey undertaken to date, centred on the Lockman Hole and ELAIS N2 regions. Here, we present new 1.4-GHz imaging of the depth and resolution necessary to reliably identify radio counterparts for 18 of 30 submm sources. Armed with this greatly improved positional information, we present and analyse new optical, near-IR and XMM-Newton X-ray imaging to identify optical/IR host galaxies to half of the submm-selected sources in those fields. As many as 15% of the submm sources detected at 1.4 GHz are resolved by the 1.4'' beam and a further 25% have more than one radio counterpart, suggesting that radio and submm emission arise from extended starbursts and that interactions are common. We note that less than a quarter of the submm-selected sample would have been recovered by targeting optically faint radio sources, underlining the selective nature of such surveys. At least 60% of the radio-confirmed optical/IR host galaxies appear to be morphologically distorted; many are composite systems - red galaxies with relatively blue companions; just over half are found to be very red or extremely red; contrary to popular belief, most are sufficiently bright to be tackled with spectrographs on 8-m telescopes. We estimate the median redshift of the S(850um) ~ 8 mJy submm galaxy population: if the radio/far-IR correlation holds at high redshift, and our sample is unbiased, we derive a conservative median z >= 2.0, or >= 2.4 using spectral templates more representative of known submm galaxies [abridged].

[ABRIDGED] We estimate the total mass distribution of the galaxy cluster CL0024+1654 from the measured source depletion due to lens magnification in the R band. Within a radius of 0.54Mpc/h, a total projected mass of (8.1+/-3.2)*10^14 M_sol/h (EdS) is measured, which corresponds to a mass- to-light ratio of M/L(B)=470+/-180. We compute the luminosity function of CL0024+1654 in order to estimate contamination of the background source counts from cluster galaxies. Three different magnification-based reconstruction methods are employed using both local and non-local techniques. We have modified the standard single power-law slope number count theory to incorporate a break and applied this to our observations. Fitting analytical magnification profiles of different cluster models to the observed number counts, we find that the cluster is best described either by a NFW model with scale radius r_s=334+/-191 kpc/h and normalisation kappa_s=0.23+/-0.08 or a power-law profile with slope xi=0.61+/-0.11, central surface mass density kappa_0=1.52+/-0.20 and assuming a core radius of r_core=35 kpc/h. The NFW model predicts that the cumulative projected mass contained within a radius R scales as M(<R)=2.9*10^14*(R/1')^[1.3-0.5lg (R/1')] M_sol/h. Finally, we have exploited the fact that flux magnification effectively enables us to probe deeper than the physical limiting magnitude of our observations in searching for a change of slope in the U band number counts. We rule out both a total flattening of the counts with a break up to U_AB<=26.6 and a change of slope, reported by some studies, from dlog N/dm=0.4->0.15 up to U_AB<=26.4 with 95% confidence.

We have obtained U and R band observations of the depletion of background galaxies due to the gravitational lensing of the galaxy cluster CL0024+1654 (z=0.39). The radial depletion curves show a significant depletion in both bands within a radius of 40''-70'' from the cluster center. This is the first time depletion is detected in the U band. This gives independent evidence for a break in the slope of the U band luminosity function at faint magnitudes. The radially averaged R band depletion curve is broader and deeper than in the U band. The differences can be attributed to the wavelength dependence of the slope of the luminosity function and to the different redshift distribution of the objects probed in the two bands. We estimate the Einstein radius of a singular isothermal sphere lens model using maximum likelihood analysis. Adopting a slope of the number counts of 0.2 and using the background density found beyond 150'' from the cluster center we find an Einstein radius of 17''+/-3'' and 25''+/-3'' in the U and R band, respectively. When combined with the redshift of the single background galaxy at z=1.675 seen as four giant arcs around 30'' from the cluster center, these values indicate a median redshift in the range 0.7 to 1.1 for the U_AB > 24 mag and R_AB > 24 mag populations.