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.

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.

We investigate the evolution of the gas mass fraction for galaxies in the COSMOS field using submillimetre emission from dust at 850$\mu$m. We use stacking methodologies on the 850$\mu$m S2COSMOS map to derive the gas mass fraction of galaxies out to high redshifts, 0 <= $z$ <= 5, for galaxies with stellar masses of $10^{9.5} < M_* (\rm M_{\odot}) < 10^{11.75}$. In comparison to previous literature studies we extend to higher redshifts, include more normal star-forming galaxies (on the main sequence), and also investigate the evolution of the gas mass fraction split by star-forming and passive galaxy populations. We find our stacking results broadly agree with scaling relations in the literature. We find tentative evidence for a peak in the gas mass fraction of galaxies at around $z$ ~ 2.5-3, just before the peak of the star formation history of the Universe. We find that passive galaxies are particularly devoid of gas, compared to the star-forming population. We find that even at high redshifts, high stellar mass galaxies still contain significant amounts of gas.

Verifying that sub-mm galaxies (SMGs) are gravitationally lensed requires time-expensive observations with over-subscribed high-resolution observatories. Here, we aim to strengthen the evidence of gravitational lensing within the Herschel Bright Sources (HerBS) by cross-comparing their positions to optical (SDSS) and near-infrared (VIKING) surveys, in order to search for the foreground lensing galaxy candidates. Resolved observations of the brightest HerBS sources have already shown that most are lensed, and a galaxy evolution model predicts that $\sim$76% of the total HerBS sources are lensed, although with the SDSS survey we are only able to identify the likely foreground lenses for 25% of the sources. With the near-infrared VIKING survey, however, we are able to identify the likely foreground lenses for 57% of the sources, and we estimate that 82% of the HerBS sources have lenses on the VIKING images even if we cannot identify the lens in every case. We find that the angular offsets between lens and Herschel source are larger than that expected if the lensing is done by individual galaxies. We also find that the fraction of HerBS sources that are lensed falls with decreasing 500-micron flux density, which is expected from the galaxy evolution model. Finally, we apply our statistical VIKING cross-identification to the entire Herschel-ATLAS catalogue, where we also find that the number of lensed sources falls with decreasing 500-micron flux density.

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<z<4.05$ with a median value of $z=2.9\pm$0.6. The sources are unresolved or barely resolved on scales of 10 kpc. In one field, two galaxies with different redshifts were detected. In two cases the sources are found to be binary galaxies with projected distances of ~140 kpc. The linewidths of the sources are large, with a mean value for the full width at half maximum of 700$\pm$300 km/s and a median of 800 km/s. We analyse the nature of the sources with currently available ancillary data to determine if they are lensed or hyper-luminous $L_{\rm FIR} > 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.

We present a CO(3-2) survey of selected regions in the M31 disc as part of the JCMT large programme, HARP and SCUBA-2 High-Resolution Terahertz Andromeda Galaxy Survey (HASHTAG). The 12 CO(3-2) fields in this survey cover a total area of 60 square arcminutes, spanning a deprojected radial range of 2 - 14 kpc across the M31 disc. Combining these observations with existing IRAM 30m CO(1-0) observations and JCMT CO(3-2) maps of the nuclear region of M31, as well as dust temperature and star formation rate surface density maps, we are able to explore the radial distribution of the CO(3-2)/CO(1-0) integrated intensity ratio (R31) and its relationship with dust temperature and star formation. We find that the value of R31 between 2 - 9 kpc galactocentric radius is 0.14, significantly lower than what is seen in the nuclear ring at ~1 kpc (R31 ~ 0.8), only to rise again to 0.27 for the fields centred on the 10 kpc star forming ring. We also found that R31 is positively correlated with dust temperature, with Spearman's rank correlation coefficient $\rho$ = 0.55. The correlation between star formation rate surface density and CO(3--2) intensity is much stronger than with CO(1-0), with $\rho$ = 0.54 compared to -0.05, suggesting that the CO(3-2) line traces warmer and denser star forming gas better. We also find that R31 correlates well with star formation rate surface density, with $\rho$ = 0.69.

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

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.

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 $2<z<3$, is consistent with that of known $z>2$ 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.

In this paper, we use the Herschel Reference Survey to make a direct test of the hypothesis that the growth of a stellar bulge leads to a reduction in the star-formation efficiency of a galaxy (or conversely a growth in the gas-depletion timescale) as a result of the stabilisation of the gaseous disk by the gravitational field of the bulge. We find a strong correlation between star-formation efficiency and specific star-formation rate in galaxies without prominent bulges and in galaxies of the same morphological type, showing that there must be some other process besides the growth of a bulge that reduces the star-formation efficiency in galaxies. However, we also find that galaxies with more prominent bulges (Hubble types E to Sab) do have significantly lower star-formation efficiencies than galaxies with later morphological types, which is at least consistent with the hypothesis that the growth of a bulge leads to the reduction in the star-formation efficiency. The answer to the question in the title is therefore, yes and no: bulges may reduce the star-formation efficiency in galaxies but there must also be some other process at work. We also find that there is a significant but small difference in the star-formation efficiencies of galaxies with and without bars, in the sense that galaxies with bars have slightly higher star-formation efficiencies.

We study the dust properties of 192 nearby galaxies from the JINGLE survey using photometric data in the 22-850micron range. We derive the total dust mass, temperature T and emissivity index beta of the galaxies through the fitting of their spectral energy distribution (SED) using a single modified black-body model (SMBB). We apply a hierarchical Bayesian approach that reduces the known degeneracy between T and beta. Applying the hierarchical approach, the strength of the T-beta anti-correlation is reduced from a Pearson correlation coefficient R=-0.79 to R=-0.52. For the JINGLE galaxies we measure dust temperatures in the range 17-30 K and dust emissivity indices beta in the range 0.6-2.2. We compare the SMBB model with the broken emissivity modified black-body (BMBB) and the two modified black-bodies (TMBB) models. The results derived with the SMBB and TMBB are in good agreement, thus applying the SMBB, which comes with fewer free parameters, does not penalize the measurement of the cold dust properties in the JINGLE sample. We investigate the relation between T and beta and other global galaxy properties in the JINGLE and Herschel Reference Survey (HRS) sample. We find that beta correlates with the stellar mass surface density (R=0.62) and anti-correlates with the HI mass fraction (M(HI)/M*, R=-0.65), whereas the dust temperature correlates strongly with the SFR normalized by the dust mass (R=0.73). These relations can be used to estimate T and beta in galaxies with insufficient photometric data available to measure them directly through SED fitting.

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

We have analysed Herschel observations of M31, using the PPMAP procedure. The resolution of PPMAP images is sufficient (31 pc on M31) that we can analyse far-IR dust emission on the scale of Giant Molecular Clouds. By comparing PPMAP estimates of the far-IR emission optical depth at 300 microns (tau_300), and the near-IR extinction optical depth at 1.1 microns (tau_1.1) obtained from the reddening of RGB stars, we show that the ratio R_OBS.tau = tau_1.1/tau_300 falls in the range 500 to 1500. Such low values are incompatible with many commonly used theoretical dust models, which predict values of R_MODEL.kappa = kappa_1.1/kappa_300 (where kappa is the dust opacity coefficient) in the range 2500 to 4000. That is, unless a large fraction, at least 60%, of the dust emitting at 300 microns is in such compact sources that they are unlikely to intercept the lines of sight to a distributed population like RGB stars. This is not a new result: variants obtained using different observations and/or different wavelengths have already been reported by other studies. We present two analytic arguments for why it is unlikely that at least 60% of the emitting dust is in sufficiently compact sources. Therefore it may be necessary to explore the possibility that the discrepancy between observed values of R_OBS.tau and theoretical values of R_MODEL.kappa is due to limitations in existing dust models. PPMAP also allows us to derive optical-depth weighted mean values for the emissivity index, beta = - dln(kappa_lambda)/dln(lambda), and the dust temperature, T, denoted betabar and Tbar. We show that, in M31, R_OBS.tau is anti-correlated with betabar according to R_OBS.tau = 2042(+/-24)-557(+/-10)betabar. If confirmed, this provides a challenging constraint on the nature of interstellar dust in M31.

We present SCUBA-2 450$\mu$m and 850$\mu$m data of the mature redshift 2 cluster CLJ1449. We combine this with archival Herschel data to explore the star forming properties of CLJ1449. Using high resolution ALMA and JVLA data we identify potentially confused galaxies, and use the Bayesian inference tool XID+ to estimate fluxes for them. Using archival optical and near infrared data with the energy-balance code CIGALE we calculate star formation rates, and stellar masses for all our cluster members, and find the star formation rate varies between 20-1600M$_{\odot}$yr$^{-1}$ over the entire 3Mpc radial range. The central 0.5Mpc region itself has a total star formation rate of 800$\pm$200M$_{\odot}$yr$^{-1}$, which corresponds to a star formation rate density of (1.2$\pm$0.3)$\times$10$^{4}$M$_{\odot}$yr$^{-1}$Mpc$^{-3}$, which is approximately five orders of magnitude greater than expected field values. When comparing this cluster to those at lower redshifts we find that there is an increase in star formation rate per unit volume towards the centre of the cluster. This indicates that there is indeed a reversal in the star formation/density relation in CLJ1449. Based on the radial star-formation rate density profile, we see evidence for an elevation in the star formation rate density, even out to radii of 3Mpc. At these radii the elevation could be an order of magnitude greater than field values, but the exact number cannot be determined due to ambiguity in the redshift associations. If this is the case it would imply that this cluster is still accreting material which is possibly interacting and undergoing vigorous star-formation.

Follow-up observations of (sub-)mm-selected gravitationally-lensed systems have allowed a more detailed study of the dust-enshrouded phase of star-formation up to very early cosmic times. Here, the case of the gravitationally lensed merger in HATLAS J142935.3-002836 (also known as H1429-0028; z_lens=0.218, z_bkg=1.027) is revisited following recent developments in the literature and new APEX observations targeting two carbon monoxide (CO) rotational transitions J_up=3 and 6. We show that the line-profiles comprise three distinct velocity components, where the fainter high-velocity one is less magnified and more compact. The modelling of the observed spectral line energy distribution of CO J_up=2 to 6 and [CI]3P_1-3P_0 assumes a large velocity gradient scenario, where the analysis is based on four statistical approaches. Since the detected gas and dust emission comes exclusively from only one of the two merging components (the one oriented North-South, NS), we are only able to determine upper-limits for the companion. The molecular gas in the NS component in H1429-0028 is found to have a temperature of ~70K, a volume density of log(n[/cm3])~3.7, to be expanding at ~10km/s/pc, and amounts to M_H2=4(-2,+3)*1e9 Msun. The CO to H2 conversion factor is estimated to be alpha_CO=0.4(-0.2,+0.3) Msun/(K.km/s.pc2). The NS galaxy is expected to have a factor of >10x more gas than its companion (M_H2<3e8 Msun). Nevertheless, the total amount of molecular gas in the system comprises only up to 15 per cent (1sigma upper-limit) of the total (dynamical) mass.

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.

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)

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.

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.

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.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

The galaxies found in optical surveys fall in two distinct regions of a diagram of optical colour versus absolute magnitude: the red sequence and the blue cloud with the green valley in between. We show that the galaxies found in a submillimetre survey have almost the opposite distribution in this diagram, forming a `green mountain'. We show that these distinctive distributions follow naturally from a single, continuous, curved Galaxy Sequence in a diagram of specific star-formation rate versus stellar mass without there being the need for a separate star-forming galaxy Main Sequence and region of passive galaxies. The cause of the red sequence and the blue cloud is the geometric mapping between stellar mass/specific star-formation rate and absolute magnitude/colour, which distorts a continuous Galaxy Sequence in the diagram of intrinsic properties into a bimodal distribution in the diagram of observed properties. The cause of the green mountain is Malmquist bias in the submillimetre waveband, with submillimetre surveys tending to select galaxies on the curve of the Galaxy Sequence, which have the highest ratios of submillimetre-to-optical luminosity. This effect, working in reverse, causes galaxies on the curve of the Galaxy Sequence to be underrepresented in optical samples, deepening the green valley. The green valley is therefore not evidence (1) for there being two distinct populations of galaxies, (2) for galaxies in this region evolving more quickly than galaxies in the blue cloud and the red sequence, (c) for rapid quenching processes in the galaxy population.

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.

The most intensively star-forming galaxies are extremely luminous at far-infrared (FIR) wavelengths, highly obscured at optical and ultraviolet wavelengths, and lie at $z\ge 1-3$. We present a programme of ${\it Herschel}$ FIR spectroscopic observations with the SPIRE FTS and photometric observations with PACS, both on board ${\it Herschel}$, towards a sample of 45 gravitationally lensed, dusty starbursts across $z\sim 1-3.6$. In total, we detected 27 individual lines down to 3-$\sigma$, including nine $[\rm C{\small II}]$ 158-$\mu$m lines with confirmed spectroscopic redshifts, five possible $[\rm C{\small II}]$ lines consistent with their far-infrared photometric redshifts, and in some individual sources a few $[\rm O{\small III}]$ 88-$\mu$m, $[\rm O{\small III}]$ 52-$\mu$m, $[\rm O{\small I}]$ 145-$\mu$m, $[\rm O{\small I}]$ 63-$\mu$m, $[\rm N{\small II}]$ 122-$\mu$m, and OH 119-$\mu$m (in absorption) lines. To derive the typical physical properties of the gas in the sample, we stack all spectra weighted by their intrinsic luminosity and by their 500-$\mu$m flux densities, with the spectra scaled to a common redshift. In the stacked spectra, we detect emission lines of $[\rm C{\small II}]$ 158-$\mu$m, $[\rm N{\small II}]$ 122-$\mu$m, $[\rm O{\small III}]$ 88-$\mu$m, $[\rm O{\small III}]$ 52-$\mu$m, $[\rm O{\small I}]$ 63-$\mu$m, and the absorption doublet of OH at 119-$\mu$m, at high fidelity. We find that the average electron densities traced by the $[\rm N{\small II}]$ and $[\rm O{\small III}]$ lines are higher than the average values in local star-forming galaxies and ULIRGs, using the same tracers. From the $[\rm N{\small II}]/[\rm C{\small II}]$ and $[\rm O{\small I}]/[\rm C{\small II}]$ ratios, we find that the $[\rm C{\small II}]$ emission is likely dominated by the photo-dominated regions (PDR), instead of by ionised gas or large-scale shocks.

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}$.

Since their discovery, submillimetre-selected galaxies (SMGs) have revolutionized the field of galaxy formation and evolution. From the hundreds of square degrees mapped at submillimetre wavelengths, only a handful of sources have been confirmed to lie at z>5 and only two at z>6. All of these SMGs are rare examples of extreme starburst galaxies with star formation rates (SFRs) of >1000 M_sun/yr and therefore are not representative of the general population of dusty star-forming galaxies. Consequently, our understanding of the nature of these sources, at the earliest epochs, is still incomplete. Here we report the spectroscopic identification of a gravitationally amplified (mu = 9.3 +/- 1.0) dusty star-forming galaxy at z=6.027. After correcting for gravitational lensing, we derive an intrinsic less-extreme SFR of 380 +/- 50 M_sun/yr for this source and find that its gas and dust properties are similar to those measured for local Ultra Luminous Infrared Galaxies (ULIRGs), extending the local trends to a poorly explored territory in the early Universe. The star-formation efficiency of this galaxy is similar to those measured in its local analogues, despite a ~12 Gyr difference in cosmic time.

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

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}$).

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.

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.

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.

We present the largest submillimeter images that have been made of the extragalactic sky. The 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{\mu}m. In this paper we present the images from our two largest fields which account for ~75% of the survey. The first field is 180.1 deg$^2$ in size centered on the North Galactic Pole (NGP) and the second field is 317.6 deg$^2$ in size centered on the South Galactic Pole. The NGP field serendipitously contains the Coma cluster. Over most (~80%) of the images, the pixel noise, including both instrumental noise and confusion noise, is approximately 3.6, and 3.5 mJy/pix at 100 and 160{\mu}m, and 11.0, 11.1 and 12.3 mJy/beam at 250, 350 and 500{\mu}m, respectively, but reaches lower values in some parts of the images. If a matched filter is applied to optimize point-source detection, our total 1{\sigma} map sensitivity is 5.7, 6.0, and 7.3 mJy at 250, 350, and 500{\mu}m, respectively. We describe the results of an investigation of the noise properties of the images. We make the most precise estimate of confusion in SPIRE maps to date finding values of 3.12+/-0.07, 4.13+/-0.02 and 4.45+/-0.04 mJy/beam at 250, 350, and 500{\mu}m in our un-convolved maps. For PACS we find an estimate of the confusion noise in our fast-parallel observations of 4.23 and 4.62 mJy/beam at 100 and 160{\mu}m. Finally, we give recipes for using these images to carry out photometry, both for unresolved and extended sources.

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.

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<z<0.8$, and a background sample of H-ATLAS galaxies with photometric redshifts $\gtrsim1.2$. It constitutes a substantial improvement over the cross-correlation measurements made by Gonzalez-Nuevo et al. (2014) with updated catalogues and wider area (with $S/N\gtrsim 5$ below 10' and reaching $S/N\sim 20$ below 30"). The better statistics allow us to split the sample in different redshift bins and to perform a tomographic analysis (with $S/N\gtrsim 3$ below 10 arcmin and reaching $S/N\sim 15$ below 30"). Moreover, we implement a halo model to extract astrophysical information about the background galaxies and the deflectors that are producing the lensing link between the foreground (lenses) and background (sources) samples. In the case of the sources, we find typical mass values in agreement with previous studies: a minimum halo mass to host a central galaxy, $M_{min}\sim 10^{12.26} M_\odot$, and a pivot halo mass to have at least one sub-halo satellite, $M_1\sim 10^{12.84} M_\odot$. However, the lenses are massive galaxies or even galaxy groups/clusters, with minimum mass of $M_{min}^{lens}\sim 10^{13.06} M_\odot$. Above a mass of $M_1^{lens}\sim 10^{14.57} M_\odot$ they contain at least one additional satellite galaxy which contributes to the lensing effect. The tomographic analysis shows that, while $M_1^{lens}$ is almost redshift independent, there is a clear evolution of increase $M_{min}^{lens}$ with redshift in agreement with theoretical estimations. Finally, the halo modeling allows us to identify a strong lensing contribution to the cross-correlation for angular scales below 30". This interpretation is supported by the results of basic but effective simulations.

The Herschel Space Observatory has revealed a very different galaxyscape from that shown by optical surveys which presents a challenge for galaxy-evolution models. The Herschel surveys reveal (1) that there was rapid galaxy evolution in the very recent past and (2) that galaxies lie on a a single Galaxy Sequence (GS) rather than a star-forming `main sequence' and a separate region of `passive' or `red-and-dead' galaxies. The form of the GS is now clearer because far-infrared surveys such as the Herschel ATLAS pick up a population of optically-red star-forming galaxies that would have been classified as passive using most optical criteria. The space-density of this population is at least as high as the traditional star-forming population. By stacking spectra of H-ATLAS galaxies over the redshift range 0.001 < z < 0.4, we show that the galaxies responsible for the rapid low-redshift evolution have high stellar masses, high star-formation rates but, even several billion years in the past, old stellar populations - they are thus likely to be relatively recent ancestors of early-type galaxies in the Universe today. The form of the GS is inconsistent with rapid quenching models and neither the analytic bathtub model nor the hydrodynamical EAGLE simulation can reproduce the rapid cosmic evolution. We propose a new gentler model of galaxy evolution that can explain the new Herschel results and other key properties of the galaxy population.

We use the energy-balance code MAGPHYS to determine stellar and dust masses, and dust corrected star-formation rates for over 200,000 GAMA galaxies, 170,000 G10-COSMOS galaxies and 200,000 3D-HST galaxies. Our values agree well with previously reported measurements and constitute a representative and homogeneous dataset spanning a broad range in stellar mass (10^8---10^12 Msol), dust mass (10^6---10^9 Msol), and star-formation rates (0.01---100 Msol per yr), and over a broad redshift range (0.0 < z < 5.0). We combine these data to measure the cosmic star-formation history (CSFH), the stellar-mass density (SMD), and the dust-mass density (DMD) over a 12 Gyr timeline. The data mostly agree with previous estimates, where they exist, and provide a quasi-homogeneous dataset using consistent mass and star-formation estimators with consistent underlying assumptions over the full time range. As a consequence our formal errors are significantly reduced when compared to the historic literature. Integrating our cosmic star-formation history we precisely reproduce the stellar-mass density with an ISM replenishment factor of 0.50 +/- 0.07, consistent with our choice of Chabrier IMF plus some modest amount of stripped stellar mass. Exploring the cosmic dust density evolution, we find a gradual increase in dust density with lookback time. We build a simple phenomenological model from the CSFH to account for the dust mass evolution, and infer two key conclusions: (1) For every unit of stellar mass which is formed 0.0065---0.004 units of dust mass is also formed; (2) Over the history of the Universe approximately 90 to 95 per cent of all dust formed has been destroyed and/or ejected.

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.

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.

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

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

A common assumption is that galaxies fall in two distinct regions on a plot of specific star-formation rate (SSFR) versus galaxy stellar mass: a star-forming Galaxy Main Sequence (GMS) and a separate region of `passive' or `red and dead galaxies'. Starting from a volume-limited sample of nearby galaxies designed to contain most of the stellar mass in this volume, and thus being a fair representation of the Universe at the end of 12 billion years of galaxy evolution, we investigate the distribution of galaxies in this diagram today. We show that galaxies follow a strongly curved extended GMS with a steep negative slope at high galaxy stellar masses. There is a gradual change in the morphologies of the galaxies along this distribution, but there is no clear break between early-type and late-type galaxies. Examining the other evidence that there are two distinct populations, we argue that the `red sequence' is the result of the colours of galaxies changing very little below a critical value of the SSFR, rather than implying a distinct population of galaxies, and that Herschel observations, which show at least half of early-type galaxies contain a cool interstellar medium, also imply continuity between early-type and late-type galaxies. This picture of a unitary population of galaxies requires more gradual evolutionary processes than the rapid quenching processes needed to to explain two distinct populations. We challenge theorists to reproduce the properties of this `Galaxy End Sequence'.

