Search

High-redshift radio(-loud) galaxies (H$z$RGs) are massive galaxies with powerful radio-loud active galactic nuclei (AGNs) and serve as beacons for protocluster identification. However, the interplay between H$z$RGs and the large-scale environment remains unclear. To understand the connection between H$z$RGs and the surrounding obscured star formation, we investigated the overdensity and spatial distribution of submillimeter-bright galaxies (SMGs) in the field of 4C\,23.56, a well-known H$z$RG at $z=2.48$. We used SCUBA-2 data ($\sigma\,{\sim}\,0.6$\,mJy) to estimate the $850\,{\rm \mu m}$ source number counts and examine the radial and azimuthal overdensities of the $850\,{\rm \mu m}$ sources in the vicinity of the H$z$RG. The angular distribution of SMGs is inhomogeneous around the H$z$RG 4C\,23.56, with fewer sources oriented along the radio jet. We also find a significant overdensity of bright SMGs (${\rm S}_{850\rm\,\mu m}\geq5\,$mJy). Faint and bright SMGs exhibit different spatial distributions. The former are concentrated in the core region, while the latter prefer the outskirts of the H$z$RG field. High-resolution observations show that the seven brightest SMGs in our sample are intrinsically bright, suggesting that the overdensity of bright SMGs is less likely due to the source multiplicity., Comment: 19 pages, 17 figures, 5 tables, accepted to A&A

Submillimetre galaxies (SMGs) are some of the most extreme star-forming systems in the Universe, whose place in the framework of galaxy evolution is as yet uncertain. It has been hypothesised that SMGs are progenitors of local early-type galaxies, requiring that SMGs generally reside in galaxy cluster progenitors at high redshift. We test this hypothesis and explore SMG environments using a narrowband VLT/HAWK-I+GRAAL study of H$\alpha$ and [OIII] emitters around an unbiased sample of three ALMA-identified and spectroscopically-confirmed SMGs at $z \sim 2.3$ and $z \sim 3.3$, where these SMGs were selected solely on spectroscopic redshift. Comparing with blank-field observations at similar epochs, we find that one of the three SMGs lies in an overdensity of emission-line sources on the $\sim4$ Mpc scale of the HAWK-I field of view, with overdensity parameter $\delta_{g} = 2.6^{+1.4}_{-1.2}$. A second SMG is significantly overdense only on $\lesssim 1.6$ Mpc scales and the final SMG is consistent with residing in a blank field environment. The total masses of the two overdensities are estimated to be $\log(M_{h}/{\rm M}_{\odot}) =$12.1--14.4, leading to present-day masses of $\log(M_{h,z=0}/{\rm M}_{\odot}) =$12.9--15.9. These results imply that SMGs occupy a range of environments, from overdense protoclusters or protogroups to the blank field, suggesting that while some SMGs are strong candidates for the progenitors of massive elliptical galaxies in clusters, this may not be their only possible evolutionary pathway., Comment: 21 pages, 12 figures. Accepted for publication in MNRAS

Measuring the environments of massive galaxies at high redshift is crucial to understanding galaxy evolution and the conditions that gave rise to the distribution of matter we see in the Universe today. While high-$z$ radio galaxies (H$z$RGs) and quasars tend to reside in protocluster-like systems, the environments of their radio-quiet counterparts are relatively unexplored, particularly in the submillimetre, which traces dust-obscured star formation. In this study we search for 850 $\mu$m-selected submillimetre galaxies in the environments of massive ($M_{\star} > 10^{11} M_{\odot}$), radio-quiet ($L_{500 {\rm MHz}} \lesssim 10^{25}$ W Hz$^{-1}$) galaxies at $z \sim 1\text{--}3$ using S2COSMOS data. By constructing number counts in circular regions of radius 1--6 arcmin and comparing with blank-field measurements, we find no significant overdensities of SMGs around massive radio-quiet galaxies at any of these scales, despite being sensitive down to overdensities of $\delta \sim 0.4$. To probe deeper than the catalogue we also examine the distribution of peaks in the SCUBA-2 SNR map, which reveals only tentative signs of any difference in the SMG densities of the radio-quiet galaxy environments compared to the blank field, and only on smaller scales (1$^{\prime}$ radii, corresponding to $\sim0.5$ Mpc) and higher SNR thresholds. We conclude that massive, radio-quiet galaxies at cosmic noon are typically in environments with $\delta\lesssim0.4$, which are either consistent with the blank field or contain only weak overdensities spanning sub-Mpc scales. The contrast between our results and studies of H$z$RGs with similar stellar masses and redshifts implies an intrinsic link between the wide-field environment and radio AGN luminosity at high redshift., Comment: 13 pages, 5 figures. Published in MNRAS

