Search

We report on deep SCUBA-2 observations at 850$\mu$m and NOEMA spectroscopic measurements at 2 mm of the environment surrounding the luminous, massive ($M_{*} \approx 2 \times 10^{11}$ M$_{\odot}$) Herschel-selected source HerBS-70. This source was revealed by previous NOEMA observations to be a binary system of dusty star-forming galaxies at $z= 2.3$, with the East component (HerBS-70E) hosting an Active Galactic Nucleus (AGN). The SCUBA-2 observations detected, in addition to the binary system, twenty-one sources at $> 3.5 \sigma$ over an area of $\sim 25$ square comoving Mpc with a sensitivity of $\sigma_{850} = 0.75$ mJy. The surface density of continuum sources around HerBS-70 is three times higher than for field galaxies. The NOEMA spectroscopic measurements confirm the protocluster membership of three of the nine brightest sources through their CO(4 - 3) line emission, yielding a volume density 36 times higher than for field galaxies. All five confirmed sub-mm galaxies in the HerBS-70 system have relatively short gas depletion times ($80 - 500$ Myr), indicating the onset of quenching for this protocluster core due to the depletion of gas. The dark matter halo mass of the HerBS-70 system is estimated around $5 \times{} 10^{13}$ M$_{\odot}$, with a projected current-day mass of $10^{15}$ M$_{\odot}$, similar to the local Virgo and Coma clusters. These observations support the claim that DSFGs, in particular the ones with observed multiplicity, can trace cosmic overdensities., Comment: 19 pages, 13 figures, accepted for publication in MNRAS

Magnetic fields are fundamental to the evolution of galaxies, playing a key role in the astrophysics of the interstellar medium and star formation. Large-scale ordered magnetic fields have been mapped in the Milky Way and nearby galaxies, but it is not known how early in the Universe such structures form. Here we report the detection of linearly polarized thermal emission from dust grains in a strongly lensed, intrinsically luminous galaxy that is forming stars at a rate more than a thousand times that of the Milky Way at redshift 2.6, within 2.5 Gyr of the Big Bang. The polarized emission arises from the alignment of dust grains with the local magnetic field. The median polarization fraction is of order one per cent, similar to nearby spiral galaxies. Our observations support the presence of a 5 kiloparsec-scale ordered magnetic field with a strength of around 500uG or lower, orientated parallel to the molecular gas disk. This confirms that such structures can be rapidly formed in galaxies, early in cosmic history., Comment: Published in Nature. Online version available at https://www.nature.com/articles/s41586-023-06346-4

The z-GAL survey observed 137 bright Herschel-selected targets with the IRAM NOrthern Extended Millimeter Array, with the aim to measure their redshift and study their properties. Several of them have been resolved into multiple sources. Consequently, robust spectroscopic redshifts have been measured for 165 individual galaxies in the range 0.8

(Abridged) We present the dust properties of 125 bright Herschel galaxies selected from the z-GAL survey. The large instantaneous bandwidth of NOEMA provides an exquisite sampling of the underlying dust continuum emission at 2 and 3 mm in the observed frame, with flux densities in at least four side bands for each source. Together with the available Herschel 250, 350, and 500 micron and SCUBA-2 850 micron flux densities, the spectral energy distribution of each source can be analyzed from the far-infrared to the millimeter, with a fine sampling of the Rayleigh-Jeans tail. This wealth of data provides a solid basis to derive robust dust properties, in particular the dust emissivity index, beta, and the dust temperature, T(dust). In order to demonstrate our ability to constrain the dust properties, we used a flux-generated mock catalog and analyzed the results under the assumption of an optically thin and optically thick modified black body emission. For the z-GAL sources, we report a range of dust emissivities with beta ~ 1.5 - 3 estimated up to high precision with relative uncertainties that vary in the range 7% - 15%, and an average of 2.2 +/- 0.3. We find dust temperatures varying from 20 to 50 K with an average of T(dust) ~ 30 K for the optically thin case and ~38 K in the optically thick case. For all the sources, we estimate the dust masses and apparent infrared luminosities (based on the optically thin approach). An inverse correlation is found between T(dust) and beta, which is similar to what is seen in the local Universe. Finally, we report an increasing trend in the dust temperature as a function of redshift at a rate of 6.5 +/- 0.5 K/z for this 500 micron-selected sample. Based on this study, future prospects are outlined to further explore the evolution of dust temperature across cosmic time., Comment: Accepted for publication on A&A; 32 pages; 26 figures

(Abridged) Using the IRAM NOEMA interferometer, we measures the redshifts of 126 bright galaxies detected in the Herschel H-ATLAS, HeLMS, and HerS surveys. We report reliable spectroscopic redshifts for a total of 124 of the Herschel-selected galaxies. The redshifts are estimated from scans of the 3 and 2-mm bands (and, in one case, the 1-mm band) and are based on the detection of at least two emission lines. Together with the Pilot Programme (Neri et al. 2020), including spectroscopic redshifts of 11 sources, our survey has derived precise redshifts for 135 bright Herschel-selected galaxies, making it the largest sample of high-z galaxies with robust redshifts to date. Most emission lines detected are from 12CO (mainly from J=2-1 to 5-4), with some sources seen in [CI] and H2O emission lines. The spectroscopic redshifts are in the range 0.8

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 present the results of a survey of CO(1-0) emission in 14 infrared luminous dusty star forming galaxies (DSFGs) at 2 < z < 4 with the NSF's Karl G. Jansky Very Large Array. All sources are detected in CO(1-0), with an ~1arcsec angular resolution. Seven sources show extended and complex structure. We measure CO luminosities of $({\mu})L'_{CO(1-0)}=0.4-2.9x10^{11}$ K km s$^{-1}$ pc$^2$, and molecular gas masses of (${\mu}$)M$_{H2}$ = 1.3 - 8.6 x 10$^{11}$ Mo, where ({\mu}) is the magnification factor. The derived molecular gas depletion times of t$_{\rm dep}$ = 40 - 460 Myr, cover the expected range of both normal star forming galaxies and starbursts. Comparing to the higher-J CO transitions previously observed for the same sources, we find CO temperature brightness ratios of r$_{32/10}$ = 0.4 - 1.4, r$_{43/10}$ = 0.4 - 1.7, and r$_{54/10}$ = 0.3 - 1.3. We find a wide range of CO spectral line energy distributions (SLEDs), in agreement with other high-z DSFGs, with the exception of three sources that are most comparable to the Cloverleaf and APM08279+5255. Based on radiative transfer modelling of the CO SLEDs we determine densities of n$_{H2}$ = 0.3 - 8.5 x 10$^3$ cm$^{-3}$ and temperatures of T$_K$ = 100 - 200 K. Lastly, four sources are detected in the continuum, three have radio emission consistent with their infrared derived star formation rates, while HerBS-70E requires an additional synchrotron radiation component from an active galactic nucleus. Overall, we find that even though the sample is similarly luminous in the infrared, by tracing the CO(1-0) emission a diversity of galaxy and excitation properties are revealed, demonstrating the importance of CO(1-0) observations in combination to higher-J transitions., Comment: 25pages , 13 figures, Accepted for publication at ApJ