We present a source-plane reconstruction of a ${\it Herschel}$ and ${\it Planck}$-detected gravitationally-lensed dusty star-forming galaxy (DSFG) at $z=1.68$ using {\it Hubble}, Sub-millimeter Array (SMA), and Keck observations. The background sub-millimeter galaxy (SMG) is strongly lensed by a foreground galaxy cluster at $z=0.997$ and appears as an arc of length $\sim 15^{\prime \prime}$ in the optical images. The continuum dust emission, as seen by SMA, is limited to a single knot within this arc. We present a lens model with source plane reconstructions at several wavelengths to show the difference in magnification between the stars and dust, and highlight the importance of a multi-wavelength lens models for studies involving lensed DSFGs. We estimate the physical properties of the galaxy by fitting the flux densities to model SEDs leading to a magnification-corrected star formation rate of $390 \pm 60$ M$_{\odot}$ yr$^{-1}$ and a stellar mass of $1.1 \pm 0.4\times 10^{11}$ M$_{\odot}$. These values are consistent with high-redshift massive galaxies that have formed most of their stars already. The estimated gas-to-baryon fraction, molecular gas surface density, and SFR surface density have values of $0.43 \pm 0.13$, $350 \pm 200$ M$_{\odot}$ pc$^{-2}$, and $\sim 12 \pm 7~$M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$, respectively. The ratio of star formation rate surface density to molecular gas surface density puts this among the most star-forming systems, similar to other measured SMGs and local ULIRGS.

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

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

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

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

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

In most studies of dust in galaxies, dust is only detected from its emission to approximately the optical radius of the galaxy. By combining the signal of 110 spiral galaxies observed as part of the Herschel Reference Survey, we are able to improve our sensitivity by an order-of-magnitude over that for a single object. Here we report the direct detection of dust from its emission that extends out to at least twice the optical radius. We find that the distribution of dust is consistent with an exponential at all radii with a gradient of ~-1.7 dex R$_{25}^{-1}$. Our dust temperature declines linearly from ~25 K in the centre to 15 K at R$_{25}$ from where it remains constant out to ~2.0 R$_{25}$. The surface-density of dust declines with radius at a similar rate to the surface-density of stars but more slowly than the surface-density of the star-formation rate. Studies based on dust extinction and reddening of high-redshift quasars have concluded that there are substantial amounts of dust in intergalactic space. By combining our results with the number counts and angular correlation function from the SDSS, we show that with Milky Way type dust we can explain the reddening of the quasars by the dust within galactic disks alone. Given the uncertainties in the properties of any intergalactic dust, we cannot rule out its existence, but our results show that statistical investigations of the dust in galactic halos that use the reddening of high-redshift objects must take account of the dust in galactic disks.