We present JWST NIRCam imaging targeting 13 $z\sim3$ infrared-luminous ($L_{\rm IR}\sim5\times10^{12}L_{\odot}$) galaxies from the ALESS survey with uniquely deep, high-resolution (0.08$''$$-$0.16$''$) ALMA 870$\mu$m imaging. The 2.0$-$4.4$\mu$m (observed frame) NIRCam imaging reveals the rest-frame near-infrared stellar emission in these submillimeter-selected galaxies (SMGs) at the same (sub-)kpc resolution as the 870$\mu$m dust continuum. The newly revealed stellar morphologies show striking similarities with the dust continuum morphologies at 870$\mu$m, with the centers and position angles agreeing for most sources, clearly illustrating that the spatial offsets reported previously between the 870$\mu$m and HST morphologies were due to strong differential dust obscuration. The F444W sizes are 78$\pm$21% larger than those measured at 870$\mu$m, in contrast to recent results from hydrodynamical simulations that predict larger 870$\mu$m sizes. We report evidence for significant dust obscuration in F444W for the highest-redshift sources, emphasizing the importance of longer-wavelength MIRI imaging. The majority of the sources show evidence that they are undergoing mergers/interactions, including tidal tails/plumes -- some of which are also detected at 870$\mu$m. We find a clear correlation between NIRCam colors and 870$\mu$m surface brightness on $\sim$1 kpc scales, indicating that the galaxies are primarily red due to dust -- not stellar age -- and we show that the dust structure on $\sim$kpc-scales is broadly similar to that in nearby galaxies. Finally, we find no strong stellar bars in the rest-frame near-infrared, suggesting the extended bar-like features seen at 870$\mu$m are highly obscured and/or gas-dominated structures that are likely early precursors to significant bulge growth., Comment: 26 pages, 10 figures. Submitted to ApJ. Comments welcome

Despite advances in observational data, theoretical models, and computational techniques to simulate key physical processes in the formation and evolution of galaxies, the stellar mass assembly of galaxies still remains an unsolved problem today. Optical spectroscopic measurements appear to show that the gas-phase metallicities of local ultra-luminous infrared galaxies (ULIRGs) are significantly lower than those of normal star-forming galaxies. This difference has resulted in the claim that ULIRGs are fueled by metal-poor gas accretion from the outskirts\cite{Mannucci10}. Here we report on a new set of gas-phase metallicity measurements making use of the far-infrared spectral lines of [O{\sc iii}]52 $\mu$m, [O{\sc iii}]88 $\mu$m, and [N{\sc iii}]57 $\mu$m instead of the usual optical lines. Photoionization models have resulted in a metallicity diagnostic based on these three lines that break the electron density degeneracy and reduce the scatter of the correlation significantly. Using new data from SOFIA and archival data from Herschel Space Observatory, we find that local ULIRGs lie on the mass-metallicity relation of star-forming galaxies and have metallicities comparable to other galaxies with similar stellar masses and star formation rates. The lack of a departure suggests that ULIRGs follow the same mass assembly mechanism as luminous star-forming galaxies and $\sim 0.3$ dex under-abundance in metallicities derived from optical lines is a result of heavily obscured metal-rich gas which has a negligible effect when using the FIR line diagnostics., Comment: Preprint submitted for publication 10/29/21

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

Previous studies measuring the clustering of submillimeter galaxies (SMGs) have based their measurements on single-dish detected sources, finding evidence for strong clustering. However, ALMA has revealed that, due to the coarse angular resolution of these instruments, single-dish sources can be comprised of multiple sources. This implies that the clustering inferred from single-dish surveys may be overestimated. Here, we measure the clustering of SMGs based on the ALESS survey, an ALMA follow-up of sources previously identified in the LABOCA ECDFS Submillimeter Survey (LESS). We present a method to measure the clustering of ALMA sources that have been previously identified using single-dish telescopes, based on forward modeling both the single-dish and the ALMA observations. We constrain upper limits for the median mass of halos hosting SMGs at $1

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

The largest Herschel extragalactic surveys, H-ATLAS and HerMES, have selected a sample of "ultrared" dusty, star-forming galaxies (DSFGs) with rising SPIRE flux densities ($S_{500} > S_{350} > S_{250}$; so-called "500 $\mu$m-risers") as an efficient way for identifying DSFGs at higher redshift ($z > 4$). In this paper, we present a large Spitzer follow-up program of 300 Herschel ultrared DSFGs. We have obtained high-resolution ALMA, NOEMA, and SMA data for 63 of them, which allow us to securely identify the Spitzer/IRAC counterparts and classify them as gravitationally lensed or unlensed. Within the 63 ultrared sources with high-resolution data, $\sim$65% appear to be unlensed, and $\sim$27% are resolved into multiple components. We focus on analyzing the unlensed sample by directly performing multi-wavelength spectral energy distribution (SED) modeling to derive their physical properties and compare with the more numerous $z \sim 2$ DSFG population. The ultrared sample has a median redshift of 3.3, stellar mass of 3.7 $\times$ 10$^{11}$ $M_{\odot}$, star formation rate (SFR) of 730 $M_{\odot}$yr$^{-1}$, total dust luminosity of 9.0 $\times$ 10$^{12}$ $L_{\odot}$, dust mass of 2.8 $\times$ 10$^9$ $M_{\odot}$, and V-band extinction of 4.0, which are all higher than those of the ALESS DSFGs. Based on the space density, SFR density, and stellar mass density estimates, we conclude that our ultrared sample cannot account for the majority of the star-forming progenitors of the massive, quiescent galaxies found in infrared surveys. Our sample contains the rarer, intrinsically most dusty, luminous and massive galaxies in the early universe that will help us understand the physical drivers of extreme star formation., Comment: 33 pages, 10 figures; accepted for publication in ApJS. The catalog tables will be available on ApJS and VizieR; authors' version is available now upon request

