We present Herschel SPIRE-FTS observations of Arp~220, a nearby ULIRG. The FTS continuously covers 190 -- 670 microns, providing a good measurement of the continuum and detection of several molecular and atomic species. We detect luminous CO (J = 4-3 to 13-12) and water ladders with comparable total luminosity; very high-J HCN absorption; OH+, H2O+, and HF in absorption; and CI and NII. Modeling of the continuum yields warm dust, with T = 66 K, and an unusually large optical depth of ~5 at 100 microns. Non-LTE modeling of the CO shows two temperature components: cold molecular gas at T ~ 50 K and warm molecular gas at T ~1350 K. The mass of the warm gas is 10% of the cold gas, but dominates the luminosity of the CO ladder. The temperature of the warm gas is in excellent agreement with H2 rotational lines. At 1350 K, H2 dominates the cooling (~20 L_sun/M_sun) in the ISM compared to CO (~0.4 L_sun/M_sun). We found that only a non-ionizing source such as the mechanical energy from supernovae and stellar winds can excite the warm gas and satisfy the energy budget of ~20 L_sun/M_sun. We detect a massive molecular outflow in Arp 220 from the analysis of strong P-Cygni line profiles observed in OH+, H2O+, and H2O. The outflow has a mass > 10^{7} M_sun and is bound to the nuclei with velocity < 250 km/s. The large column densities observed for these molecular ions strongly favor the existence of an X-ray luminous AGN (10^{44} ergs/s) in Arp 220.

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.

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.

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.

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.

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.

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

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

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.

We present a catalog of optical spectroscopic identifications of sources detected by Spitzer at 3.6 or 24 micron down to 10 and 280 microJy, respectively, in the SWIRE/XMM-Newton/ELAIS-S1 field and classified via line width analysis and diagnostic diagrams. A total of 1376 sources down to R~24.2 mag have been identified (1362 detected at 3.6 micron, 419 at 24 micron, and 405 at both) by low-resolution optical spectroscopy carried out with FORS2, VIMOS, and EFOSC2 at the Very Large Telescope and 3.6m ESO telescopes. The spectroscopic campaigns have been carried out over the central 0.6 square degrees area of ELAIS-S1 which, in particular, has also been observed by XMM-Newton and Chandra. We find the first direct optical spectroscopic evidence that the fraction of active galactic nuclei (AGN; mostly AGN2) increases with increasing F(24 micron)/F(R) ratio, reaching values of 70(+/-20)% in the range 316<F(24 micron)F(R)<1000. We present an IRAC-MIPS color-color diagram able to separate AGN1 from obscured AGN2 candidates. After having corrected for the spectroscopic incompleteness of our sample, it results that the AGN fraction at F(24 micron)=0.8 mJy is ~22(+/-7)% and decreases slowly to ~19(+/-5)% down to F(24 micron)=0.3 mJy.

We present optical identifications and a multi-band catalogue of a sample of 478 X-ray sources in the XMM and Chandra surveys of the central 0.6 deg^2 of the ELAIS-S1 field. The optical/infrared counterpart of each X-ray source was identified using R and IRAC 3.6 um bands. This method was complemented by the precise positions obtained through Chandra observations. Approximately 94% of the counterparts are detected in the R band, while the remaining are blank fields in the optical down to R~24.5, but have a near-infrared counterpart detected by IRAC within 6 arcsec from the XMM centroid. The multi-band catalogue contains photometry in ten photometric bands (B to the MIPS 24 um). We determined redshift and classification for 237 sources (~50% of the sample) brighter than R=24. We classified 47% of the sources with spectroscopic redshift as broad-line active galactic nuclei (BL AGNs) with z=0.1-3.5, while sources without broad-lines are about 46% of the spectroscopic sample and are found up to z=2.6. We identified 11 type 2 QSOs among the sources with X/O>8, with z=0.9-2.6, high 2-10 keV luminosity (log(L2-10keV)>=43.8 erg/s) and hard X-ray colors suggesting large absorbing columns at the rest frame (logN_H up to 23.6 cm^-2). BL AGNs show on average blue optical-to-near-infrared colors, softer X-ray colors and X-ray-to-optical colors typical of optically selected AGNs. Conversely, narrow-line sources show redder optical colors, harder X-ray flux ratio and span a wider range of X-ray-to-optical colors. On average the SEDs of high-luminosity BL AGNs resemble the power-law typical of unobscured AGNs. The SEDs of NOT BL AGNs are dominated by the galaxy emission in the optical/near-infrared, and show a rise in the mid-infrared which suggests the presence of an obscured active nucleus.

