Search

The amount of evolution in the dust content of galaxies over the past five billion years of cosmic history is contested in the literature. Here we present a far-infrared census of dust based on a sample of 29,241 galaxies with redshifts ranging from 0 < z < 0.5 using data from the Herschel Astrophysical Terahertz Survey (H-ATLAS). We use the spectral energy distribution fitting tool MAGPHYS and a stacking analysis to investigate the evolution of dust mass and temperature of far-infrared-selected galaxies as a function of both luminosity and redshift. At low redshifts, we find that the mass-weighted and luminosity-weighted dust temperatures from the stacking analysis both exhibit a trend for brighter galaxies to have warmer dust. In higher redshift bins, we see some evolution in both mass-weighted and luminosity-weighted dust temperatures with redshift, but the effect is strongest for luminosity-weighted temperature. The measure of dust content in galaxies at z<0.1 (the Dust Mass Function) has a different shape to that derived using optically-selected galaxies from the same region of sky. We revise the local dust mass density (z<0.1) to $\rho_{\rm d} =(1.37\pm0.08)\times 10^5 {\rm\,M_{\odot}\,Mpc^{-3}}\,h_{70}^{-1}$; corresponding to an overall fraction of baryons (by mass) stored in dust of $f_{\rm mb} {(\rm dust)} = (2.22\pm 0.13) \times 10^{-5}$. We confirm evolution in both the luminosity density and dust mass density over the past few billion years ($\rho_{\rm d} \propto (1+z)^{2.6 \pm 0.6}$), with a flatter evolution than observed in previous FIR-selected studies. We attribute the evolution in $\rho_L$ and $\rho_m$ to an evolution in the dust mass., Comment: Accepted for publication in MNRAS

We use the Atacama Large sub/Millimetre Array (ALMA) to efficiently observe spectral lines across Bands 3, 4, 5, 6, 7, and 8 at high-resolution (0.5" - 0.1") for 16 bright southern Herschel sources at $1.5 < z < 4.2$. With only six and a half hours of observations, we reveal 66 spectral lines in 17 galaxies. These observations detect emission from CO (3-2) to CO(18-17), as well as atomic ([CI](1-0), (2-1), [OI] 145 $\mu$m and [NII] 205 $\mu$m) lines. Additional molecular lines are seen in emission (${\rm H_2O}$ and ${\rm H_2O^+}$) and absorption (OH$^+$ and CH$^+$). The morphologies based on dust continuum ranges from extended sources to strong lensed galaxies with magnifications between 2 and 30. CO line transitions indicate a diverse set of excitation conditions with a fraction of the sources ($\sim 35$%) showcasing dense, warm gas. The resolved gas to star-formation surface densities vary strongly per source, and suggest that the observed diversity of dusty star-forming galaxies could be a combination of lensed, compact dusty starbursts and extended, potentially-merging galaxies. The predicted gas depletion timescales are consistent with 100 Myr to 1 Gyr, but require efficient fueling from the extended gas reservoirs onto the more central starbursts, in line with the Doppler-shifted absorption lines that indicate inflowing gas for two out of six sources. This pilot paper explores a successful new method of observing spectral lines in large samples of galaxies, supports future studies of larger samples, and finds that the efficiency of this new observational method will be further improved with the planned ALMA Wideband Sensitivity Upgrade., Comment: Accepted for publication in MNRAS. 42 pages, including 20 figures + spectra. Comments and discussion are warmly welcomed

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

Variations in the dust emissivity index, $\beta$, within and between galaxies, are evidence that the chemistry and physics of dust must vary on large scales, although the nature of the physical and/or chemical variations is still unknown. In this paper we estimate values of $\beta$ and dust temperature for a sample of 109 dusty star-forming galaxies (DSFGs) over the range, $2 < z < 6$. We compare the results obtained with both an optically-thin model and a general opacity model, finding that our estimates of $\beta$ are similar between the models but our estimates of dust temperature are not. We find no evidence of a change in $\beta$ with redshift, with a median value of $\beta = 1.96$ for the optically-thin model with a confidence interval (16 - 84%) of 1.67 to 2.35 for the population. Using simulations, we estimate the measurement errors from our procedure and show that the variation of $\beta$ in the population results from intrinsic variations in the properties of the dust in DSFGs. At a fixed far-infrared luminosity, we find no evidence for a change in dust temperature, $T_\textrm{dust}$, with redshift. After allowing for the effects of correlated measurement errors, we find an inverse correlation between $\beta$ and $T_\textrm{dust}$ in DSFGs, for which there is also evidence in low-redshift galaxies., Comment: 16 pages, accepted for publication in MNRAS

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 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

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 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