We present CO(J= 1-0; 3-2; 5-4; 10-9) and 1.2-kpc resolution [CII] line observations of the dusty star-forming galaxy (SFG) HXMM05 -- carried out with the Karl G. Jansky Very Large Array, the Combined Array for Research in Millimeter-wave Astronomy, the Plateau de Bure Interferometer, and the Atacama Large Millimeter/submillimeter Array, measuring an unambiguous redshift of z = 2.9850+/-0.0009. We find that HXMM05 is a hyper-luminous infrared galaxy (LIR=(4+/-1)x10^13 Lsun) with a total molecular gas mass of (2.1+/-0.7)x10^11 (alpha_CO/0.8) Msun. The CO(J=1-0) and [CII] emission are extended over ~9 kpc in diameter, and the CO line FWHM exceeds 1100 km s^-1. The [CII] emission shows a monotonic velocity gradient consistent with a disk, with a maximum rotation velocity of v_c = 616+/-100 km s^-1 and a dynamical mass of (7.7+/-3.1)x10^11 Msun. We find a star formation rate (SFR) of 2900^750_-595 Msun yr^-1. HXMM05 is thus among the most intensely star-forming galaxies known at high redshift. Photo-dissociation region modeling suggests physical conditions similar to nearby SFGs, showing extended star formation, which is consistent with our finding that the gas and dust emission are co-spatial. Its molecular gas excitation resembles the local major merger Arp 220. The broad CO and [CII] lines and a pair of compact dust nuclei suggest the presence of a late-stage major merger at the center of the extended disk, again reminiscent of Arp 220. The observed gas kinematics and conditions together with the presence of a companion and the pair of nuclei suggest that HXMM05 is experiencing multiple mergers as a part of the evolution., Comment: 32 pages, 19 figures, 7 tables. Accepted to the ApJ

In this paper, we investigate the physical associations between blended far-infrared (FIR)-emitting galaxies, in order to identify the level of line-of-sight projection contamination in the single-dish Herschel data. Building on previous work, and as part of the Herschel Extragalactic Legacy Project (HELP), we identify a sample of galaxies in the COSMOS field which are found to be both FIR-bright (typically $\sim 15$ mJy) and blended within the Herschel 250 $\mu$m beam. We identify a spectroscopic or photometric redshift for each FIR-bright source. We conduct a joint probability distribution analysis on the redshift probability density functions to determine the fraction of the FIR sources with multiple FIR-bright counterparts which are likely to be found at consistent ($\Delta z$ $< 0.01$) redshifts. We find that only 3 (0.4 per cent) of the pair permutations between counterparts are $>50$ per cent likely to be at consistent redshifts. A majority of counterparts (72 per cent) have no overlap in their redshift probability distributions whatsoever. This is in good agreement with the results of recent simulations, which indicate that single-dish observations of the FIR sky should be strongly contaminated by line of sight projection effects. We conclude that for our sample of 3.6- and 24-$\mu$m selected, FIR-bright objects in the COSMOS field, the overwhelming majority of multi-component FIR systems are line of sight projections within the 18.1 arcsec Herschel beam, rather than physical associations., Comment: Accepted for publication in MNRAS. 12 pages excluding Appendix; 16 pages total. 6 figures, 3 tables

We present the detection at 89 $\mu$m (observed frame) of the {\it Herschel}-selected gravitationally lensed starburst galaxy HATLASJ1429-0028 (also known as G15v2.19) in 15 minutes with the High-resolution Airborne Wideband Camera-plus (HAWC+) onboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The spectacular lensing system consists of an edge-on foreground disk galaxy at $z$ = 0.22 and a nearly complete Einstein ring of an intrinsic ultra-luminous infrared galaxy at $z$ = 1.03. Is this high IR luminosity powered by pure star formation (SF) or also an active galactic nucleus (AGN)? Previous nebular line diagnostics indicate that it is star-formation dominated. We perform a 27-band multi-wavelength spectral energy distribution modeling (SED) including the new SOFIA/HAWC+ data to constrain the fractional AGN contribution to the total IR luminosity. The AGN fraction in the IR turns out to be negligible. In addition, J1429-0028 serves as a testbed for comparing SED results from different models/templates and SED codes (MAGPHYS, SED3FIT, and CIGALE). We stress that star formation history is the dominant source of uncertainty in the derived stellar mass (as high as a factor of $\sim$ 10) even in the case of extensive photometric coverage. Furthermore, the detection of a source at $z$ $\sim$ 1 with SOFIA/HAWC+ demonstrates the potential of utilizing this facility for distant galaxy studies including the decomposition of SF/AGN components, which cannot be accomplished with other current facilities., Comment: 8 pages, 3 figures; accepted for publication in ApJ

