Search

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 $\mu$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., Comment: 16 pages, 13 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

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

We present ALMA observations of a small but statistically complete sample of twelve 250 micron selected galaxies at $z=0.35$ designed to measure their dust submillimeter continuum emission as well as their CO(1-0) and atomic carbon [CI](3P1-3P0) spectral lines. This is the first sample of galaxies with global measures of all three $H_2$-mass tracers and which show star formation rates (4-26 Msun yr$^{-1}$) and infra-red luminosities ($1-6\times10^{11}$ Lsun) typical of star forming galaxies in their era. We find a surprising diversity of morphology and kinematic structure; one-third of the sample have evidence for interaction with nearby smaller galaxies, several sources have disjoint dust and gas morphology. Moreover two galaxies have very high $L_{CI}/L_{CO}$ ratios for their global molecular gas reservoirs; if confirmed, such extreme intensity ratios in a sample of dust selected, massive star forming galaxies presents a challenge to our understanding of ISM. Finally, we use the emission of the three molecular gas tracers, to determine the carbon abundance, $X_{ci}$, and CO-$\rm{H_2}$ conversion $\alpha_{co}$ in our sample, using a weak prior that the gas-to-dust ratio is similar to that of the Milky Way for these massive and metal rich galaxies. Using a likelihood method which simultaneously uses all three gas tracer measurements, we find mean values and errors on the mean of $\alpha_{co}=3.0\pm0.5\,\rm{Msun\,(K\,kms^{-1}\,pc^2)^{-1}}$ and $X_{ci}=1.6\pm0.1\times 10^{-5}$ (or $\alpha_{ci}=18.8\,K kms^{-1}\,pc^2 (Msun)^{-1}$) and $\delta_{GDR}=128\pm16$ (or $\alpha_{850}=5.9\times10^{12}\,\rm{W\,Hz^{-1}\, Msun^{-1}}$), where our starting assumption is that these metal rich galaxies have an average gas-to-dust ratio similar to that of the Milky Way centered on $\delta_{GDR}=135$., Comment: 34 pages, published in MNRAS

