We used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02<z<0.5$, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500$\,\mu$m. SPIRE flux densities were also combined with Spitzer photometry and multi-wavelength archival data to perform a complete SED fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (8-1000$\,\mu$m) luminosity functions and their low-$z$ evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density (LLD) and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The luminosity functions show significant and rapid luminosity evolution already at low redshifts, $0.02<z<0.2$, with L$_{IR}^* \propto (1+z)^{6.0\pm0.4}$ and $\Phi_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $\Phi_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,\mu$m LFs respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter IMF and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate a SFRD scaling of SFRD$_0+0.08 z$, where SFRD$_0\simeq (1.9\pm 0.03)\times 10^{-2} [\mathrm{M}_\odot\,\mathrm{Mpc}^{-3}]$ is our total SFRD estimate at $z\sim0.02$.

We present three independent catalogs of point-sources extracted from SPIRE images at 250, 350, and 500 micron as a part of the Herschel Virgo Cluster Survey (HeViCS). The source positions are determined by estimating the likelihood to be a real source for each peak on the maps and the flux densities are estimated using the sourceExtractorTimeline, a timeline-based point source fitter. Afterwards, each source is subtracted from the maps, removing a Gaussian function in every position with the full width half maximum equal to that estimated in sourceExtractorTimeline. This procedure improves the robustness of our algorithm in terms of source identification. The HeViCS catalogs contain about 52000, 42200, and 18700 sources selected at 250, 350, and 500 micron above 3sigma and are ~ 75%, 62%, and 50% complete at flux densities of 20 mJy at 250, 350, 500 micron, respectively. We then measured source number counts at 250, 350, and 500 micron and we also cross-correlated the catalogs with the Sloan Digital Sky Survey to investigate the redshift distribution of the nearby sources. From this cross-correlation, we select ~2000 sources with reliable fluxes and a high signal-to-noise ratio, finding an average redshift z~0.3+/-0.22. The number counts at 250, 350, and 500 micron show an increase in the slope below 200 mJy, indicating a strong evolution in number of density for galaxies at these fluxes. In general, models tend to overpredict the counts at brighter flux densities, underlying the importance of studying the Rayleigh-Jeans part of the spectral energy distribution to refine the theoretical recipes of the models. Our iterative method for source identification allowed the detection of a family of 500 micron sources that are not foreground objects belonging to Virgo and not found in other catalogs.

Photometric instruments operating at far infrared to millimetre wavelengths often have broad spectral passbands (central wavelength/bandwidth ~ 3 or less), especially those operating in space. A broad passband can result in significant variation of the beam profile and aperture efficiency across the passband, effects which thus far have not generally been taken into account in the flux calibration of such instruments. With absolute calibration uncertainties associated with the brightness of primary calibration standards now in the region of 5% or less, variation of the beam properties across the passband can be a significant contributor to the overall calibration accuracy for extended emission. We present a calibration framework which takes such variations into account for both antenna-coupled and absorber-coupled focal plane architectures. The scheme covers point source and extended source cases, and also the intermediate case of a semi-extended source profile. We apply the new method to the Herschel-SPIRE space-borne photometer.

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 present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.

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.

We present new photometric data from our Herschel Key Programme, the Dwarf Galaxy Survey (DGS), dedicated to the observation of the gas and dust in 48 low-metallicity environments. They were observed with PACS and SPIRE onboard Herschel at 70,100,160,250,350, and 500 microns. We focus on a systematic comparison of the derived FIR properties (FIR luminosity, dust mass, dust temperature and emissivity index) with more metal-rich galaxies and investigate the detection of a potential submm excess. The data reduction method is adapted for each galaxy to derive the most reliable photometry from the final maps. PACS flux densities are compared with the MIPS 70 and 160 microns bands. We use colour-colour diagrams and modified blackbody fitting procedures to determine the dust properties of the DGS galaxies. We also include galaxies from the Herschel KINGFISH sample, containing more metal-rich environments, totalling 109 galaxies. The location of the DGS galaxies on Herschel colour-colour diagrams highlights the differences in global environments of low-metallicity galaxies. The dust in DGS galaxies is generally warmer than in KINGFISH galaxies (T_DGS~32 K, T_KINGFISH~23 K). The emissivity index, beta, is ~1.7 in the DGS, but metallicity does not make a strong effect on beta. The dust-to-stellar mass ratio is lower in low-metallicity galaxies: M_dust/M_star~0.02% for the DGS vs 0.1% for KINGFISH. Per unit dust mass, dwarf galaxies emit ~6 times more in the FIR than higher metallicity galaxies. Out of the 22 DGS galaxies detected at 500 micron, 41% present an excess in the submm not explained by our dust SED model. The excess mainly appears in lower metallicity galaxies (12+log(O/H) < 8.3), and the strongest excesses are detected in the most metal-poor galaxies. We stress the need for observations longwards of the Herschel wavelengths to detect any submm excess appearing beyond 500 micron.

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.

The Dwarf Galaxy Survey (DGS) program is studying low-metallicity galaxies using 230h of far-infrared (FIR) and submillimetre (submm) photometric and spectroscopic observations of the Herschel Space Observatory and draws to this a rich database of a wide range of wavelengths tracing the dust, gas and stars. This sample of 50 galaxies includes the largest metallicity range achievable in the local Universe including the lowest metallicity (Z) galaxies, 1/50 Zsun, and spans 4 orders of magnitude in star formation rates. The survey is designed to get a handle on the physics of the interstellar medium (ISM) of low metallicity dwarf galaxies, especially on their dust and gas properties and the ISM heating and cooling processes. The DGS produces PACS and SPIRE maps of low-metallicity galaxies observed at 70, 100, 160, 250, 350, and 500 mic with the highest sensitivity achievable to date in the FIR and submm. The FIR fine-structure lines, [CII] 158 mic, [OI] 63 mic, [OI] 145 mic, [OIII] 88 mic, [NIII] 57 mic and [NII] 122 and 205 mic have also been observed with the aim of studying the gas cooling in the neutral and ionized phases. The SPIRE FTS observations include many CO lines (J=4-3 to J=13-12), [NII] 205 mic and [CI] lines at 370 and 609 mic. This paper describes the sample selection and global properties of the galaxies, the observing strategy as well as the vast ancillary database available to complement the Herschel observations. The scientific potential of the full DGS survey is described with some example results included.

