Search

Wide, deep, blind continuum surveys at submillimetre/millimetre (submm/mm) wavelengths are required to provide a full inventory of the dusty, distant Universe. However, conducting such surveys to the necessary depth, with sub-arcsec angular resolution, is prohibitively time-consuming, even for the most advanced submm/mm telescopes. Here, we report the most recent results from the ALMACAL project, which exploits the 'free' calibration data from the Atacama Large Millimetre/submillimetre Array (ALMA) to map the lines of sight towards and beyond the ALMA calibrators. ALMACAL has now covered 1,001 calibrators, with a total sky coverage around 0.3 deg2, distributed across the sky accessible from the Atacama desert, and has accumulated more than 1,000h of integration. The depth reached by combining multiple visits to each field makes ALMACAL capable of searching for faint, dusty, star-forming galaxies (DSFGs), with detections at multiple frequencies to constrain the emission mechanism. Based on the most up-to-date ALMACAL database, we report the detection of 186 DSFGs with flux densities down to S870um ~ 0.2mJy, comparable with existing ALMA large surveys but less susceptible to cosmic variance. We report the number counts at five wavelengths between 870um and 3mm, in ALMA bands 3, 4, 5, 6 and 7, providing a benchmark for models of galaxy formation and evolution. By integrating the observed number counts and the best-fitting functions, we also present the resolved fraction of the cosmic infrared background (CIB) and the CIB spectral shape. Combining existing surveys, ALMA has currently resolved about half of the CIB in the submm/mm regime., Comment: 22 pages, 16 figures, accepted for publication in MNRAS

We report the detection of the ground state rotational emission of ammonia, ortho-NH$_3$ $(J_K=1_0\rightarrow0_0)$ in a gravitationally lensed, intrinsically hyperluminous, star-bursting galaxy at $z=2.6$. The integrated line profile is consistent with other molecular and atomic emission lines which have resolved kinematics well-modelled by a 5 kpc-diametre rotating disc. This implies that the gas responsible for NH$_3$ emission is broadly tracing the global molecular reservoir, but likely distributed in pockets of high density ($n\gtrsim5\times10^4$ cm$^{-3}$). With a luminosity of $2.8\times10^{6}$ $L_\odot$, the NH$_3$ emission represents $2.5\times10^{-7}$ of the total infrared luminosity of the galaxy, comparable to the ratio observed in the Kleinmann-Low nebula in Orion and consistent with sites of massive star formation in the Milky Way. If $L_{\rm NH_3}/L_{\rm IR}$ serves as a proxy for the 'mode' of star formation, this hints that the nature of star formation in extreme starbursts in the early Universe is similar to that of Galactic star-forming regions, with a large fraction of the cold interstellar medium in this state, plausibly driven by a storm of violent disc instabilities in the gas-dominated disc. This supports the 'full of Orions' picture of star formation in the most extreme galaxies seen close to the peak epoch of stellar mass assembly., Comment: 5 pages, 3 figures, accepted for publication in MNRAS Letters

We analyse the molecular and atomic emission lines of 71 bright Herschel-selected galaxies between redshifts 1.4 to 4.6 detected by the Atacama Large Millimetre/submillimetre Array. These lines include a total of 156 CO, [C I], and H2O emission lines. For 46 galaxies, we detect two transitions of CO lines, and for these galaxies we find gas properties similar to those of other dusty star-forming galaxy (DSFG) samples. A comparison to photo-dissociation models suggests that most of Herschel-selected galaxies have similar interstellar medium conditions as local infrared-luminous galaxies and high-redshift DSFGs, although with denser gas and more intense far-ultraviolet radiation fields than normal star-forming galaxies. The line luminosities agree with the luminosity scaling relations across five orders of magnitude, although the star-formation and gas surface density distributions (i.e., Schmidt-Kennicutt relation) suggest a different star-formation phase in our galaxies (and other DSFGs) compared to local and low-redshift gas-rich, normal star-forming systems. The gas-to-dust ratios of these galaxies are similar to Milky Way values, with no apparent redshift evolution. Four of 46 sources appear to have CO line ratios in excess of the expected maximum (thermalized) profile, suggesting a rare phase in the evolution of DSFGs. Finally, we create a deep stacked spectrum over a wide rest-frame frequency (220-890 GHz) that reveals faint transitions from HCN and CH, in line with previous stacking experiments., Comment: 30 pages, 17 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal. Comments are warmly welcomed