We present 101 and 151 GHz ALMA continuum images for 85 fields selected from Herschel observations that have 500 micron flux densities >80 mJy and 250-500 micron colours consistent with z > 2, most of which are expected to be gravitationally lensed or hyperluminous infrared galaxies. Approximately half of the Herschel 500 micron sources were resolved into multiple ALMA sources, but 11 of the 15 brightest 500 micron Herschel sources correspond to individual ALMA sources. For the 37 fields containing either a single source with a spectroscopic redshift or two sources with the same spectroscopic redshift, we examined the colour temperatures and dust emissivity indices. The colour temperatures only vary weakly with redshift and are statistically consistent with no redshift-dependent temperature variations, which generally corresponds to results from other samples selected in far-infrared, submillimetre, or millimetre bands but not to results from samples selected in optical or near-infrared bands. The dust emissivity indices, with very few exceptions, are largely consistent with a value of 2. We also compared spectroscopic redshifts to photometric redshifts based on spectral energy distribution templates designed for infrared-bright high-redshift galaxies. While the templates systematically underestimate the redshifts by ~15%, the inclusion of ALMA data decreases the scatter in the predicted redshifts by a factor of ~2, illustrating the potential usefulness of these millimetre data for estimating photometric redshifts., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society

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 present spectroscopic measurements for 71 galaxies associated with 62 of the brightest high-redshift submillimeter sources from the Southern fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), while targeting 85 sources which resolved into 142. We have obtained robust redshift measurements for all sources using the 12-m Array and an efficient tuning of ALMA to optimise its use as a redshift hunter, with 73 per cent of the sources having a robust redshift identification. Nine of these redshift identifications also rely on observations from the Atacama Compact Array. The spectroscopic redshifts span a range $1.41

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., Comment: This article has been accepted for publication in MNRAS, published by Oxford University Press on behalf of the Royal Astronomical Society. 24 pages

We present new observations with the Atacama Large Millimeter/sub-millimeter Array of the 122um and 205um fine-structure line emission of singly-ionised nitrogen in a strongly lensed starburst galaxy at z=2.6. The 122/205um [NII] line ratio is sensitive to electron density, n_e, in the ionised interstellar medium, and we use this to measure n_e~300cm^-3 averaged across the galaxy. This is over an order of magnitude higher than the Milky Way average, but comparable to localised Galactic star-forming regions. Combined with observations of the atomic carbon (CI(1-0)) and carbon monoxide (CO(4-3)) in the same system, we reveal the conditions in this intensely star-forming system. The majority of the molecular interstellar medium has been driven to high density, and the resultant conflagration of star formation produces a correspondingly dense ionised phase, presumably co-located with myriad HII regions that litter the gas-rich disk., Comment: 9 pages, 4 figures. Accepted for publication in ApJ

Using the IRAM NOrthern Extended Millimeter Array (NOEMA), we conducted a program to measure redshifts for 13 bright galaxies detected in the Herschel Astrophysical Large Area Survey (H-ATLAS) with $S_{500{\mu}\rm m}\ge$80 mJy. We report reliable spectroscopic redshifts for 12 individual sources, which are derived from scans of the 3 and 2 mm bands, covering up to 31 GHz in each band, and are based on the detection of at least two emission lines. The spectroscopic redshifts are in the range $2.08 10^{13}\,L_\odot$ galaxies. We also present a reanalysis of the spectral energy distributions including the continuum flux densities measured at 3 and 2 mm to derive the overall properties of the sources. Future prospects based on these efficient measurements of redshifts of high-z galaxies using NOEMA are outlined, including a comprehensive survey of all the brightest Herschel galaxies., Comment: 18 pages, 21 figures

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

We present Atacama Large Millimeter/submillimeter Array observations of a radio-loud and millimeter-bright galaxy at z=2.6. Gravitational lensing by a foreground galaxy at z~0.2 provides access to physical scales of approximately 360 pc, and we resolve a 2.5 kpc-radius ring of star-forming molecular gas, traced by atomic carbon CI(1-0) and carbon monoxide CO(4-3). We also detect emission from the cyanide radical, CN(4-3). With a velocity width of 680 km/s, this traces dense molecular gas travelling at velocities nearly a factor of two larger than the rotation speed of the molecular ring. While this could indicate the presence of a dynamical and photochemical interaction between the active galactic nucleus and molecular interstellar medium on scales of a few 100 pc, on-going feedback is unlikely to have a significant impact on the assembly of stellar mass in the molecular ring, given the ~10s Myr depletion timescale due to star formation., Comment: Accepted in ApJ Letters: The Astrophysical Journal Letters, Volume 866, Number 1

With the advent of wide-area submillimeter surveys, a large number of high-redshift gravitationally lensed dusty star-forming galaxies (DSFGs) has been revealed. Due to the simplicity of the selection criteria for candidate lensed sources in such surveys, identified as those with $S_{500\mu m} > 100$ mJy, uncertainties associated with the modelling of the selection function are expunged. The combination of these attributes makes submillimeter surveys ideal for the study of strong lens statistics. We carried out a pilot study of the lensing statistics of submillimetre-selected sources by making observations with the Atacama Large Millimetre Array (ALMA) of a sample of strongly-lensed sources selected from surveys carried out with the Herschel Space Observatory. We attempted to reproduce the distribution of image separations for the lensed sources using a halo mass function taken from a numerical simulation which contains both dark matter and baryons. We used three different density distributions, one based on analytical fits to the halos formed in the EAGLE simulation and two density distributions (Singular Isothermal Sphere (SIS) and SISSA) that have been used before in lensing studies. We found that we could reproduce the observed distribution with all three density distributions, as long as we imposed an upper mass transition of $\sim$$10^{13} M_{\odot}$ for the SIS and SISSA models, above which we assumed that the density distribution could be represented by an NFW profile. We show that we would need a sample of $\sim$500 lensed sources to distinguish between the density distributions, which is practical given the predicted number of lensed sources in the Herschel surveys., Comment: 14 pages, 7 figures, accepted for publication in MNRAS