We present LOFAR High-Band Array (HBA) observations of the Herschel-ATLAS North Galactic Pole survey area. The survey we have carried out, consisting of four pointings covering around 142 square degrees of sky in the frequency range 126--173 MHz, does not provide uniform noise coverage but otherwise is representative of the quality of data to be expected in the planned LOFAR wide-area surveys, and has been reduced using recently developed `facet calibration' methods at a resolution approaching the full resolution of the datasets ($\sim 10 \times 6$ arcsec) and an rms off-source noise that ranges from 100 $\mu$Jy beam$^{-1}$ in the centre of the best fields to around 2 mJy beam$^{-1}$ at the furthest extent of our imaging. We describe the imaging, cataloguing and source identification processes, and present some initial science results based on a 5-$\sigma$ source catalogue. These include (i) an initial look at the radio/far-infrared correlation at 150 MHz, showing that many Herschel sources are not yet detected by LOFAR; (ii) number counts at 150 MHz, including, for the first time, observational constraints on the numbers of star-forming galaxies; (iii) the 150-MHz luminosity functions for active and star-forming galaxies, which agree well with determinations at higher frequencies at low redshift, and show strong redshift evolution of the star-forming population; and (iv) some discussion of the implications of our observations for studies of radio galaxy life cycles.

We present an improved and extended analysis of the cross-correlation between the map of the Cosmic Microwave Background (CMB) lensing potential derived from the \emph{Planck} mission data and the high-redshift galaxies detected by the \emph{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) in the photometric redshift range $z_{\rm ph} \ge 1.5$. We compare the results based on the 2013 and 2015 \textit{Planck} datasets, and investigate the impact of different selections of the H-ATLAS galaxy samples. Significant improvements over our previous analysis have been achieved thanks to the higher signal-to-noise ratio of the new CMB lensing map recently released by the \textit{Planck} collaboration. The effective galaxy bias parameter, $b$, for the full galaxy sample, derived from a joint analysis of the cross-power spectrum and of the galaxy auto-power spectrum is found to be $b = 3.54^{+0.15}_{-0.14}$. Furthermore, a first tomographic analysis of the cross-correlation signal is implemented, by splitting the galaxy sample into two redshift intervals: $1.5 \le z_{\rm ph} < 2.1$ and $z_{\rm ph}\ge 2.1$. A statistically significant signal was found for both bins, indicating a substantial increase with redshift of the bias parameter: $b=2.89\pm0.23$ for the lower and $b=4.75^{+0.24}_{-0.25}$ for the higher redshift bin. Consistently with our previous analysis we find that the amplitude of the cross correlation signal is a factor of $1.45^{+0.14}_{-0.13}$ higher than expected from the standard $\Lambda$CDM model for the assumed redshift distribution. The robustness of our results against possible systematic effects has been extensively discussed although the tension is mitigated by passing from 4 to 3$\sigma$.

We investigate the region around the Planck-detected z=3.26 gravitationally lensed galaxy HATLAS J114637.9-001132 (hereinafter HATLAS12-00) using both archival Herschel data from the H-ATLAS survey and using submm data obtained with both LABOCA and SCUBA2. The lensed source is found to be surrounded by a strong overdensity of both Herschel-SPIRE sources and submm sources. We detect 17 bright (S_870 >~7 mJy) sources at >4sigma closer than 5 arcmin to the lensed object at 850/870 microns. Ten of these sources have good cross-identifications with objects detected by Herschel-SPIRE which have redder colours than other sources in the field, with 350 micron flux > 250 micron flux, suggesting that they lie at high redshift. Submillimeter Array (SMA) observations localise one of these companions to ~1 arcsecond, allowing unambiguous cross identification with a 3.6 and 4.5 micron Spitzer source. The optical/near-IR spectral energy distribution (SED) of this source is measured by further observations and found to be consistent with z>2, but incompatible with lower redshifts. We conclude that this system may be a galaxy cluster/protocluster or larger scale structure that contains a number of galaxies undergoing starbursts at the same time.

We compare the absolute gain photometric calibration of the Planck/HFI and Herschel/SPIRE instruments on diffuse emission. The absolute calibration of HFI and SPIRE each relies on planet flux measurements and comparison with theoretical far-infrared emission models of planetary atmospheres. We measure the photometric cross calibration between the instruments at two overlapping bands, 545 GHz / 500 $\mu$m and 857 GHz / 350 $\mu$m. The SPIRE maps used have been processed in the Herschel Interactive Processing Environment (Version 12) and the HFI data are from the 2015 Public Data Release 2. For our study we used 15 large fields observed with SPIRE, which cover a total of about 120 deg^2. We have selected these fields carefully to provide high signal-to-noise ratio, avoid residual systematics in the SPIRE maps, and span a wide range of surface brightness. The HFI maps are bandpass-corrected to match the emission observed by the SPIRE bandpasses. The SPIRE maps are convolved to match the HFI beam and put on a common pixel grid. We measure the cross-calibration relative gain between the instruments using two methods in each field, pixel-to-pixel correlation and angular power spectrum measurements. The SPIRE / HFI relative gains are 1.047 ($\pm$ 0.0069) and 1.003 ($\pm$ 0.0080) at 545 and 857 GHz, respectively, indicating very good agreement between the instruments. These relative gains deviate from unity by much less than the uncertainty of the absolute extended emission calibration, which is about 6.4% and 9.5% for HFI and SPIRE, respectively, but the deviations are comparable to the values 1.4% and 5.5% for HFI and SPIRE if the uncertainty from models of the common calibrator can be discounted. Of the 5.5% uncertainty for SPIRE, 4% arises from the uncertainty of the effective beam solid angle, which impacts the adopted SPIRE point source to extended source unit conversion factor (Abridged)

We present a list of candidate gravitationally lensed dusty star-forming galaxies (DSFGs) from the HerMES Large Mode Survey (HeLMS) and the Herschel Stripe 82 Survey (HerS). Together, these partially overlapping surveys cover 372 deg$^{2}$ on the sky. After removing local spiral galaxies and known radio-loud blazars, our candidate list of lensed DSFGs is composed of 77 sources with 500 $\mu$m flux densities ($S_{500}$) greater than 100 mJy. Such sources are dusty starburst galaxies similar to the first bright Sub Millimeter Galaxies (SMGs) discovered with SCUBA. We expect a large fraction of this list to be strongly lensed, with a small fraction made up of bright SMG-SMG mergers that appear as Hyper-Luminous Infrared Galaxies (HyLIRGs; $\rm L_{IR}>10^{13} L_{\odot}$). Thirteen of the 77 candidates have spectroscopic redshifts from CO spectroscopy with ground-based interferometers, putting them at $z>1$ and well above the redshift of the foreground lensing galaxies. The surface density of our sample is 0.21 $\pm$ 0.03 deg$^{-2}$. We present follow-up imaging of a few of the candidates confirming their lensing nature. The sample presented here is an ideal tool for higher resolution imaging and spectroscopic observations to understand detailed properties of starburst phenomena in distant galaxies.

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

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

The B-mode Foreground Experiment (BFORE) is a proposed NASA balloon project designed to make optimal use of the sub-orbital platform by concentrating on three dust foreground bands (270, 350, and 600 GHz) that complement ground-based cosmic microwave background (CMB) programs. BFORE will survey ~1/4 of the sky with 1.7 - 3.7 arcminute resolution, enabling precise characterization of the Galactic dust that now limits constraints on inflation from CMB B-mode polarization measurements. In addition, BFORE's combination of frequency coverage, large survey area, and angular resolution enables science far beyond the critical goal of measuring foregrounds. BFORE will constrain the velocities of thousands of galaxy clusters, provide a new window on the cosmic infrared background, and probe magnetic fields in the interstellar medium. We review the BFORE science case, timeline, and instrument design, which is based on a compact off-axis telescope coupled to >10,000 superconducting detectors.

The dust in nearby galaxies absorbs a fraction of the UV-optical-near-infrared radiation produced by stars. This energy is consequently re-emitted in the infrared. We investigate the portion of the stellar radiation absorbed by spiral galaxies from the HRS by modelling their UV-to-submillimetre spectral energy distributions. Our models provide an attenuated and intrinsic SED from which we find that on average 32 % of all starlight is absorbed by dust. We define the UV heating fraction as the percentage of dust luminosity that comes from absorbed UV photons and find that this is 56 %, on average. This percentage varies with morphological type, with later types having significantly higher UV heating fractions. We find a strong correlation between the UV heating fraction and specific star formation rate and provide a power-law fit. Our models allow us to revisit the IRX-AFUV relations, and derive these quantities directly within a self-consistent framework. We calibrate this relation for different bins of NUV-r colour and provide simple relations to relate these parameters. We investigated the robustness of our method and we conclude that the derived parameters are reliable within the uncertainties which are inherent to the adopted SED model. This calls for a deeper investigation on how well extinction and attenuation can be determined through panchromatic SED modelling.

We combine Herschel/SPIRE sub-millimeter (submm) observations with existing multi-wavelength data to investigate the characteristics of low redshift, optically red galaxies detected in submm bands. We select a sample of galaxies in the redshift range 0.01$\leq$z$\leq$0.2, having >5$\sigma$ detections in the SPIRE 250 micron submm waveband. Sources are then divided into two sub-samples of $red$ and $blue$ galaxies, based on their UV-optical colours. Galaxies in the $red$ sample account for $\approx$4.2 per cent of the total number of sources with stellar masses M$_{*}\gtrsim$10$^{10}$ Solar-mass. Following visual classification of the $red$ galaxies, we find that $\gtrsim$30 per cent of them are early-type galaxies and $\gtrsim$40 per cent are spirals. The colour of the $red$-spiral galaxies could be the result of their highly inclined orientation and/or a strong contribution of the old stellar population. It is found that irrespective of their morphological types, $red$ and $blue$ sources occupy environments with more or less similar densities (i.e., the $\Sigma_5$ parameter). From the analysis of the spectral energy distributions (SEDs) of galaxies in our samples based on MAGPHYS, we find that galaxies in the $red$ sample (of any morphological type) have dust masses similar to those in the $blue$ sample (i.e. normal spiral/star-forming systems). However, in comparison to the $red$-spirals and in particular $blue$ systems, $red$-ellipticals have lower mean dust-to-stellar mass ratios. Besides galaxies in the $red$-elliptical sample have much lower mean star-formation/specific-star-formation rates in contrast to their counterparts in the $blue$ sample. Our results support a scenario where dust in early-type systems is likely to be of an external origin.

We present the GAMA Panchromatic Data Release (PDR) constituting over 230deg$^2$ of imaging with photometry in 21 bands extending from the far-UV to the far-IR. These data complement our spectroscopic campaign of over 300k galaxies, and are compiled from observations with a variety of facilities including: GALEX, SDSS, VISTA, WISE, and Herschel, with the GAMA regions currently being surveyed by VST and scheduled for observations by ASKAP. These data are processed to a common astrometric solution, from which photometry is derived for 221,373 galaxies with r<19.8 mag. Online tools are provided to access and download data cutouts, or the full mosaics of the GAMA regions in each band. We focus, in particular, on the reduction and analysis of the VISTA VIKING data, and compare to earlier datasets (i.e., 2MASS and UKIDSS) before combining the data and examining its integrity. Having derived the 21-band photometric catalogue we proceed to fit the data using the energy balance code MAGPHYS. These measurements are then used to obtain the first fully empirical measurement of the 0.1-500$\mu$m energy output of the Universe. Exploring the Cosmic Spectral Energy Distribution (CSED) across three time-intervals (0.3-1.1Gyr, 1.1-1.8~Gyr and 1.8---2.4~Gyr), we find that the Universe is currently generating $(1.5 \pm 0.3) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$, down from $(2.5 \pm 0.2) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$ 2.3~Gyr ago. More importantly, we identify significant and smooth evolution in the integrated photon escape fraction at all wavelengths, with the UV escape fraction increasing from 27(18)% at z=0.18 in NUV(FUV) to 34(23)% at z=0.06. The GAMA PDR will allow for detailed studies of the energy production and outputs of individual systems, sub-populations, and representative galaxy samples at $z<0.5$. The GAMA PDR can be found at: http://gama-psi.icrar.org/

In this work we present IRAM-30m telescope observations of a sample of bulge-dominated galaxies with large dust lanes, which have had a recent minor merger. We find these galaxies are very gas rich, with H2 masses between 4x10^8 and 2x10^10 Msun. We use these molecular gas masses, combined with atomic gas masses from an accompanying paper, to calculate gas-to-dust and gas-to-stellar mass ratios. The gas-to-dust ratios of our sample objects vary widely (between ~50 and 750), suggesting many objects have low gas-phase metallicities, and thus that the gas has been accreted through a recent merger with a lower mass companion. We calculate the implied minor companion masses and gas fractions, finding a median predicted stellar mass ratio of ~40:1. The minor companion likely had masses between ~10^7 - 10^10 Msun. The implied merger mass ratios are consistent with the expectation for low redshift gas-rich mergers from simulations. We then go on to present evidence that (no matter which star-formation rate indicator is used) our sample objects have very low star-formation efficiencies (star-formation rate per unit gas mass), lower even than the early-type galaxies from ATLAS3D which already show a suppression. This suggests that minor mergers can actually suppress star-formation activity. We discuss mechanisms that could cause such a suppression, include dynamical effects induced by the minor merger.

Using results from the Herschel Astrophysical Terrahertz Large-Area Survey and the Galaxy and Mass Assembly project, we show that, for galaxy masses above approximately 1.0e8 solar masses, 51% of the stellar mass-density in the local Universe is in early-type galaxies (ETGs: Sersic n > 2.5) while 89% of the rate of production of stellar mass-density is occurring in late-type galaxies (LTGs: Sersic n < 2.5). From this zero-redshift benchmark, we have used a calorimetric technique to quantify the importance of the morphological transformation of galaxies over the history of the Universe. The extragalactic background radiation contains all the energy generated by nuclear fusion in stars since the Big Bang. By resolving this background radiation into individual galaxies using the deepest far-infrared survey with the Herschel Space Observatory and a deep near-infrared/optical survey with the Hubble Space Telescope (HST), and using measurements of the Sersic index of these galaxies derived from the HST images, we estimate that approximately 83% of the stellar mass-density formed over the history of the Universe occurred in LTGs. The difference between this and the fraction of the stellar mass-density that is in LTGs today implies there must have been a major transformation of LTGs into ETGs after the formation of most of the stars.

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

We present the properties of the first 250 $\mu$m blind sample of nearby galaxies (15 < D < 46 Mpc) containing 42 objects from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Herschel's sensitivity probes the faint end of the dust luminosity function for the first time, spanning a range of stellar mass (7.4 < log$_{10}$ M$_{\star}$ < 11.3 M$_{\odot}$), star formation activity (-11.8 < log$_{10}$ SSFR < -8.9 yr$^{-1}$), gas fraction (3-96 per cent), and colour (0.6 < FUV-Ks < 7.0 mag). The median cold dust temperature is 14.6 K, colder than in the Herschel Reference Survey (18.5 K) and Planck Early Release Compact Source Catalogue 17.7 K. The mean dust-to-stellar mass ratio in our sample is higher than these surveys by factors of 3.7 and 1.8, with a dust mass volume density of (3.7 $\pm$ 0.7) x 10$^{5}$ M$_{\odot}$ Mpc$^{-3}$. Counter-intuitively, we find that the more dust rich a galaxy, the lower its UV attenuation. Over half of our dust-selected sample are very blue in FUV-Ks colour, with irregular and/or highly flocculent morphology, these galaxies account for only 6 per cent of the sample's stellar mass but contain over 35 per cent of the dust mass. They are the most actively star forming galaxies in the sample, with the highest gas fractions and lowest UV attenuation. They also appear to be in an early stage of converting their gas into stars, providing valuable insights into the chemical evolution of young galaxies.

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

We exploit long-baseline ALMA sub-mm observations of the lensed star-forming galaxy SDP 81 at z=3.042 to investigate the properties of inter-stellar medium on scales of 50-100pc. The kinematics of the CO gas within this system are well described by a rotationally-supported disk with an inclination-corrected rotation speed, v=320+/-20km/s and a dynamical mass of M=(3.5+/-1.0)x10^10Mo within a radius of 1.5 kpc. The disk is gas rich and unstable, with a Toomre parameter, Q=0.30+/-0.10 and so should collapse in to star-forming regions with Jeans length L_J~130pc. We identify five star-forming regions within the ISM on these scales and show that their scaling relations between luminosity, line-widths and sizes are significantly offset from those typical of molecular clouds in local Galaxies (Larson's relations). These offsets are likely to be caused by the high external hydrostatic pressure for the interstellar medium (ISM), P/kB=(40+/-20)x10^7K/cm3, which is ~10,000x higher than the typical ISM pressure in the Milky Way. The physical conditions of the star-forming ISM and giant molecular clouds appears to be similar to the those found in the densest environments in the local Universe, such as those in the Galactic center.

Interstellar dust in galaxies can be traced either through its extinction effects on the star light, or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of three. We investigate the dust energy balance for IC4225 and NGC5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC4225 underestimates the observed fluxes by a factor of about three. For NGC5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy.

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

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

We explore the behaviour of [CII]-157.74um forbidden fine-structure line observed in a sample of 28 galaxies selected from ~50deg^2 of the H-ATLAS survey. The sample is restricted to galaxies with flux densities higher than S_160um>150mJy and optical spectra from the GAMA survey at 0.02<z<0.2. Far-IR spectra centred on this redshifted line were taken with the PACS instrument on-board the Herschel Space Observatory. The galaxies span 10<log(L_IR/Lo)<12 (where L_IR=L_IR[8-1000um]) and 7.3<log(L_[CII]/Lo)<9.3, covering a variety of optical galaxy morphologies. The sample exhibits the so-called [CII] deficit at high IR luminosities, i.e. L_[CII]/L_IR (hereafter [CII]/IR) decreases at high L_IR. We find significant differences between those galaxies presenting [CII]/IR>2.5x10^-3 with respect to those showing lower ratios. In particular, those with high ratios tend to have: (1) L_IR<10^11Lo; (2) cold dust temperatures, T_d<30K; (3) disk-like morphologies in r-band images; (4) a WISE colour 0.5<S_12um/S_22um<1.0; (5) low surface brightness Sigma_IR~10^8-9 Lo kpc^-2, (6) and specific star-formation rates of sSFR~0.05-3 Gyr^-1. We suggest that the strength of the far-UV radiation fields (<G_O>) is main parameter responsible for controlling the [CII]/IR ratio. It is possible that relatively high <G_O> creates a positively charged dust grain distribution, impeding an efficient photo-electric extraction of electrons from these grains to then collisionally excite carbon atoms. Within the brighter IR population, 11<log(L_IR/Lo)<12, the low [CII]/IR ratio is unlikely to be modified by [CII] self absorption or controlled by the presence of a moderately luminous AGN (identified via the BPT diagram).

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Gamma-ray bursts (GRBs) are the most energetic phenomena in the Universe; believed to result from the collapse and subsequent explosion of massive stars. Even though it has profound consequences for our understanding of their nature and selection biases, little is known about the dust properties of the galaxies hosting GRBs. We present analysis of the far-infrared properties of an unbiased sample of 20 \textit{BeppoSAX} and \textit{Swift} GRB host galaxies (at an average redshift of $z\,=\,3.1$) located in the {\it Herschel} Astrophysical Terahertz Large Area Survey, the {\it Herschel} Virgo Cluster Survey, the {\it Herschel} Fornax Cluster Survey, the {\it Herschel} Stripe 82 Survey and the {\it Herschel} Multi-tiered Extragalactic Survey, totalling $880$ deg$^2$, or $\sim 3$\% of the sky in total. Our sample selection is serendipitous, based only on whether the X-ray position of a GRB lies within a large-scale {\it Herschel} survey -- therefore our sample can be considered completely unbiased. Using deep data at wavelengths of 100\,--\,500$\,\mu$m, we tentatively detected 1 out of 20 GRB hosts located in these fields. We constrain their dust masses and star formation rates (SFRs), and discuss these in the context of recent measurements of submillimetre galaxies and ultraluminous infrared galaxies. The average far-infrared flux of our sample gives an upper limit on SFR of $<114\,{\rm M}\odot\,\mbox{yr}^{-1}$. The detection rate of GRB hosts is consistent with that predicted assuming that GRBs trace the cosmic SFR density in an unbiased way, i.e. that the fraction of GRB hosts with $\mbox{SFR}>500\,{\rm M}\odot\,\mbox{yr}^{-1}$ is consistent with the contribution of such luminous galaxies to the cosmic star formation density.

Surveys with submillimetre telescopes are revealing large numbers of gravitationally lensed high-redshift sources. I describe how, in practice, these lensed systems could be simultaneously used to estimate the values of cosmological parameters, test models for the evolution of the distribution of dark-matter halos and investigate the properties of the source population. Even the existing sample of lenses found with the Hershcel Space Observatory is enough to formally rule out the standard models of the evolving population of dark-matter halos, with the likely explanation a combination of baryon physics and the perturbation by infalling baryons of the density distribution of dark matter at the centres of the halos. Independently of the evolution of the halos, observations of a sample of 100 lensed systems would be enough to estimate Omega_Lambda with a precision of 5% and observations of 1000 lenses would be enough to estimate w, the parameter in the equation-of-state of dark energy, with a precision similar to that obtained from the Planck observations of the cosmic microwave background. While the fraction of submillimetre sources that are lensed depends weakly on the specific halo mass function that is used in the model, it depends very strongly on the evolution of the submillimetre luminosity function of the source population. Therefore measurements of the lensing fraction could be used to investigate galaxy evolution in a way that is independent of the properties of the intervening halos.

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

We use 10,387 galaxies from the Herschel Astrophysical TeraHertz Large Area Survey (H-ATLAS) to probe the far-infrared radio correlation (FIRC) of star forming galaxies as a function of redshift, wavelength, and effective dust temperature. All of the sources in our 250 {\mu}m-selected sample have spectroscopic redshifts, as well as 1.4 GHz flux density estimates measured from the Faint Images of the Radio Sky at Twenty centimetres (FIRST) survey. This enables us to study not only individual sources, but also the average properties of the 250 {\mu}m selected population using median stacking techniques. We find that individual sources detected at $\geq 5\sigma$ in both the H-ATLAS and FIRST data have logarithmic flux ratios (i.e. FIRC $q_\lambda$ parameters) consistent with previous studies of the FIRC. In contrast, the stacked values show larger $q_\lambda$, suggesting excess far-IR flux density/luminosity in 250{\mu}m selected sources above what has been seen in previous analyses. In addition, we find evidence that 250 {\mu}m sources with warm dust SEDs have a larger 1.4 GHz luminosity than the cooler sources in our sample. Though we find no evidence for redshift evolution of the monochromatic FIRC, our analysis reveals significant temperature dependence. Whilst the FIRC is reasonably constant with temperature at 100 {\mu}m, we find increasing inverse correlation with temperature as we probe longer PACS and SPIRE wavelengths. These results may have important implications for the use of monochromatic dust luminosity as a star formation rate indicator in star-forming galaxies, and in the future, for using radio data to determine galaxy star formation rates.

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

In this paper we present a catalogue of Giant Molecular Clouds (GMCs) in the Andromeda (M31) galaxy extracted from the Hershel Exploitation of Local Galaxy Andromeda (HELGA) dataset. GMCs are identified from the Herschel maps using a hierarchical source extraction algorithm. We present the results of this new catalogue and characterise the spatial distribution and spectral energy properties of its clouds based on the radial dust/gas properties found by Smith et al (2012). 236 GMCs in the mass range 10^4-10^7 M_sol are identified, their cumulative mass distribution is found to be proportional to M^-1.45 in agreement with earlier studies. The GMCs appear to follow the same cloud mass to L_CO correlation observed in the Milky Way. However, comparison between this catalogue and interferometry studies also shows that the GMCs are substructured below the Herschel resolution limit suggesting that we are observing associations of GMCs. Following Gordon et al. (2006), we study the spatial structure of M31 by splitting the observed structure into a set of spiral arms and offset rings. We fit radii of 10.5 and 15.5 kpc to the two most prominent rings. We then fit a logarithmic spiral with a pitch angle of 8.9 deg to the GMCs not associated with either ring. Lastly, we comment upon the effects of deprojection on our results and investigate the effect different models for M31's inclination will have upon the projection of an unperturbed spiral arm system.

We use the Herschel ATLAS (H-ATLAS) Phase I data to study the conditional luminosity function of far-IR (250 um) selected galaxies in optically-selected galaxy groups from the GAMA spectroscopic survey, as well as environmental effects on the far-IR-to-optical colour. We applied two methods, which gave consistent results for the far-IR conditional luminosity functions. The direct matching method matches H-ATLAS sources to GAMA/SDSS galaxies, then links the optical counterparts to GAMA groups. The stacking method counts the number of far-IR sources within the projected radii of GAMA groups, subtracting the local background. We investigated the dependence of the far-IR (250 um) luminosity function on group mass in the range 10d12 < M_h < 10d14Msun/h and on redshift in the range 0 < z < 0.4, using a sample of 3000 groups containing H-ATLAS sources with GAMA redshifts over an area of 126 sq. deg. We find that the characteristic 250 um luminosity, L^*(250), increases with group mass up to M_h ~ 10d13Msun/h, but is roughly constant above this, while it increases with redshift at high group masses, but less so at low masses. We also find that the group far-IR luminosity-to-mass ratio L(250)/M_h increases with redshift and is higher in low-mass groups. We estimate that around 70% of the 250 um luminosity density in the local universe is contributed by groups with M_h > 10d12Msun/h. We also find that the far-IR-to-optical colours of H-ATLAS galaxies are independent of group mass over the range 10d12 < M_h <10d14Msun/h in the local universe. We also compare our observational results with recent semi-analytical models, and find that none of these galaxy formation model can reproduce the conditional far-IR luminosity functions of galaxy groups.

[Abridged] Aims: This work focuses on one lensed system, HATLAS J142935.3-002836 (H1429-0028), selected in the Herschel-ATLAS field. Gathering a rich, multi-wavelength dataset, we aim to confirm the lensing hypothesis and model the background source's morphology and dynamics, as well as to provide a full physical characterisation. Methods: Multi-wavelength high-resolution data is utilised to assess the nature of the system. A lensing-analysis algorithm which simultaneously fits different wavebands is adopted to characterise the lens. The background galaxy dynamical information is studied by reconstructing the 3-D source-plane of the ALMA CO(J:4-3) transition. Near-IR imaging from HST and Keck-AO allows to constrain rest-frame optical photometry independently for the foreground and background systems. Physical parameters (such as stellar and dust masses) are estimated via modelling of the spectral energy distribution taking into account source blending, foreground obscuration, and differential magnification. Results: The system comprises a foreground edge-on disk galaxy (at z_sp=0.218) with an almost complete Einstein ring around it. The background source (at z_sp=1.027) is magnified by a factor of ~8-10 depending on wavelength. It is comprised of two components and a tens of kpc long tidal tail resembling the Antennae merger. As a whole, the system is a massive stellar system (1.32[-0.41,+0.63] x1E11 Mo) forming stars at a rate of 394+-90 Mo/yr, and has a significant gas reservoir M_ISM = 4.6+-1.7 x1E10 Mo. Its depletion time due to star formation alone is thus expected to be tau_SF=M_ISM/SFR=117+-51 Myr. The dynamical mass of one of the components is estimated to be 5.8+-1.7 x1E10 Mo, and, together with the photometric total mass estimate, it implies that H1429-0028 is a major merger system (1:2.8[-1.5,+1.8]).

We report a highly significant ($>10\sigma$) spatial correlation between galaxies with $S_{350\mu\rm m}\ge 30\,$mJy detected in the equatorial fields of the \textsl{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) with estimated redshifts $\gtrsim 1.5$, and SDSS or GAMA galaxies at $0.2\le z\le 0.6$. The significance of the cross-correlation is much higher than those reported so far for samples with non-overlapping redshift distributions selected in other wavebands. Extensive, realistic simulations of clustered sub-mm galaxies amplified by foreground structures confirm that the cross-correlation is explained by weak gravitational lensing ($\mu<2$). The simulations also show that the measured amplitude and range of angular scales of the signal are larger than can be accounted for by galaxy-galaxy weak lensing. However, for scales $\lesssim 2\,$arcmin, the signal can be reproduced if SDSS/GAMA galaxies act as signposts of galaxy groups/clusters with halo masses in the range $10^{13.2}$--$10^{14.5} M_{\odot}$. The signal detected on larger scales appears to reflect the clustering of such halos.

We present a comparison of the physical properties of a rest-frame $250\mu$m selected sample of massive, dusty galaxies from $0<z<5.3$. Our sample comprises 29 high-redshift submillimetre galaxies (SMGs) from the literature, and 843 dusty galaxies at $z<0.5$ from the Herschel-ATLAS, selected to have a similar stellar mass to the SMGs. The $z>1$ SMGs have an average SFR of $390^{+80}_{-70}\,$M$_\odot$yr$^{-1}$ which is 120 times that of the low-redshift sample matched in stellar mass to the SMGs (SFR$=3.3\pm{0.2}$ M$_\odot$yr$^{-1}$). The SMGs harbour a substantial mass of dust ($1.2^{+0.3}_{-0.2}\times{10}^9\,$M$_\odot$), compared to $(1.6\pm0.1)\times{10}^8\,$M$_\odot$ for low-redshift dusty galaxies. At low redshifts the dust luminosity is dominated by the diffuse ISM, whereas a large fraction of the dust luminosity in SMGs originates from star-forming regions. At the same dust mass SMGs are offset towards a higher SFR compared to the low-redshift H-ATLAS galaxies. This is not only due to the higher gas fraction in SMGs but also because they are undergoing a more efficient mode of star formation, which is consistent with their bursty star-formation histories. The offset in SFR between SMGs and low-redshift galaxies is similar to that found in CO studies, suggesting that dust mass is as good a tracer of molecular gas as CO.

We present new deep images of the Coma Cluster from the ESA Herschel Space Observatory at wavelengths of 70, 100 and 160 microns, covering an area of 1.75 x 1.0 square degrees encompassing the core and southwest infall region. Our data display an excess of sources at flux densities above 100 mJy compared to blank-field surveys, as expected. We use extensive optical spectroscopy of this region to identify cluster members and hence produce cluster luminosity functions in all three photometric bands. We compare our results to the local field galaxy luminosity function, and the luminosity functions from the Herschel Virgo Cluster Survey (HeViCS). We find consistency between the shapes of the Coma and field galaxy luminosity functions at all three wavelengths, however we do not find the same level of agreement with that of the Virgo Cluster.

We have constructed a sample of radio-loud and radio-quiet quasars from the Faint Im- ages Radio Sky at Twenty-one centimetres (FIRST) and the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), over the H-ATLAS Phase 1 Area (9h, 12h and 14.5h). Using a stacking analysis we find a significant correlation between the far-infrared luminosity and 1.4-GHz luminosity for radio-loud quasars. Partial correlation analysis confirms the intrinsic correlation after removing the redshift contribution while for radio-quiet quasars no partial correlation is found. Using a single-temperature grey-body model we find a general trend of lower dust temperatures in the case of radio-loud quasars comparing to radio-quiet quasars. Also, radio-loud quasars are found to have almost constant mean values of dust mass along redshift and optical luminosity bins. In addition, we find that radio-loud quasars at lower optical luminosities tend to have on average higher FIR and 250-micron luminosity with respect to radio-quiet quasars with the same optical luminosites. Even if we use a two-temperature grey-body model to describe the FIR data, the FIR luminosity excess remains at lower optical luminosities. These results suggest that powerful radio jets are associated with star formation especially at lower accretion rates.

We have determined the mass-density radial profiles of the first five strong gravitational lens systems discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We present an enhancement of the semi-linear lens inversion method of Warren & Dye which allows simultaneous reconstruction of several different wavebands and apply this to dual-band imaging of the lenses acquired with the Hubble Space Telescope. The five systems analysed here have lens redshifts which span a range, 0.22<z<0.94. Our findings are consistent with other studies by concluding that: 1) the logarithmic slope of the total mass density profile steepens with decreasing redshift; 2) the slope is positively correlated with the average total projected mass density of the lens contained within half the effective radius and negatively correlated with the effective radius; 3) the fraction of dark matter contained within half the effective radius increases with increasing effective radius and increases with redshift.

We present Herschel/PACS 100 and 160 micron integrated photometry for the 323 galaxies in the Herschel Reference Survey (HRS), a K-band-, volume-limited sample of galaxies in the local Universe. Once combined with the Herschel/SPIRE observations already available, these data make the HRS the largest representative sample of nearby galaxies with homogeneous coverage across the 100-500 micron wavelength range. In this paper, we take advantage of this unique dataset to investigate the properties and shape of the far-infrared/sub-millimeter spectral energy distribution in nearby galaxies. We show that, in the stellar mass range covered by the HRS (8<log(M*/Msun)<12), the far-infrared/sub-millimeter colours are inconsistent with a single modified black-body having the same dust emissivity index beta for all galaxies. In particular, either beta decreases, or multiple temperature components are needed, when moving from metal-rich/gas-poor to metal-poor/gas-rich galaxies. We thus investigate how the dust temperature and mass obtained from a single modified black-body depend on the assumptions made on beta. We show that, while the correlations between dust temperature, galaxy structure and star formation rate are strongly model dependent, the dust mass scaling relations are much more reliable, and variations of beta only change the strength of the observed trends.

We gather infrared (IR) photometric data from 8 to 500 microns (Spitzer, WISE, IRAS and Herschel) for all of the HRS galaxies. Draine & Li (2007) models are fit to the data from which the stellar contribution has been carefully removed. We find that our photometric coverage is sufficient to constrain all of the models parameters and that a strong constraint on the 20-60 microns range is mandatory to estimate the relative contribution of the photo-dissociation regions to the IR SED. The SED models tend to systematically under-estimate the observed 500 microns flux densities, especially for low mass systems. We provide the output parameters for all of the galaxies: the minimum intensity of the interstellar radiation field (ISRF), the fraction of PAH, the relative contribution of PDR and evolved stellar population to the dust heating, the $M_{dust}$ and the $L_{IR}$. For a subsample of gas-rich galaxies, we analyze the relations between these parameters and the integrated properties of galaxies, such as $M_*$, SFR, metallicity, H$\alpha$ and H-band surface brightness, and the FUV attenuation. A good correlation between the fraction of PAH and the metallicity is found implying a weakening of the PAH emission in galaxies with low metallicities. The intensity of the IRSF and the H-band and H$\alpha$ surface brightnesses are correlated, suggesting that the diffuse dust component is heated by both the young stars in star forming regions and the diffuse evolved population. We use these results to provide a new set of IR templates calibrated with Herschel observations on nearby galaxies and a mean SED template to provide the z=0 reference for cosmological studies. For the same purpose, we put our sample on the SFR-$M_*$ diagram. The templates are compared to the most popular IR SED libraries, enlightening a large discrepancy between all of them in the 20-100 microns range.

We report on deep near-infrared observations obtained with the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) of the first five confirmed gravitational lensing events discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We succeed in disentangling the background galaxy from the lens to gain separate photometry of the two components. The HST data allow us to significantly improve on previous constraints of the mass in stars of the lensed galaxy and to perform accurate lens modelling of these systems, as described in the accompanying paper by Dye et al. We fit the spectral energy distributions of the background sources from near-IR to millimetre wavelengths and use the magnification factors estimated by Dye et al. to derive the intrinsic properties of the lensed galaxies. We find these galaxies to have star-formation rates of approximately 400 to 2000 M_sol/yr, with approximately (6-25)x10^10 M_sol of their baryonic mass already turned into stars. At these rates of star formation, all remaining molecular gas will be exhausted in less than 100 Myr, reaching a final mass in stars of a few 10^11 M_sol. These galaxies are thus proto-ellipticals caught during their major episode of star formation, and observed at the peak epoch z=1.5-3 of the cosmic star formation history of the Universe.

We use deep Herschel PACS and SPIRE observations in GOODSS, GOODSN and COSMOS to estimate the average dust mass (Mdust) of galaxies on a redshift-stellar mass (Mstar)-SFR grid. We study the scaling relations between Mdust, Mstar and SFR at z<=2.5. No clear evolution of Mdust is observed at fixed SFR and Mstar. We find a tight correlation between SFR and Mdust, likely a consequence of the Schmidt-Kennicutt (S-K) law. The Mstar-Mdust correlation observed by previous works flattens or sometimes disappears when fixing the SFR. Most of it likely derives from the combination of the Mdust-SFR and Mstar-SFR correlations. We then investigate the gas content as inferred by converting Mdust by assuming that the dust/gas ratio scales linearly with the gas metallicity. All galaxies in the sample follow, within uncertainties, the same SFR-Mgas relation (integrated S-K law), which broadly agrees with CO-based results for the bulk of the population, despite the completely different approaches. The majority of galaxies at z~2 form stars with an efficiency (SFE=SFR/Mgas) ~5 times higher than at z~0. It is not clear what fraction of such variation is an intrinsic redshift evolution and what fraction arises from selection effects. The gas fraction (fgas) decreases with Mstar and increases with SFR, and does not evolve with z at fixed Mstar and SFR. We explain these trends by introducing a universal relation between fgas, Mstar and SFR, non-evolving out to z~2.5. Galaxies move across this relation as their gas content evolves in time. We use the 3D fundamental fgas-Mstar-SFR relation and the redshift evolution of the Main Sequence to estimate the evolution of fgas in the average population of galaxies as a function of z and Mstar, and we find evidence a downsizing scenario.

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

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

We use galaxies from the Herschel-ATLAS survey, and a suite of ancillary simulations based on an isothermal dust model, to study our ability to determine the effective dust temperature, luminosity and emissivity index of 250um selected galaxies in the local Universe (z < 0.5). As well as simple far-infrared SED fitting of individual galaxies based on chi^2 minimisation, we attempt to derive the best global isothermal properties of 13,826 galaxies with reliable optical counterparts and spectroscopic redshifts. Using our simulations, we highlight the fact that applying traditional SED fitting techniques to noisy observational data in the Herschel Space Observatory bands introduces artificial anti-correlation between derived values of dust temperature and emissivity index. This is true even for galaxies with the most robust detections in our sample, making the results hard to interpret. We apply a method to determine the best-fit global values of isothermal effective temperature and emissivity index for z < 0.5 galaxies in H-ATLAS, deriving T = 22.3 +/- 0.1K and Beta = 1.98 +/- 0.02 (or T = 23.5 +/- 0.1K and Beta = 1.82 +/- 0.02 if we attempt to correct for bias by assuming that T and Beta are independent and normally distributed). We use our technique to test for an evolving emissivity index, finding only weak evidence. The median dust luminosity of our sample is log(Ldust/Lsolar) = 10.72 +/- 0.05 which (unlike T) shows little dependence on the choice of Beta used in our analysis, including whether it is variable or fixed. We use a further suite of simulations to emphasise the importance of the H-ATLAS PACS data for deriving dust temperatures at these redshifts, even though they are less sensitive than the SPIRE data. The majority of galaxies detected by H-ATLAS are normal star-forming galaxies, though a substantial minority (~31%) fall in the Luminous Infrared Galaxy category. (Abridged)

Upon its completion the Herschel ATLAS (H-ATLAS) will be the largest submillimetre survey to date, detecting close to half-a-million sources. It will only be possible to measure spectroscopic redshifts for a small fraction of these sources. However, if the rest-frame spectral energy distribution (SED) of a typical H-ATLAS source is known, this SED and the observed Herschel fluxes can be used to estimate the redshifts of the H-ATLAS sources without spectroscopic redshifts. In this paper, we use a subset of 40 H-ATLAS sources with previously measured redshifts in the range 0.5<z<4.2 to derive a suitable average template for high redshift H-ATLAS sources. We find that a template with two dust components T_c = 23.9 K, T_h = 46.9 K and ratio of mass of cold dust to mass of warm dust of 30.1) provides a good fit to the rest-frame fluxes of the sources in our calibration sample. We use a jackknife technique to estimate the accuracy of the redshifts estimated with this template, finding a root mean square of Delta z/(1+z) = 0.26. For sources for which there is prior information that they lie at z > 1 we estimate that the rms of Delta z/(1+z) = 0.12. We have used this template to estimate the redshift distribution for the sources detected in the H-ATLAS equatorial fields, finding a bimodal distribution with a mean redshift of 1.2, 1.9 and 2.5 for 250, 350 and 500 um selected sources respectively. \end{abstract}

We report the redshift of HATLAS J132427.0+284452 (hereafter HATLAS J132427), a gravitationally lensed starburst galaxy, the first determined 'blind' by the Herschel Space Observatory. This is achieved via the detection of [C II] consistent with z = 1.68 in a far-infrared spectrum taken with the SPIRE Fourier Transform Spectrometer. We demonstrate that the [C II] redshift is secure via detections of CO J = 2 - 1 and 3 - 2 using the Combined Array for Research in Millimeter-wave Astronomy and the Institut de Radioastronomie Millimetrique's Plateau de Bure Interferometer. The intrinsic properties appear typical of high-redshift starbursts despite the high lensing-amplified fluxes, proving the ability of the FTS to probe this population with the aid of lensing. The blind detection of [C II] demonstrates the potential of the SAFARI imaging spectrometer, proposed for the much more sensitive SPICA mission, to determine redshifts of multiple dusty galaxies simultaneously without the benefit of lensing.

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

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

Panchromatic observations of the best candidate HyLIRG from the widest Herschel extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z=2.41 galaxies across a ~100-kpc region - a cluster of starbursting proto-ellipticals. Via sub-arcsecond interferometric imaging we have measured accurate gas and star-formation surface densities. The two brightest galaxies span ~3 kpc FWHM in submm/radio continuum and CO J=4-3, and double that in CO J=1-0. The broad CO line is due partly to the multitude of constituent galaxies and partly to large rotational velocities in two counter-rotating gas disks -- a scenario predicted to lead to the most intense starbursts, which will therefore come in pairs. The disks have M(dyn) of several x 10^11 Msun, and gas fractions of ~40%. Velocity dispersions are modest so the disks are unstable, potentially on scales commensurate with their radii: these galaxies are undergoing extreme bursts of star formation, not confined to their nuclei, at close to the Eddington limit. Their specific star-formation rates place them ~>5x above the main sequence, which supposedly comprises large gas disks like these. Their high star-formation efficiencies are difficult to reconcile with a simple volumetric star-formation law. N-body and dark matter simulations suggest this system is the progenitor of a B(inary)-type ~10^14.6-Msun cluster.

We compare the environmental and star formation properties of far-infrared detected and non--far-infrared detected galaxies out to $z \sim0.5$. Using optical spectroscopy and photometry from the Galaxy And Mass Assembly (GAMA) and Sloan Digital Sky Survey (SDSS), with far-infrared observations from the {\em Herschel}-ATLAS Science Demonstration Phase (SDP), we apply the technique of Voronoi Tessellations to analyse the environmental densities of individual galaxies. Applying statistical analyses to colour, $r-$band magnitude and redshift-matched samples, we show there to be a significant difference at the 3.5$\sigma$ level between the normalized environmental densities of these two populations. This is such that infrared emission (a tracer of star formation activity) favours underdense regions compared to those inhabited by exclusively optically observed galaxies selected to be of the same $r-$band magnitude, colour and redshift. Thus more highly star-forming galaxies are found to reside in the most underdense environments, confirming previous studies that have proposed such a correlation. However, the degeneracy between redshift and far-infrared luminosity in our flux-density limited sample means that we are unable to make a stronger statement in this respect. We then apply our method to synthetic light cones generated from semi-analytic models, finding that over the whole redshift distribution the same correlations between star-formation rate and environmental density are found.

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.

We present a detailed study of how the Star Formation Rate (SFR) relates to the interstellar medium (ISM) of M31 at ~140pc scales. The SFR is calculated using the far-ultraviolet and 24um emission, corrected for the old stellar population in M31. We find a global value for the SFR of 0.25+/-0.05Msun/yr and compare this with the SFR found using the total far-infrared (FIR) luminosity. There is general agreement in regions where young stars dominate the dust heating. Atomic hydrogen (HI) and molecular gas (traced by carbon monoxide, CO) or the dust mass is used to trace the total gas in the ISM. We show that the global surface densities of SFR and gas mass place M31 amongst a set of low-SFR galaxies in the plot of Kennicutt (1998b). The relationship between SFR and gas surface density is tested in six radial annuli across M31, assuming a power law relationship with index, N. The star formation law using total gas traced by HI and CO gives a global index of N=2.03+/-0.04, with a significant variation with radius; the highest values are observed in the 10kpc ring. We suggest that this slope is due to HI turning molecular at ~10Msun/pc2. When looking at H2 regions, we measure a higher mean SFR suggesting a better spatial correlation between H2 and SF. We find N~0.6 with consistent results throughout the disk - this is at the low end of values found in previous work and argues against a superlinear SF law on small scales.

We perform a stacking analysis of H-ATLAS data in order to obtain isothermal dust temperatures and rest-frame luminosities at 250um (L250), for 1599 radio sources over the H-ATLAS P1 GAMA area. The radio sample is generated using a combination of NVSS data and K-band UKIDSS-LAS data, over 0.01<z<0.8. The FIR properties of the sample are investigated as a function of L_1.4 GHz, z, projected radio-source size and radio spectral index. In order to search for stellar-mass dependent relations, we split the parent sample into those sources which are below and above 1.5 L_{K}^{*}. Correcting for stellar mass and redshift, we find no relation between the L250 and L_1.4GHz of radio AGN. This implies that a galaxy's nominal radio luminosity has little or no bearing on the SFR and/or dust mass content of the host system. The L250 of both the radio detected and non radio-detected galaxies (defined as those sources not detected at 1.4GHz but detected in the SDSS with r'<22) rises with increasing z. For mass and colour matched samples, sub-1.5 L_{K}^{*} and super-1.5 L_{K}^{*} radio-detected galaxies have 0.89\pm0.18 and 0.49\pm0.12 times the L250 of their non radio-detected counterparts. We explain these results in terms of the hotter, denser and richer halo environments massive radio-galaxies maintain and are embedded in. These results indicate that all massive radio galaxies (>1.5 L_{K}^{*}) may have systematically lower FIR luminosities (~25%) than their colour-matched non radio-detected counterparts. Compact radio sources (<30kpc) are associated with higher L250s and dust temperatures than their more extended (>30kpc) counterparts. The higher dust temperature suggests that this may be attributed to enhanced SFRs, but whether this is directly or indirectly due to radio activity (e.g. jet induced or merger-driven SF) is as yet unknown.

[abridged] Aims. The aim of the present paper is to provide new and more detailed relations at the kpc scale between the gas surface density and the face-on optical depth directly calibrated on galaxies, in order to compute the attenuation not only for semi-analytic models but also observationally as new and upcoming radio observatories are able to trace gas ever farther in the Universe. Methods. We have selected a sample of 4 nearby resolved galaxies and a sample of 27 unresolved galaxies from the Herschel Reference Survey and the Very Nearby Galaxies Survey, for which we have a large set of multi-wavelength data from the FUV to the FIR including metallicity gradients for resolved galaxies, along with radio HI and CO observations. For each pixel in resolved galaxies and for each galaxy in the unresolved sample, we compute the face-on optical depth from the attenuation determined with the CIGALE SED fitting code and an assumed geometry. We determine the gas surface density from HI and CO observations with a metallicity-dependent XCO factor. Results. We provide new, simple to use, relations to determine the face-on optical depth from the gas surface density, taking the metallicity into account, which proves to be crucial for a proper estimate. The method used to determine the gas surface density or the face-on optical depth has little impact on the relations except for galaxies that have an inclination over 50d. Finally, we provide detailed instructions on how to compute the attenuation practically from the gas surface density taking into account possible information on the metallicity. Conclusions. Examination of the influence of these new relations on simulated FUV and IR luminosity functions shows a clear impact compared to older oft-used relations, which in turn could affect the conclusions drawn from studies based on large scale cosmological simulations.

We present two large, nearby (0.013$\le$z$\le$0.06) samples of Early-Type Galaxies (ETGs): a visually classified sample of 220 ETGs, created using source-matched data from the Galaxy and Mass Assembly (GAMA) database with FIR/sub-mm detections from $Herschel$-ATLAS; and a visually classified sample of 551 ETGs which are undetected with $Herschel$-ATLAS. Active galactic nuclei (AGN) are removed from our samples using optical emission line diagnostics. These samples are scrutinised to determine characteristics of sub-mm detected versus undetected ETGs. We find similarities in the stellar mass distributions of the two ETG samples but testing other properties uncovers significant differences. The sub-mm detected sample is shown to have lower concentration and S\'ersic indices than those in the undetected sample - a result which may be linked to the presence of dust in the former. Optical and UV-optical colours are also shown to be much bluer, indicating that the dust is linked with recent star formation. The intrinsic effective radii are on average 1.5 times larger for the sub-mm detected ETGs. Surface densities and groups data from the GAMA database are examined for the two samples, leading to the conclusion that dusty ETGs inhabit sparser environments than non-dusty ETGs in the nearby universe, although environments of the brightest ETGs are shown to differ the least. Modified Planck functions are fit to the H-ATLAS detected PACS and SPIRE fluxes for ETGs with sub-mm flux densities of at least 3$\sigma$ in the 350$\mu$m SPIRE band, giving a resultant mean cold dust temperature of T$_{d}$=22.1K, with a range of 9-30K. The corresponding mean dust mass is 1.8$\times10^{7}$M$_{\odot}$, with a range of (0.08-35.0)$\times10^{7}$M$_{\odot}$. The dust masses calculated from these fits, normalised by stellar mass, are shown to increase with decreasing stellar mass and bluer colours.

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.

Using IRAM PdBI we report the detection of H2O in six new lensed ultra-luminous starburst galaxies at high redshift, discovered in the Herschel H-ATLAS survey. The sources are detected either in the 2_{02}-1_{11} or 2_{11}-2_{02} H_2O emission lines with integrated line fluxes ranging from 1.8 to 14 Jy.km/s. The corresponding apparent luminosities are mu x L_H2O ~ 3-12 x 10^8 Lo, where mu is the lensing magnification factor (3 < mu < 12). These results confirm that H2O lines are among the strongest molecular lines in such galaxies, with intensities almost comparable to those of the high-J CO lines, and same profiles and line widths (200-900 km/s) as the latter. With the current sensitivity of PdBI, H2O can therefore easily be detected in high-z lensed galaxies (with F(500um) > 100 mJy) discovered in the Herschel surveys. Correcting the luminosities for lensing amplification, L_H2O is found to have a strong dependence on the IR luminosity, varying as ~L_IR^{1.2}. This relation which needs to be confirmed with better statistics, may indicate a role of radiative (IR) excitation of the H2O lines, and implies that high-z galaxies with L_IR >~ 10^13 Lo tend to be very strong emitters in H2O, that have no equivalent in the local universe.

Dust emission at sub-millimetre wavelengths allows us to trace the early phases of star formation in the Universe. In order to understand the physical processes involved in this mode of star formation, it is essential to gain knowledge about the dark matter structures - most importantly their masses - that sub-millimetre galaxies live in. Here we use the magnification effect of gravitational lensing to determine the average mass and dust content of sub-millimetre galaxies with 250mu flux densities of S_250>15mJy selected using data from the Herschel Multi-tiered Extragalactic Survey. The positions of hundreds of sub-millimetre foreground lenses are cross-correlated with the positions of background Lyman-break galaxies at z~3-5 selected using optical data from the Canada-France Hawaii Telescope Legacy Survey. We detect a cross-correlation signal at the 7-sigma level over a sky area of one square degree, with ~80% of this signal being due to magnification, whereas the remaining ~20% comes from dust extinction. Adopting some simple assumptions for the dark matter and dust profiles and the redshift distribution enables us to estimate the average mass of the halos hosting the sub-millimetre galaxies to be log(M_200/M_sun)=13.17+0.05-0.08(stat.) and their average dust mass fraction (at radii of >10kpc) to be M_dust/M_200~6x10^-5. This supports the picture that sub-millimetre galaxies are dusty, forming stars at a high rate, reside in massive group-sized halos, and are a crucial phase in the assembly and evolution of structure in the Universe.

We present a measurement of the angular power spectrum of the cosmic far-infrared background (CFIRB) anisotropies in one of the extragalactic fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) at 250, 350 and 500 \mu m bands. Consistent with recent measurements of the CFIRB power spectrum in Herschel-SPIRE maps, we confirm the existence of a clear one-halo term of galaxy clustering on arcminute angular scales with large-scale two-halo term of clustering at 30 arcminutes to angular scales of a few degrees. The power spectrum at the largest angular scales, especially at 250 \mu m, is contaminated by the Galactic cirrus. The angular power spectrum is modeled using a conditional luminosity function approach to describe the spatial distribution of unresolved galaxies that make up the bulk of the CFIRB. Integrating over the dusty galaxy population responsible for the background anisotropies, we find that the cosmic abundance of dust, relative to the critical density, to be between \Omega_dust=10^{-6} and 8 x 10^{-6} in the redshift range z ~ 0-3. This dust abundance is consistent with estimates of the dust content in the Universe using quasar reddening and magnification measurements in the SDSS.

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

We have constructed a sample of radio-loud objects with optical spectroscopy from the Galaxy and Mass Assembly (GAMA) project over the Herschel-ATLAS Phase 1 fields. Classifying the radio sources in terms of their optical spectra, we find that strong-emission-line sources (`high-excitation radio galaxies') have, on average, a factor ~4 higher 250-micron Herschel luminosity than weak-line (`low-excitation') radio galaxies and are also more luminous than magnitude-matched radio-quiet galaxies at the same redshift. Using all five H-ATLAS bands, we show that this difference in luminosity between the emission-line classes arises mostly from a difference in the average dust temperature; strong-emission-line sources tend to have comparable dust masses to, but higher dust temperatures than, radio galaxies with weak emission lines. We interpret this as showing that radio galaxies with strong nuclear emission lines are much more likely to be associated with star formation in their host galaxy, although there is certainly not a one-to-one relationship between star formation and strong-line AGN activity. The strong-line sources are estimated to have star-formation rates at least a factor 3-4 higher than those in the weak-line objects. Our conclusion is consistent with earlier work, generally carried out using much smaller samples, and reinforces the general picture of high-excitation radio galaxies as being located in lower-mass, less evolved host galaxies than their low-excitation counterparts.

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