The submillimeter galaxy (SMG) HERMES J105751.1+573027 (hereafter HLock01) at z = 2.9574 +/- 0.0001 is one of the brightest gravitationally lensed sources discovered in the Herschel Multi-tiered Extragalactic Survey. Apart from the high flux densities in the far-infrared, it is also extremely bright in the rest-frame ultraviolet (UV), with a total apparent magnitude m_UV = 19.7 mag. We report here deep spectroscopic observations with the Gran Telescopio Canarias of the optically bright lensed images of HLock01. Our results suggest that HLock01 is a merger system composed of the Herschel-selected SMG and an optically bright Lyman break-like galaxy (LBG), separated by only 3.3 kpc in projection. While the SMG appears very massive (M* = 5x10^11 Msun), with a highly extinguished stellar component (A_V = 4.3), the LBG is a young, lower-mass (M* = 1x10^10 Msun), but still luminous (10xL*_UV) satellite galaxy. Detailed analysis of the high signal-to-noise (S/N) rest-frame UV spectrum of the LBG shows complex kinematics of the gas, exhibiting both blueshifted and redshifted absorption components. While the blueshifted component is associated with strong galactic outflows from the massive stars in the LBG, as is common in most star-forming galaxies, the redshifted component may be associated with gas inflow seen along a favorable sightline to the LBG. We also find evidence of an extended gas reservoir around HLock01 at an impact parameter of 110 kpc, through the detection of C II 1334$\AA$ absorption in the red wing of a bright Ly-alpha emitter at z = 3.327. The data presented here highlight the power of gravitational lensing in high S/N studies to probe deeply into the physics of high-z star forming galaxies., Comment: 24 pages, 11 figures. Accepted for publication in ApJ

The SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) is a three-year JCMT Large Program aiming at reaching the 450 $\mu$m confusion limit in the COSMOS-CANDELS region, to study a representative sample of the high-redshift far-infrared galaxy population that gives rise to the bulk of the far-infrared background. We present the first-year data from STUDIES. We have reached a 450 $\mu$m noise level of 0.91~mJy for point sources at the map center, covered an area of 151 arcmin$^2$, and detected 98 and 141 sources at 4.0 and 3.5 $\sigma$, respectively. Our derived counts are best constrained in the 3.5-25 mJy regime using directly detected sources. Below the detection limits, our fluctuation analysis further constrains the slope of the counts down to 1 mJy. The resulting counts at 1-25 mJy are consistent with a power law having a slope of $-2.59$ ($\pm0.10$ for 3.5-25 mJy, and $^{+0.4}_{-0.7}$ for 1-3.5 mJy). There is no evidence of a faint-end termination or turn-over of the counts in this flux density range. Our counts are also consistent with previous SCUBA-2 blank-field and lensing cluster surveys. The integrated surface brightness from our counts down to 1 mJy is $90.0\pm17.2$ Jy deg$^{-2}$, which can account for up to $83^{+15}_{-16}\%$ of the COBE 450 $\mu$m background. We show that Herschel counts at 350 and 500 $\mu$m are significantly higher than our 450 $\mu$m counts, likely caused by its large beam and source clustering. High-angular resolution instruments like SCUBA-2 at 450 $\mu$m are therefore highly beneficial for measuring the luminosity and spatial density of high-redshift dusty galaxies., Comment: ApJ accepted. revised version uploaded

We present spectroscopic redshifts of S(870)>2mJy submillimetre galaxies (SMGs) which have been identified from the ALMA follow-up observations of 870um detected sources in the Extended Chandra Deep Field South (the ALMA-LESS survey). We derive spectroscopic redshifts for 52 SMGs, with a median of z=2.4+/-0.1. However, the distribution features a high redshift tail, with ~25% of the SMGs at z>3. Spectral diagnostics suggest that the SMGs are young starbursts, and the velocity offsets between the nebular emission and UV ISM absorption lines suggest that many are driving winds, with velocity offsets up to 2000km/s. Using the spectroscopic redshifts and the extensive UV-to-radio photometry in this field, we produce optimised spectral energy distributions (SEDs) using Magphys, and use the SEDs to infer a median stellar mass of M*=(6+/-1)x10^{10}Msol for our SMGs with spectroscopic redshifts. By combining these stellar masses with the star-formation rates (measured from the far-infrared SEDs), we show that SMGs (on average) lie a factor ~5 above the main-sequence at z~2. We provide this library of 52 template fits with robust and well-sampled SEDs available as a resource for future studies of SMGs, and also release the spectroscopic catalog of ~2000 (mostly infrared-selected) galaxies targeted as part of the spectroscopic campaign., Comment: 26 pages; 15 figures. Accepted for publication in ApJ. Redshift catalogs, spectra and SED templates available from: http://astro.dur.ac.uk/~ams/zLESS/