We perform lens modelling and source reconstruction of Submillimeter Array (SMA) data for a sample of 12 strongly lensed galaxies selected at 500$\mu$m in the Herschel Astrophysical Terahertz Large Area Survey H-ATLAS. A previous analysis of the same dataset used a single S\`ersic profile to model the light distribution of each background galaxy. Here we model the source brightness distribution with an adaptive pixel scale scheme, extended to work in the Fourier visibility space of interferometry. We also present new SMA observations for seven other candidate lensed galaxies from the H-ATLAS sample. Our derived lens model parameters are in general consistent with previous findings. However, our estimated magnification factors, ranging from 3 to 10, are lower. The discrepancies are observed in particular where the reconstructed source hints at the presence of multiple knots of emission. We define an effective radius of the reconstructed sources based on the area in the source plane where emission is detected above 5$\sigma$. We also fit the reconstructed source surface brightness with an elliptical Gaussian model. We derive a median value $r_{eff}\,\sim 1.77\,$kpc and a median Gaussian full width at half maximum $\sim1.47\,$kpc. After correction for magnification, our sources have intrinsic star formation rates SFR$\,\sim900-3500\,M_{\odot}yr^{-1}$, resulting in a median star formation rate surface density $\Sigma_{SFR}\sim132\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ (or $\sim 218\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ for the Gaussian fit). This is consistent with what observed for other star forming galaxies at similar redshifts, and is significantly below the Eddington limit for a radiation pressure regulated starburst., Comment: 19 pages, 8 figures, accepted for publication in MNRAS

The {\it Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 660 deg$^2$ with the PACS and SPIRE cameras in five photometric bands: 100, 160, 250, 350 and 500\mic. This is the second of three papers describing the data release for the large fields at the south and north Galactic poles (NGP and SGP). In this paper we describe the catalogues of far-infrared and submillimetre sources for the NGP and SGP, which cover 177 deg$^2$ and 303 deg$^2$, respectively. The catalogues contain 153,367 sources for the NGP field and 193,527 sources for the SGP field detected at more than 4$\sigma$ significance in any of the 250, 350 or 500\mic\ bands. The source detection is based on the 250\mic\ map, and we present photometry in all five bands for each source, including aperture photometry for sources known to be extended. The rms positional accuracy for the faintest sources is about 2.4 arc seconds in both right ascension and declination. We present a statistical analysis of the catalogues and discuss the practical issues -- completeness, reliability, flux boosting, accuracy of positions, accuracy of flux measurements -- necessary to use the catalogues for astronomical projects.

By determining the nature of all the Planck compact sources within 808.4 deg^2 of large Herschel surveys, we have identified 27 candidate proto-clusters of dusty star forming galaxies (DSFGs) that are at least 3{\sigma} overdense in either 250, 350 or 500 $\mu$mm sources. We find roughly half of all the Planck compact sources are resolved by Herschel into multiple discrete objects, with the other half remaining unresolved by Herschel. We find a significant difference between versions of the Planck catalogues, with earlier releases hosting a larger fraction of candidate proto-clusters and Galactic Cirrus than later releases, which we ascribe to a difference in the filters used in the creation of the three catalogues. We find a surface density of DSFG candidate proto-clusters of (3.3 $\pm$ 0.7) x 10^(-2) sources deg^(-2), in good agreement with previous similar studies. We find that a Planck colour selection of S_{857}/S_{545} < 2 works well to select candidate proto-clusters, but can miss proto-clusters at z < 2. The Herschel colours of individual candidate proto-cluster members indicate our candidate proto-clusters all likely all lie at z > 1. Our candidate proto-clusters are a factor of 5 times brighter at 353 GHz than expected from simulations, even in the most conservative estimates. Further observations are needed to confirm whether these candidate proto-clusters are physical clusters, multiple proto-clusters along the line of sight, or chance alignments of unassociated sources., Comment: Resubmitted to MNRAS after request for minor revision

This paper forms part of the second major public data release of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). In this work, we describe the identification of optical and near-infrared counterparts to the submillimetre detected sources in the $177$ deg$^2$ North Galactic Plane (NGP) field. We used the likelihood ratio method to identify counterparts in the Sloan Digital Sky Survey and in the UKIRT Imaging Deep Sky Survey within a search radius of $10$ arcsec of the H-ATLAS sources with a $4\sigma$ detection at $250 \, \mu$m. We obtained reliable ($R \ge 0.8 $) optical counterparts with $r< 22.4$ for 42429 H-ATLAS sources ($37.8$ per cent), with an estimated completeness of $71.7$ per cent and a false identification rate of $4.7$ per cent. We also identified counterparts in the near-infrared using deeper $K$-band data which covers a smaller $\sim25$ deg$^2$. We found reliable near-infrared counterparts to $61.8$ per cent of the $250$-$\mu$m-selected sources within that area. We assessed the performance of the likelihood ratio method to identify optical and near-infrared counterparts taking into account the depth and area of both input catalogues. Using catalogues with the same surface density of objects in the overlapping $\sim25$ deg$^2$ area, we obtained that the reliable fraction in the near-infrared ($54.8$ per cent) is significantly higher than in the optical ($36.4$ per cent). Finally, using deep radio data which covers a small region of the NGP field, we found that $80 - 90$ per cent of our reliable identifications are correct., Comment: 19 pages, 24 figures, submitted to MNRAS

We present the Herschel Bright Sources (HerBS) sample, a sample of bright, high-redshift Herschel sources detected in the 616.4 square degree H-ATLAS survey. The HerBS sample contains 209 galaxies, selected with a 500 {\mu}m flux density greater than 80 mJy and an estimated redshift greater than 2. The sample consists of a combination of HyLIRGs and lensed ULIRGs during the epoch of peak cosmic star formation. In this paper, we present SCUBA-2 observations at 850 ${\mu}$m of 189 galaxies of the HerBS sample, 152 of these sources were detected. We fit a spectral template to the Herschel-SPIRE and 850 ${\mu}$m SCUBA-2 flux densities of 22 sources with spectroscopically determined redshifts, using a two-component modified blackbody spectrum as a template. We find a cold- and hot-dust temperature of 21.29 K and 45.80 K, a cold-to-hot dust mass ratio of 26.62 and a $\beta$ of 1.83. The poor quality of the fit suggests that the sample of galaxies is too diverse to be explained by our simple model. Comparison of our sample to a galaxy evolution model indicates that the fraction of lenses is high. Out of the 152 SCUBA-2 detected galaxies, the model predicts 128.4 $\pm$ 2.1 of those galaxies to be lensed (84.5%). The SPIRE 500 ${\mu}$m flux suggests that out of all 209 HerBS sources, we expect 158.1 $\pm$ 1.7 lensed sources, giving a total lensing fraction of 76 per cent.

This paper defines the UK Infra-Red Telescope (UKIRT) Hemisphere Survey (UHS) and release of the remaining ~12,700 sq.deg of J-band survey data products. The UHS will provide continuous J and K-band coverage in the northern hemisphere from a declination of 0 deg to 60 deg by combining the existing Large Area Survey, Galactic Plane Survey and Galactic Clusters Survey conducted under the UKIRT Infra-red Deep Sky Survey (UKIDSS) programme with this new additional area not covered by UKIDSS. The released data includes J-band imaging and source catalogues over the new area, which, together with UKIDSS, completes the J-band UHS coverage over the full ~17,900 sq.deg area. 98 per cent of the data in this release have passed quality control criteria, the remaining 2 per cent being scheduled for re-observation. The median 5-sigma point source sensitivity of the released data is 19.6 mag (Vega). The median full width at half-maximum of the point spread function across the dataset is 0.75 arcsec. In this paper, we outline the survey management, data acquisition, processing and calibration, quality control and archiving as well as summarising the characteristics of the released data products. The data are initially available to a limited consortium with a world-wide release scheduled for August 2018., Comment: 14 pages, 9 figures, MNRAS in press. v2 includes some changes suggested by referee

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

In this work we measure and study the cross-correlation signal between a foreground sample of GAMA galaxies with spectroscopic redshifts in the range $0.2

Combining large-area optical quasar surveys with the new far-infrared Herschel-ATLAS Data Release 1, we search for an observational signature associated with the minority of quasars possessing bright far-infrared (FIR) luminosities. We find that FIR-bright quasars show broad CIV emission line blueshifts in excess of that expected from the optical luminosity alone, indicating particularly powerful nuclear outflows. The quasars show no signs of having redder optical colours than the general ensemble of optically-selected quasars, ruling out differences in line-of-sight dust within the host galaxies. We postulate that these objects may be caught in a special evolutionary phase, with unobscured, high black hole accretion rates and correspondingly strong nuclear outflows. The high FIR emission found in these objects is then either a result of star formation related to the outflow, or is due to dust within the host galaxy illuminated by the quasar. We are thus directly witnessing coincident small-scale nuclear processes and galaxy-wide activity, commonly invoked in galaxy simulations which rely on feedback from quasars to influence galaxy evolution., Comment: Comments: 7 pages, 3 figures, accepted for publication in MNRAS

We present an extragalactic survey using observations from the Atacama Large Millimeter/submillimeter Array (ALMA) to characterise galaxy populations up to $z=0.35$: the Valpara\'iso ALMA Line Emission Survey (VALES). We use ALMA Band-3 CO(1--0) observations to study the molecular gas content in a sample of 67 dusty normal star-forming galaxies selected from the $Herschel$ Astrophysical Terahertz Large Area Survey ($H$-ATLAS). We have spectrally detected 49 galaxies at $>5\sigma$ significance and 12 others are seen at low significance in stacked spectra. CO luminosities are in the range of $(0.03-1.31)\times10^{10}$ K km s$^{-1}$ pc$^2$, equivalent to $\log({\rm M_{gas}/M_{\odot}}) =8.9-10.9$ assuming an $\alpha_{\rm CO}$=4.6(K km s$^{-1}$ pc$^{2}$)$^{-1}$, which perfectly complements the parameter space previously explored with local and high-z normal galaxies. We compute the optical to CO size ratio for 21 galaxies resolved by ALMA at $\sim 3$."$5$ resolution (6.5 kpc), finding that the molecular gas is on average $\sim$ 0.6 times more compact than the stellar component. We obtain a global Schmidt-Kennicutt relation, given by $\log [\Sigma_{\rm SFR}/({\rm M_{\odot} yr^{-1}kpc^{-2}})]=(1.26 \pm 0.02) \times \log [\Sigma_{\rm M_{H2}}/({\rm M_{\odot}\,pc^{-2}})]-(3.6 \pm 0.2)$. We find a significant fraction of galaxies lying at `intermediate efficiencies' between a long-standing mode of star-formation activity and a starburst, specially at $\rm L_{IR}=10^{11-12} L_{\odot}$. Combining our observations with data taken from the literature, we propose that star formation efficiencies can be parameterised by $\log [{\rm SFR/M_{H2}}]=0.19 \times {\rm (\log {L_{IR}}-11.45)}-8.26-0.41 \times \arctan[-4.84 (\log {\rm L_{IR}}-11.45) ]$. Within the redshift range we explore ($z<0.35$), we identify a rapid increase of the gas content as a function of redshift., Comment: Accepted for publication in MNRAS. 31 pages, including 9 figures and 2 tables