We describe the procedure used to flux calibrate the three-band submillimetre photometer in the Spectral and Photometric Imaging REceiver (SPIRE) instrument on the Herschel Space Observatory. This includes the equations describing the calibration scheme, a justification for using Neptune as the primary calibration source, a description of the observations and data processing procedures used to derive flux calibration parameters (for converting from voltage to flux density) for every bolometer in each array, an analysis of the error budget in the flux calibration for the individual bolometers, and tests of the flux calibration on observations of primary and secondary calibrators. The procedure for deriving the flux calibration parameters is divided into two parts. In the first part, we use observations of astronomical sources in conjunction with the operation of the photometer internal calibration source to derive the unscaled derivatives of the flux calibration curves. To scale the calibration curves in Jy/beam/V, we then use observations of Neptune in which the beam of each bolometer is mapped using Neptune observed in a very fine scan pattern. The total instrumental uncertainties in the flux calibration for the individual bolometers is ~0.5% for most bolometers, although a few bolometers have uncertainties of ~1-5% because of issues with the Neptune observations. Based on application of the flux calibration parameters to Neptune observations performed using typical scan map observing modes, we determined that measurements from each array as a whole have instrumental uncertainties of 1.5%. This is considerably less than the absolute calibration uncertainty associated with the model of Neptune, which is estimated at 4%.

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