We have measured the clustering properties of low-redshift (z < 0.3) sub-mm galaxies detected at 250 micron in the Herschel-ATLAS Science Demonstration Phase (SDP) field. We selected a sample for which we have high-quality spectroscopic redshifts, obtained from reliably matching the 250-micron sources to a complete (for r < 19.4) sample of galaxies from the GAMA database. Both the angular and spatial clustering strength are measured for all z < 0.3 sources as well as for five redshift slices with thickness delta z=0.05 in the range 0.05 < z < 0.3. Our measured spatial clustering length r_0 is comparable to that of optically-selected, moderately star-forming (blue) galaxies: we find values around 5 Mpc. One of the redshift bins contains an interesting structure, at z = 0.164.

We combine new dust continuum observations of the edge-on spiral galaxy NGC 4565 in all Herschel/SPIRE (250, 350, 500 micron) wavebands, obtained as part of the Herschel Reference Survey, and a large set of ancillary data (Spitzer, SDSS, GALEX) to analyze its dust energy balance. We fit a radiative transfer model for the stars and dust to the optical maps with the fitting algorithm FitSKIRT. To account for the observed UV and mid-infrared emission, this initial model was supplemented with both obscured and unobscured star-forming regions. Even though these star-forming complexes provide an additional heating source for the dust, the far-infrared/submillimeter emission long wards of 100 micron is underestimated by a factor of 3-4. This inconsistency in the dust energy budget of NGC 4565 suggests that a sizable fraction (two-thirds) of the total dust reservoir (Mdust ~ 2.9e+8 Msun) consists of a clumpy distribution with no associated young stellar sources. The distribution of those dense dust clouds would be in such a way that they remain unresolved in current far-infrared/submillimeter observations and hardly comtribute to the attenuation at optical wavelengths. More than two-thirds of the dust heating in NGC 4565 is powered by the old stellar population, with localized embedded sources supplying the remaining dust heating in NGC 4565. The results from this detailed dust energy balance study in NGC 4565 is consistent with that of similar analyses of other edge-on spirals.