We have modelled high resolution ALMA imaging of six strong gravitationally lensed galaxies detected by the Herschel Space Observatory. Our modelling recovers mass properties of the lensing galaxies and, by determining magnification factors, intrinsic properties of the lensed sub-millimetre sources. We find that the lensed galaxies all have high ratios of star formation rate to dust mass, consistent with or higher than the mean ratio for high redshift sub-millimetre galaxies and low redshift ultra-luminous infra-red galaxies. Source reconstruction reveals that most galaxies exhibit disturbed morphologies. Both the cleaned image plane data and the directly observed interferometric visibilities have been modelled, enabling comparison of both approaches. In the majority of cases, the recovered lens models are consistent between methods, all six having mass density profiles that are close to isothermal. However, one system with poor signal to noise shows mildly significant differences., Comment: 13 pages, 4 figures. v3 includes minor alterations following further refereeing of v2

We present the calibration between the dust continuum luminosity and interstellar gas content obtained from the Valpara\'{i}so ALMA Line Emission Survey (VALES) sample of 67 main-sequence star-forming galaxies at 0.02<$z$<0.35. We use CO(1-0) observations from the Atacama Large Millimetre/submillimetre Array (ALMA) to trace the molecular gas mass, $M_{\mathrm{H}_{2}}$, and estimate the rest-frame monochromatic luminosity at 850 $\mu$m, $L_{\nu_{850}}$, by extrapolating the dust continuum from MAGPHYS modelling of the far-ultraviolet to submillimetre spectral energy distribution sampled by the Galaxy And Mass Assembly (GAMA) survey. Adopting $\alpha_{\rm CO}$ = 6.5 (K km s$^{-1}$ pc$^{2}$)$^{-1}$, the average ratio of $L_{\nu_{850}}/M_{\mathrm{H}_{2}}$ = (6.4$\pm$1.4)$\times10^{19}$ erg s$^{-1}$ Hz$^{-1}$ $\mathrm{M}_{\odot}^{-1}$, in excellent agreement with literature values. We obtain a linear fit of $\log_{10}$ ($M_{\mathrm{H}_{2}}/\mathrm{M}_{\odot}$) = (0.92$\pm$0.02) $\log_{10}$ ($L_{\nu_{850}}$/erg s$^{-1}$ Hz$^{-1}$)-(17.31$\pm$0.59). We provide relations between $L_{\nu_{850}}$, $M_{\mathrm{H}_{2}}$ and $M_{\mathrm{ISM}}$ when combining the VALES and literature samples, and adopting a Galactic $\alpha_{\rm CO}$ value., Comment: Accepted for publication in MNRAS. 5 pages, including 2 figures and 1 table