We present stacked average far-infrared spectra of a sample of 197 dusty, star-forming galaxies (DSFGs) at $0.005 < z < 4$ using close to 90% of the SPIRE Fourier Transform Spectrometer (FTS) extragalactic data archive from the Herschel Space Observatory based on 3.5 years of science operations. These spectra explore an observed-frame $\rm 447\,GHz-1568\,GHz$ ($\rm 191\,\mu m-671\,\mu m$) frequency (wavelength) range allowing us to observe the main atomic and molecular lines emitted by gas in the interstellar medium. The sample is sub-divided into five redshift bins at $0.005 < z < 0.05$, $0.05 < z < 0.2$, $0.2 < z < 0.5$, $0.8 < z <2$, and $2 < z < 4$. To study the dependence of observed spectral lines on total infrared luminosity, the sources in a subset of the redshift bins are stacked in luminosity bins. These stacked spectra are used to determine the average properties of the interstellar medium and dense molecular gas properties of DSFGs, in particular, the fine-structure line ([CII] 158 $\mu$m and [OI] 63 $\mu$m) luminosity ratios, and the line to far-IR luminosity ratios are used to model the gas density and radiation field strength in the photodissociation regions (PDRs). For the low-redshift sample, we additionally present the average spectral line energy distributions (SLED) of CO and $\rm{H_2O}$ rotational transitions and also consider PDR conditions based on a combination of [CI] 370 $\mu$m and 609 $\mu$m and $\rm CO (7-6)$ lines. For the high-z ($0.8 < z < 4$) sample PDR models suggest a molecular gas distribution in the presence of a radiation field that is at least a factor of 10$^3$ larger than the Milky-Way and with a neutral gas density of roughly 10$^3$ to 10$^5$ cm$^{-3}$. The corresponding PDR models for the low-z sample suggest a UV radiation field and gas density comparable to those at high-z., Comment: 30 pages, 19 figures

Submillimeter galaxies (SMGs) at $z\gtrsim1$ are luminous in the far-infrared and have star-formation rates, SFR, of hundreds to thousands of solar masses per year. However, it is unclear whether they are true analogs of local ULIRGs or whether the mode of their star formation is more similar to that in local disk galaxies. We target these questions by using Herschel-PACS to examine the conditions in the interstellar medium (ISM) in far-infrared luminous SMGs at z~1-4. We present 70-160 micron photometry and spectroscopy of the [OIV]26 micron, [FeII]26 micron, [SIII]33 micron, [SiII]34 micron, [OIII]52 micron, [NIII]57 micron, and [OI]63 micron fine-structure lines and the S(0) and S(1) hydrogen rotational lines in 13 lensed SMGs identified by their brightness in early Herschel data. Most of the 13 targets are not individually spectroscopically detected and we instead focus on stacking these spectra with observations of an additional 32 SMGs from the \herschel\ archive -- representing a complete compilation of PACS spectroscopy of SMGs. We detect [OI]63 micron, [SiII]34 micron, and [NIII]57 micron at >3sigma in the stacked spectra, determining that the average strengths of these lines relative to the far-IR continuum are $(0.36\pm0.12)\times10^{-3}$, $(0.84\pm0.17)\times10^{-3}$, and $(0.27\pm0.10)\times10^{-3}$, respectively. Using the [OIII]52/[NIII]57 emission line ratio we show that SMGs have average gas-phase metallicities $\gtrsim Z_{\rm sun}$. By using PDR modelling and combining the new spectral measurements with integrated far-infrared fluxes and existing [CII]158 micron data we show that SMGs have average gas densities, n, of $\sim10^{1-3}{\rm cm^{-3}}$ and FUV field strengths, $G_0\sim10^{2.2-4.5}$ (in Habing units: $1.6\times10^{-3}{\rm erg~cm^{-2}~s^{-1}}$), consistent with both local ULIRGs and lower luminosity star-forming galaxies., Comment: Accepted for publication in ApJ