(abridged) The ESIS survey is the optical follow up of the SWIRE/Spitzer in the ELAIS-S1 region of the sky. In the era of observational cosmology, the main efforts are focused on the study of galaxy evolution and its environmental dependence. Wide area, multiwavelength, extragalactic surveys are needed in order to probe sufficiently large volumes, minimize cosmic variance and find significant numbers of rare objects. We present VIMOS I and z band imaging belonging to the ESIS survey. A total of ~4 deg2 were targeted in I and ~1 deg2 in z. More than 300000 galaxies have been detected in the I band and ~50000 in the z band. Object coordinates are defined within an uncertainty of ~0.2 arcsec r.m.s., with respect to GSC 2.2. We reach a 90% average completeness at 23.1 and 22.5 mag (Vega) in the I and z bands, respectively. On the basis of IRAC colors, we identified galaxies having the 1.6 um stellar peak shifted to z=1-3. The new I, z band data provide reliable constraints to avoid low-redshift interlopers and reinforce this selection. Roughly 1000 galaxies between z=2-3 were identified over the ESIS ~4 deg2, at the SWIRE 5.8 um depth (25.8 uJy at 3sigma). These are the best galaxy candidates to dominate the massive tail (M>1e11 Msun) of the z>2 mass function.

We present the broad-band Spectral Energy Distributions (SEDs) of the largest available highly (72%) complete spectroscopic sample of mid-infrared (MIR) selected galaxies and AGN at intermediate redshift. The sample contains 203 extragalactic sources from the 15-micron survey in the ELAIS-SWIRE field S1, all with measured spectroscopic redshift. Most of these sources have full multi-wavelength coverage from the far-UV to the far-infrared and lie in the redshift range 0.1<z<1.3. Due to its size, this sample allows us for the first time to characterise the spectral properties of the sources responsible for the strong evolution observed in the MIR. Based on SED-fitting technique we have classified the MIR sources, identifying AGN signatures in about 50% of them. This fraction is significantly higher than that derived from optical spectroscopy (~29%) and is due in particular to the identification of AGN activity in objects spectroscopically classified as galaxies. It is likely that in most of our objects, the AGN is either obscured or of low-luminosity, and thus it does not dominate the energetic output at any wavelength, except in the MIR, showing up just in the range where the host galaxy SED has a minimum. The fraction of AGN strongly depends on the flux density, with that derived through the SED-fitting being about 20% at S(15)~0.5-1 mJy and gradually increasing up to 100% at S(15)>10 mJy, while that obtained from optical spectroscopy never being >30%, even at the higher flux densities. The results of this work will be very useful for updating all the models aimed at interpreting the deep infrared survey data and, in particular, for constraining the nature and the role of dust-obscured systems in the intermediate/high-redshift Universe.

(Abridged) We selected high-z massive galaxies at 5.8 microns, in the SWIRE ELAIS-S1 field (1 sq. deg.). Galaxies with the 1.6 microns stellar peak redshifted into the IRAC bands (z~1-3, called ``IR-peakers'') were identified. Stellar masses were derived by means of spectro-photometric fitting and used to compute the stellar mass function (MF) at z=1-2 and 2-3. A parametric fit to the MF was performed, based on a Bayesian formalism, and the stellar mass density of massive galaxies above z=2 determined. We present the first systematic study of the very-massive tail of the galaxy stellar mass function at high redshift. A total of 326 sources were selected. The majority of these galaxies have stellar masses in excess of 1e11 Msun and lie at z>1.5. The availability of mid-IR data turned out to be a valuable tool to constrain the contribution of young stars to galaxy SEDs, and thus their M(stars)/L ratio. The influence of near-IR data and of the chosen stellar library on the SED fitting are also discussed. A significant evolution is found not only for galaxies with M~1e11 Msun, but also in the highest mass bins considered. The comoving number density of these galaxies was lower by more than a factor of 10 at z=2-3, with respect to the local estimate. SWIRE 5.8 micron peakers more massive than 1.6x1e11 Msun provide 30-50% of the total stellar mass density in galaxies at z=2-3.