We have used data from the Herschel-ATLAS at 250, 350 and 500 \mu m to determine the far-infrared (FIR) properties of 50 Broad Absorption Line Quasars (BAL QSOs). Our sample contains 49 high-ionization BAL QSOs (HiBALs) and 1 low-ionization BAL QSO (LoBAL) which are compared against a sample of 329 non-BAL QSOs. These samples are matched over the redshift range 1.5 \leq z < 2.3 and in absolute i-band magnitude over the range -28 \leq M_{i} \leq -24. Of these, 3 BAL QSOs (HiBALs) and 27 non-BAL QSOs are detected at the > 5 sigma level. We calculate star-formation rates (SFR) for our individually detected HiBAL QSOs and the non-detected LoBAL QSO as well as average SFRs for the BAL and non-BAL QSO samples based on stacking the Herschel data. We find no difference between the HiBAL and non-BAL QSO samples in the FIR, even when separated based on differing BAL QSO classifications. Using Mrk 231 as a template, the weighted mean SFR is estimated to be \approx240\pm21 M_{\odot} yr^{-1} for the full sample, although this figure should be treated as an upper limit if AGN-heated dust makes a contribution to the FIR emission. Despite tentative claims in the literature, we do not find a dependence of {\sc C\,iv} equivalent width on FIR emission, suggesting that the strength of any outflow in these objects is not linked to their FIR output. These results strongly suggest that BAL QSOs (more specifically HiBALs) can be accommodated within a simple AGN unified scheme in which our line-of-sight to the nucleus intersects outflowing material. Models in which HiBALs are caught towards the end of a period of enhanced spheroid and black-hole growth, during which a wind terminates the star-formation activity, are not supported by the observed FIR properties.

We have obtained Herschel images at five wavelengths from 100 to 500 micron of a ~5.5x2.5 degree area centred on the local galaxy M31 (Andromeda), our nearest neighbour spiral galaxy, as part of the Herschel guaranteed time project "HELGA". The main goals of HELGA are to study the characteristics of the extended dust emission, focusing on larger scales than studied in previous observations of Andromeda at an increased spatial resolution, and the obscured star formation. In this paper we present data reduction and Herschel maps, and provide a description of the far-infrared morphology, comparing it with features seen at other wavelengths. We use high--resolution maps of the atomic hydrogen, fully covering our fields, to identify dust emission features that can be associated to M31 with confidence, distinguishing them from emission coming from the foreground Galactic cirrus. Thanks to the very large extension of our maps we detect, for the first time at far-infrared wavelengths, three arc-like structures extending out to ~21, ~26 and ~31 kpc respectively, in the south-western part of M31. The presence of these features, hosting ~2.2e6 Msol of dust, is safely confirmed by their detection in HI maps. Overall, we estimate a total dust mass of ~5.8e7 Msol, about 78% of which is contained in the two main ring-like structures at 10 and 15 kpc, at an average temperature of 16.5 K. We find that the gas-to-dust ratio declines exponentially as a function of the galacto-centric distance, in agreement with the known metallicity gradient, with values ranging from 66 in the nucleus to ~275 in the outermost region. [Abridged]

We present new observations from Z-Spec, a broadband 185-305 GHz spectrometer, of five sub-millimeter bright lensed sources selected from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration phase (SDP) catalog. We construct a redshift finding algorithm using combinations of the signal-to-noise of all the lines falling in the Z-Spec bandpass to determine redshifts with high confidence, even in cases where the signal-to-noise in individual lines is low. We measure the dust continuum in all sources and secure CO redshifts for four out of five (z~1.5-3). In one source, SDP.17, we tentatively identify two independent redshifts and a water line, confirmed at z=2.308. Our sources have properties characteristic of dusty starburst galaxies, with magnification-corrected star formation rates of 10^(2-3) M_sun/yr. Lower limits for the dust masses (~a few 10^8 M_sun) and spatial extents (~1 kpc equivalent radius) are derived from the continuum spectral energy distributions, corresponding to dust temperatures between 54 and 69K. In the LTE approximation, we derive relatively low CO excitation temperatures (< 100 K) and optical depths (tau<1). Performing a non-LTE excitation analysis using RADEX, we find that the CO lines measured by Z-Spec (from J=4->3 to 10->9, depending on the galaxy) localize the best solutions to either a high-temperature / low-density region, or a low-temperature / high-density region near the LTE solution, with the optical depth varying accordingly. Observations of additional CO lines, CO(1-0) in particular, are needed to constrain the non-LTE models.