Over the past decade, studies of dust in the Andromeda galaxy (M31) have shown radial variations in the dust emissivity index ($\beta$). Understanding the astrophysical reasons behind these radial variations may give clues about the chemical composition of dust grains, their physical structure, and the evolution of dust. We use $^{12}$CO(J=1-0) observations taken by the Combined Array for Research in Millimeter Astronomy (CARMA) and dust maps derived from \textit{Herschel} images, both with an angular resolution of 8" and spatial resolution of 30 pc, to study variations in $\beta$ across an area of $\approx$ 18.6 kpc$^2$ in M31. We extract sources, which we identify as molecular clouds, by applying the astrodendro algorithm to the $^{12}$CO and dust maps, which as a byproduct allows us to compare continuum emission from dust and CO emission as alternative ways of finding molecular clouds. We then use these catalogues to investigate whether there is evidence that $\beta$ is different inside and outside molecular clouds. Our results confirm the radial variations of $\beta$ seen in previous studies. However, we find little difference between the average $\beta$ inside molecular clouds compared to outside molecular clouds, in disagreement with models which predict an increase of $\beta$ in dense environments. Finally, we find some clouds traced by dust with very little CO which may be either clouds dominated by atomic gas or clouds of molecular gas that contain little CO., Comment: 14 pages, 11 figures, accepted for publication in MNRAS

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 the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenises, and creates derived data products for most of the premium multi-wavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the Herschel Atlas survey (H-ATLAS). Here, we describe the motivation and principal elements in the design of the project. Guiding principles are transparent or "open" methodologies with care for reproducibility and identification of provenance. A key element of the design focuses around the homogenisation of calibration, meta data and the provision of information required to define the selection of the data for statistical analysis. We apply probabilistic methods that extract information directly from the images at long wavelengths, exploiting the prior information available at shorter wavelengths and providing full posterior distributions rather than maximum likelihood estimates and associated uncertainties as in traditional catalogues. With this project definition paper we provide full access to the first data release of HELP; Data Release 1 (DR1), including a monolithic map of the largest SPIRE extragalactic field at 385 deg2 and 18 million measurements of PACS and SPIRE fluxes. We also provide tools to access and analyse the full HELP database. This new data set includes far-infrared photometry, photometric redshifts, and derived physical properties estimated from modelling the spectral energy distributions.

(Abridged) Exploiting the sensitivity and broad band width of NOEMA, we have studied the molecular gas and dust in the galaxy HerBS-89a, at z=2.95. High angular resolution images reveal a partial 1.0" diameter Einstein ring in the dust continuum emission and the molecular emission lines of 12CO(9-8) and H2O(2_02-1_11). We report the detection of the three fundamental transitions of the molecular ion OH+, seen in absorption; the molecular ion CH+(1-0) seen in absorption (and tentatively in emission); two transitions of amidogen (NH2), seen in emission; and HCN(11-10) and/or NH(1_2-0_1) seen in absorption. The NOEMA data are complemented with VLA data tracing the 12CO(1-0) emission line, which provides a measurement of the total mass of molecular gas and an anchor for a CO excitation analysis. In addition, we present HST imaging that reveals the foreground lensing galaxy in the near-infrared. Together with data from the GTC, we derive a photometric redshift of z(phot)~0.9 for the foreground lensing galaxy. Modelling the lensing of HerBS-89a, we reconstruct the dust continuum and molecular emission lines (magnified by a factor ~4-5) in the source plane. The 12CO(9-8) and H2O emission lines have comparable spatial and kinematic distributions; the source-plane reconstructions do not clearly distinguish between a one-component and a two-component scenario, but the latter accounts for the observed broad line widths. HerBS-89a is a powerful star forming galaxy with a dust-to-gas ratio delta(GDR)~80, a SFR = 614 +/- 59 Msun/yr and a depletion timescale tau(depl) = (3.4 +/- 1.0) 1e8 years. The OH+ and CH+ absorption lines, all have their main velocity component red-shifted by \Delta(V)~100 km/s relative to the global CO reservoir. We argue that these absorption lines trace a rare example of gas inflow towards the center of the galaxy., Comment: 33 pages, Accepted for publication on A&A