We have carried out the first spatially-resolved investigation of the multi-phase interstellar medium (ISM) at high redshift, using the z=4.24 strongly-lensed sub-millimetre galaxy H-ATLASJ142413.9+022303 (ID141). We present high-resolution (down to ~350 pc) ALMA observations in dust continuum emission and in the CO(7-6), H_2O (2_{1,1} - 2_{0,2}), CI(1-0) and CI(2-1) lines, the latter two allowing us to spatially resolve the cool phase of the ISM for the first time. Our modelling of the kinematics reveals that the system appears to be dominated by a rotationally-supported gas disk with evidence of a nearby perturber. We find that the CI(1-0) line has a very different distribution to the other lines, showing the existence of a reservoir of cool gas that might have been missed in studies of other galaxies. We have estimated the mass of the ISM using four different tracers, always obtaining an estimate in the range (3.2-3.8) x 10^{11} M_sol, significantly higher than our dynamical mass estimate of (0.8-1.3) x 10^{11} M_sol. We suggest that this conflict and other similar conflicts reported in the literature is because the gas-to-tracer ratios are ~4 times lower than the Galactic values used to calibrate the ISM in high-redshift galaxies. We demonstrate that this could result from a top-heavy initial mass function and strong chemical evolution. Using a variety of quantitative indicators, we show that, extreme though it is at z=4.24, ID141 will likely join the population of quiescent galaxies that appears in the Universe at z~3., This article has been accepted for publication in MNRAS, published by Oxford University Press on behalf of the Royal Astronomical Society. 24 pages

In this paper we present the third data release (DR3) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We identify likely near-infrared counterparts to submillimetre sources in the South Galactic Pole (SGP) field using the VISTA VIKING survey. We search for the most probable counterparts within 15 arcsec of each Herschel source using a probability measure based on the ratio between the likelihood the true counterpart is found close to the submillimetre source and the likelihood that an unrelated object is found in the same location. For 110 374 (57.0$\%$) sources we find galaxies on the near-infrared images where the probability that the galaxy is associated to the source is greater than 0.8. We estimate the false identification rate to be 4.8$\%$, with a probability that the source has an associated counterpart on the VIKING images of 0.835$\pm$0.009. We investigate the effects of gravitational lensing and present 41 (0.14 deg$^{-2}$) candidate lensed systems with observed flux densities > 100 mJy at 500 $��$m. We include in the data release a probability that each source is gravitationally lensed and discover an additional 5 923 sources below 100 mJy that have a probability greater than 0.94 of being gravitationally lensed. We estimate that $\sim$ 400 - 1 000 sources have multiple true identifications in VIKING based on the similarity of redshift estimates for multiple counterparts close to a Herschel source. The data described in this paper can be found at the H-ATLAS website., 16 pages, 13 figures, accepted for publication in MNRAS