We present a pan-chromatic analysis of an unprecedented sample of 1402 250 micron-selected galaxies at z < 0.5 (mean z = 0.24) from the Herschel-ATLAS survey. We complement our Herschel 100-500 micron data with UV-K-band photometry from the Galaxy And Mass Assembly (GAMA) survey and apply the MAGPHYS energy-balance technique to produce pan-chromatic SEDs for a representative sample of 250 micron selected galaxies spanning the most recent 5 Gyr of cosmic history. We derive estimates of physical parameters, including star formation rates, stellar masses, dust masses and infrared luminosities. The typical H-ATLAS galaxy at z < 0.5 has a far-infrared luminosity in the range 10^10 - 10^12 Lsolar (SFR: 1-50 Msolar/yr) thus is broadly representative of normal star forming galaxies over this redshift range. We show that 250 micron-selected galaxies contain a larger mass of dust at a given infra-red luminosity or star formation rate than previous samples selected at 60 micron from IRAS. We derive typical SEDs for H-ATLAS galaxies, and show that the emergent SED shape is most sensitive to specific star formation rate. The optical-UV SEDs also become more reddened due to dust at higher redshifts. Our template SEDs are significantly cooler than existing infra-red templates. They may therefore be most appropriate for inferring total IR luminosities from moderate redshift submillimetre selected samples and for inclusion in models of the lower redshift submillimetre galaxy populations.

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

We present an analysis of the dust and gas in Andromeda, using Herschel images sampling the entire far-infrared peak. We fit a modified-blackbody model to ~4000 quasi-independent pixels with spatial resolution of ~140pc and find that a variable dust-emissivity index (beta) is required to fit the data. We find no significant long-wavelength excess above this model suggesting there is no cold dust component. We show that the gas-to-dust ratio varies radially, increasing from ~20 in the center to ~70 in the star-forming ring at 10kpc, consistent with the metallicity gradient. In the 10kpc ring the average beta is ~1.9, in good agreement with values determined for the Milky Way (MW). However, in contrast to the MW, we find significant radial variations in beta, which increases from 1.9 at 10kpc to ~2.5 at a radius of 3.1kpc and then decreases to 1.7 in the center. The dust temperature is fairly constant in the 10kpc ring (ranging from 17-20K), but increases strongly in the bulge to ~30K. Within 3.1kpc we find the dust temperature is highly correlated with the 3.6 micron flux, suggesting the general stellar population in the bulge is the dominant source of dust heating there. At larger radii, there is a weak correlation between the star formation rate and dust temperature. We find no evidence for 'dark gas' in M31 in contrast to recent results for the MW. Finally, we obtained an estimate of the CO X-factor by minimising the dispersion in the gas-to-dust ratio, obtaining a value of (1.9+/-0.4)x10^20 cm^-2 [K kms^-1]^-1.

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

The identification of strong gravitational lenses in large surveys has historically been a rather time consuming exercise. Early data from the Herschel Astrophysical Terahertz Large Area Survey (Herschel-ATLAS) demonstrate that lenses can be identified efficiently at submillimetre wavelengths using a simple flux criteria. Motivated by that development, this work considers the statistical properties of strong gravitational lens systems which have been, and will be, found by the Herschel-ATLAS. Analytical models of lens statistics are tested with the current best estimates for the various model ingredients. These include the cosmological parameters, the mass function and the lens density profile, for which we consider the singular isothermal sphere (SIS) and the Navarro, Frenk & White (NFW) approximations. The five lenses identified in the Herschel-ATLAS Science Demonstration Phase suggest a SIS density profile is preferred, but cannot yet constrain \Omega_\Lambda to an accuracy comparable with other methods. The complete Herschel-ATLAS data set should be sufficient to yield competitive constraints on \Omega_\Lambda. Whilst this huge number of lenses has great potential for constraining cosmological parameters, they will be most powerful in constraining uncertainty in astrophysical processes. Further investigation is needed to fully exploit this unprecedented data set.

Using infrared imaging from the Herschel Space Observatory, observed as part of the VNGS, we investigate the spatially resolved dust properties of the interacting Whirlpool galaxy system (NGC 5194 and NGC 5195), on physical scales of 1 kpc. Spectral energy distribution modelling of the new infrared images in combination with archival optical, near- through mid-infrared images confirms that both galaxies underwent a burst of star formation ~370-480 Myr ago and provides spatially resolved maps of the stellar and dust mass surface densities. The resulting average dust-to-stellar mass ratios are comparable to other spiral and spheroidal galaxies studied with Herschel, with NGC 5194 at log M(dust)/M(star)= -2.5+/-0.2 and NGC 5195 at log M(dust)/M(star)= -3.5+/-0.3. The dust-to-stellar mass ratio is constant across NGC 5194 suggesting the stellar and dust components are coupled. In contrast, the mass ratio increases with radius in NGC 5195 with decreasing stellar mass density. Archival mass surface density maps of the neutral and molecular hydrogen gas are also folded into our analysis. The gas-to-dust mass ratio, 94+/-17, is relatively constant across NGC 5194. Somewhat surprisingly, we find the dust in NGC 5195 is heated by a strong interstellar radiation field, over 20 times that of the ISRF in the Milky Way, resulting in relatively high characteristic dust temperatures (~30 K). This post-starburst galaxy contains a substantial amount of low-density molecular gas and displays a gas-to-dust ratio (73+/-35) similar to spiral galaxies. It is unclear why the dust in NGC 5195 is heated to such high temperatures as there is no star formation in the galaxy and its active galactic nucleus is 5-10 times less luminous than the one in NGC 5194, which exhibits only a modest enhancement in the amplitude of its ISRF.

The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions in the galaxies in this sample and explore the correlation of the far-infrared luminosity, which traces star formation, with the CO luminosity, which traces the molecular gas mass. By comparing the NGLS data with merging galaxies at low and high redshift which have also been observed in the CO J=3-2 line, we show that the correlation of far-infrared and CO luminosity shows a significant trend with luminosity. This trend is consistent with a molecular gas depletion time which is more than an order of magnitude faster in the merger galaxies than in nearby normal galaxies. We also find a strong correlation of the L(FIR)/L(CO3-2) ratio with the atomic to molecular gas mass ratio. This correlation suggests that some of the far-infrared emission originates from dust associated with atomic gas and that its contribution is particularly important in galaxies where most of the gas is in the atomic phase.

We present GALEX far-ultraviolet (FUV) and near-ultraviolet (NUV) as well as SDSS g, r, i photometry and structural parameters for the Herschel Reference Survey, a magnitude-, volume-limited sample of nearby galaxies in different environments. We use this unique dataset to investigate the ultraviolet (UV) structural scaling relations of nearby galaxies and to determine how the properties of the UV disk vary with atomic hydrogen content and environment. We find a clear change of slope in the stellar mass vs. effective surface brightness relation when moving from the optical to the UV, with more massive galaxies having brighter optical but fainter UV surface brightnesses than smaller systems. A similar change of slope is also seen in the radius vs. surface brightness relation. By comparing our observations with the predictions of a simple multi-zone chemical model of galaxy evolution, we show that these findings are a natural consequence of a much more efficient inside-out growth of the stellar disk in massive galaxies. We confirm that isophotal radii are always a better proxy for the size of the stellar/star-forming disk than effective quantities and we show that the extent of the UV disk (normalized to the optical size) is strongly correlated to the integrated HI gas fraction. This relation still holds even when cluster spirals are considered, with HI-deficient systems having less extended star-forming disks than HI-normal galaxies. Interestingly, the star formation in the inner part of HI-deficient galaxies is significantly less affected by the removal of the atomic hydrogen, as expected in a simple ram-pressure stripping scenario. These results suggest that it is the amount of HI that regulates the growth of the star-forming disk in the outskirts of galaxies.

We identify near-infrared Ks band counterparts to Herschel-ATLAS sub-mm sources, using a preliminary object catalogue from the VISTA VIKING survey. The sub-mm sources are selected from the H-ATLAS Phase 1 catalogue of the GAMA 9h field, which includes all objects detected at 250, 350 or 500 um with the SPIRE instrument. We apply and discuss a likelihood ratio (LR) method for VIKING candidates within a search radius of 10" of the 22,000 SPIRE sources with a 5 sigma detection at 250 um. We find that 11,294(51%) of the SPIRE sources have a best VIKING counterpart with a reliability $R\ge 0.8$, and the false identification rate of these is estimated to be 4.2%. We expect to miss ~5% of true VIKING counterparts. There is evidence from Z-J and J-Ks colours that the reliable counterparts to SPIRE galaxies are marginally redder than the field population. We obtain photometric redshifts for ~68% of all (non-stellar) VIKING candidates with a median redshift of 0.405. Comparing to the results of the optical identifications supplied with the Phase I catalogue, we find that the use of medium-deep near-infrared data improves the identification rate of reliable counterparts from 36% to 51%.

The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) provides an unprecedented opportunity to search for blazars at sub-mm wavelengths. We cross-matched the FIRST radio source catalogue with the 11655 sources brighter than 35 mJy at 500{\mu}m in the \sim 135 square degrees of the sky covered by the H-ATLAS equatorial fields at 9 h and 15 h, plus half of the field at 12 h. We found that 379 of the H-ATLAS sources have a FIRST counterpart within 10 arcsec, including 8 catalogued blazars (plus one known blazar that was found at the edge of one the H-ATLAS maps). To search for additional blazar candidates we have devised new diagnostic diagrams and found that known blazars occupy a region of the log(S500{\mu}m/S350{\mu}m) vs. log(S500{\mu}m/S1.4GHz) plane separated from that of the other sub-mm sources with radio counterparts. Using this diagnostic we have selected 12 further candidates that turn out to be scattered in the (r-z) vs. (u-r) plane or in the WISE colour-colour diagram proposed by Massaro et al. (2012), where known blazars are concentrated in well defined strips. This suggests that the majority of them either are not blazars or have spectral energy distributions contaminated by their host galaxies. A significant fraction of true blazars are found to be hosted by star-forming galaxies. This finding, supported by an analysis of blazars detected in Planck 545 and 857 GHz bands, is at odds with the notion that blazar hosts are passive ellipticals and indicates that the sub-mm selection is providing a novel prospect on blazar properties. Based on an inspection of the available photometric data, including the WISE all-sky survey, the unpublished VIKING survey and new radio observations, we tentatively estimate that there are 11 blazars with synchrotron flux density S500{\mu}m > 35mJy over the considered area. This result already allows us to constrain blazar evolution models.

We present high-resolution maps of stars, dust, and molecular gas in a strongly lensed submillimeter galaxy (SMG) at z = 3.259. HATLAS12--00 is selected from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) as a strong lens candidate mainly based on its unusually high 500um flux density (~300 mJy). It is the only high-redshift Planck detection in the 130 deg^2 H-ATLAS Phase 1 area. Keck Adaptive Optics images reveal a quadruply imaged galaxy in the K-band while the Submillimeter Array and the Extended Very Large Array show doubly imaged 880um and CO(1-0) sources, indicating differentiated distributions of the various components in the galaxy. In the source plane, the stars reside in three major kpc-scale clumps extended over ~1.6 kpc, the dust in a compact (~1 kpc) region ~3 kpc north of the stars, and the cold molecular gas in an extended (~7 kpc) disk ~5 kpc northeast of the stars. The emission from the stars, dust, and gas are magnified by ~17, 8, and 7 times, respectively, by four lensing galaxies at z ~ 1. Intrinsically, the lensed galaxy is a warm (T_dust ~ 40-65 K), hyper-luminous (L_IR ~ 1.7e13 Lsun; SFR ~ 2000 Msun/yr), gas-rich (M_gas/M_baryon ~ 70%), young (M_stellar/SFR ~ 20 Myr), and short-lived (M_gas/SFR ~ 40 Myr) starburst, without a significant active galactic nucleus. With physical properties similar to unlensed z > 2 SMGs, HATLAS12--00 offers a detailed view of a typical SMG through a powerful cosmic microscope.

We report measurements of the carbon monoxide ground state rotational transition (12C16O J = 1--0) with the Zpectrometer ultra-wideband spectrometer on the 100-m diameter Green Bank Telescope. The sample comprises 11 galaxies with redshifts between z = 2.1 and 3.5 from a total sample of 24 targets identified by Herschel-ATLAS photometric colors from the SPIRE instrument. Nine of the CO measurements are new redshift determinations, substantially adding to the number of detections of galaxies with rest-frame peak submillimeter emission near 100um. The CO detections confirm the existence of massive gas reservoirs within these luminous dusty star-forming galaxies (DSFGs). The CO redshift distribution of the 350um-selected galaxies is strikingly similar to the optical redshifts of 850um-selected submillimeter galaxies (SMGs) in 2.1 < z < 3.5. Spectroscopic redshifts break a temperature-redshift degeneracy; optically thin dust models fit to the far-infrared photometry indicate characteristic dust temperatures near 34 K for most of the galaxies we detect in CO. Detections of two warmer galaxies and statistically significant nondetections hint at warmer or molecule-poor DSFGs with redshifts difficult determine from from Herschel-SPIRE photometric colors alone. Many of the galaxies identified by H-ATLAS photometry are expected to be amplified by foreground gravitational lenses. Analysis of CO linewidths and luminosities provides a method for finding approximate gravitational lens magnifications mu from spectroscopic data alone, yielding mu ~ 3--20. Corrected for magnification, most galaxy luminosities are consistent with an ultra-luminous infrared galaxy (ULIRG) classification, but three are candidate hyper-LIRGs with luminosities greater than 10^13 L_sun.

The Herschel Reference Survey (HRS) is a guaranteed time Herschel key project aimed at studying the physical properties of the interstellar medium in galaxies of the nearby universe. This volume limited, K-band selected sample is composed of galaxies spanning the whole range of morphological types (from ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster). We present flux density measurements of the whole sample of 323 galaxies of the HRS in the three bands of the Spectral and Photometric Imaging Receiver (SPIRE), at 250, 350 and 500 microns. Aperture photometry is performed on extended galaxies and point spread function (PSF) fitting on timeline data for unresolved objects; we carefully estimate errors and upper limits. The flux densities are found to be in good agreement with those of the HeViCS and KINGFISH key projects in all SPIRE bands, and of the Planck consortium at 350 and 550 microns, for the galaxies in common. This submillimetre catalogue of nearby galaxies is a benchmark for the study of the dust properties in the local universe, giving the zero redshift reference for any cosmological survey.

We present the results of a cross-correlation of the Planck Early Release Compact Source Catalog (ERCSC) with the catalog of Herschel-ATLAS sources detected in the Phase 1 fields, covering 134.55 deg2. There are 28 ERCSC sources detected by Planck at 857 GHz in this area. As many as 16 of them are probably high Galactic latitude cirrus; 10 additional sources can be clearly identified as bright, low-z galaxies; one further source is resolved by Herschel as two relatively bright sources; and the last is resolved into an unusual condensation of low-flux, probably high-redshift point sources, around a strongly lensed Herschel-ATLAS source at z = 3.26. Our results demonstrate that the higher sensitivity and higher angular resolution H-ATLAS maps provide essential information for the interpretation of candidate sources extracted from Planck sub-mm maps.

We describe the infrared properties of sources detected over ~36 deg^2 of sky in the GAMA 15-hr equatorial field, using data from both the Herschel Astrophysical Terahertz Large-Area Survey (H-ATLAS) and Wide-field Infrared Survey (WISE). With 5-sigma point-source depths of 34 and 0.048 mJy at 250 micron and 3.4 micron, respectively, we are able to identify 50.6% of the H-ATLAS sources in the WISE survey, corresponding to a surface density of ~630 deg^{-2}. Approximately two-thirds of these sources have measured spectroscopic or optical/near-IR photometric redshifts of z<1. For sources with spectroscopic redshifts at z<0.3, we find a linear correlation between the infrared luminosity at 3.4 micron and that at 250 micron, with +-50% scatter over ~1.5 orders of magnitude in luminosity, ~10^9 - 10^{10.5} L_sun. By contrast, the matched sources without previously measured redshifts (r>~20.5) have 250-350 micron flux density ratios that suggest either high-redshift galaxies (z>~1.5) or optically faint low-redshift galaxies with unusually low temperatures (T<~20). Their small 3.4-250 micron flux ratios favor a high-redshift galaxy population, as only the most actively star-forming galaxies at low redshift (e.g., Arp 220) exhibit comparable flux density ratios. Furthermore, we find a relatively large AGN fraction (~30%) in a 12 micron flux-limited subsample of H-ATLAS sources, also consistent with there being a significant population of high-redshift sources in the no-redshift sample.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