We use new Band-3 CO(1-0) observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) to study the physical conditions in the interstellar gas of a sample of 27 dusty main-sequence star-forming galaxies at 0.03<$z$<0.2 present in the Valpara\'iso ALMA Line Emission Survey (VALES). The sample is drawn from far-IR bright galaxies over $\sim$160 deg$^{2}$ in the Herschel Astrophysical Terahertz Large Area Survey (HATLAS), which is covered by Herschel [CII] 158 $\mu$m spectroscopy and far-infrared (FIR) photometry. The [CII] and CO lines are both detected at >5$\sigma$ in 26 sources. We find an average [CII] to CO(1-0) luminosity ratio of 3500$\pm$1200 for our sample that is consistent with previous studies. Using the [CII], CO and FIR measurements as diagnostics of the physical conditions of the interstellar medium, we compare these observations to the predictions of a photodissociation region (PDR) model to determine the gas density, surface temperature, pressure, and the strength of the incident far-ultraviolet (FUV) radiation field, $G_{0}$, normalised to the Habing Field. The majority of our sample exhibit hydrogen densities of 4 < $\log n/\mathrm{cm}^{3}$ < 5.5 and experience an incident FUV radiation field with strengths of 2 < $\log G_0$ < 3 when adopting standard adjustments. A comparison to galaxy samples at different redshifts indicates that the average strength of the FUV radiation field appears constant up to redshift $z\sim$6.4, yet the neutral gas density increases with redshift by a factor of $\sim$100, that persists regardless of various adjustments to our observable quantities. This evolution could provide an explanation for the observed evolution of the star formation rate density with cosmic time, yet could arise from a combination of observational biases when using different suites of emission lines as diagnostic tracers of PDR gas., Comment: Accepted for publication in A&A. Abstract abridged for arXiv. 13 pages, including 5 figures and 2 tables

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500{\mu}m extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500{\mu}m and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies., Comment: 25 pages, 17 figures, accepted for publication in MNRAS

We present a study of the dust, stars and atomic gas (HI) in an HI-selected sample of local galaxies (z<0.035) in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) fields. This HI-selected sample reveals a population of very high gas fraction (>80 per cent), low stellar mass sources that appear to be in the earliest stages of their evolution. We compare this sample with dust and stellar mass selected samples to study the dust and gas scaling relations over a wide range of gas fraction (proxy for evolutionary state of a galaxy). The most robust scaling relations for gas and dust are those linked to NUV-r (SSFR) and gas fraction, these do not depend on sample selection or environment. At the highest gas fractions, our additional sample shows the dust content is well below expectations from extrapolating scaling relations for more evolved sources, and dust is not a good tracer of the gas content. The specific dust mass for local galaxies peaks at a gas fraction of ~75 per cent. The atomic gas depletion time is also longer for high gas fraction galaxies, opposite to the trend found for molecular gas depletion timescale. We link this trend to the changing efficiency of conversion of HI to H2 as galaxies increase in stellar mass surface density as they evolve. Finally, we show that galaxies start out barely obscured and increase in obscuration as they evolve, yet there is no clear and simple link between obscuration and global galaxy properties., Comment: 16 pages (+11 pages appendix), 12 figures, Published in 2017 in MNRAS, 464, 4680; This version 2 has corrected a small error in the table heading of Table 1, but is otherwise unchanged

(abridged) We report rest-frame submillimeter H2O emission line observations of 11 HyLIRGs/ULIRGs at z~2-4 selected among the brightest lensed galaxies discovered in the Herschel-ATLAS. Using the IRAM NOEMA, we have detected 14 new H2O emission lines. The apparent luminosities of the H2O emission lines are $\mu L_{\rm{H_2O}} \sim 6-21 \times 10^8 L_\odot$, with velocity-integrated line fluxes ranging from 4-15 Jy km s$^{-1}$. We have also observed CO emission lines using EMIR on the IRAM 30m telescope in seven sources. The velocity widths for CO and H2O lines are found to be similar. With almost comparable integrated flux densities to those of the high-J CO line, H2O is found to be among the strongest molecular emitters in high-z Hy/ULIRGs. We also confirm our previously found correlation between luminosity of H2O ($L_{\rm{H_2O}}$) and infrared ($L_{\rm{IR}}$) that $L_{\rm{H_2O}} \sim L_{\rm{IR}}^{1.1-1.2}$, with our new detections. This correlation could be explained by a dominant role of far-infrared (FIR) pumping in the H2O excitation. Modelling reveals the FIR radiation fields have warm dust temperature $T_\rm{warm}$~45-75 K, H2O column density per unit velocity interval $N_{\rm{H_2O}}/\Delta V \gtrsim 0.3 \times 10^{15}$ cm$^{-2}$ km$^{-1}$ s and 100 $\mu$m continuum opacity $\tau_{100} > 1$ (optically thick), indicating that H2O is likely to trace highly obscured warm dense gas. However, further observations of $J\geq4$ H2O lines are needed to better constrain the continuum optical depth and other physical conditions of the molecular gas and dust. We have also detected H2O+ emission in three sources. A tight correlation between $L_{\rm{H_2O}}$ and $L_{\rm{H_2O^+}}$ has been found in galaxies from low to high redshift. The velocity-integrated flux density ratio between H2O+ and H2O suggests that cosmic rays generated by strong star formation are possibly driving the H2O+ formation., Comment: Minor changes to match the published version in A&A

Aims. We aim to study the 250 micron luminosity function (LF) down to much fainter luminosities than achieved by previous efforts. Methods. We developed a modified stacking method to reconstruct the 250 micron LF using optically selected galaxies from the SDSS survey and Herschel maps of the GAMA equatorial fields and Stripe 82. Our stacking method not only recovers the mean 250 micron luminosities of galaxies that are too faint to be individually detected, but also their underlying distribution functions. Results. We find very good agreement with previous measurements in the overlapping luminosity range. More importantly, we are able to derive the LF down to much fainter luminosities (around 25 times fainter) than achieved by previous studies. We find strong positive luminosity evolution \propto (1 + z)^4.89\pm1.07 and moderate negative density evolution \propto (1 + z)^-1.02\pm0.54 over the redshift range z=[0.02, 0.5]., Comment: 5 pages including appendix, 4 figures, accepted for publication as a letter in A&A