Using the EMIR instrument on the IRAM 30m telescope, we conducted a spectroscopic redshift search of seven z$_{\rm phot}$ $\sim$ 4 sub-millimetre bright galaxies selected from the Herschel Bright Sources (HerBS) sample with fluxes at 500 $\mu$m greater than 80 mJy. For four sources, we obtained spectroscopic redshifts between 3.4 < z < 4.1 through the detection of multiple CO-spectral lines with J $\leq$ 3. Later, we detected low-J transitions for two of these sources with the GBT including the CO(1-0) transition. For the remaining three sources, more data are needed to determine the spectroscopic redshift unambiguously. The measured CO luminosities and line widths suggest that all these sources are gravitationally lensed. These observations demonstrate that the 2 mm window is indispensable to confirm robust spectroscopic redshifts for z < 4 sources. Finally, we present an efficient graphical method to correctly identify spectroscopic redshifts., Comment: 20 Pages, 8 Figures + 7 Appendix Figures + 1 animated figure. Resubmitted in order to reflect the version as printed in MNRAS

We measure the 850-$\mu$m source densities of 46 candidate protoclusters selected from the Planck High-z catalogue (PHz) and the Planck Catalogue of Compact Sources (PCCS) that were followed up with Herschel-SPIRE and SCUBA-2. This paper aims to search for overdensities of 850-$\mu$m sources in order to select the fields that are most likely to be genuine protoclusters. Of the 46 candidate protoclusters, 25 have significant overdensities ($>$5 times the field counts), 11 have intermediate overdensities (3--5 times the field counts) and 10 have no overdensity ($<$3 times the field counts) of 850-$\mu$m sources. We find that the enhanced number densities are unlikely to be the result of sample variance. Compared with the number counts of another sample selected from Planck's compact source catalogues, this [PHz+PCCS]-selected sample has a higher fraction of candidate protoclusters with significant overdensities, though both samples show overdensities of 850-$\mu$m sources above intermediate level. Based on the estimated star-formation rate densities (SFRDs), we suggest that both samples can efficiently select protoclusters with starbursting galaxies near the redshift at which the global field SFRD peaks ($2 < z < 3$). Based on the confirmation of overdensities found here, future follow-up observations on other PHz targets may greatly increase the number of genuine DSFG-rich clusters/protoclusters., Comment: 8 pages, 1 figure, 1 table, accepted for publication in MNRAS

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

Dusty star-forming galaxies (DSFGs) detected at $z > 4$ provide important examples of the first generations of massive galaxies. However, few examples with spectroscopic confirmation are currently known, with Hershel struggling to detect significant numbers of $z > 6$ DSFGs. NGP6_D1 is a bright 850 $\mu m$ source (12.3 $\pm$ 2.5 mJy) with no counterparts at shorter wavelengths (a SPIRE dropout). Interferometric observations confirm it is a single source, with no evidence for any optical or NIR emission, or nearby likely foreground lensing sources. No $>3\sigma$ detected lines are seen in both LMT RSR and IRAM 30m EMIR spectra of NGP6_D1 across 32 $GHz$ of bandwidth despite reaching detection limits of $\sim 1 mJy/500 km~s^{-1}$, so the redshift remains unknown. Template fitting suggests that NGP6_D1 is most likely between $z = 5.8$ and 8.3. SED analysis finds that NGP6_D1 is a ULIRG, with a dust mass $\sim 10^8$ - $10^9$ $M_{\odot}$ and a SFR of $\sim$ 500 $M_{\odot}~yr^{-1}$. We place upper limits on the gas mass of NGP6_D1 of $M_{H2}$ $ < (1.1~\pm~3.5) \times 10^{11}$ $M_{\odot}$, consistent with a gas-to-dust ratio of $\sim$ 100 - 1000. We discuss the nature of NGP6_D1 in the context of the broader submm population, and find that comparable SPIRE dropouts account for $\sim$ 20% of all SCUBA-2 detected sources, but with a similar flux density distribution to the general population., Comment: Accepoted for publication in MNRAS