We examine the far-IR properties of a sample of 5391 optically selected QSOs in the 0.544.7, using SPIRE data from Herschel-ATLAS. We split the sample in a grid of 74 luminosity-redshift bins and compute the average optical-infrared spectral energy distribution (SED) in each bin. By normalising an intrinsic AGN template to the AGN optical power (at 5100A) we decompose the total infrared emission (L_IR; 8-1000um) into an AGN (L_IR,AGN) and star-forming component (L_IR,SF). We find that the AGN contribution to L_IR increases as a function of AGN power which manifests as a reduction of the `far-IR bump' in the average QSO SEDs. We note that L_IR,SF does not correlate with AGN power; the mean star formation rates (SFRs) of AGN host galaxies are a function of redshift only and they range from ~6 Msun/yr at z~0 to a plateau of, Comment: 20 pages, 10 figures, accepted for publication in MNRAS

We present a self-consistent cross-calibration of the three main molecular gas mass tracers in galaxies, the $\rm ^{12}CO$(1-0), [CI]($^3P_1$-$^3P_0$) lines, and the submm dust continuum emission, using a sample of 407 galaxies, ranging from local disks to submillimetre-selected galaxies (SMGs) up to $z \approx 6$. A Bayesian method is used to produce galaxy-scale universal calibrations of these molecular gas indicators, that hold over 3-4 orders of magnitude in infrared luminosity, $L_{\rm IR}$. Regarding the dust continuum, we use a mass-weighted dust temperature, $T_{\rm mw}$, determined using new empirical relations between temperature and luminosity. We find the average $L/M_{\rm mol}$ gas mass conversion factors to be $\alpha_{850}= 6.9\times10^{12}\,\rm W\,Hz^{-1}\,M_{\odot}^{-1}$, $\alpha_{\rm CO} = \rm 4\,M_{\odot} (K\,km\,s^{-1}\,pc^2)^{-1}$ and $\alpha_{\rm CI} = \rm 17.0 \,M_{\odot} (K\,km\,s^{-1}\,pc^2)^{-1}$, based on the assumption that the mean dust properties of the sample ($\kappa_H$ = gas-to-dust ratio/dust emissivity) will be similar to those of local metal rich galaxies and the MW. The tracer with the least intrinsic scatter is [CI](1-0), while CO(1-0) has the highest. The conversion factors show a weak but significant correlation with $L_{\rm IR}$. Assuming dust properties typical of metal-rich galaxies, we infer a neutral carbon abundance $X_{\rm CI} = [C^0/\rm mol]=1.6\times 10^{-5}$, similar to that in the MW. We find no evidence for bimodality of $\alpha_{\rm CO}$ between main-sequence (MS) galaxies and those with extreme star-formation intensity, i.e. ULIRGs and SMGs. The means of the three conversion factors are found to be similar between MS galaxies and ULIRGs/SMGs, to within 10-20%. We show that for metal-rich galaxies, near-universal average values for $\alpha_{\rm CO}$, $X_{\rm CI}$ and $\kappa_H$ are adequate for global molecular gas estimates., Comment: 42 pages, accepted for publication in MNRAS. Abstract abridged. Data tables will be available on CDS imminently. Minor changes from proof stage including swapping order of Tables 8 and 9 and small corrections to the bibliography

While unobscured and radio-quiet active galactic nuclei are regularly being found at redshifts $z > 6$, their obscured and radio-loud counterparts remain elusive. We build upon our successful pilot study, presenting a new sample of low-frequency-selected candidate high-redshift radio galaxies (HzRGs) over a sky area twenty times larger. We have refined our selection technique, in which we select sources with curved radio spectra between 72-231 MHz from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. In combination with the requirements that our GLEAM-selected HzRG candidates have compact radio morphologies and be undetected in near-infrared $K_{\rm s}$-band imaging from the Visible and Infrared Survey Telescope for Astronomy Kilo-degree Infrared Galaxy (VIKING) survey, we find 51 new candidate HzRGs over a sky area of approximately 1200 deg$^2$. Our sample also includes two sources from the pilot study: the second-most distant radio galaxy currently known, at $z=5.55$, with another source potentially at $z \sim 8$. We present our refined selection technique and analyse the properties of the sample. We model the broadband radio spectra between 74 MHz and 9 GHz by supplementing the GLEAM data with both publicly available data and new observations from the Australia Telescope Compact Array at 5.5 and 9 GHz. In addition, deep $K_{\rm s}$-band imaging from the High-Acuity Widefield $K$-band Imager (HAWK-I) on the Very Large Telescope and from the Southern Herschel Astrophysical Terahertz Large Area Survey Regions $K_{\rm s}$-band Survey (SHARKS) is presented for five sources. We discuss the prospects of finding very distant radio galaxies in our sample, potentially within the epoch of reionisation at $z \gtrsim 6.5$., 49 pages, 3 figures (one of which is a multi-page figure with 102 separate panels), 9 tables, accepted for publication in PASA

This is a correction to: Monthly Notices of the Royal Astronomical Society, Volume 473, Issue 2, January 2018, Pages 1751–1773, https://doi.org/10.1093/mnras/stx2267

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., Comment: 16 pages, 11 figures, includes revisions recommended by ApJ referee

We present a multi-wavelength analysis of 52 sub-millimeter galaxies (SMGs), identified using ALMA 870$\mu$m continuum imaging in a pilot program to precisely locate bright SCUBA2-selected sub-mm sources in the UKIDSS Ultra Deep Survey (UDS) field. Using the available deep (especially near-infrared), panoramic imaging of the UDS field at optical-to-radio wavelengths we characterize key properties of the SMG population. The median photometric redshift of the bright ALMA/SCUBA-2 UDS (AS2UDS) SMGs that are detected in a sufficient number of wavebands to derive a robust photometric redshift is $z$=2.65$\pm$0.13. However, similar to previous studies, 27% of the SMGs are too faint at optical-to-near-infrared wavelengths to derive a reliable photometric redshift. Assuming that these SMGs lie at z$\gtrsim$3 raises the median redshift of the full sample to $z$=2.9$\pm$0.2. A subset of 23, unlensed, bright AS2UDS SMGs have sizes measured from resolved imaging of their rest-frame far-infrared emission. We show that the extent and luminosity of the far-infrared emission are consistent with the dust emission arising from regions that are optically thick, on average, at a wavelength of $\lambda_0$$\ge$75$\mu$m (1-$\sigma$ dispersion of 55-90$\mu$m). Using the dust masses derived from our optically-thick spectral energy distribution models we determine that these galaxies have a median hydrogen column density of $N_{H}$=9.8$_{-0.7}^{+1.4}$$\times$10$^{23}$cm$^{-2}$, or a corresponding median $V$-band obscuration of $A_\mathrm{v}$=540$^{+80}_{-40}$mag, averaged along the line of sight to the source of their restframe $\sim$200$\mu$m emission. We discuss the implications of this extreme attenuation by dust for the multiwavelength study of dusty starbursts and reddening-sensitive tracers of star formation., Comment: Accepted for publication in ApJ

Exploiting the sensitivity and spatial resolution of the Atacama Large Millimeter/submillimeter Array (ALMA), we have studied the morphology and the physical scale of the interstellar medium - both gas and dust - in SGP38326, an unlensed pair of interacting starbursts at $z= 4.425$. SGP38326 is the most luminous star bursting system known at $z > 4$ with an IR-derived ${\rm SFR \sim 4300 \,} M_\odot \, {\rm yr}^{-1}$. SGP38326 also contains a molecular gas reservoir among the most massive ever found in the early Universe, and it is the likely progenitor of a massive, red-and-dead elliptical galaxy at $z \sim 3$. Probing scales of $\sim 0.1"$ or $\sim 800 \, {\rm pc}$ we find that the smooth distribution of the continuum emission from cool dust grains contrasts with the more irregular morphology of the gas, as traced by the [CII] fine structure emission. The gas is also extended over larger physical scales than the dust. The velocity information provided by the resolved [CII] emission reveals that the dynamics of the two components of SGP38326 are compatible with disk-like, ordered rotation, but also reveals an ISM which is turbulent and unstable. Our observations support a scenario where at least a subset of the most distant extreme starbursts are highly dissipative mergers of gas-rich galaxies., Submitted to ApJ

We have carried out a pilot study for the SCUBA-2 'All-Sky' Survey, SASSy, a wide and shallow mapping project at 850 microns, designed to find rare objects, both Galactic and extragalactic. Two distinct sets of exploratory observations were undertaken, and used to test the SASSy approach and data reduction pipeline. The first was a 0.5 by 0.5 degrees map around the nearby galaxy NGC 2559. The galaxy was easily detected at 156 mJy, but no other convincing sources are present in the map. Comparison with other galaxies with similar wavelength coverage indicates that NGC 2559 has relatively warm dust. The second observations cover 1 square degree around the W5-E HII region. As well as diffuse structure in the map, a filtering approach was able to extract 27 compact sources with signal-to-noise greater than 6. By matching with data at other wavelengths we can see that the SCUBA-2 data can be used to discriminate the colder cores. Together these observations show that the SASSy project will be able to meet its original goals of detecting new bright sources which will be ideal for follow-up observations with other facilities.

We have reduced the data taken with the Spectral and Photometric Imaging Receiver (SPIRE) photometer on board the Herschel Space Observatory in the Science Demonstration Phase (SDP) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We describe the data reduction, which poses specific challenges, both because of the sheer size of the data, and because only two scans are made for each region. We implement effective solutions to process the bolometric timelines into maps, and show that correlations among detectors are negligible, and that the photometer is stable on time scales up to 250 s. This is longer than the time the telescope takes to cross the observed sky region, and it allows us to use naive binning methods for an optimal reconstruction of the sky emission. The maps have equal contribution of confusion and white instrumental noise, and the latter is estimated to 5.3, 6.4, and 6.7 mJy/beam (1-{\sigma}), at 250, 350, and 500 \mu{m}, respectively. This pipeline is used to reduce other H-ATLAS observations, as they became available, and we discuss how it can be used with the optimal map maker implemented in the Herschel Interactive Processing Environment (HIPE), to improve computational efficiency and stability. The SDP dataset is available from this http URL

We use data from the Herschel-ATLAS to investigate the evolution of the far-infrared--radio correlation over the redshift range 0 5sigma sources in the Herschel-ATLAS Science Demonstration Field and cross-matching these data with radio data from the Faint Images of the Radio Sky at Twenty-Centimetres (FIRST) survey and the NRAO VLA Northern Sky Survey (NVSS), we obtain 104 radio counterparts to the Herschel sources. With these data we find no evidence for evolution in the far-infrared--radio correlation over the redshift range 0

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux = 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

We present spectral energy distributions (SEDs) for 68 Herschel sources detected at 5-sigma at 250, 350 and 500 mu in the HerMES SWIRE-Lockman field. We explore whether existing models for starbursts, quiescent star-forming galaxies and for AGN dust tori are able to model the full range of SEDs measured with Herschel. We find that while many galaxies (~ 56 %) are well fitted with the templates used to fit IRAS, ISO and Spitzer sources, for about half the galaxies two new templates are required: quiescent ('cirrus') models with colder (10-20 K) dust, and a young starburst model with higher optical depth than Arp 220. Predictions of submillimetre fluxes based on model fits to 4.5-24 mu data agree rather poorly with the observed fluxes, but the agreement is better for fits to 4.5-70 mu data. Herschel galaxies detected at 500 mu tend to be those with the very highest dust masses.

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, =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., Comment: 8 pages, 3 figures, accepted in ApJ

We present the scientific motivation and observing plan for an upcoming detection survey for debris disks using the James Clerk Maxwell Telescope. The SCUBA-2 Unbiased Nearby Stars (SUNS) Survey will observe 500 nearby main sequence and sub-giant stars (100 of each of the A, F, G, K and M spectral classes) to the 850 micron extragalactic confusion limit to search for evidence of submillimeter excess, an indication of circumstellar material. The survey distance boundaries are 8.6, 16.5, 22, 25 and 45 pc for M, K, G, F and A stars, respectively, and all targets lie between the declinations of -40 deg to 80 deg. In this survey, no star will be rejected based on its inherent properties: binarity, presence of planetary companions, spectral type or age. This will be the first unbiased survey for debris disks since IRAS. We expect to detect ~125 debris disks, including ~50 cold disks not detectable in current shorter wavelength surveys. A substantial amount of complementary data will be required to constrain the temperatures and masses of discovered disks. High resolution studies will likely be required to resolve many of the disks. Therefore, these systems will be the focus of future observational studies using a variety of observatories to characterize their physical properties. For non-detected systems, this survey will set constraints (upper limits) on the amount of circumstellar dust, of typically 200 times the Kuiper Belt mass, but as low as 10 times the Kuiper Belt mass for the nearest stars in the sample (approximately 2 pc)., Comment: 11 pages, 7 figures (3 color), accepted by the Publications of the Astronomical Society of the Pacific

We present Keck-Adaptive Optics and Hubble Space Telescope high resolution near-infrared (IR) imaging for 500 um-bright candidate lensing systems identified by the Herschel Multi-tiered Extra-galactic Survey (HerMES) and Herschel Astrophysical Terahertz Survey (H-ATLAS). Out of 87 candidates with near-IR imaging, 15 (~17%) display clear near-IR lensing morphologies. We present near-IR lens models to reconstruct and recover basic rest-frame optical morphological properties of the background galaxies from 12 new systems. Sources with the largest near-IR magnification factors also tend to be the most compact, consistent with the size bias predicted from simulations and pre- vious lensing models for sub-millimeter galaxies. For four new sources that also have high-resolution sub-mm maps, we test for differential lensing between the stellar and dust components and find that the 880 um magnification factor (u_880) is ~1.5 times higher than the near-IR magnification factor (u_NIR), on average. We also find that the stellar emission is ~2 times more extended in size than dust. The rest-frame optical properties of our sample of Herschel-selected lensed SMGs are consistent with those of unlensed SMGs, which suggests that the two populations are similar., Comment: 27 pages, 13 figures, 7 tables, ApJ accepted

Despite extensive observational efforts, the brightest sub-mm source in the Hubble Deep Field, HDF850.1, has failed to yield a convincing optical/infrared identifica- tion almost 4 years after its discovery. This failure is all the more notable given the availability of supporting multi-frequency data of unparalleled depth, and sub-arcsec positional accuracy for the sub-mm/mm source. Consequently, HDF850.1 has become a test case of the possibility that the most violently star-forming objects in the universe are too red and/or distant to be seen in the deepest optical images. Here we report the discovery of the host galaxy of HDF850.1. This object has been revealed by careful analysis of a new, deep Kimage of the HDF obtained with the Subaru 8.2-m telescope. Its reality is confirmed by a similar analysis of the HST NICMOS F160W image of the same region. This object is extremely faint (K ≃ 23.5), clumpy (on sub-arcsec scales) and very red (I − K > 5.2; H − K = 1.4 ± 0.35). The likelihood that it is the correct identification is strongly reinforced by a reanalysis of the combined MERLIN+VLA 1.4-GHz map of the field which provides a new radio detection of HDF850.1 only 0.1arcsec from the new near-infrared counterpart, and with sufficient positional accuracy to exclude all previously considered alternative optical candidates. We have calculated new confidence limits on the estimated redshift of HDF850.1 in the light of the new radio detection, and find z = 4.1±0.5. We have also determined the scale-length, and hence estimated the mass of the apparently nearby (0.5arcsec distant) z ≃ 1 elliptical galaxy 3-586.0. From this we calculate that the flux density of HDF850.1 has been boosted by a factor of ≃ 3 through gravitational lensing by this intervening elliptical, consistent with predictions that a small but significant fraction of blank-field sub-mm sources are lensed by foreground galaxies. We discuss the wider implications of these results for the sub-mm population and cosmic star-formation history.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z1., A&A Herschel Special Issue, in press as a Letter. 5 pages

We examine the physical conditions of the molecular gas in the extremely red object (ERO), HR10 at z=1.44, as constrained by the recent detection of CO J+1 -> J, J=1, 4 lines. The line ratios of such widely spaced CO transitions are extremely sensitive probes of gas excitation; in the case of the HR10, they provide a rare opportunity to study the state of the molecular gas in a gas-rich starburst at a cosmologically significant distance. Contrary to earlier claims, we find the CO J=5-4 transition to be of low excitation. The brightness temperature ratio of (5-4)/(2-1) ~ 0.15. This has serious consequences for estimates of H2 gas mass and the detectability of high-J CO lines from similar objects at high redshifts (z>~3).

We quantify the fraction of the cosmic infrared background (CIB) that originates from galaxies identified in the UV/optical/near-infrared by stacking 81,250 (~35.7 arcmin^2) K-selected sources (K_AB < 24.0), split according to their rest-frame U - V vs. V - J colors into 72,216 star-forming and 9,034 quiescent galaxies, on maps from Spitzer/MIPS (24um), Herschel/SPIRE (250, 350, 500um), Herschel/PACS (100, 160um), and AzTEC (1100um). The fraction of the CIB resolved by our catalog is (69 $\pm$ 15)% at 24um, (78 $\pm$ 17)% at 70um, (58 $\pm$ 13)% at 100um, (78 $\pm$ 18)% at 160um, (80 $\pm$ 17)% at 250um, (69 $\pm$ 14)% at 350um, (65 $\pm$ 12)% at 500um, and (45 $\pm$ 8)% at 1100um. Of that total, about 95% originates from star-forming galaxies, while the remaining 5% is from apparently quiescent galaxies. The CIB at $��$ < 200um is sourced predominantly from galaxies at z < 1, while at $��$ > 200um the bulk originates from 1 < z < 2. Galaxies with stellar masses log(M/ M_sun)=9.5-11 are responsible for the majority of the CIB, with those in the log(M/ M_sun)=9.5-10 contributing mostly at $��$ < 250um, and those in the log(M/ M_sun)=10.5-11 bin dominating at $��$ > 350um. The contribution from galaxies in the log(M/ M_sun)=9.0-9.5 and log(M/ M_sun)=11.0-12.0 stellar mass bins contribute the least, both of order 5%, although the highest stellar-mass bin is a significant contributor to the luminosity density at z > 2. The luminosities of the galaxies responsible for the CIB shifts from a combination of "normal" and luminous infrared galaxies (LIRGs) at $��$ < 160um, to LIRGs at 160um < $��$ < 500um, to finally LIRGs and ultra-luminous infrared galaxies (ULIRGs) at $��$ > 500um. Stacking analyses were performed with SIMSTACK (available at http://www.astro.caltech.edu/~viero/viero_homepage/toolbox.html) which accounts for possible biases due to clustering., 24 pages, 12 Figures, and 8 Tables. Accepted for publication in The Astrophysical Journal. Tabulated data and Software available at http://www.astro.caltech.edu/~viero/viero_homepage/toolbox.html

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 03x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs., Comment: Astronomy and Astrophysics, vol.533, A119 (version with updated figures 2 and 3). [26 pages, 25 figures]

We report the detection of a significant excess in the surface density of far-infrared sources from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) within ∼ 1 Mpc of the centres of 66 optically-selected clusters of galaxie s in the SDSS with hzi ∼0.25. From the analysis of the multiwavelength properties of their cou nterparts we conclude that the farinfrared emission is associated with dust-obscured star fo rmation and/or active galactic nuclei within galaxies in the clusters themselves. The excess reaches a maximum at a radius of ∼ 0.8 Mpc, where we find 1.0±0.3 S250 >34 mJy sources on average per cluster above what would be expected for random field locations. If the far-infr ared emission is dominated by star formation (as opposed to AGN) then this corresponds to an average star formation rate of ∼7 M⊙ yr −1 per cluster in sources with LIR > 5×10 10 L⊙. Although lensed sources make a negligible contribution to the excess signal, a fraction of the sources around the clusters could be gravitationally lensed, and we have identified a sample of potential cases of cluster-lensed Herschel sources that could be targeted in follow-up studies.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2., Comment: Accepted for publication in Astronomy & Astrophyics (Herschel first results special issue)