We compare common star-formation rate (SFR) indicators in the local Universe in the GAMA equatorial fields (around 160 sq. deg.), using ultraviolet (UV) photometry from GALEX, far-infrared (FIR) and sub-millimetre (sub-mm) photometry from H-ATLAS, and Halpha spectroscopy from the GAMA survey. With a high-quality sample of 745 galaxies (median redshift 0.08), we consider three SFR tracers: UV luminosity corrected for dust attenuation using the UV spectral slope beta (SFRUV,corr), Halpha line luminosity corrected for dust using the Balmer decrement (BD) (SFRHalpha,corr), and the combination of UV and IR emission (SFRUV+IR). We demonstrate that SFRUV,corr can be reconciled with the other two tracers after applying attenuation corrections by calibrating IRX (i.e. the IR to UV luminosity ratio) and attenuation in the Halpha (derived from BD) against beta. However, beta on its own is very unlikely to be a reliable attenuation indicator. We find that attenuation correction factors depend on parameters such as stellar mass, z and dust temperature (Tdust), but not on Halpha equivalent width (EW) or Sersic index. Due to the large scatter in the IRX vs beta correlation, when compared to SFRUV+IR, the beta-corrected SFRUV,corr exhibits systematic deviations as a function of IRX, BD and Tdust., Comment: 21 pages. 19 figures. Accepted for publication in MNRAS

We present the first major data release of the largest single key-project in area carried out in open time with the Herschel Space Observatory. The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 600 deg^2 in five photometric bands - 100, 160, 250, 350 and 500 um - with the PACS and SPIRE cameras. In this paper and a companion paper (Bourne et al. 2016) we present the survey of three fields on the celestial equator, covering a total area of 161.6 deg^2 and previously observed in the Galaxy and Mass Assembly (GAMA) spectroscopic survey. This paper describes the Herschel images and catalogues of the sources detected on the SPIRE 250 um images. The 1-sigma noise for source detection, including both confusion and instrumental noise, is 7.4, 9.4 and 10.2 mJy at 250, 350 and 500 um. Our catalogue includes 120230 sources in total, with 113995, 46209 and 11011 sources detected at >4-sigma at 250, 350 and 500 um. The catalogue contains detections at >3-sigma at 100 and 160 um for 4650 and 5685 sources, and the typical noise at these wavelengths is 44 and 49 mJy. We include estimates of the completeness of the survey and of the effects of flux bias and also describe a novel method for determining the true source counts. The H-ATLAS source counts are very similar to the source counts from the deeper HerMES survey at 250 and 350 um, with a small difference at 500 um. Appendix A provides a quick start in using the released datasets, including instructions and cautions on how to use them., Comment: Accepted by MNRAS; data available at www.h-atlas.org