We present SCUBA-2 850-$\mu$m observations of 13 candidate starbursting protoclusters selected using Planck and Herschel data. The cumulative number counts of the 850-$\mu$m sources in 9/13 of these candidate protoclusters show significant overdensities compared to the field, with the probability $<$10$^{-2}$ assuming the sources are randomly distributed in the sky. Using the 250-, 350-, 500- and 850-$\mu$m flux densities, we estimate the photometric redshifts of individual SCUBA-2 sources by fitting spectral energy distribution (SED) templates with an MCMC method. The photometric redshift distribution, peaking at $22$ protoclusters and the peak of the cosmic star-formation rate density (SFRD). We find that the 850-$\mu$m sources in our candidate protoclusters have infrared luminosities of $L_{\mathrm{IR}}\gtrsim$10$^{12}L_{\odot}$ and star-formation rates of SFR=(500-1,500)$M_{\odot}$yr$^{-1}$. By comparing with results in the literature considering only Herschel photometry, we conclude that our 13 candidate protoclusters can be categorised into four groups: six of them being high-redshift starbursting protoclusters, one being a lower-redshift cluster/protocluster, three being protoclusters that contain lensed DSFG(s) or are rich in 850-$\mu$m sources, and three regions without significant Herschel or SCUBA-2 source overdensities. The total SFRs of the candidate protoclusters are found to be comparable or higher than those of known protoclusters, suggesting our sample contains some of the most extreme protocluster population. We infer that cross-matching Planck and Herschel data is a robust method for selecting candidate protoclusters with overdensities of 850-$\mu$m sources., Comment: 30 pages, 8 figures, 5 tables, accepted for publication in MNRAS

We explore the constraints that can be placed on the evolutionary timescales for typical low redshift galaxies evolving from the blue cloud through the green valley and onto the red sequence. We utilise galaxies from the GAMA survey with 0.1 < z < 0.2 and classify them according to the intrinsic (u-r?) colours of their stellar populations, as determined by fits to their multi-wavelength spectral energy distributions. Using these fits to also determine stellar population ages and star formation timescales, we argue that our results are consistent with a green valley population dominated by galaxies that are simply decreasing their star formation (running out of gas) over a timescale of 2-4 Gyr which are seen at a specific epoch in their evolution (approximately 1.6 e-folding times after their peak in star formation). If their fitted star formation histories are extrapolated forwards, the green galaxies will further redden over time, until they attain the colours of a passive population. In this picture, no specific quenching event which cuts-off their star formation is required, though it remains possible that the decline in star formation in green galaxies may be expedited by internal or external forces. However, there is no evidence that green galaxies have recently changed their star formation timescales relative to their previous longer term star formation histories., Comment: MNRAS, accepted

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

Unresolved gas and dust observations show a surprising diversity in the amount of interstellar matter in early-type galaxies. Using ALMA observations we resolve the ISM in z$\sim$0.05 early-type galaxies. From a large sample of early-type galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) we selected five of the dustiest cases, with dust masses M$_d\sim$several$\times10^7$M$_\odot$, with the aim of mapping their submillimetre continuum and $^{12}$CO(2-1) line emission distributions. These observations reveal molecular gas disks. There is a lack of associated, extended continuum emission in these ALMA observations, most likely because it is resolved out or surface brightness limited, if the dust distribution is as extended as the CO gas. However, two galaxies have central continuum ALMA detections. An additional, slightly offset, continuum source is revealed in one case, which may have contributed to confusion in the Herschel fluxes. Serendipitous continuum detections further away in the ALMA field are found in another case. Large and massive rotating molecular gas disks are mapped in three of our targets, reaching a few$\times10^{9}$M$_\odot$. One of these shows evidence of kinematic deviations from a pure rotating disc. The fields of our two remaining targets contain only smaller, weak CO sources, slightly offset from the optical galaxy centres. These may be companion galaxies seen in ALMA observations, or background objects. These heterogeneous findings in a small sample of dusty early-type galaxies reveal the need for more such high spatial resolution studies, to understand statistically how dust and gas are related in early-type galaxies., Comment: 15 pages, 10 figures, 3 tables, 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

