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.

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.

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.

The $Herschel$ Fornax Cluster Survey (HeFoCS) is a deep, far-infrared (FIR) survey of the Fornax cluster. The survey is in 5 $Herschel$ bands (100 - 500 $\mu$m) and covers an area of 16 deg$^2$ centred on NGC1399. This paper presents photometry, detection rates, dust masses and temperatures using an optically selected sample from the Fornax Cluster Catalogue (FCC). Our results are compared with those previously obtained using data from the $Herschel$ Virgo Cluster Survey (HeViCS). In Fornax, we detect 30 of the 237 (13%) optically selected galaxies in at least one $Herschel$ band. The global detection rates are significantly lower than Virgo, reflecting the morphological make up of each cluster - Fornax has a lower fraction of late-type galaxies. For galaxies detected in at least 3 bands we fit a modified blackbody with a $\beta = 2$ emissivity. Detected early-type galaxies (E/S0) have a mean dust mass, temperature, and dust-to-stars ratio of $\log_{10}(<M_{dust}>/\mathrm{M_{\odot}}) = 5.82 \pm 0.20$, $<T_{dust}> = 20.82 \pm 1.77$K, and $\log_{10}(M_{dust}/M_{stars}) = -3.87 \pm 0.28$, respectively. Late-type galaxies (Sa to Sd) have a mean dust mass, temperature, and dust-to-stars ratio of $\log_{10}(<M_{dust}>/\mathrm{M_{\odot}}) = 6.54 \pm 0.19$, $<T_{dust}> = 17.47 \pm 0.97$K, and $\log_{10}(M_{dust}/M_{stars}) = -2.93 \pm 0.09$, respectively. The different cluster environments seem to have had little effect on the FIR properties of the galaxies and so we conclude that any environment dependent evolution, has taken place before the cluster was assembled.

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.

Herschel FIR observations are used to construct Virgo cluster galaxy luminosity functions and to show that the cluster lacks the very bright and the numerous faint sources detected in field galaxy surveys. The far-infrared SEDs are fitted to obtain dust masses and temperatures and the dust mass function. The cluster is over dense in dust by about a factor of 100 compared to the field. The same emissivity (beta) temperature relation applies for different galaxies as that found for different regions of M31. We use optical and HI data to show that Virgo is over dense in stars and atomic gas by about a factor of 100 and 20 respectively. Metallicity values are used to measure the mass of metals in the gas phase. The mean metallicity is about 0.7 solar and 50% of the metals are in the dust. For the cluster as a whole the mass density of stars in galaxies is 8 times that of the gas and the gas mass density is 130 times that of the metals. We use our data to consider the chemical evolution of the individual galaxies, inferring that the measured variations in effective yield are due to galaxies having different ages, being affected to varying degrees by gas loss. Four galaxy scaling relations are considered: mass-metallicity, mass-velocity, mass-star formation rate and mass-radius - we suggest that initial galaxy mass is the prime driver of a galaxy's ultimate destiny. Finally, we use X-ray observations and galaxy dynamics to assess the dark and baryonic matter content compared to the cosmological model.

We present Herschel observations of the Fornax cluster at 100, 160, 250, 350 and 500u with a spatial resolution of 7 - 36 arc sec (10" = 1 kpc at d_Fornax=17.9 Mpc). We define a sample of 11 bright galaxies, selected at 500u, directly comparable with our past work on Virgo. We find good agreement with previous observations made by IRAS and Planck. The FIR luminosity density is higher (factor of three) in Fornax compared to Virgo. The 100u (42.5-122.5u) luminosity is two orders of magnitude larger in Fornax than in the local field as measured by IRAS. Using stellar (L_{0.4-2.5}) and FIR (L_{100-500}) luminosities we estimate a mean optical depth of tau=0.4+/-0.1 - the same value as Virgo. For 10 of the 11 galaxies (NGC1399 excepted) we fit a modified blackbody curve (beta=2.0) to the SEDs to derive dust masses and temperatures of 10^{6.54-8.35} M_0 and T=14.6-24.2K respectively, comparable to Virgo. The derived stars-to-gas(atomic) and gas(atomic)-to-dust ratios vary from 1.1-67.6 and 9.8-436.5 respectively, again consistent with Virgo. Fornax is a mass overdensity in stars and dust of about 120 compared to the local field (30 for Virgo). Fornax and Virgo are both a factor of 6 lower over densities in gas(atomic) than in stars and dust indicating loss of gas, but not dust and stars, in the cluster environment. As the brightest FIR source in either Fornax and Virgo, NGC1365 is detected by Planck. The Planck data fit the PACS/SPIRE SED out to 1382u with no evidence of other sources of emission ('spinning dust', free-free, synchrotron). At the opposite end of the scale NGC1399 is detected only at 500$\mu$m with the emission probably arising from the nuclear radio source rather than inter-stellar dust.