This paper is the second in a pair of articles presenting data release 1 (DR1) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), the largest single open-time key project carried out with the Herschel Space Observatory. The H-ATLAS is a wide-area imaging survey carried out in five photometric bands at 100, 160, 250, 350 and 500$\mu$m covering a total area of 600deg$^2$. In this paper we describe the identification of optical counterparts to submillimetre sources in DR1, comprising an area of 161 deg$^2$ over three equatorial fields of roughly 12$^\circ$x4.5$^\circ$ centred at 9$^h$, 12$^h$ and 14.5$^h$ respectively. Of all the H-ATLAS fields, the equatorial regions benefit from the greatest overlap with current multi-wavelength surveys spanning ultraviolet (UV) to mid-infrared regimes, as well as extensive spectroscopic coverage. We use a likelihood-ratio technique to identify SDSS counterparts at r<22.4 for 250-$\mu$m-selected sources detected at $\geq$ 4$\sigma$ ($\approx$28mJy). We find `reliable' counterparts (reliability R$\geq$0.8) for 44,835 sources (39 per cent), with an estimated completeness of 73.0 per cent and contamination rate of 4.7 per cent. Using redshifts and multi-wavelength photometry from GAMA and other public catalogues, we show that H-ATLAS-selected galaxies at $z<0.5$ span a wide range of optical colours, total infrared (IR) luminosities, and IR/UV ratios, with no strong disposition towards mid-IR-classified AGN in comparison with optical selection. The data described herein, together with all maps and catalogues described in the companion paper (Valiante et al. 2016), are available from the H-ATLAS website at www.h-atlas.org., Comment: MNRAS accepted

We describe a far infrared survey of the Coma cluster and the galaxy filament it resides within. Our survey covers an area of $\sim$150 deg$^2$ observed by $Herschel$ H-ATLAS in five bands at 100, 160, 250, 350 and 500$\mu$m. The SDSS spectroscopic survey ($m_{r} \le 17.8)$ is used to define an area (within the Virial radius) and redshift selected ($4268 < v < 9700$ km s$^{-1}$) sample of 744 Coma cluster galaxies - the Coma Cluster Catalogue (CCC). For comparison we also define a sample of 951 galaxies in the connecting filament - the Coma Filament Catalogue (CFC). The optical positions and parameters are used to define appropriate apertures to measure each galaxy's far infrared emission. We have detected 99 of 744 (13\%) and 422 of 951 (44\%) of the cluster and filament galaxies in the SPIRE 250$\mu$m band. We consider the relative detection rates of galaxies of different morphological types finding that it is only the S0/Sa population that shows clear differences between the cluster and filament. We find no differences between the dust masses and temperatures of cluster and filament galaxies with the exception of early type galaxy dust temperatures, which are significantly hotter in the cluster than in the filament (X-ray heating?). From a chemical evolution model we find no evidence for different evolutionary processes (gas loss or infall) between galaxies in the cluster and filament., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society

The Herschel Space Observatory has had a tremendous impact on the study of extragalactic dust. Specifically, early-type galaxies (ETG) have been the focus of several studies. In this paper we combine results from two Herschel studies - a Virgo cluster study HeViCS and a broader, low-redshift H-ATLAS/GAMA study - and contrast the dust and associated properties for similar mass galaxies. This comparison is motivated by differences in results exhibited between multiple Herschel studies of early-type galaxies. A comparison between consistent modified blackbody derived dust mass is carried out, revealing strong differences between the two samples in both dust mass and dust-to-stellar mass ratio. In particular, the HeViCS sample lacks massive ETG with as high a specific dust content as found in H-ATLAS. This is most likely connected with the difference in environment for the two samples. We calculate nearest neighbour environment densities in a consistent way, showing that H-ATLAS ETG occupy sparser regions of the local Universe, whereas HeViCS ETG occupy dense regions. This is also true for ETG that are not Herschel-detected but are in the Virgo and GAMA parent samples. Spectral energy distributions are fit to the panchromatic data. From these we find that in H-ATLAS the specific star formation rate anticorrelates with stellar mass and reaches values as high as in our Galaxy. On the other hand HeViCS ETG appear to have little star formation. Based on the trends found here, H-ATLAS ETG are thought to have more extended star formation histories and a younger stellar population than HeViCS ETG., Comment: 23 pages, 16 figures, 2 tables, accepted for publication in MNRAS

We have modelled Atacama Large Millimeter/sub-millimeter Array (ALMA) long baseline imaging of the strong gravitational lens system H-ATLAS J090311.6+003906 (SDP.81). We have reconstructed the distribution of band 6 and 7 continuum emission in the z=3.042 source and we have determined its kinematic properties by reconstructing CO(5-4) and CO(8-7) line emission in bands 4 and 6. The continuum imaging reveals a highly non-uniform distribution of dust with clumps on scales of 200 pc. In contrast, the CO line emission shows a relatively smooth, disk-like velocity field which is well fit by a rotating disk model with an inclination angle of 40+/-5 degrees and an asymptotic rotation velocity of 320 km/s. The inferred dynamical mass within 1.5kpc is (3.5+/-0.5)x10^{10} M_sol which is comparable to the total molecular gas masses of (2.7+/-0.5)x10^{10} M_sol and (3.5+/-0.6)x10^{10} M_sol from the dust continuum emission and CO emission respectively. Our new reconstruction of the lensed HST near-infrared emission shows two objects which appear to be interacting, with the rotating disk of gas and dust revealed by ALMA distinctly offset from the near-infrared emission. The clumpy nature of the dust and a low value of the Toomre parameter of Q=0.3 suggest that the disk is in a state of collapse. We estimate a star formation rate in the disk of 470+/-80 M_sol/yr with an efficiency 65 times greater than typical low-redshift galaxies. Our findings add to the growing body of evidence that the most infra-red luminous, dust obscured galaxies in the high redshift Universe represent a population of merger induced starbursts., Comment: 11 pages, 7 figures, resubmitted to MNRAS following referee's comments. Includes a correction to the molecular gas mass from CO

We explore the behaviour of [CII]-157.74um forbidden fine-structure line observed in a sample of 28 galaxies selected from ~50deg^2 of the H-ATLAS survey. The sample is restricted to galaxies with flux densities higher than S_160um>150mJy and optical spectra from the GAMA survey at 0.022.5x10^-3 with respect to those showing lower ratios. In particular, those with high ratios tend to have: (1) L_IR<10^11Lo; (2) cold dust temperatures, T_d<30K; (3) disk-like morphologies in r-band images; (4) a WISE colour 0.5) is main parameter responsible for controlling the [CII]/IR ratio. It is possible that relatively high creates a positively charged dust grain distribution, impeding an efficient photo-electric extraction of electrons from these grains to then collisionally excite carbon atoms. Within the brighter IR population, 11

We present the first measurement of the correlation between the map of the cosmic microwave background (CMB) lensing potential derived from the \emph{Planck} nominal mission data and $z\gtrsim 1.5$ galaxies detected by the \emph{Herschel}-ATLAS (H-ATLAS) survey covering about $600\,\hbox{deg}^2$, i.e. about 1.4\% of the sky. We reject the hypothesis that there is no correlation between CMB lensing and galaxy detection at a $20\,\sigma$ significance, checking the result by performing a number of null tests. The significance of the detection of the theoretically expected cross-correlation signal is found to be $10\,\sigma$. The galaxy bias parameter, $b$, derived from a joint analysis of the cross-power spectrum and of the auto-power spectrum of the galaxy density contrast is found to be $b=2.80^{+0.12}_{-0.11}$, consistent with earlier estimates for H-ATLAS galaxies at similar redshifts. On the other hand, the amplitude of the cross-correlation is found to be a factor $1.62 \pm 0.16$ higher than expected from the standard model and also found by cross-correlation analyses with other tracers of the large-scale structure. The enhancement due to lensing magnification can account for only a fraction of the excess cross-correlation signal. We suggest that part of it may be due to an incomplete removal of the contamination of the CIB, that includes the H-ATLAS sources we are cross-correlating with. In any case, the highly significant detection reported here using a catalog covering only 1.4\% of the sky demonstrates the potential of CMB lensing correlations with submillimeter surveys., Comment: 14 pages, 17 figures, 3 tables, accepted for pubblication in ApJ, minor changes to match the published version

We report an unexpected variation in the positional offset distributions between Herschel-ATLAS sub-millimetre (submm) sources and their optical associations, depending on both 250-{\mu}m signal-to-noise ratio and 250/350-{\mu}m colour. We show that redder and brighter submm sources have optical associations with a broader distribution of positional offsets than would be expected if these offsets were due to random positional errors in the source extraction. The observation can be explained by two possible effects: either red submm sources trace a more clustered population than blue ones, and their positional errors are increased by confusion; or red submm sources are generally at high redshifts and are frequently associated with low-redshift lensing structures which are identified as false counterparts. We perform various analyses of the data, including the multiplicity of optical associations, the redshift and magnitude distributions in H-ATLAS in comparison to HerMES, and simulations of weak lensing, and we conclude that the effects are most likely to be explained by widespread weak lensing of Herschel-SPIRE sources by foreground structures. This has important consequences for counterpart identification and derived redshift distributions and luminosity functions of submm surveys., Comment: 10 pages, 6 figures. Accepted for publication in MNRAS

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

We present a comparison of the physical properties of a rest-frame $250\mu$m selected sample of massive, dusty galaxies from $01$ SMGs have an average SFR of $390^{+80}_{-70}\,$M$_\odot$yr$^{-1}$ which is 120 times that of the low-redshift sample matched in stellar mass to the SMGs (SFR$=3.3\pm{0.2}$ M$_\odot$yr$^{-1}$). The SMGs harbour a substantial mass of dust ($1.2^{+0.3}_{-0.2}\times{10}^9\,$M$_\odot$), compared to $(1.6\pm0.1)\times{10}^8\,$M$_\odot$ for low-redshift dusty galaxies. At low redshifts the dust luminosity is dominated by the diffuse ISM, whereas a large fraction of the dust luminosity in SMGs originates from star-forming regions. At the same dust mass SMGs are offset towards a higher SFR compared to the low-redshift H-ATLAS galaxies. This is not only due to the higher gas fraction in SMGs but also because they are undergoing a more efficient mode of star formation, which is consistent with their bursty star-formation histories. The offset in SFR between SMGs and low-redshift galaxies is similar to that found in CO studies, suggesting that dust mass is as good a tracer of molecular gas as CO., Comment: 17 pages, 8 pages in Appendix, 12 figures, 2 tables. Accepted for publication in MNRAS. Minor typos corrected to match published version

We apply a chemical evolution model to investigate the sources and evolution of dust in a sample of 26 high-redshift ($z>1$) submillimetre galaxies (SMGs) from the literature, with complete photometry from ultraviolet to the submillimetre. We show that dust produced only by low-intermediate mass stars falls a factor 240 short of the observed dust masses of SMGs, the well-known `dust-budget crisis'. Adding an extra source of dust from supernovae can account for the dust mass in 19 per cent of the SMG sample. Even after accounting for dust produced by supernovae the remaining deficit in the dust mass budget provides support for higher supernova yields, substantial grain growth in the interstellar medium or a top-heavy IMF. Including efficient destruction of dust by supernova shocks increases the tension between our model and observed SMG dust masses. The models which best reproduce the physical properties of SMGs have a rapid build-up of dust from both stellar and interstellar sources and minimal dust destruction. Alternatively, invoking a top-heavy IMF or significant changes in the dust grain properties can solve the dust budget crisis only if dust is produced by both low mass stars and supernovae and is not efficiently destroyed by supernova shocks., Comment: 21 pages including 1 page appendix, 5 figures, 6 tables. Accepted for publication in MNRAS

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S_500>100 mJy, 21 are strongly lensed (multiply imaged), 4 are moderately lensed (singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r_half) and far-infrared luminosities (L_FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z_lens>0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500um flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L_FIR (median L_FIR=7.9x10^12 L_sun) and two decades in FIR luminosity surface density (median Sigma_FIR=6.0x10^11 L_sun kpc^-2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift., Comment: 28 pages, 7 tables, 9 figures. Accepted to ApJ; revised following minor comments from the referee and the community at large

We present an analysis of CO molecular gas tracers in a sample of 500{\mu}m-selected Herschel-ATLAS galaxies at z<0.05 (cz<14990km/s). Using 22-500{\mu}m photometry from WISE, IRAS and Herschel, with HI data from the literature, we investigate correlations between warm and cold dust, and tracers of the gas in different phases. The correlation between global CO(3-2) line fluxes and FIR-submillimetre fluxes weakens with increasing IR wavelength ({\lambda}>60{\mu}m), as a result of colder dust being less strongly associated with dense gas. Conversely, CO(2-1) and HI line fluxes both appear to be better correlated with longer wavelengths, suggesting that cold dust is more strongly associated with diffuse atomic and molecular gas phases, consistent with it being at least partially heated by radiation from old stellar populations. The increased scatter at long wavelengths implies that submillimetre fluxes are a poorer tracer of SFR. Fluxes at 22 and 60{\mu}m are also better correlated with diffuse gas tracers than dense CO(3-2), probably due to very-small-grain emission in the diffuse interstellar medium, which is not correlated with SFR. The FIR/CO luminosity ratio and the dust mass/CO luminosity ratio both decrease with increasing luminosity, as a result of either correlations between mass and metallicity (changing CO/H2) or between CO luminosity and excitation [changing CO(3-2)/CO(1-0)]., Comment: 27 pages, 11 figures. MNRAS accepted

We present high-spatial resolution imaging obtained with the Submillimeter Array (SMA) at 880um and the Keck Adaptive Optics (AO) system at Ks-band of a gravitationally lensed sub-millimeter galaxy (SMG) at z=4.243 discovered in the Herschel-Astrophysical Terahertz Large Area Survey. The SMA data (angular resolution ~0.6") resolve the dust emission into multiple lensed images, while the Keck AO Ks-band data (angular resolution ~0.1") resolve the lens into a pair of galaxies separated by 0.3". We present an optical spectrum of the foreground lens obtained with the Gemini-South telescope that provides a lens redshift of z_lens = 0.595 +/- 0.005. We develop and apply a new lens modeling technique in the visibility plane that shows that the SMG is magnified by a factor of mu = 4.1 +/- 0.2 and has an intrinsic infrared (IR) luminosity of L_IR = (2.1 +/- 0.2) x 10^13 Lsun. We measure a half-light radius of the background source of r_s = 4.4 +/- 0.5 kpc which implies an IR luminosity surface density of Sigma_IR = (3.4 +/- 0.9) x 10^11 Lsun kpc^-2, a value that is typical of z > 2 SMGs but significantly lower than IR luminous galaxies at z~0. The two lens galaxies are compact (r_lens ~ 0.9 kpc) early-types with Einstein radii of theta_E1 = 0.57 +/- 0.01 and theta_E2 = 0.40 +/- 0.01 that imply masses of M_lens1 = (7.4 +/- 0.5) x 10^10 Msun and M_lens2 = (3.7 +/- 0.3) x 10^10 Msun. The two lensing galaxies are likely about to undergo a dissipationless merger, and the mass and size of the resultant system should be similar to other early-type galaxies at z~0.6. This work highlights the importance of high spatial resolution imaging in developing models of strongly lensed galaxies discovered by Herschel., Comment: Accepted to the ApJ. 12 pages, 6 figures, 3 tables

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 study the Spectral Energy Distribution (SED) and the power spectrum of Galactic cirrus emission observed in the 14 deg^2 Science Demonstration Phase field of the Herschel-ATLAS using Herschel and IRAS data from 100 to 500 um. We compare the SPIRE 250, 350 and 500um maps with IRAS 100um emission, binned in 6' pixels. We assume a modified black-body SED with dust emissivity parameter beta (F ~ lambda^(-beta)) and a single dust temperature T_d, and find that the dust temperature and emissivity index varies over the science demonstration field as 10< T_rm < 25 K and 1 < beta< 4. The latter values are somewhat higher than the range of beta often quoted in the literature (1< beta< 2). We estimate the mean values of these parameters to be T_d=19.0 +/- 2.4 K and beta = 1.4 +/- 0.4. In regions of bright cirrus emission, we find that the dust has similar temperatures with T_d = 18.0 +/- 2.5 K, and similar values of beta, ranging from 1.4 +- 0.5 to 1.9+/- 0.5. We show that T_d and beta associated with diffuse cirrus emission are anti-correlated and can be described by the relationship: beta(T_d) = NT_d^alpha with [N=116+/-38, alpha=-1.4+/1 0.1]. The strong correlation found in this analysis is not just limited to high density clumps of cirrus emission as seen in previous studies, but is also seen in diffuse cirrus in low density regions. To provide an independent measure of $T_{\rm d}$ and $\beta$, we obtain the angular power spectrum of the cirrus emission in the {\it IRAS} and SPIRE maps, which is consistent with a power spectrum of the form P(k)=P_0(k/k_0)^gamma where gamma = \^a H R2.6+/-m 0.2 for scales of 50-200' in the SPIRE maps. The cirrus rms fluctuation amplitude at angular scales of 100' is consistent with a modified blackbody SED with T_d = 20.1+/- 0.9 K and beta = 1.3+/- 0.2, in agreement with the values obtained above., Comment: 11 Pages; accepted for publication in MNRAS Herschel Special Issue