Context: The counts of galaxies and AGN in the mid infra-red (MIR) bands are important instruments for studying their cosmological evolution. However, the classic spectral line ratios techniques can become misleading when trying to properly separate AGN from starbursts or even from apparently normal galaxies. Aims: We use X-ray band observations to discriminate AGN activity in previously classified MIR-selected starburst galaxies and to derive updated AGN1 and (Compton thin) AGN2 counts at 15 um. Methods: XMM observations of the ELAIS-S1 15um sample down to flux limits ~2x10^-15 erg cm^-2 s^-1 (2-10 keV band) were used. We classified as AGN all those MIR sources with a unabsorbed 2-10 keV X-ray luminosity higher that ~10^42 erg/s. Results: We find that at least about 13(+/-6) per cent of the previously classified starburst galaxies harbor an AGN. According to these figures, we provide an updated estimate of the counts of AGN1 and (Compton thin) AGN2 at 15 um. It turns out that at least 24% of the extragalactic sources brighter than 0.6 my at 15 um are AGN (~13% contribution to the extragalactic background produced at fluxes brighter than 0.6 mJy).

We have surveyed with XMM-Newton the central ~0.6 deg2 region of the ELAIS-S1 field down to flux limits of ~5.5X10-16 cgs (0.5-2 keV, S band), ~2X10-15 cgs (2-10 keV, H band), and ~4X10-15 cgs (5-10 keV, HH band). We detect a total of 478 sources, 395 and 205 of which detected in the S and H bands respectively. We identified 7 clearly extended sources and estimated their redshift through X-ray spectral fits with thermal models. In four cases the redshift is consistent with z=0.4. We have computed the angular correlation function of the sources in the S and H bands, finding best fit correlation angles theta_0=5.2+/-3.8 arcsec and theta_0=12.8+/-7.8 arcsec respectively. A rough estimate of the present-day correlation length r_0 can be obtained inverting the Limber equation and assuming an appropriate redshift distribution dN/dz. The results range between 12.8 and 9.8 h-1 Mpc in the S band and between 17.9 and 13.4 h-1 Mpc in the H band, with 30-40% statistical errors, assuming either smooth redshift distributions or redshift distributions with spikes accounting for the presence of a structure at z=0.4. The relative density of the S band sources is higher near the clusters and groups at z~0.4 and extends toward East and toward South/West. This suggests that the structure is complex, with a size comparable to the full XMM-Newton field. Conversely, the highest relative source densities of the H band sources are located in the central-west region of the field.

We have determined the cosmological evolution of the density of AGN and of their Nh distribution as a function of the un-absorbed 2-10 keV luminosity up to redshift 4. We used the HELLAS2XMM sample combined with other published catalogs, yielding a total of 508 AGN. Our best fit is obtained with a luminosity-dependent density evolution (LDDE) model where low luminosity (Lx~10^43 erg s^-1) AGN peak at z~0.7, while high luminosity AGN (Lx>10^45 erg s^-1) peak at z~2.0. A pure luminosity evolution model (PLE) can instead be rejected. There is evidence that the fraction of absorbed (Nh>10^22 cm^-2) AGN decreases with the intrinsic X-ray luminosity, and increases with the redshift. Our best fit solution provides a good fit to the observed counts, the cosmic X-ray background, and to the observed fraction of absorbed AGN as a function of the flux in the 10^-15<S(2-10)<10^-10 erg s^-1 cm^-2 range. We find that the absorbed, high luminosity (Lx>10^44 erg s^-1) AGN have a density of 267 deg^-2 at fluxes S(2-10)>10^-15 erg s^-1 cm^-2. Using these results, we estimate a density of supermassive black holes in the local Universe of 3.2 (h^2_70) x 10^5 M_sol Mpc^-3, which is consistent with the recent measurements of the black hole mass function in the local galaxies.

The spectroscopic analysis of 117 serendipitous sources in the HELLAS2XMM 1df (1 degree field) survey is described. Of these, 106 sources, of which 86% have a spectroscopic redshift, are used to evaluate the fraction of X-ray absorbed (log NH>22) Active Galactic Nuclei (AGN) in the 2--10 keV flux range 0.8-20E-14 erg/cm2/s. This fraction turns out lower than what is predicted by two well known Cosmic X-Ray Background synthesis models, and the discrepancy is significant at the 99.999% level. This result consolidates the findings recently obtained by other authors. In the flux interval explored, the data are consistent with an intrinsic distribution of the absorbing columns (flat per decade above logNH>21) independent of luminosity and redshift, together with an AGN luminosity function evolving purely in luminosity. It is shown that, on the other hand, extrapolation to lower fluxes fails to reproduce the results inferred from the Chandra Deep Field North survey. It is found that about 40% of the high luminosity sources in our sample have best fit logNH>22, and the surface density of these X-ray obscured QSOs can then be estimated at about 48 per square degree, at the flux limit of ~1E-14} erg/cm2/s of the HELLAS2XMM 1df survey. As a side issue, 5 or 6 out of 60 sources, that is about 10%, identified with broad line AGN, turn out to be affected by logNH>22 absorption.