We build a rigorous statistical framework to provide constraints on the chemical and dust evolution parameters for nearby late-type galaxies with a wide range of gas fractions ($3\%

We present LMT/AzTEC 1.1mm observations of $\sim100$ luminous high-redshift dusty star-forming galaxy candidates from the $\sim600\,$sq.deg $Herschel$-ATLAS survey, selected on the basis of their SPIRE red far-infrared colours and with $S_{500\mu\rm m}=35-80$ mJy. With an effective $\theta_{\rm FWHM}\approx9.5\,$ arcsec angular resolution, our observations reveal that at least 9 per cent of the targets break into multiple systems with SNR $\geq 4$ members. The fraction of multiple systems increases to $\sim23\,$ per cent (or more) if some non-detected targets are considered multiples, as suggested by the data. Combining the new AzTEC and deblended $Herschel$ photometry we derive photometric redshifts, IR luminosities, and star formation rates. While the median redshifts of the multiple and single systems are similar $(z_{\rm med}\approx3.6)$, the redshift distribution of the latter is skewed towards higher redshifts. Of the AzTEC sources $\sim85\,$ per cent lie at $z_{\rm phot}>3$ while $\sim33\,$ per cent are at $z_{\rm phot}>4$. This corresponds to a lower limit on the space density of ultra-red sources at $44$ galaxies., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 23 pages, 18 figures

We present ALMA observations of a small but statistically complete sample of twelve 250 micron selected galaxies at $z=0.35$ designed to measure their dust submillimeter continuum emission as well as their CO(1-0) and atomic carbon [CI](3P1-3P0) spectral lines. This is the first sample of galaxies with global measures of all three $H_2$-mass tracers and which show star formation rates (4-26 Msun yr$^{-1}$) and infra-red luminosities ($1-6\times10^{11}$ Lsun) typical of star forming galaxies in their era. We find a surprising diversity of morphology and kinematic structure; one-third of the sample have evidence for interaction with nearby smaller galaxies, several sources have disjoint dust and gas morphology. Moreover two galaxies have very high $L_{CI}/L_{CO}$ ratios for their global molecular gas reservoirs; if confirmed, such extreme intensity ratios in a sample of dust selected, massive star forming galaxies presents a challenge to our understanding of ISM. Finally, we use the emission of the three molecular gas tracers, to determine the carbon abundance, $X_{ci}$, and CO-$\rm{H_2}$ conversion $\alpha_{co}$ in our sample, using a weak prior that the gas-to-dust ratio is similar to that of the Milky Way for these massive and metal rich galaxies. Using a likelihood method which simultaneously uses all three gas tracer measurements, we find mean values and errors on the mean of $\alpha_{co}=3.0\pm0.5\,\rm{Msun\,(K\,kms^{-1}\,pc^2)^{-1}}$ and $X_{ci}=1.6\pm0.1\times 10^{-5}$ (or $\alpha_{ci}=18.8\,K kms^{-1}\,pc^2 (Msun)^{-1}$) and $\delta_{GDR}=128\pm16$ (or $\alpha_{850}=5.9\times10^{12}\,\rm{W\,Hz^{-1}\, Msun^{-1}}$), where our starting assumption is that these metal rich galaxies have an average gas-to-dust ratio similar to that of the Milky Way centered on $\delta_{GDR}=135$., Comment: 34 pages, published in MNRAS

We build a rigorous statistical framework to provide constraints on the chemical and dust evolution parameters for nearby late-type galaxies with a wide range of gas fractions ($3\%

We present LMT/AzTEC 1.1mm observations of $\sim100$ luminous high-redshift dusty star-forming galaxy candidates from the $\sim600\,$sq.deg $Herschel$-ATLAS survey, selected on the basis of their SPIRE red far-infrared colours and with $S_{500\mu\rm m}=35-80$ mJy. With an effective $\theta_{\rm FWHM}\approx9.5\,$ arcsec angular resolution, our observations reveal that at least 9 per cent of the targets break into multiple systems with SNR $\geq 4$ members. The fraction of multiple systems increases to $\sim23\,$ per cent (or more) if some non-detected targets are considered multiples, as suggested by the data. Combining the new AzTEC and deblended $Herschel$ photometry we derive photometric redshifts, IR luminosities, and star formation rates. While the median redshifts of the multiple and single systems are similar $(z_{\rm med}\approx3.6)$, the redshift distribution of the latter is skewed towards higher redshifts. Of the AzTEC sources $\sim85\,$ per cent lie at $z_{\rm phot}>3$ while $\sim33\,$ per cent are at $z_{\rm phot}>4$. This corresponds to a lower limit on the space density of ultra-red sources at $44$ galaxies., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 23 pages, 18 figures

Context. As recently demonstrated, high-z submillimetre galaxies (SMGs) are the perfect background sample for tracing the mass density profiles of galaxies and clusters (baryonic and dark matter) and their time-evolution through gravitational lensing. Their magnification bias, a weak gravitational lensing effect, is a powerful tool for constraining the free parameters of a halo occupation distribution (HOD) model and potentially also some of the main cosmological parameters. Aims. The aim of this work is to test the capability of the magnification bias produced on high-z SMGs as a cosmological probe. We exploit cross-correlation data to constrain not only astrophysical parameters ($M_{min}$, $M_1$, and $\alpha$), but also some of the cosmological ones ($\Omega_m$, $\sigma_8$, and $H_0$) for this proof of concept. Methods. The measured cross-correlation function between a foreground sample of GAMA galaxies with spectroscopic redshifts in the range 0.2 < z < 0.8 and a background sample of H-ATLAS galaxies with photometric redshifts >1.2 is modelled using the traditional halo model description that depends on HOD and cosmological parameters. These parameters are then estimated by performing a Markov chain Monte Carlo analysis using different sets of priors to test the robustness of the results and to study the performance of this novel observable with the current set of data Results. With our current results, $\Omega_m$ and $H_0$ cannot be well constrained. However, we can set a lower limit of >0.24 at 95\% confidence level (CL) on $\Omega_m$ and we see a slight trend towards $H_0>70$ values. For our constraints on $\sigma_8$ we obtain only a tentative peak around 0.75, but an interesting upper limit of $\sigma_8\lesssim 1$ at 95\% CL. We also study the possibility to derive better constraints by imposing more restrictive priors on the astrophysical parameters., Comment: accepted in Astronomy & Astrophysics

We describe the method used to detect sources for the Herschel-ATLAS survey. The method is to filter the individual bands using a matched filter, based on the point-spread function (PSF) and confusion noise, and then form the inverse variance weighted sum of the individual bands, including weights determined by a chosen spectral energy distribution. Peaks in this combined image are used to estimate the source positions. The fluxes for each source are estimated from the filtered single-band images, interpolated to the exact sub-pixel position. We test the method by creating simulated maps in three bands with PSFs, pixel sizes and Gaussian instrumental noise that match the 250, 350 and 500 micron bands of Herschel-ATLAS. We use our method to detect sources and compare the measured positions and fluxes to the input sources. The multi-band approach allows reliable source detection a factor 1.2 to 3 lower in flux compared to single-band source detection, depending on the source colours. The false detection rate is reduced by a factor between 4 and 10, and the variance of the source position errors is reduced by about a factor 1.5. We also consider the effect of confusion noise and find that the appropriate matched filter gives a further improvement in completeness and noise over the standard PSF filter approach. Overall the two modifications give a factor of 1.5 to 3 improvement in the depth of the recovered catalogues compared to a single-band PSF filter approach., Comment: 10 pages, 11 figures

Context. As recently demonstrated, high-z submillimetre galaxies (SMGs) are the perfect background sample for tracing the mass density profiles of galaxies and clusters (baryonic and dark matter) and their time-evolution through gravitational lensing. Their magnification bias, a weak gravitational lensing effect, is a powerful tool for constraining the free parameters of a halo occupation distribution (HOD) model and potentially also some of the main cosmological parameters. Aims. The aim of this work is to test the capability of the magnification bias produced on high-z SMGs as a cosmological probe. We exploit cross-correlation data to constrain not only astrophysical parameters ($M_{min}$, $M_1$, and $\alpha$), but also some of the cosmological ones ($\Omega_m$, $\sigma_8$, and $H_0$) for this proof of concept. Methods. The measured cross-correlation function between a foreground sample of GAMA galaxies with spectroscopic redshifts in the range 0.2 < z < 0.8 and a background sample of H-ATLAS galaxies with photometric redshifts >1.2 is modelled using the traditional halo model description that depends on HOD and cosmological parameters. These parameters are then estimated by performing a Markov chain Monte Carlo analysis using different sets of priors to test the robustness of the results and to study the performance of this novel observable with the current set of data Results. With our current results, $\Omega_m$ and $H_0$ cannot be well constrained. However, we can set a lower limit of >0.24 at 95\% confidence level (CL) on $\Omega_m$ and we see a slight trend towards $H_0>70$ values. For our constraints on $\sigma_8$ we obtain only a tentative peak around 0.75, but an interesting upper limit of $\sigma_8\lesssim 1$ at 95\% CL. We also study the possibility to derive better constraints by imposing more restrictive priors on the astrophysical parameters., Comment: accepted in Astronomy & Astrophysics

We describe the method used to detect sources for the Herschel-ATLAS survey. The method is to filter the individual bands using a matched filter, based on the point-spread function (PSF) and confusion noise, and then form the inverse variance weighted sum of the individual bands, including weights determined by a chosen spectral energy distribution. Peaks in this combined image are used to estimate the source positions. The fluxes for each source are estimated from the filtered single-band images, interpolated to the exact sub-pixel position. We test the method by creating simulated maps in three bands with PSFs, pixel sizes and Gaussian instrumental noise that match the 250, 350 and 500 micron bands of Herschel-ATLAS. We use our method to detect sources and compare the measured positions and fluxes to the input sources. The multi-band approach allows reliable source detection a factor 1.2 to 3 lower in flux compared to single-band source detection, depending on the source colours. The false detection rate is reduced by a factor between 4 and 10, and the variance of the source position errors is reduced by about a factor 1.5. We also consider the effect of confusion noise and find that the appropriate matched filter gives a further improvement in completeness and noise over the standard PSF filter approach. Overall the two modifications give a factor of 1.5 to 3 improvement in the depth of the recovered catalogues compared to a single-band PSF filter approach., Comment: 10 pages, 11 figures

The Herschel-ATLAS unbiased survey of cold dust in the local Universe is dominated by a surprising population of very blue (FUV-K < 3.5), dust-rich galaxies with high gas fractions (f_HI = M_HI/(M*+M_HI)>0.5)). Dubbed `Blue and Dusty Gas Rich Sources' (BADGRS) they have cold diffuse dust temperatures, and the highest dust-to-stellar mass ratios of any galaxies in the local Universe. Here, we explore the molecular ISM in a representative sample of BADGRS, using very deep CO(J_up=1,2,3) observations across the central and outer disk regions. We find very low CO brightnesses (Tp=15-30 mK), despite the bright far-infrared emission and metallicities in the range 0.5

The Herschel-ATLAS unbiased survey of cold dust in the local Universe is dominated by a surprising population of very blue (FUV-K < 3.5), dust-rich galaxies with high gas fractions (f_HI = M_HI/(M*+M_HI)>0.5)). Dubbed `Blue and Dusty Gas Rich Sources' (BADGRS) they have cold diffuse dust temperatures, and the highest dust-to-stellar mass ratios of any galaxies in the local Universe. Here, we explore the molecular ISM in a representative sample of BADGRS, using very deep CO(J_up=1,2,3) observations across the central and outer disk regions. We find very low CO brightnesses (Tp=15-30 mK), despite the bright far-infrared emission and metallicities in the range 0.5

We present measurements of the angular correlation function of sub-millimeter (sub-mm) galaxies (SMGs) identified in four out of the five fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) - GAMA-9h, GAMA-12h, GAMA-15h and NGP - with flux densities $S_{250\mu m}$>30 mJy at 250 {\mu}m. We show that galaxies selected at this wavelength trace the underlying matter distribution differently at low and high redshifts. We study the evolution of the clustering finding that at low redshifts sub-mm galaxies exhibit clustering strengths of $r_0$ $\sim$ 2 - 3 $h^{-1}$ Mpc, below z < 0.3. At high redshifts, on the other hand, we find that sub-mm galaxies are more strongly clustered with correlation lengths $r_0$ = 8.1 $\pm$ 0.5, 8.8 $\pm$ 0.8 and 13.9 $\pm$ 3.9 $h^{-1}$Mpc at z = 1 - 2, 2 - 3 and 3 - 5, respectively. We show that sub-mm galaxies across the redshift range 1 < z < 5, typically reside in dark-matter halos of mass of the order of ~ $10^{12.5}$ - $10^{13.0}$ $h^{-1} \, M_{\odot}$ and are consistent with being the progenitors of local massive elliptical galaxies that we see in the local Universe., Comment: 17 pages, 11 figures, 3 tables. Accepted for publication in MNRAS

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 present measurements of the angular correlation function of sub-millimeter (sub-mm) galaxies (SMGs) identified in four out of the five fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) - GAMA-9h, GAMA-12h, GAMA-15h and NGP - with flux densities $S_{250\mu m}$>30 mJy at 250 {\mu}m. We show that galaxies selected at this wavelength trace the underlying matter distribution differently at low and high redshifts. We study the evolution of the clustering finding that at low redshifts sub-mm galaxies exhibit clustering strengths of $r_0$ $\sim$ 2 - 3 $h^{-1}$ Mpc, below z < 0.3. At high redshifts, on the other hand, we find that sub-mm galaxies are more strongly clustered with correlation lengths $r_0$ = 8.1 $\pm$ 0.5, 8.8 $\pm$ 0.8 and 13.9 $\pm$ 3.9 $h^{-1}$Mpc at z = 1 - 2, 2 - 3 and 3 - 5, respectively. We show that sub-mm galaxies across the redshift range 1 < z < 5, typically reside in dark-matter halos of mass of the order of ~ $10^{12.5}$ - $10^{13.0}$ $h^{-1} \, M_{\odot}$ and are consistent with being the progenitors of local massive elliptical galaxies that we see in the local Universe., Comment: 17 pages, 11 figures, 3 tables. Accepted for publication in MNRAS

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 present an extragalactic survey using observations from the Atacama Large Millime- ter/submillimeter Array (ALMA) to characterize galaxy populations up to z = 0.35: the Valpara´ıso ALMA Line Emission Survey (VALES). We use ALMA Band-3 CO(1–0) obser- vations 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σ significance and 12 others are seen at low significance in stacked spectra. CO luminosities are in the range of (0.03–1.31) × 1010 K km s−1 pc2, equivalent to log(Mgas/M) = 8.9–10.9 assuming an αCO= 4.6 (K km s−1 pc2)−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 ∼3.5 arcsec resolution (6.5 kpc), finding that the molecular gas is on average ∼ 0.6 times more compact than the stellar component. We obtain a global Schmidt–Kennicutt relation, given by log[:ESFR/(Myr−1 kpc−2)] = (1.26 ± 0.02) × log[:EMH2/(Mpc−2)] − (3.6 ± 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 LIR = 1011–12 L. Combining our observations with data taken from the literature, we propose that star formation efficiencies can be parametrized by log [SFR/MH2 ] = 0.19 × (log LIR − 11.45) − 8.26 −0.41 × arctan[−4.84 (log LIR − 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.

We present the dust mass function (DMF) of 15,750 galaxies with redshift $z< 0.1$, drawn from the overlapping area of the GAMA and {\it H-}ATLAS surveys. The DMF is derived using the density corrected $V_{\rm max}$ method, where we estimate $V_{\rm max}$ using: (i) the normal photometric selection limit ($pV_{\rm max}$) and (ii) a bivariate brightness distribution (BBD) technique, which accounts for two selection effects. We fit the data with a Schechter function, and find $M^{*}=(4.65\pm0.18)\times 10^{7}\,h^2_{70}\, M_{\odot}$, $\alpha=(1.22\pm 0.01)$, $\phi^{*}=(6.26\pm 0.28)\times 10^{-3}\,h^3_{70}\,\rm Mpc^{-3}\,dex^{-1}$. The resulting dust mass density parameter integrated down to $10^4\,M_{\odot}$ is $\Omega_{\rm d}=(1.11 \pm0.02)\times 10^{-6}$ which implies the mass fraction of baryons in dust is $f_{m_b}=(2.40\pm0.04)\times 10^{-5}$; cosmic variance adds an extra 7-17\,per\,cent uncertainty to the quoted statistical errors. Our measurements have fewer galaxies with high dust mass than predicted by semi-analytic models. This is because the models include too much dust in high stellar mass galaxies. Conversely, our measurements find more galaxies with high dust mass than predicted by hydrodynamical cosmological simulations. This is likely to be from the long timescales for grain growth assumed in the models. We calculate DMFs split by galaxy type and find dust mass densities of $\Omega_{\rm d}=(0.88\pm0.03)\times 10^{-6}$ and $\Omega_{\rm d}=(0.060\pm0.005)\times 10^{-6}$ for late-types and early-types respectively. Comparing to the equivalent galaxy stellar mass functions (GSMF) we find that the DMF for late-types is well matched by the GMSF scaled by $(8.07\pm0.35) \times 10^{-4}$., Comment: 25 Pages, 18 Figures. Submitted December 2017

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.

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.

We present images obtained with LABOCA on the APEX telescope of a sample of 22 galaxies selected via their red Herschel SPIRE 250-, 350- and $500\textrm{-}\mu\textrm{m}$ colors. We aim to see if these luminous, rare and distant galaxies are signposting dense regions in the early Universe. Our $870\textrm{-}\mu\textrm{m}$ survey covers an area of $\approx0.8\,\textrm{deg}^2$ down to an average r.m.s. of $3.9\,\textrm{mJy beam}^{-1}$, with our five deepest maps going $\approx2\times$ deeper still. We catalog 86 DSFGs around our 'signposts', detected above a significance of $3.5\sigma$. This implies a $100\pm30\%$ over-density of $S_{870}>8.5\,\textrm{mJy}$ DSFGs, excluding our signposts, when comparing our number counts to those in 'blank fields'. Thus, we are $99.93\%$ confident that our signposts are pinpointing over-dense regions in the Universe, and $\approx95\%$ confident that these regions are over-dense by a factor of at least $\ge1.5\times$. Using template SEDs and SPIRE/LABOCA photometry we derive a median photometric redshift of $z=3.2\pm0.2$ for our signposts, with an interquartile range of $z=2.8\textrm{-}3.6$. We constrain the DSFGs likely responsible for this over-density to within $|\Delta z|\le0.65$ of their respective signposts. These 'associated' DSFGs are radially distributed within $1.6\pm0.5\,\textrm{Mpc}$ of their signposts, have median SFRs of $\approx(1.0\pm0.2)\times10^3\,M_{\odot}\,\textrm{yr}^{-1}$ (for a Salpeter stellar IMF) and median gas reservoirs of $\sim1.7\times10^{11}\,M_{\odot}$. These candidate proto-clusters have average total SFRs of at least $\approx (2.3\pm0.5)\times10^3\,M_{\odot}\,\textrm{yr}^{-1}$ and space densities of $\sim9\times10^{-7}\,\textrm{Mpc}^{-3}$, consistent with the idea that their constituents may evolve to become massive ETGs in the centers of the rich galaxy clusters we see today., Comment: 24 pages, 13 figures, 5 tables, submitted for publication in ApJ

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 report the identification of an extreme proto-cluster of galaxies in the early Universe whose core (nicknamed Distant Red Core, DRC) is formed by at least ten dusty star-forming galaxies (DSFGs), confirmed to lie at $z_{\rm spec} = 4.002$ via detection of [CI](1-0), $^{12}$CO(6-5), $^{12}$CO(4-3), $^{12}$CO(2-1) and ${\rm H_2O} (2_{11} - 2_{02})$ emission lines, detected using ALMA and ATCA. The spectroscopically-confirmed components of the proto-cluster are distributed over a ${\rm 260\, kpc \times 310\, kpc}$ region and have a collective obscured star-formation rate (SFR) of $\sim 6500 \, M_\odot \, {\rm yr}^{-1}$, considerably higher than has been seen before in any proto-cluster of galaxies or over-densities of DSFGs at $z \gtrsim 4$. Most of the star formation is taking place in luminous DSFGs since no Ly$\alpha$ emitters are detected in the proto-cluster core, apart from a Ly$\alpha$ blob located next to one of the DRC dusty components and extending over $60\,{\rm kpc}$. The total obscured SFR of the proto-cluster could rise to ${\rm SFR} \sim 14,400 \, M_\odot \, {\rm yr}^{-1}$ if all the members of an over-density of bright DSFGs discovered around DRC in a wide-field LABOCA 870-$\mu$m image are part of the same structure. The total halo mass of DRC could be as high as $\sim 4.4 \times 10^{13}\,M_\odot$ and could be the progenitor of a Coma-like cluster at $z = 0$. The relatively short gas-depletion times of the DRC components suggest either the presence of a mechanism able to trigger extreme star formation simultaneously in galaxies spread over a few hundred kpc or the presence of gas flows from the cosmic web able to sustain star formation over several hundred million years., Comment: Submitted to ApJ. Minor updates added, including a change of the source name. Comments welcome

We combine samples of nearby galaxies with Herschel photometry selected on their dust, metal, HI, and stellar mass content, and compare these to chemical evolution models in order to discriminate between different dust sources. In a companion paper, we used a HI-selected sample of nearby galaxies to reveal a sub-sample of very gas rich (gas fraction > 80 per cent) sources with dust masses significantly below predictions from simple chemical evolution models, and well below $M_d/M_*$ and $M_d/M_{gas}$ scaling relations seen in dust and stellar-selected samples of local galaxies. We use a chemical evolution model to explain these dust-poor, but gas-rich, sources as well as the observed star formation rates (SFRs) and dust-to-gas ratios. We find that (i) a delayed star formation history is required to model the observed SFRs; (ii) inflows and outflows are required to model the observed metallicities at low gas fractions; (iii) a reduced contribution of dust from supernovae (SNe) is needed to explain the dust-poor sources with high gas fractions. These dust-poor, low stellar mass galaxies require a typical core-collapse SN to produce 0.01 - 0.16 $M_{\odot}$ of dust. To match the observed dust masses at lower gas fractions, significant grain growth is required to counteract the reduced contribution from dust in SNe and dust destruction from SN shocks. These findings are statistically robust, though due to intrinsic scatter it is not always possible to find one single model that successfully describes all the data. We also show that the dust-to-metals ratio decreases towards lower metallicity., Comment: 15 pages (+10 pages appendix), 8 figures, Accepted for publication in MNRAS

The process of planet formation offers a rich source of dust production via grain growth in protostellar discs, and via grinding of larger bodies in debris disc systems. Chemical evolution models, designed to follow the build up of metals and dust in galaxies, do not currently account for planet formation. We consider the possibility that the apparent under-prediction of dust mass in high redshift galaxies by chemical evolution models could be in part, due to these models neglecting this process, specifically due to their assumption that a large fraction of the dust mass is removed from the interstellar medium during star formation (so-called astration). By adding a planet formation phase into galaxy chemical evolution, we demonstrate that the dust budget crisis can be partially ameliorated by a factor of 1.3-1.5 only if a) circumstellar discs prevent a large fraction of the dust mass entering the star during its birth, and b) that dust mass is preferentially liberated via jets, winds and outflows rather than accreted into planetary-mass bodies., Comment: 9 pages, 11 figures, accepted for publication in MNRAS

We present 1.3- and/or 3-mm continuum images and 3-mm spectral scans, obtained using NOEMA and ALMA, of 21 distant, dusty, star-forming galaxies (DSFGs). Our sample is a subset of the galaxies selected by Ivison et al. (2016) on the basis of their extremely red far-infrared (far-IR) colours and low {\it Herschel} flux densities; most are thus expected to be unlensed, extraordinarily luminous starbursts at $z \gtrsim 4$, modulo the considerable cross-section to gravitational lensing implied by their redshift. We observed 17 of these galaxies with NOEMA and four with ALMA, scanning through the 3-mm atmospheric window. We have obtained secure redshifts for seven galaxies via detection of multiple CO lines, one of them a lensed system at $z=6.027$ (two others are also found to be lensed); a single emission line was detected in another four galaxies, one of which has been shown elsewhere to lie at $z=4.002$. Where we find no spectroscopic redshifts, the galaxies are generally less luminous by 0.3-0.4 dex, which goes some way to explaining our failure to detect line emission. We show that this sample contains amongst the most luminous known star-forming galaxies. Due to their extreme star-formation activity, these galaxies will consume their molecular gas in $\lesssim 100$ Myr, despite their high molecular gas masses, and are therefore plausible progenitors of the massive, `red-and-dead' elliptical galaxies at $z \approx 3$., Comment: 15 pages, 10 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 derive the low redshift galaxy stellar mass function (GSMF), inclusive of dust corrections, for the equatorial Galaxy And Mass Assembly (GAMA) dataset covering 180 deg$^2$. We construct the mass function using a density-corrected maximum volume method, using masses corrected for the impact of optically thick and thin dust. We explore the galactic bivariate brightness plane ($M_\star-\mu$), demonstrating that surface brightness effects do not systematically bias our mass function measurement above 10$^{7.5}$ M$_{\odot}$. The galaxy distribution in the $M-\mu$-plane appears well bounded, indicating that no substantial population of massive but diffuse or highly compact galaxies are systematically missed due to the GAMA selection criteria. The GSMF is {fit with} a double Schechter function, with $\mathcal M^\star=10^{10.78\pm0.01\pm0.20}M_\odot$, $\phi^\star_1=(2.93\pm0.40)\times10^{-3}h_{70}^3$Mpc$^{-3}$, $\alpha_1=-0.62\pm0.03\pm0.15$, $\phi^\star_2=(0.63\pm0.10)\times10^{-3}h_{70}^3$Mpc$^{-3}$, and $\alpha_2=-1.50\pm0.01\pm0.15$. We find the equivalent faint end slope as previously estimated using the GAMA-I sample, although we find a higher value of $\mathcal M^\star$. Using the full GAMA-II sample, we are able to fit the mass function to masses as low as $10^{7.5}$ $M_\odot$, and assess limits to $10^{6.5}$ $M_\odot$. Combining GAMA-II with data from G10-COSMOS we are able to comment qualitatively on the shape of the GSMF down to masses as low as $10^{6}$ $M_\odot$. Beyond the well known upturn seen in the GSMF at $10^{9.5}$ the distribution appears to maintain a single power-law slope from $10^9$ to $10^{6.5}$. We calculate the stellar mass density parameter given our best-estimate GSMF, finding $\Omega_\star= 1.66^{+0.24}_{-0.23}\pm0.97 h^{-1}_{70} \times 10^{-3}$, inclusive of random and systematic uncertainties., Comment: 22 pages; 15 figures; Accepted for publication in MNRAS

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

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.

We present the dust mass function (DMF) of 15,750 galaxies with redshift $z< 0.1$, drawn from the overlapping area of the GAMA and {\it H-}ATLAS surveys. The DMF is derived using the density corrected $V_{\rm max}$ method, where we estimate $V_{\rm max}$ using: (i) the normal photometric selection limit ($pV_{\rm max}$) and (ii) a bivariate brightness distribution (BBD) technique, which accounts for two selection effects. We fit the data with a Schechter function, and find $M^{*}=(4.65\pm0.18)\times 10^{7}\,h^2_{70}\, M_{\odot}$, $\alpha=(1.22\pm 0.01)$, $\phi^{*}=(6.26\pm 0.28)\times 10^{-3}\,h^3_{70}\,\rm Mpc^{-3}\,dex^{-1}$. The resulting dust mass density parameter integrated down to $10^4\,M_{\odot}$ is $\Omega_{\rm d}=(1.11 \pm0.02)\times 10^{-6}$ which implies the mass fraction of baryons in dust is $f_{m_b}=(2.40\pm0.04)\times 10^{-5}$; cosmic variance adds an extra 7-17\,per\,cent uncertainty to the quoted statistical errors. Our measurements have fewer galaxies with high dust mass than predicted by semi-analytic models. This is because the models include too much dust in high stellar mass galaxies. Conversely, our measurements find more galaxies with high dust mass than predicted by hydrodynamical cosmological simulations. This is likely to be from the long timescales for grain growth assumed in the models. We calculate DMFs split by galaxy type and find dust mass densities of $\Omega_{\rm d}=(0.88\pm0.03)\times 10^{-6}$ and $\Omega_{\rm d}=(0.060\pm0.005)\times 10^{-6}$ for late-types and early-types respectively. Comparing to the equivalent galaxy stellar mass functions (GSMF) we find that the DMF for late-types is well matched by the GMSF scaled by $(8.07\pm0.35) \times 10^{-4}$., Comment: 25 Pages, 18 Figures. Submitted December 2017

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 images obtained with LABOCA on the APEX telescope of a sample of 22 galaxies selected via their red Herschel SPIRE 250-, 350- and $500\textrm{-}\mu\textrm{m}$ colors. We aim to see if these luminous, rare and distant galaxies are signposting dense regions in the early Universe. Our $870\textrm{-}\mu\textrm{m}$ survey covers an area of $\approx0.8\,\textrm{deg}^2$ down to an average r.m.s. of $3.9\,\textrm{mJy beam}^{-1}$, with our five deepest maps going $\approx2\times$ deeper still. We catalog 86 DSFGs around our 'signposts', detected above a significance of $3.5\sigma$. This implies a $100\pm30\%$ over-density of $S_{870}>8.5\,\textrm{mJy}$ DSFGs, excluding our signposts, when comparing our number counts to those in 'blank fields'. Thus, we are $99.93\%$ confident that our signposts are pinpointing over-dense regions in the Universe, and $\approx95\%$ confident that these regions are over-dense by a factor of at least $\ge1.5\times$. Using template SEDs and SPIRE/LABOCA photometry we derive a median photometric redshift of $z=3.2\pm0.2$ for our signposts, with an interquartile range of $z=2.8\textrm{-}3.6$. We constrain the DSFGs likely responsible for this over-density to within $|\Delta z|\le0.65$ of their respective signposts. These 'associated' DSFGs are radially distributed within $1.6\pm0.5\,\textrm{Mpc}$ of their signposts, have median SFRs of $\approx(1.0\pm0.2)\times10^3\,M_{\odot}\,\textrm{yr}^{-1}$ (for a Salpeter stellar IMF) and median gas reservoirs of $\sim1.7\times10^{11}\,M_{\odot}$. These candidate proto-clusters have average total SFRs of at least $\approx (2.3\pm0.5)\times10^3\,M_{\odot}\,\textrm{yr}^{-1}$ and space densities of $\sim9\times10^{-7}\,\textrm{Mpc}^{-3}$, consistent with the idea that their constituents may evolve to become massive ETGs in the centers of the rich galaxy clusters we see today., Comment: 24 pages, 13 figures, 5 tables, submitted for publication in ApJ

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.

We report the identification of an extreme proto-cluster of galaxies in the early Universe whose core (nicknamed Distant Red Core, DRC) is formed by at least ten dusty star-forming galaxies (DSFGs), confirmed to lie at $z_{\rm spec} = 4.002$ via detection of [CI](1-0), $^{12}$CO(6-5), $^{12}$CO(4-3), $^{12}$CO(2-1) and ${\rm H_2O} (2_{11} - 2_{02})$ emission lines, detected using ALMA and ATCA. The spectroscopically-confirmed components of the proto-cluster are distributed over a ${\rm 260\, kpc \times 310\, kpc}$ region and have a collective obscured star-formation rate (SFR) of $\sim 6500 \, M_\odot \, {\rm yr}^{-1}$, considerably higher than has been seen before in any proto-cluster of galaxies or over-densities of DSFGs at $z \gtrsim 4$. Most of the star formation is taking place in luminous DSFGs since no Ly$\alpha$ emitters are detected in the proto-cluster core, apart from a Ly$\alpha$ blob located next to one of the DRC dusty components and extending over $60\,{\rm kpc}$. The total obscured SFR of the proto-cluster could rise to ${\rm SFR} \sim 14,400 \, M_\odot \, {\rm yr}^{-1}$ if all the members of an over-density of bright DSFGs discovered around DRC in a wide-field LABOCA 870-$\mu$m image are part of the same structure. The total halo mass of DRC could be as high as $\sim 4.4 \times 10^{13}\,M_\odot$ and could be the progenitor of a Coma-like cluster at $z = 0$. The relatively short gas-depletion times of the DRC components suggest either the presence of a mechanism able to trigger extreme star formation simultaneously in galaxies spread over a few hundred kpc or the presence of gas flows from the cosmic web able to sustain star formation over several hundred million years., Comment: Submitted to ApJ. Minor updates added, including a change of the source name. Comments welcome

The process of planet formation offers a rich source of dust production via grain growth in protostellar discs, and via grinding of larger bodies in debris disc systems. Chemical evolution models, designed to follow the build up of metals and dust in galaxies, do not currently account for planet formation. We consider the possibility that the apparent under-prediction of dust mass in high redshift galaxies by chemical evolution models could be in part, due to these models neglecting this process, specifically due to their assumption that a large fraction of the dust mass is removed from the interstellar medium during star formation (so-called astration). By adding a planet formation phase into galaxy chemical evolution, we demonstrate that the dust budget crisis can be partially ameliorated by a factor of 1.3-1.5 only if a) circumstellar discs prevent a large fraction of the dust mass entering the star during its birth, and b) that dust mass is preferentially liberated via jets, winds and outflows rather than accreted into planetary-mass bodies., Comment: 9 pages, 11 figures, accepted for publication in MNRAS

We present 1.3- and/or 3-mm continuum images and 3-mm spectral scans, obtained using NOEMA and ALMA, of 21 distant, dusty, star-forming galaxies (DSFGs). Our sample is a subset of the galaxies selected by Ivison et al. (2016) on the basis of their extremely red far-infrared (far-IR) colours and low {\it Herschel} flux densities; most are thus expected to be unlensed, extraordinarily luminous starbursts at $z \gtrsim 4$, modulo the considerable cross-section to gravitational lensing implied by their redshift. We observed 17 of these galaxies with NOEMA and four with ALMA, scanning through the 3-mm atmospheric window. We have obtained secure redshifts for seven galaxies via detection of multiple CO lines, one of them a lensed system at $z=6.027$ (two others are also found to be lensed); a single emission line was detected in another four galaxies, one of which has been shown elsewhere to lie at $z=4.002$. Where we find no spectroscopic redshifts, the galaxies are generally less luminous by 0.3-0.4 dex, which goes some way to explaining our failure to detect line emission. We show that this sample contains amongst the most luminous known star-forming galaxies. Due to their extreme star-formation activity, these galaxies will consume their molecular gas in $\lesssim 100$ Myr, despite their high molecular gas masses, and are therefore plausible progenitors of the massive, `red-and-dead' elliptical galaxies at $z \approx 3$., Comment: 15 pages, 10 figures; accepted for publication in MNRAS

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 derive the low redshift galaxy stellar mass function (GSMF), inclusive of dust corrections, for the equatorial Galaxy And Mass Assembly (GAMA) dataset covering 180 deg$^2$. We construct the mass function using a density-corrected maximum volume method, using masses corrected for the impact of optically thick and thin dust. We explore the galactic bivariate brightness plane ($M_\star-\mu$), demonstrating that surface brightness effects do not systematically bias our mass function measurement above 10$^{7.5}$ M$_{\odot}$. The galaxy distribution in the $M-\mu$-plane appears well bounded, indicating that no substantial population of massive but diffuse or highly compact galaxies are systematically missed due to the GAMA selection criteria. The GSMF is {fit with} a double Schechter function, with $\mathcal M^\star=10^{10.78\pm0.01\pm0.20}M_\odot$, $\phi^\star_1=(2.93\pm0.40)\times10^{-3}h_{70}^3$Mpc$^{-3}$, $\alpha_1=-0.62\pm0.03\pm0.15$, $\phi^\star_2=(0.63\pm0.10)\times10^{-3}h_{70}^3$Mpc$^{-3}$, and $\alpha_2=-1.50\pm0.01\pm0.15$. We find the equivalent faint end slope as previously estimated using the GAMA-I sample, although we find a higher value of $\mathcal M^\star$. Using the full GAMA-II sample, we are able to fit the mass function to masses as low as $10^{7.5}$ $M_\odot$, and assess limits to $10^{6.5}$ $M_\odot$. Combining GAMA-II with data from G10-COSMOS we are able to comment qualitatively on the shape of the GSMF down to masses as low as $10^{6}$ $M_\odot$. Beyond the well known upturn seen in the GSMF at $10^{9.5}$ the distribution appears to maintain a single power-law slope from $10^9$ to $10^{6.5}$. We calculate the stellar mass density parameter given our best-estimate GSMF, finding $\Omega_\star= 1.66^{+0.24}_{-0.23}\pm0.97 h^{-1}_{70} \times 10^{-3}$, inclusive of random and systematic uncertainties., Comment: 22 pages; 15 figures; Accepted for publication in MNRAS

We combine samples of nearby galaxies with Herschel photometry selected on their dust, metal, HI, and stellar mass content, and compare these to chemical evolution models in order to discriminate between different dust sources. In a companion paper, we used a HI-selected sample of nearby galaxies to reveal a sub-sample of very gas rich (gas fraction > 80 per cent) sources with dust masses significantly below predictions from simple chemical evolution models, and well below $M_d/M_*$ and $M_d/M_{gas}$ scaling relations seen in dust and stellar-selected samples of local galaxies. We use a chemical evolution model to explain these dust-poor, but gas-rich, sources as well as the observed star formation rates (SFRs) and dust-to-gas ratios. We find that (i) a delayed star formation history is required to model the observed SFRs; (ii) inflows and outflows are required to model the observed metallicities at low gas fractions; (iii) a reduced contribution of dust from supernovae (SNe) is needed to explain the dust-poor sources with high gas fractions. These dust-poor, low stellar mass galaxies require a typical core-collapse SN to produce 0.01 - 0.16 $M_{\odot}$ of dust. To match the observed dust masses at lower gas fractions, significant grain growth is required to counteract the reduced contribution from dust in SNe and dust destruction from SN shocks. These findings are statistically robust, though due to intrinsic scatter it is not always possible to find one single model that successfully describes all the data. We also show that the dust-to-metals ratio decreases towards lower metallicity., Comment: 15 pages (+10 pages appendix), 8 figures, Accepted for publication in MNRAS

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