We present Atacama Large Millimeter/submillimeter Array (ALMA) resolved observations of molecular gas in galaxies up to $z=0.35$ to characterise the role of global galactic dynamics on the global interstellar medium (ISM) properties. These observations consist of a sub-sample of 39 galaxies taken from the Valpara\'iso ALMA Line Emission Survey (VALES). From the CO($J=1-0)$ emission line, we quantify the kinematic parameters by modelling the velocity fields. We find that the IR luminosity increases with the rotational to dispersion velocity ratio ($V_{\rm rot}/\sigma_v$, corrected for inclination). We find a dependence between $V_{\rm rot}/\sigma_v$ and the [CII]/IR ratio, suggesting that the so-called `[CII] deficit' is related to the dynamical state of the galaxies. We find that global pressure support is needed to reconcile the dynamical mass estimates with the stellar masses in our systems with low $V_{\rm rot}/\sigma_v$ values. The star formation rate (SFR) is weakly correlated with the molecular gas fraction ($f_{\rm H_2}$) in our sample, suggesting that the release of gravitational energy from cold gas may not be the main energy source of the turbulent motions seen in the VALES galaxies. By defining a proxy of the `star formation efficiency' parameter as the SFR divided by the CO luminosity (SFE$'\equiv$ SFR/L$'_{\rm CO}$), we find a constant SFE$'$ per crossing time ($t_{\rm cross}$). We suggest that $t_{\rm cross}$ may be the controlling timescale in which the star formation occurs in dusty $z\sim0.03-0.35$ galaxies., Comment: Accepted for publication in MNRAS. 28 pages, 15 figures, 2 Tables

We present new Herschel-derived images of warm dust in the Andromeda Galaxy, M31, with unprecedented spatial resolution (~ 30 pc), column density accuracy, and constraints on the three-dimensional distributions of dust temperature and dust opacity index (hence grain size and composition), based on the new PPMAP Bayesian analysis procedure. We confirm the overall radial variation of dust opacity index reported by other recent studies, including the central decrease within ~ 3 kpc of the nucleus. We also investigate the detailed distribution of dust in the nuclear region, a prominent feature of which is a ~ 500 pc bar-like structure seen previously in H{\alpha}. The nature of this feature has been the subject of some debate. Our maps show it to be the site of the warmest dust, with a mean line-of-sight temperature ~ 30 K. A comparison with the stellar distribution, based on 2MASS data, provides strong evidence that it is a gravitationally induced bar. A comparison with radial velocity maps suggests the presence of an inflow towards the nucleus from opposite directions along this bar, fed presumably by the nuclear spiral with which it appears to connect. Such behaviour is common in large-scale bars in spiral galaxies, as is the phenomenon of nested bars whereby a subkiloparsec nuclear bar exists within a large-scale primary bar. We suggest that M31 represents an example of such nesting., Comment: 11 pages, 8 figures. 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

High-redshift, luminous, dusty star forming galaxies (DSFGs) constrain the extremity of galaxy formation theories. The most extreme are discovered through follow-up on candidates in large area surveys. Here we present 850 $\mu$m SCUBA-2 follow-up observations of 188 red DSFG candidates from the \textit{Herschel} Multi-tiered Extragalactic Survey (HerMES) Large Mode Survey, covering 274 deg$^2$. We detected 87 per cent with a signal-to-noise ratio $>$ 3 at 850~$\mu$m. We introduce a new method for incorporating the confusion noise in our spectral energy distribution fitting by sampling correlated flux density fluctuations from a confusion limited map. The new 850~$\mu$m data provide a better constraint on the photometric redshifts of the candidates, with photometric redshift errors decreasing from $\sigma_z/(1+z)\approx0.21$ to $0.15$. Comparison spectroscopic redshifts also found little bias ($\langle (z-z_{\rm spec})/(1+z_{\rm spec})\rangle = 0.08 $). The mean photometric redshift is found to be 3.6 with a dispersion of $0.4$ and we identify 21 DSFGs with a high probability of lying at $z > 4$. After simulating our selection effects we find number counts are consistent with phenomenological galaxy evolution models. There is a statistically significant excess of WISE-1 and SDSS sources near our red galaxies, giving a strong indication that lensing may explain some of the apparently extreme objects. Nevertheless, our sample should include examples of galaxies with the highest star formation rates in the Universe ($\gg10^3$ M$_\odot$yr$^{-1}$)., Comment: 22 pages, 14 figures and 3 tables

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 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

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 high-spatial-resolution ($\sim 0.12''$ or $\approx 800 \, {\rm pc}$ at $z = 4.5$) ALMA $870\,\mu$m dust continuum observations of a sample of 44 ultrared dusty star-forming galaxies (DSFGs) selected from the H-ATLAS and HerMES far-infrared surveys because of their red colors from 250 to 500 $\mu$m: $S_{500} / S_{250} > 1.5$ and $S_{500} / S_{350} > 1.0$. With photometric redshifts in the range $z \sim 4$-6, our sample includes the most luminous starbursting systems in the early Universe known so far, with total obscured star-formation rates (SFRs) of up to $\sim 4,500 \, M_\odot \, {\rm yr}^{-1}$, as well as a population of lensed, less intrinsically luminous sources. The lower limit on the number of ultrared DSFGs at 870 $\mu$m (with flux densities measured from the ALMA maps and thus not affected by source confusion) derived in this work is in reasonable agreement with models of galaxy evolution, whereas there have been reports of conflicts at 500 $\mu$m (where flux densities are derived from SPIRE). Ultrared DSFGs have a variety of morphologies (from relatively extended disks with smooth radial profiles, to compact sources, both isolated and interacting) and an average size, $\theta_{\rm FWHM}$, of $1.46 \pm 0.41\, {\rm kpc}$, considerably smaller than the values reported in previous work for less-luminous DSFGs at lower redshifts. The size and the estimated gas-depletion times of our sources are compatible with their being the progenitors of the most massive, compact, red-and-dead galaxies at $z \sim 2$-3, and ultimately of local ultra-massive elliptical galaxies or massive galaxy clusters. We are witnessing the birth of the high-mass tail of the red sequence of galaxies., Comment: Submitted to ApJ. Comments welcome

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 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 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 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 present ALMA and VLA detections of the dense molecular gas tracers HCN, HCO$^+$ and HNC in two lensed, high-redshift starbursts selected from the {\it Herschel}-ATLAS survey: {\it H}-ATLAS\,J090740.0$-$004200 (SDP.9, $z \sim 1.6$) and {\it H}-ATLAS\,J091043.1$-$000321 (SDP.11, $z \sim 1.8$). ALMA observed the $J = 3-2$ transitions in both sources, while the VLA observed the $J = 1-0$ transitions in SDP.9. We have detected all observed HCN and HCO$^+$ lines in SDP.9 and SDP.11, and also HNC(3--2) in SDP.9. The amplification factors for both galaxies have been determined from sub-arcsec resolution CO and dust emission observations carried out with NOEMA and the SMA. The HNC(1--0)/HCN(1--0) line ratio in SDP.9 suggests the presence of photon-dominated regions, as it happens to most local (U)LIRGs. The CO, HCN and HCO$^+$ SLEDs of SDP.9 are compatible to those found for many local, infrared (IR) bright galaxies, indicating that the molecular gas in local and high-redshift dusty starbursts can have similar excitation conditions. We obtain that the correlation between total IR ($L_{\rm IR}$) and dense line ($L_{\rm dense}$) luminosity in SDP.9 and SDP.11 and local star-forming galaxies can be represented by a single relation. The scatter of the $L_{\rm IR} - L_{\rm dense}$ correlation, together with the lack of sensitive dense molecular gas tracer observations for a homogeneous sample of high-redshift galaxies, prevents us from distinguishing differential trends with redshift. Our results suggest that the intense star formation found in some high-redshift dusty, luminous starbursts is associated with more massive dense molecular gas reservoirs and higher dense molecular gas fractions., Comment: Submitted to ApJ. Comments most welcome

We present the results of combined deep Keck/NIRC2, HST/WFC3 near-infrared and Herschel far infrared observations of an extremely star forming dusty lensed galaxy identified from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS J133542.9+300401). The galaxy is gravitationally lensed by a massive WISE identified galaxy cluster at $z\sim1$. The lensed galaxy is spectroscopically confirmed at $z=2.685$ from detection of $\rm {CO (1 \rightarrow 0)}$ by GBT and from detection of $\rm {CO (3 \rightarrow 2)}$ obtained with CARMA. We use the combined spectroscopic and imaging observations to construct a detailed lens model of the background dusty star-forming galaxy (DSFG) which allows us to study the source plane properties of the target. The best-fit lens model provide magnification of $\mu_{\rm star}=2.10\pm0.11$ and $\mu_{\rm dust}=2.02\pm0.06$ for the stellar and dust components respectively. Multi-band data yields a magnification corrected star formation rate of $1900(\pm200)\,M_{\odot}{\rm yr^{-1}}$ and stellar mass of $6.8_{-2.7}^{+0.9}\times10^{11}\,M_{\odot}$ consistent with a main sequence of star formation at $z\sim2.6$. The CO observations yield a molecular gas mass of $8.3(\pm1.0)\times10^{10}\,M_{\odot}$, similar to the most massive star-forming galaxies, which together with the high star-formation efficiency are responsible for the intense observed star formation rates. The lensed DSFG has a very short gas depletion time scale of $\sim40$ Myr. The high stellar mass and small gas fractions observed indicate that the lensed DSFG likely has already formed most of its stellar mass and could be a progenitor of the most massive elliptical galaxies found in the local Universe., Comment: 15 Pages, 10 Figures, 2 Tables. Accepted to the Astrophysical Journal

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

Until recently, only a handful of dusty, star-forming galaxies (DSFGs) were known at $z>4$, most of them significantly amplified by gravitational lensing. Here, we have increased the number of such DSFGs substantially, selecting galaxies from the uniquely wide 250-, 350- and 500-$\mu$m Herschel-ATLAS imaging survey on the basis of their extremely red far-infrared colors and faint 350- and 500-$\mu$m flux densities - ergo they are expected to be largely unlensed, luminous, rare and very distant. The addition of ground-based continuum photometry at longer wavelengths from the JCMT and APEX allows us to identify the dust peak in their SEDs, better constraining their redshifts. We select the SED templates best able to determine photometric redshifts using a sample of 69 high-redshift, lensed DSFGs, then perform checks to assess the impact of the CMB on our technique, and to quantify the systematic uncertainty associated with our photometric redshifts, $\sigma=0.14\,(1+z)$, using a sample of 25 galaxies with spectroscopic redshifts, each consistent with our color selection. For Herschel-selected ultrared galaxies with typical colors of $S_{500}/S_{250}\sim 2.2$ and $S_{500}/S_{350}\sim 1.3$ and flux densities, $S_{500}\sim 50\,$mJy, we determine a median redshift, $\hat{z}_{\rm phot}=3.66$, an interquartile redshift range, 3.30$-$4.27, with a median rest-frame 8$-$1000-$\mu$m luminosity, $\hat{L}_{\rm IR}$, of $1.3\times 10^{13}\,$L$_\odot$. A third lie at $z>4$, suggesting a space density, $\rho_{z>4}$, of $\approx 6 \times 10^{-7}\,$Mpc$^{-3}$. Our sample contains the most luminous known star-forming galaxies, and the most over-dense cluster of starbursting proto-ellipticals yet found., Comment: 22 pages; in press at The Astrophysical Journal; full-resolution PDF version available at http://www.roe.ac.uk/~rji/ivison2016.pdf