We investigate the multiplicity of extragalactic sources detected by the Herschel Space Observatory in the COSMOS field. Using 3.6- and 24-$\mu$m catalogues, in conjunction with 250-$\mu$m data from Herschel, we seek to determine if a significant fraction of Herschel sources are composed of multiple components emitting at 250 $\mu$m. We use the XID+ code, using Bayesian inference methods to produce probability distributions of the possible contributions to the observed 250-$\mu$m flux for each potential component. The fraction of Herschel flux assigned to the brightest component is highest for sources with total 250-$\mu$m fluxes < 45 mJy; however, the flux in the brightest component is still highest in the brightest Herschel sources. The faintest 250-$\mu$m sources (30-45 mJy) have the majority of their flux assigned to a single bright component; the second brightest component is typically significantly weaker, and contains the remainder of the 250-$\mu$m source flux. At the highest 250-$\mu$m fluxes (45-110 mJy), the brightest and second brightest components are assigned roughly equal fluxes, and together are insufficient to reach 100 per cent of the 250-$\mu$m source flux. This indicates that additional components are required, beyond the brightest two components, to reproduce the observed flux. 95 per cent of the sources in our sample have a second component that contains more than 10 per cent of the total source flux. Particularly for the brightest Herschel sources, assigning the total flux to a single source may overestimate the flux contributed by around 150 per cent., Comment: 12 pages, 7 figures. MNRAS accepted