We present photometry of the nearby galaxy NGC 5128 (Centaurus A) observed with the PACS and SPIRE instruments on board the Herschel Space Observatory, at 70, 160, 250, 350 and 500 {\mu}m, as well as new CO J = 3-2 observations taken with the HARP-B instrument on the JCMT. Using a single component modified blackbody, we model the dust spectral energy distribution within the disk of the galaxy using all five Herschel wavebands, and find dust temperatures of ~30 K towards the centre of the disk and a smoothly decreasing trend to ~20 K with increasing radius. We find a total dust mass of (1.59 \pm 0.05) \times 10^7 M\odot, and a total gas mass of (2.7 \pm 0.2) \times 10^9 M\odot. The average gas-to-dust mass ratio is 103 \pm 8 but we find an interesting increase in this ratio to approximately 275 toward the centre of Cen A. We discuss several possible physical processes that may be causing this effect, including dust sputtering, jet entrainment and systematic variables such as the XCO factor. Dust sputtering by X-rays originating in the AGN or the removal of dust by the jets are our most favoured explanations.

We combine new Herschel/SPIRE sub-millimeter observations with existing multiwavelength data to investigate the dust scaling relations of the Herschel Reference Survey, a magnitude-, volume-limited sample of ~300 nearby galaxies in different environments. We show that the dust-to-stellar mass ratio anti-correlates with stellar mass, stellar mass surface density and NUV-r colour across the whole range of parameters covered by our sample. Moreover, the dust-to-stellar mass ratio decreases significantly when moving from late- to early-type galaxies. These scaling relations are similar to those observed for the HI gas-fraction, supporting the idea that the cold dust is tightly coupled to the cold atomic gas component in the interstellar medium. We also find a weak increase of the dust-to-HI mass ratio with stellar mass and colour but no trend is seen with stellar mass surface density. By comparing galaxies in different environments we show that, although these scaling relations are followed by both cluster and field galaxies, HI-deficient systems have, at fixed stellar mass, stellar mass surface density and morphological type systematically lower dust-to-stellar mass and higher dust-to-HI mass ratios than HI-normal/field galaxies. This provides clear evidence that dust is removed from the star-forming disk of cluster galaxies but the effect of the environment is less strong than what is observed in the case of the HI disk. Such effects naturally arise if the dust disk is less extended than the HI and follows more closely the distribution of the molecular gas phase, i.e., if the dust-to-atomic gas ratio monotonically decreases with distance from the galactic center.

We use the Herschel-ATLAS survey to conduct the first large-scale statistical study of the submm properties of optically selected galaxies. Using ~80,000 r-band selected galaxies from 126 deg^2 of the GAMA survey, we stack into submm imaging at 250, 350 and 500{\mu}m to gain unprecedented statistics on the dust emission from galaxies at z < 0.35. We find that low redshift galaxies account for 5% of the cosmic 250{\mu}m background (4% at 350{\mu}m; 3% at 500{\mu}m), of which approximately 60% comes from 'blue' and 20% from 'red' galaxies (rest-frame g - r). We compare the dust properties of different galaxy populations by dividing the sample into bins of optical luminosity, stellar mass, colour and redshift. In blue galaxies we find that dust temperature and luminosity correlate strongly with stellar mass at a fixed redshift, but red galaxies do not follow these correlations and overall have lower luminosities and temperatures. We make reasonable assumptions to account for the contaminating flux from lensing by red sequence galaxies and conclude that galaxies with different optical colours have fundamentally different dust emission properties. Results indicate that while blue galaxies are more luminous than red galaxies due to higher temperatures, the dust masses of the two samples are relatively similar. Dust mass is shown to correlate with stellar mass, although the dust/stellar mass ratio is much higher for low stellar mass galaxies, consistent with the lowest mass galaxies having the highest specific star formation rates. We stack the 250{\mu}m/NUV luminosity ratio, finding results consistent with greater obscuration of star formation at lower stellar mass and higher redshift. Submm luminosities and dust masses of all galaxies are shown to evolve strongly with redshift, indicating a fall in the amount of obscured star formation in ordinary galaxies over the last four billion years.

We present Herschel observations of 62 Early-Type Galaxies (ETGs), including 39 galaxies morphologically classified as S0+S0a and 23 galaxies classified as ellipticals using SPIRE at 250, 350 and 500 microns (and PACS 100 and 160 microns for 19 sources) as part of the volume-limited Herschel Reference Survey. We detect dust emission in 24% of the ellipticals and 62% of the S0s. The mean temperature of the dust is 23.9+/-0.8 K, warmer than that found for late-type galaxies in the Virgo Cluster. Including the non-detections, the mean dust mass is log(Mdust) = 5.9+/-0.1 and 5.2+/-0.1 Msun for the S0s and elliptical galaxies respectively. The mean dust-to-stellar mass is log(Mdust/Mstar) = -4.4+/-0.1 (S0s) and -5.8+/-0.1 (ellipticals). Virtually all the galaxies lie close to the red sequence yet the large number of detections of cool dust, the gas-to-dust ratios and the ratios of far-infrared to radio emission all suggest that many ETGs contain a cool interstellar medium similar to that in late-type galaxies. The mean dust-to-stellar mass ratio for S0s is approximatly a factor of ten less than for early-type spirals and the sizes of the dust sources in the S0s are also much smaller. We show that the difference cannot be explained by either the different bulge-to-disk ratios or environmental effects such as ram-pressure stripping. The wide range in the dust-to-stellar mass ratio for ETGs and the lack of a correlation between dust mass and optical luminosity suggest that much of the dust in the ETGs detected by Herschel has been acquired as the result of gravitational interactions; these interactions are unlikely to have had a major effect on the stellar masses of the ETGs. The Herschel observations tentatively suggest that in the most massive ETGs, the mass of the interstellar medium is unconnected to the evolution of the stellar populations.

While the selection of strongly lensed galaxies with 500{\mu}m flux density S_500>100 mJy has proven to be rather straightforward (Negrello et al. 2010), for many applications it is important to analyze samples larger than the ones obtained when confining ourselves to such a bright limit. Moreover, only by probing to fainter flux densities is possible to exploit strong lensing to investigate the bulk of the high-z star-forming galaxy population. We describe HALOS (the Herschel-ATLAS Lensed Objects Selection), a method for efficiently selecting fainter candidate strongly lensed galaxies, reaching a surface density of ~1.5-2 deg^-2, i.e. a factor of about 4 to 6 higher than that at the 100 mJy flux limit. HALOS will allow the selection of up to ~1000 candidate strongly lensed galaxies (with amplifications \mu>2) over the full H-ATLAS survey area. Applying HALOS to the H-ATLAS Science Demonstration Phase field (~14.4 deg^2) we find 31 candidate strongly lensed galaxies, whose candidate lenses are identified in the VIKING near-infrared catalog. Using the available information on candidate sources and candidate lenses we tentatively estimate a ~72% purity of the sample. The redshift distribution of the candidate lensed sources is close to that reported for most previous surveys for lensed galaxies, while that of candidate lenses extends to substantially higher redshifts than found in the other surveys. The counts of candidate strongly lensed galaxies are also in good agreement with model predictions (Lapi et al. 2011). Even though a key ingredient of the method is the deep near-infrared VIKING photometry, we show that H-ATLAS data alone allow the selection of a similarly deep sample of candidate strongly lensed galaxies with an efficiency close to 50%; a slightly lower surface density (~1.45 deg^-2) can be reached with a ~70% efficiency.

It has often been suggested that an alternative to the standard CO/21-cm method for estimating the mass of the interstellar medium (ISM) in a galaxy might be to estimate the mass of the ISM from the continuum dust emission. In this paper, we investigate the potential of this technique using Herschel observations of ten galaxies in the Herschel Reference Survey and in the Herschel Virgo Cluster Survey. We show that the emission detected by Herschel is mostly from dust that has a temperature and emissivity index similar to that of dust in the local ISM in our galaxy, with the temperature generally increasing towards the centre of each galaxy. We calibrate the dust method using the CO and 21-cm observations to provide an independent estimate of the mass of hydrogen in each galaxy, solving the problem of the uncertain `X factor' for the molecular gas by minimizing the dispersion in the ratio of the masses estimated using the two methods. With the calibration for the dust method and the estimate of the X-factor produced in this way, the dispersion in the ratio of the two gas masses is 30%, which gives an upper limit on the fundamental accuracy of the dust method. The calibration we obtain for the dust method is very similar to an independent Herschel measurement for M31 and to the calibration for the Milky Way from Planck measurements.

UV and optical surveys are essential to gain insight into the processes driving galaxy formation and evolution. The rest-frame UV emission is key to measure the cosmic SFR. However, UV light is strongly reddened by dust. In starburst galaxies, the UV colour and the attenuation are linked, allowing to correct for dust extinction. Unfortunately, evidence has been accumulating that the relation between UV colour and attenuation is different for normal star-forming galaxies when compared to starburst galaxies. It is still not understood why star-forming galaxies deviate from the UV colour-attenuation relation of starburst galaxies. Previous work and models hint that the role of the shape of the attenuation curve and the age of stellar populations have an important role. In this paper we aim at understanding the fundamental reasons to explain this deviation. We have used the CIGALE SED fitting code to model the far UV to the far IR emission of a set of 7 reasonably face-on spiral galaxies from the HRS. We have explored the influence of a wide range of physical parameters to quantify their influence and impact on the accurate determination of the attenuation from the UV colour, and why normal galaxies do not follow the same relation as starburst galaxies. We have found that the deviation can be best explained by intrinsic UV colour differences between different regions in galaxies. Variations in the shape of the attenuation curve can also play a secondary role. Standard age estimators of the stellar populations prove to be poor predictors of the intrinsic UV colour. These results are also retrieved on a sample of 58 galaxies when considering their integrated fluxes. When correcting the emission of normal star-forming galaxies for the attenuation, it is crucial to take into account possible variations in the intrinsic UV colour as well as variations of the shape of the attenuation curve.

We study the far infrared (60-500 $\mu$m) colours of late-type galaxies in the $Herschel$ Reference Survey, a K-band selected, volume limited sample of nearby galaxies. The far infrared colours are correlated with each other, with tighter correlations for the indices that are closer in wavelength. We also compare the different colour indices to various tracers of the physical properties of the target galaxies, such as the surface brightness of the ionising and non-ionising stellar radiation, the dust attenuation and the metallicity. The emission properties of the cold dust dominating the far infrared spectral domain are regulated by the properties of the interstellar radiation field. Consistent with that observed in nearby, resolved galaxies, our analysis shows that the ionising and the non-ionising stellar radiation, including that emitted by the most evolved, cold stars, both contribute to the heating of the cold dust component. This work also shows that metallicity is another key parameter characterising the cold dust emission of normal, late-type galaxies. A single modified black body with a grain emissivity index $\beta$=1.5 better fits the observed SPIRE flux density ratios $S250/S350$ vs. $S350/S500$ than $\beta$=2, although values of $\beta$ $\simeq$ 2 are possible in metal rich, high surface brightness galaxies. Values of $\beta$ $\lesssim$ 1.5 better represent metal poor, low surface brightness objects. This observational evidence provides strong constraints for dust emission models of normal, late type galaxies.

We examine the dust and gas properties of the nearby, barred galaxy M83, which is part of the Very Nearby Galaxy Survey. Using images from the PACS and SPIRE instruments of Herschel, we examine the dust temperature and dust mass surface density distribution. We find that the nuclear, bar and spiral arm regions exhibit higher dust temperatures and masses compared to interarm regions. However, the distribution of dust temperature and mass are not spatially coincident. Assuming a trailing spiral structure, the dust temperature peaks in the spiral arms lie ahead of the dust surface density peaks. The dust mass surface density correlates well with the distribution of molecular gas as traced by CO (J=3-2) images (JCMT) and the star formation rate as traced by H?2 with a correction for obscured star formation using 24 micron emission. Using HI images from THINGS to trace the atomic gas component, we make total gas mass surface density maps and calculate the gas-to-dust ratio. We find a mean gas-to-dust ratio of 84 \pm 4 with higher values in the inner region assuming a constant CO-to-H2 conversion factor. We also examine the gas-to-dust ratio using CO-to-H2 conversion factor that varies with metallicity.

The FIR/submm window is one of the least-studied regions of the electromagnetic spectrum, yet this wavelength range is absolutely crucial for understanding the physical processes and properties of the ISM in galaxies. The advent of the Herschel Space Observatory has opened up the entire FIR/submm window for spectroscopic studies. We present the first FIR/submm spectroscopic results on both nearby and distant galaxies obtained in the frame of two Herschel key programs: the Very Nearby Galaxies Survey and the Herschel ATLAS.

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

We present the dust properties and star-formation histories of local submillimetre-selected galaxies in Herschel-ATLAS, classified by optical morphology. The early-type galaxies (ETGs) that are detected contain as much dust as typical spirals, and form a unique sample that has been blindly selected at submillimetre wavelengths. Comparing H-ATLAS galaxies to a control sample of optically selected galaxies, we find 5.5% of luminous ETGs are detected in H-ATLAS. The H-ATLAS ETGs contain a significant mass of cold dust: the mean dust mass is 5.5x10^7 Msun, with individual galaxies ranging from 9x10^5-4x10^8 Msun. This is comparable to that of spirals in our sample, and is an order of magnitude more dust than that found for the control ETGs, which have a median dust mass inferred from stacking of (0.8-4.0)x10^6 Msun. The ETGs detected in H-ATLAS have bluer NUV-r colours, higher specific star-formation rates and younger stellar populations than ETGs which are optically selected, and may be transitioning from the blue cloud to the red sequence. We also find that H-ATLAS and control ETGs inhabit similar low-density environments. We conclude that the dust in H-ATLAS and control ETGs cannot be solely from stellar sources, and a large contribution from dust formed in the ISM or external sources is required. Alternatively, dust destruction may not be as efficient as predicted. We also explore the properties of the most passive spiral galaxies in our sample with SSFR<10^-11/yr. We find these passive spirals have lower dust-to-stellar mass ratios, higher stellar masses and older stellar population ages than normal spirals. The passive spirals inhabit low density environments similar to those of the normal spiral galaxies in our sample. This shows that the processes which turn spirals passive do not occur solely in the intermediate density environments of group and cluster outskirts. (Abridged)

We present the first results of a multi-wavelength survey, incoporating Herschel-SPIRE, Spitzer, GALEX and ATCA observations, of a 1 deg x 1 deg field centred on Centaurus A. As well as detecting the inner lobes of the active galactic nucleus (AGN) jet and counterjet, we have found two clouds, bright at sub-mm wavelengths, ~15 kpc from the centre of Cen A that are co-aligned with the jets. Flux measurements at Herschel wavelengths have proved vital in constraining fits to the Spectral Energy Distributions (SEDs). The clouds are well fit by a single-temperature, modified blackbody spectrum (beta=2) indicating that we are looking at two cold dust clouds on the outskirts of Cen A. The temperature and masses of the clouds are: T_{north} = 12.6^{+1.1}_{-1.2} K, T_{south} = 15.1^{+1.7}_{-1.6} K; log(M_{north} / M_o) = 5.8^{+0.2}_{-0.2}, log(M_{south} / M_o) = 5.6^{+0.2}_{-0.2} and the gas-dust ratio for both clouds is ~100. The measured values for the northern dust cloud are consistent with previous measurements from ISO while the southern cloud is a new sub-mm detection. The two dust clouds are located at the termini of the partial HI ring that surrounds Cen A which is also where the gas column density peaks... abridged

We describe the generation of single-band point source catalogues from submillimetre Herschel-SPIRE observations taken as part of the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found by means of peak-finding and the fitting of a Gaussian point-response function. With highly-confused images, careful checks must be made on the completeness and flux density accuracy of the detected sources. This is done by injecting artificial sources into the images and analysing the resulting catalogues. Measured flux densities at which 50 per cent of injected sources result in good detections at (250, 350, 500) {\mu}m range from (11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depending on the depth of the observation (where a `good' detection is taken to be one with positional offset less than one full-width half-maximum of the point-response function, and with the measured flux density within a factor of 2 of the flux density of the injected source). This paper acts as a reference for the 2010 July HerMES public data release.

Exploiting the Herschel-ATLAS Science Demonstration Phase (SDP) survey data, we have determined the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 micron and at several redshifts z>1, for bright sub-mm galaxies with star formation rates (SFR) >100 M_sun/yr. We find that the evolution of the comoving LF is strong up to z~2.5, and slows down at higher redshifts. From the LFs and the information on halo masses inferred from clustering analysis, we derived an average relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main episode of dust-enshrouded star formation in massive halos (M_H>3*10^12 M_sun) amounts to ~7*10^8 yr. Given the SFRs, which are in the range 10^2-10^3 M_sun/yr, this timescale implies final stellar masses of order of 10^11-10^12 M_sun. The corresponding stellar mass function matches the observed mass function of passively evolving galaxies at z>1. The comparison of the statistics for sub-mm and UV selected galaxies suggests that the dust-free, UV bright phase, is >10^2 times shorter than the sub-mm bright phase, implying that the dust must form soon after the onset of star formation. Using a single reference Spectral Energy Distribution (SED; the one of the z~2.3 galaxy SMM J2135-0102), our simple physical model is able to reproduce not only the LFs at different redshifts > 1 but also the counts at wavelengths ranging from 250 micron to ~1 mm. Owing to the steepness of the counts and their relatively broad frequency range, this result suggests that the dispersion of sub-mm SEDs of z>1 galaxies around the reference one is rather small.