[Abridged] We present an integrated photometric spectral energy distribution (SED) of the Magellanic-type galaxy NGC 4449 from the far-ultraviolet (UV) to the submillimetre, including new observations acquired by the Herschel Space Observatory. We include integrated UV photometry from the Swift Ultraviolet and Optical Telescope using a measurement technique which is appropriate for extended sources with coincidence loss. In this paper, we examine the available multiwavelength data to infer a range of ages, metallicities and star formation rates for the underlying stellar populations, as well as the composition and the total mass of dust in NGC 4449. We present an iterative scheme, which allows us to build an in-depth and multicomponent representation of NGC 4449 `bottom-up', taking advantage of the broad capabilities of the photoionization and radiative transfer code MOCASSIN (MOnte CArlo SimulationS of Ionized Nebulae). We fit the observed SED, the global ionization structure and the emission line intensities, and infer a recent SFR of 0.4 Msolar/yr and a total stellar mass of approximately 1e9 Msolar emitting with a bolometric luminosity of 5.7e9 Lsolar. Our fits yield a total dust mass of 2.9e6 Msolar including 2 per cent attributed to polycyclic aromatic hydrocarbons. We deduce a dust to gas mass ratio of 1/190 within the modelled region. While we do not consider possible additional contributions from even colder dust, we note that including the extended HI envelope and the molecular gas is likely to bring the ratio down to as low as ~ 1/800.

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.

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.1<z<2 and determine how these evolve with cosmic time. The unique angle of this work is the rigorous analysis of survey selection effects, making this the first study of the star-formation-dominated, IR-luminous population within a framework almost entirely free of selection biases. We find that IR-luminous galaxies have SEDs with broad far-IR peaks characterised by cool/extended dust emission and average dust temperatures in the 25-45K range. Hot (T>45K) SEDs and cold (T<25K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity-temperature (L-T) relation, where the average dust temperature of log[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase.

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.

The Herschel Virgo Cluster Survey (HeViCS) is the deepest, confusion-limited survey of the Virgo Cluster at far-infrared (FIR) wavelengths. The entire survey at full depth covers $\sim$55 sq. deg. in 5 bands (100-500 \micron), encompassing the areas around the central dominant elliptical galaxies (M87, M86 & M49) and extends as far as the NW cloud, the W cloud and the Southern extension. The survey extends beyond this region with lower sensitivity so that the total area covered is 84 sq. deg. In this paper we describe the data, the data acquisition techniques and present the detection rates of the optically selected Virgo Cluster Catalogue (VCC). We detect 254 (34%) of 750 VCC galaxies found within the survey boundary in at least one band and 171 galaxies are detected in all five bands. For the remainder of the galaxies we have measured strict upper limits for their FIR emission. The population of detected galaxies contains early- as well as late-types although the latter dominate the detection statistics. We have modelled 168 galaxies, showing no evidence of a strong synchrotron component in their FIR spectra, using a single-temperature modified blackbody spectrum with a fixed emissivity index ($\beta = 2$). A study of the $\chi^2$ distribution indicates that this model is not appropriate in all cases, and this is supported by the FIR colours which indicate a spread in $\beta$=1--2. Statistical comparison of the dust mass and temperature distributions from 140 galaxies with $\chi^2_{dof=3} < 7.8$ (95% confidence level) shows that late-types have typically colder, more massive dust reservoirs; the early-type dust masses have a mean of ${\rm log}(<M> / M_{\sun}) = 6.3 \pm 0.3 $, while for late-types ${\rm log}(<M> / M_{\sun}) =7.1 \pm 0.1$... (abridged)

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 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 investigate the dust-to-gas mass ratio and the environmental effects on the various components of the interstellar medium for a spatially resolved sample of Virgo spirals. We have used the IRAM-30m telescope to map over their full extent NGC 4189, NGC 4298, NGC 4388, and NGC 4299 in the 12CO(1-0) and the 12CO(2-1) lines. We observed the same lines in selected regions of NGC 4351, NGC 4294, and NGC 4424. The CO observations are combined with Herschel maps in 5 bands between 100-500 {\mu}m from the HeViCS survey, and with HI data from the VIVA survey, to obtain spatially resolved dust and gas distributions. We studied the environmental dependencies by adding to our sample eight galaxies with 12CO(1-0) maps from the literature. We estimate the integrated mass of molecular hydrogen for the galaxies observed in the CO lines. We find molecular-to-total gas mass fractions between 0.04 \leq fmol \leq 0.65, with the lowest values for the dimmest galaxy in the B-band. The integrated dust-to-gas ratio ranges between 0.011 and 0.004. For the 12 mapped galaxies we derive the radial distributions of the atomic gas, molecular gas, and dust. We also study the effect of different CO-to-H2 conversion factors. Both the molecular gas and the dust distributions show steeper radial profiles for HI-deficient galaxies and the average dust-to-gas ratio for these galaxies increases or stays radially constant. On scales of \sim 3 kpc, we find a strong correlation between the molecular gas and the 250 micron surface brightness that is tighter than average for non-deficient galaxies. The correlation becomes linear if we consider the total gas surface mass density. However, the inclusion of atomic hydrogen does not improve the statistical significance of the correlation. The environment can modify the distributions of molecules and dust within a galaxy, although these components are more tightly bound than the atomic gas.

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

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

Using the far-infrared emission, as observed by the Herschel Virgo Cluster Survey (HeViCS), and the integrated HI and CO brightness, we infer the dust and total gas mass for a magnitude limited sample of 35 metal rich spiral galaxies in Virgo. The CO flux correlates tightly and linearly with far-infrared fluxes observed by Herschel. Molecules in these galaxies are more closely related to cold dust rather than to dust heated by star formation or to optical/NIR brightness. We show that dust mass establishes a stronger correlation with the total gas mass than with the atomic or molecular component alone. The dust-to-gas ratio increases as the HI deficiency increases, but in highly HI deficient galaxies it stays constant. Dust is in fact less affected than atomic gas by weak cluster interactions, which remove most of the HI gas from outer and high latitudes regions. Highly disturbed galaxies, in a dense cluster environment, can instead loose a considerable fraction of gas and dust from the inner regions of the disk keeping constant the dust-to-gas ratio. There is evidence that the molecular phase is also quenched. This quencing becomes evident by considering the molecular gas mass per unit stellar mass. Its amplitude, if confirmed by future studies, highlights that molecules are missing in Virgo HI deficient spirals, but to a somewhat lesser extent than dust.

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.

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.

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.

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.

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.

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 describe the Herschel Virgo Cluster Survey (HeViCS) and the first data that cover the complete survey area (four 4 x 4 deg2 regions). We use these data to measure and compare the global far infrared properties of 78 optically bright galaxies that are selected at 500 \mum and detected in all five far-infrared bands. We show that our measurements and calibration are broadly consistent with previous data obtained by IRAS, ISO, Spitzer and Planck. We use SPIRE and PACS photometry data to produce 100, 160, 250, 350 and 500 \mum cluster luminosity distributions. These luminosity distributions are not power laws, but peaked, with small numbers of both faint and bright galaxies. We measure a cluster 100-500 micron far-infrared luminosity density of 1.6(7.0) \pm 0.2 x 10^9 Lsun/Mpc3. This compares to a cluster 0.4-2.5 \mum optical luminosity density of 5.0(20.0) x 10^9 Lsun/Mpc3, some 3.2(2.9) times larger than the far-infrared. A typical photon originates from an optical depth of 0.4\pm0.1. Most of our sample galaxies are well fitted by a single modified blackbody (beta=2), leading to a mean dust mass of log Mdust = 7.31 Msun and temperature of 20.0 K. We also derive both stellar and atomic hydrogen masses from which we calculate mean values for the stars:gas(atomic) and gas(atomic): dust mass ratios of 15.1 and 58.2 respectively. Using our derived dust, atomic gas and stellar masses we estimate cluster mass densities of 8.6(27.8) x 10^6, 4.6(13.9) x 10^8, 7.8(29.7) x 10^9 Msun/Mpc3, respectively for dust, atomic gas and stars. These values are higher than those derived for field galaxies by factors of 39(126), 6(18) and 34(129) respectively. In the above luminosity/mass densities are given using the whole sample with values in brackets using just those galaxies that lie between 17 and 23 Mpc. We provide a data table of flux densities in all the Herschel bands for all 78 bright Virgo cluster galaxies.

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.

We use Spitzer Space Telescope and Herschel Space Observatory far-infrared data along with ground-based optical and near-infrared data to understand how dust heating in the nearby face-on spiral galaxies M81, M83, and NGC 2403 is affected by the starlight from all stars and by the radiation from star forming regions. We find that 70/160 micron surface brightness ratios tend to be more strongly influenced by star forming regions. However, the 250/350 micron and 350/500 micron surface brightness ratios are more strongly affected by the light from the total stellar populations, suggesting that the dust emission at >250 microns originates predominantly from a component that is colder than the dust seen at <160 microns and that is relatively unaffected by star formation activity. We conclude by discussing the implications of this for modelling the spectral energy distributions of both nearby and more distant galaxies and for using far-infrared dust emission to trace star formation.

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.

Using Herschel data from the Open Time Key Project the Herschel Virgo Cluster Survey (HeViCS), we investigated the relationship between the metallicity gradients expressed by metal abundances in the gas phase as traced by the chemical composition of HII regions, and in the solid phase, as traced by the dust-to-gas mass ratio. We derived the radial gradient of the dust-to-gas mass ratio for all galaxies observed by HeViCS whose metallicity gradients are available in the literature. They are all late type Sbc galaxies, namely NGC4254, NGC4303, NGC4321, and NGC4501. We examined different dependencies on metallicity of the CO-to-H$_2$ conversion factor (\xco), used to transform the $^{12}$CO observations into the amount of molecular hydrogen. We found that in these galaxies the dust-to-gas mass ratio radial profile is extremely sensitive to choice of the \xco\ value, since the molecular gas is the dominant component in the inner parts. We found that for three galaxies of our sample, namely NGC4254, NGC4321, and NGC4501, the slopes of the oxygen and of the dust-to-gas radial gradients agree up to $\sim$0.6-0.7R$_{25}$ using \xco\ values in the range 1/3-1/2 Galactic \xco. For NGC4303 a lower value of \xco$\sim0.1\times$ 10$^{20}$ is necessary. We suggest that such low \xco\ values might be due to a metallicity dependence of \xco (from close to linear for NGC4254, NGC4321, and NGC4501 to superlinear for NGC4303), especially in the radial regions R$_G<$0.6-0.7R$_{25}$ where the molecular gas dominates. On the other hand, the outer regions, where the atomic gas component is dominant, are less affected by the choice of \xco, and thus we cannot put constraints on its value.

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.

The Herschel Dwarf Galaxy Survey investigates the interplay of star formation activity and the the metal-poor gas and dust of dwarf galaxies using FIR and submillimetre imaging spectroscopic and photometric observations in the 50 to 550mu window of the Herschel Space Observatory. The dust SEDs are well constrained with the new Herschel and MIR Spitzer data. A submillimetre excess is often found in low metallicity galaxies, which,if tracing very cold dust, would highlight large dust masses not easily reconciled in some cases, given the low metallicities and expected gas-to-dust mass ratios. The galaxies are also mapped in the FIR fine-structure lines (63 and 145mu OI, 158mu CII, 122 and 205mu NII, 88mu OIII) probing the low density ionised gas, the HII regions and photodissociation regions. While still early in the Herschel mission we can already see, along with earlier studies, that line ratios in the metal-poor ISM differ remarkably from those in the metal-rich starburst environments. In dwarf galaxies, L[CII]/L(CO) (>10^4) is at least an order of magnitude greater than in the most metal-rich starburst galaxies. The enhanced [CII] arises from the larger photodissociation region where H2, not traced by the CO, can exist. The 88mu [OIII] line usually dominates the FIR line emission over galaxy-wide scales in dwarf galaxies, not the 158mu [CII] line which is the dominant FIR cooling line in metal-rich galaxies. All of the FIR lines together can contribute 1% to 2% of the L(TIR). The Herschel Dwarf Galaxy survey will provide statistical information on the nature of the dust and gas in low metallicity galaxies, elucidating the origin of the submm excess in dwarf galaxies, and help determine a ([CII] +CO) to H2 conversion factor, thus providing observational constraints on chemical evolution models of galaxies.

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

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

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

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24{\mu}m data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24{\mu}m and 15/45 (all MIPS-24{\mu}m detections) are detected at SPIRE-250{\mu}m, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).

We have reduced the data taken with the Spectral and Photometric Imaging Receiver (SPIRE) photometer on board the Herschel Space Observatory in the Science Demonstration Phase (SDP) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We describe the data reduction, which poses specific challenges, both because of the sheer size of the data, and because only two scans are made for each region. We implement effective solutions to process the bolometric timelines into maps, and show that correlations among detectors are negligible, and that the photometer is stable on time scales up to 250 s. This is longer than the time the telescope takes to cross the observed sky region, and it allows us to use naive binning methods for an optimal reconstruction of the sky emission. The maps have equal contribution of confusion and white instrumental noise, and the latter is estimated to 5.3, 6.4, and 6.7 mJy/beam (1-{\sigma}), at 250, 350, and 500 \mu{m}, respectively. This pipeline is used to reduce other H-ATLAS observations, as they became available, and we discuss how it can be used with the optimal map maker implemented in the Herschel Interactive Processing Environment (HIPE), to improve computational efficiency and stability. The SDP dataset is available from http://www.h-atlas.org/.

Galactic cirrus emission at far-infrared wavelengths affects many extragalactic observations. Separating this emission from that associated with extragalactic objects is both important and difficult. In this paper we discuss a particular case, the M81 group, and the identification of diffuse structures prominent in the infrared, but also detected at optical wavelengths. The origin of these structures has previously been controversial, ranging from them being the result of a past interaction between M81 and M82 or due to more local Galactic emission. We show that over of order a few arcminute scales the far-infrared (Herschel 250 &\mu&m) emission correlates spatially very well with a particular narrow velocity (2-3 km/s) component of the Galactic HI. We find no evidence that any of the far-infrared emission associated with these features actually originates in the M81 group. Thus we infer that the associated diffuse optical emission must be due to galactic light back scattered off dust in our galaxy. Ultra-violet observations pick out young stellar associations around M81, but no detectable far-infrared emission. We consider in detail one of the Galactic cirrus features, finding that the far-infrared HI relation breaks down below arc minute scales and that at smaller scales there can be quite large dust temperature variations.

The origin of the far-infrared emission from the nearby radio galaxy M87 remains a matter of debate. Some studies find evidence of a far-infrared excess due to thermal dust emission, whereas others propose that the far-infrared emission can be explained by synchrotron emission without the need for an additional dust emission component. We observed M87 with PACS and SPIRE as part of the Herschel Virgo Cluster Survey (HeViCS). We compare the new Herschel data with a synchrotron model based on infrared, submm and radio data to investigate the origin of the far-infrared emission. We find that both the integrated SED and the Herschel surface brightness maps are adequately explained by synchrotron emission. At odds with previous claims, we find no evidence of a diffuse dust component in M87.

We describe the production and verification of sky maps of the five SPIRE fields observed as part of the Herschel Multi-tiered Extragalactic Survey (HerMES) during the Science Demonstration Phase (SDP) of the Herschel mission. We have implemented an iterative map-making algorithm (SHIM; The SPIRE-HerMES Iterative Mapper) to produce high fidelity maps that preserve extended diffuse emission on the sky while exploiting the repeated observations of the same region of the sky with many detectors in multiple scan directions to minimize residual instrument noise. We specify here the SHIM algorithm and outline the various tests that were performed to determine and characterize the quality of the maps and verify that the astrometry, point source flux and power on all relevant angular scales meets the needs of the HerMES science goals. These include multiple jackknife tests, determination of the map transfer function and detailed examination of the power spectra of both sky and jackknife maps. The map transfer function is approximately unity on scales from one arcminute to one degree. Final maps (v1.0), including multiple jackknives, as well as the SHIM pipeline, have been used by the HerMES team for the production of SDP papers.

We use the Herschel-ATLAS science demonstration data to investigate the star-formation properties of radio-selected galaxies in the GAMA-9h field as a function of radio luminosity and redshift. Radio selection at the lowest radio luminosities, as expected, selects mostly starburst galaxies. At higher radio luminosities, where the population is dominated by AGN, we find that some individual objects are associated with high far-infrared luminosities. However, the far-infrared properties of the radio-loud population are statistically indistinguishable from those of a comparison population of radio-quiet galaxies matched in redshift and K-band absolute magnitude. There is thus no evidence that the host galaxies of these largely low-luminosity (Fanaroff-Riley class I), and presumably low-excitation, AGN, as a population, have particularly unusual star-formation histories. Models in which the AGN activity in higher-luminosity, high-excitation radio galaxies is triggered by major mergers would predict a luminosity-dependent effect that is not seen in our data (which only span a limited range in radio luminosity) but which may well be detectable with the full Herschel-ATLAS dataset.

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

We use data from the Herschel-ATLAS to investigate the evolution of the far-infrared--radio correlation over the redshift range 0<z<0.5. Using the total far-infrared luminosity of all >5sigma sources in the Herschel-ATLAS Science Demonstration Field and cross-matching these data with radio data from the Faint Images of the Radio Sky at Twenty-Centimetres (FIRST) survey and the NRAO VLA Northern Sky Survey (NVSS), we obtain 104 radio counterparts to the Herschel sources. With these data we find no evidence for evolution in the far-infrared--radio correlation over the redshift range 0<z<0.5, where the median value for the ratio between far-infrared and radio luminosity, $q_{\rm IR}$, over this range is $q_{\rm IR} = 2.40\pm 0.12$ (and a mean of $q_{\rm IR}=2.52 \pm 0.03$ accounting for the lower limits), consistent with both the local value determined from {\em IRAS} and values derived from surveys targeting the high-redshift Universe. By comparing the radio fluxes of our sample measured from both FIRST and NVSS we show that previous results suggesting an increase in the value of $q_{\rm IR}$ from high to low redshift may be the result of resolving out extended emission of the low-redshift sources with relatively high-resolution interferometric data, although AGN contamination could still play a significant role. We also find tentative evidence that the longer wavelength, cooler dust is heated by an evolved stellar population which does not trace the star-formation rate as closely as the shorter wavelength $\ltsim 250~\mu$m emission or the radio emission, supporting suggestions based on detailed models of individual galaxies.

(abridged) We present the first study of the farIR properties of high redshift, radio-selected ULIRGs using deep observations obtained with SPIRE from the Herschel Multi-tiered Extragalactic Survey (HerMES). These galaxies span a large range of 850um fluxes from submillimetre-luminous ~10mJy "SCUBA galaxies -- SMGs" to ~1.5mJy from stacked SCUBA non-detections, thus likely representing a complete distribution of ULIRG spectral energy distributions. From Keck spectroscopic surveys in the Lockman-North field we identified a sample of 31 SMGs and 37 submillimetre-faint, optically-faint radio galaxies (OFRGs), all with radio-inferred IR luminosities >10^12 Lsun. These galaxies were cross-identified with SPIRE 250, 350 and 500um catalogs based on fluxes extracted at 24um positions in the SWIRE survey, yielding a sample of more than half of the galaxies well detected in at least two of the SPIRE bandpasses. By fitting greybody dust models to the SPIRE photometry together with SCUBA 850um measurements, we infer dust temperatures and far-infrared luminosities. The OFRGs detected by SPIRE have median <T_d>= 41+-5 K and the SMGs have <T_d>= 34+-5 K, both in reasonable agreement with previous (pre-Herschel) estimates, reaffirming that the local FIR/radio correlation holds (at least for this subset of high-z ULIRGs) at high redshift. Our observations firstly confirm that a substantial fraction of OFRGs exhibit large infrared luminosities corresponding to SFRs of ~400 Msun/yr. The SPIRE observations secondly confirm the higher dust temperatures for these OFRGs than similarly selected SMGs, consistent with early predictions of the submm-faint radio populations. Our observations also clearly confirm the large infrared luminosities of most SMGs selected with S850um>5 mJy and radio and strong 24um detections, corresponding to SFRs of ~700 Msun/yr.

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

We describe the reduction of data taken with the PACS instrument on board the Herschel Space Observatory in the Science Demonstration Phase of the Herschel-ATLAS (H-ATLAS) survey, specifically data obtained for a 4x4-deg^2 region using Herschel's fast-scan (60 arcsec/s) parallel mode. We describe in detail a pipeline for data reduction using customised procedures within HIPE from data retrieval to the production of science-quality images. We found that the standard procedure for removing Cosmic-Ray glitches also removed parts of bright sources and so implemented an effective two-stage process to minimise these problems. The pronounced 1/f noise is removed from the timelines using 3.4- and 2.5-arcmin boxcar high-pass filters at 100 and 160-um. Empirical measurements of the point-spread function (PSF) are used to determine the encircled energy fraction as a function of aperture size. For the 100- and 160-um bands, the effective PSFs are ~9 and ~13 arcsec (FWHM), and the 90-per-cent encircled energy radii are 13 and 18 arcsec. Astrometric accuracy is good to ~<2 arcsec. The noise in the final maps is correlated between neighbouring pixels and rather higher than advertised prior to launch. For a pair of cross-scans, the 5-sigma point-source sensitivities are 125-165 mJy for 9-13-arcsec radius apertures at 100-um and 150-240 mJy for 13-18-arcsec radius apertures at 160-um.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .

By combining Herschel-SPIRE data with archival Spitzer, HI, and CO maps, we investigate the spatial distribution of gas and dust in the two famous grand-design spirals M99 and M100 in the Virgo cluster. Thanks to the unique resolution and sensitivity of the Herschel-SPIRE photometer, we are for the first time able to measure the distribution and extent of cool, submillimetre (submm)-emitting dust inside and beyond the optical radius. We compare this with the radial variation in both the gas mass and the metallicity. Although we adopt a model-independent, phenomenological approach, our analysis provides important insights. We find the dust extending to at least the optical radius of the galaxy and showing breaks in its radial profiles at similar positions as the stellar distribution. The colour indices f350/f500 and f250/f350 decrease radially consistent with the temperature decreasing with radius. We also find evidence of an increasing gas to dust ratio with radius in the outer regions of both galaxies.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

SPIRE, the Spectral and Photometric Imaging Receiver, is the Herschel Space Observatory's submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 {\mu}m, and an imaging Fourier transform spectrometer (FTS) covering 194-671 {\mu}m (447-1550 GHz). In this paper we describe the initial approach taken to the absolute calibration of the SPIRE instrument using a combination of the emission from the Herschel telescope itself and the modelled continuum emission from solar system objects and other astronomical targets. We present the photometric, spectroscopic and spatial accuracy that is obtainable in data processed through the "standard" pipelines. The overall photometric accuracy at this stage of the mission is estimated as 15% for the photometer and between 15 and 50% for the spectrometer. However, there remain issues with the photometric accuracy of the spectra of low flux sources in the longest wavelength part of the SPIRE spectrometer band. The spectrometer wavelength accuracy is determined to be better than 1/10th of the line FWHM. The astrometric accuracy in SPIRE maps is found to be 2 arcsec when the latest calibration data are used. The photometric calibration of the SPIRE instrument is currently determined by a combination of uncertainties in the model spectra of the astronomical standards and the data processing methods employed for map and spectrum calibration. Improvements in processing techniques and a better understanding of the instrument performance will lead to the final calibration accuracy of SPIRE being determined only by uncertainties in the models of astronomical standards.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

We present a resolved dust analysis of three of the largest angular size spiral galaxies, NGC 4501 and NGC 4567/8, in the Herschel Virgo Cluster Survey (HeViCS) Science Demonstration field. Herschel has unprecedented spatial resolution at far-infrared wavelengths and with the PACS and SPIRE instruments samples both sides of the peak in the far infrared spectral energy distribution (SED).We present maps of dust temperature, dust mass, and gas-to-dust ratio, produced by fitting modified black bodies to the SED for each pixel. We find that the distribution of dust temperature in both systems is in the range ~19 - 22 K and peaks away from the centres of the galaxies. The distribution of dust mass in both systems is symmetrical and exhibits a single peak coincident with the galaxy centres. This Letter provides a first insight into the future analysis possible with a large sample of resolved galaxies to be observed by Herschel.

We use the Science Demonstration Phase data of the Herschel Virgo Cluster Survey to search for dust emission of early-type dwarf galaxies in the central regions of the Virgo Cluster as an alternative way of identifying the interstellar medium.We present the first possible far-infrared detection of cluster early-type dwarf galaxies: VCC781 and VCC951 are detected at the 10 sigma level in the SPIRE 250 micron image. Both detected galaxies have dust masses of the order of 10^5 Msun and average dust temperatures ~20K. The detection rate (less than 1%) is quite high compared to the 1.7% detection rate for Hi emission, considering that dwarfs in the central regions are more Hi deficient. We conclude that the removal of interstellar dust from dwarf galaxies resulting from ram pressure stripping, harassment, or tidal effects must be as efficient as the removal of interstellar gas.

The origin of the far-infrared emission from the nearby radio galaxy M87 remains a matter of debate. Some studies find evidence of a far-infrared excess due to thermal dust emission, whereas others propose that the far-infrared emission can be explained by synchrotron emission without the need for an additional dust emission component. We present Herschel PACS and SPIRE observations of M87, taken as part of the science demonstration phase observations of the Herschel Virgo Cluster Survey. We compare these data with a synchrotron model based on mid-infrared, far-infrared, submm and radio data from the literature to investigate the origin of the far-infrared emission. Both the integrated SED and the Herschel surface brightness maps are adequately explained by synchrotron emission. At odds with previous claims, we find no evidence of a diffuse dust component in M87, which is not unexpected in the harsh X-ray environment of this radio galaxy sitting at the core of the Virgo Cluster.

We present the dust properties of a small sample of Virgo cluster dwarf galaxies drawn from the science demonstration phase data set of the Herschel Virgo Cluster Survey. These galaxies have low metallicities (7.8 < 12 + log(O/H) < 8.3) and star-formation rates < 10^{-1} M_{sun}/yr. We measure the spectral energy distribution (SED) from 100 to 500 um and derive dust temperatures and dust masses. The SEDs are fitted by a cool component of temperature T < 20 K, implying dust masses around 10^{5} M_{sun} and dust-to-gas ratios D within the range 10^{-3}-10^{-2}. The completion of the full survey will yield a larger set of galaxies, which will provide more stringent constraints on the dust content of star-forming dwarf galaxies.

Passive early-type galaxies (ETGs) provide an ideal laboratory for studying the interplay between dust formation around evolved stars and its subsequent destruction in a hot gas. Using Spitzer-IRS and Herschel data we compare the dust production rate in the envelopes of evolved AGB stars with a constraint on the total dust mass. Early-type galaxies which appear to be truly passively evolving are not detected by Herschel. We thus derive a distance independent upper limit to the dust grain survival time in the hostile environment of ETGs of < 46 +/- 25 Myr for amorphous silicate grains. This implies that ETGs which are detected at far-infrared wavelengths have acquired a cool dusty medium via interaction. Given likely time-scales for ram-pressure stripping, this also implies that only galaxies with dust in a cool (atomic) medium can release dust into the intra-cluster medium.

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

By combining Herschel-SPIRE observations obtained as part of the Herschel Virgo Cluster Survey with 21 cm HI data from the literature, we investigate the role of the cluster environment on the dust content of Virgo spiral galaxies.We show for the first time that the extent of the dust disk is significantly reduced in HI-deficient galaxies, following remarkably well the observed 'truncation' of the HI disk. The ratio of the submillimetre-to- optical diameter correlates with the HI-deficiency, suggesting that the cluster environment is able to strip dust as well as gas. These results provide important insights not only into the evolution of cluster galaxies but also into the metal enrichment of the intra-cluster medium.

We describe the Herschel Virgo Cluster Survey (HeViCS) and the first data obtained as part of the Science Demonstration Phase (SDP). The data cover a central 4x4 sq deg region of the cluster. We use SPIRE and PACS photometry data to produce 100, 160, 250, 350 and 500 micron luminosity functions (LFs) for optically bright galaxies that are selected at 500 micron and detected in all bands. We compare these LFs with those previously derived using IRAS, BLAST and Herschel-ATLAS data. The Virgo Cluster LFs do not have the large numbers of faint galaxies or examples of very luminous galaxies seen previously in surveys covering less dense environments.

Aims: We aim to demonstrate that the Herschel ATLAS (H-ATLAS) is suitable for a blind and unbiased survey for debris disks by identifying candidate debris disks associated with main sequence stars in the initial science demonstration field of the survey. We show that H-ATLAS reveals a population of far-infrared/sub-mm sources that are associated with stars or star-like objects on the SDSS main-sequence locus. We validate our approach by comparing the properties of the most likely candidate disks to those of the known population. Methods: We use a photometric selection technique to identify main sequence stars in the SDSS DR7 catalogue and a Bayesian Likelihood Ratio method to identify H-ATLAS catalogue sources associated with these main sequence stars. Following this photometric selection we apply distance cuts to identify the most likely candidate debris disks and rule out the presence of contaminating galaxies using UKIDSS LAS K-band images. Results: We identify 78 H-ATLAS sources associated with SDSS point sources on the main-sequence locus, of which two are the most likely debris disk candidates: H-ATLAS J090315.8 and H-ATLAS J090240.2. We show that they are plausible candidates by comparing their properties to the known population of debris disks. Our initial results indicate that bright debris disks are rare, with only 2 candidates identified in a search sample of 851 stars. We also show that H-ATLAS can derive useful upper limits for debris disks associated with Hipparcos stars in the field and outline the future prospects for our debris disk search programme.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

We present measurements of the angular correlation function of galaxies selected from the first field of the H-ATLAS survey. Careful removal of the background from galactic cirrus is essential, and currently dominates the uncertainty in our measurements. For our 250 micron-selected sample we detect no significant clustering, consistent with the expectation that the 250 micron-selected sources are mostly normal galaxies at z<~ 1. For our 350 micron and 500 micron-selected samples we detect relatively strong clustering with correlation amplitudes A of 0.2 and 1.2 at 1', but with relatively large uncertainties. For samples which preferentially select high redshift galaxies at z~2-3 we detect significant strong clustering, leading to an estimate of r_0 ~ 7-11 h^{-1} Mpc. The slope of our clustering measurements is very steep, delta~2. The measurements are consistent with the idea that sub-mm sources consist of a low redshift population of normal galaxies and a high redshift population of highly clustered star-bursting galaxies.

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

We present Herschel-SPIRE observations of the perturbed galaxy NGC4438 in the Virgo cluster. These images reveal the presence of extra-planar dust up to ~4-5 kpc away from the galaxy's disk. The dust closely follows the distribution of the stripped atomic and molecular hydrogen, supporting the idea that gas and dust are perturbed in a similar fashion by the cluster environment. Interestingly, the extra-planar dust lacks a warm temperature component when compared to the material still present in the disk, explaining why it was missed by previous far-infrared investigations. Our study provides evidence for dust stripping in clusters of galaxies and illustrates the potential of Herschel data for our understanding of environmental effects on galaxy evolution.

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.

We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 $\sigma$, we find an average dust temperature of $(28 \pm 8)$K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3$\sigma$, we find an average redshift of $2.2 \pm 0.6$.

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500 microns when completed, reaching flux limits (5 sigma) from 32 to 145mJy. We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500 microns. Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to model catalogs and then applied to the raw observational counts. Results. We find a steep rise in the number counts at flux levels of 100-200mJy in all three SPIRE bands, consistent with results from BLAST. The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift.

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$\mu$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $\gamma$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

We present Herschel SPIRE and PACS photometeric observations of the low metallicity (Z ~ 0.35 solar) nearby dwarf galaxy, NGC 1705, in six wavelength bands as part of the Dwarf Galaxy Survey guaranteed time Herschel Key Program. We confirm the presence of two dominant circumnuclear IR-bright regions surrounding the central super star cluster that had been previously noted at mid-IR wavelengths and in the sub-mm by LABOCA. On constructing a global spectral energy distribution using the SPIRE and PACS photometry, in conjunction with archival IR measurements, we note the presence of an excess at sub-mm wavelengths. This excess suggests the presence of a significant cold dust component within NGC 1705 and was modeled as an additional cold component in the SED. Although alternative explanations for the sub-mm excess beyond 350 microns, such as changes to the dust emissivity cannot be ruled out, the most likely explanation for the observed submillimetre excess is that of an additional cold dust component.

The standard method of mapping the interstellar medium in a galaxy, by observing the molecular gas in the CO 1-0 line and the atomic gas in the 21-cm line, is largely limited with current telescopes to galaxies in the nearby universe. In this letter, we use SPIRE observations of the galaxies M99 and M100 to explore the alternative approach of mapping the interstellar medium using the continuum emission from the dust. We have compared the methods by measuring the relationship between the star-formation rate and the surface density of gas in the galaxies. We find the two methods give relationships with a similar dispersion, confirming that observing the continuum emission from the dust is a promising method of mapping the interstellar medium in galaxies.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

With appropriate spatial resolution, images of spiral galaxies in thermal infrared (~10 micron and beyond) often reveal a bright central component, distinct from the stellar bulge, superimposed on a disk with prominent spiral arms. ISO and Spitzer studies have shown that much of the scatter in the mid-infrared colors of spiral galaxies is related to changes in the relative importance of these two components, rather than to other modifications, such as the morphological type or star formation rate, that affect the properties of the galaxy as a whole. With the Herschel imaging capability from 70 to 500 micron, we revisit this two-component approach at longer wavelengths, to see if it still provides a working description of the brightness distribution of galaxies, and to determine its implications on the interpretation of global far-infrared properties of galaxies.

We use Herschel PACS and SPIRE observations of the edge-on spiral galaxy UGC 4754, taken as part of the H-ATLAS SDP observations, to investigate the dust energy balance in this galaxy. We build detailed SKIRT radiative models based on SDSS and UKIDSS maps and use these models to predict the far-infrared emission. We find that our radiative transfer model underestimates the observed FIR emission by a factor of two to three. Similar discrepancies have been found for other edge-on spiral galaxies based on IRAS, ISO, and SCUBA data. Thanks to the good sampling of the SED at FIR wavelengths, we can rule out an underestimation of the FIR emissivity as the cause for this discrepancy. Instead we support highly obscured star formation that contributes little to the optical extinction as a more probable explanation.

We present the first Herschel PACS and SPIRE images of the low-metallicity galaxy NGC6822 observed from 70 to 500 mu and clearly resolve the HII regions with PACS and SPIRE. We find that the ratio 250/500 is dependent on the 24 mu surface brightness in NGC6822, which would locally link the heating processes of the coldest phases of dust in the ISM to the star formation activity. We model the SEDs of some regions HII regions and less active regions across the galaxy and find that the SEDs of HII regions show warmer ranges of dust temperatures. We derive very high dust masses when graphite is used in our model to describe carbon dust. Using amorphous carbon, instead, requires less dust mass to account for submm emission due to its lower emissivity properties. This indicates that SED models including Herschel constraints may require different dust properties than commonly used.

We present the observations of the starburst galaxy M82 taken with the Herschel SPIRE Fourier Transform Spectrometer. The spectrum (194-671 {\mu}m) shows a prominent CO rotational ladder from J = 4-3 to 13-12 emitted by the central region of M82. The fundamental properties of the gas are well constrained by the high J lines observed for the first time. Radiative transfer modeling of these high-S/N 12CO and 13CO lines strongly indicates a very warm molecular gas component at ~500 K and pressure of ~3x10^6 K cm^-3, in good agreement with the H_2 rotational lines measurements from Spitzer and ISO. We suggest that this warm gas is heated by dissipation of turbulence in the interstellar medium (ISM) rather than X-rays or UV flux from the straburst. This paper illustrates the promise of the SPIRE FTS for the study of the ISM of nearby galaxies.

We use Herschel Space Observatory data to place observational constraints on the peak and Rayleigh-Jeans slope of dust emission observed at 70-500 microns in the nearby spiral galaxy M81. We find that the ratios of wave bands between 160 and 500 microns are primarily dependent on radius but that the ratio of 70 to 160 micron emission shows no clear dependence on surface brightness or radius. These results along with analyses of the spectral energy distributions imply that the 160-500 micron emission traces 15-30 K dust heated by evolved stars in the bulge and disc whereas the 70 micron emission includes dust heated by the active galactic nucleus and young stars in star forming regions.

We present infrared colours (in the 25-500 mic spectral range) and UV to radio continuum spectral energy distributions of a sample of 51 nearby galaxies observed with SPIRE on Herschel. The observed sample includes all morphological classes, from quiescent ellipticals to active starbursts. Active galaxies have warmer colour temperatures than normal spirals. In ellipticals hosting a radio galaxy, the far-infrared (FIR) emission is dominated bynthe synchrotron nuclear emission. The colour temperature of the cold dust is higher in quiescent E-S0a than in star-forming systems probably because of the different nature of their dust heating sources (evolved stellar populations, X-ray, fast electrons) and dust grain properties. In contrast to the colour temperature of the warm dust, the f350/f500 index sensitive to the cold dust decreases with star formation and increases with metallicity, suggesting an overabundance of cold dust or an emissivity parameter beta<2 in low metallicity, active systems.

M82 is a unique representative of a whole class of galaxies, starbursts with superwinds, in the Very Nearby Galaxy Survey with Herschel. In addition, its interaction with the M81 group has stripped a significant portion of its interstellar medium from its disk. SPIRE maps now afford better characterization of the far-infrared emission from cool dust outside the disk, and sketch a far more complete picture of its mass distribution and energetics than previously possible. They show emission coincident in projection with the starburst wind and in a large halo, much more extended than the PAH band emission seen with Spitzer. Some complex substructures coincide with the brightest PAH filaments, and others with tidal streams seen in atomic hydrogen. We subtract the far-infrared emission of the starburst and underlying disk from the maps, and derive spatially-resolved far-infrared colors for the wind and halo. We interpret the results in terms of dust mass, dust temperature, and global physical conditions. In particular, we examine variations in the dust physical properties as a function of distance from the center and the wind polar axis, and conclude that more than two thirds of the extraplanar dust has been removed by tidal interaction, and not entrained by the starburst wind.

We present Herschel-SPIRE observations at 250-500um of the giant elliptical galaxy M86 and examine the distribution of the resolved cold dust emission and its relation with other galactic tracers. The SPIRE images reveal three dust components: emission from the central region; a dust lane extending north-south; and a bright emission feature 10kpc to the south-east. We estimate that approximately 10^6 solar masses of dust is spatially coincident with atomic and ionized hydrogen, originating from stripped material from the nearby spiral NGC4438 due to recent tidal interactions with M86. The gas-to-dust ratio of the cold gas component ranges from ~20-80. We discuss the different heating mechanisms for the dust features.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

We present our recent results on the properties of the outskirts of disk galaxies. In particular, we focus on spiral galaxies with stellar disk truncations in their radial surface brightness profiles. Using SDSS, UDF and GOODS data we show how the position of the break (i.e., a direct estimator of the size of the stellar disk) evolves with time since z~1. Our findings agree with an evolution on the radial position of the break by a factor of 1.3+/-0.1 in the last 8 Gyr for galaxies with similar stellar masses. We also present radial color gradients and how they evolve with time. At all redshift we find a radial inside-out bluing reaching a minimum at the position of the break radius, this minimum is followed by a reddening outwards. Our results constrain several galaxy disk formation models and favour a scenario where stars are formed inside the break radius and are relocated in the outskirts of galaxies through secular processes.

The Herschel Reference Survey is a guaranteed time Herschel key project and will be a benchmark study of dust in the nearby universe. The survey will complement a number of other Herschel key projects including large cosmological surveys that trace dust in the distant universe. We will use Herschel to produce images of a statistically-complete sample of 323 galaxies at 250, 350 and 500 micron. The sample is volume-limited, containing sources with distances between 15 and 25 Mpc and flux limits in the K-band to minimize the selection effects associated with dust and with young high-mass stars and to introduce a selection in stellar mass. The sample spans the whole range of morphological types (ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster) and as such will be useful for other purposes than our own. We plan to use the survey to investigate (i) the dust content of galaxies as a function of Hubble type, stellar mass and environment, (ii) the connection between the dust content and composition and the other phases of the interstellar medium and (iii) the origin and evolution of dust in galaxies. In this paper, we describe the goals of the survey, the details of the sample and some of the auxiliary observing programs that we have started to collect complementary data. We also use the available multi-frequency data to carry out an analysis of the statistical properties of the sample.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

We have conducted a deep multi-wavelength analysis (0.15-160 mum) to study the outer region of the nearby galaxy M94. We show that the non-optical data support the idea that the outskirts of this galaxy is not formed by a closed stellar ring (as traditionally claimed in the literature) but by a spiral arm structure. In this sense, M94 is a good example of a Type III (anti-truncated) disk galaxy having a very bright outer disk. The outer disk of this galaxy contains ~23% of the total stellar mass budget of the galaxy and contributes ~10% of the new stars created showing that this region of the galaxy is active. In fact, the specific star formation rate of the outer disk (~0.012 Gyr^{-1}) is a factor of ~2 larger (i.e. the star formation is more efficient per unit stellar mass) than in the inner disk. We have explored different scenarios to explain the enhanced star formation in the outer disk. We find that the inner disk (if considered as an oval distortion) can dynamically create a spiral arm structure in the outer disk which triggers the observed relatively high star formation rate as well as an inner ring similar to what is found in this galaxy.

Within the hierarchical framework for galaxy formation, merging and tidal interactions are expected to shape large galaxies to this day. While major mergers are quite rare at present, minor mergers and satellite disruptions - which result in stellar streams - should be common, and are indeed seen in both the Milky Way and the Andromeda Galaxy. As a pilot study, we have carried out ultra-deep, wide-field imaging of some spiral galaxies in the Local Volume, which has revealed external views of such stellar tidal streams at unprecedented detail, with data taken at small robotic telescopes (0.1-0.5-meter) that provide exquisite surface brightness sensitivity. The goal of this project is to undertake the first systematic and comprehensive imaging survey of stellar tidal streams, from a sample of ~50 nearby Milky-Way-like spiral galaxies within 15 Mpc, that features a surface brightness sensitivity of ~ 30 mag/arcsec^2 The survey will result in estimates of the incidence, size/geometry and stellar luminosity/mass distribution of such streams. This will not only put our Milky Way and M31 in context but, for the first time, also provide an extensive statistical basis for comparison with state-of-the-art, self-consistent cosmological simulations of this phenomenon.

We report the discovery of a giant, loop-like stellar structure around the edge-on spiral galaxy NGC 4013. This arcing feature extends 6 arcmin (~26 kpc in projected distance) northeast from the center and 3 arcmin (~=12 kpc) from the disk plane; likely related features are also apparent on the southwest side of the disk, extending to 4 arcmin (~17 kpc). The detection of this low surface-brightness muR= 27.0+0.3-0.2 mag/sqarcsec) structure is independently confirmed in three separate datasets from three different telescopes. Although its true three dimensional geometry is unknown, the sky- projected morphology of this structure displays a match with the theoretical predictions for the edge-on, projected view of a stellar tidal streams of a dwarf satellite moving in a low inclined (~25deg), nearly circular orbit. Using the recent model of the Monoceros tidal stream in the Milky Way by Penarrubia et al. as template, we find that the progenitor system may have been a galaxy with an initial mass 6*10^8 Msun, of which current position and final fate is unknown. According to this simulation, the tidal stream may be approximately ~2.8 Gyr of age. Our results demonstrate that NGC 4013, previously considered a prototypical isolated disk galaxy in spite of having one of the most prominent HI warps detected thus far, may have in fact suffered a recent minor merger. This discovery highlights that undisturbed disks at high surface brightness levels in the optical but warped in HI maps may in fact reveal complex signatures of recent accretion events in deep photometric surveys.

We have explored radial color and stellar surface mass density profiles for a sample of 85 late-type spiral galaxies with deep (down to ~27 mag arcsec^-2) SDSS g'- and r'-band surface brightness profiles. About 90% of the light profiles have been classified as broken exponentials, exhibiting either truncations (Type II galaxies) or antitruncations (Type III galaxies). The color profiles of Type II galaxies show a "U shape" with a minimum of (g' - r') = 0.47 +- 0.02 mag at the break radius. Around the break radius, Type III galaxies have a plateau region with a color of (g' - r') = 0.57 +- 0.02. Using the color to calculate the stellar surface mass density profiles reveals a surprising result. The breaks, well established in the light profiles of the truncated galaxies, are almost gone, and the mass profiles resemble now those of the pure exponential (Type I) galaxies. This result suggests that the origin of the break in Type II galaxies is more likely due to a radial change in stellar population than being associated to an actual drop in the distribution of mass. Type III galaxies, however, seem to preserve their shape in the stellar mass density profiles. We find that the stellar surface mass density at the break for truncated galaxies is 13.6 +- 1.6 Msun pc^-2 and for the antitruncated ones is 9.9 +- 1.3 Msun pc^-2 . We estimate that the fraction of stellar mass outside the break radius is ~15% for truncated galaxies and ~9% for antitruncated galaxies.

We present an extragalactic perspective of an extended stellar tidal stream wrapping around the edge-on, spiral galaxy NGC 5907. Our deep images reveal for the first time a large scale complex of arcing loops that is an excellent example of how a low-mass satellite accretion can produce an interwoven, rosette-like structure of debris dispersed in the halo of its host galaxy. The existence of this structure, which has probably formed and survived for several Gigayears, confirms that halos of spiral galaxies in the Local Universe may still contain a significant number of galactic fossils from their hierarchical formation. To examine the validity of the external accretion scenario, we present N-body simulations of the tidal disruption of a dwarf galaxy-like system in a disk galaxy plus dark halo potential that demonstrate that most of the observed tidal features observed in NGC 5907 can be explained by a single accretion event. Unfortunately, with no kinematic data and only the projected geometry of the stream as constraint, the parameters of our model are considerably degenerate and, for now, must be considered illustrative only. Interestingly, NGC 5907 has long been considered a prototypical example of a warped spiral in relative isolation. The presence of an extended tidal stream challenges this picture and suggests that the gravitational perturbations induced by the stream progenitor may be the cause for the warp. The detection of an old, complex tidal stream in a nearby galaxy with rather modest instrumentation points to the viability of surveys to find extragalactic tidal substructures around spiral galaxies in the Local Volume (< 15 Mpc) -- with the prospect of obtaining a census with enough statistical significance to be compared with cosmological simulations.