It is thought that a class of broad absorption line (BAL) QSOs, characterised by Fe absorption features in their UV spectra (called `FeLoBALs'), could mark a transition stage between the end of an obscured starburst event and a youthful QSO beginning to shed its dust cocoon, where Fe has been injected into the interstellar medium by the starburst. To test this hypothesis we have undertaken deep SCUBA-2 850 $\mu$m observations of a sample of 17 FeLoBAL QSOs with 0.89 $\leq$ z $\leq$ 2.78 and -23.31 $\leq$ M$_{B}$ $\leq$-28.50 to directly detect an excess in the thermal emission of the dust which would probe enhanced star-formation activity. We find that FeLoBALs are not luminous sources in the submillimetre, none of them are individually detected at 850 $\mu$m, nor as a population through stacking ($F_{s}=1.14\pm0.58$ mJy). Statistical and survival analyses reveal that FeLoBALs have sub-mm properties consistent with BAL and non-BAL QSOs with matched redshifts and magnitudes. An SED fitting analysis shows that the FIR emission is dominated by AGN activity, and a starburst component is required only in 6/17 sources of our sample; moreover the integrated total luminosity of 16/17 sources is L$\geq$10$^{12}$L$_{\odot}$, high enough to classify FeLoBALs as infrared luminous. In conclusion, we do not find any evidence in support of FeLoBAL QSOs being a transition population between a ULIRG and an unobscured QSO; in particular, FeLoBALs are not characterized by a cold starburst which would support this hypothesis., Comment: 15 pages, 7 figures, MNRAS in press

We show that the use of red colour as the basis for selecting candidate high redshift dusty galaxies from surveys made with Herschel has proved highly successful. The highest redshift such object, HFLS3, lies at z=6.34 and numerous other sources have been found. Spectroscopic followup confirms that most of these lie at z>4. These sources are found in such numbers that they represent a challenge to current models of galaxy evolution. We also examine the prospects for finding dusty galaxies at still higher redshifts. These would not appear in the SPIRE surveys from Herschel but would be detected in longer wavelength, submm, surveys. Several such `SPIRE-dropouts' have been found and are now subject to followup observations., Comment: 4 pages, 2 figures. Accepted for publication in Proceedings of IAU Symposium 319 - Galaxies at High Redshift and Their Evolution over Cosmic Time

Measurements of the UV-continuum slopes provide valuable information on the physical properties of galaxies forming in the early universe, probing the dust reddening, age, metal content, and even the escape fraction. While constraints on these slopes generally become more challenging at higher redshifts as the UV continuum shifts out of the Hubble Space Telescope bands (particularly at z>7), such a characterisation actually becomes abruptly easier for galaxies in the redshift window z=9.5-10.5 due to the Spitzer/IRAC 3.6um-band probing the rest-UV continuum and the long wavelength baseline between this Spitzer band and the Hubble H-band. Higher S/N constraints on the UV slope are possible at z~10 than at z=8. Here we take advantage of this opportunity and five recently discovered bright z=9.5-10.5 galaxies to present the first measurements of the mean slope for a multi-object sample of galaxy candidates at z~10. We find the measured observed slopes of these candidates are $-2.1\pm0.3\pm0.2$ (random and systematic), only slightly bluer than the measured slopes at 3.5

We present Hubble Space Telescope (HST) WFC3 imaging and grism spectroscopy observations of the {\it Herschel}-selected gravitationally-lensed starburst galaxy HATLASJ1429-0028. The lensing system consists of an edge-on foreground disk galaxy at $z=0.218$ with a nearly complete Einstein ring of the infrared luminous galaxy at $z=1.027$. The WFC3 spectroscopy with G102 and G141 grisms, covering the wavelength range of 0.8 to 1.7 $\mu$m, resulted in detections of H$\alpha$+[NII], H$\beta$, [SII], and [OIII] for the background galaxy from which we measure line fluxes and ratios. The Balmer line ratio H$\alpha$/H$\beta$ of $7.5 \pm 4.4$, when corrected for [NII], results in an extinction for the starburst galaxy of $E(B-V)=0.8 \pm 0.5$. The $H\alpha$ based star-formation rate, when corrected for extinction, is $60 \pm 50$ M$_{\odot}$ yr$^{-1}$, lower than the instantaneous star-formation rate of 390 $\pm$ 90 M$_{\odot}$ yr$^{-1}$ from the total IR luminosity. We also compare the nebular line ratios of HATLASJ1429-0028 with other star-forming and sub-mm bright galaxies. The nebular line ratios are consistent with an intrinsic ultra-luminous infrared galaxy with no evidence for excitation by an active galactic nucleus (AGN). We estimate the metallicity, $12 + log(O/H)$, of HATLASJ1429-0028 to be 8.49 $\pm$ 0.16. Such a low value is below the average relations for stellar mass vs. metallicity of galaxies at $z \sim 1$ for a galaxy with stellar mass of $\sim 2 \times 10^{11}$ M$_{\odot}$. The combination of high stellar mass, lack of AGN indicators, low metallicity, and the high star-formation rate of HATLASJ1429-0028 suggest that this galaxy is currently undergoing a rapid formation., Comment: 7 pages, 8 figures, Accepted for publication in the Astrophysical Journal

We have studied in detail a sample of 967 SPIRE sources with 5-sigma detections at 350 and 500 micron and associations with Spitzer-SWIRE 24 micron galaxies in the HerMES-Lockman survey area, fitting their mid- and far-infrared, and submillimetre, SEDs in an automatic search with a set of six infrared templates. For almost 300 galaxies we have modelled their SEDs individually to ensure the physicality of the fits. We confirm the need for the new cool and cold cirrus templates, and also of the young starburst template, introduced in earlier work. We also identify 109 lensing candidates via their anomalous SEDs and provide a set of colour-redshift constraints which allow lensing candidates to be identified from combined Herschel and Spitzer data. The picture that emerges of the submillimetre galaxy population is complex, comprising ultraluminous and hyperluminous starbursts, lower luminosity galaxies dominated by interstellar dust emission, lensed galaxies and galaxies with surprisingly cold (10-13K) dust. 11 % of 500$\mu$m selected sources are lensing candidates. 70 % of the unlensed sources are ultraluminous infrared galaxies and 26 % are hyperluminous. 34 % are dominated by optically thin interstellar dust ('cirrus') emission, but most of these are due to cooler dust than is characteristic of our Galaxy. At the highest infrared luminosities we see SEDs dominated by M82, Arp220 and young starburst types, in roughly equal proportions., Comment: 18 pages, 46 figures, accepted for publication MNRAS. arXiv admin note: text overlap with arXiv:1006.0120

We discuss the restframe UV emission from the starbursting galaxy HFLS3 at z=6.34, discovered in Herschel/SPIRE data due to its red color in the submm wavelengths from 250-500 um. The apparent inst. SFR of HFLS3 inferred from the total FIR luminosity measured with over 15 photometric data points between 100 to 1000 um is 2900 Msun/yr. Keck/NIRC2 Ks band adaptive optics imaging data showed two potential NIR counterparts near HFLS3. Previously, the northern galaxy was taken to be in the foreground at z=2.1 while the southern galaxy was assumed to HFLS3's NIR counterpart. New HST/WFC3 and ACS imaging data show both optically bright galaxies are in the foreground at z<6. A new lensing model based on HST data and mm-wave continuum emission yields a magnification of 2.2+/-0.3. The lack of multiple imaging constrains the lensing magnification to be lower than either 2.7 or 3.5 at the 95% confidence level for the two scenarios, which attribute one or two components to HFLS3 in the source plane. Correcting for gravitational lensing, the inst. SFR is 1320 Msun/yr with the 95% confidence lower limit around 830 Msun/yr. Using models for the restframe UV to FIR SED, the ave. SFR over the last 100 Myr is around 660 Msun/yr. The dust and stellar masses of HFLS3 from the same SED models are 3x10^8 Msun and ~5x10^10 Msun, respectively, with large systematic uncertainties on assumptions related to the SED model. With HST/WFC3 images we also find diffuse NIR emission about 0.5" (~3 kpc) SW of HFLS3 that remains undetected in the ACS data. The emission has a photometric redshift consistent with either z~6 or a dusty galaxy template at z~2. If at the same redshift as HFLS3 the detected diffuse emission could be part of the complex merger system that could be triggering the starburst. Alternatively, it could be part of the foreground structure at z~2.1 that is responsible for lensing of HFLS3., Comment: Submitted to ApJ, 10 pages, 6 figures

We study the far-infrared properties of 498 Lyman Alpha Emitters (LAEs) at z=2.8, 3.1 and 4.5 in the Extended Chandra Deep Field-South, using 250, 350 and 500 micron data from the Herschel Multi-tiered Extragalactic Survey (HerMES) and 870 micron data from the LABOCA ECDFS Submillimeter Survey (LESS). None of the 126, 280 or 92 LAEs at z=2.8, 3.1 and 4.5, respectively, are individually detected in the far-infrared data. We use stacking to probe the average emission to deeper flux limits, reaching $1\sigma$ depths of ~0.1 to 0.4 mJy. The LAEs are also undetected at $\ge3\sigma$ in the stacks, although a $2.5\sigma$ signal is observed at 870 micron for the z=2.8 sources. We consider a wide range of far-infrared spectral energy distributions (SEDs), including a M82 and an Sd galaxy template, to determine upper limits on the far-infrared luminosities and far-infrared-derived star-formation rates of the LAEs. These star-formation rates are then combined with those inferred from the Ly$\alpha$ and UV emission to determine lower limits on the LAEs Ly$\alpha$ escape fraction ($f_{\rm esc}($Ly$\alpha$)). For the Sd SED template, the inferred LAEs $f_{\rm esc}($Ly$\alpha$) are $\gtrsim30%$ ($1\sigma$) at z=2.8, 3.1 and 4.5, which are all significantly higher than the global $f_{\rm esc}($Ly$\alpha$) at these redshifts. Thus, if the LAEs $f_{\rm esc}($Ly$\alpha$) follows the global evolution then they have warmer far-infrared SEDs than the Sd galaxy template. The average and M82 SEDs produce lower limits on the LAE $f_{\rm esc}($Ly$\alpha$) of ~10 to 20% ($1\sigma$), all of which are slightly higher than the global evolution of $f_{\rm esc}($Ly$\alpha$) but consistent with it at the 2 to 3$\sigma$ level., Comment: 9 pages, 4 figures. Accepted for publication in ApJ

MS$\,$0451.6$-$0305 is a rich galaxy cluster whose strong lensing is particularly prominent at submm wavelengths. We combine new SCUBA-2 data with imaging from Herschel SPIRE and PACS and HST in order to try to understand the nature of the sources being lensed. In the region of the "giant submm arc," we uncover seven multiply imaged galaxies (up from the previously known three), of which six are found to be at a redshift of $z\sim2.9$, and possibly constitute an interacting system. Using a novel forward-modelling approach, we are able to simultaneously deblend and fit SEDs to the individual galaxies that contribute to the giant submm arc, constraining their dust temperatures, far infrared luminosities and star formation rates. The submm arc first identified by SCUBA can now be seen to be composed of at least five distinct sources, four of these within the galaxy group at $z\sim2.9$. The total unlensed luminosity for this galaxy group is $(3.1\pm0.3) \times 10^{12}\,\mathrm{L}_\odot$, which gives an unlensed star formation rate of $(450\pm50)$ M$_\odot$ yr$^{-1}$. From the properties of this system, we see no evidence of evolution towards lower temperatures in the dust temperature versus far-infrared luminosity relation for high redshift galaxies.

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., Comment: Accepted by ApJ. 45 pages, 16 figures

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., Comment: Accepted for publication in MNRAS. 24 pages, 12 figures

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., Comment: 9 pages, 3 figures, 2 tables; revised to match version accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS. 10 pages, 5 figures, 2 tables

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

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane., Comment: ApJ in press. 31 pages, 16 figures, 5 tables. This version updated to match accepted version

We use spitzer-IRAC data to identify near-infrared counterparts to submillimeter galaxies detected with Herschel-SPIRE at 250um in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Using a likelihood ratio analysis we identify 146 reliable IRAC counterparts to 123 SPIRE sources out of the 159. We find that, compared to the field population, the SPIRE counterparts occupy a distinct region of 3.6 and 4.5um color-magnitude space, and we use this property to identify a further 23 counterparts to 13 SPIRE sources. The IRAC identification rate of 86% is significantly higher than those that have been demonstrated with wide-field ground-based optical and near-IR imaging of Herschel fields. We estimate a false identification rate of 3.6%, corresponding to 4 to 5 sources. Among the 73 counterparts that are undetected in SDSS, 57 have both 3.6 and 4.5um coverage. Of these 43 have [3.6] - [4.5]> 0 indicating that they are likely to be at z > 1.4. Thus, ~ 40% of identified SPIRE galaxies are likely to be high redshift (z > 1.4) sources. We discuss the statistical properties of the IRAC-identified SPIRE galaxy sample including far-IR luminosities, dust temperatures, star-formation rates, and stellar masses. The majority of our detected galaxies have 10^10 to 10^11 L_sun total IR luminosities and are not intense starbursting galaxies as those found at z ~ 2, but they have a factor of 2 to 3 above average specific star-formation rates compared to near-IR selected galaxy samples., Comment: 16 pages, 12 figures, two tables. Accepted for publication in ApJ

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs., Comment: White paper supplement to the proposal submitted by the HSLS science team to ESA for Herschel open-time programs

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., Comment: Accepted for publication in MNRAS Letters July 5, 2010: 5 pages, 3 figures

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., Comment: 10 pages, 4 figures, 1 table, Accepted for publication in MNRAS Feb 23