The Herschel survey, H-ATLAS, with its large areal coverage, has recently discovered a number of bright, strongly lensed high-z submillimeter galaxies. The strong magnification makes it possible to study molecular species other than CO, which are otherwise difficult to observe in high-z galaxies. Among the lensed galaxies already identified by H-ATLAS, the source J090302.9-014127B (SDP.17b) at z = 2.305 is remarkable due to its excitation conditions and a tentative detection of the H2O 202-111 emission line (Lupu et al. 2010). We report observations of this line in SDP.17b using the IRAM interferometer equipped with its new 277- 371GHz receivers. The H2O line is detected at a redshift of z = 2.3049+/-0.0006, with a flux of 7.8+/-0.5 Jy km s-1 and a FWHM of 250+/-60 km s-1. The new flux is 2.4 times weaker than the previous tentative detection, although both remain marginally consistent within 1.6-sigma. The intrinsic line luminosity and ratio of H2O(202-111)/CO8-7 seem comparable with those of the nearby starburst/enshrouded-AGN Mrk 231, suggesting that SDP.17b could also host a luminous AGN. The detection of a strong H2O 202-111 line in SDP.17b implies an efficient excitation mechanism of the water levels that must occur in very dense and warm interstellar gas.

We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) and radio follow-up observations of two Herschel-ATLAS (H-ATLAS) detected strongly lensed distant galaxies. In one of the targeted galaxies H-ATLAS J090311.6+003906 (SDP.81) we detect [OIII] 88\mum and [CII] 158\mum lines at a signal-to-noise ratio of ~5. We do not have any positive line identification in the other fainter target H-ATLAS J091305.0-005343 (SDP.130). Currently SDP.81 is the faintest sub-mm galaxy with positive line detections with the FTS, with continuum flux just below 200 mJy in the 200-600 \mum wavelength range. The derived redshift of SDP.81 from the two detections is z=3.043 +/-0.012, in agreement with ground-based CO measurements. This is the first detection by Herschel of the [OIII] 88\mum line in a galaxy at redshift higher than 0.05. Comparing the observed lines and line ratios with a grid of photo-dissociation region (PDR) models with different physical conditions, we derive the PDR cloud density n ~ 2000 cm^{-3} and the far-UV ionizing radiation field G_0 ~ 200 (in units of the Habing field -- the local Galactic interstellar radiation field of 1.6x10^{-6} W/m^2). Using the CO derived molecular mass and the PDR properties we estimate the effective radius of the emitting region to be 500-700 pc. These characteristics are typical for star-forming, high redshift galaxies. The radio observations indicate that SDP.81 deviates significantly from the local FIR/radio correlation, which hints that some fraction of the radio emission is coming from an AGN. The constraints on the source size from millimiter-wave observations put a very conservative upper limit of the possible AGN contribution to less than 33%. These indications, together with the high [OIII]/FIR ratio and the upper limit of [OI] 63\mum/[CII] 158\mum suggest that some fraction of the ionizing radiation is likely to originate from an AGN.

We present the first direct and unbiased measurement of the evolution of the dust mass function of galaxies over the past 5 billion years of cosmic history using data from the Science Demonstration Phase of the Herschel-ATLAS. The sample consists of galaxies selected at 250{\mu}m which have reliable counterparts from SDSS at z < 0.5, and contains 1867 sources. Dust masses are calculated using both a single temperature grey-body model for the spectral energy distribution and also using a model with multiple temperature components. The dust temperature for either model shows no trend with redshift. Splitting the sample into bins of redshift reveals a strong evolution in the dust properties of the most massive galaxies. At z = 0.4 - 0.5, massive galaxies had dust masses about five times larger than in the local Universe. At the same time, the dust-to-stellar mass ratio was about 3-4 times larger, and the optical depth derived from fitting the UV-sub-mm data with an energy balance model was also higher. This increase in the dust content of massive galaxies at high redshift is difficult to explain using standard dust evolution models and requires a rapid gas consumption timescale together with either a more top-heavy IMF, efficient mantle growth, less dust destruction or combinations of all three. This evolution in dust mass is likely to be associated with a change in overall ISM mass, and points to an enhanced supply of fuel for star formation at earlier cosmic epochs.

We report ground-based follow-up observations of the exceptional source, ID141, one the brightest sources detected so far in the H-ATLAS cosmological survey. ID141 was observed using the IRAM 30-meter telescope and Plateau de Bure interferometer (PdBI), the Submillimeter Array (SMA) and the Atacama Pathfinder Experiment (APEX) submillimeter telescope to measure the dust continuum and emission lines of the main isotope of carbon monoxide and carbon ([C I] and [C II]). The detection of strong CO emission lines with the PdBI confirms that ID141 is at high redshift (z=4.243 +/- 0.001). The strength of the continuum and emission lines suggests that ID141 is gravitationally lensed. The width (Delta V (FWHM) ~ 800 km/s}) and asymmetric profiles of the CO and carbon lines indicate orbital motion in a disc or a merger. The properties derived for ID141 are compatible with a ultraluminous (L_FIR ~ 8.5 +/- 0.3 x 10^13/mu_L Lsun, where mu_L is the amplification factor, dense (n ~ 10^4 cm^-3) and warm (T_kin ~ 40K) starburst galaxy, with an estimated star-formation rate of (0.7 to 1.7) x 10^4/mu_L Msun/yr. The carbon emission lines indicate a dense (n ~ 10^4 cm^-3) Photo-Dominated Region, illuminated by a far-UV radiation field a few thousand times more intense than that in our Galaxy. In conclusion, the physical properties of the high-z galaxy, ID141, are remarkably similar to those of local ultraluminous infrared galaxies.

We present a technique to identify optical counterparts of 250 um-selected sources from the Herschel-ATLAS survey. Of the 6621 250 um > 32 mJy sources in our science demonstration catalogue we find that ~60 percent have counterparts brighter than r=22.4 mag in the Sloan Digital Sky Survey. Applying a likelihood ratio technique we are able to identify 2423 of the counterparts with a reliability R > 0.8. This is approximately 37 percent of the full 250 micron catalogue. We have estimated photometric redshifts for each of these 2423 reliable counterparts, while 1099 also have spectroscopic redshifts collated from several different sources, including the GAMA survey. We estimate the completeness of identifying counterparts as a function of redshift, and present evidence that 250 um-selected Herschel-ATLAS galaxies have a bimodal redshift distribution. Those with reliable optical identifications have a redshift distribution peaking at z ~ 0.25 +/- 0.05, while sub-mm colours suggest that a significant fraction with no counterpart above the r-band limit have z > 1. We also suggest a method for selecting populations of strongly-lensed high redshift galaxies. Our identifications are matched to UV--NIR photometry from the GAMA survey, and these data are available as part of the Herschel-ATLAS public data release.

For the first time, we investigate the X-ray/infrared (IR) correlation for star-forming galaxies at z~1, using SPIRE submm data from the recently-launched Herschel Space Observatory and deep X-ray data from the 2Ms Chandra deep field north (CDFN) survey. We examine the X-ray/IR correlation in the soft X-ray (SX, 0.5-2 keV) and hard X-ray (HX, 2-10 keV) bands by comparing our z~1 SPIRE-detected star-forming galaxies (SFGs) to equivalently IR-luminous (L_IR >10^10 L_sun) samples in the local/low redshift Universe. Our results suggest that the X-ray/IR properties of the SPIRE SFGs are on average similar to those of their local counterparts, as we find no evidence for evolution in the L_SX/L_IR and L_HX/L_IR ratios with redshift. We note however, that at all redshifts, both L_SX/L_IR and L_HX/L_IR are strongly dependent on IR luminosity, with luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs,L_IR >10^11 L_sun) having up to an order of magnitude lower values than normal infrared galaxies (L_IR <10^11 L_sun). We derive a L_SX-L_IR relation and confirm the applicability of an existing L_HX-L_IR relation for both local and distant LIRGs and ULIRGs, consistent with a scenario where X-ray luminosity is correlated with the star-formation rate (SFR).

We investigate the ultraviolet and optical properties and environment of low redshift galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration data. We use the Sloan Digital Sky Survey seventh release and the Galaxy And Mass Assembly database to select galaxies with r_Petro < 19.0 mag in the redshift range 0.02 < z < 0.2 and look for their submillimeter counterparts in H-ATLAS. Our results show that at low redshift, H-ATLAS detects mainly blue/star-forming galaxies with a minor contribution from red systems which are highly obscured by dust. In addition we find that the colour of a galaxy rather than the local density of its environment determines whether it is detectable by H-ATLAS. The average dust temperature of galaxies that are simultaneously detected by both PACS and SPIRE is 25K \pm 4K, independent of environment. This analysis provides a glimpse of the potential of the H-ATLAS data to investigate the submillimeter properties of galaxies in the local universe.

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

We report the detection of a significant excess in the surface density of far-infrared sources from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) within ~1 Mpc of the centres of 66 optically-selected clusters of galaxies in the SDSS with <z>~0.25. From the analysis of the multiwavelength properties of their counterparts we conclude that the far-infrared emission is associated with dust-obscured star formation and/or active galactic nuclei within galaxies in the clusters themselves. The excess reaches a maximum at a radius of ~0.8 Mpc, where we find 1.0\pm0.3 S_250um>34 mJy sources on average per cluster above what would be expected for random field locations. If the far-infrared emission is dominated by star formation (as opposed to AGN) then this corresponds to an average star formation rate of ~7 M_sun/yr per cluster in sources with L_IR>5d10 L_sun. Although lensed sources make a negligible contribution to the excess signal, a fraction of the sources around the clusters could be gravitationally lensed, and we have identified a sample of potential cases of cluster-lensed Herschel sources that could be targeted in follow-up studies.

As part of the Herschel Multi-tiered Extragalactic Survey we have investigated the rest-frame far-infrared (FIR) properties of a sample of more than 4800 Lyman Break Galaxies (LBGs) in the Great Observatories Origins Deep Survey North field. Most LBGs are not detected individually, but we do detect a sub-sample of 12 objects at 0.7 < z < 1.6 and one object at z ~ 2.0. The ones detected by Herschel SPIRE have redder observed NUV-U and U-R colors than the others, while the undetected ones have colors consistent with average LBGs at z > 2.5. The UV-to-FIR spectral energy distributions of the objects detected in the rest-frame FIR are investigated using the code CIGALE to estimate physical parameters. We find that LBGs detected by SPIRE are high mass, luminous infrared galaxies. It appears that LBGs are located in a triangle-shaped region in the A_FUV vs. Log L_FUV diagram limited by A_FUV=0 at the bottom and by a diagonal following the temporal evolution of the most massive galaxies from the bottom-right to the top-left of the diagram. This upper envelop can be used as upper limits for the UV dust attenuation as a function of L_FUV}. The limits of this region are well explained using a closed-box model, where the chemical evolution of galaxies produces metals, which in turn lead to higher dust attenuation when the galaxies age.

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

We report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

The Herschel-ATLAS is a survey of 550 square degrees with the Herschel Space Observatory in five far--infrared and submillimetre bands. The first data for the survey, observations of a field 4x4 sq. degrees in size, were taken during the Science Demonstration Phase, and reach a 5 sigma noise level of 33 mJy/beam at 250 microns. This paper describes the source extraction methods used to create the corresponding Science Demonstration Phase catalogue, which contains 6876 sources, selected at 250 microns, within ~14 sq. degrees. SPIRE sources are extracted using a new method specifically developed for Herschel data; PACS counterparts of these sources are identified using circular apertures placed at the SPIRE positions. Aperture flux densities are measured for sources identified as extended after matching to optical wavelengths. The reliability of this catalogue is also discussed, using full simulated maps at the three SPIRE bands. These show that a significant number of sources at 350 and 500 microns have undergone flux density enhancements of up to a factor of ~2, due mainly to source confusion. Correction factors are determined for these effects. The SDP dataset and corresponding catalogue will be available from http://www.h-atlas.org/.

We use the science demonstration field data of the Herschel-ATLAS to study how star formation, traced by the far-infrared Herschel data, is related to both the accretion luminosity and redshift of quasars selected from the Sloan Digital Sky Survey and the 2SLAQ survey. By developing a maximum likelihood estimator to investigate the presence of correlations between the far-infrared and optical luminosities we find evidence that the star-formation in quasar hosts is correlated with both redshift and quasar accretion luminosity. Assuming a relationship of the form L_IR \propto L_QSO^{\theta} (1 + z)^{\zeta}, we find {\theta} = 0.22 +/- 0.08 and {\zeta} = 1.6 +/- 0.4, although there is substantial additional uncertainty in {\zeta} of order +/- 1, due to uncertainties in the host galaxy dust temperature. We find evidence for a large intrinsic dispersion in the redshift dependence, but no evidence for intrinsic dispersion in the correlation between L_QSO and L_IR, suggesting that the latter may be due to a direct physical connection between star formation and black hole accretion. This is consistent with the idea that both the quasar activity and star formation are dependent on the same reservoir of cold gas, so that they are both affected by the influx of cold gas during mergers or heating of gas via feedback processes.

We use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, $\beta$, the Balmer decrement, and the far infrared (IR) to $150\,$nm far ultraviolet (FUV) luminosity ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order to identify possible systematics in these various measures. We reiterate the finding of other authors that there is a large scatter between the Balmer decrement and the $\beta$ parameter, and that $\beta$ may be poorly constrained when derived from only two broad passbands in the UV. We also emphasise that FUV derived SFRs, corrected for dust obscuration using $\beta$, will be overestimated unless a modified relation between $\beta$ and the attenuation factor is used. Even in the optimum case, the resulting SFRs have a significant scatter, well over an order of magnitude. While there is a stronger correlation between the IR to FUV luminosity ratio and $\beta$ parameter than with the Balmer decrement, neither of these correlations are particularly tight, and dust corrections based on $\beta$ for high redshift galaxy SFRs must be treated with caution. We conclude with a description of the extent to which the different obscuration measures are consistent with each other as well as the effects of including other galactic properties on these correlations.

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

We present physical properties of two submillimeter selected gravitationally lensed sources, identified in the Herschel Astrophysical Terahertz Large Area Survey. These submillimeter galaxies (SMGs) have flux densities > 100 mJy at 500 um, but are not visible in existing optical imaging. We fit light profiles to each component of the lensing systems in Spitzer IRAC 3.6 and 4.5 um data and successfully disentangle the foreground lens from the background source in each case, providing important constraints on the spectral energy distributions (SEDs) of the background SMG at rest-frame optical-near-infrared wavelengths. The SED fits show that these two SMGs have high dust obscuration with Av ~4 to 5 and star formation rates of ~100 M_sun/yr. They have low gas fractions and low dynamical masses compared to 850 um selected galaxies.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) has uncovered a population of strongly-lensed submillimeter galaxies (SMGs). The Zpectrometer instrument on the Green Bank Telescope (GBT) was used to measure the redshifts and constrain the masses of the cold molecular gas reservoirs for two candidate high-redshift lensed sources. We derive CO(1-0) redshifts of z=3.042+/-0.001 and z=2.625+/-0.001, and measure molecular gas masses of (1--3)x10^{10}Msun, corrected for lens amplification and assuming a conversion factor of alhpa=0.8 Msun(K km/s pc^2)^{-1}. We find typical L(IR)/L'(CO) ratios of 120+/-40 and 140+/-50 Lsun (K km/s pc^2)^{-1}, which are consistent with those found for local ULIRGs and other high-redshift SMGs. From analysis of published data, we find no evidence for enhanced L(IR)/L'(CO(1-0)) ratios for the SMG population in comparison to local ULIRGs. The GBT results highlight the power of using the CO lines to derive blind redshifts, which is challenging for the SMG population at optical wavelengths given their high obscuration.

We carry out a multi-wavelength study of individual galaxies detected by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) and identified at other wavelengths, using data spanning the radio to the ultraviolet (UV). We develop a Monte Carlo method to account for flux boosting, source blending, and correlations among bands, which we use to derive deboosted far-infrared (FIR) luminosities for our sample. We estimate total star-formation rates for BLAST counterparts with z < 0.9 by combining their FIR and UV luminosities. Star formation is heavily obscured at L_FIR > 10^11 L_sun, z > 0.5, but the contribution from unobscured starlight cannot be neglected at L_FIR < 10^11 L_sun, z < 0.25. We assess that about 20% of the galaxies in our sample show indication of a type-1 active galactic nucleus (AGN), but their submillimeter emission is mainly due to star formation in the host galaxy. We compute stellar masses for a subset of 92 BLAST counterparts; these are relatively massive objects, with a median mass of ~10^11 M_sun, which seem to link the 24um and SCUBA populations, in terms of both stellar mass and star-formation activity. The bulk of the BLAST counterparts at z<1 appear to be run-of-the-mill star-forming galaxies, typically spiral in shape, with intermediate stellar masses and practically constant specific star-formation rates. On the other hand, the high-z tail of the BLAST counterparts significantly overlaps with the SCUBA population, in terms of both star-formation rates and stellar masses, with observed trends of specific star-formation rate that support strong evolution and downsizing.

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity \nu L_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint (\nu L_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.